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<< "Reading transportProperties\n"<< endl;immiscibleIncompressibleTwoPhaseMixture mixture(U, phi);volScalarField &alpha1(mixture.alpha1());volScalarField &alpha2(mixture.alpha2());const dimensionedScalar &rho1=mixture.rho1();const dimensionedScalar &rho2=mixture.rho2();volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), alpha1 *rho1+alpha2 *rho2);rho.oldTime();volScalarField mu(IOobject("mu", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT), mixture.mu(), calculatedFvPatchScalarField::typeName);surfaceScalarField rhoPhi(IOobject("rhoPhi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), fvc::interpolate(rho) *phi);volScalarField p(IOobject("p", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), p_rgh+rho *gh);label pRefCell=0;scalar pRefValue=0.0;setRefCell(p, p_rgh, mesh.solutionDict().subDict("PIMPLE"), pRefCell, pRefValue);if(p_rgh.needReference()){ p+=dimensionedScalar("p", p.dimensions(), pRefValue - getRefCellValue(p, pRefCell));p_rgh=p - rho *gh;}mesh.setFluxRequired(p_rgh.name());mesh.setFluxRequired(alpha1.name());volScalarField alphac(IOobject("alphac", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), mesh, dimensionedScalar(dimless, Zero), zeroGradientFvPatchScalarField::typeName);alphac.oldTime();volScalarField alphacRho(alphac *rho);alphacRho.oldTime();Info<< "Constructing kinematicCloud "<< endl;basicKinematicCloud kinematicCloud("kinematicCloud", rho, U, mu, g);scalar alphacMin(1.0 -(kinematicCloud.particleProperties().subDict("constantProperties") .get< scalar >("alphaMax")));alphac=max(1.0 - kinematicCloud.theta(), alphacMin);alphac.correctBoundaryConditions();surfaceScalarField alphacf("alphacf", fvc::interpolate(alphac));surfaceScalarField alphaRhoPhic("alphaRhoPhic", alphacf *rhoPhi);surfaceScalarField alphaPhic("alphaPhic", alphacf *phi);autoPtr< PhaseCompressibleTurbulenceModel< immiscibleIncompressibleTwoPhaseMixture > turbulence (PhaseCompressibleTurbulenceModel< immiscibleIncompressibleTwoPhaseMixture >::New(alphac, rho, U, alphaRhoPhic, rhoPhi, mixture))
 

Function Documentation

◆ turbulence()

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;immiscibleIncompressibleTwoPhaseMixture mixture(U, phi);volScalarField &alpha1(mixture.alpha1());volScalarField &alpha2(mixture.alpha2());const dimensionedScalar &rho1=mixture.rho1();const dimensionedScalar &rho2=mixture.rho2();volScalarField rho( IOobject( "rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), alpha1 *rho1+alpha2 *rho2);rho.oldTime();volScalarField mu( IOobject( "mu", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT), mixture.mu(), calculatedFvPatchScalarField::typeName);surfaceScalarField rhoPhi( IOobject( "rhoPhi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), fvc::interpolate(rho) *phi);volScalarField p( IOobject( "p", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), p_rgh+rho *gh);label pRefCell=0;scalar pRefValue=0.0;setRefCell( p, p_rgh, mesh.solutionDict().subDict("PIMPLE"), pRefCell, pRefValue);if(p_rgh.needReference()){ p+=dimensionedScalar( "p", p.dimensions(), pRefValue - getRefCellValue(p, pRefCell)); p_rgh=p - rho *gh;}mesh.setFluxRequired(p_rgh.name());mesh.setFluxRequired(alpha1.name());volScalarField alphac( IOobject( "alphac", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), mesh, dimensionedScalar(dimless, Zero), zeroGradientFvPatchScalarField::typeName);alphac.oldTime();volScalarField alphacRho(alphac *rho);alphacRho.oldTime();Info<< "Constructing kinematicCloud "<< endl;basicKinematicCloud kinematicCloud( "kinematicCloud", rho, U, mu, g);scalar alphacMin( 1.0 -( kinematicCloud.particleProperties().subDict("constantProperties") .get< scalar >("alphaMax")));alphac=max(1.0 - kinematicCloud.theta(), alphacMin);alphac.correctBoundaryConditions();surfaceScalarField alphacf("alphacf", fvc::interpolate(alphac));surfaceScalarField alphaRhoPhic("alphaRhoPhic", alphacf *rhoPhi);surfaceScalarField alphaPhic("alphaPhic", alphacf *phi);autoPtr< PhaseCompressibleTurbulenceModel< immiscibleIncompressibleTwoPhaseMixture > turbulence ( PhaseCompressibleTurbulenceModel< immiscibleIncompressibleTwoPhaseMixture >  ::Newalphac, rho, U, alphaRhoPhic, rhoPhi, mixture)