Related papers: Computing the solutions of the van der Pol equatio…
We develop a validated numerical procedure for continuation of local stable/unstable manifold patches attached to equilibrium solutions of ordinary differential equations. The procedure has two steps. First we compute an accurate high order…
The use of proper orthogonal decomposition (POD) to explore the complex fluid flows that are common in engineering applications is increasing and has yielded new physical insights. However, for most engineering systems the dimension of the…
Probabilistic numerical solvers for ordinary differential equations (ODEs) treat the numerical simulation of dynamical systems as problems of Bayesian state estimation. Aside from producing posterior distributions over ODE solutions and…
This paper deals with the asymptotic study of the so-called canard solutions, which arise in the study of real singularly perturbed ODEs. Starting near an attracting branch of the "slow curve", those solutions are crossing a turning point…
Two-step predictor/corrector methods are provided to solve three classes of problems that present themselves as systems of ordinary differential equations (ODEs). In the first class, velocities are given from which displacements are to be…
An adaptation of the arbitrary high order ADER-DG numerical method with local DG predictor for solving the IVP for a first-order non-linear ODE system is proposed. The proposed numerical method is a completely one-step ODE solver with…
We consider a Hamiltonian system of coupled Van der Pol-Duffing(VdPD) oscillators with balanced loss and gain. The system is analyzed perturbatively by using Renormalization Group(RG) techniques as well as Multiple Scale Analysis(MSA). Both…
The E(xplicit)I(implicit)N(null) method was developed recently to remove numerical instability from PDEs, adding and subtracting an operator $\mathcal{D}$ of arbitrary structure, treating the operator implicitly in one case, and explicitly…
Neural ordinary differential equations (Neural ODEs) define continuous time dynamical systems with neural networks. The interest in their application for modelling has sparked recently, spanning hybrid system identification problems and…
The standard series expansion for the period of a finite amplitude pendulum as a function of energy (and hence amplitude) provides a lower limit on the period when the series is truncated. An adjustment to the last term in the truncated…
In this paper we define a small variation of the Taylor method and a formula for the global error of this new numerical method that allows us to keep track of the round-off error and does not require previous knowledge of the exact…
Numerical solutions to wave-type PDEs utilizing method-of-lines require the ODE solver's stability domain to include a large stretch of the imaginary axis surrounding the origin. We show here that extrapolation based solvers of…
In a previous paper the \textit{real} evolution of the system of ODEs \ddot{z}_{n} + z_{n}=\sum\limits_{m = 1, m \ne n}^{N} g_{nm}{(z_{n} - z_{m})} ^{- 3}, z_{n} \equiv z_{n}(t), \qquad \dot {z}_{n} \equiv \frac{d z_{n}(t)}{dt}, \qquad n =…
We present a novel method for calculating Pad\'e approximants that is capable of eliminating spurious poles placed at the point of development and of identifying and eliminating spurious poles created by precision limitations and/or noisy…
In this article, we introduce a novel parallel-in-time solver for nonlinear ordinary differential equations (ODEs). We state the numerical solution of an ODE as a root-finding problem that we solve using Newton's method. The affine…
We propose a fast and scalable optimization method to solve chance or probabilistic constrained optimization problems governed by partial differential equations (PDEs) with high-dimensional random parameters. To address the critical…
Virial expansions are the series in powers of density assumed to be small. However, the equations of state require to consider finite densities for which virial expansions, as a rule, diverge. In order to extrapolate a virial expansion to…
We present general series solutions to the Tolman-Oppenheimer-Volkoff equations for compact stellar objects. We develop an algorithm to compute the coefficients of the power series in terms of the equation of state and its derivatives with…
The solution of systems of non-autonomous linear ordinary differential equations is crucial in a variety of applications, such us nuclear magnetic resonance spectroscopy. A new method with spectral accuracy has been recently introduced in…
The study of strongly out-of-equilibrium states and their time evolution towards thermalization is critical to the understanding of an ever widening range of physical processes. We present a numerical method that for the first time allows…