#include "src/finiteVolume/cfdTools/incompressible/createPhi.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 << "Reading field alpha2\n" << endl;volScalarField alpha2(IOobject("alpha2", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info << "Reading field alpha3\n" << endl;volScalarField alpha3(IOobject("alpha3", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);alpha3==1.0-alpha1-alpha2;Info << "Reading field U\n"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);threePhaseMixture threePhaseProperties(U, phi);const dimensionedScalar &rho1=threePhaseProperties.rho1();const dimensionedScalar &rho2=threePhaseProperties.rho2();const dimensionedScalar &rho3=threePhaseProperties.rho3();dimensionedScalar D23(threePhaseProperties.lookup("D23"));volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT), alpha1 *rho1+alpha2 *rho2+alpha3 *rho3, alpha1.boundaryField().types());rho.oldTime();surfaceScalarField rhoPhi(IOobject("rho*phi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho1 *phi);threePhaseInterfaceProperties interface(threePhaseProperties);autoPtr < incompressible::turbulenceModel | turbulence )(incompressible::turbulenceModel::New(U, phi, threePhaseProperties)) |
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<< "Reading field alpha2\n" << endl; volScalarField alpha2 ( IOobject ( "alpha2", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); Info<< "Reading field alpha3\n" << endl; volScalarField alpha3 ( IOobject ( "alpha3", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); alpha3 == 1.0 - alpha1 - alpha2; Info<< "Reading field U\n" << endl; volVectorField U ( IOobject ( "U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); threePhaseMixture threePhaseProperties(U, phi); const dimensionedScalar& rho1 = threePhaseProperties.rho1(); const dimensionedScalar& rho2 = threePhaseProperties.rho2(); const dimensionedScalar& rho3 = threePhaseProperties.rho3(); dimensionedScalar D23(threePhaseProperties.lookup("D23")); volScalarField rho ( IOobject ( "rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT ), alpha1*rho1 + alpha2*rho2 + alpha3*rho3, alpha1.boundaryField().types() ); rho.oldTime(); surfaceScalarField rhoPhi ( IOobject ( "rho*phi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE ), rho1*phi ); threePhaseInterfaceProperties interface(threePhaseProperties); autoPtr<incompressible::turbulenceModel turbulence)(incompressible::turbulenceModel::New(U, phi, threePhaseProperties)) |
Definition at line 125 of file createFields.H.