
[Sponsors] 
October 13, 2012, 14:39 
Production term in Finite Volumes Schemes

#1 
New Member
Ivan
Join Date: Aug 2012
Posts: 22
Rep Power: 7 
Sponsored Links
I have a working scheme for LES for a transport equation of a scalar, let's say Phi, done with Finite Volumes schemes in the code I am using. I need to modify this code by adding a production term which consists substantially of two nonlinear terms, let's say k1*Phi^2 and k2*grad(Phi)*grad(Phi), where k1 and k2 are constants. I want to implement it as an implicit scheme, which should be like: A_d * Phi = A_nd * Phi + b where A_d and A_nd are the diagonal matrix and notdiagonal matrix of coefficients and b is the vector of known terms. I am confused on how to implement these new terms. Does someone have an idea or can suggest me some book where to look for? Thank you! 

Sponsored Links 
October 13, 2012, 17:53 

#2  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,416
Rep Power: 39 
Quote:
First, you have to consider the accuracy order of your FV discretization. At second order of accuracy the volume integral can be discretized simply by the mean value formula but I do not suggest that approach since it is not effective in terms of tophat filtering. Better is the multidimensional extension of the trapezoidal rule. If you have an accuracy more than two then your volume integral must be discretized congruently, for example by using shape functions of high degree so that you can integrate them analitically 

October 13, 2012, 18:31 

#3 
New Member
Ivan
Join Date: Aug 2012
Posts: 22
Rep Power: 7 
It is said that the terms in the coefficient matrix are calculated as upwind and then a second order correction is added on the RHS in a deferred way (so that they can use a blending factor I guess)...so the accuracy is second order. Therefore, I'm thinking:
1) On one side, wether using an upwindlike discretization for those terms (plus a second order correction which I would not be sure how to apply) or directly a second order scheme (ex. a central scheme). 2) I'm confused, chosen the scheme, what I should write. For example, since for each time step, a number of subiterations are performed, can I discretize the term, i.e., k1*Phi^2 in cell i as k1*Vol_i*Phi(n1)_i * Phi(n)_i, where n is the current subiteration? And what about the term grad(Phi)^2 

October 13, 2012, 18:35 

#4  
New Member
Ivan
Join Date: Aug 2012
Posts: 22
Rep Power: 7 
Quote:


October 14, 2012, 04:25 

#5  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,416
Rep Power: 39 
Quote:
Furthermore, the volume integral discretized with the mean value formula is second order of accuracy but it is not effective as filter, as you can simply see in the wavenumber space 

October 14, 2012, 08:54 

#6  
New Member
Ivan
Join Date: Aug 2012
Posts: 22
Rep Power: 7 
Quote:
Do you know some book/paper which explain this issue of the Upwind when using sgs model and the issue of the "noneffctiveness" of the filter you were mentioning (I will be required to justify my choices)? 

October 17, 2012, 17:07 

#7  
Senior Member
Join Date: Aug 2011
Posts: 271
Rep Power: 9 
Quote:
k1(Phy_p)^2*Volp +k2*(grad(Phi)_p*grad(Phi)_p)*Volp which is second order too where Phy_p and grad(Phi)_p are evaluated at the cell center P and Volp is the volume of cell P. these terms of course are located in rhight hand side, i.e in b You could also implict the term k1(Phy_p)^2 while adding the quantity k1(Phy_p)*Volp in the main diagonal A_d. It can add stability if large values of Phy occurs. for k2*grad(Phi)*grad(Phi) I do not see any other way instead of expliciting this term. If you have a kind of deferred correction used in the code, it should work fine. 

October 17, 2012, 17:19 

#8  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,416
Rep Power: 39 
Quote:
Hi leflix, consider that you are correct in term of accuracy order in physical space but for the aim of an LES, the volume average should be like a tophat filter. That means, for example in 1D, that the transfer function is in the form sin(kh/2)/kh/2. If you use simply the mean value you will find that the transfer function is the unity function...thus is not effective as filter 

October 17, 2012, 17:24 

#9  
Senior Member
Join Date: Aug 2011
Posts: 271
Rep Power: 9 
Quote:
Hi Filippo, ok my contribution here was just about discretization principles. From LES concerns I do trust you ! 

October 17, 2012, 17:36 

#10  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,416
Rep Power: 39 
Quote:
Regards Filippo 

October 17, 2012, 17:39 

#11 
Senior Member
Join Date: Aug 2011
Posts: 271
Rep Power: 9 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Plotting a contour of the turbulence production term in CFXPost  kellio9  CFX  4  March 19, 2013 12:34 
production term on staggered grid  behzadmaha  Main CFD Forum  0  April 24, 2011 16:30 
ATTENTION! Reliability problems in CFX 5.7  Joseph  CFX  14  April 20, 2010 15:45 
Finite volumes vs. finite elements: applications  Lionel S.  Main CFD Forum  10  January 22, 2007 11:51 
Numerical method for the kepsilon production term  carlslack  Main CFD Forum  0  January 23, 2003 17:13 
Sponsored Links 