CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

Cartesian grid specific solverIterationpresure solver

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

LinkBack Thread Tools Search this Thread Display Modes
Old   May 26, 2006, 06:04
Default Hello, I would like to know
anne dejoan
Join Date: Mar 2009
Location: madrid, spain
Posts: 66
Rep Power: 17
anne is on a distinguished road

I would like to know if there is any flag option
in OpenFoam that could specify that the
solver has to be applied to a cartesian grid.

I asked this because (I think ?) it could
avoid some CPU time when running on a cartesian grid (it would avoid all the surfave normal vectors computations and so on).

Another question relates to the number
of iterations required to solve the
pressure. In a developped channel flow, when using
the IGCCg solver it requires about 30 iterations for each pressure resolution and when using AMG it requires about 10-15. Is there aby way to
optimize the iteration number ?

Finally, I would like to know if OpenFoam makes use of Rhie-Chow interpolation.

Thank You for your help,

anne is offline   Reply With Quote

Old   May 26, 2006, 06:21
Default Cartesian Meshes You are ri
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,905
Rep Power: 33
hjasak will become famous soon enough
Cartesian Meshes

You are right, on Cartesian meshes you could use special tricks to make the discretisation run faster by using the properties of the mesh. In order to do that, some methods (e.g. gradient) would have to be written again specifically for Cartesian meshes. This also leads to the issues with signalling (you would need 2 mesh classes) etc. and Cartesian meshes are not too useful. On balance, I would say this is not worth it - unless you wish to do things that are only possible on a Cartesian mesh (e.g. sixth order gradients). For this reason, OpenFOAM treats all meshes as unstructured polyhedral.

The code has a natural slot for this kind of development: you could copy the finiteVolume library and re-implement the discretisation for the new mesh type. All other stuff, from pre-processing to matrices and fields can be re-used as appropriate. An example of something simular would be the finiteArea library, solving FVM on a curved moving 2-D surface mesh.


Sounds like for your case ICCG works much better - just pick that. This is why we've got several solvers in the first place :-)

Rhie-Chow Interpolation

Yes, of course, but not "by the book". Have a look at my Thesis or try to catch Rasmus Hemph - he has recently written a note about this.


Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting:
hjasak is offline   Reply With Quote

Old   May 26, 2006, 09:59
Default Dear Hrv, Thank you for you
anne dejoan
Join Date: Mar 2009
Location: madrid, spain
Posts: 66
Rep Power: 17
anne is on a distinguished road
Dear Hrv,

Thank you for your fast answer.

In our project we would like to perform LES for the study of contaminant dispersion in a urban context.
For this we are thinking in using OpenFoam because
it is an open source code.
I have just started with DNS and LES of fully channel flow for a first go (I don't know anything about c++). The second step I think about is to perform a surface mounted cube in a channel. As it is still
cartesian I was wondering about the possibility of
switch off the unstructuired-related part to gain
in CPU. But in viw of my poor knowledge
in c++ I am not sure I will re-write
the discrete operators.

However, could you let me know where I can access your thesis ? It would be very helpfull for me to have it.

I have two more quetions ....

The first one is on turbulent inlet conditions issued from a precursor channel flow computation
that I had posted. I still not have any answer,
I know that you must be very busy but
some indications on how to do it it
would be great.

The second is: is there any good tutorial to
start with tracking particules and scalar passive

Thanks again,

anne is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

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

Similar Threads
Thread Thread Starter Forum Replies Last Post
cartesian vs cylindrical solver giles Main CFD Forum 0 August 26, 2008 19:00
Cartesian grid generation bobbyzhu Main CFD Forum 0 January 10, 2008 22:33
Cartesian Grid Generation FASA Main CFD Forum 1 June 19, 2006 05:13
Cartesian grid apple Main CFD Forum 0 April 9, 2006 09:09
Cartesian Grid Harish FLUENT 0 January 3, 2006 00:09

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