Functions
createFields.H File Reference
Include dependency graph for createFields.H:

Go to the source code of this file.

Functions

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<< "Constructing twoPhaseMixtureThermo\n"<< endl;twoPhaseMixtureThermo mixture(mesh);volScalarField &alpha1(mixture.alpha1());volScalarField &alpha2(mixture.alpha2());Info<< "Reading thermophysical properties\n"<< endl;volScalarField &p=mixture.p();volScalarField &T=mixture.T();volScalarField &rho1=mixture.thermo1().rho();const volScalarField &psi1=mixture.thermo1().psi();volScalarField &rho2=mixture.thermo2().rho();const volScalarField &psi2=mixture.thermo2().psi();volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), alpha1 *rho1+alpha2 *rho2);dimensionedScalar pMin("pMin", dimPressure, mixture);mesh.setFluxRequired(p_rgh.name());surfaceScalarField rhoPhi(IOobject("rhoPhi", 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, mixture);autoPtr< compressible::turbulenceModel > turbulence (compressible::turbulenceModel::New(rho, U, rhoPhi, mixture))
 

Function Documentation

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<< "Constructing twoPhaseMixtureThermo\n" << endl;twoPhaseMixtureThermo mixture(mesh);volScalarField& alpha1(mixture.alpha1());volScalarField& alpha2(mixture.alpha2());Info<< "Reading thermophysical properties\n" << endl;volScalarField& p = mixture.p();volScalarField& T = mixture.T();volScalarField& rho1 = mixture.thermo1().rho();const volScalarField& psi1 = mixture.thermo1().psi();volScalarField& rho2 = mixture.thermo2().rho();const volScalarField& psi2 = mixture.thermo2().psi();volScalarField rho( IOobject ( "rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE ), alpha1*rho1 + alpha2*rho2);dimensionedScalar pMin( "pMin", dimPressure, mixture);mesh.setFluxRequired(p_rgh.name());surfaceScalarField rhoPhi( IOobject ( "rhoPhi", 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, mixture);autoPtr<compressible::turbulenceModel> turbulence ( compressible::turbulenceModel::New(rho, U, rhoPhi, mixture )