CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Non equilibrium wall function : pressure acquisition problem

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 6, 2010, 18:52
Question Non equilibrium wall function : pressure acquisition problem
  #1
New Member
 
Mattia Maffioletti
Join Date: Oct 2009
Posts: 1
Rep Power: 0
MafioTia is on a distinguished road
Hi everyone.

I'm trying to implement a fluent-type wall function sensitized to pressure gradient, i.e. "non-equilibrium WF". There is no problem in epsilon/omega WF when I get the pressure gradient, needed for my calculation, in this way:

const volScalarField& P = db().lookupObject<volScalarField>("p");
volVectorField grdntP = fvc::grad(P);

If I apply the same procedure in nutSpalartAllmarasWF ( that follows the general Spalding law-of-the-wall but it has to be sentivized to pressure gradient, too!) the solution explode in an unexplainable way. I've already tried to debug my code on every line that I added and I've concluded that the problematic line is the second I wrote before.
My idea is that this procedure applied to nut-WF result in an unintentional modify of the pressure field, but I thin it's a weak idea, if I had understood OF works.

Any idea?
MafioTia is offline   Reply With Quote

Old   October 10, 2010, 04:45
Default
  #2
New Member
 
Maolong LIU
Join Date: Apr 2010
Location: Japan
Posts: 28
Rep Power: 7
maolongliu is on a distinguished road
Hi, have you solved this problem?
I am now also trying to write a new wall function including pressure gradient.

Quote:
Originally Posted by MafioTia View Post
Hi everyone.

I'm trying to implement a fluent-type wall function sensitized to pressure gradient, i.e. "non-equilibrium WF". There is no problem in epsilon/omega WF when I get the pressure gradient, needed for my calculation, in this way:

const volScalarField& P = db().lookupObject<volScalarField>("p");
volVectorField grdntP = fvc::grad(P);

If I apply the same procedure in nutSpalartAllmarasWF ( that follows the general Spalding law-of-the-wall but it has to be sentivized to pressure gradient, too!) the solution explode in an unexplainable way. I've already tried to debug my code on every line that I added and I've concluded that the problematic line is the second I wrote before.
My idea is that this procedure applied to nut-WF result in an unintentional modify of the pressure field, but I thin it's a weak idea, if I had understood OF works.

Any idea?
maolongliu is offline   Reply With Quote

Reply

Tags
pressure gradient, wall functions

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Fluent3DMeshToFoam simvun OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 48 May 14, 2012 05:20
Compilation errors in ThirdPartymallochoard feng_w OpenFOAM Installation 1 January 25, 2009 07:59
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51
a problem in calculating pressure drop in Fluent? yu chun FLUENT 1 May 18, 2004 03:40
X-Y plot of Yplus in Fluent 5.3 Luo FLUENT 24 April 11, 2000 06:07


All times are GMT -4. The time now is 14:15.