| NlTransientTransportProblem | nsteps(in) # |
| [renumber(in) #] | |
| deltaT(rn) # | |
| alpha(rn) # | |
| [lumpedcapa() #] | |
| [nsmax(in) #] | |
| rtol(rn) # | |
| [manrmsteps(in) #] | |
| [sparselinsolverparams() #] | |
| [exportfields() #] | |
| [atomicfields(in) #] |
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. Parameter maxiter determines the maximum
number of iterations allowed to reach equilibrium (default is 30).
The convergence is reached, when norms of both residual fluxes and iterative
change of solution vector is less than the value given by rtol.
If manrmsteps parameter is nonzero, then the modified
N-R scheme is used, with the left-hand side matrix updated after
manrmsteps steps.
If lumpedcapa is set, then the stabilization of numerical
algorithm using lumped capacity matrix will be used, reducing the
initial oscillations.
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.
See section 6.17 for explanation of
exportfields and atomicfields parameters.
Note: This problem type requires transport module and it can be used only when this module is included in your oofem configuration.
Borek Patzak 2009-08-24