6 Info<<
"Reading field p_rgh\n" << endl;
21 Info<<
"Reading field U\n" << endl;
35 #include "createPhi.H" 37 Info<<
"Constructing twoPhaseMixtureThermo\n" << endl;
38 twoPhaseMixtureThermo mixture(mesh);
40 volScalarField& alpha1(mixture.alpha1());
41 volScalarField& alpha2(mixture.alpha2());
43 Info<<
"Reading thermophysical properties\n" << endl;
45 volScalarField& p = mixture.p();
46 volScalarField& T = mixture.T();
48 Info<<
"Reading field rho_foam \n" << endl;
49 volScalarField rho_foam
62 Info<<
"Reading field rho_gas \n" << endl;
63 volScalarField rho_gas
82 IOobject::READ_IF_PRESENT,
85 alpha1*rho_gas + (1-alpha1)*rho_foam
94 IOobject::READ_IF_PRESENT,
98 dimensionedScalar(
"Psi1", dimensionSet(0,-2,2,0,0,0,0), 1e-5)
108 IOobject::READ_IF_PRESENT,
112 dimensionedScalar(
"Psi2", dimensionSet(0,-2,2,0,0,0,0), 1e-2)
115 volScalarField Psi_bubble
122 IOobject::READ_IF_PRESENT,
126 dimensionedScalar(
"Psi_bubble", dimensionSet(0,-2,2,0,0,0,0), 1e-2)
129 volScalarField Psi_foam
136 IOobject::READ_IF_PRESENT,
140 dimensionedScalar(
"Psi_foam", dimensionSet(0,-2,2,0,0,0,0), 1e-2)
143 dimensionedScalar pMin(mixture.lookup(
"pMin"));
145 Info<<
"Calculating field g.h\n" << endl;
146 volScalarField gh(
"gh", g & mesh.C());
147 surfaceScalarField ghf(
"ghf", g & mesh.Cf());
152 surfaceScalarField rhoPhi
162 fvc::interpolate(rho)*phi
167 pos(alpha2)*fvc::div(phi)/max(alpha2, scalar(0.0001))
172 interfaceProperties interface(alpha1, U, mixture);
179 autoPtr<compressible::turbulenceModel> turbulence
181 compressible::turbulenceModel::New(rho, U, rhoPhi, mixture)
184 Info<<
"Creating field kinetic energy K\n" << endl;
185 volScalarField K(
"K", 0.5*magSqr(U));
187 dimensionedScalar densityDimension (
"densityDimension", dimensionSet(1,-3,0,0,0,0,0), 1.0);
188 dimensionedScalar initialFoamMass =
190 fvc::domainIntegrate(rho_foam*alpha2)
191 + densityDimension*
L0*
rhoPoly*fvc::domainIntegrate(alpha2)
193 dimensionedScalar cumulativeContinuityError = 0.0;
creates the field variables related to the PU kinetics
double rhoPoly
density of the liquid polymer, kg/m3
defines mold-filling and validation as two simulation modes
double L0
Initial weight fraction of blowing agent in the liquid, -.
creates the field variables related to the moments
creates the field variables for the temperature equation
creates the field variables for the rheology