Commit d50c5ce8 authored by Andrew Heather's avatar Andrew Heather
Browse files

ENH: DEShybrid - added robustness for case that nuEff goes negative (typically...

ENH: DEShybrid - added robustness for case that nuEff goes negative (typically from negative denisty in compressible cals).  See #1224
parent f330921c
......@@ -197,13 +197,26 @@ class DEShybrid
CH3_*Omega*max(S, Omega)
/max(0.5*(sqr(S) + sqr(Omega)), sqr(OmegaLim_/tau0_))
);
const volScalarField K
(
max(Foam::sqrt(0.5*(sqr(S) + sqr(Omega))), 0.1/tau0_)
);
const volScalarField lTurb(Foam::sqrt(nuEff/(pow(0.09, 1.5)*K)));
const volScalarField lTurb
(
Foam::sqrt
(
max
(
nuEff/(pow(0.09, 1.5)*K),
dimensionedScalar("l0", sqr(dimLength), 0)
)
)
);
const volScalarField g(tanh(pow4(B)));
const volScalarField A
(
CH2_*max
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment