56 coeffsDict_(dict.
subDict(typeName +
"Coeffs")),
58 Comega_(
readScalar(coeffsDict_.lookup(
"Comega"))),
61 WeCrit_(
readScalar(coeffsDict_.lookup(
"WeCrit"))),
62 WeTransition_(
readScalar(coeffsDict_.lookup(
"WeTransition"))),
66 AWe_ = (k21*
sqrt(WeTransition_) - 1.0)/
pow(WeTransition_, 4.0);
88 scalar
pc = spray_.p()[p.
cell()];
93 scalar
rho = fuels.
rho(pc, T, p.
X());
94 scalar sigma = fuels.
sigma(pc, T, p.
X());
95 scalar
mu = fuels.
mu(pc, T, p.
X());
98 scalar rtd = 0.5*Cmu_*mu/(rho*r2);
101 scalar omega2 = Comega_*sigma/(rho*r3) - rtd*rtd;
105 scalar omega =
sqrt(omega2);
106 scalar romega = 1.0/omega;
108 scalar rhog = spray_.rho()[p.
cell()];
109 scalar We = p.
We(Ug, rhog, sigma);
110 scalar Wetmp = We/WeCrit_;
112 scalar
y1 = p.
dev() - Wetmp;
113 scalar y2 = p.
ddev()*romega;
115 scalar a =
sqrt(y1*y1 + y2*y2);
123 phic =
max(
min(phic, 1), -1);
125 scalar phit =
acos(phic);
137 scalar theta =
acos((1.0 - Wetmp)/a);
147 tb = (theta-
phi)*romega;
153 p.
ddev() = -a*omega*
sin(omega*tb + phi);
160 scalar sqrtWe = AWe_*
pow(We, 4.0) + 1.0;
161 scalar Kbr = k1_*omega*sqrtWe;
163 if (We > WeTransition_)
166 Kbr =k2_*omega*sqrtWe;
169 scalar rWetmp = 1.0/Wetmp;
170 scalar cosdtbu =
max(-1.0,
min(1.0, 1.0-rWetmp));
171 scalar dtbu = romega*
acos(cosdtbu);
172 scalar decay =
exp(-Kbr*dtbu);
174 scalar rNew = decay*r;