compressibleAlphaEqnSubCycle.H
Go to the documentation of this file.
1 tmp<surfaceScalarField> talphaPhi1(alphaPhi10);
2 
4 {
5  dimensionedScalar totalDeltaT = runTime.deltaT();
6 
8  (
9  IOobject
10  (
11  "alphaPhi1",
12  runTime.timeName(),
13  mesh
14  ),
15  mesh,
16  dimensionedScalar("0", alphaPhi10.dimensions(), 0)
17  );
18 
19  surfaceScalarField rhoPhiSum
20  (
21  IOobject
22  (
23  "rhoPhiSum",
24  runTime.timeName(),
25  mesh
26  ),
27  mesh,
28  dimensionedScalar("0", rhoPhi.dimensions(), 0)
29  );
30 
31  tmp<volScalarField> trSubDeltaT;
32 
33  if (LTS)
34  {
35  trSubDeltaT =
36  fv::localEulerDdt::localRSubDeltaT(mesh, nAlphaSubCycles);
37  }
38 
39  for
40  (
41  subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
42  !(++alphaSubCycle).end();
43  )
44  {
45  #include "alphaEqn.H"
46  talphaPhi1.ref() += (runTime.deltaT()/totalDeltaT)*alphaPhi10;
47  rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
48  }
49 
50  rhoPhi = rhoPhiSum;
51 }
52 else
53 {
54  #include "alphaEqn.H"
55 }
56 
58 
59 const surfaceScalarField& alphaPhi1 = talphaPhi1();
60 surfaceScalarField alphaPhi2("alphaPhi2", phi - alphaPhi1);
61 
62 volScalarField::Internal contErr
63 (
64  (
66  - (fvOptions(alpha1, mixture.thermo1().rho())&rho1)
67  - (fvOptions(alpha2, mixture.thermo2().rho())&rho2)
68  )()
69 );
fv::options & fvOptions
surfaceScalarField & phi
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
engineTime & runTime
surfaceScalarField alphaPhi10(alphaPhi10Header, phi *fvc::interpolate(alpha1))
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45
volScalarField::Internal contErr((fvc::ddt(rho)+fvc::div(rhoPhi) -(fvOptions(alpha1, mixture.thermo1().rho())&rho1) -(fvOptions(alpha2, mixture.thermo2().rho())&rho2))())
dynamicFvMesh & mesh
const volScalarField & alpha1
alpha2
Definition: alphaEqn.H:115
bool LTS
Definition: createRDeltaT.H:1
Info<< "Reading field p_rgh\"<< endl;volScalarField p_rgh(IOobject("p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Creating phaseChangeTwoPhaseMixture\"<< endl;autoPtr< phaseChangeTwoPhaseMixture > mixture
Definition: createFields.H:33
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
rhoPhi
Definition: rhoEqn.H:10
tmp< surfaceScalarField > talphaPhi1(alphaPhi10)
surfaceScalarField alphaPhi2("alphaPhi2", phi - alphaPhi1)
const dimensionedScalar & rho2
Definition: createFields.H:40
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
const dimensionedScalar & rho1
Definition: createFields.H:39
label nAlphaSubCycles(readLabel(alphaControls.lookup("nAlphaSubCycles")))