13,694 research outputs found
Discretization of Homogeneous Systems Using Euler Method with a State-Dependent Step
International audienceNumeric approximations to the solutions of asymptotically stable homogeneous systems by Euler method, with a step of discretization scaled by the state norm, are investigated (for the explicit and implicit integration schemes). It is proven that for a sufficiently small discretization step the convergence of the approximating solutions to zero can be guaranteed globally in a finite or a fixed time depending on the degree of homogeneity of the system, but in an infinite number of discretization iterations. The maximal admissible step can be estimated by analyzing the system properties on the sphere. It is shown that the absolute and relative errors of the discretizations are globally bounded functions, thus the approximations approaching the solutions with the step converging to zero. In addition, it is established that the proposed discretization approach preserves robustness with respect to exogenous perturbations. Efficiency of the designed discretization algorithms is demonstrated in simulations
Proofs for "Discretization of Homogeneous Systems Using Euler Method with a State-Dependent Step"
This note contains some proofs for the paper "Discretization of Homogeneous Systems Using Euler Method with a State-Dependent Step" of the same authors
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
Degenerate anisotropic elliptic problems and magnetized plasma simulations
This paper is devoted to the numerical approximation of a degenerate
anisotropic elliptic problem. The numerical method is designed for arbitrary
space-dependent anisotropy directions and does not require any specially
adapted coordinate system. It is also designed to be equally accurate in the
strongly and the mildly anisotropic cases. The method is applied to the
Euler-Lorentz system, in the drift-fluid limit. This system provides a model
for magnetized plasmas
Geometric integration on spheres and some interesting applications
Geometric integration theory can be employed when numerically solving ODEs or
PDEs with constraints. In this paper, we present several one-step algorithms of
various orders for ODEs on a collection of spheres. To demonstrate the
versatility of these algorithms, we present representative calculations for
reduced free rigid body motion (a conservative ODE) and a discretization of
micromagnetics (a dissipative PDE). We emphasize the role of isotropy in
geometric integration and link numerical integration schemes to modern
differential geometry through the use of partial connection forms; this
theoretical framework generalizes moving frames and connections on principal
bundles to manifolds with nonfree actions.Comment: This paper appeared in prin
New numerical solver for flows at various Mach numbers
Many problems in stellar astrophysics feature flows at low Mach numbers.
Conventional compressible hydrodynamics schemes frequently used in the field
have been developed for the transonic regime and exhibit excessive numerical
dissipation for these flows. While schemes were proposed that solve
hydrodynamics strictly in the low Mach regime and thus restrict their
applicability, we aim at developing a scheme that correctly operates in a wide
range of Mach numbers. Based on an analysis of the asymptotic behavior of the
Euler equations in the low Mach limit we propose a novel scheme that is able to
maintain a low Mach number flow setup while retaining all effects of
compressibility. This is achieved by a suitable modification of the well-known
Roe solver. Numerical tests demonstrate the capability of this new scheme to
reproduce slow flow structures even in moderate numerical resolution. Our
scheme provides a promising approach to a consistent multidimensional
hydrodynamical treatment of astrophysical low Mach number problems such as
convection, instabilities, and mixing in stellar evolution.Comment: 16 pages, 8 figures, accepted for publication by A&
A class of Galerkin schemes for time-dependent radiative transfer
The numerical solution of time-dependent radiative transfer problems is
challenging, both, due to the high dimension as well as the anisotropic
structure of the underlying integro-partial differential equation. In this
paper we propose a general framework for designing numerical methods for
time-dependent radiative transfer based on a Galerkin discretization in space
and angle combined with appropriate time stepping schemes. This allows us to
systematically incorporate boundary conditions and to preserve basic properties
like exponential stability and decay to equilibrium also on the discrete level.
We present the basic a-priori error analysis and provide abstract error
estimates that cover a wide class of methods. The starting point for our
considerations is to rewrite the radiative transfer problem as a system of
evolution equations which has a similar structure like first order hyperbolic
systems in acoustics or electrodynamics. This analogy allows us to generalize
the main arguments of the numerical analysis for such applications to the
radiative transfer problem under investigation. We also discuss a particular
discretization scheme based on a truncated spherical harmonic expansion in
angle, a finite element discretization in space, and the implicit Euler method
in time. The performance of the resulting mixed PN-finite element time stepping
scheme is demonstrated by computational results
- …