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

Total time for simulation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 5, 2011, 06:48
Red face Total time for simulation
  #1
New Member
 
Join Date: Sep 2011
Posts: 1
Rep Power: 0
Meerkat is on a distinguished road
Hello,

I am running a hydrodyamics simulation that is using a 3rd order TVD Runge-Kutta time stepping method along with appropriate 3rd order reconstruction routines to calculate primitive variables at cell boundaries for use in a variety of Riemann solvers to calculate the fluxes. However, my question is quite simple:

After each completed loop of the three steps in the RK algorithm, is the total amount of time elapsed increased by a single time step, three time steps or a third of the time step?

I have tried to compare my results against test problems in the literature and while I can reproduce the solution it is usually at the wrong time. I am just not sure how to control the time step in the calculation of the overall time elapsed. Everything else is working fine!

Thoughts anyone?
Meerkat is offline   Reply With Quote

Old   September 5, 2011, 07:17
Default
  #2
Senior Member
 
cfdnewbie
Join Date: Mar 2010
Posts: 551
Rep Power: 11
cfdnewbie is on a distinguished road
Hello Meerkat, this link might help you:
http://en.wikipedia.org/wiki/List_of...d-order_method

note that the step (timestep) here is h, and the coefficients are the b's....So in the three step method, your bs are 1/6, 2/3, 1/6 (note they add up to 1!). This means that if you have applied your method correctly, the solution is advanced by h (in space or time).

Might I suggest that you turn to the explicit Euler: du/dt = (u(n+1)-u(n))/delta t... that's the easiest and most simple, I'm quite sure there won't be stability issues with an O(3) method....

So pick the Euler, and try a simple testcase, i.e. a density bubble transported by a constant flow velocity.... thereby, you know the physical speed of the problem and can debug your code!

cheers!
cfdnewbie is offline   Reply With Quote

Reply

Tags
time integrators

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
Hydrostatic Pressure and Gravity miliante OpenFOAM Running, Solving & CFD 132 October 7, 2012 22:50
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03
VOF özgür FLUENT 8 January 6, 2004 09:23
Physical Reason for stability of Implicit Schemes? radhakrishnan Main CFD Forum 25 December 7, 2000 16:25
unsteady calcs in FLUENT Sanjay Padhiar Main CFD Forum 1 March 31, 1999 12:32


All times are GMT -4. The time now is 03:09.