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

How to deal with source term when using RKschemes?

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

LinkBack Thread Tools Display Modes
Old   May 10, 2006, 07:49
Default How to deal with source term when using RKschemes?
Posts: n/a
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?
  Reply With Quote

Old   May 10, 2006, 09:46
Default Re: How to deal with source term when using RKsche
Jeff Moder
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.

  Reply With Quote

Old   May 11, 2006, 10:34
Default Re: How to deal with source term when using RKsche
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.

  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
swak4foam building problem GGerber OpenFOAM Installation 54 April 24, 2015 16:02
momentum source term zwdi FLUENT 13 December 5, 2013 18:58
DxFoam reader update hjasak OpenFOAM Post-Processing 69 April 24, 2008 01:24
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51
UDF Scalar Code: HT 1 Greg Perkins FLUENT 8 October 20, 2000 12:40

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