63int main(
int argc,
char *argv[])
67 "Solver for two incompressible, isothermal immiscible fluids using"
68 " VOF phase-fraction based interface capturing\n"
69 "With optional mesh motion and mesh topology changes including"
70 " adaptive re-meshing."
78 #include "initContinuityErrs.H"
81 #include "createDyMControls.H"
82 #include "createFields.H"
93 IOobject::READ_IF_PRESENT,
102 #include "correctPhi.H"
113 #include "CourantNo.H"
114 #include "setInitialDeltaT.H"
118 Info<<
"\nStarting time loop\n" <<
endl;
122 #include "readControls.H"
126 #include "setRDeltaT.H"
130 #include "CourantNo.H"
131 #include "alphaCourantNo.H"
132 #include "setDeltaT.H"
144 scalar timeBeforeMeshUpdate =
runTime.elapsedCpuTime();
150 Info<<
"Execution time for mesh.update() = "
151 <<
runTime.elapsedCpuTime() - timeBeforeMeshUpdate
156 if (
mesh.topoChanging())
178 Uf += (1-faceMaskOld)*Uint;
183 Uf.boundaryFieldRef()[patchI] =
184 Uint.boundaryField()[patchI];
192 #include "correctPhi.H"
206 fvc::makeRelative(
phi,
U);
217 #include "alphaControls.H"
218 #include "alphaEqnSubCycle.H"
244 runTime.printExecutionTime(Info);
CMULES: Multidimensional universal limiter for explicit corrected implicit solution.
const uniformDimensionedVectorField & g
const surfaceScalarField & ghf
const volScalarField & gh
surfaceScalarField faceMask(localMin< scalar >(mesh).interpolate(cellMask))
tmp< surfaceScalarField > talphaPhi1Corr0
Read the control parameters used by setDeltaT.
autoPtr< surfaceVectorField > Uf
Creates and initialises the velocity velocity field Uf.
compressible::turbulenceModel & turbulence
Provides functions smooth spread and sweep which use the FaceCellWave algorithm to smooth and redistr...
tmp< volScalarField > rAU
Calculates and outputs the mean and maximum Courant Numbers.
GeometricField< vector, fvsPatchField, surfaceMesh > surfaceVectorField
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
messageStream Info
Information stream (stdout output on master, null elsewhere)
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Ostream & endl(Ostream &os)
Add newline and flush stream.
constexpr char nl
The newline '\n' character (0x0a)
Adjust the balance of fluxes on the faces between interpolated and calculated to obey continuity.
Execute application functionObjects to post-process existing results.
Sets blocked cells mask field.
Sets blocked cells mask field.
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
#define forAll(list, i)
Loop across all elements in list.