59int main(
int argc,
char *argv[])
63 "Solver for two incompressible, isothermal immiscible fluids using"
64 " VOF phase-fraction based interface capturing.\n"
65 "Includes MRF and an MPPIC cloud."
74 #include "createControl.H"
76 #include "initContinuityErrs.H"
77 #include "createFields.H"
80 #include "correctPhi.H"
87 #include "CourantNo.H"
88 #include "setInitialDeltaT.H"
93 Info<<
"\nStarting time loop\n" <<
endl;
101 #include "setRDeltaT.H"
105 #include "CourantNo.H"
106 #include "alphaCourantNo.H"
107 #include "setDeltaT.H"
114 Info<<
"Evolving " << kinematicCloud.name() <<
endl;
116 kinematicCloud.evolve();
119 alphac =
max(1.0 - kinematicCloud.theta(), alphacMin);
120 alphac.correctBoundaryConditions();
122 Info<<
"Continuous phase-1 volume fraction = "
123 << alphac.weightedAverage(
mesh.Vsc()).value()
124 <<
" Min(alphac) = " <<
min(alphac).value()
125 <<
" Max(alphac) = " <<
max(alphac).value()
128 alphacf = fvc::interpolate(alphac);
129 alphaRhoPhic = alphacf*
rhoPhi;
130 alphaPhic = alphacf*
phi;
131 alphacRho = alphac*
rho;
147 cloudVolSUSu.primitiveFieldRef() = -cloudSU.source()/
mesh.V();
148 cloudVolSUSu.correctBoundaryConditions();
150 cloudSU.source() = vector::zero;
155 #include "alphaControls.H"
156 #include "alphaEqnSubCycle.H"
176 runTime.printExecutionTime(Info);
CMULES: Multidimensional universal limiter for explicit corrected implicit solution.
Read the control parameters used by setDeltaT.
compressible::turbulenceModel & turbulence
Provides functions smooth spread and sweep which use the FaceCellWave algorithm to smooth and redistr...
GeometricField< vector, fvPatchField, volMesh > volVectorField
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
fvMatrix< vector > fvVectorMatrix
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
dimensioned< vector > dimensionedVector
Dimensioned vector obtained from generic dimensioned type.
constexpr char nl
The newline '\n' character (0x0a)
Execute application functionObjects to post-process existing results.
Read the control parameters used by setDeltaT.
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
static const char *const typeName
The type name used in ensight case files.