General linear methods: connection to one step methods and invariant curves.
General theorems for numerical approximation of stochastic processes on the Hilbert space .
Geometric integrators for piecewise smooth Hamiltonian systems
In this paper, we consider C1,1 Hamiltonian systems. We prove the existence of a first derivative of the flow with respect to initial values and show that it satisfies the symplecticity condition almost everywhere in the phase-space. In a second step, we present a geometric integrator for such systems (called the SDH method) based on B-splines interpolation and a splitting method introduced by McLachlan and Quispel [Appl. Numer. Math. 45 (2003) 411–418], and we prove it is convergent, and that...
Half-explicit Runge-Kutta methods for semi-explicit differential-algebraic equations of index 1.
High order explicit Runge-Kutta pairs for ephemerides of the solar system and the Moon.
Higher-index differential-algebraic equations: analysis and numerical treatment
Highly anisotropic nonlinear temperature balance equation and its numerical solution using asymptotic-preserving schemes of second order in time
This paper deals with the numerical study of a nonlinear, strongly anisotropic heat equation. The use of standard schemes in this situation leads to poor results, due to the high anisotropy. An Asymptotic-Preserving method is introduced in this paper, which is second-order accurate in both, temporal and spacial variables. The discretization in time is done using an L-stable Runge−Kutta scheme. The convergence of the method is shown to be independent of the anisotropy parameter , and this for fixed...
Implementing Adams methods with preassigned stepsize ratios.
Implicit Runge-Kutta methods for transferable differential-algebraic equations
The numerical solution of transferable differential-algebraic equations (DAE’s) by implicit Runge-Kutta methods (IRK) is studied. If the matrix of coefficients of an IRK is non-singular then the arising systems of nonlinear equations are uniquely solvable. These methods are proved to be stable if an additional contractivity condition is satisfied. For transferable DAE’s with smooth solution we get convergence of order , where is the classical order of the IRK and is the stage order. For transferable...
Implicit-explicit Runge–Kutta schemes and finite elements with symmetric stabilization for advection-diffusion equations
We analyze a two-stage implicit-explicit Runge–Kutta scheme for time discretization of advection-diffusion equations. Space discretization uses continuous, piecewise affine finite elements with interelement gradient jump penalty; discontinuous Galerkin methods can be considered as well. The advective and stabilization operators are treated explicitly, whereas the diffusion operator is treated implicitly. Our analysis hinges on L2-energy estimates on discrete functions in physical space. Our main...
Implicit-explicit Runge–Kutta schemes and finite elements with symmetric stabilization for advection-diffusion equations
We analyze a two-stage implicit-explicit Runge–Kutta scheme for time discretization of advection-diffusion equations. Space discretization uses continuous, piecewise affine finite elements with interelement gradient jump penalty; discontinuous Galerkin methods can be considered as well. The advective and stabilization operators are treated explicitly, whereas the diffusion operator is treated implicitly. Our analysis hinges on L2-energy estimates on discrete functions in physical space. Our main...
Implicit-explicit Runge–Kutta schemes and finite elements with symmetric stabilization for advection-diffusion equations
We analyze a two-stage implicit-explicit Runge–Kutta scheme for time discretization of advection-diffusion equations. Space discretization uses continuous, piecewise affine finite elements with interelement gradient jump penalty; discontinuous Galerkin methods can be considered as well. The advective and stabilization operators are treated explicitly, whereas the diffusion operator is treated implicitly. Our analysis hinges on L2-energy estimates on discrete functions in physical space. Our main...
-error analysis for a system of quasivariational inequalities with noncoercive operators.
La résolution numérique des problèmes différentiels aux conditions initiales par des méthodes de collocation
Linear Multistep Methods for a Class of Functional Differential Equations.
Neurons: a numerical approach.
New Resolution Strategy for Multi-scale Reaction Waves using Time Operator Splitting and Space Adaptive Multiresolution: Application to Human Ischemic Stroke*
We tackle the numerical simulation of reaction-diffusion equations modeling multi-scale reaction waves. This type of problems induces peculiar difficulties and potentially large stiffness which stem from the broad spectrum of temporal scales in the nonlinear chemical source term as well as from the presence of large spatial gradients in the reactive fronts, spatially very localized. A new resolution strategy was recently introduced ? that combines...
Nine-stage Multi-derivative Runge-Kutta method of order 12
Numerical integration of differential equations in the presence of first integrals: observer method
We introduce a simple and powerful procedure-the observer method-in order to obtain a reliable method of numerical integration over an arbitrary long interval of time for systems of ordinary differential equations having first integrals. This aim is achieved by a modification of the original system such that the level manifold of the first integrals becomes a local attractor. We provide a theoretical justification of this procedure. We report many tests and examples dealing with a large spectrum...
Numerical modelling of steady and unsteady flows of generalized Newtonian fluids
This work presents the numerical solution of laminar incompressible viscous flow in a three dimensional branching channel with circular cross section for generalized Newtonian fluids. This model can be generalized by cross model in shear thinning meaning. The governing system of equations is based on the system of balance laws for mass and momentum. Numerical tests are performed on a three dimensional geometry, the branching channel with one entrance and two outlet parts. Numerical solution of the...