Elliptic-blending k - ε model

Formulation
All-y + wall-treatment
Those terms needs to be defined at the first cell away from the
wall:
Remarks:
Pk =
(U + uk) 2 νt 0.95(1−νt)
2 classical u∗
u4 ∗ dU
ν
+
dy +
1− dU+
dy +
iterative u∗
ε = Γ νk
+(1−Γ)
y2 E = − u2 k
µ 0.0875exp
Cµu∗k
2κy
−20.66
◮ U + is taken from Reichardt’s law:
U + = 1
κ ln(1+κy+ )+C
ln(y + )
ln(10) −1.01
2
1−e
with C = 1
κ ln
E 1 y +
κ and b = mκ 1
2 C +
y + m
−y +
y + m + y+
y + m
.
e −by+