8,624 research outputs found
Square-root filtering via covariance SVD factors in the accurate continuous-discrete extended-cubature Kalman filter
This paper continues our research devoted to an accurate nonlinear Bayesian
filters' design. Our solution implies numerical methods for solving ordinary
differential equations (ODE) when propagating the mean and error covariance of
the dynamic state. The key idea is that an accurate implementation strategy
implies the methods with a discretization error control involved. This means
that the filters' moment differential equations are to be solved accurately,
i.e. with negligible error. In this paper, we explore the continuous-discrete
extended-cubature Kalman filter that is a hybrid method between Extended and
Cubature Kalman filters (CKF). Motivated by recent results obtained for the
continuous-discrete CKF in Bayesian filtering realm, we propose the numerically
stable (to roundoff) square-root approach within a singular value decomposition
(SVD) for the hybrid filter. The new method is extensively tested on a few
application examples including stiff systems
Multiple Right-Hand Side Techniques in Semi-Explicit Time Integration Methods for Transient Eddy Current Problems
The spatially discretized magnetic vector potential formulation of
magnetoquasistatic field problems is transformed from an infinitely stiff
differential algebraic equation system into a finitely stiff ordinary
differential equation (ODE) system by application of a generalized Schur
complement for nonconducting parts. The ODE can be integrated in time using
explicit time integration schemes, e.g. the explicit Euler method. This
requires the repeated evaluation of a pseudo-inverse of the discrete curl-curl
matrix in nonconducting material by the preconditioned conjugate gradient (PCG)
method which forms a multiple right-hand side problem. The subspace projection
extrapolation method and proper orthogonal decomposition are compared for the
computation of suitable start vectors in each time step for the PCG method
which reduce the number of iterations and the overall computational costs.Comment: 4 pages, 5 figure
Exponential Runge-Kutta methods for stiff kinetic equations
We introduce a class of exponential Runge-Kutta integration methods for
kinetic equations. The methods are based on a decomposition of the collision
operator into an equilibrium and a non equilibrium part and are exact for
relaxation operators of BGK type. For Boltzmann type kinetic equations they
work uniformly for a wide range of relaxation times and avoid the solution of
nonlinear systems of equations even in stiff regimes. We give sufficient
conditions in order that such methods are unconditionally asymptotically stable
and asymptotic preserving. Such stability properties are essential to guarantee
the correct asymptotic behavior for small relaxation times. The methods also
offer favorable properties such as nonnegativity of the solution and entropy
inequality. For this reason, as we will show, the methods are suitable both for
deterministic as well as probabilistic numerical techniques
Flux Splitting for stiff equations: A notion on stability
For low Mach number flows, there is a strong recent interest in the
development and analysis of IMEX (implicit/explicit) schemes, which rely on a
splitting of the convective flux into stiff and nonstiff parts. A key
ingredient of the analysis is the so-called Asymptotic Preserving (AP)
property, which guarantees uniform consistency and stability as the Mach number
goes to zero. While many authors have focussed on asymptotic consistency, we
study asymptotic stability in this paper: does an IMEX scheme allow for a CFL
number which is independent of the Mach number? We derive a stability criterion
for a general linear hyperbolic system. In the decisive eigenvalue analysis,
the advective term, the upwind diffusion and a quadratic term stemming from the
truncation in time all interact in a subtle way. As an application, we show
that a new class of splittings based on characteristic decomposition, for which
the commutator vanishes, avoids the deterioration of the time step which has
sometimes been observed in the literature
Physical and numerical sources of computational inefficiency in integration of chemical kinetic rate equations: Etiology, treatment and prognosis
The design of a very fast, automatic black-box code for homogeneous, gas-phase chemical kinetics problems requires an understanding of the physical and numerical sources of computational inefficiency. Some major sources reviewed in this report are stiffness of the governing ordinary differential equations (ODE's) and its detection, choice of appropriate method (i.e., integration algorithm plus step-size control strategy), nonphysical initial conditions, and too frequent evaluation of thermochemical and kinetic properties. Specific techniques are recommended (and some advised against) for improving or overcoming the identified problem areas. It is argued that, because reactive species increase exponentially with time during induction, and all species exhibit asymptotic, exponential decay with time during equilibration, exponential-fitted integration algorithms are inherently more accurate for kinetics modeling than classical, polynomial-interpolant methods for the same computational work. But current codes using the exponential-fitted method lack the sophisticated stepsize-control logic of existing black-box ODE solver codes, such as EPISODE and LSODE. The ultimate chemical kinetics code does not exist yet, but the general characteristics of such a code are becoming apparent
Algorithmic aspects of transient heat transfer problems in structures
It is noted that the application of finite element or finite difference techniques to the solution of transient heat transfer problems in structures often results in a stiff system of ordinary differential equations. Such systems are usually handled most efficiently by implicit integration techniques which require the solution of large and sparse systems of algebraic equations. The assembly and solution of these systems using the incomplete Cholesky conjugate gradient algorithm is examined. Several examples are used to demonstrate the advantage of the algorithm over other techniques
- …