Related papers: Transforming Stiffness and Chaos
In this introductory work I will present the Finite Difference method for hyperbolic equations, focusing on a method which has second order precision both in time and space (the so-called staggered leapfrog method) and applying it to the…
Locally refined meshes impose severe stability constraints on explicit time-stepping methods for the numerical simulation of time dependent wave phenomena. Local time-stepping methods overcome that bottleneck by using smaller time-steps…
Applied to the master equation, the usual numerical integration methods, such as Runge-Kutta, become inefficient when the rates associated with various transitions differ by several orders of magnitude. We introduce an integration scheme…
A fundamental requirement for the emergence of classical behavior from an underlying quantum description is that certain observed quantum systems make a transition to chaotic dynamics as their action is increased relative to $\hbar$. While…
This paper develops an $\alpha$-parametrized framework for analyzing the strong convergence of the stochastic theta (ST) method for stochastic differential equations driven by time-changed L\'evy noise (TCSDEwLNs) with time-space-dependent…
We show that any solution of a smooth first order dynamic equation can be made Lyapunov stable at will by the choice of an appropriate time-dependent Riemannian metric.
We show that Newton methods for generalized equations are input-to-state stable with respect to disturbances such as due to inexact computations. We then use this result to obtain convergence and robustness of a multistep Newton-type method…
Lyapunov exponents measure the average exponential growth rate of typical linear perturbations in a chaotic system, and the inverse of the largest exponent is a measure of the time horizon over which the evolution of the system can be…
Time integration of Fourier pseudo-spectral DNS is usually performed using the classical fourth-order accurate Runge--Kutta method, or other methods of second or third order, with a fixed step size. We investigate the use of higher-order…
This paper is part of a program to combine a staggered time and staggered spatial discretization of continuum wave equations so that important properties of the continuum that are proved using vector calculus can be proven in an analogous…
We decide the stability and compute the Lyapunov exponent of continuous-time linear switching systems with a guaranteed dwell time. The main result asserts that the discretization method with step size~$h$ approximates the Lyapunov exponent…
We prove the theorem of linearized asymptotic stability for fractional differential equations. More precisely, we show that an equilibrium of a nonlinear Caputo fractional differential equation is asymptotically stable if its linearization…
A stabilized finite element method is introduced for the simulation of time-periodic creeping flows, such as those found in the cardiorespiratory systems. The new technique, which is formulated in the frequency rather than time domain,…
Approximation theory for Lyapunov and Sacker-Sell spectra based upon QR techniques is used to analyze the stability of a one-step method solving a time-dependent, linear, ordinary differential equation (ODE) initial value problem in terms…
We show that for continuous time dynamical systems described by polynomial differential equations of modest degree (typically equal to three), the following decision problems which arise in numerous areas of systems and control theory…
In problems where the temporal evolution of a nonlinear system cannot be followed, a method for studying the fluctuations of spatial patterns has been developed. That method is applied to well-known problems in deterministic chaos (the…
It is known in \cite{beccari} that the standard explicit Euler-type scheme (such as the exponential Euler and the linear-implicit Euler schemes) with a uniform timestep, though computationally efficient, may diverge for the stochastic…
In ordinary turbulence research it has been a long standing tradition to solve the equations in spectral space giving the best possible accuracy. This is indeed a natural choice for incompressible problems with periodic boundaries, but it…
In this paper we address the question of whether it is possible to integrate time-dependent high-dimensional PDEs with hierarchical tensor methods and explicit time stepping schemes. To this end, we develop sufficient conditions for…
We present fifth order Runge-Kutta-Nystr\"om methods, where we allow the timestep coefficients to assume complex values. Among the methods with complex timesteps, we focus on the ones with the coefficients that have positive real parts.…