In this paper, we study a model for the magnetization in thin ferromagnetic films. It comes as a variational problem for $S^1$-valued maps $m^{\text{'}}$ (the magnetization) of two variables $x^{\text{'}}$: $E_{\epsilon }\left(m^{\text{'}}\right)=\epsilon \int {|{\nabla }^{\text{'}}·m^{\text{'}}|}^{2}d{x}^{\text{'}}+\frac{1}{2}\int {\left||{\nabla }^{\text{'}}{|}^{-1/2}{\nabla }^{\text{'}}·m^{\text{'}}\right|}^{2}d{x}^{\text{'}}$. We are interested in the behavior of minimizers as $\epsilon \to 0$. They are expected to be $S^1$-valued maps $m^{\text{'}}$ of vanishing distributional divergence ${\nabla }^{\text{'}}·m^{\text{'}}=0$, so that appropriate boundary conditions enforce line discontinuities. For finite $\epsilon >0$, these line discontinuities are approximated by smooth transition layers, the so-called Néel walls. Néel...

We consider the numerical approximation of a first order stationary hyperbolic equation by the method of characteristics with pseudo time step k using discontinuous finite elements on a mesh ${𝒯}_h$. For this method, we exhibit a “natural” norm || ||h,k for which we show that the discrete variational problem $P_h^k$ is well posed and we obtain an error estimate. We show that when k goes to zero problem $\left(P_h^k\right)$ (resp. the || ||h,k norm) has as a limit problem (Ph) (resp. the || ||h norm) associated to the...

We discuss the formulation of a simulator in three spatial dimensions for a multicomponent, two phase (air, water) system of groundwater flow and transport with biodegradation kinetics and wells with multiple screens. The simulator has been developed for parallel, distributed memory, message passing machines. The numerical procedures employed are a fully implicit expanded mixed finite element method for flow and either a characteristics-mixed method or a Godunov method for transport and reactions...

In this article we introduce an adaptive multi-level method in space and time for convection diffusion problems. The scheme is based on a multi-level spatial splitting and the use of different time-steps. The temporal discretization relies on the characteristics method. We derive an a posteriori error estimate and design a corresponding adaptive algorithm. The efficiency of the multi-level method is illustrated by numerical experiments, in particular for a convection-dominated problem.

We describe a numerical method for the equation $u_t+p{u}_x-\epsilon u_{xx}=f$ in $(0,1)\times (0,T)$ with Dirichlet boundary and initial conditions which is a combination of the method of characteristics and the finite-difference method. We prove both an a priori local error-estimate of a high order and stability. Example 3.3 indicates that our approximate solutions are disturbed only by a minimal amount of the artificial diffusion.

In this paper, a sub-optimal boundary control strategy for a free boundary problem is investigated. The model is described by a non-smooth convection-diffusion equation. The control problem is addressed by an instantaneous strategy based on the characteristics method. The resulting time independent control problems are formulated as function space optimization problems with complementarity constraints. At each time step, the existence of an optimal solution is proved and first-order optimality conditions...

We present the fiber-spring elastic model of the arterial wall with atherosclerotic plaque composed of a lipid pool and a fibrous cap. This model allows us to reproduce pressure to cross-sectional area relationship along the diseased vessel which is used in the network model of global blood circulation. Atherosclerosis attacks a region of systemic arterial network. Our approach allows us to examine the impact of the diseased region onto global haemodynamics....

The modelling and the numerical resolution of the electrical charging of a spacecraft in interaction with the Earth magnetosphere is considered. It involves the Vlasov-Poisson system, endowed with non standard boundary conditions. We discuss the pros and cons of several numerical methods for solving this system, using as benchmark a simple 1D model which exhibits the main difficulties of the original models.

Numerical approximation schemes are discussed for the solution of contaminant transport with adsorption in dual-well flow. The method is based on time stepping and operator splitting for the transport with adsorption and diffusion. The nonlinear transport is solved by Godunov’s method. The nonlinear diffusion is solved by a finite volume method and by Newton’s type of linearization. The efficiency of the method is discussed.

We propose a numerical scheme to compute the motion of a two-dimensional rigid body in a viscous fluid. Our method combines the method of characteristics with a finite element approximation to solve an ALE formulation of the problem. We derive error estimates implying the convergence of the scheme.

The Maxwell type elastic-plastic solids are characterized by decaying the absolute values of the principal components of the deviatoric part of the stress tensor during the plastic relaxation step. We propose a mathematical formulation of such a model which is compatible with the von Mises criterion of plasticity. Numerical examples show the ability of the model to deal with complex physical phenomena.