Elliptic Blending Model

UNDER CONSTRUCTION

The model implemented here is the one described in Thielen et al. (Int Journal of Heat and Mass Transfer, 48, 2005, pp.1583-1598).

The model solves the full set Reynolds stresses %$\overline{u_i u_j}$%, a dissipation equation %$\varepsilon$% and an elliptic equation for %$\alpha$%.

The Reynolds stress equation can be written as:

%BEGINLATEX% \begin{align} \frac{D\overline{u_i u_j}}{Dt} = &P_{ij} + {\Phi^*}_{ij} - \varepsilon_{ij}+\frac{\partial}{\partial x_k}\left[ (\nu \delta_{kl} + C_s\overline{u_ku_l}\tau) \frac{\partial \overline{u_i u_j}}{\partial x_l} \right] \\ \frac{D\varepsilon}{Dt} = & \frac{C'_{\varepsilon 1}P - C_{\varepsilon 2} \varepsilon}{\tau}+\left[ (\nu \delta_{kl} + C_{\varepsilon}\overline{u_ku_l}\tau) \frac{\partial \varepsilon }{\partial x_l} \right] \end{align} %ENDLATEX%

The main idea of the model is to blend the near wall part with the farfiled one via %$\alpha$%, which is defined as:

%BEGINLATEX% \begin{equation} \alpha-L^2\nabla^2\alpha = 1 \end{equation} %ENDLATEX%

With this parameter, the equation for the velocity-pressure-gradient correlation %$\Phi^*_{ij}$% can be written as:

%BEGINLATEX% \begin{equation} \Phi^*_{ij} = (1-\alpha^2)\Phi^w_{ij}+\alpha^2\Phi^h_{ij} \end{equation} %ENDLATEX%

The "homogeneous" part is taken from the SSG (Speziale et al.) %BEGINLATEX% \begin{align*} \Phi^h_{ij} =& -\left( C_1 +C_2\frac{P}{\varepsilon}\right) \varepsilon a_{ij} + C_3kS_{ij}\&+C_4k \left(a_{ik}S_{jk} +a_{jk}S_{ik} -\frac{2}{3}\delta_{ij}a_{lk}S_{kl} \right)\&+C_5k \left( a_{ik}\Omega_{jk} + a_{jk}\Omega_{ik} \right) \end{align*} %ENDLATEX%

And the near wall part is derived to satisfy the near wall asymptotic behaviour and stress budgets:

%BEGINLATEX% \begin{equation*} \Phi^w_{ij} = -5\frac{\varepsilon}{k} \left( \overline{u_i u_k} n_j n_k + \overline{u_j u_k} n_i n_k - \frac{1}{2} \overline{u_k u_l}n_k n_l (n_i n_j + \delta_{ij} ) \right) \end{equation*} %ENDLATEX%

with the following definitions:

%BEGINLATEX% \begin{align*} &a_{ij} = \frac{\overline{u_i u_j}}{k}-\frac{2}{3}\delta_{ij}; \quad S_{ij} = \frac{1}{2}\left(\frac{\partial U_i}{\partial x_j} + \frac{\partial U_j}{\partial x_i} \right); \quad \Omega_{ij} = \frac{1}{2}\left(\frac{\partial U_i}{\partial x_j} - \frac{\partial U_j}{\partial x_i} \right) \\ &\mathbf{n} = \frac{\nabla \alpha}{\parallel\nabla \alpha\parallel}; \quad \tau = \max \left( \frac{k}{\varepsilon}; \, C_{\tau} \sqrt[2]{\frac{\nu}{\varepsilon}} \right); \quad L=C_L\max \left( \frac{k^{3/2}}{\varepsilon}; \, C_{\eta} \sqrt[4]{\frac{\nu^3}{\varepsilon}} \right) \end{align*} %ENDLATEX%

The constants used are:

%BEGINLATEX% \begin{tabular}{|c c c c c c c c c c c|} \hline $C_s$ & $C_{\varepsilon}$ & $C_{\varepsilon}^0$ & $C_{\varepsilon2}$ & $C_1$ & $C_2$ & $C_4$ & $C_5$ & $C_{\tau}$ & $C_L$ & $C_{\eta}$ \\ \hline 0.21&0.18&1.44&1.83&1.7&0.9&0.625&0.2&6.0&0.161&80 \\ \hline \end{tabular} \\ \centering $C_{\varepsilon 1} = C_{\varepsilon 1}^0 \left(1+0.03(1-\alpha^2)\sqrt{\frac{k}{\overline{u_i u_j} n_i n_j}}\right)$ %ENDLATEX%

-- JuanUribe - 2009-11-05


Current Tags:
create new tag
, view all tags
Topic attachments
I Attachment Action Size Date Who CommentSorted ascending
elsegz ebm_channel.tar.gz manage 283.5 K 2009-11-05 - 12:41 JuanUribe File for the EBM for a channel flow calculation (v1.4.0)
Topic revision: r3 - 2009-11-06 - 10:17:22 - FlavienBillard
Saturne Web
10 Dec 2019
Links
Home
Download
FAQ
News
Test Cases
Developments
Publications

tip Saturne Tip of the Day
No graphical interface
To not use the graphical interface, create the cases with the noihm flag, i.e. cree sat STUDY CASE ... Read on Read more

 

Code_Saturne collaborative website
@ the University of Manchester
Copyright & by the contributing authors. Unless noted otherwise, all material on this web site is the property of the contributing authors.