19 (
b*c*mgb)().weightedAverage(
mesh.V())
20 /((
b*c)().weightedAverage(
mesh.V()) + SMALL)
28 nfVec /= (
mag(nfVec) + dMgb);
41 mesh.surfaceInterpolation::deltaCoeffs()
43 ).value()*runTime.deltaT().value();
45 Info<<
"Max St-Courant Number = " << StCoNum <<
endl;
108 )*(Xi + scalar(1))/(2*Xi);
120 if (SuModel ==
"unstrained")
124 else if (SuModel ==
"equilibrium")
128 else if (SuModel ==
"transport")
134 (sigmas*SuInf*(Su0 - SuInf) +
sqr(SuMin)*sigmaExt)
135 /(
sqr(Su0 - SuInf) +
sqr(SuMin));
166 if (XiModel ==
"fixed")
170 else if (XiModel ==
"algebraic")
176 (scalar(1) + (2*XiShapeCoef)*(scalar(0.5) -
b))
177 *XiCoef*
sqrt(up/(
Su + SuMin))*Reta;
179 else if (XiModel ==
"transport")
187 scalar(1.001) + XiCoef*
sqrt(up/(
Su + SuMin))*Reta;
191 + (scalar(1) + (2*XiShapeCoef)*(scalar(0.5) -
b))
192 *(XiEqStar - scalar(1.001));
228 Info<<
"max(XiEq) = " <<
max(XiEq).value() <<
endl;
237 Info<<
"Combustion progress = "
238 << 100*(scalar(1) -
b)().weightedAverage(
mesh.V()).value() <<
"%"