Commit graph

210 commits

Author SHA1 Message Date
Henry Weller
cf9ed8cdbb LTSReactingParcelFoam: Minor improvement 2015-06-27 12:30:10 +01:00
Henry Weller
869af5bf17 reactingTwoPhaseEulerFoam: Removed superfluous trSubDeltaT field 2015-06-26 22:05:39 +01:00
Henry Weller
d314ed3b4a LTS, MULES and sub-cycling: Improved the handling of the rSubDeltaT field 2015-06-26 21:58:27 +01:00
Henry Weller
3025df4e48 interFoam: Change construction order 2015-06-26 18:54:37 +01:00
Henry Weller
e0d7809ca5 reactingTwoPhaseEulerFoam: Added experimental run-time selectable LTS support
Select LTS via the ddtScheme:

    ddtSchemes
    {
        default         localEuler rDeltaT;
    }

The LTS algorithm is currently controlled with the standard settings in
controlDict, e.g.:

    maxCo           0.5;
    maxDeltaT       2e-8;

with the addition of the optional rDeltaT smoothing coefficient:

    rDeltaTSmoothingCoeff 0.02;

which defaults to 0.02.

    ddtSchemes
    {
        default         localEuler rDeltaT;
    }
2015-06-26 18:52:23 +01:00
Henry Weller
2ac4a4e84c interFoam family: Added run-time selectable LTS support
LTS is selected by the ddt scheme e.g. in the
tutorials/multiphase/interFoam/ras/DTCHull case:

ddtSchemes
{
    default         localEuler rDeltaT;
}

LTSInterFoam is no longer needed now that interFoam includes LTS
support.
2015-06-26 18:32:20 +01:00
Henry Weller
f67c9ab9bd reactingTwoPhaseEulerFoam: Change the implicit handling of phase-pressure and dispersion
to support any number of phases
2015-06-26 15:15:10 +01:00
Henry Weller
0f94c22141 reactingTwoPhaseEulerFoam: Minor reorganization 2015-06-25 22:43:22 +01:00
Henry Weller
acfe44eac4 reactingTwoPhaseEulerFoam: Add fvOption handling to the continuity error
correction in MovingPhaseModel<BasePhaseModel>::correct()
2015-06-25 22:29:08 +01:00
Henry Weller
64d8d90f03 sampledPlane: Correct handling of coordinate system specification
Patches provided by Timm Severin
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1764
2015-06-25 22:04:00 +01:00
Henry Weller
5feec38881 twoPhaseEulerFoam, reactingTwoPhaseEulerFoam: Removed duplicate files 2015-06-25 20:27:09 +01:00
Henry Weller
ed402f7a0c smoluchowskiJumpTFvPatchScalarField: Updated use of basicThermo 2015-06-25 19:28:16 +01:00
Henry Weller
27d956d278 reactingTwoPhaseEulerFoam: Construct MRF and fvOptions in phaseSystem 2015-06-25 19:22:47 +01:00
Henry Weller
a3c6e30bc7 reactingTwoPhaseEulerFoam: Change the sign of the compressibility correction field
to correspond to the phase dilatation rate.
2015-06-25 16:37:53 +01:00
Henry Weller
d83df30b98 reactingTwoPhaseEulerFoam: Generalize the handling of the dilatation rate
to support any number of phases
2015-06-25 16:10:07 +01:00
Henry Weller
fa6902fde0 twoPhaseEulerFoam: Change the implicit particle-pressure and turbulence dispersion
to be phase-symmetric so that the results are independent of which
phase-fraction is solved.
2015-06-25 16:08:21 +01:00
Henry Weller
4d1372c358 Update headers 2015-06-24 10:45:42 +01:00
Henry Weller
61ba5f0268 Use basicThermo::dictName rather than hard-coding "thermophysicalProperties" 2015-06-24 10:44:57 +01:00
Henry Weller
7a624425cb utilities/postProcessing: Automate the selection of incompressible/compressible modes 2015-06-24 10:44:02 +01:00
Henry Weller
97b0882ef4 timeSelector: Add support for -newTimes option
//- Return the set of times selected based on the argList options
    //  including support for \b -newTimes in which times are selected
    //  if the file <fName> does not exist in the time directory.
    //  Also set the runTime to the first instance or the
    //  \c constant/ directory if no instances are specified or available
    static instantList select
    (
        Time& runTime,
        const argList& args,
        const word& fName
    );

This is experimental functionality and currently on test in the yPlus
post-processing utility.
2015-06-24 10:04:29 +01:00
Henry Weller
d9db912bdd Update headers 2015-06-23 11:27:05 +01:00
Henry Weller
c01b22e038 Removed support for the legacy SiCortex64 architecture 2015-06-23 11:26:45 +01:00
Henry Weller
9ecf5abd53 simpleFoam: Construct a turbulenceModel rather than RASModel
to support laminar simulations in the new framework.

Unfortunately this allows LES models to be instantiated although they
are not compatible with steady-state simulations.
2015-06-23 10:52:45 +01:00
Henry Weller
52025b8776 rhoCentralFoam: Added experimental LTS support
Select LTS via the ddtScheme:

ddtSchemes
{
    default         localEuler rDeltaT;
}

The LTS algorithm is controlled with the standard settings in
controlDict, e.g.:

maxCo           0.5;
maxDeltaT       2e-8;

with the addition of the optional rDeltaT smoothing coefficient:

rDeltaTSmoothingCoeff 0.02;

which defaults to 0.02.

For cases with reasonably uniform meshes like the forwardStep tutorial
LTS does not provide much benefit but for cases with large variation in
cell-size like the biconic25-55Run35 tutorial LTS provides significant
speed-up to convergence particularly if started from uniform conditions.
2015-06-19 11:52:48 +01:00
Henry Weller
b0e7c33736 rhoCentralFoam: Reverted correction to BCs
Reverted changes proposed in
http://openfoam.org/mantisbt/view.php?id=1548 as it adversely affects
fixed-value BCs and is formulated to fix an issue with an unphysical
case.  Further analysis of the handling of fixed pressure outlet
conditions as the Mach number approaches 1 is required.
2015-06-18 15:10:48 +01:00
Henry Weller
d28a72f14c rhoCentralFoam: Updated Courant number calculation to be cell-based rather than face-based
Now consistent with the way the Courant number is calculated for other solvers
2015-06-17 20:20:49 +01:00
Henry Weller
0d3bfc0a96 rhoCentralFoam: Ensure fixed value boundary conditions are preserved
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1748
2015-06-17 08:46:46 +01:00
Henry Weller
04824fb8fc createTurbulenceFields: Construct fields and turbulence model for every time step processed
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1749
2015-06-16 12:49:02 +01:00
Henry Weller
6c66ef9c88 Rationalized the LTS solvers 2015-06-16 12:42:37 +01:00
Henry Weller
c0cb1d6e8f LTSCoalChemistryFoam: LTS version of coalChemistryFoam 2015-06-15 22:36:27 +01:00
Henry Weller
3ed90ae73d reactingTwoPhaseEulerFoam: New twoPhaseEulerFoam supporting mass-transfer and reactions
Multi-species, mass-transfer and reaction support and multi-phase
structure provided by William Bainbridge.

Integration of the latest p-U and face-p_U algorithms with William's
multi-phase structure is not quite complete due to design
incompatibilities which needs further development.  However the
integration of the functionality is complete.

The results of the tutorials are not exactly the same for the
twoPhaseEulerFoam and reactingTwoPhaseEulerFoam solvers but are very
similar.  Further analysis in needed to ensure these differences are
physical or to resolve them; in the meantime the twoPhaseEulerFoam
solver will be maintained.
2015-06-12 09:52:17 +01:00
Henry Weller
b5edc2c9c7 twoPhaseEulerFoam: Minor reorganization 2015-06-12 09:42:19 +01:00
Henry
d617e9b0e7 twoPhaseEulerFoam/interfacialModels/heatTransferModels/sphericalHeatTransfer: new heat-transfer model
Model which applies an analytical solution for heat transfer from the
surface of a sphere to the fluid within the sphere.

Provided by William Bainbridge
2015-06-07 18:57:13 +01:00
Henry
c937c4c9ba twoPhaseEulerFoam: Move the residualAlpha used for drag into the phaseModel
This is necessary to guarantee consistency between the residualAlpha
used for drag and buoyancy in a multi-phase system
2015-06-07 18:55:24 +01:00
Henry
2c61c22312 rhoCentralFoam, sonicFoam: Remove superfluous includes
Resolves bug-report http://openfoam.org/mantisbt/view.php?id=1730
2015-06-03 09:40:25 +01:00
Henry
2b9a2adf8c MRF: Separate MRF from fvOptions
fvOptions does not have the appropriate structure to support MRF as it
is based on option selection by user-specified fields whereas MRF MUST
be applied to all velocity fields in the particular solver.  A
consequence of the particular design choices in fvOptions made it
difficult to support MRF for multiphase and it is easier to support
frame-related and field related options separately.

Currently the MRF functionality provided supports only rotations but
the structure will be generalized to support other frame motions
including linear acceleration, SRF rotation and 6DoF which will be
run-time selectable.
2015-05-29 23:35:43 +01:00
Henry
8b7e02167f MRFZone: rationalize to allow support for general frame acceleration 2015-05-29 15:31:53 +01:00
Henry
1b6f46fae0 simpleFoam/SRFSimpleFoam: Added support for SIMPLEC
SIMPLEC (SIMPLE-consistent) is selected by setting "consistent" option true/yes:

SIMPLE
{
    nNonOrthogonalCorrectors 0;
    consistent yes;
}

which relaxes the pressure in a "consistent" manner and additional
relaxation of the pressure is not generally necessary.  In addition
convergence of the p-U system is better and reliable with less
aggressive relaxation of the momentum equation, e.g. for the motorbike
tutorial:

relaxationFactors
{
    equations
    {
        U               0.9;
        k               0.7;
        omega           0.7;
    }
}

The cost per iteration is marginally higher but the convergence rate is
better so the number of iterations can be reduced.

The SIMPLEC algorithm also provides benefit for cases with large
body-forces, e.g. SRF, see tutorials/incompressible/SRFSimpleFoam/mixer
and feature request http://www.openfoam.org/mantisbt/view.php?id=1714
2015-05-29 11:30:40 +01:00
Henry
e8c7aac911 DPMFoam, MPPICFoam: Corrected contiuousPhaseName -> continuousPhaseName
Resolves second part of http://www.openfoam.org/mantisbt/view.php?id=1717
2015-05-28 17:01:18 +01:00
Henry
59134ff9f8 viewFactorsGen: Correct container resizing
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1712
2015-05-28 16:13:29 +01:00
Henry
71c6af8fe9 mapFields: Reinstated mapFields from OpenFOAM-2.2.x and renamed the current mapFields -> mapFieldsPar
This required the addition of the meshToMesh class in the sampling
library from OpenFOAM-2.2.x which is now named meshToMesh0.
2015-05-26 11:32:46 +01:00
Henry
f60808b4b8 streamFunction: Evaluate which coordinate plan the 2D geometry is in filter-out the normal component
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1703
2015-05-22 22:38:39 +01:00
Henry
658af8b5b2 driftFluxFoam/relativeVelocityModels/Make/options: add newline 2015-05-19 14:34:37 +01:00
Henry
11a0aa03ba extrudeMesh: Add wedge patches after edge collapse 2015-05-18 11:54:03 +01:00
Henry
262ad3f005 functionObjects: Upgrade to use new turbulence modeling library 2015-05-17 14:56:40 +01:00
Henry
bfb726aba3 twoPhaseEulerFoam: rationalize handling of fixed-flux BC update 2015-05-09 23:45:15 +01:00
Henry
227f3a33f5 foamListTimes: Added '-rm' option to remove the selected time directories
If -rm is specified with the -processor option the selected time
directories are removed from all the processor directories.
2015-05-08 22:47:31 +01:00
Henry
5c2fa016e4 Update header 2015-05-08 10:10:47 +01:00
Henry
cd245e9a3f twoPhaseEulerFoam: Update only the fixed-value phi patch fields before constructing the pressure eqn
Avoids small continuity error in parallel
2015-05-08 09:51:36 +01:00
Henry
57b25f22b3 planeExtrusion.[HC]: Renamed from plane.[HC] to avoid name clash 2015-05-04 21:16:44 +01:00