#include "src/finiteVolume/cfdTools/incompressible/createPhi.H"
#include "src/finiteVolume/cfdTools/general/include/readGravitationalAcceleration.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();volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT), alpha1 *rho1+(scalar(1)-alpha1)*rho2, alpha1.boundaryField().types());rho.oldTime();surfaceScalarField rhoPhi(IOobject("rho*phi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho1 *phi);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<< "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(); volScalarField rho ( IOobject ( "rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT ), alpha1*rho1 + (scalar(1) - alpha1)*rho2, alpha1.boundaryField().types() ); rho.oldTime(); surfaceScalarField rhoPhi ( IOobject ( "rho*phi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE ), rho1*phi ); interfaceProperties interface(alpha1, U, twoPhaseProperties); autoPtr<incompressible::turbulenceModel turbulence)(incompressible::turbulenceModel::New(U, phi, twoPhaseProperties)) |
Definition at line 92 of file createFields.H.