Definition in file createFields.H.
#include "src/finiteVolume/cfdTools/incompressible/createPhi.H"
Include dependency graph for createFields.H:Go to the source code of this file.
Variables | |
| Info<< "Reading field p_rgh\n" << endl;volScalarField p_rgh(IOobject("p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info << "Reading field alpha1\n" << endl;volScalarField alpha1(IOobject("alpha1", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info << "Calculating field alpha1\n" << endl;volScalarField alpha2("alpha2", scalar(1)-alpha1);Info << "Reading field U\n"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info << "Reading transportProperties\n"<< endl;twoPhaseMixture twoPhaseProperties(U, phi);dimensionedScalar rho10(twoPhaseProperties.subDict(twoPhaseProperties.phase1Name()).lookup("rho0"));dimensionedScalar rho20(twoPhaseProperties.subDict(twoPhaseProperties.phase2Name()).lookup("rho0"));dimensionedScalar psi1(twoPhaseProperties.subDict(twoPhaseProperties.phase1Name()).lookup("psi"));dimensionedScalar psi2(twoPhaseProperties.subDict(twoPhaseProperties.phase2Name()).lookup("psi"));dimensionedScalar pMin(twoPhaseProperties.lookup("pMin"));Info << "Calculating field g.h\n" << endl;volScalarField gh("gh", g &mesh.C());surfaceScalarField ghf("ghf", g &mesh.Cf());volScalarField p(IOobject("p", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), max((p_rgh+gh *(alpha1 *rho10+alpha2 *rho20))/(1.0-gh *(alpha1 *psi1+alpha2 *psi2)), pMin));volScalarField rho1=rho10+psi1 *p;volScalarField rho2=rho20+psi2 *p;volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), alpha1 *rho1+alpha2 *rho2);surfaceScalarField rhoPhi(IOobject("rho*phi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), fvc::interpolate(rho)*phi);volScalarField dgdt=pos(alpha2)*fvc::div(phi)/max(alpha2, scalar(0.0001));interfaceProperties interface(alpha1, U, twoPhaseProperties);autoPtr < incompressible::turbulenceModel | turbulence )(incompressible::turbulenceModel::New(U, phi, twoPhaseProperties)) |
| Info<< "Reading field p_rgh\n" << endl; volScalarField p_rgh ( IOobject ( "p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); Info<< "Reading field alpha1\n" << endl; volScalarField alpha1 ( IOobject ( "alpha1", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); Info<< "Calculating field alpha1\n" << endl; volScalarField alpha2("alpha2", scalar(1) - alpha1); Info<< "Reading field U\n" << endl; volVectorField U ( IOobject ( "U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); Info<< "Reading transportProperties\n" << endl; twoPhaseMixture twoPhaseProperties(U, phi); dimensionedScalar rho10 ( twoPhaseProperties.subDict ( twoPhaseProperties.phase1Name() ).lookup("rho0") ); dimensionedScalar rho20 ( twoPhaseProperties.subDict ( twoPhaseProperties.phase2Name() ).lookup("rho0") ); dimensionedScalar psi1 ( twoPhaseProperties.subDict ( twoPhaseProperties.phase1Name() ).lookup("psi") ); dimensionedScalar psi2 ( twoPhaseProperties.subDict ( twoPhaseProperties.phase2Name() ).lookup("psi") ); dimensionedScalar pMin(twoPhaseProperties.lookup("pMin")); Info<< "Calculating field g.h\n" << endl; volScalarField gh("gh", g & mesh.C()); surfaceScalarField ghf("ghf", g & mesh.Cf()); volScalarField p ( IOobject ( "p", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), max ( (p_rgh + gh*(alpha1*rho10 + alpha2*rho20)) /(1.0 - gh*(alpha1*psi1 + alpha2*psi2)), pMin ) ); volScalarField rho1 = rho10 + psi1*p; volScalarField rho2 = rho20 + psi2*p; volScalarField rho ( IOobject ( "rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE ), alpha1*rho1 + alpha2*rho2 ); surfaceScalarField rhoPhi ( IOobject ( "rho*phi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE ), fvc::interpolate(rho)*phi ); volScalarField dgdt = pos(alpha2)*fvc::div(phi)/max(alpha2, scalar(0.0001)); interfaceProperties interface(alpha1, U, twoPhaseProperties); autoPtr<incompressible::turbulenceModel turbulence)(incompressible::turbulenceModel::New(U, phi, twoPhaseProperties)) |
Definition at line 149 of file createFields.H.
