AStableDirk4

Fourth-order diagonally implicit Runge Kutta method (Dirk) with three stages plus an update.

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

var element = document.getElementById("moose-equation-dcc46439-2901-4405-8fa0-31d8ed028621");katex.render("\\begin{array}{c|ccc} \\gamma & \\gamma & 0 & 0 \\\\ 1/2 & 1/2 - \\gamma & \\gamma & 0 \\\\ 1 - \\gamma & 2 \\gamma & 1 - 4 \\gamma & \\gamma \\\\ \\hline & \\dfrac{1}{24(1/2-\\gamma)} & 1 - \\dfrac{1}{12(1/2-\\gamma)^2} & \\dfrac{1}{24(1/2-\\gamma)^2} \\end{array}", element, {displayMode:true,throwOnError:false});

where $var element = document.getElementById("moose-equation-ad62d4eb-9060-4e06-b670-987044d5c8f2");katex.render("\\gamma = 1/2 + \\dfrac{\\sqrt{3}}{3} \\cos( \\dfrac{\\pi}{18}) \\approx 1.06857902130162881", element, {displayMode:false,throwOnError:false});$

The stability function for this method is:

var element = document.getElementById("moose-equation-749646c4-5df4-4eae-9791-6c08daea39c7");katex.render("R(z) = \\dfrac{-0.76921266689461 z^3 - 0.719846311954034 z^2 + 2.20573706179581 z - 1.0}{ 1.22016884572716 z^3 - 3.42558337748497 z^2 + 3.20573706551682 z - 1.0}", element, {displayMode:true,throwOnError:false});

The method is not L-stable; it is only A-stable:

var element = document.getElementById("moose-equation-ff9c0e1c-48a8-4edb-a921-f744ed1ec13b");katex.render("\\lim_{z->\\infty} R(z) = -0.630414937726258", element, {displayMode:true,throwOnError:false});

Notes

Method is originally due to Crouzeix (1975)
Since $var element = document.getElementById("moose-equation-b28a3bef-ed0f-42fd-bef6-8fa842d0747c");katex.render("\\gamma", element, {displayMode:false,throwOnError:false});$ is slightly larger than $var element = document.getElementById("moose-equation-17be7884-d5e6-4c07-a604-c487e7bb2945");katex.render("1", element, {displayMode:false,throwOnError:false});$ , the first stage involves evaluation of the non-time residuals for $var element = document.getElementById("moose-equation-4a7d69dc-8bc4-41a7-93fb-d03471528d15");katex.render("t > t_n+\\Delta t", element, {displayMode:false,throwOnError:false});$ , while the third stage involves evaluation of the non-time residual for $var element = document.getElementById("moose-equation-626ec226-173f-4192-9e0e-b57fbec49d37");katex.render("t < t_n", element, {displayMode:false,throwOnError:false});$ , which may present an issue for the first timestep (if e.g. material properties or forcing functions are not defined for $var element = document.getElementById("moose-equation-17a120a1-4f48-4bd8-b252-1a3809707f17");katex.render("t<0", element, {displayMode:false,throwOnError:false});$ . We could handle this by using an alternate (more expensive) method in the first timestep, or by using a lower-order method for the first timestep and then switching to this method for subsequent timesteps.

Input Parameters

safe_startTrueIf true, use LStableDirk4 to bootstrap this method.
Default:True
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:If true, use LStableDirk4 to bootstrap this method.
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

References

M. Crouzeix. Sur l'approximation des equations differentielles operationelles lineaires par des methodes de Runge Kutta. PhD thesis, Universite Paris VI, Paris, 1975.

@phdthesis{crouzeix1975,
    author = "Crouzeix, M.",
    title = "Sur l'approximation des equations differentielles operationelles lineaires par des methodes de Runge Kutta",
    school = "Universite Paris VI, Paris",
    year = "1975"
}

Notes
Input Parameters
References