PLUTO
4.4-patch2
|
Runge-Kutta-Legendre driver for integration of diffusion terms. More...
#include "pluto.h"
Functions | |
void | RKL (const Data *d, double dt, timeStep *Dts, Grid *grid) |
Take one step in the solution of the diffusion equation using Runge-Kutta-Legendre (RKL) method. Equations are evolved in time by taking
s
sub-steps including only the diffusion (parabolic) terms. The full step tau
equals the current integration time step (from hydro or a fixed one).
This driver can be used with a 1st or 2nd order algorithm (default is 2).
For the 1st order method (RKL_ORDER
= 1):
where
For the 2nd order method (RKL_ORDER
= 2):
where
References