1 Info<< "Reading field p\n" << endl;
15 Info<< "Reading field alpha1\n" << endl;
29 Info<< "Reading field U\n" << endl;
43 # include "createPhi.H"
46 Info<< "Reading transportProperties\n" << endl;
47 twoPhaseMixture twoPhaseProperties(U, phi);
49 const dimensionedScalar& rho1 = twoPhaseProperties.rho1();
50 const dimensionedScalar& rho2 = twoPhaseProperties.rho2();
53 // Need to store rho for ddt(rho, U)
61 IOobject::READ_IF_PRESENT
63 alpha1*rho1 + (scalar(1) - alpha1)*rho2,
64 alpha1.boundaryField().types()
70 // Initialisation does not matter because rhoPhi is reset after the
71 // alpha1 solution before it is used in the U equation.
72 surfaceScalarField rhoPhi
86 // Construct interface from alpha1 distribution
87 interfaceProperties interface(alpha1, U, twoPhaseProperties);
89 // Construct incompressible turbulence model
90 autoPtr<incompressible::turbulenceModel> turbulence
92 incompressible::turbulenceModel::New(U, phi, twoPhaseProperties)
97 p.boundaryField().size(),
98 zeroGradientFvPatchScalarField::typeName
101 for (label i=0; i<p.boundaryField().size(); i++)
103 if (p.boundaryField()[i].fixesValue())
105 pcorrTypes[i] = fixedValueFvPatchScalarField::typeName;
111 scalar pRefValue = 0.0;
112 setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);