Number of equation by using eddy viscosity turbulence models
1 Attachment(s)
Hello,
i'm trying to understand the theory behind the eddy viscosity turbulence modelling. First i'm doing the RANS approach on the conservation equations and i get the Reynolds Stress tensor like in the picture (Reynold_Stress_Tensor) below. Each line characterize the additonal terms that i get in the three momentum equation. By choosing a eddy viscosity turbulence model, i use the approach of Boussinesq for modelling the expressions of the Reynolds Stress Tensor. There is an unknown variable, the turbulent viscosity. The definition of the turbulent viscosity depends on the choosen turbulence model. Let's choose the k-epsilon model for example. https://www.cfd-online.com/Wiki/Stan...-epsilon_model In the equations of k and epsilon, you can find the indices i and j which go from 1 to 3 each (direction x,y,z). If i=1, j could be 1,2,3, so i have 3 equations in total for k in x-direction ? All in all i have 9 equations for k and 9 equations for epsilon? Than i have lots of opportunities to calculate the turbulent viscosity and the entries of the Reynolds Stress Tensor. Or do i have just 9 additonal equations (one for each entrie of the Reynolds Stress Tensor). For example i want to calculate the first entry. Than, i and j are 1 and it "goes down" to the equation of the turbulence model? So the main question is, how many equations do i get in total. I hope that i describe my problem in a good way. It would be great if someone could explain it to me. Thanks a lot, bumper |
You get 3 scalar momentum equations because it's a 3-dim vector. Really you could equivalently say you have 1 momentum equation.
You get 1 equation for k and 1 equation for epsilon because they're both scalars. The fact that you see an i,j, or k in the k and epsilon equations doesn't really mean you get 3 new equations each. If you write out the entire equations, you'll find there is just 1 equation each. You have to follow the summation rules for the notation being used. Even though the Stress Tensor has 9 entries in the matrix, you still get 3 momentum equations. You have to pay attention to the products of vectors and tensors. What you need to get the turbulent viscosity is k and epsilon. Once the turbulent viscosity is known, you get all the components of the reynolds stress tensor from the boussinesq hypothesis. You get k and epsilon by solving their transport equations, respectively. |
Okay, i have three momentum equations, that's clear. 3 entries (one line) of the Reynolds Stress Tensor belong to one momentum equation.
So the i,j or k in the equation for the energy or dissipation doesn"t mean, that i get a 3dimensional equation. There is only one equation for the energy k and one for the dissipation rate. Which summation rules do you mean? I see just a partial derivation. That's the reason why i was expacting 3 equation each. Can you show me how the entire equation looks like? Thanks a lot for your help. |
Quote:
Can you suggest a book with numerical schemes solving k and epsilon equations? Many of CFD books describe for example SIMPLE, SIMPLER, PISO etc algorithms but do not mention numerical methods for solving for example k-epsilon equations. |
Quote:
As an example, the momentum equations look like this when you expand all the summations. Yes this is the navier-stokes & not the RANS. The stress tensor and reynolds stress tensor are similar anyway and the difference is mainly in the effective viscosity. Despite there being 9 combinations of indices, you get only an x,y,z equations. There's no magic here, you just have to follow the rules for the notation being used. Quote:
Other equations like energy, k, epsilon, etc. don't have this issue. All you need is some type of linear solver (i.e. solve Ax=b). You just need a way to solve this linear system. The dumb way is to invert A. Gauss-Seidel algorithm is another archetype. Even when you execute the SIMPLE algorithm, you will use similar solvers to solve the individual pieces of the SIMPLE algorithm because the first step in SIMPLE is to solve the momentum equation with an assumed pressure field. |
I think that a brief observation is required in case of incompressible flows.
The eddy viscosity model is written for the deviatoric part of the tensor, the isotropic part (the trace) being incorporated into the "pressure". From that you can write the three scalar equations |
1 Attachment(s)
Quote:
Thanks for your help. |
I am not sure to understand the key of these doubts...In the eddy viscosity model, you substitute the full unknown tensor with the modelled tensor. One does not pretend to substitute each entry of the Reynolds tensor with a modelled entry.
At the end you have a term like this: Div(2*mu_total*S) being S the symmetric part of the velocity gradient with zero trace. Now, to get the contribution to one of the scalar momentum equations, just project the term along a direction ei |
Quote:
Expressions in equation for x mometum: First line Expressions in equation for y mometum: Second line And so on... |
Quote:
So, you are not interested to the eddy viscosity modelling but only on the unresolved Reynolds stress? In such case just consider that this tensor is nothing but as same as you have in the classical convective term Div (vv). Along the x-direction you write i.Div (v'v'), along y j.Div (v'v') and so on (I did not write the averaging). for example, the x component: d/dx(u'u')+d/dy(v'u')+d/dz(w'u') |
Sure, i"m interested to the eddy viscosity modelling but i want to understand the Reynolds Stress tensor and the equations behind.
If the tensor is nothing but as same as you have in the classical convective term Div, i'm right with the lines of the Reynolds stress tensor and the equations of momentum, right? |
Quote:
yes, the expression I wrote above is the same in the figure you posted |
To sum it up and to check if i unstand it in the right way, the Reynolds Stress Tensor is the summary of additonal entries which i get through the Reynolds-approach in the equations of momentum right?
|
Quote:
yes, is exactly as same expression as when you write the convective terms, that is d/dxi(u'i u'j). |
Perfect, than thanks a lot for your help
|
Quote:
Thank you. So lets take the k equation in k-epsilon model. Can we discretize this equation with standard schemes applied for convective-diffusive terms like ADI and solve it for example with Gauss-Seidel or Jacobi? |
Quote:
Yes. Usually in CFD the advection terms are the trouble maker in terms of accuracy and stability, but specifically the advection term in the momentum equation (it's non-linear) is where all the attention goes. The advection terms in k and epsilon equations are linear with respect to k and epsilon and much easier to solve and give much fewer headaches. That's why we don't really talk about them much. If you can solve the much tougher momentum equation, the rest are easier. |
Quote:
Well, actually there is a problem in the transport of k as it is assumed to be a non-negative variable and a monotonicity-preserving scheme should be properly used. |
Quote:
Thank you. |
All times are GMT -4. The time now is 22:44. |