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<< "Constructing twoPhaseMixtureThermo\"<< endl;twoPhaseMixtureThermo mixture(U, phi);volScalarField &alpha1(mixture.alpha1());volScalarField &alpha2(mixture.alpha2());Info<< "Reading thermophysical properties\"<< endl;const volScalarField &rho1=mixture.thermo1().rho();const volScalarField &rho2=mixture.thermo2().rho();volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), alpha1 *rho1+alpha2 *rho2);dimensionedScalar pMin("pMin", dimPressure, mixture);mesh.setFluxRequired(p_rgh.name());mesh.setFluxRequired(alpha1.name());surfaceScalarField rhoPhi(IOobject("rhoPhi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), fvc::interpolate(rho) *phi);volScalarField dgdt(alpha1 *fvc::div(phi));compressibleInterPhaseTransportModel turbulence(rho, U, phi, rhoPhi, alphaPhi10, mixture);{ wordHashSet &nonInt=const_cast< wordHashSet & > |
( |
Stencil::New(mesh).nonInterpolatedFields() |
| ) |
|