1 Info<< "Reading field pd\n" << endl;
15 Info<< "Reading field gamma\n" << endl;
29 Info<< "Reading field U\n" << endl;
43 # include "createPhi.H"
46 Info<< "Reading transportProperties\n" << endl;
47 twoPhaseMixture twoPhaseProperties(U, phi, "gamma");
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 gamma*rho1 + (scalar(1) - gamma)*rho2,
64 gamma.boundaryField().types()
70 // Initialisation does not matter because rhoPhi is reset after the
71 // gamma solution before it is used in the U equation.
72 surfaceScalarField rhoPhi
86 Info<< "Calculating field g.h\n" << endl;
87 volScalarField gh("gh", g & mesh.C());
88 surfaceScalarField ghf("gh", g & mesh.Cf());
106 scalar pdRefValue = 0.0;
107 setRefCell(pd, mesh.solutionDict().subDict("PISO"), pdRefCell, pdRefValue);
110 // Construct interface from gamma distribution
111 interfaceProperties interface(gamma, U, twoPhaseProperties);