52 includeOscillation_(dict_.lookup(
"includeOscillation")),
53 TABcoeffsDict_(dict.
subDict(
"TABCoeffs")),
60 if (includeOscillation_)
63 yDot0_ =
readScalar(TABcoeffsDict_.lookup(
"yDot0"));
64 TABComega_ =
readScalar(TABcoeffsDict_.lookup(
"Comega"));
65 TABCmu_ =
readScalar(TABcoeffsDict_.lookup(
"Cmu"));
66 TABWeCrit_ =
readScalar(TABcoeffsDict_.lookup(
"WeCrit"));
87 if(includeOscillation_)
91 scalar
pc = spray_.p()[p.
cell()];
92 scalar r = 0.5 * p.
d();
96 scalar
rho = fuels.
rho(pc, T, p.
X());
97 scalar sigma = fuels.
sigma(pc, T, p.
X());
98 scalar
mu = fuels.
mu(pc, T, p.
X());
101 scalar rtd = 0.5*TABCmu_*mu/(rho*r2);
104 scalar omega2 = TABComega_ * sigma /(rho*r3) - rtd*rtd;
109 scalar omega =
sqrt(omega2);
110 scalar rhog = spray_.rho()[p.
cell()];
111 scalar We = p.
We(Ug, rhog, sigma);
112 scalar Wetmp = We/TABWeCrit_;
114 scalar
y1 = p.
dev() - Wetmp;
115 scalar y2 = p.
ddev()/omega;
118 scalar c =
cos(omega*deltaT);
119 scalar s =
sin(omega*deltaT);
120 scalar
e =
exp(-rtd*deltaT);
121 y2 = (p.
ddev() + y1*rtd)/omega;
123 p.
dev() = Wetmp + e*(y1*c + y2*s);
131 p.
ddev() = (Wetmp-p.
dev())*rtd + e*omega*(y2*c - y1*s);