Publications

2011

• The role of electrode direction during axonal bipolar electrical stimulation : a bidomain computational model study
  
  • Pantz Olivier
  • Mandonnet Emmanuel
  Acta Neurochirurgica, Springer Verlag, 2011.
• Large Time-Step Numerical Scheme for the Seven-Equation Model of Compressible Two-Phase Flows
  
  • Chalons Christophe
  • Coquel Frédéric
  • Kokh Samuel
  • Spillane Nicole
  , 2011, pp.pp. 225-233. We consider the seven-equation model for compressible two-phase flows and propose a large time-step numerical scheme based on a time implicit-explicit Lagrange-Projection strategy introduced by Coquel et al. for Euler equations. The main objective is to get a Courant-Friedrichs-Lewy (CFL) condition driven by (slow) contact waves instead of (fast) acoustic waves. (10.1007/978-3-642-20671-9_24)
  DOI : 10.1007/978-3-642-20671-9_24
• Damage and fracture evolution in brittle materials by shape optimization methods
  
  • Allaire Grégoire
  • Jouve François
  • van Goethem Nicolas
  Journal of Computational Physics, Elsevier, 2011, 230 (12), pp.5010--5044. This paper is devoted to a numerical implementation of the Francfort-Marigo model of damage evolution in brittle materials. This quasi-static model is based, at each time step, on the minimization of a total energy which is the sum of an elastic energy and a Griffith-type dissipated energy. Such a minimization is carried over all geometric mixtures of the two, healthy and damaged, elastic phases, respecting an irreversibility constraint. Numerically, we consider a situation where two well-separated phases coexist, and model their interface by a level set function that is transported according to the shape derivative of the minimized total energy. In the context of interface variations (Hadamard method) and using a steepest descent algorithm, we compute local minimizers of this quasi-static damage model. Initially, the damaged zone is nucleated by using the so-called topological derivative. We show that, when the damaged phase is very weak, our numerical method is able to predict crack propagation , including kinking and branching. Several numerical examples in 2d and 3d are discussed.
• Optimal Control of the Atmospheric Reentry of a Space Shuttle by an Homotopy Method
  
  • Hermant Audrey
  Optimal Control Applications and Methods, Wiley, 2011, 32 (6), pp.627-646. This paper deals with the optimal control problem of the atmospheric reentry of a space shuttle with a second-order state constraint on the thermal flux. We solve the problem using the shooting algorithm combined with an homotopy method which automatically determines the structure of the optimal trajectory (composed of one boundary arc and one touch point). (10.1002/oca.961)
  DOI : 10.1002/oca.961
• Stable reconstruction of generalized impedance boundary conditions
  
  • Bourgeois Laurent
  • Chaulet Nicolas
  • Haddar Houssem
  Inverse Problems, IOP Publishing, 2011, 27 (9), pp.095002. (10.1088/0266-5611/27/9/095002)
  DOI : 10.1088/0266-5611/27/9/095002
• A discontinuous Galerkin solver for front propagation
  
  • Bokanowski Olivier
  • Cheng Yingda
  • Shu Chi-Wang
  SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2011, 33 (2), pp.923-938. We propose a new discontinuous Galerkin (DG) method based on [Cheng and Shu, JCP, 2007] to solve a class of Hamilton-Jacobi equations that arises from optimal control problems. These equations are connected to front propagation problems or minimal time problems with non isotropic dynamics. Several numerical experiments show the relevance of our method, in particular for front propagation.
• Swing Options Valuation:a BSDE with Constrained Jumps Approach
  
  • Bernhart Marie
  • Pham Huyên
  • Tankov Peter
  • Warin Xavier
  , 2011. We introduce a new probabilistic method for solving a class of impulse control problems based on their representations as Backward Stochastic Differential Equations (BSDEs for short) with constrained jumps. As an example, our method is used for pricing Swing options. We deal with the jump constraint by a penalization procedure and apply a discrete-time backward scheme to the resulting penalized BSDE with jumps. We study the convergence of this numerical method, with respect to the main approximation parameters: the jump intensity $\lambda$, the penalization parameter $p > 0$ and the time step. In particular, we obtain a convergence rate of the error due to penalization of order $(\lambda p)^{\alpha - \frac{1}{2}}, \forall \alpha \in \left(0, \frac{1}{2}\right)$. Combining this approach with Monte Carlo techniques, we then work out the valuation problem of (normalized) Swing options in the Black and Scholes framework. We present numerical tests and compare our results with a classical iteration method.
• Characterization of the value function of final state constrained control problems with BV trajectories
  
  • Briani Ariela
  • Zidani Hasnaa
  Communications on Pure and Applied Mathematics, Wiley, 2011, 10 (6), pp.1567-1587. This paper aims to investigate a control problem governed by differential equations with Radon measure as data and with final state constraints. By using a known reparametrization method (by Dal Maso and Rampazzo [18]), we obtain that the value function can be characterized by means of an auxiliary control problem of absolutely continuous trajectories, involving time-measurable Hamiltonian. We study the characterization of the value function of this auxiliary problem and discuss its numerical approximations. (10.3934/cpaa.2011.10.1567)
  DOI : 10.3934/cpaa.2011.10.1567
• The Electromagnetic Interior Transmission Problem for Regions with Cavities
  
  • Cossonnière Anne
  • Haddar Houssem
  SIAM Journal on Mathematical Analysis, Society for Industrial and Applied Mathematics, 2011, 43 (4), pp.1698-1715. We consider the electromagnetic interior transmission problem in the case when the medium has cavities, i.e. regions in which the index of refraction is the same as in the host medium. We address the configuration where the electromagnetic permeability is constant while the electric permittivity is variable and may be anisotropic. In this case, using appropriate reformulation of the problem into a fourth order pde, we establish the Fredholm property for this problem and show that transmission eigenvalues exist and form a discrete set. Monotonicity properties of the first eigenvalue in terms of the permittivity and the size of the cavity are established. (10.1137/100813890)
  DOI : 10.1137/100813890
• Application of convex lexicographical optimization to the balance of GRTgaz gas grid
  
  • Bonnans Joseph Frederic
  • Paraisy Ruben
  • Veyrat Sébastien
  • Adam Soizic
  Journal of Global Optimization, Springer Verlag, 2011, 49 (3), pp.415--423.