A finite difference equivalent circuit approach to secondary current modelling in annular porous electrodes.
An unconditionally stable parallel difference scheme for telegraph equation.
Mathematical and numerical modelling of piezoelectric sensors
The present work aims at proposing a rigorous analysis of the mathematical and numerical modelling of ultrasonic piezoelectric sensors. This includes the well-posedness of the final model, the rigorous justification of the underlying approximation and the design and analysis of numerical methods. More precisely, we first justify mathematically the classical quasi-static approximation that reduces the electric unknowns to a scalar electric potential. We next justify the reduction of the computation...
Mathematical and numerical modelling of piezoelectric sensors
The present work aims at proposing a rigorous analysis of the mathematical and numerical modelling of ultrasonic piezoelectric sensors. This includes the well-posedness of the final model, the rigorous justification of the underlying approximation and the design and analysis of numerical methods. More precisely, we first justify mathematically the classical quasi-static approximation that reduces the electric unknowns to a scalar electric potential. We next justify the reduction of the computation...
Mathematical and numerical modelling of piezoelectric sensors
The present work aims at proposing a rigorous analysis of the mathematical and numerical modelling of ultrasonic piezoelectric sensors. This includes the well-posedness of the final model, the rigorous justification of the underlying approximation and the design and analysis of numerical methods. More precisely, we first justify mathematically the classical quasi-static approximation that reduces the electric unknowns to a scalar electric potential. We next justify the reduction of the computation...
Mathematical and numerical studies of non linear ferromagnetic materials
Mathematical and numerical studies of non linear ferromagnetic materials
In this paper we are interested in the numerical modeling of absorbing ferromagnetic materials obeying the non-linear Landau-Lifchitz-Gilbert law with respect to the propagation and scattering of electromagnetic waves. In this work we consider the 1D problem. We first show that the corresponding Cauchy problem has a unique global solution. We then derive a numerical scheme based on an appropriate modification of Yee's scheme, that we show to preserve some important properties of the continuous...
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 solution of a scattering problem in a wave guide using a grid generation technique.
Numerical solution of the Maxwell equations in time-varying media using Magnus expansion
For the Maxwell equations in time-dependent media only finite difference schemes with time-dependent conductivity are known. In this paper we present a numerical scheme based on the Magnus expansion and operator splitting that can handle time-dependent permeability and permittivity too. We demonstrate our results with numerical tests.
On the analysis of Bérenger’s perfectly matched layers for Maxwell’s equations
In this work, we investigate the Perfectly Matched Layers (PML) introduced by Bérenger [3] for designing efficient numerical absorbing layers in electromagnetism. We make a mathematical analysis of this model, first via a modal analysis with standard Fourier techniques, then via energy techniques. We obtain uniform in time stability results (that make precise some results known in the literature) and state some energy decay results that illustrate the absorbing properties of the model. This last...
On the analysis of Bérenger's Perfectly Matched Layers for Maxwell's equations
In this work, we investigate the Perfectly Matched Layers (PML) introduced by Bérenger [3] for designing efficient numerical absorbing layers in electromagnetism. We make a mathematical analysis of this model, first via a modal analysis with standard Fourier techniques, then via energy techniques. We obtain uniform in time stability results (that make precise some results known in the literature) and state some energy decay results that illustrate the absorbing properties of the model. This...
On the importance of solid deformations in convection-dominated liquid/solid phase change of pure materials
We analyse the effect of the mechanical response of the solid phase during liquid/solid phase change by numerical simulation of a benchmark test based on the well-known and debated experiment of melting of a pure gallium slab counducted by Gau & Viskanta in 1986. The adopted mathematical model includes the description of the melt flow and of the solid phase deformations. Surprisingly the conclusion reached is that, even in this case of pure material, the contribution of the solid phase to the...
Solutions to Maxwell equations for an electromagnetic beam injected into a plasma.
Stable upwind schemes for the magnetic induction equation
We consider the magnetic induction equation for the evolution of a magnetic field in a plasma where the velocity is given. The aim is to design a numerical scheme which also handles the divergence constraint in a suitable manner. We design and analyze an upwind scheme based on the symmetrized version of the equations in the non-conservative form. The scheme is shown to converge to a weak solution of the equations. Furthermore, the discrete divergence produced by the scheme is shown to be...