Generalized interpolation and some applications in ordinary differential equations
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...
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...
The purpose of this article is the analysis and the development of Eulerian multi-fluid models to describe the evolution of the mass density of evaporating liquid sprays. First, the classical multi-fluid model developed in [Laurent and Massot, Combust. Theor. Model.5 (2001) 537–572] is analyzed in the framework of an unsteady configuration without dynamical nor heating effects, where the evaporation process is isolated, since it is a key issue. The classical multi-fluid method consists then in...
This paper presents a mathematical model for photo-excited carrier decay in a semiconductor. Due to the carrier trapping states and recombination centers in the bandgap, the carrier decay process is defined by the system of nonlinear differential equations. The system of nonlinear ordinary differential equations is approximated by linearized backward Euler scheme. Some a priori estimates of the discrete solution are obtained and the convergence of the linearized backward Euler method is proved....
We propose here a model and a numerical scheme to compute the motion of rigid particles interacting through the lubrication force. In the case of a particle approaching a plane, we propose an algorithm and prove its convergence towards the solutions to the gluey particle model described in [B. Maury, ESAIM: Proceedings 18 (2007) 133–142]. We propose a multi-particle version of this gluey model which is based on the projection of the velocities onto a set of admissible velocities. Then, we describe...
We propose here a model and a numerical scheme to compute the motion of rigid particles interacting through the lubrication force. In the case of a particle approaching a plane, we propose an algorithm and prove its convergence towards the solutions to the gluey particle model described in [B. Maury, ESAIM: Proceedings18 (2007) 133–142]. We propose a multi-particle version of this gluey model which is based on the projection of the velocities onto a set of admissible velocities. Then, we describe...