CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Local stepping for steady state

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

LinkBack Thread Tools Display Modes
Old   June 16, 2003, 17:57
Default Local stepping for steady state
Posts: n/a
I am trying to compte a 2-D MHD steady state solution using an upwind scheme and a multigrid space domain.

I don't really understand how to compte the steady state with the multi-stage runge-kutta scheme using a local time step for each cell.

When we compute the solution for a cell at a local time t we use its neighbors to evaluate the cell residual. But the parameters of these neighbors cells are obtained for another local times different of the local time t (since we are using diffrent time stepping in function of the volume of each cell). I tried to do a linear extrapolation of these parameters from the local times of their cells to the local time t (using the last known residual for each neighbor cell). This method didn't work because the residual is not linear.

Is there any other method to do that ?

Thanks Mehdi
  Reply With Quote

Old   June 17, 2003, 05:00
Default Re: Local stepping for steady state
Posts: n/a
I am not very familiar with MHD. But it seems that the problem is more related to numerics.

A Runge-Kutta scheme is an explicit scheme. For stability reasons, the time-step has to be limited by a local or global criteria.

Computing a steady-state solution is equivalent to compute the solution at large physical times. Using a local time-step for each cell is equivalent to take the larger numerically-stable time-step for each cell, thus "bringing cell as close as possible to the steady-state independantly of the others".

In theory, you don't have to extrapolate in time the variables in the neighbouring cells, eg construct the RHS as if all cells were at the same physical time. This should work since you are only interested by the steady-state solution, when the variables are not evolving anymore.

For pure fluid dynamic, implicit methods are prefered since larger "time-steps" can be used and steady-state solution is reached quicker.

  Reply With Quote

Old   June 17, 2003, 09:43
Default Re: Local stepping for steady state
Posts: n/a
Thank you for you answer JDC, I think that you understood perfectly my problem.

I tried first you method : computing the Residual using the parameters of the neighbors cell without taking into account their local time. This method gave me a lot of problem especially during the reconstruction process (negative pressure or negative mass denstities).

I will try it again and see what can happen and let you know if it works.

Thank you anyway for you help JDC
  Reply With Quote


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
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03
Constant velocity of the material Sas CFX 15 July 13, 2010 08:56
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07

All times are GMT -4. The time now is 05:38.