60 int main(
int argc,
char *argv[])
64 "Solver for two incompressible, isothermal immiscible fluids using" 65 " VOF phase-fraction based interface capturing\n" 66 "With optional mesh motion and mesh topology changes including" 67 " adaptive re-meshing." 75 #include "initContinuityErrs.H" 78 #include "createFields.H" 89 IOobject::READ_IF_PRESENT,
97 #include "createControls.H" 106 #include "CourantNo.H" 107 #include "setInitialDeltaT.H" 111 Info<<
"\nStarting time loop\n" <<
endl;
120 #include "setRDeltaT.H" 124 #include "CourantNo.H" 125 #include "alphaCourantNo.H" 126 #include "setDeltaT.H" 138 scalar timeBeforeMeshUpdate =
runTime.elapsedCpuTime();
144 Info<<
"Execution time for mesh.update() = " 145 <<
runTime.elapsedCpuTime() - timeBeforeMeshUpdate
150 if (
mesh.topoChanging())
182 #include "alphaControls.H" 183 #include "alphaEqnSubCycle.H"
Sets blocked cells mask field.
constexpr char nl
The newline '\n' character (0x0a)
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Ostream & endl(Ostream &os)
Add newline and flush stream.
tmp< volScalarField > rAU
Sets blocked cells mask field.
const surfaceScalarField & ghf
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Creates and initialises the velocity velocity field Uf.
bool oversetAdjustPhi(surfaceScalarField &phi, const volVectorField &U, const label zoneId=-1)
Adjust the balance of fluxes to obey continuity.
const uniformDimensionedVectorField & g
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar("pos", dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar("neg", dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
Correct phi on new faces C-I faces.
Calculates and outputs the mean and maximum Courant Numbers.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const volScalarField & gh
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
messageStream Info
Information stream (stdout output on master, null elsewhere)
Execute application functionObjects to post-process existing results.
void makeRelative(surfaceScalarField &phi, const volVectorField &U)
Make the given flux relative.
Provides functions smooth spread and sweep which use the FaceCellWave algorithm to smooth and redistr...
CMULES: Multidimensional universal limiter for explicit corrected implicit solution.
Adjust the balance of fluxes on the faces between interpolated and calculated to obey continuity...
tmp< surfaceScalarField > talphaPhi1Corr0