The turbulent kinetic energy transport equation, \(k\), [[41], Eq. 3]:
\[ \partial_t k + u_i \partial_{x_i} k = P - \epsilon + \partial_{x_j} \left\{ \left( \nu + \frac{\nu_t}{\sigma_k} \right) \partial_{x_j} k \right\} \]
where
| \(k \) | = | Turbulent kinetic energy [ \(\text{m}^2 \text{s}^{-3}\)] |
| \(P \) | = | Turbulent kinetic energy production rate [ \(\text{m}^2 \text{s}^{-3}\)] |
| \(\epsilon \) | = | Turbulent kinetic energy dissipation rate [ \(\text{m}^2 \text{s}^{-3}\)] |
| \(\nu \) | = | Kinematic viscosity of fluid [ \(\text{m}^2 \text{s}^{-1}\)] |
| \(\nu_t \) | = | Turbulent viscosity [ \(\text{m}^2 \text{s}^{-1}\)] |
| \(\sigma_k \) | = | Turbulent Prandtl number for \(k\) [-] |
The turbulent kinetic energy dissipation rate transport equation, \(\epsilon\), [[41], Eq. 4]:
\[ \partial_t \epsilon + u_i \partial_{x_i} \epsilon = \frac{C_{\epsilon_1} P}{T} - \frac{C_{\epsilon_2} \epsilon}{T} + \partial_{x_j} \left\{ \left( \nu + \frac{\nu_t}{\sigma_\epsilon} \right) \partial_{x_j} \epsilon \right\} \]
where
| \(T \) | = | Modelled turbulent time scale [ \(s\)] |
| \(C_{\epsilon_*} \) | = | Empirical model constants [-] |
| \(\sigma_\epsilon \) | = | Turbulent Prandtl number for \(\epsilon\) [-] |
The normalised wall-normal fluctuating velocity scale transport equation, \(\phi\), [[41], Eq. 17]:
\[ \partial_t \phi + u_i \partial_{x_i} \phi = f - P \frac{\phi}{k} + \frac{2 \nu_t}{k \sigma_k} \partial_{x_j} \phi \partial_{x_j} k + \partial_{x_j} \left\{ \left( \frac{\nu_t}{\sigma_k} \right) \partial_{x_j} \phi \right\} \]
where
| \(f \) | = | Elliptic relaxation factor [ \(s^{-1}\)] |
| \(\phi \) | = | Normalised wall-normal fluctuating velocity scale [-] |
| \(\sigma_\phi \) | = | Turbulent Prandtl number for \(\phi\) [-] |
The elliptic relaxation factor equation, \(f\), [[41], Eq. 18]:
\[ L^2 \partial^2_{x_j} f - f = \frac{1}{T} (C_1 - 1) \left[ \phi - \frac{2}{3} \right] - C_2 \frac{P}{k} - 2 \frac{\nu }{k} \partial_{x_j} \phi \partial_{x_j} k - \nu \partial^2_{x_j} \phi \]
where
| \(L \) | = | Modelled turbulent length scale [ \(m\)] |
| \(C_{f_*} \) | = | Empirical model constants [-] |
The turbulent time scale equation, \(T\), [[41], Eq. 7]:
\[ T = \max \left[ \frac{k}{\epsilon}, 6.0 \sqrt{\frac{\nu}{\epsilon}} \right] \]
The turbulent length scale equation, \(L\), [[41], Eq. 7]:
\[ L = C_L \max \left[ \frac{k^{1.5}}{\epsilon}, C_\eta \left( \frac{\nu^3}{\epsilon} \right)^{0.25} \right] \]
where
| \(C_L \) | = | Empirical model constant [-] |
| \(C_\eta \) | = | Empirical model constant [-] |
The turbulent viscosity equation, \(\nu_t\), [[41], p. 173]:
\[ \nu_t = C_\mu \phi k T \]
where
| \(C_{\mu} \) | = | Model coefficient for the turbulent viscosity [-] |
| \(\nu_t \) | = | Turbulent viscosity [ \(\text{m}^2 \text{s}^{-1}\)] |
The main differences between the implemented and original governing equations are as follows:
nu into DphitEff even though it is not present in [[41], Eq. 17] (This change has provided higher level of resemblance to benchmarks for the tests considered, particularly for the peak skin friction, yet pressure-related predictions were unaffected. Users can switch off nu in DphitEff by using includeNu entry in kEpsilonPhitFCoeffs as shown below in order to follow the reference paper thereat. includeNu is left true by default.)The implementation of the \(\epsilon\) transport equation:
\[ \ddt{\alpha \rho \epsilon} + \div \left(\alpha \rho \u \epsilon \right) - \laplacian \left(\alpha \rho D_\epsilon \epsilon \right) = \alpha \rho C_{\epsilon_1} \frac{G}{T} - \left( \left\{ \frac{2}{3} C_{\epsilon_1} \right\} \left\{ \alpha \rho \div \u \right\} \epsilon \right) - \left( \alpha \rho \frac{C_{\epsilon_2}}{T} \epsilon \right) + S_{\text{fvOptions}} \]
where
| \(\alpha \) | = | Phase fraction of the given phase [-] |
| \(\rho \) | = | Density of the fluid [ \( \text{kg} m^{-3} \)] |
| \(S_{\text{fvOptions}} \) | = | Source terms introduced by fvOptions dictionary |
The implementation of the \(k\) transport equation:
\[ \ddt{ \alpha \rho k } + \div \left( \alpha \rho \u k \right) - \laplacian \left( \alpha \rho D_k k \right) = \alpha \rho G - \left( \left[ \frac{2}{3} \alpha \rho \div \u k \right] \right) - \left( \frac{\alpha \rho}{T} k \right) + S_{\text{fvOptions}} \]
where
| \(G \) | = | Anisotropic contribution part of the turbulent kinetic energy production rate, i.e. \( P = G - 2/3 k \div \u \) |
The implementation of the \(f\) elliptic relaxation equation:
\[ - \laplacian f = - \left( \frac{1}{L^2} f \right) - \left( (C_{f_1} - 1) \frac{\phi - \frac{2}{3}}{T} - \frac{C_{f_2} G}{k} + C_{f_2} \frac{2}{3} \div \u - \frac{2 \nu (\nabla \phi \cdot \nabla k )}{k} - \nu \laplacian \phi \right)\frac{1}{L^2} \]
The implementation of the \(\phi\) transport equation:
\[ \ddt{\alpha \rho \phi} + \div \left(\alpha \rho \u \phi \right) - \laplacian \left(\alpha \rho D_\phi \, \phi \right) = \alpha \rho f - \left( \alpha \rho \left\{ \frac{G}{k} - \frac{2}{3} \div \u - \frac{2 \nu (\nabla \phi \cdot \nabla k )}{k \sigma_k \phi } \right\} \phi \right) + S_{\text{fvOptions}} \]
The implementation of the \( T \) equation:
\[ T = \max \left( \frac{k}{\epsilon}, C_T \frac{\sqrt{\max\{\nu, 0\}}}{\epsilon} \right) \]
The implementation of the \(L\) equation:
\[ L = C_L \max \left( \frac{k^{1.5}}{\epsilon}, C_\eta \left[ \frac{(\max\{\nu, 0\})^3}{\epsilon} \right]^{0.25} \right) \]
The model coefficients are [[41], Eqs. 19-20]:
\[ C_\mu = 0.22; \quad C_{\epsilon_1} = 1.4(1.0 + 0.05\sqrt{1.0/\phi} ); \quad C_{\epsilon_2} = 1.9; \quad C_T = 6.0; \quad C_L = 0.25;\\ C_{f_1} = 1.4; \quad C_{f_2} = 0.3; \quad C_\eta = 110.0; \quad \sigma_k = 1.0; \quad \sigma_\epsilon = 1.3; \quad \sigma_\phi = 1.0. \]
Please note that [[41], p. 176] stated that the model constants above have been calibrated to match the flow statistics of the smooth-wall plane channel flow at \(\text{Re}_\tau = 395\).
The model can be enabled by using constant/turbulenceProperties dictionary:
simulationType RAS;
RAS
{
// Mandatory entries
RASModel kEpsilonPhitF;
// Optional entries
turbulence on;
printCoeffs on;
kEpsilonPhitFCoeffs
{
includeNu true; // include nu in DphitEff()
Cmu 0.22; // Turbulent viscosity constant
Ceps1a 1.4; // Model constant for epsilon
Ceps1b 1.0; // Model constant for epsilon
Ceps1c 0.05; // Model constant for epsilon
Ceps2 1.9; // Model constant for epsilon
Cf1 1.4; // Model constant for f
Cf2 0.3; // Model constant for f
CL 0.25; // Model constant for L
Ceta 110.0; // Model constant for L
CT 6.0; // Model constant for T
sigmaK 1.0; // Turbulent Prandtl number for k
sigmaEps 1.3; // Turbulent Prandtl number for epsilon
sigmaPhit 1.0; // Turbulent Prandtl number for phit = sigmaK
}
Please note that \( \phi \) is a symmetric, and \( f \) is an asymmetric matrix.
For \( k \) and \( \epsilon \) fields, the initial conditions can be estimated by using the recommendations made for the kEpsilon model.
For \( \phi \):
object phit; dimensions [0 0 0 0 0 0 0]; internalField uniform 0.66;
For \( f \):
object f; dimensions [0 0 -1 0 0 0 0]; internalField uniform 1.0;
For \( k \) and \( \epsilon \) fields, please use the boundary conditions that are available for the kEpsilon model.
For \( \phi \) and \( f \):
Inlet/Outlet:
Walls:
For example, \( \phi \):
object phit;
dimensions [0 0 0 0 0 0 0];
wall
{
type fixedValue;
value uniform 1e-10;
}
"(inlet|outlet)"
{
type zeroGradient;
}
For example, \( f \):
object f;
dimensions [0 0 -1 0 0 0 0];
wall
{
type fixedValue;
value uniform 0;
}
"(inlet|outlet)"
{
type zeroGradient;
}
Source code:
References:
| Would you like to suggest an improvement to this page? | Create an issue |
Copyright © 2019-2020 OpenCFD Ltd.
Licensed under the Creative Commons License BY-NC-ND 
1.8.17