Go to the source code of this file.
|
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 U\n"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading transportProperties\n"<< endl;incompressibleTwoPhaseMixture mixture(U, phi);volScalarField &alpha1(mixture.alpha1());volScalarField &alpha2(mixture.alpha2());const dimensionedScalar &rho1=mixture.rho1();const dimensionedScalar &rho2=mixture.rho2();dimensionedScalar Dab("Dab", dimViscosity, mixture);dimensionedScalar alphatab("alphatab", dimless, mixture);volScalarField rho("rho", alpha1 *rho1+alpha2 *rho2);rho.oldTime();surfaceScalarField rhoPhi(IOobject("rhoPhi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho1 *phi);autoPtr< incompressible::turbulenceModel > | turbulence (incompressible::turbulenceModel::New(U, phi, mixture)) |
|
label | pRefCell = 0 |
|
scalar | pRefValue = 0.0 |
|
volScalarField p |
( |
IOobject("p", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE) |
, |
|
|
p_rgh+rho * |
gh |
|
) |
| |
if |
( |
p_rgh. |
needReference() | ) |
|
mesh setFluxRequired |
( |
p_rgh. |
name() | ) |
|
mesh setFluxRequired |
( |
alpha1. |
name() | ) |
|
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 U\n" << endl;volVectorField U( IOobject ( "U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh);Info<< "Reading transportProperties\n" << endl;incompressibleTwoPhaseMixture mixture(U, phi);volScalarField& alpha1(mixture.alpha1());volScalarField& alpha2(mixture.alpha2());const dimensionedScalar& rho1 = mixture.rho1();const dimensionedScalar& rho2 = mixture.rho2();dimensionedScalar Dab("Dab", dimViscosity, mixture);dimensionedScalar alphatab("alphatab", dimless, mixture);volScalarField rho("rho", alpha1*rho1 + alpha2*rho2);rho.oldTime();surfaceScalarField rhoPhi( IOobject ( "rhoPhi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE ), rho1*phi);autoPtr<incompressible::turbulenceModel> turbulence(incompressible::turbulenceModel::New(U, phi, mixture)) |