| NonStationaryProblem | nsteps(in) # |
| [renumber(in) #] | |
| deltaT(rn) # | |
| alpha(rn) # | |
| [lumpedcapa() #] | |
| [sparselinsolverparams() #] | |
| [exportfields(ia) #] | |
| [changingProblemSize() #] |
Linear implicit
integration scheme for transient transport problems. The generalized
midpoint rule (sometimes called
-method) is used for time discretization, with alpha parameter,
which has limits
. For
explicit Euler
forward method is obtained, for
implicit trapezoidal rule
is recovered, which is unconditionally stable, second-order accurate
in
, and
yields implicit Euler backward method,
which is unconditionally stable, and first-order accurate
in
. deltaT is time step length used for
integration, nsteps parameter specifies
number of time steps to be solved.
Nonzero value of optional parameter renumber turns on the equation renumbering to optimize the profile of characteristic matrix (uses Sloan algorithm). By default, profile optimization is not performed. It will not work in parallel mode. If lumpedcapa is set, then the stabilization of numerical algorithm using lumped capacity matrix will be used, reducing the initial oscillations. See section 6.11 for an explanation on exportfields.
This linear transport problem supports changes in number of equations. It is possible to impose/remove Dirichlet boundary conditions during solution. This feature is enabled with changingProblemSize, which ensures storing solution values on nodes (DoFs) directly. If the problem does not grow/decrease during solution, it is more efficient to use conventional solution strategy and the parameter should not be mentioned.
Note: This problem type requires transport module and it can be used only when this module is included in your oofem configuration.
Borek Patzak 2011-12-29