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

Go to the source code of this file.

Functions

Info<< "Reading thermophysical properties\"<< endl;autoPtr< psiuReactionThermo > pThermo (psiuReactionThermo::New(mesh))
 
thermo validate (args.executable(), "ha", "ea")
 
volScalarField rho (IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho())
 
Info<< "min(b) = "<< min(b).value()<< endl;Info<< "\Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);mesh.schemes().setFluxRequired(p.name());Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::RASModel > turbulence (compressible::New< compressible::RASModel >(rho, U, phi, thermo))
 
Info<< "Creating thermophysical transport model\"<< endl;autoPtr< RASFluidThermophysicalTransportModel > thermophysicalTransport (RASFluidThermophysicalTransportModel::New(turbulence(), thermo))
 
 if (composition.contains("ft"))
 
fields add (b)
 
fields add (thermo.he())
 
flameWrinkling addXi (fields)
 

Variables

psiuReactionThermo & thermo = pThermo()
 
basicCombustionMixture & composition = thermo.composition()
 
volScalarField & p = thermo.p()
 
volScalarField & b = composition.Y("b")
 Wien displacement law constant: default SI units: [m K]. More...
 
Info<< "Creating field dpdt\"<< endl;volScalarField dpdt(IOobject("dpdt", runTime.timeName(), mesh), mesh, dimensionedScalar(p.dimensions()/dimTime, 0));Info<< "Creating field kinetic energy K\"<< endl;volScalarField K("K", 0.5 *magSqr(U));Info<< "Creating the unstrained laminar flame speed\"<< endl;autoPtr< laminarFlameSpeed > unstrainedLaminarFlameSpeed (laminarFlameSpeed::New(thermo))
 
Info<< "Reading strained laminar flame speed field Su\"<< endl;volScalarField Su(IOobject("Su", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field betav\"<< endl;volScalarField betav(IOobject("betav", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field Lobs\"<< endl;volScalarField Lobs(IOobject("Lobs", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field CT\"<< endl;volSymmTensorField CT(IOobject("CT", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field Nv\"<< endl;volScalarField Nv(IOobject("Nv", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field nsv\"<< endl;volSymmTensorField nsv(IOobject("nsv", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);IOdictionary PDRProperties(IOobject("PDRProperties", runTime.constant(), mesh, IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE));autoPtr< PDRDragModel > drag
 
autoPtr< XiModel > flameWrinkling
 Create the flame-wrinkling model. More...
 
Info<< "Calculating turbulent flame speed field St\"<< endl;volScalarField St(IOobject("St", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), flameWrinkling->Xi() *Su);multivariateSurfaceInterpolationScheme< scalar >::fieldTable fields
 

Function Documentation

◆ pThermo()

Info<< "Reading thermophysical properties\n" << endl;autoPtr<psiuReactionThermo> pThermo ( psiuReactionThermo::New(mesh )

◆ validate()

thermo validate ( args.  executable(),
"ha"  ,
"ea"   
)

◆ rho()

volScalarField rho ( IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE)  ,
thermo.  rho() 
)

◆ turbulence()

Info<< "min(b) = " << min(b).value() << endl;Info<< "\nReading field U\n" << endl;volVectorField U( IOobject ( "U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh);mesh.schemes().setFluxRequired(p.name());Info<< "Creating turbulence model\n" << endl;autoPtr<compressible::RASModel> turbulence ( compressible::New< compressible::RASModel >  rho, U, phi, thermo)

◆ thermophysicalTransport()

Info<< "Creating thermophysical transport model\n" << endl;autoPtr<RASFluidThermophysicalTransportModel> thermophysicalTransport ( RASFluidThermophysicalTransportModel::New(turbulence(), thermo )

◆ if()

if ( composition.  contains"ft")

Definition at line 230 of file createFields.H.

◆ add() [1/2]

fields add ( b  )

◆ add() [2/2]

fields add ( thermo.  he())

◆ addXi()

flameWrinkling addXi ( fields  )

Variable Documentation

◆ thermo

psiuReactionThermo& thermo = pThermo()

Definition at line 7 of file createFields.H.

◆ composition

basicCombustionMixture& composition = thermo.composition()

Definition at line 10 of file createFields.H.

◆ p

volScalarField& p = thermo.p()

Definition at line 25 of file createFields.H.

◆ b

volScalarField& b = composition.Y("b")

Wien displacement law constant: default SI units: [m K].

Definition at line 27 of file createFields.H.

Referenced by greyMean::aCont(), wideBand::aCont(), bufferedAccumulator< scalar >::addToBuffers(), tetrahedron< Point, PointRef >::bounds(), Boussinesq< Specie >::Boussinesq(), lineFace::calcSamples(), points::calcSamples(), tetrahedron< Point, PointRef >::circumCentre(), tetrahedron< Point, PointRef >::circumRadius(), triSurfaceTools::collapseEdges(), volFieldValue::compareScalars(), SSG< BasicMomentumTransportModel >::correct(), Foam::correctContactAngle(), Foam::vectorTools::cosPhi(), cyclicPolyPatch::coupledEdges(), cyclicPolyPatch::coupledPoints(), PengRobinsonGas< Specie >::Cp(), PengRobinsonGas< Specie >::CpMCv(), PengRobinsonGas< Specie >::Cv(), cyclicTransform::cyclicTransform(), LangmuirHinshelwoodReactionRate::ddc(), LangmuirHinshelwoodReactionRate::ddT(), dmat_uniform(), edge::edge(), Foam::eigenValues(), globalIndexAndTransform::encodeTransformIndex(), treeBoundBox::extend(), searchableSphere::findLine(), searchableCylinder::findLine(), searchableSphere::findLineAny(), searchableCylinder::findLineAny(), fieldsExpression::foldAB(), forAll(), coordSet::gather(), Foam::HasherDual(), Foam::HasherInt(), if(), Foam::interpolate(), jenkins_hashbig(), jenkins_hashlittle(), Smagorinsky< BasicMomentumTransportModel >::k(), structuredRenumber::layerLess::layerLess(), Foam::max(), Foam::min(), Foam::minimaxs(), homogeneousMixture< ThermoType >::mixture(), inhomogeneousMixture< ThermoType >::mixture(), egrMixture< ThermoType >::mixture(), slidingInterface::modifyMotionPoints(), Moment< Type >::Moment(), line< Point, PointRef >::nearestDist(), Foam::triIntersect::operator!=(), Foam::operator!=(), procLess::operator()(), minusOpAuto< A, B, R >::operator()(), multiplyOpAuto< A, B, R >::operator()(), plusOpAuto< A, B, R >::operator()(), divideOpAuto< A, B, R >::operator()(), lessProcPatches::operator()(), SortableListDRGEP< Type >::less::operator()(), SortableListEFA< Type >::more::operator()(), structuredRenumber::layerLess::operator()(), globalIndexAndTransform::less::operator()(), LangmuirHinshelwoodReactionRate::operator()(), UList< T >::less::operator()(), UList< T >::greater::operator()(), ListOp< BinaryOp< Type > >::operator()(), ListEqOp< BinaryEqOp< Type > >::operator()(), Boussinesq< Specie >::operator*=(), Foam::operator+(), Foam::operator-(), Foam::operator/(), Foam::operator<(), Foam::operator<<(), Foam::operator<=(), word::operator=(), Foam::operator==(), Foam::operator>=(), Pair< word >::other(), Foam::parseNASCoord(), plane::planePlaneIntersect(), searchableCylinder::points(), Foam::pow(), Foam::triIntersect::protectedDivide(), Foam::vectorTools::radAngleBetween(), Foam::readBool(), Foam::readLabel(), Foam::readULabel(), normal::sample(), PointHit< point >::setDistance(), bufferedAccumulator< scalar >::setSizes(), lduAddressing::size(), Foam::Swap(), swapec(), tetCell::tetCell(), particle::trackToMovingTri(), particle::trackToStationaryTri(), globalIndexAndTransform::transformIndicesForPatches(), triFace::triFace(), vbedg(), XiReactionRate::write(), Boussinesq< Specie >::write(), PengRobinsonGas< Specie >::Z(), and FSD::~FSD().

◆ unstrainedLaminarFlameSpeed

Info<< "Creating field dpdt\n" << endl;volScalarField dpdt( IOobject ( "dpdt", runTime.timeName(), mesh ), mesh, dimensionedScalar(p.dimensions()/dimTime, 0));Info<< "Creating field kinetic energy K\n" << endl;volScalarField K("K", 0.5*magSqr(U));Info<< "Creating the unstrained laminar flame speed\n" << endl;autoPtr<laminarFlameSpeed> unstrainedLaminarFlameSpeed(laminarFlameSpeed::New(thermo))

Definition at line 90 of file createFields.H.

◆ drag

Info<< "Reading strained laminar flame speed field Su\n" << endl;volScalarField Su( IOobject ( "Su", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh);Info<< "Reading field betav\n" << endl;volScalarField betav( IOobject ( "betav", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE ), mesh);Info<< "Reading field Lobs\n" << endl;volScalarField Lobs( IOobject ( "Lobs", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE ), mesh);Info<< "Reading field CT\n" << endl;volSymmTensorField CT( IOobject ( "CT", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE ), mesh);Info<< "Reading field Nv\n" << endl;volScalarField Nv( IOobject ( "Nv", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE ), mesh);Info<< "Reading field nsv\n" << endl;volSymmTensorField nsv( IOobject ( "nsv", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE ), mesh);IOdictionary PDRProperties( IOobject ( "PDRProperties", runTime.constant(), mesh, IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE ));autoPtr<PDRDragModel> drag
Initial value:
(
PDRProperties,
rho,
U,
)
autoPtr< CompressibleMomentumTransportModel > New(const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi, const viscosity &viscosity)
rho
Definition: createFields.H:81
phi
Definition: correctPhi.H:3
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::momentumTransportModel > turbulence(compressible::momentumTransportModel::New(rho, U, phi, thermo))
Definition: createFields.H:94
volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mesh, dimensionedVector(dimVelocity, Zero))

Definition at line 192 of file createFields.H.

Referenced by NicenoKEqn< BasicMomentumTransportModel >::bubbleG(), LaheyKEpsilon< BasicMomentumTransportModel >::bubbleG(), mixtureKEpsilon< BasicMomentumTransportModel >::bubbleG(), and mixtureKEpsilon< BasicMomentumTransportModel >::Ct2().

◆ flameWrinkling

autoPtr<XiModel> flameWrinkling
Initial value:
(
PDRProperties,
Su,
rho,
b,
)
fluidReactionThermo & thermo
Definition: createFields.H:28
autoPtr< CompressibleMomentumTransportModel > New(const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi, const viscosity &viscosity)
volScalarField & b
Wien displacement law constant: default SI units: [m K].
Definition: createFields.H:27
rho
Definition: createFields.H:81
phi
Definition: correctPhi.H:3
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::momentumTransportModel > turbulence(compressible::momentumTransportModel::New(rho, U, phi, thermo))
Definition: createFields.H:94
const tmp< volScalarField::Internal > & Su
Definition: alphaSuSp.H:6

Create the flame-wrinkling model.

Definition at line 202 of file createFields.H.

◆ fields

Info<< "Calculating turbulent flame speed field St\n" << endl;volScalarField St( IOobject ( "St", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), flameWrinkling->Xi()*Su);multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields

Definition at line 228 of file createFields.H.