equilibrium wall function

jyothishkumar · July 5, 2011, 06: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, 13: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, 06:17	equilibrium wall function	#1
jyothishkumar Member jk Join Date: Jun 2009 Posts: 64 Rep Power: 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(*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	Search this Thread
Show Printable Version Email this Page	Search this Thread: Advanced Search
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 14:06
Non equilibrium wall function : pressure acquisition problem	MafioTia	OpenFOAM	1	October 10, 2010 05:45
[blockMesh] BlockMesh FOAM warning	gaottino	OpenFOAM Meshing & Mesh Conversion	7	July 19, 2010 15:11
Problem with rhoSimpleFoam	matteo_gautero	OpenFOAM Running, Solving & CFD	0	February 28, 2008 07:51
Question about wall function	Gary	Main CFD Forum	1	December 3, 2007 11:54

August 3, 2011, 13:16		#2
phoenics_cfd Member michael malin Join Date: Mar 2009 Location: london Posts: 35 Rep Power: 18	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+)