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

Fourth order surface normal gradient interpolation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 16, 2007, 09:33
Default I have been having troubles im
  #1
Member
 
Adam Donaldson
Join Date: Mar 2009
Location: Ottawa, Ontario, Canada
Posts: 37
Rep Power: 8
adona058 is on a distinguished road
I have been having troubles implementing the fourth order surface normal gradient interpolation scheme. Each time I try, I receive an error message stating that only corrected, uncorrected and limited can be chosen.

I have checked the config dictionary files for FoamX, and they appear to contain the Fourth entry. Is there a different location where the solvers look for valid schemes?

Has anyone else experienced this issue, and if so, would you have any suggestions on how I can implement the fourth order scheme.

Thank you,

Adam
adona058 is offline   Reply With Quote

Old   September 6, 2007, 10:16
Default I just stumbled on the same pr
  #2
Member
 
Ola Widlund
Join Date: Mar 2009
Location: Sweden
Posts: 87
Rep Power: 8
olwi is on a distinguished road
I just stumbled on the same problem. It has nothing to do with FoamX, cause I've never used it. Must be a list somewhere...

What version are you using? I'm still in 1.3.

/Ola
olwi is offline   Reply With Quote

Old   September 24, 2007, 07:10
Default Even though it is listed in th
  #3
pbo
Member
 
Patrick Bourdin
Join Date: Mar 2009
Posts: 40
Rep Power: 8
pbo is on a distinguished road
Even though it is listed in the user's manual,
the fourth order surface normal gradient interpolation scheme simply does not exist anymore!!
(as well as the bounded scheme)

have a look at: $FOAM_SRC/finiteVolume/finiteVolume/snGradSchemes

This is for OF 1.3 and above.
However in OF 1.2 if you look there:
$FOAM_SRC/OpenFOAM/finiteVolume/snGradSchemes,
tada! boundedSnGrad and fourthSnGrad are available.

Why did those schemes get removed after OF 1.2?
Hrv, Henry, any comments?

Another question related to surface normal gradient schemes:
the current snGradSchemes (based on a central finite difference between neighbour and owner cell centres) limit the application of OF to orthogonal or midly non-orthogonal meshes. Wouldn't a linear interpolation of the cell gradients at the face provide an estimate as accurate, yet able to tackle severly non-orthogonal meshes (such as the ones with anisotropic tetrahedral elements in the boundary layer of a wall-bounded high Re flow -- see picture below for an example)?


pbo is offline   Reply With Quote

Old   September 24, 2007, 08:19
Default Fourth order snGrad still pres
  #4
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Fourth order snGrad still present and correct in the dev line - you should be able to simply pick up the files and add to the compilation...

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   September 24, 2007, 08:34
Default It's actually present in my OF
  #5
Member
 
Ola Widlund
Join Date: Mar 2009
Location: Sweden
Posts: 87
Rep Power: 8
olwi is on a distinguished road
It's actually present in my OF 1.3 installation (from Hrv) as well, but it's still not a valid runtime selection. Should it appear somewhere else as well?

/Ola
olwi is offline   Reply With Quote

Old   September 24, 2007, 08:40
Default ? Just tested it: Both thin
  #6
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
? Just tested it:

Both things like:

laplacian(nu,U) Gauss linear fourth;

and


snGradSchemes
{
default fourth;

}

work just fine. I have found that it was missing in the Make/files for the finiteVolume library, added it there (checked in the change into SVN) and all is well.

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   September 24, 2007, 09:37
Default Hrv, Is there any reason wh
  #7
pbo
Member
 
Patrick Bourdin
Join Date: Mar 2009
Posts: 40
Rep Power: 8
pbo is on a distinguished road
Hrv,

Is there any reason why surface normal gradients are always approximated by:
{phi_N - phi_P}/|d| + nonOrthoCorrectors.

Why is the approach based on interpolating the cell gradients at the face intersection (like in the skew correction) discarded?
How can one accurately compute the normal gradients in the case depicted in my previous message? The nonOrthogonal correction seems only to work if the non-orthogonality is moderate...
pbo is offline   Reply With Quote

Old   September 24, 2007, 11:55
Default Yes, this is fundamental - if
  #8
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Yes, this is fundamental - if you do not have the phi_N - phi_P/delta bit, the implicit Laplacian will not be formed correctly on the compact molecule. Everything else is treated as a correction, including any non-orthogonal or skewness terms you may want to have. For details, see my Thesis.

Enjoy,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   September 24, 2007, 15:12
Default ok, I got that part. My qu
  #9
pbo
Member
 
Patrick Bourdin
Join Date: Mar 2009
Posts: 40
Rep Power: 8
pbo is on a distinguished road
ok,

I got that part.
My question was in fact about the alternative way to compute the surface normal gradients:
Instead of using {phi_N - phi_P}/|d| + ..., why not evaluating the cell-centred gradient for the 2 cells adjacent to the face at hand (using, say, Gauss linear), and interpolate it to the face. This involves a bigger computational molecule but on the other hand, it does not require any non-orthoganility correction and therefore should be more robust wrt non-orthogonality.
Has it ever been implemented in FOAM?

Patrick
pbo is offline   Reply With Quote

Reply

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
Fourth Order Runge Kutta time integration lr103476 OpenFOAM Running, Solving & CFD 10 November 11, 2014 09:27
How can I get the pressure gradient along the direction normal to the local wall surface qtian OpenFOAM Running, Solving & CFD 3 April 4, 2014 05:05
Help for fourth order accurate convective schemes Z.C.Wang Main CFD Forum 0 January 15, 2009 07:53
4thorder surface normal gradient interpolation smehdi609 OpenFOAM Running, Solving & CFD 0 June 13, 2008 19:52
surface tension as function of normal gradient tom FLUENT 0 April 30, 2008 14:47


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