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.setFluxRequired(p.name());Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::RASModel > turbulence (compressible::New< compressible::RASModel >(rho, U, phi, thermo))
 
 if (composition.contains("ft"))
 
fields add (b)
 
fields add (thermo.he())
 
flameWrinkling addXi (fields)
 

Variables

psiuReactionThermo & thermo = pThermo()
 
basicSpecieMixture & 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.setFluxRequired(p.name());Info<< "Creating turbulence model\n" << endl;autoPtr<compressible::RASModel> turbulence ( compressible::New< compressible::RASModel >  rho, U, phi, thermo)

◆ if()

if ( composition.  contains"ft")

Definition at line 220 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

basicSpecieMixture& 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.

◆ 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 80 of file createFields.H.

Referenced by if().

◆ 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,
)
surfaceScalarField & phi
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::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
Definition: createFields.H:94
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions)
rho
Definition: createFields.H:81
volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mesh, dimensionedVector(dimVelocity, Zero))

Definition at line 182 of file createFields.H.

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

◆ flameWrinkling

autoPtr<XiModel> flameWrinkling
Initial value:
(
PDRProperties,
Su,
rho,
b,
)
zeroField Su
Definition: alphaSuSp.H:1
surfaceScalarField & phi
volScalarField & b
Wien displacement law constant: default SI units: [m K].
Definition: createFields.H:27
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::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
Definition: createFields.H:94
rhoReactionThermo & thermo
Definition: createFields.H:28
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions)
rho
Definition: createFields.H:81

Create the flame-wrinkling model.

Definition at line 192 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 218 of file createFields.H.