# Problem with wall function

 Register Blogs Members List Search Today's Posts Mark Forums Read

August 23, 2011, 03:01
Problem with wall function
#1
Member

Gregoire Junqua
Join Date: Jun 2011
Location: China
Posts: 58
Blog Entries: 1
Rep Power: 7
Dear Openfoamer.

I am a new in modelling, I want to study turbulence over a plane plate, but now
I am dealing a problem with the wall function.

Normally i should have a beautiful log law function but I got the following wind profile showed in the following picture.

My boundary are as follow

Inlet
pressure zerogradient U fixed value k fixedvalue epsilon fixedvalue

Outlet

Wall
pressure zerogradient U fixedvalue k kqrwallfunction epsilon epsilonwallfunctionn

The case consist of a box
inlet is in the left
outlet is on the top and on the right
wall is on the bottom
the case is in 2D.

If someone could give me an advice it would be nice
Attached Images
 windProfile.jpg (14.5 KB, 45 views)

 August 23, 2011, 15:11 #2 Member   Join Date: Mar 2010 Posts: 42 Rep Power: 8 Are the front and back walls with the "empty" boundary condition? And what is the size of your mesh? Can you post the values of y+ (i think the command is yPlusRas. Standart k-epsilon turbulence mode? What about the convergence? Did it pass the 1e-5?

August 23, 2011, 21:15
#3
Member

Gregoire Junqua
Join Date: Jun 2011
Location: China
Posts: 58
Blog Entries: 1
Rep Power: 7
Hi
The front and back got empty conditions. and yes i use a k-epsilon model
The profile of yplus is in the atached file.
I had done several test with several model, it show that it's the boundary conditions are the problem
I tried Pisofoam, simplefoam, with the both condition LES RAS
but the given answer is the same.
For the convergeance
i don't knbow how to get it. I give you a sample of a time step perhaps it s inside.
Quote:
 Time = 10.9215 Courant Number mean: 0.729365 max: 1.54273 DILUPBiCG: Solving for Ux, Initial residual = 3.69415e-05, Final residual = 6.31557e-10, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.000546025, Final residual = 9.49463e-09, No Iterations 1 DICPCG: Solving for p, Initial residual = 5.72492e-05, Final residual = 4.1219e-06, No Iterations 2 time step continuity errors : sum local = 7.61108e-13, global = 3.79935e-14, cumulative = 1.00392e-07 DICPCG: Solving for p, Initial residual = 2.02581e-05, Final residual = 9.95151e-07, No Iterations 56 time step continuity errors : sum local = 1.83755e-13, global = -5.35126e-15, cumulative = 1.00392e-07 DILUPBiCG: Solving for k, Initial residual = 9.97577e-06, Final residual = 9.97577e-06, No Iterations 0
I still searching why i got that problem
Attached Images
 windprofileyplus.png (6.7 KB, 27 views)

 August 24, 2011, 07:24 #4 Member   Join Date: Mar 2010 Posts: 42 Rep Power: 8 Hi, is the left y axis the y+ value? Try increasing the size of the mesh near the walls and get a value of y+ between 30 and 300, your mesh seems to small for the k-epsilon model.

 August 24, 2011, 20:57 #5 Member   Gregoire Junqua Join Date: Jun 2011 Location: China Posts: 58 Blog Entries: 1 Rep Power: 7 Hi, No the wind is on the left the input is an uniform 10m/s wind. on the right is the yplus. I will today increase the mesh size near the wall... I hope it would work. Thanks for the advice

 August 25, 2011, 23:32 #6 Member   Gregoire Junqua Join Date: Jun 2011 Location: China Posts: 58 Blog Entries: 1 Rep Power: 7 Hi, I found that it the wall function or some misunderstood of this experiment because i tried different mesh size, I try the dnsFoam And all give the same answer. So i think the configuration I gave is wrong If somebody know how to have a beautiful loglaw please tell me... perhaps the input should be a loglaw too?

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Yr0gErG FLUENT 3 June 12, 2013 02:12 Attesz CFX 7 January 5, 2013 04:32 alxasli FLUENT 3 January 10, 2011 07:22 asherah STAR-CCM+ 0 February 19, 2010 18:45 hills1 CFX 2 October 12, 2009 05:36

All times are GMT -4. The time now is 10:56.