65L80 Methods for differential-algebraic equations
Refine
Language
- English (12)
Keywords
- strangeness index (5)
- differential-algebraic equation (4)
- Lyapunov exponent (3)
- Sacker-Sell spectrum (3)
- exponential dichotomy (3)
- spectral interval (3)
- Bohl exponent (2)
- Hamiltonian matrix (2)
- Steklov function (2)
- canonical form (2)
This paper is devoted to the numerical approximation of Lyapunov and Sacker-Sell spectral intervals for linear differential-algebraic equations (DAEs). The spectral analysis for DAEs is improved and the concepts of leading directions and solution subspaces associated with spectral intervals are extended to DAEs. Numerical methods
based on smooth singular value decompositions are introduced for computing all or only some spectral intervals and their associated leading directions. The numerical algorithms as well as implementation issues are discussed in detail and numerical examples are presented to illustrate the theoretical results.
In this paper, we propose and investigate numerical methods based on QR factorization for computing all or some Lyapunov or Sacker-Sell spectral intervals for
linear differential-algebraic equations.
Furthermore, a perturbation and error analysis for these methods is presented. We
investigate how errors in the data and in the numerical integration affect the
accuracy of the approximate spectral intervals. Although we need to integrate
numerically some differential-algebraic systems on usually very long
time-intervals, under certain assumptions, it is shown that the error of the
computed spectral intervals can be controlled by the local error of numerical
integration and the error in solving the algebraic constraint.
Some numerical examples are presented to illustrate the theoretical results.
We discuss the solution of linear second order differential-algebraic equations with variable coefficients.
Since index reduction and order reduction for higher order higher index differential-algebraic systems do not commute, appropriate index reduction methods for higher order DAEs are required.
We present an index reduction method based on derivative arrays that allows to determine
an equivalent second order system of lower index in a numerical computable way.
For such an equivalent second order system
an appropriate order reduction method allows to formulate a suitable first order DAE system of low index
that has the same solution components as the original second order system.
We consider hybrid systems of differential-algebraic equations and present
a general framework for general nonlinear over- and underdetermined hybrid
systems that allows the
analysis of existence and uniqueness and the application of index reduction
methods for hybrid differential-algebraic systems.
A particular difficulty in the numerical simulation of hybrid systems is
(numerical) chattering, i.e., fast oscillations between modes of operations.
A regularization technique using sliding modes allows to regularize the
system behavior in the case of chattering.
Further, we show how chattering behavior during the numerical solution can
be prevented using sliding mode simulation. The advantage of the sliding mode
simulation is illustrated by numerical examples.
Canonical forms for matrix triples $(A,G,\hat G)$, where
$A$ is arbitrary rectangular and $G$, $\hat G$ are either real symmetric
or skew symmetric, or complex Hermitian or skew Hermitian, are derived.
These forms generalize classical product Schur forms as well as
singular value decompositions.
An new proof for the complex case is given, where there is no need to
distinguish whether $G$ and $\hat G$ are Hermitian or skew Hermitian.
This proof is independent from the results in Bolschakov/Reichstein 1995, where
a similar canonical form has been obtained for the complex case,
and it allows generalization to the real case. Here,
the three cases, i.e., that
$G$ and $\hat G$ are both symmetric, both skew symmetric or one each,
are treated separately.
The classical singular value decomposition for a matrix $A\in\Cmn$ is a
canonical form for $A$ that also displays the eigenvalues
of the Hermitian matrices $AA^\ast$ and $A^\ast A$. In this paper, we develop
a corresponding decomposition for $A$ that provides the Jordan canonical forms
for the complex symmetric matrices $AA^T$ and $A^TA$. More generally, we consider
the matrix triple $(A,G_1,G_2)$, where $G_1\in\CC{m}, G_2\in\CC{n}$
are invertible and either complex symmetric and complex skew-symmetric, and we
provide a canonical form under transformations of the form
$(A,G_1,G_2)\mapsto(X^T A Y, X^T G_1X, Y^T G_2Y)$, where $X,Y$ are nonsingular.
In this paper we present the new numerical algorithm GEOMS for the numerical integration of the most general form of the equations of motion of multibody systems, including nonholonomic constraints and possible redundancies in the constraints, as they may appear in industrial applications. Besides the numerical integration it offers some additional features like stabilization of the model equations, use of different decomposition strategies, or checking and correction of the initial values with respect to their consistency. Furthermore, GEOMS preserves hidden constraints and (possibly) existing solution invariants if they are provided as equations.
We will also demonstrate the performance and the applicability of GEOMS for two mechanical examples of different degrees of complexity.
Lyapunov and exponential dichotomy spectral theory is extended
from ordinary differential equations (ODEs) to nonautonomous
differential-algebraic equations (DAEs). By using orthogonal
changes of variables, the original DAE system is transformed into
appropriate condensed forms, for which concepts such as Lyapunov
exponents, Bohl exponents, exponential dichotomy and spectral
intervals of various kinds can be analyzed via the resulting
underlying ODE. Some essential differences between the spectral
theory for ODEs and that for DAEs are pointed out. Numerical
methods for computing the spectral intervals associated with
Lyapunov and Sacker-Sell (exponential dichotomy) spectra are
derived by modifying and extending those methods proposed for ODEs. Perturbation theory and error analysis are discussed, as
well. Finally, some numerical examples are presented to illustrate
the theoretical results and the properties of the numerical
methods.
We discuss the eigenvalue problem for
general and structured matrix polynomials which may
be singular and may have eigenvalues at infinity.
We derive staircase
condensed forms that allow deflation of the infinite eigenvalue and
singular structure of the matrix polynomial.
The remaining reduced order staircase form leads to
new types of linearizations which determine the finite eigenvalues and
and corresponding eigenvectors. The new linearizations
also simplify the construction of structure preserving linearizations.
We review some known results for POD model reduction applied to ODEs. Then, these
results are generalized to several types of DAEs. We provide algorithms for the
model reduction and error bounds for the reduced order models. Some limits of
the approach are pointed out and alternative methods for reduced order subspace approximation
are presented. The POD approach is tested and evaluated for a medium sized DAE example
from multibody dynamics.