#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 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);const dimensionedScalar &rho1=twoPhaseProperties.rho1();const dimensionedScalar &rho2=twoPhaseProperties.rho2();dimensionedScalar Dab(twoPhaseProperties.lookup("Dab"));dimensionedScalar alphatab(twoPhaseProperties.lookup("alphatab"));volScalarField rho("rho", alpha1 *rho1+(scalar(1)-alpha1)*rho2);rho.oldTime();surfaceScalarField rhoPhi(IOobject("rho*phi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho1 *phi);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<< "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); const dimensionedScalar& rho1 = twoPhaseProperties.rho1(); const dimensionedScalar& rho2 = twoPhaseProperties.rho2(); dimensionedScalar Dab(twoPhaseProperties.lookup("Dab")); dimensionedScalar alphatab(twoPhaseProperties.lookup("alphatab")); volScalarField rho("rho", alpha1*rho1 + (scalar(1) - alpha1)*rho2); rho.oldTime(); surfaceScalarField rhoPhi ( IOobject ( "rho*phi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE ), rho1*phi ); autoPtr<incompressible::turbulenceModel turbulence)(incompressible::turbulenceModel::New(U, phi, twoPhaseProperties)) |
Definition at line 79 of file createFields.H.