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 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)) |
| 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)) |
Definition at line 188 of file createFields.H.
1.8.17 