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"
45 Info<< "Creating phaseChangeTwoPhaseMixture\n" << endl;
46 autoPtr<phaseChangeTwoPhaseMixture> twoPhaseProperties =
47 phaseChangeTwoPhaseMixture::New(U, phi);
49 const dimensionedScalar& rho1 = twoPhaseProperties->rho1();
50 const dimensionedScalar& rho2 = twoPhaseProperties->rho2();
51 const dimensionedScalar& pSat = twoPhaseProperties->pSat();
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 scalar pRefValue = 0.0;
71 setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);
74 // Construct interface from alpha1 distribution
75 interfaceProperties interface(alpha1, U, twoPhaseProperties());
77 // Construct incompressible turbulence model
78 autoPtr<incompressible::turbulenceModel> turbulence
80 incompressible::turbulenceModel::New(U, phi, twoPhaseProperties())