# How to deal with source term when using RKschemes?

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

 May 10, 2006, 07:49 How to deal with source term when using RKschemes? #1 leaf Guest   Posts: n/a Sponsored Links There are a lot of numerical schemes to deal with the derivatives in the NS equation. My question is how to deal with the chemical raction source term? For example: in Rounge-Kutta schemes there are several intermediate 'times steps' between n to n+1 i.e. v(n+1)=a*v1(n+1/3)+b*v2(n+1/2)+c*v(n) my question is how to deal with the energy and species source terms between those small "fraction time steps" such as from "n to n+1/3" or form "n+1/3 to n+1/2", let them to be constant or let them advance with the dependent variabes such as temperatures and so on?

 May 10, 2006, 09:46 Re: How to deal with source term when using RKsche #2 Jeff Moder Guest   Posts: n/a If you are using the R-K scheme to march in "pseudo-time" to a steady-state solution, then you should try updating source terms (and possibly viscous fluxes) on only a "few" stages (maybe first and last, for example) to save on CPU time for your simulation. You can experiment yourself with updating source terms every R-K stage, versus only on some of the R-K stages. This is a very typical procedure and since you are just using the R-K procedure (in this case) as an iteration technique to find the steady-state solution, you can play around with what gets updated on each stage -- all you care about is that the simulation converges and converges in some "reasonable" number of iterations. If you are using the R-K scheme to march explicity in physical time (ie, a time-accurate simulation), then you should update your source terms every R-K stage. Of course, you could evaluate the loss of accuracy in not doing so and make your own choice in accuracy loss versus CPU time reduction for a time-accurate simulation using R-K to march in physical time.

 May 11, 2006, 10:34 Re: How to deal with source term when using RKsche #3 Harish Guest   Posts: n/a Though the intermediate t* are not clearly known evaluating the source term at all the intermediate steps using the same t seems to work for me. -H

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post zwdi FLUENT 14 June 27, 2017 15:40 GGerber OpenFOAM Installation 54 April 24, 2015 16:02 hjasak OpenFOAM Post-Processing 69 April 24, 2008 01:24 jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51 Greg Perkins FLUENT 8 October 20, 2000 12:40