Commit Graph

16006 Commits

Author SHA1 Message Date
Henry Weller
91ee15e3ac 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
80e13f8002 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
d1c18321a7 LimitedScheme/LimitFuncs: Make the limiter function for symmTensor consistent with tensor
Use the trace as the limiter function for both symmTensor and tensor to
bound the normal stresses rather than the shear stresses.

Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1751
2015-06-18 11:45:40 +01:00
Henry Weller
164448f8b1 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
bc493180d9 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
280f0fbf31 PatchInjection::parcelsToInject: Ensure parallel consistency
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1746
2015-06-16 15:14:15 +01:00
Henry Weller
35c53a0ee8 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
b1a1a0c60d filmViscosityModel/thixotropicViscosity: Corrected the handling of the effect of droplet deposition
on the lambda of the film
2015-06-16 12:43:40 +01:00
Henry Weller
8293227964 Rationalized the LTS solvers 2015-06-16 12:42:37 +01:00
Henry Weller
03adb37ab4 LTSCoalChemistryFoam: LTS version of coalChemistryFoam 2015-06-15 22:36:27 +01:00
Henry Weller
fc8d1c4c14 externalWallHeatFluxTemperature: Correct doc
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1745
2015-06-15 19:13:31 +01:00
Henry Weller
08191ee828 ReactingMultiphaseParcel,ReactingParcel: Use T0 for consistent radiation transfer
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1636
2015-06-15 17:40:48 +01:00
Henry Weller
f9358d33a0 pairPatchAgglomeration: Stop agglomeration loop if no further agglomeration is achieved
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1157
2015-06-15 17:14:24 +01:00
Henry Weller
57e5a11041 Corrected typo
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1741
2015-06-15 11:20:56 +01:00
Henry Weller
406f0ab830 pairPatchAgglomeration: Do not use the result agglomeratePatch if the return indicates it is not valid
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1157
2015-06-13 23:05:45 +01:00
Henry Weller
e9f35a9d19 Lagrangian radiation models and fvDOM: Apply updates
Applying patches provided by Timo Niemi
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1636
2015-06-12 15:35:43 +01:00
Henry Weller
eb53f9bdf0 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
ee6660794f Updated header 2015-06-12 09:44:48 +01:00
Henry Weller
573021c9ba TurbulenceModels/phaseCompressible: Do not call correctNut in constructor
if it depends on the availability of the physical properties of both
phases as the construction of the phase-system may not be complete.
2015-06-12 09:43:07 +01:00
Henry Weller
360604b104 twoPhaseEulerFoam: Minor reorganization 2015-06-12 09:42:19 +01:00
Henry
d71b839320 polyMesh::findCell: Call tetBasePtIs for CELL_TETS to ensure parallel consistency
Resolves additional bug-report http://www.openfoam.org/mantisbt/view.php?id=1544
2015-06-11 17:45:25 +01:00
Henry
a5b1708e63 tetDecomposer: Rename enumeration elements for consistency with polyMesh 2015-06-11 16:08:41 +01:00
Henry
7366e04ac4 calculateMeshToMesh0Addressing: Change decomposition method to CELL_TETS
Resolves additional request in http://www.openfoam.org/mantisbt/view.php?id=1544
2015-06-11 15:29:34 +01:00
Henry
8eda4b8632 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
61e52b2cb4 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
3f62ba192f Update headers 2015-06-07 18:54:19 +01:00
Henry
661f5aa50c TurbulenceModels/phaseCompressible: Up-cast the phase-system to twoPhaseSystem 2015-06-07 18:53:22 +01:00
Henry
57d75a691f combustionModel: Change the base turbulenceModel to compressibleTurbulenceModel
to allow combustionModel to be used with more general forms of multi-phase
2015-06-07 18:51:36 +01:00
Henry
320b8708c2 combustionModels: Added multiphase-phase support
by naming fields including the phase name is needed
2015-06-05 08:19:13 +01:00
Henry
28f58cb679 IOobject: Added member function "member()" which returns the member name of the object
i.e. without the group suffix.

Patch provided by William Bainbridge
2015-06-04 17:09:42 +01:00
Henry
76f7fd18f9 filmViscosityModel/ArrheniusViscosity: temperature-dependent viscosity function
Applied as a correction pre-factor to the viscosity obtained from a
run-time selectable base-model.
2015-06-04 11:34:21 +01:00
Henry
f5ce55c181 filmViscosityModel/thixotropicViscosity/thixotropicViscosity: Remove the spurious mask and correct the blending 2015-06-03 17:29:18 +01:00
Henry
af51444760 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
ca82a1078a solidificationMeltingSource: Removed all need for the "alwaysApply" hack
Also removed the underlying structure for "alwaysApply" from fvOption
2015-06-01 21:53:13 +01:00
Henry
35e2d24dce Updated header 2015-06-01 15:47:50 +01:00
Henry
7608bb5bc6 fixedTemperatureConstraint,effectivenessHeatExchangerSource: Set the field name to that of the energy
field from which the temperature is obtained

Removes the need for the 'alwaysApply' hack.
2015-06-01 15:46:17 +01:00
Henry
25944abbc8 limitTemperature: Set the field name to that of the energy field from which the temperature is obtained
Removes the need for the 'alwaysApply' hack.
2015-06-01 14:21:44 +01:00
Henry
257d71ebf1 fvOption: Rationalize corrections 2015-06-01 11:18:45 +01:00
Henry
bbbad9cfc0 fvOptions: Rationalize naming convention for constraints 2015-06-01 09:58:33 +01:00
Henry
6ffc90ad29 Update headers 2015-05-31 17:15:44 +01:00
Henry
b4a94650c9 fvOption: Rationalize read to support re-reading
writeData is not yet completely consistent with read.
2015-05-31 17:14:32 +01:00
Henry
3a004fda10 fvOptions: Separate options for all cells, cellSets and inter-region coupling
by introducing rational base-classes rather than using the hideous
'switch' statement.  Further rationalization of the cell-selection
mechanism will be implemented via an appropriate class hierarchy to
replace the remaining 'switch' statement.

Mesh-motion is currently handled very inefficiently for cellSets and not
at all for inter-region coupling.  The former will be improved when the
cell-selection classes are written and the latter by making the
meshToMesh class a MeshObject after it has been corrected for mapFields.
2015-05-31 16:38:01 +01:00
Henry
c3ee2348a6 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
ec6ad11cbd tutorials/incompressible/simpleFoam/pitzDaily/0/omega: Correct value to be consistent with epsilon
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1720
2015-05-29 21:26:49 +01:00
Henry
38171e0c7a MRFZone: rationalize to allow support for general frame acceleration 2015-05-29 15:31:53 +01:00
Henry
38177526ff 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
1971c1eb88 fvPatchFields: Support explicitly named U, phi and rho in all BCs
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1717
2015-05-28 20:42:47 +01:00
Henry
0d69b108a2 tutorials/lagrangian/MPPICFoam: Remove spurious phi specifications in p BCs 2015-05-28 17:09:55 +01:00
Henry
1568641986 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
35fd6f4894 cyclicAMIPolyPatch: Correct report of sweep angle
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1716
2015-05-28 16:29:06 +01:00