Convectional term fvm::Sp(fvc::div(phi_), epsilon_) in RNGkEpsilon OF 1.7.1 missing
Hi!
I have a quistion: why this terms are missing in RNGkEpsilon model in epsilon and k equations? Code:
fvm::ddt(epsilon_) the Code:
+ fvm::SuSp(fvc::div(phi_), epsilon_) Only one remark: why there is no such treatment for kEpsilon model in OF1.6ext (in kEpsilon.C)? BTW the residuals normalization used now in OF is not appropriated, because it is connected with the source term and for the turbulent transport equations produces underestimated residuals level. Regards, Alexander 
Quote:
It seems you are trying to treat the convective term in nonconservative form (why?). If you think to: div(kU) = U . grad(k) + k div(U), you are discretizing U . grad(k), since you have the full convective term in conservative form at the line fvm::div(phi_, epsilon_) Best, 
... I ask because IT IS, for example, in realizableKE (incompressible) in OF 1.7.1:
Code:
// Dissipation equation Heh..., it is funny, the flux phi is divergence free... So why is it included? div(UK) k = U grad(k) + k div(U) = U grad(k) for incompressible liquid. That is well known transition between conservative and nonconservatives forms of the convectional term. Alberto, I am sorry, but would you be kind to clarify: is convectional term in RNG somehow different from those in standard and realizable kepsilon? And are you tallking about 1.7.x, which is not 1.7.1? Do you mean some development version? 
So I agree with Alberto,
that U*grad(k) or U*grad(epsilon) or U*grad(omega) are discretized. Hmm..., why? Is it some error? For example, standard kepsilon, transport of kinetic energy is: dk/dt + div(U K) = div(D_k grad(k)) + G  epsilon There is no terms, from which we can get U*grad(k). Any comments? To avoid missunderstanding: I am not trying to put "div(U*k)  k*div(U) = U*grad(k)" in the OF source code, it is all already there))) I only noticed the difference between approaches how convectional term is discretized in RNGkEpsilon.C and in other incompressible RAS models. I recognized it while modifying RAS models for my solidification part. 
Hi,
Quote:
As a consequence, they are implementing U . grad(k) = div(U, k)  k div(U), to avoid using U & grad(k), and use the conservative flux. What is not clear to me is why not implementing directly div(phi, k), as in the original equation from the literature? 
Sometimes for the discretization (e.g. for the finite differences method) it is easier to use nonconservative form of the convectional term
div(U*k) = k*div(U)+U & grad(k) But here we get the wrong term. So I will try direct conservative form and compare the results. 
Quote:
Quote:
div(U*k) = U & grad(k) if div(U) = 0. Now, since numerically it is not exactly zero, they probably enforce the "incompressible form" by using the term U & grad(k) = div(U*k)  k*div(U) This in OpenFOAM becomes fmv::div(phi, k)  fvm::Sp(fvc::div(phi), k), since using phi and not U does not introduce errors, being phi corrected to enforce conservation after solving for p. 
Alberto, thanks! I got it:
Code:
U = U* + delta Code:
div(U*.k) = U*&grad(k) Code:
Err1 = div(U.k)  U*&grad(k) Code:
Err1 = k.div(delta) + delta&grad(k) = k.div(U) + delta&grad(k) div(U.k)  k.div(U) we reduce the calculation error to the delta&grad(k) ))) I am curious, does it worth to use such approach for other scalar equations? E.g. for heat transfer... It would be really nice to here some comment from developers, why they decided to use this method for the turbulence modeling? 
I have the same questions, I feel very confused:
1. Is using "fvm::SuSp(fvc::div(phi_), epsilon_)" better using " fvm::Sp(fvc::div(phi_), epsilon_)"? 2. For other transport equations, should we use this approach? Why it is commented out from RNG kepsilon? 
Hi Daniel!
1. Using SuSp operator is to increase matrix diagonal influence. With Sp operator it goes always to the diagonal and can reduce diagonal elements which makes its solution more problematic for the linear solver. 2. If you check OpenFOAMext such approach is used in all RANS model (*reffering to the source code.... yap! it's there!). I am just using it for RANS models and for energy equation as well. But I did not make any test if it is better for the scalar equation in general. Let's do it! What are your suggestions for the benchmark? It should be some simple case... 
pitzDaily maybe or other BFS flow?

Hmm, Two more questions,
1. When to use source term explicitly and when shall one use it implicitly using Sp or Susp? 2. In ke formulation, why epsilon is solved first instead of k? Thanks 
All times are GMT 4. The time now is 15:25. 