Heun

Heun's (aka improved Euler) time integration method.

Heun's (or "improved Euler") time integration method is second-order accurate in time. It is a two-step explicit method and a special case of the 2nd-order Runge-Kutta method.

Description

With UU, the vector of nonlinear variables, and AA, a nonlinear operator, we write the PDE of interest as:

Ut=A(t,U(t))\dfrac{\partial U}{\partial t} = A(t, U(t))

Using t+Δtt+\Delta t for the current time step, and tt for the previous step, Heun's method can be written:

U(t+Δt)=U(t)+Δt2(A(t,U(t))+A(t,U(t)+ΔtA(t,U(t))))U(t+\Delta t) = U(t) + \dfrac{\Delta t}{2} \left(A(t, U(t)) + A(t,U(t) + \Delta t A(t, U(t)) ) \right)

This method can be expressed as a Runge-Kutta method with the following Butcher Tableau:

001101/21/2\begin{array}{c|cc} 0 & 0 \\ 1 & 1 & 0 \\ \hline & 1/2 & 1/2 \end{array}
warningwarning

All kernels except time-(derivative)-kernels should have the parameter implicit=false to use this time integrator.

warningwarning

ExplicitRK2-derived TimeIntegrators (ExplicitMidpoint, Heun, Ralston) and other multistage TimeIntegrators are known not to work with Materials/AuxKernels that accumulate 'state' and should be used with caution.

Input Parameters

  • variablesA subset of the variables that this time integrator should be applied to

    C++ Type:std::vector<VariableName>

    Unit:(no unit assumed)

    Controllable:No

    Description:A subset of the variables that this time integrator should be applied to

Optional Parameters

  • control_tagsAdds user-defined labels for accessing object parameters via control logic.

    C++ Type:std::vector<std::string>

    Unit:(no unit assumed)

    Controllable:No

    Description:Adds user-defined labels for accessing object parameters via control logic.

  • enableTrueSet the enabled status of the MooseObject.

    Default:True

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:No

    Description:Set the enabled status of the MooseObject.

Advanced Parameters