1 dimensionedScalar
StCorr(
"StCorr", dimless, 1.0);
6 dimensionedScalar Vk(
"Vk", dimVolume, gSum(c*
mesh.V().field()));
7 dimensionedScalar Ak(
"Ak", dimArea, Zero);
9 if (Vk.value() > SMALL)
14 switch(
mesh.nGeometricD())
21 combustionProperties.get<scalar>
23 "ignitionSphereFraction"
27 Ak = sphereFraction*4.0*constant::mathematical::pi
31 /(sphereFraction*4.0*constant::mathematical::pi),
40 dimensionedScalar thickness
42 "ignitionThickness", dimLength, combustionProperties
47 combustionProperties.get<scalar>
49 "ignitionCircleFraction"
53 Ak = circleFraction*constant::mathematical::pi*thickness
60 *constant::mathematical::pi
68 Ak = dimensionedScalar
70 "ignitionKernelArea", dimArea, combustionProperties
77 const volScalarField mgb
79 fvc::div(nf,
b,
"div(phiSt,b)") -
b*fvc::div(nf) + dMgb
81 dimensionedScalar AkEst = gSum(mgb*
mesh.V().field());
83 StCorr.value() =
max(min((Ak/AkEst).value(), 10.0), 1.0);
85 Info<<
"StCorr = " <<
StCorr.value() << endl;
dimensionedScalar StCorr("StCorr", dimless, 1.0)