American Institute of Mathematical Sciences

In this paper we investigate splitting methods in the presence of non-homogeneous boundary conditions. In particular, we consider the corrections that have been described and analyzed in Einkemmer, Ostermann 2015 and Alonso-Mallo, Cano, Reguera 2016. The latter method is extended to the non-linear case, and a rigorous convergence analysis is provided. We perform numerical simulations for diffusion-reaction, advection-reaction, and dispersion-reaction equations in order to evaluate the relative performance of these two corrections. Furthermore, we introduce an extension of both methods to obtain order three locally and evaluate under what circumstances this is beneficial.

We propose novel positive numerical integrators for approximating predator-prey models. The schemes are based on suitable symplectic procedures applied to the dynamical system written in terms of the log transformation of the original variables. Even if this approach is not new when dealing with Hamiltonian systems, it is of particular interest in population dynamics since the positivity of the approximation is ensured without any restriction on the temporal step size. When applied to separable M-systems, the resulting schemes are proved to be explicit, positive, Poisson maps. The approach is generalized to predator-prey dynamics which do not exhibit an M-system structure and successively to reaction-diffusion equations describing spatially extended dynamics. A classical polynomial Krylov approximation for the diffusive term joint with the proposed schemes for the reaction, allows us to propose numerical schemes which are explicit when applied to well established ecological models for predator-prey dynamics. Numerical simulations show that the considered approach provides results which outperform the numerical approximations found in recent literature.

The equations describing engineering and real-life models are usually derived in an approximated way. Thus, in most cases it is necessary to deal with equations containing some kind of perturbation. In this paper we consider fractional differential equations and study the effects on the continuous and numerical solution, of perturbations on the given function, over long-time intervals. Some bounds on the global error are also determined.

The paper is focused on the analysis of stability properties of a family of numerical methods designed for the numerical solution of stochastic Volterra integral equations. Stability properties are provided with respect to the basic test equation, as well as to the convolution test equation. For each equation, stability properties are intended both in the mean-square and in the asymptotic sense. Stability regions are also provided for a selection of methods. Numerical experiments confirming the theoretical study are also given.

We propose two-step collocation methods for the numerical solution of fractional differential equations. These methods increase the order of convergence of one-step collocation methods, with the same number of collocation points. Moreover, they are continuous methods, i.e. they furnish an approximation of the solution at each point of the time interval. We describe the derivation of two-step collocation methods and analyse convergence. Some numerical experiments confirm theoretical expectations.

A recent pseudospectral collocation is used to reduce a nonlinear delay differential equation to a system of ordinary differential equations. Standard methods are then applied to compute Lyapunov exponents. The validity of this simple approach is shown experimentally. Matlab codes are also included.

We consider a model planar system with discontinuous right-hand side possessing an attracting periodic orbit, and we investigate what happens to a Euler discretization with stepsize \begin{document}$ τ$\end{document} of this system. We show that, in general, the resulting discrete dynamical system does not possess an invariant curve, in sharp contrast to what happens for smooth problems. In our context, we show that the numerical trajectories are forced to remain inside a band, whose width is proportional to the discretization stepsize \begin{document}$ τ$\end{document}. We further show that if we consider an event-driven discretization of the model problem, whereby the solution is forced to step exactly on the discontinuity line, then there is a discrete periodic solution near the one of the original problem (for sufficiently small \begin{document}$ τ$\end{document}). Finally, we consider what happens to the Euler discretization of the regularized system rewritten in polar coordinates, and give numerical evidence that the discrete solution now undergoes a period doubling cascade with respect to the regularization parameter \begin{document}$ \epsilon$\end{document}, for fixed \begin{document}$ τ$\end{document}.

The paper aims to explore the long-term behaviour of stochastic two-step methods applied to a class of second order stochastic differential equations. In particular, the treatment focuses on preserving long-term statistics related to the dynamics of a linear stochastic damped oscillator whose velocity, in the stationary regime, is distributed as a Gaussian variable and uncorrelated with the position. By computing the solution of a very simple matrix equality, we a-priori determine the long-term statistics characterizing the numerical dynamics and analyze the behaviour of a selection of methods.

In this article we present, for the first time, domain-growth induced pattern formation for reaction-diffusion systems with linear cross-diffusion on evolving domains and surfaces. Our major contribution is that by selecting parameter values from spaces induced by domain and surface evolution, patterns emerge only when domain growth is present. Such patterns do not exist in the absence of domain and surface evolution. In order to compute these domain-induced parameter spaces, linear stability theory is employed to establish the necessary conditions for domain-growth induced cross-diffusion-driven instability for reaction-diffusion systems with linear cross-diffusion. Model reaction-kinetic parameter values are then identified from parameter spaces induced by domain-growth only; these exist outside the classical standard Turing space on stationary domains and surfaces. To exhibit these patterns we employ the finite element method for solving reaction-diffusion systems with cross-diffusion on continuously evolving domains and surfaces.

Magnetohydrodynamics (MHD) studies the dynamics of electrically conducting fluids, involving Navier-Stokes equations coupled with Maxwell equations via Lorentz force and Ohm's law. Monolithic methods, which solve fully coupled MHD systems, are computationally expensive. Partitioned methods, on the other hand, decouple the full system and solve subproblems in parallel, and thus reduce the computational cost.

This paper is devoted to the design and analysis of a partitioned method for the MHD system in the Elsässer variables. The stability analysis shows that for magnetic Prandtl number of order unity, the method is unconditionally stable. We prove the error estimates and present computational tests that support the theory.

In this work, we introduce a method based on piecewise polynomial interpolation to enclose rigorously solutions of nonlinear ODEs. Using a technique which we call a priori bootstrap, we transform the problem of solving the ODE into one of looking for a fixed point of a high order smoothing Picard-like operator. We then develop a rigorous computational method based on a Newton-Kantorovich type argument (the radii polynomial approach) to prove existence of a fixed point of the Picard-like operator. We present all necessary estimates in full generality and for any nonlinearities. With our approach, we study two systems of nonlinear equations: the Lorenz system and the ABC flow. For the Lorenz system, we solve Cauchy problems and prove existence of periodic and connecting orbits at the classical parameters, and for ABC flows, we prove existence of ballistic spiral orbits.

We generalize the theory of underlying one-step methods to strictly stable general linear methods (GLMs) solving nonautonomous ordinary differential equations (ODEs) that satisfy a global Lipschitz condition. We combine this theory with the Lyapunov and Sacker-Sell spectral stability theory for one-step methods developed in [34,35,36] to analyze the stability of a strictly stable GLM solving a nonautonomous linear ODE. These results are applied to develop a stability diagnostic for the solution of nonautonomous linear ODEs by strictly stable GLMs.

In this paper, we study the relative error propagation in the solution of linear autonomous ordinary differential equations with respect to perturbations in the initial value. We also consider equations with a constant forcing term and a nonzero equilibrium. The study is carried out for equations defined by normal matrices.

Several models in the applied sciences are characterized by instantaneous changes in the solutions or discontinuities in the vector field. Knowledge of the geometry of interaction of the flow with the discontinuities can give new insights on the behaviour of such systems. Here, we focus on the class of the piecewise smooth systems of Filippov type. We describe some numerical techniques to locate crossing and sliding regions on the discontinuity surface, to compute the sets of attraction of these regions together with the mathematical form of the separatrices of such sets. Some numerical tests will illustrate our approach.

We consider a smooth planar system having slow-fast motion, where the slow motion takes place near a curve γ. We explore the idea of replacing the original smooth system with a system with discontinuous right-hand side (DRHS system for short), whereby the DRHS system coincides with the smooth one away from a neighborhood of γ. After this reformulation, in the region of phase-space where γ is attracting for the DRHS system, we will obtain sliding motion on γ and numerical methods apt at integrating for sliding motion can be applied. Moreover, we further bypass resolving the sliding motion and monitor entries (transversal) and exits (tangential) on the curve γ, a fact that can be done independently of resolving for the motion itself. The end result is a method free from the need to adopt stiff integrators or to worry about resolving sliding motion for the DRHS system. We illustrate the performance of our method on a few problems, highlighting the feasibility of using simple explicit Runge-Kutta schemes, and that we obtain much the same orbits of the original smooth system.