Publications

2015

• A remark on accelerated block coordinate descent for computing the proximity operators of a sum of convex functions.
  
  • Chambolle Antonin
  • Pock Thomas
  SMAI Journal of Computational Mathematics, Société de Mathématiques Appliquées et Industrielles (SMAI), 2015, 1, pp.29 - 54. We analyze alternating descent algorithms for minimizing the sum of a quadratic function and block separable non-smooth functions. In case the quadratic interactions between the blocks are pairwise, we show that the schemes can be accelerated, leading to improved convergence rates with respect to related accelerated parallel proximal descent. As an application we obtain very fast algorithms for computing the proximity operator of the 2D and 3D total variation. (10.5802/smai-jcm.3)
  DOI : 10.5802/smai-jcm.3
• Numerical analysis of the nonlinear Schrödinger equation with white noise dispersion
  
  • Belaouar Radoin
  • de Bouard Anne
  • Debussche Arnaud
  Stochastics and Partial Differential Equations: Analysis and Computations, Springer US, 2015, 3 (1), pp.103-132. This article is devoted to the numerical study of a nonlinear Schrödinger equation in which the coefficient in front of the group velocity dispersion is multiplied by a real valued Gaussian white noise. We first perform the numerical analysis of a semi-discrete Crank-Nicolson scheme in the case when the continuous equation possesses a unique global solution. We prove that the strong order of convergence in probability is equal to one in this case. In a second step, we numerically investigate, in space dimension one, the behavior of the solutions of the equation for different power nonlinearities, corresponding to subcritical, critical or supercritical nonlinearities in the deterministic case. Numerical evidence of a change in the critical power due to the presence of the noise is pointed out. (10.1007/s40072-015-0044-z)
  DOI : 10.1007/s40072-015-0044-z
• Boundary Integral Equations for the Transmission Eigenvalue Problem for Maxwell’s Equations
  
  • Cakoni Fioralba
  • Haddar Houssem
  • Meng Shixu
  Journal of Integral Equations and Applications, Rocky Mountain Mathematics Consortium, 2015, pp.26. In this paper we consider the transmission eigenvalue problem for Maxwell’s equations corresponding to non-magnetic inhomogeneities with contrast in electric permittivity that changes sign inside its support. We formulate the transmission eigenvalue problem as an equivalent homogeneous system of boundary integral equa- tion, and assuming that the contrast is constant near the boundary of the support of the inhomogeneity, we prove that the operator associated with this system is Fredholm of index zero and depends analytically on the wave number. Then we show the existence of wave numbers that are not transmission eigenvalues which by an application of the analytic Fredholm theory implies that the set of transmission eigenvalues is discrete with positive infinity as the only accumulation point.
• Axisymmetric eddy current inspection of highly conducting thin layers via asymptotic models
  
  • Haddar Houssem
  • Jiang Zixian
  Inverse Problems, IOP Publishing, 2015. Thin copper deposits covering the steam generator tubes can blind eddy current probes in non-destructive testings of problematic faults and are therefore important to be identified. Existing methods based on shape reconstruction using eddy current signals encounter difficulties of high numerical costs due to the layer's small thickness and high conductivity. In this article, we approximate the axisymmetric eddy current problem with some appropriate asymptotic models using effective transmission conditions representing the thin deposits. In these models, the geometrical information related to the deposit is transformed into parameter coefficients on a fictitious interface. Standard iterative inversion algorithm is then applied to the asymptotic models in order to reconstruct the thickness of the thin copper layers. Numerical tests both validating the asymptotic model and benefit of the inversion procedure are provided.
• On the scaling limits of Galton–Watson processes in varying environments
  
  • Bansaye Vincent
  • Simatos Florian
  Electronic Journal of Probability, Institute of Mathematical Statistics (IMS), 2015, 20 (75), pp.1-36. We establish a general sufficient condition for a sequence of Galton–Watson branching processes in varying environments to converge weakly. This condition extends previ- ous results by allowing offspring distributions to have infinite variance. Our assumptions are stated in terms of pointwise convergence of a triplet of two real- valued functions and a measure. The limiting process is characterized by a backwards integro-differential equation satisfied by its Laplace exponent, which generalizes the branching equation satisfied by continuous state branching processes. Several examples are discussed, namely branching processes in random environment, Feller diffusion in varying environments and branching processes with catastrophes. (10.1214/EJP.v20-3812)
  DOI : 10.1214/EJP.v20-3812
• Converse Lyapunov-Krasovskii theorems for uncertain retarded differential equations
  
  • Haidar Ihab
  • Mason Paolo
  • Sigalotti Mario
  Automatica, Elsevier, 2015, 62, pp.263-273. In this article we give a collection of converse Lyapunov–Krasovskii theorems for uncertain retarded differential equations. We show that the existence of a weakly-degenerate Lyapunov–Krasovskii functional is a necessary and sufficient condition for the global exponential stability of linear retarded functional differential equations. This is carried out using a switched system representation approach. (10.1016/j.automatica.2015.09.034)
  DOI : 10.1016/j.automatica.2015.09.034
• Definable Zero-Sum Stochastic Games
  
  • Bolte Jérôme
  • Gaubert Stéphane
  • Vigeral Guillaume
  Mathematics of Operations Research, INFORMS, 2015, 40 (1), pp.171-191. Definable zero-sum stochastic games involve a finite number of states and action sets, reward and transition functions that are definable in an o-minimal structure. Prominent examples of such games are finite, semi-algebraic or globally subanalytic stochastic games. We prove that the Shapley operator of any definable stochastic game with separable transition and reward functions is definable in the same structure. Definability in the same structure does not hold systematically: we provide a counterexample of a stochastic game with semi-algebraic data yielding a non semi-algebraic but globally subanalytic Shapley operator. %Showing the definability of the Shapley operator in full generality appears thus as a complex and challenging issue. } Our definability results on Shapley operators are used to prove that any separable definable game has a uniform value; in the case of polynomially bounded structures we also provide convergence rates. Using an approximation procedure, we actually establish that general zero-sum games with separable definable transition functions have a uniform value. These results highlight the key role played by the tame structure of transition functions. As particular cases of our main results, we obtain that stochastic games with polynomial transitions, definable games with finite actions on one side, definable games with perfect information or switching controls have a uniform value. Applications to nonlinear maps arising in risk sensitive control and Perron-Frobenius theory are also given. (10.1287/moor.2014.0666)
  DOI : 10.1287/moor.2014.0666
• Dispersal is a major driver of the latitudinal diversity gradient of Carnivora
  
  • Rolland Jonathan
  • Condamine Fabien L.
  • Beeravolu Reddy Champak
  • Jiguet Frédéric
  • Morlon Hélène
  Global Ecology and Biogeography, Wiley, 2015, 24 (9), pp.1059 - 1071. <strong>Aim</strong> Understanding the relative contribution of diversification rates (speciation and extinction) and dispersal in the formation of the latitudinal diversity gradient - the decrease in species richness with increasing latitude - is a main goal of biogeography. The mammalian order Carnivora, which comprises 286 species, displays the traditional latitudinal diversity gradient seen in almost all mammalian orders. Yet the processes driving high species richness in the tropics may be fundamentally different in this group from that in other mammalian groups. Indeed, a recent study suggested that in Carnivora, unlike in all other major mammalian orders, net diversification rates are not higher in the tropics than in temperate regions. Our goal was thus to understand the reasons why there are more species of Carnivora in the tropics. <strong>Location</strong> World-wide. <strong>Methods</strong> We reconstructed the biogeographical history of Carnivora using a time-calibrated phylogeny of the clade comprising all terrestrial species and dispersal-extinction-cladogenesis models. We also analysed a fossil dataset of carnivoran genera to examine how the latitudinal distribution of Carnivora varied through time. <strong>Results</strong> Our biogeographical analyses suggest that Carnivora originated in the East Palaearctic (i.e. Central Asia, China) in the early Palaeogene. Multiple independent lineages dispersed to low latitudes following three main paths: toward Africa, toward India/Southeast Asia and toward South America via the Bering Strait. These dispersal events were probably associated with local extinctions at high latitudes. Fossil data corroborate a high-latitude origin of the group, followed by late dispersal events toward lower latitudes in the Neogene. <strong>Main conclusions</strong> Unlike most other mammalian orders, which originated and diversified at low latitudes and dispersed out of the tropics', Carnivora originated at high latitudes, and subsequently dispersed southward. Our study provides an example of combining phylogenetic and fossil data to understand the generation and maintenance of global-scale geographical variations in species richness. (10.1111/geb.12354)
  DOI : 10.1111/geb.12354
• Controllability of spin-boson systems
  
  • Boscain Ugo
  • Mason Paolo
  • Panati Gianluca
  • Sigalotti Mario
  Journal of Mathematical Physics, American Institute of Physics (AIP), 2015, 56. In this paper we study the so-called spin-boson system, namely {a two-level system} in interaction with a distinguished mode of a quantized bosonic field. We give a brief description of the controlled Rabi and Jaynes--Cummings models and we discuss their appearance in the mathematics and physics literature. We then study the controllability of the Rabi model when the control is an external field acting on the bosonic part. Applying geometric control techniques to the Galerkin approximation and using perturbation theory to guarantee non-resonance of the spectrum of the drift operator, we prove approximate controllability of the system, for almost every value of the interaction parameter.
• Configuration Tracking for Mechanical Systems by Kinematic Reduction and Fast Oscillating Controls
  
  • Barbero-Liñán M.
  • Sigalotti Mario
  , 2015.
• Second order mean field games with degenerate diffusion and local coupling
  
  • Cardaliaguet Pierre
  • Graber J.
  • Porretta Alessio
  • Tonon Daniela
  Nonlinear Differential Equations and Applications, Springer Verlag, 2015, 22 (5), pp.1287-1317. We analyze a (possibly degenerate) second order mean field games system of partial differential equations. The distinguishing features of the model considered are (1) that it is not uniformly parabolic, including the first order case as a possibility, and (2) the coupling is a local operator on the density. As a result we look for weak, not smooth, solutions. Our main result is the existence and uniqueness of suitably defined weak solutions, which are characterized as minimizers of two optimal control problems. We also show that such solutions are stable with respect to the data, so that in particular the degenerate case can be approximated by a uniformly parabolic (viscous) perturbation.
• A Priori Error Estimate of a Multiscale Finite Element Method for Transport Modeling
  
  • Ouaki Franck
  • Allaire Grégoire
  • Desroziers Sylvain
  • Enchéry Guillaume
  SeMA Journal: Boletin de la Sociedad Española de Matemática Aplicada, Springer, 2015, 67 (1), pp.1-37. This work proposes an \textit{a priori} error estimate of a multiscale finite element method to solve convection-diffusion problems where both velocity and diffusion coefficient exhibit strong variations at a scale which is much smaller than the domain of resolution. In that case, classical discretization methods, used at the scale of the heterogeneities, turn out to be too costly. Our method, introduced in~\cite{ECCOMAS}, aims at solving this kind of problems on coarser grids with respect to the size of the heterogeneities by means of particular basis functions. These basis functions are defined using cell problems and are designed to reproduce the variations of the solution on an underlying fine grid. Since all cell problems are independent from each other, these problems can be solved in parallel, which makes the method very efficient when used on parallel architectures. This article focuses on the proof of an \textit{a priori} error estimate of this method.
• Mathematical Methods in Elasticity Imaging
  
  • Ammari Habib
  • Bretin Elie
  • Garnier Josselin
  • Kang Hyeonbae
  • Lee Hyundae
  • Wahab Abdul
  , 2015, pp.240 p..
• Hawkes Processes in Finance
  
  • Bacry Emmanuel
  • Mastromatteo Iacopo
  • Muzy Jean-François
  Market microstructure and liquidity, World scientific publishing company, 2015, 01 (01), pp.1550005. no abstract (10.1142/S2382626615500057)
  DOI : 10.1142/S2382626615500057
• Singular limits for reaction-diffusion equations with fractional Laplacian and local or nonlocal nonlinearity
  
  • Méléard Sylvie
  • Mirrahimi Sepideh
  Communications in Partial Differential Equations, Taylor & Francis, 2015, 40 (5), pp.957-993. We perform an asymptotic analysis of models of population dynamics with a fractional Laplacian and local or nonlocal reaction terms. The first part of the paper is devoted to the long time/long range rescaling of the fractional Fisher-KPP equation. This rescaling is based on the exponential speed of propagation of the population. In particular we show that the only role of the fractional Laplacian in determining this speed is at the initial layer where it determines the thickness of the tails of the solutions. Next, we show that such rescaling is also possible for models with non-local reaction terms, as selection-mutation models. However, to obtain a more relevant qualitative behavior for this second case, we introduce, in the second part of the paper, a second rescaling where we assume that the diffusion steps are small. In this way, using a WKB ansatz, we obtain a Hamilton-Jacobi equation in the limit which describes the asymptotic dynamics of the solutions, similarly to the case of selection-mutation models with a classical Laplace term or an integral kernel with thin tails. However, the rescaling introduced here is very different from the latter cases. We extend these results to the multidimensional case. (10.1080/03605302.2014.963606)
  DOI : 10.1080/03605302.2014.963606
• MatVPC: A User-Friendly MATLAB-Based Tool for the Simulation and Evaluation of Systems Pharmacology Models
  
  • Biliouris Kostas
  • Lavielle Marc
  • Trame Mirjam
  CPT: Pharmacometrics and Systems Pharmacology, American Society for Clinical Pharmacology and Therapeutics ; International Society of Pharmacometrics, 2015. Quantitative systems pharmacology (QSP) models are progressively entering the arena of contemporary pharmacology. The efficient implementation and evaluation of complex QSP models necessitates the development of flexible computational tools that are built into QSP mainstream software. To this end, we present MatVPC, a versatile MATLAB-based tool that accommodates QSP models of any complexity level. MatVPC executes Monte Carlo simulations as well as automatic construction of visual predictive checks (VPCs) and quantified VPCs (QVPCs). VPC is a model diagnostic tool that facilitates the evaluation of both the structural and the stochastic part of a model. It is constructed by superimposing the observations over the model simulations while accounting for both the interindivid-ual variability as well as the residual variability. 1 Once underutilized, 2 the VPC now is recognized as one of the most valuable model diagnostics in pharmacological model evaluation. 3–5 Its superiority over comparable diagnostic tools has been established 6 and reflected by the fact that regulatory agencies recommend it as one of the central model diagnostics. 7 (10.1002/psp4.12011)
  DOI : 10.1002/psp4.12011
• Ensuring robustness of domain decomposition methods by block strategies
  
  • Gosselet Pierre
  • Rixen Daniel
  • Spillane Nicole
  • Roux François-Xavier
  , 2015. no abstract
• Coupling techniques for nonlinear hyperbolic equations. IV. Well-balanced schemes for scalar multi-dimensional and multi-component laws
  
  • Boutin Benjamin
  • Coquel Frédéric
  • LeFloch Philippe G.
  Mathematics of Computation, American Mathematical Society, 2015, 84 (294), pp.1663-1702. This series of papers is devoted to the formulation and the approximation of coupling problems for nonlinear hyperbolic equations. The coupling across an interface in the physical space is formulated in term of an augmented system of partial differential equations. In an earlier work, this strategy allowed us to develop a regularization method based on a thick interface model in one space variable for coupling scalar equations. In the present paper, we significantly extend this framework and, in addition, encompass equations in several space variables. This new formulation includes the coupling of several distinct scalar conservation laws and allows for a possible covering in space. Our main contributions are, on one hand, the design and analysis of a well–balanced finite volume method on general triangulations and, on the other hand, a proof of convergence of this method toward entropy solutions, extending Coquel, Cockburn, and LeFloch's theory (restricted to a single conservation law without coupling). The core of our analysis is, first, the derivation of entropy inequalities as well as a discrete entropy dissipation estimate and, second, a proof of convergence toward the entropy solution of the coupling problem. (10.1090/S0025-5718-2015-02933-0)
  DOI : 10.1090/S0025-5718-2015-02933-0
• Developmental Partial Differential Equations
  
  • Pouradier Duteil Nastassia
  • Rossi Francesco
  • Boscain Ugo
  • Piccoli Benedetto
  , 2015. In this paper, we introduce the concept of Developmental Partial Differential Equation (DPDE), which consists of a Partial Differential Equation (PDE) on a time-varying manifold with complete coupling between the PDE and the manifold’s evolution. In other words, the manifold’s evolution depends on the solution to the PDE, and vice versa the differential operator of the PDE depends on the manifold’s geometry. DPDE is used to study a diffusion equation with source on a growing surface whose growth depends on the intensity of the diffused quantity. The surface may, for instance, represent the membrane of an egg chamber and the diffused quantity a protein activating a signaling pathway leading to growth. Our main objective is to show controllability of the surface shape using a fixed source with variable intensity for the diffusion. More specifically, we look for a control driving a symmetric manifold shape to any other symmetric shape in a given time interval. For the diffusion we take directly the Laplace-Beltrami operator of the surface, while the surface growth is assumed to be equal to the value of the diffused quantity. We introduce a theoretical framework, provide approximate controllability and show numerical results. Future applications include a specific model for the oogenesis of Drosophila melanogaster.
• Phaseless inverse scattering in the one-dimensional case
  
  • Novikov Roman
  Eurasian Journal of Mathematical and Computer Applications, Eurasian National University, Kazakhstan (Nur-Sultan), 2015, 3 (1), pp.64-70. We consider the one-dimensional Schrödinger equation with a potential satisfying the standard assumptions of the inverse scattering theory and supported on the half-line x ≥ 0. For this equation at fixed positive energy we give explicit formulas for finding the full complex valued reflection coefficient to the left from appropriate phaseless scattering data measured on the left, i.e. for x < 0. Using these formulas and known inverse scattering results we obtain global uniqueness and reconstruction results for phaseless inverse scattering in dimension d = 1.
• Artificial boundary conditions for axisymmetric eddy current probe problems
  
  • Haddar Houssem
  • Jiang Zixian
  • Lechleiter Armin
  Computers & Mathematics with Applications, Elsevier, 2015, 68 (12, Part A,), pp.1844–1870. We study different strategies for the truncation of computational domains in the simulation of eddy current probes of elongated axisymmetric tubes. For axial fictitious boundaries, an exact Dirichlet-to-Neumann map is proposed and mathematically analyzed via a non-selfadjoint spectral problem: under general assumptions we show convergence of the solution to an eddy current problem involving a truncated Dirichlet-to-Neumann map to the solution on the entire, unbounded axisymmetric domain as the truncation parameter tends to infinity. Under stronger assumptions on the physical parameters of the eddy current problem, convergence rates are shown. We further validate our theoretical results through numerical experiments for a realistic physical setting inspired by eddy current probes of nuclear reactor core tubes. (10.1016/j.camwa.2014.10.008)
  DOI : 10.1016/j.camwa.2014.10.008
• Lyapunov and Minimum-Time Path Planning for Drones
  
  • Maillot Thibault
  • Boscain Ugo
  • Gauthier Jean-Paul
  • Serres Ulysse
  Journal of Dynamical and Control Systems, Springer Verlag, 2015, 21 (1), pp.1-34. (10.1007/s10883-014-9222-y)
  DOI : 10.1007/s10883-014-9222-y
• Inverse scattering without phase information
  
  • Novikov Roman
  Séminaire Laurent Schwartz - EDP et applications, Centre de mathématiques Laurent Schwartz, 2015, 2014-2015, pp.Exposé no. 16, 13 pp. We report on nonuniqueness, uniqueness and reconstruction results in quantum mechanical and acoustic inverse scattering without phase information. We are motivated by recent and very essential progress in this domain. This paper is an extended version of the talk given at Séminaire Laurent Schwartz on March 31, 2015. (10.5802/slsedp.74)
  DOI : 10.5802/slsedp.74
• The topological derivative of stress-based cost functionals in anisotropic elasticity
  
  • Delgado Gabriel
  • Bonnet Marc
  Computers & Mathematics with Applications, Elsevier, 2015, 69, pp.1144-1166. The topological derivative of cost functionals J that depend on the stress (through the displacement gradient, assuming a linearly elastic material behavior) is considered in a quite general 3D setting where both the background and the inhomogeneity may have arbitrary anisotropic elastic properties. The topological derivative dJ(z) of J quantifies the asymptotic behavior of J under the nucleation in the background elastic medium of a small anisotropic inhomogeneity of characteristic radius a at a specified location z. The fact that the strain perturbation inside an elastic inhomogeneity remains finite for arbitrarily small a makes the small-inhomogeneity asymptotics of stress-based cost functionals quite different than that of the more usual displacement-based functionals. The asymptotic perturbation of J is shown to be of order O(a^3) for a wide class of stress-based cost functionals having smooth densities. The topological derivative of J, i.e. the coefficient of the O(a^3) perturbation, is established, and computational procedures then discussed. The resulting small-inhomogeneity expansion of J is mathematically justified (i.e. its remainder is proved to be of order o(a^3)). Several 2D and 3D numerical examples are presented, in particular demonstrating the proposed formulation of \dJ on cases involving anisotropic elasticity and non-quadratic cost functionals. (10.1016/j.camwa.2015.03.010)
  DOI : 10.1016/j.camwa.2015.03.010
• Energy release rate for non smooth cracks in planar elasticity
  
  • Babadjian Jean-François
  • Chambolle Antonin
  • Lemenant Antoine
  Journal de l'École polytechnique — Mathématiques, École polytechnique, 2015, 2, pp.117-152. This paper is devoted to the characterization of the energy release rate of a crack which is merely closed, connected, and with density $1/2$ at the tip. First, the blow-up limit of the displacement is analyzed, and the convergence to the corresponding positively $1/2$-homogenous function in the cracked plane is established. Then, the energy release rate is obtained as the derivative of the elastic energy with respect to an infinitesimal additional crack increment.