Stability analysis of a singularly perturbed coupled ODE-PDE system

International audienceThis paper is concerned with a coupled ODE-PDE system with two time scales modeled by a perturbation parameter. Firstly, the perturbation parameter is introduced into the PDE system. We show that the stability of the full system is guaranteed by the stability of the reduced and the boundary-layer subsystems. A numerical simulation on a gas flow transport model is used to illustrate the first result. Secondly, an example is used to show that the full system can be unstable even though both subsystems are stable when the perturbation parameter is introduced into the ODE system

Unfolding symmetric Bogdanov-Takens bifurcations for front dynamics in a reaction-diffusion system

This manuscript extends the analysis of a much studied singularly perturbed three-component reaction-diffusion system for front dynamics in the regime where the essential spectrum is close to the origin. We confirm a conjecture from a preceding paper by proving that the triple multiplicity of the zero eigenvalue gives a Jordan chain of length three. Moreover, we simplify the center manifold reduction and computation of the normal form coefficients by using the Evans function for the eigenvalues. Finally, we prove the unfolding of a Bogdanov-Takens bifurcation with symmetry in the model. This leads to stable periodic front motion, including stable traveling breathers, and these results are illustrated by numerical computations.Comment: 39 pages, 7 figure

Spatial Manifestations of Order Reduction in Runge-Kutta Methods for Initial Boundary Value Problems

This paper studies the spatial manifestations of order reduction that occur when time-stepping initial-boundary-value problems (IBVPs) with high-order Runge-Kutta methods. For such IBVPs, geometric structures arise that do not have an analog in ODE IVPs: boundary layers appear, induced by a mismatch between the approximation error in the interior and at the boundaries. To understand those boundary layers, an analysis of the modes of the numerical scheme is conducted, which explains under which circumstances boundary layers persist over many time steps. Based on this, two remedies to order reduction are studied: first, a new condition on the Butcher tableau, called weak stage order, that is compatible with diagonally implicit Runge-Kutta schemes; and second, the impact of modified boundary conditions on the boundary layer theory is analyzed.Comment: 41 pages, 9 figure

Differential-Algebraic Equations and Beyond: From Smooth to Nonsmooth Constrained Dynamical Systems

The present article presents a summarizing view at differential-algebraic equations (DAEs) and analyzes how new application fields and corresponding mathematical models lead to innovations both in theory and in numerical analysis for this problem class. Recent numerical methods for nonsmooth dynamical systems subject to unilateral contact and friction illustrate the topicality of this development.Comment: Preprint of Book Chapte

Emergence of steady and oscillatory localized structures in a phytoplankton-nutrient model

Co-limitation of marine phytoplankton growth by light and nutrient, both of which are essential for phytoplankton, leads to complex dynamic behavior and a wide array of coherent patterns. The building blocks of this array can be considered to be deep chlorophyll maxima, or DCMs, which are structures localized in a finite depth interior to the water column. From an ecological point of view, DCMs are evocative of a balance between the inflow of light from the water surface and of nutrients from the sediment. From a (linear) bifurcational point of view, they appear through a transcritical bifurcation in which the trivial, no-plankton steady state is destabilized. This article is devoted to the analytic investigation of the weakly nonlinear dynamics of these DCM patterns, and it has two overarching themes. The first of these concerns the fate of the destabilizing stationary DCM mode beyond the center manifold regime. Exploiting the natural singularly perturbed nature of the model, we derive an explicit reduced model of asymptotically high dimension which fully captures these dynamics. Our subsequent and fully detailed study of this model - which involves a subtle asymptotic analysis necessarily transgressing the boundaries of a local center manifold reduction - establishes that a stable DCM pattern indeed appears from a transcritical bifurcation. However, we also deduce that asymptotically close to the original destabilization, the DCM looses its stability in a secondary bifurcation of Hopf type. This is in agreement with indications from numerical simulations available in the literature. Employing the same methods, we also identify a much larger DCM pattern. The development of the method underpinning this work - which, we expect, shall prove useful for a larger class of models - forms the second theme of this article