5 molecules.potential().nIds(),
11 molecules.potential().nIds(),
17 molecules.potential().nIds(),
23 molecules.potential().nIds(),
37 molecules.potential().nIds()
42 Info<<
" Creating number density field for " 43 << molecules.potential().idList()[rN] <<
endl;
52 "rhoN_" + molecules.potential().idList()[rN],
67 Info <<
" Creating total number density field" <<
endl;
83 totalRhoN.primitiveFieldRef() = 0;
84 totalRhoN.correctBoundaryConditions();
92 molecules.potential().nIds()
97 Info<<
" Creating mass density field for " 98 << molecules.potential().idList()[rM] <<
endl;
107 "rhoM_" + molecules.potential().idList()[rM],
122 Info <<
" Creating total mass density field" <<
endl;
138 totalRhoM.primitiveFieldRef() = 0;
139 totalRhoM.correctBoundaryConditions();
147 molecules.potential().nIds()
152 Info<<
" Creating velocity field for " 153 << molecules.potential().idList()[v] <<
endl;
162 "velocity_" + molecules.potential().idList()[v],
177 Info <<
" Creating total velocity field" <<
endl;
220 molecules.potential().nIds()
225 Info<<
" Creating temperature field for " 226 << molecules.potential().idList()[t] <<
endl;
235 "temperature_" + molecules.potential().idList()[t],
250 Info <<
" Creating total temperature field" <<
endl;
266 totalTemperature.primitiveFieldRef() = 0;
267 totalTemperature.correctBoundaryConditions();
276 molecules.potential().nIds()
281 Info<<
" Creating mean kinetic energy field for " 282 << molecules.potential().idList()[mKE] <<
endl;
291 "meanKE_" + molecules.potential().idList()[mKE],
298 dimensionSet(1, 2, -2, 0, 0, 0, 0),
306 Info <<
" Creating total mean kinetic energy field" <<
endl;
319 dimensionSet(1, 2, -2, 0, 0, 0, 0),
323 totalMeanKE.correctBoundaryConditions();
Info<< " Creating total mass density field"<< endl;volScalarField totalRhoM(IOobject("rhoM_total", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimMass/dimVolume, "zeroGradient");totalRhoM.primitiveFieldRef()=0;totalRhoM.correctBoundaryConditions();PtrList< volVectorField > allSpeciesVelocity(molecules.potential().nIds())
List< scalarField > allSpeciesM_RU(molecules.potential().nIds(), scalarField(mesh.nCells(), 0.0))
dimensioned< vector > dimensionedVector
Dimensioned vector obtained from generic dimensioned type.
Ostream & endl(Ostream &os)
Add newline and flush stream.
forAll(allSpeciesRhoN, rN)
volVectorField vectorField(fieldObject, mesh)
List< scalarField > allSpeciesVelocityMagSquaredSum_RU(molecules.potential().nIds(), scalarField(mesh.nCells(), 0.0))
List< vectorField > allSpeciesVelocitySum_RU(molecules.potential().nIds(), vectorField(mesh.nCells(), Zero))
GeometricField< vector, fvPatchField, volMesh > volVectorField
Info<< " Creating total temperature field"<< endl;volScalarField totalTemperature(IOobject("temperature_total", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimTemperature, "zeroGradient");totalTemperature.primitiveFieldRef()=0;totalTemperature.correctBoundaryConditions();PtrList< volScalarField > allSpeciesMeanKE(molecules.potential().nIds())
Info<< " Creating total velocity field"<< endl;volVectorField totalVelocity(IOobject("velocity_total", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimensionedVector("zero", dimVelocity, Zero));PtrList< volScalarField > allSpeciesTemperature(molecules.potential().nIds())
GeometricField< scalar, fvPatchField, volMesh > volScalarField
const dimensionSet dimVolume(pow3(dimLength))
const dimensionSet dimTemperature(0, 0, 0, 1, 0, 0, 0)
List< scalarField > allSpeciesN_RU(molecules.potential().nIds(), scalarField(mesh.nCells(), 0.0))
volScalarField scalarField(fieldObject, mesh)
Internal::FieldType & primitiveFieldRef()
Return a reference to the internal field.
Info<< " Creating total number density field"<< endl;volScalarField totalRhoN(IOobject("rhoN_total", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimless/dimVolume, "zeroGradient");totalRhoN.primitiveFieldRef()=0;totalRhoN.correctBoundaryConditions();PtrList< volScalarField > allSpeciesRhoM(molecules.potential().nIds())
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
Info<< nl<< "Creating fields."<< endl;PtrList< volScalarField > allSpeciesRhoN(molecules.potential().nIds())
const dimensionSet dimMass(1, 0, 0, 0, 0, 0, 0)
const dimensionSet dimVelocity