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

@Developers: External RK4 time stepping

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 4, 2014, 11:24
Default @Developers: External RK4 time stepping
  #1
New Member
 
nilesh
Join Date: Mar 2014
Location: Kanpur / Mumbai, India
Posts: 27
Rep Power: 4
nilesh is on a distinguished road
Hi,
I need to implement RK4 for the external time stepping (i.e. instead of first order dual time stepping).
From what I understood, line no 6125 of solver_direct_mean.cpp has to be modified (I do not mind entering the coefficients directly as constants at this point of time instead of taking them from the config file).
However, I do not understand how to go about finding the solution (space derivatives) at the intermediate stages of the Runge kutta iteration. Please help me with this.
Thanks.
nilesh is offline   Reply With Quote

Old   June 4, 2014, 16:44
Default
  #2
New Member
 
nilesh
Join Date: Mar 2014
Location: Kanpur / Mumbai, India
Posts: 27
Rep Power: 4
nilesh is on a distinguished road
Sorry, apparently that line is only for calculating the residue. I am confused by the object oriented structure of the code. Which is the function to be called for calculating the space derivative part? Does ->GetSolution() which returns SOLUTION only return the solution or also calculated it?
I am very eager (desperate) to use RK3 (or RK4) for the external loop. If any user is interested, we can collaborate and work on it.

Thanks.
nilesh is offline   Reply With Quote

Old   June 5, 2014, 15:22
Default
  #3
New Member
 
David Manosalvas
Join Date: Feb 2014
Posts: 24
Rep Power: 4
demanosalvas is on a distinguished road
Nilesh,

RK time integration is already implemented in the code, all you need to do is select this option in your config file.

David
demanosalvas is offline   Reply With Quote

Old   June 7, 2014, 10:22
Default
  #4
New Member
 
nilesh
Join Date: Mar 2014
Location: Kanpur / Mumbai, India
Posts: 27
Rep Power: 4
nilesh is on a distinguished road
Quote:
Originally Posted by demanosalvas View Post
Nilesh,

RK time integration is already implemented in the code, all you need to do is select this option in your config file.

David
That is for internal iterations, I wanted RK for external iterations when using dual-time stepping.
nilesh is offline   Reply With Quote

Old   June 13, 2014, 13:55
Default
  #5
New Member
 
David Manosalvas
Join Date: Feb 2014
Posts: 24
Rep Power: 4
demanosalvas is on a distinguished road
Nilesh,

Can you elaborate a bit more on why would you like to do that?

David
demanosalvas is offline   Reply With Quote

Old   June 13, 2014, 14:39
Default
  #6
New Member
 
nilesh
Join Date: Mar 2014
Location: Kanpur / Mumbai, India
Posts: 27
Rep Power: 4
nilesh is on a distinguished road
I wanted to use RK3 for the external time marching which is presently being done by euler explicit or 3-time level BDF method. I have done it for the the SU2_CFD module as of now (without support for mesh deformation or turbulence models) and the results look great.
Would be nice if you include it throughout the SU2 suite.

Thanks.
nilesh is offline   Reply With Quote

Old   June 24, 2014, 16:53
Default
  #7
Super Moderator
 
Thomas D. Economon
Join Date: Jan 2013
Location: Stanford, CA
Posts: 270
Rep Power: 6
economon is on a distinguished road
Hi,

Just to make sure everything is clear.. you mean that you would like to perform time-accurate simulations with an explicit RK scheme rather than the current dual-time stepping method (1st- or 2nd-order BDF), right?

We do have a "time-stepping" mode for unsteady calculations (check the config file options), which will essentially use a global time step for the entire mesh and solve the flow in a time-accurate manner with the chosen scheme. You might try experimenting with this to see if it meets your needs. Otherwise, perhaps we can chat a little more about your implementation.

All the best,
Tom
economon is offline   Reply With Quote

Old   June 25, 2014, 02:58
Default
  #8
New Member
 
nilesh
Join Date: Mar 2014
Location: Kanpur / Mumbai, India
Posts: 27
Rep Power: 4
nilesh is on a distinguished road
Quote:
Originally Posted by economon View Post
Hi,

Just to make sure everything is clear.. you mean that you would like to perform time-accurate simulations with an explicit RK scheme rather than the current dual-time stepping method (1st- or 2nd-order BDF), right?

We do have a "time-stepping" mode for unsteady calculations (check the config file options), which will essentially use a global time step for the entire mesh and solve the flow in a time-accurate manner with the chosen scheme. You might try experimenting with this to see if it meets your needs. Otherwise, perhaps we can chat a little more about your implementation.

All the best,
Tom
Thanks. Last time I tried it, I was unable to write files. Actually the filename of the output files did not increment with the time step and consequently the output files kept overwriting the previous ones. I will give it another try and keep you posted.
nilesh is offline   Reply With Quote

Old   June 25, 2014, 20:30
Default
  #9
Super Moderator
 
Thomas D. Economon
Join Date: Jan 2013
Location: Stanford, CA
Posts: 270
Rep Power: 6
economon is on a distinguished road
I see.. yep, please let us know your experience, and if you are still running into issues, we should be able to modify the code quickly to resolve it.

T
economon is offline   Reply With Quote

Old   March 23, 2015, 05:58
Default Issue with "TIME_STEPPING"
  #10
New Member
 
Deep Ray
Join Date: Jan 2014
Posts: 6
Rep Power: 4
d_ray is on a distinguished road
Quote:
Originally Posted by economon View Post
I see.. yep, please let us know your experience, and if you are still running into issues, we should be able to modify the code quickly to resolve it.

T
I have been trying to run an unsteady problem with RK3 (w/o dual stepping), as mentioned previously by Nilesh. I am facing the following issues:
1) The solution files are not being saved in an incremental manner. The solutions are overwritten into the previous file. However, the incremental storage is working for the dual time stepping scenario.
2) It is unclear how to specify the final time for the simulation using TIME_STEPPING.

I have been trying to look into the cfg files of various unsteady examples. Unfortunately all of them seem to use the dual time stepping strategy.

I am attaching my cfg file for your reference.

Thanks.
Attached Files
File Type: gz forward_hllc.cfg.gz (3.0 KB, 3 views)
d_ray is offline   Reply With Quote

Reply

Tags
dual time stepping method, explicit, rk4, su2

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
Micro Scale Pore, icoFoam gooya_kabir OpenFOAM Running, Solving & CFD 2 November 2, 2013 14:58
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
number of iterations per time step chouki FLUENT 1 August 13, 2013 00:11
pisoFoam with k-epsilon turb blows up - Some questions Heroic OpenFOAM Running, Solving & CFD 26 December 17, 2012 04:34
Transient simulation not converging skabilan OpenFOAM Running, Solving & CFD 12 September 17, 2007 17:48


All times are GMT -4. The time now is 22:40.