Numerical approximation of singular boundary value problems for a nonlinear differential equation
Numerical Approximations of the Relative Rearrangement: The piecewise linear case. Application to some Nonlocal Problems
We first prove an abstract result for a class of nonlocal problems using fixed point method. We apply this result to equations revelant from plasma physic problems. These equations contain terms like monotone or relative rearrangement of functions. So, we start the approximation study by using finite element to discretize this nonstandard quantities. We end the paper by giving a numerical resolution of a model containing those terms.
Numerical Computation of Branch Points in Ordinary Differential Equations.
Numerical Computation of Forced Oscillations in Coupled Duffing Equations.
Numerical computation of solitons for optical systems
In this paper, we present numerical methods for the determination of solitons, that consist in spatially localized stationary states of nonlinear scalar equations or coupled systems arising in nonlinear optics. We first use the well-known shooting method in order to find excited states (characterized by the number of nodes) for the classical nonlinear Schrödinger equation. Asymptotics can then be derived in the limits of either large are large nonlinear exponents . In a second part, we compute...
Numerical computation of solitons for optical systems
In this paper, we present numerical methods for the determination of solitons, that consist in spatially localized stationary states of nonlinear scalar equations or coupled systems arising in nonlinear optics. We first use the well-known shooting method in order to find excited states (characterized by the number k of nodes) for the classical nonlinear Schrödinger equation. Asymptotics can then be derived in the limits of either large k are large nonlinear exponents σ. In a second part, we compute...
Numerical computation of the eigenvalues for the spheroidal wave equation with accurate error estimation by matrix method.
Numerical experiments with different schemes for a singularly perturbed problem.
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 methods for approximating eigenvalues of boundary value problems.
Numerical model of a pine in a wind
Steady-state nonlinear differential equations govering the stem curve of a wind-loaded pine are derived and solved numerically. Comparison is made between the results computed and the data from photographs of a pine stem during strong wind. The pine breaking is solved at the end.
Numerical modelling of a bridge subjected to simultaneous effect of a moving load and a vertical seismic ground excitation
A simple beam subjected to a row of regularly distributed moving forces and simultaneous vertical motions of its supports is described using a simplified theoretical model and a finite differences approach. Several levels of simplification of the structure and input data are supposed. Numerical results confirm legitimacy of the assumptions.
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...
Numerical optimization of parameters in systems of differential equations
We present results on the estimation of unknown parameters in systems of ordinary differential equations in order to fit the output of models to real data. The numerical method is based on the nonlinear least squares problem along with the solution of sensitivity equations corresponding to the differential equations. We will present the performance of the method on the problem of fitting the output of basic compartmental epidemic models to data from the Covid-19 epidemic. This allows us to draw...
Numerical precision for differential inclusions with uniqueness
In this article, we show the convergence of a class of numerical schemes for certain maximal monotone evolution systems; a by-product of this results is the existence of solutions in cases which had not been previously treated. The order of these schemes is in general and when the only non Lipschitz continuous term is the subdifferential of the indicatrix of a closed convex set. In the case of Prandtl’s rheological model, our estimates in maximum norm do not depend on spatial dimension.
Numerical precision for differential inclusions with uniqueness
In this article, we show the convergence of a class of numerical schemes for certain maximal monotone evolution systems; a by-product of this results is the existence of solutions in cases which had not been previously treated. The order of these schemes is 1/2 in general and 1 when the only non Lipschitz continuous term is the subdifferential of the indicatrix of a closed convex set. In the case of Prandtl's rheological model, our estimates in maximum norm do not depend on spatial dimension. ...
Numerical simulation of generalized Newtonian fluids flow in bypass geometry
The aim of this work is to present numerical results of non-Newtonian fluid flow in a model of bypass. Different angle of a connection between narrowed channel and the bypass graft is considered. Several rheology viscosity models were used for the non-Newtonian fluid, namely the modified Cross model and the Carreau-Yasuda model. The results of non-Newtonian fluid flow are compared to the results of Newtonian fluid. The fundamental system of equations is the generalized system of Navier-Stokes equations...
Numerical simulation of gluey particles
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...
Numerical simulation of gluey particles
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...
Numerical solution of a 1-d elastohydrodynamic problem in magnetic storage devices
In this work we present new numerical methods to simulate the mechanics of head-tape magnetic storage devices. The elastohydrodynamic problem is formulated in terms of a coupled system which is governed by a nonlinear compressible Reynolds equation for the air pressure over the head, and a rod model for the tape displacement. A fixed point algorithm between the solutions of the elastic and hydrodynamic problems is proposed. For the nonlinear Reynolds equation, a characteristics method and a...