createFields.H File Reference
Include dependency graph for createFields.H:
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions

Info<< "Reading thermophysical properties\n"<< endl;autoPtr< psiuMulticomponentThermo > pThermo (psiuMulticomponentThermo::New(mesh))
 
thermo validate (args.executable(), "ha", "ea")
 
volScalarField rho (IOobject("rho", runTime.name(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.renameRho())
 
Info<< "min(b) = "<< min(b).value()<< endl;Info<< "\nReading field U\n"<< endl;volVectorField U(IOobject("U", runTime.name(), 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))
 
Info<< "Creating thermophysical transport model\n"<< endl;turbulenceThermophysicalTransportModels::unityLewisEddyDiffusivity< RASThermophysicalTransportModel< ThermophysicalTransportModel< compressibleMomentumTransportModel, fluidThermo > >> thermophysicalTransport (turbulence(), thermo, true)
 
Info<< "Creating field dpdt\n"<< endl;volScalarField dpdt(IOobject("dpdt", runTime.name(), 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))
 
 if (composition.contains("ft"))
 
fields add (b)
 
fields add (thermo.he())
 
flameWrinkling addXi (fields)
 

Variables

psiuMulticomponentThermo & thermo = pThermo()
 
basicCombustionMixture & composition = thermo.composition()
 
volScalarField & p = thermo.p()
 
volScalarField & b = composition.Y("b")
 
Info<< "Reading strained laminar flame speed field Su\n"<< endl;volScalarField Su(IOobject("Su", runTime.name(), 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
 
autoPtr< XiModel > flameWrinkling
 Create the flame-wrinkling model. More...
 
Info<< "Calculating turbulent flame speed field St\n"<< endl;volScalarField St(IOobject("St", runTime.name(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), flameWrinkling->Xi() *Su);multivariateSurfaceInterpolationScheme< scalar >::fieldTable fields
 

Function Documentation

◆ pThermo()

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

◆ validate()

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

◆ rho()

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

◆ turbulence()

Info<< "min(b) = " << min(b).value() << endl;Info<< "\nReading field U\n" << endl;volVectorField U( IOobject ( "U", runTime.name(), 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;turbulenceThermophysicalTransportModels::unityLewisEddyDiffusivity< RASThermophysicalTransportModel < ThermophysicalTransportModel < compressibleMomentumTransportModel, fluidThermo > >> thermophysicalTransport ( turbulence()  ,
thermo  ,
true   
)

◆ unstrainedLaminarFlameSpeed()

Info<< "Creating field dpdt\n" << endl;volScalarField dpdt( IOobject ( "dpdt", runTime.name(), 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 )

Referenced by XiFluid::bSolve(), and if().

Here is the caller graph for this function:

◆ if()

if ( composition.  contains"ft")

Definition at line 232 of file createFields.H.

References composition, and fields.

◆ add() [1/2]

fields add ( b  )

◆ add() [2/2]

fields add ( thermo.  he())

◆ addXi()

flameWrinkling addXi ( fields  )

Variable Documentation

◆ thermo

psiuMulticomponentThermo& thermo = pThermo()

Definition at line 7 of file createFields.H.

◆ composition

basicCombustionMixture& composition = thermo.composition()

Definition at line 10 of file createFields.H.

Referenced by if().

◆ p

volScalarField& p = thermo.p()

Definition at line 25 of file createFields.H.

Referenced by if().

◆ b

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

Definition at line 27 of file createFields.H.

Referenced by greyMean::aCont(), wideBand::aCont(), LuoSvendsen::addToBinaryBreakupRate(), bufferedAccumulator< scalar >::addToBuffers(), cubicEqn::b(), linearEqn::b(), quadraticEqn::b(), BarycentricTensor< Cmpt >::BarycentricTensor(), BarycentricTensor2D< Cmpt >::BarycentricTensor2D(), tetrahedron< Point, PointRef >::bounds(), XiFluid::bSolve(), points::calcSamples(), lineFace::calcSamples(), tetrahedron< Point, PointRef >::circumSphere(), triSurfaceTools::collapseEdges(), Foam::combine(), edge::compare(), face::compare(), triFace::compare(), volFieldValue::compareScalars(), SSG< BasicMomentumTransportModel >::correct(), Foam::correctContactAngle(), cyclicPolyPatch::coupledEdges(), cyclicPolyPatch::coupledPoints(), PengRobinsonGas< Specie >::Cp(), PengRobinsonGas< Specie >::CpMCv(), cubicEqn::cubicEqn(), PengRobinsonGas< Specie >::Cv(), cyclicTransform::cyclicTransform(), LangmuirHinshelwoodReactionRate::ddc(), LangmuirHinshelwoodReactionRate::ddT(), cubicEqn::derivative(), quadraticEqn::derivative(), dmat_uniform(), pairPotentialList::dr(), edge::edge(), Foam::eigenValues(), globalIndexAndTransform::encodeTransformIndex(), pairPotentialList::energy(), instant::equal(), cubicEqn::error(), linearEqn::error(), quadraticEqn::error(), treeBoundBox::extend(), tetCell::face(), searchableCylinder::findLine(), searchableSphere::findLine(), searchableCylinder::findLineAny(), searchableSphere::findLineAny(), fieldsExpression::foldAB(), forAll(), pairPotentialList::force(), coordSet::gather(), Foam::HasherDual(), Foam::HasherInt(), if(), jenkins_hashbig(), jenkins_hashlittle(), Smagorinsky< BasicMomentumTransportModel >::k(), linearEqn::linearEqn(), Foam::max(), Foam::min(), Foam::minimaxs(), homogeneousMixture< ThermoType >::mixture(), inhomogeneousMixture< ThermoType >::mixture(), egrMixture< ThermoType >::mixture(), slidingInterface::modifyMotionPoints(), Moment< Type >::Moment(), line< Point, PointRef >::nearestDist(), XiModel::New(), fieldsExpression::opAndStore(), Foam::triIntersect::operator!=(), Foam::operator!=(), Foam::operator&(), plusOpAuto< A, B, R >::operator()(), divideOpAuto< A, B, R >::operator()(), multiplyOpAuto< A, B, R >::operator()(), minusOpAuto< A, B, R >::operator()(), instant::less::operator()(), lessProcPatches::operator()(), UList< T >::less::operator()(), UList< T >::greater::operator()(), procLess::operator()(), structuredRenumber::layerLess::operator()(), SortableListDRGEP< Type >::less::operator()(), SortableListEFA< Type >::more::operator()(), normalLess::operator()(), globalIndexAndTransform::less::operator()(), labelRange::less::operator()(), ListOp< BinaryOp< Type > >::operator()(), LangmuirHinshelwoodReactionRate::operator()(), BinaryOpAdd::operator()(), ListEqOp< BinaryEqOp< Type > >::operator()(), Foam::operator+(), Foam::operator-(), Foam::operator/(), Foam::operator<(), Foam::operator<<(), Foam::operator<=(), Switch::operator=(), bufferedAccumulator< scalar >::operator=(), Foam::triIntersect::operator==(), Foam::operator==(), Foam::operator>(), Foam::operator>=(), Foam::operator>>(), pairPotentialList::pairPotentialFunction(), Vector2D< Cmpt >::perp(), plane::plane(), plane::planePlaneIntersect(), Foam::pow(), Foam::triIntersect::protectedDivide(), quadraticEqn::quadraticEqn(), pairPotentialList::rCut(), pairPotentialList::rCutSqr(), Foam::readBool(), Foam::repeat(), pairPotentialList::rMin(), linearEqn::roots(), face::sameVertices(), Random::sampleAB(), unintegrable::sampleInterval(), Random::scalarAB(), bufferedAccumulator< scalar >::setSizes(), Foam::sqrtSumSqr(), XiFluid::StCorr(), Foam::Swap(), swapec(), tetCell::tetCell(), particle::trackToMovingTri(), particle::trackToStationaryTri(), globalIndexAndTransform::transformIndicesForPatches(), triFace::triFace(), lduAddressing::triIndex(), regIOobject::upToDate(), basicThermo::validate(), cubicEqn::value(), linearEqn::value(), quadraticEqn::value(), vbedg(), XiReactionRate::write(), basicSubGrid::XiEq(), PengRobinsonGas< Specie >::Z(), and bufferedAccumulator< scalar >::~bufferedAccumulator().

◆ drag

Info<< "Reading strained laminar flame speed field Su\n" << endl;volScalarField Su( IOobject ( "Su", runTime.name(), 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,
phi
)
autoPtr< CompressibleMomentumTransportModel > New(const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi, const viscosity &viscosity)
autoPtr< incompressible::momentumTransportModel > turbulence(incompressible::momentumTransportModel::New(U, phi, viscosity))
volScalarField rho(IOobject("rho", runTime.name(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.renameRho())
volVectorField U(IOobject("U", runTime.name(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mesh, dimensionedVector(dimVelocity, Zero))

Definition at line 194 of file createFields.H.

Referenced by NicenoKEqn< BasicMomentumTransportModel >::bubbleG(), LaheyKEpsilon< BasicMomentumTransportModel >::bubbleG(), mixtureKEpsilon< BasicMomentumTransportModel >::bubbleG(), PDRkEpsilon::correct(), mixtureKEpsilon< BasicMomentumTransportModel >::Ct2(), Burns::D(), and Gosman::D().

◆ flameWrinkling

autoPtr<XiModel> flameWrinkling
Initial value:
(
PDRProperties,
Su,
rho,
b,
phi
)
volScalarField & b
Definition: createFields.H:27
tmp< VolField< Type > > Su(const VolField< Type > &su, const VolField< Type > &vf)
Definition: fvcSup.C:44
fluidMulticomponentThermo & thermo
Definition: createFields.H:31

Create the flame-wrinkling model.

Definition at line 204 of file createFields.H.

Referenced by if().

◆ fields