56 coeffsDict_(dict.
subDict(typeName +
"Coeffs")),
58 Comega_(
readScalar(coeffsDict_.lookup(
"Comega"))),
59 WeCrit_(
readScalar(coeffsDict_.lookup(
"WeCrit")))
63 const scalar xx0 = 12.0;
64 const scalar rrd100 = 1.0/(1.0-
exp(-xx0)*(1+xx0+
pow(xx0, 2)/2+
pow(xx0, 3)/6));
66 for(label n=0; n<100; n++)
68 scalar
xx = 0.12*(n+1);
69 rrd_[n] = (1-
exp(-xx)*(1 + xx +
pow(xx, 2)/2 +
pow(xx, 3)/6))*rrd100;
92 scalar
pc = spray_.p()[p.
cell()];
97 scalar
rho = fuels.
rho(pc, T, p.
X());
98 scalar sigma = fuels.
sigma(pc, T, p.
X());
99 scalar
mu = fuels.
mu(pc, T, p.
X());
102 scalar rtd = 0.5*Cmu_*mu/(rho*r2);
105 scalar omega2 = Comega_*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/WeCrit_;
114 scalar
y1 = p.
dev() - Wetmp;
115 scalar y2 = p.
ddev()/omega;
117 scalar a =
sqrt(y1*y1 + y2*y2);
125 phic =
max(
min(phic, 1), -1);
127 scalar phit =
acos(phic);
149 scalar theta =
acos((coste-Wetmp)/a);
159 tb = (theta-
phi)/omega;
165 p.
ddev() = -a*omega*
sin(omega*tb + phi);
176 + (4.0/3.0)*
pow(p.
dev(), 2)
177 + rho*r3/(8*sigma)*
pow(p.
ddev(), 2)
182 scalar random = rndGen_.scalar01();
183 while (!found && (n<99))
192 scalar rNew = 0.04*n*rs;