CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Different Meshes for Different Physics

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

LinkBack Thread Tools Search this Thread Display Modes
Old   August 26, 2018, 17:59
Default Different Meshes for Different Physics
New Member
physicsunited's Avatar
Join Date: Feb 2018
Posts: 10
Rep Power: 4
physicsunited is on a distinguished road
Dear fellow cfd-online members,

I am using a modified icoFoam setup, that couples transient passive scalar transport to the momentum equations.

My issue is that my momentum equations don't require a fine mesh, while my scalar transport problem requires a relatively very fine mesh (due to high Peclet numbers). This further constrains my time step choice (due to the Courant number and viscous stability constraints) and adds a significant number of time steps to my problem. These time steps do not add significant physical information to my problem, but are rather there due to numerical stability requirements. Basically in one of the equations I have to use very fine mesh while the other one works OK with a coarser mesh.

I was wondering whether it would be possible to solve my coupled system on different meshes (one for each set of equations), and project the velocity results from the coarse mesh to the finer mesh to be used in solving my scalar transport equation.

Is there anyway that I can assign different meshes for each of these equations and thus save computational time?

physicsunited is offline   Reply With Quote

Old   August 27, 2018, 01:40
Senior Member
piu58's Avatar
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 590
Rep Power: 8
piu58 is on a distinguished road
You may always use a script which runs both simulations, one after the other. You need a tiny program in the programming language of your choice, which reads the velocity results of the coarse mesh and transfers it (with interpolation) to the finer one.
Uwe Pilz
Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950)
piu58 is offline   Reply With Quote

Old   August 27, 2018, 02:43
Senior Member
Peter Baskovich
Join Date: Jul 2014
Posts: 120
Rep Power: 7
pete20r2 is on a distinguished road
I'd be curious to see how you go, I suspect you will spend more time interpolating and IO operations than just solving on the fine mesh. Unless you do it inside the one solver with two meshes and keep the fields in memory while you solve I don't think you'll see the benefit.
pete20r2 is offline   Reply With Quote

Old   August 27, 2018, 05:42
Senior Member
Gerhard Holzinger
Join Date: Feb 2012
Location: Austria
Posts: 267
Rep Power: 23
GerhardHolzinger will become famous soon enoughGerhardHolzinger will become famous soon enough
It can be done, see this post of mine, which contains a screenshot.

I took inspiration from the multi-region solvers, and the mapFields utility.

I have no idea, when the additional effort (the second mesh, the mapping of fields each time step) is worth it.
GerhardHolzinger is offline   Reply With Quote


coupled equations, multiphysics, multiple meshes

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
how to set periodic boundary conditions Ganesh FLUENT 14 November 26, 2018 12:26
Different Physics Contnuum for Same Region itsanisant STAR-CCM+ 0 October 10, 2014 01:30
Hex and Tet meshes - simplefoam comparison danvica OpenFOAM Running, Solving & CFD 10 January 4, 2013 02:18
Multiple Physics Continua boathead STAR-CCM+ 6 December 4, 2011 13:33
physics phD to CFD? jck87 Main CFD Forum 13 October 9, 2010 17:15

All times are GMT -4. The time now is 02:33.