
[Sponsors] 
January 17, 2002, 09:20 
Setting PRT as a function of stratification

#1 
Guest
Posts: n/a

Hi!
I would like to set the turbulent Prandtl number as a function of the stratification in the ke model. One way of doing this would be to set the following in Q1. 1. To calculate ENUT in GROUND from KE and EP with my own formulas for PRT: ENUT=GRND 2. To set the buoyancy term in the equations for KE and EP in GROUND using my own values for PRT: PATCH(patchname,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP) COVAL(patchname,KE,FIXFLU,GRND1) COVAL(patchname,EP,FIXLFU,GRND2) Assuming the above is correct, my remaining problem is the following. Is the section in GROUND where I set my own formula for ENUT called for each variable, or only once per slab/timestep/whatever? Can I test on INDVAR so I get different PRTformulas for the viscosity and the diffusivity? I run 3.2 on SGI IRIX. Thanks, Olof 

January 22, 2002, 05:57 
Re: Setting PRT as a function of stratification

#2 
Guest
Posts: n/a

Hi OlOf,
I may be wrong, but if my recollection is correct, I don't think you can do it this way.  for a start, I am doubting that ENUT can be set indvidually and differently for a given variable. ENUT is ENUT, in the same way as RHO is RHO.  there is no provision for a user's coding of PRT (turbulent) in ground (but please check). Therefore, the only way left to alter the PRT value in ground seems to be to set PRT(indvar in question) to a large number in order to cancel out the builtin turbulent diffusive coefficient. Then, you can rewrite your own coefficient via the PRNDTL(indvar in question), for which access is definitively provided in ground. This coding should hence reflect both the laminar and turbulent contribution. Well, at least, sthing like that. Good luck, Hervé. 

January 22, 2002, 12:17 
Re: Setting PRT as a function of stratification

#3 
Guest
Posts: n/a

Thanks, Hervé!
I got an answer saying more or less the same thing from Michael Malin at CHAM, though a bit more elaborated. I'll try it right away. Olof 

January 29, 2002, 12:49 
Re: Setting PRT as a function of stratification

#4 
Guest
Posts: n/a

This is to record the detailed advice for the benefit of other interested users. The recommended approach for introducing stabilitydependent turbulent viscosity and scalar diffusivity is as follows:
STORE(CMUM,CMUS) in the Q1 file, so as create 3d stores for the parameters CMU and CMU', as these will be a function of local stability parameters. This storage will facilitate GROUND computations and permit printout of these parameters to the RESULT and PHI files. In Group 9 Section 5 of GROUND compute ENUT for the current IZ slab from CMUM*(KE**0.5)*EL1 (note that ENUT is NOT computed for each INDVAR, but rather once for use in all transport equations). Note that the user will have to create the do loop and arrange Farray access to ENUT, CMUM, EL1, etc. For those scalars that require variable turbulent Prandtl number, set PRT(PHI)=1.E10 so as to eliminate the turbulent contribution to the effective exchange coefficient for scalar diffusive transport. Then set PRNDTL(PHI)=GRND so as to inform EARTH that you will be setting the laminar diffusivity for that PHI in Group 9 Section 7 of GROUND. However, in practice you will be setting the effective diffusivity (laminar+turbulent) because we have already suppressed the turbulent contribution by setting PRT(PHI)=1.E10. This procedure is necessary because there is no facility which permits the user to set a variable PRT directly in GROUND. In Group 9 Section 7 of GROUND compute for the current IZ slab the effective diffusivity LAMPR for the diffusive transport of PHI, i.e. (CMUS*(KE**0.5)*EL1 + ENUL/laminar Prandtl number for PHI).Note that again the user will have to create the do loop and arrange Farray access to CMUS, CMUM, EL1, ENUL etc. The use of variable CMU in the ke model means that you will have to be very careful what sourceterm linearisation (KELIN) you employ for the builtin ke model. This is because CMU appears explicity in the ke source terms for some KELIN options (see attached documentation of KELIN). I recommend that the usen of KELIN=3 because this will avoid the problem. The last area to consider is the buoyant production/destruction of k and e. You can modify do loop 30 as required in subroutine gxkegb of the grex source file gxgenk. This is the builtin ke buoyancy term option, which uses a linearisation procedure that prevents negative values of k and e in the event of strong stable stratification. If you use your own coding with explicit source terms, there is a possibility of numerical instability under these conditions. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Gas stratification  Zoltan Turzo  CFX  2  October 8, 2008 01:55 
Multicomponent gas stratification  Zoltan Turzo  CFX  0  October 7, 2008 03:36 
Help setting up pressure function across opening  Jenny  CFX  3  October 18, 2007 00:16 
Setting Prt as a function of stratification?  varghese  FLUENT  0  October 14, 2004 10:27 
stratification of lpg  Ricardo Zuniga  Main CFD Forum  1  October 28, 2003 13:29 