Total time for simulation

Meerkat · September 5, 2011, 07:48

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?

cfdnewbie · September 5, 2011, 08:17

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!

September 5, 2011, 07:48	Total time for simulation	#1
Meerkat New Member Join Date: Sep 2011 Posts: 1 Rep Power: 0	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?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Physical Reason for stability of Implicit Schemes?	radhakrishnan	Main CFD Forum	26	October 3, 2023 23:05
Hydrostatic Pressure and Gravity	miliante	OpenFOAM Running, Solving & CFD	132	October 7, 2012 23:50
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3	bookie56	OpenFOAM Installation	8	August 13, 2011 05:03
VOF	özgür	FLUENT	8	January 6, 2004 09:23
unsteady calcs in FLUENT	Sanjay Padhiar	Main CFD Forum	1	March 31, 1999 13:32

September 5, 2011, 08:17		#2
cfdnewbie Senior Member cfdnewbie Join Date: Mar 2010 Posts: 557 Rep Power: 20	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!