I need to have a look to LES with wall functions. Is there anyone that already done some work about?. I know that in Chatou someone already used LES with WF (if I'm right was the topic of David phd, and also Thomas developped some dynamic WF) . Is it possible to have the subroutines and a quick explanation how to use. If something exist under which version is implemented? Thanks very much for the help.

-- StefanoRolfo - 05 Jun 2008

Hi Stefano, by default, if the centre of your wall cell is further from the wall that y+=10,88, Code_Saturne will automatically use a logarithmic wall function (default), or a power-law wall function (with ILOGPO=0). It's been used quite often in Chatou on industrial applications where we do not have a choice (mesh size limitations) with sometimes good results ... but not always! Indeed, both David and Thomas have developped alternative more precise (but more specific) wall functions. I'll let them comment on their specificities.

-- MarcSakiz - 10 Jun 2008

You can also try using the Reichardt function which goes smoothly between the laminar and logarithmic regions.

I have a test version: causta.F where you should set YPLULI(IPHAS) = 1D-6 in the usini1.F (in SUBROUTINE USIPSU) so that the Reichardt function is called all the time.

-- RichardHoward - 25 Jun 2009