equilibrium wall function

jyothishkumar · July 5, 2011, 05:17

Dear Friends,

Can anyone help me in the expression given in phoenics code for equilibrium wall function.

What is given in the Phoenics documentation is

sturb = (k/LN(E*Re*SQRT(sturb))**2 (brackets are not actually closed)

Is it (k/LN(E*Re*SQRT(sturb)))**2 (sturb^2)

or

sturb = k/LN(E*Re*SQRT(sturb))**2 (is it LN(E*Re*SQRT(sturb))**2 ?)

Someone Please help me in this regard.

Advance thanks

Jyothish

phoenics_cfd · August 3, 2011, 12:16

The following expression can be derived

sturb=kappa/[ ln (E*Re*sqrt(sturb) ]

by re-arranging the log law

u+ = ln(E*y+)/kappa

and noting that

sturb=1/(u+)^2

Re = (u+)*(y+)

July 5, 2011, 05:17	equilibrium wall function	#1
jyothishkumar Member jk Join Date: Jun 2009 Posts: 43 Rep Power: 4	Dear Friends, Can anyone help me in the expression given in phoenics code for equilibrium wall function. What is given in the Phoenics documentation is sturb = (k/LN(*EReSQRT(sturb))2 (brackets are not actually closed)* Is it (k/LN(*EReSQRT(sturb)))2 (sturb^2)* or sturb = k/LN(*EReSQRT(sturb))2 (is it LN(EReSQRT(sturb))2 ?)* Someone Please help me in this regard. Advance thanks Jyothish

Thread Tools
Show Printable Version Email this Page
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
channelFoam for a 3D pipe	AlmostSurelyRob	OpenFOAM	3	June 24, 2011 13:06
Non equilibrium wall function : pressure acquisition problem	MafioTia	OpenFOAM	1	October 10, 2010 04:45
BlockMesh FOAM warning	gaottino	OpenFOAM Native Meshers: blockMesh	7	July 19, 2010 14:11
Problem with rhoSimpleFoam	matteo_gautero	OpenFOAM Running, Solving & CFD	0	February 28, 2008 06:51
Question about wall function	Gary	Main CFD Forum	1	December 3, 2007 10:54

August 3, 2011, 12:16		#2
phoenics_cfd New Member michael malin Join Date: Mar 2009 Location: london Posts: 21 Rep Power: 6	The following expression can be derived sturb=kappa/[ ln (EResqrt(sturb) ] by re-arranging the log law u+ = ln(Ey+)/kappa and noting that sturb=1/(u+)^2 Re = (u+)(y+)