7,542 research outputs found
Implicit-Explicit Runge-Kutta schemes for numerical discretization of optimal control problems
Implicit-explicit (IMEX) Runge-Kutta methods play a major rule in the
numerical treatment of differential systems governed by stiff and non-stiff
terms. This paper discusses order conditions and symplecticity properties of a
class of IMEX Runge-Kutta methods in the context of optimal control problems.
The analysis of the schemes is based on the continuous optimality system. Using
suitable transformations of the adjoint equation, order conditions up to order
three are proven as well as the relation between adjoint schemes obtained
through different transformations is investigated. Conditions for the IMEX
Runge-Kutta methods to be symplectic are also derived. A numerical example
illustrating the theoretical properties is presented
Feedback Stabilization Methods for the Numerical Solution of Systems of Ordinary Differential Equations
In this work we study the problem of step size selection for numerical
schemes, which guarantees that the numerical solution presents the same
qualitative behavior as the original system of ordinary differential equations,
by means of tools from nonlinear control theory. Lyapunov-based and Small-Gain
feedback stabilization methods are exploited and numerous illustrating
applications are presented for systems with a globally asymptotically stable
equilibrium point. The obtained results can be used for the control of the
global discretization error as well.Comment: 33 pages, 9 figures. Submitted for possible publication to BIT
Numerical Mathematic
Solving periodic semilinear stiff PDEs in 1D, 2D and 3D with exponential integrators
Dozens of exponential integration formulas have been proposed for the
high-accuracy solution of stiff PDEs such as the Allen-Cahn, Korteweg-de Vries
and Ginzburg-Landau equations. We report the results of extensive comparisons
in MATLAB and Chebfun of such formulas in 1D, 2D and 3D, focusing on fourth and
higher order methods, and periodic semilinear stiff PDEs with constant
coefficients. Our conclusion is that it is hard to do much better than one of
the simplest of these formulas, the ETDRK4 scheme of Cox and Matthews
- …