Publications

2018

• Résolution des équations de Navier-Stokes linéarisées pour l'aéroélasticité, l’optimisation de forme et l’aéroacoustique
  
  • Bissuel Aloïs
  , 2018. Les équations de Navier-Stokes linéarisées sont utilisées dans l’industrie aéronautique pour l’optimisation de forme aérodynamique, l’aéroélasticité et l’aéroacoustique. Deux axes ont été suivis pour accélérer et rendre plus robuste la résolution de ces équations. Le premier est l’amélioration de la méthode itérative de résolution de systèmes linéaires utilisée, et le deuxième la formulation du schéma numérique conduisant à ce système linéaire. Dans cette première partie, l’extension de l’algorithme GMRES avec déflation spectrale à des systèmes à plusieurs seconds membres a été testée sur des cas tests industriels. L’amélioration du préconditionnement de la méthode GMRES par l’utilisation d'une méthode de Schwarz additive avec préconditionneur ILU(k) a permis une accélération du temps de résolution allant jusqu’à un facteur dix, ainsi que la convergence de cas jusqu’alors impossibles à résoudre. La deuxième partie présente d’abord un travail sur la stabilisation SUPG du schéma élément fini utilisé. La forme proposée de la matrice de stabilisation, dite complète, a donné des résultats encourageants en non-linéaire qui ne se sont pas transposés en linéarisé. Une étude sur les conditions aux limites de Dirichlet clôt cette partie. Une méthode algébrique d’imposition de conditions non homogènes sur des variables non triviales du calcul, qui a permis l’application industrielle à l’aéroacoustique, y est détaillée. De plus, la preuve est apportée que le caractère transparent d’une condition de Dirichlet homogène sur toutes les variables s’explique par le schéma SUPG.
• Principal-Agent Problem with Common Agency without Communication
  
  • Mastrolia Thibaut
  • Ren Zhenjie
  , 2018. In this paper, we consider a problem of contract theory in which several Principals hire a common Agent and we study the model in the continuous time setting. We show that optimal contracts should satisfy some equilibrium conditions and we reduce the optimisation problem of the Principals to a system of coupled Hamilton-Jacobi-Bellman (HJB) equations. We provide conditions ensuring that for risk-neutral Principals, the system of coupled HJB equations admits a solution. Further, we apply our study in a more specific linear-quadratic model where two interacting Principals hire one common Agent. In this continuous time model, we extend the result of Bernheim and Whinston (1986) in which the authors compare the optimal effort of the Agent in a non-cooperative Principals model and that in the aggregate model, by showing that these two optimisations coincide only in the first best case. We also study the sensibility of the optimal effort and the optimal remunerations with respect to appetence parameters and the correlation between the projects. (10.1137/17M1133609)
  DOI : 10.1137/17M1133609
• Development and performance of npde for the evaluation of time-to-event models
  
  • Cerou Marc
  • Lavielle Marc
  • Brendel Karl
  • Chenel Marylore
  • Comets Emmanuelle
  Pharmaceutical Research, American Association of Pharmaceutical Scientists, 2018, 35 (2), pp.30. Purpose - Normalised prediction distribution errors (npde) are used to graphically and statistically evaluate mixed-effect models for continuous responses. In this study, our aim was to extend npde to time-to-event (TTE) models and evaluate their performance. Methods - Let V denote a dataset with censored TTE observations. The null hypothesis (H) is that observations in V can be described by model M. We extended npde to TTE models using imputations to take into account censoring. We then evaluated their performance in terms of type I error and power to detect model misspecifications for TTE data by means of a simulation study with different sample sizes. Results - Type I error was found to be close to the expected 5% significance level for all sample sizes tested. The npde were able to detect misspecifications in the baseline hazard as well as in the link between the longitudinal variable and the survival function. The ability to detect model misspecifications increased as the difference in the shape of the survival function became more apparent. As expected, the power also increased as the sample size increased. Imputing the censored events tended to decrease the percentage of rejections. Conclusions - We have shown that npde can be readily extended to TTE data and that they perform well with an adequate type I error. (10.1007/s11095-017-2291-3)
  DOI : 10.1007/s11095-017-2291-3
• Quantitative estimates for the flux of TASEP with dilute site disorder
  
  • Bahadoran Christophe
  • Bodineau Thierry
  Electronic Journal of Probability, Institute of Mathematical Statistics (IMS), 2018, 23. We prove that the flux function of the totally asymmetric simple exclusion process (TASEP) with site disorder exhibits a flat segment for sufficiently dilute disorder. For high dilution, we obtain an accurate description of the flux. The result is established under a decay assumption of the maximum current in finite boxes, which is implied in particular by a sufficiently slow power tail assumption on the disorder distribution near its minimum. To circumvent the absence of explicit invariant measures, we use an original renormalization procedure and some ideas inspired by homogenization. (10.1214/18-EJP137)
  DOI : 10.1214/18-EJP137
• Avis en réponse à la saisine HCB - dossier EFSA-GMO-RX009. Paris, le 4 juin 2018
  
  • Comité Scientifique Du Haut Conseil Des Biotechnologies .
  • Angevin Frédérique
  • Bagnis Claude
  • Bar-Hen Avner
  • Barny Marie-Anne
  • Boireau Pascal
  • Brévault Thierry
  • Chauvel Bruno B.
  • Collonnier Cécile
  • Couvet Denis
  • Dassa Elie
  • de Verneuil Hubert
  • Demeneix Barbara
  • Franche Claudine
  • Guerche Philippe
  • Guillemain Joël
  • Hernandez Raquet Guillermina
  • Khalife Jamal
  • Klonjkowski Bernard
  • Lavielle Marc
  • Le Corre Valérie
  • Lefèvre François
  • Lemaire Olivier
  • Lereclus Didier D.
  • Maximilien Rémy
  • Meurs Eliane
  • Naffakh Nadia
  • Négre Didier
  • Noyer Jean-Louis
  • Ochatt Sergio
  • Pages Jean-Christophe
  • Raynaud Xavier
  • Regnault-Roger Catherine
  • Renard Michel M.
  • Renault Tristan
  • Saindrenan Patrick
  • Simonet Pascal
  • Troadec Marie-Bérengère
  • Vaissière Bernard
  • Vilotte Jean-Luc
  , 2018.
• Log-barrier interior point methods are not strongly polynomial
  
  • Allamigeon Xavier
  • Benchimol Pascal
  • Gaubert Stéphane
  • Joswig Michael
  SIAM Journal on Applied Algebra and Geometry, Society for Industrial and Applied Mathematics, 2018, 2 (1), pp.140-178. We prove that primal-dual log-barrier interior point methods are not strongly polynomial, by constructing a family of linear programs with $3r+1$ inequalities in dimension $2r$ for which the number of iterations performed is in $\Omega(2^r)$. The total curvature of the central path of these linear programs is also exponential in $r$, disproving a continuous analogue of the Hirsch conjecture proposed by Deza, Terlaky and Zinchenko. Our method is to tropicalize the central path in linear programming. The tropical central path is the piecewise-linear limit of the central paths of parameterized families of classical linear programs viewed through logarithmic glasses. This allows us to provide combinatorial lower bounds for the number of iterations and the total curvature, in a general setting. (10.1137/17M1142132)
  DOI : 10.1137/17M1142132
• SEME 2017 : identification de véhicules en utilisant le numéro VIN
  
  • Besson Rémi
  • Etchegaray Christèle
  • Ferrari Luca
  • Nordmann Samuel
  , 2018.
• Darboux–Moutard transformations and Poincare–Steklov operators
  
  • Novikov Roman
  • Taimanov Iskander
  Proceedings of the Steklov Institute of Mathematics, MAIK Nauka/Interperiodica, 2018, 302, pp.315–324. Formulas relating Poincare–Steklov operators for Schrödinger equations related by Darboux–Moutard transformations are derived. They can be used for testing algorithms of reconstruction of the potential from measurements at the boundary. (10.1134/S0081543818060160)
  DOI : 10.1134/S0081543818060160
• FEM and BEM simulations with the Gypsilab framework
  
  • Alouges François
  • Aussal Matthieu
  SMAI Journal of Computational Mathematics, Société de Mathématiques Appliquées et Industrielles (SMAI), 2018, 4, pp.297-318. (10.5802/smai-jcm.36)
  DOI : 10.5802/smai-jcm.36
• Infinite Horizon Stochastic Optimal Control Problems with Running Maximum Cost
  
  • Kröner Axel
  • Picarelli Athena
  • Zidani Hasnaa
  SIAM Journal on Control and Optimization, Society for Industrial and Applied Mathematics, 2018, 56 (5), pp.3296-3319. An infinite horizon stochastic optimal control problem with running maximum cost is considered. The value function is characterized as the viscosity solution of a second-order Hamilton-Jacobi-Bellman (HJB) equation with mixed boundary condition. A general numerical scheme is proposed and convergence is established under the assumptions of consistency, monotonicity and stability of the scheme. These properties are verified for a specific semi-Lagrangian scheme. (10.1137/17M115253X)
  DOI : 10.1137/17M115253X
• Derivation of an ornstein-uhlenbeck process for a massive particle in a rarified gas of particles
  
  • Bodineau Thierry
  • Gallagher Isabelle
  • Saint-Raymond Laure
  Annales de l'Institut Henri Poincaré (A). Physique Theorique, Birkhäuser, 2018, 19 (6). We consider the statistical motion of a convex rigid body in a gas of N smaller (spherical) atoms close to thermodynamic equilibrium. Because the rigid body is much bigger and heavier, it undergoes a lot of collisions leading to small deflections. We prove that its velocity is described, in a suitable limit, by an Ornstein-Uhlenbeck process. The strategy of proof relies on Lanford's arguments [17] together with the pruning procedure from [3] to reach diffusive times, much larger than the mean free time. Furthermore, we need to introduce a modified dynamics to avoid pathological collisions of atoms with the rigid body: these collisions, due to the geometry of the rigid body, require developing a new type of trajectory analysis. (10.1007/s00023-018-0674-6)
  DOI : 10.1007/s00023-018-0674-6
• An integrate-and-fire model to generate spike trains with long memory
  
  • Richard Alexandre
  • Orio Patricio
  • Tanré Etienne
  Journal of Computational Neuroscience, Springer Verlag, 2018. Long-range dependence (LRD) has been observed in a variety of phenomena in nature, and for several years also in the spiking activity of neurons. Often, this is interpreted as originating from a non-Markovian system. Here we show that a purely Markovian integrate-and-re (IF) model, with a noisy slow adaptation term, can generate data that appears as having LRD with a Hurst exponent (H) greater than 0.5. A proper analysis shows that the asymptotic value of H is 0.5 if a long enough sequence of events is taken into account. For comparison, we also consider a new model of individual IF neuron with fractional noise. The correlations of its spike trains are studied and proved to have long memory, unlike classical IF models. On the other hand, to correctly measure long-range dependence, it is usually necessary to know if the data are stationary. Thus, a methodology to evaluate stationarity of the interspike intervals (ISIs) is presented and applied to the various IF models. In conclusion, the spike trains of our fractional model have the long-range dependence property, while those from classical Markovian models do not. However, Markovian IF models may seem to have it because of apparent non-stationarities. (10.1007/s10827-018-0680-1)
  DOI : 10.1007/s10827-018-0680-1
• Optimizing supports for additive manufacturing
  
  • Allaire Grégoire
  • Bogosel Beniamin
  Structural and Multidisciplinary Optimization, Springer Verlag, 2018, 58 (6), pp.2493-2515. In additive manufacturing process support structures are often required to ensure the quality of the final built part. In this article we present mathematical models and their numerical implementations in an optimization loop, which allow us to design optimal support structures. Our models are derived with the requirement that they should be as simple as possible, computationally cheap and yet based on a realistic physical modeling. Supports are optimized with respect to two different physical properties. First, they must support overhanging regions of the structure for improving the stiffness of the supported structure during the building process. Second, supports can help in channeling the heat flux produced by the source term (typically a laser beam) and thus improving the cooling down of the structure during the fabrication process. Of course, more involved constraints or manufacturability conditions could be taken into account, most notably removal of supports. Our work is just a first step, proposing a general framework for support optimization. Our optimization algorithm is based on the level set method and on the computation of shape derivatives by the Hadamard method. In a first approach, only the shape and topology of the supports are optimized, for a given and fixed structure. In second and more elaborated strategy, both the supports and the structure are optimized, which amounts to a specific multiphase optimization problem. Numerical examples are given in 2-d and 3-d.
• New interior transmission problem applied to a single Floquet–Bloch mode imaging of local perturbations in periodic media
  
  • Cakoni Fioralba
  • Haddar Houssem
  • Nguyen Thi-Phong
  Inverse Problems, IOP Publishing, 2018, 35 (1), pp.015009.
• Optimization of dispersive coefficients in the homogenization of the wave equation in periodic structures
  
  • Allaire Grégoire
  • Yamada T
  Numerische Mathematik, Springer Verlag, 2018, 140 (2), pp.265-326. We study dispersive effects of wave propagation in periodic media, which can be modelled by adding a fourth-order term in the homogenized equation. The corresponding fourth-order dispersive tensor is called Burnett tensor and we numerically optimize its values in order to minimize or maximize dispersion. More precisely, we consider the case of a two-phase composite medium with an 8-fold symmetry assumption of the periodicity cell in two space dimensions. We obtain upper and lower bound for the dispersive properties, along with optimal microgeometries.
• Study of new rare event simulation schemes and their application to extreme scenario generation
  
  • Agarwal Ankush
  • de Marco Stefano
  • Gobet Emmanuel
  • Liu Gang
  Mathematics and Computers in Simulation, Elsevier, 2018, 143, pp.89-98. This is a companion paper based on our previous work [ADGL15] on rare event simulation methods. In this paper, we provide an alternative proof for the ergodicity of shaking transformation in the Gaussian case and propose two variants of the existing methods with comparisons of numerical performance. In numerical tests, we also illustrate the idea of extreme scenario generation based on the convergence of marginal distributions of the underlying Markov chains and show the impact of the discretization of continuous time models on rare event probability estimation. (10.1016/j.matcom.2017.05.004)
  DOI : 10.1016/j.matcom.2017.05.004
• Intrinsic random walks in Riemannian and sub-Riemannian geometry via volume sampling
  
  • Agrachev Andrei
  • Boscain Ugo
  • Neel Robert
  • Rizzi Luca
  ESAIM: Control, Optimisation and Calculus of Variations, EDP Sciences, 2018, 24 (3), pp.1075–1105. We relate some basic constructions of stochastic analysis to differential geometry , via random walk approximations. We consider walks on both Riemannian and sub-Riemannian manifolds in which the steps consist of travel along either geodesics or integral curves associated to orthonormal frames, and we give particular attention to walks where the choice of step is influenced by a volume on the manifold. A primary motivation is to explore how one can pass, in the parabolic scaling limit, from geodesics, orthonormal frames, and/or volumes to diffusions, and hence their infinitesimal generators , on sub-Riemannian manifolds, which is interesting in light of the fact that there is no completely canonical notion of sub-Laplacian on a general sub-Riemannian mani-fold. However, even in the Riemannian case, this random walk approach illuminates the geometric significance of Ito and Stratonovich stochastic differential equations as well as the role played by the volume. (10.1051/cocv/2017037)
  DOI : 10.1051/cocv/2017037
• Random planar maps and growth-fragmentations
  
  • Bertoin Jean
  • Curien Nicolas
  • Kortchemski Igor
  The Annals of Probability, Institute of Mathematical Statistics, 2018, 46 (1), pp.207-260. (10.1214/17-AOP1183)
  DOI : 10.1214/17-AOP1183
• Fluctuations and Temperature Effects in Bose-Einstein Condensation
  
  • de Bouard Anne
  • Debussche Arnaud
  • Fukuizumi Reika
  • Poncet Romain
  ESAIM: Proceedings and Surveys, EDP Sciences, 2018, 61, pp.55-67. The modeling of cold atoms systems has known an increasing interest in the theoretical physics community, after the first experimental realizations of Bose Einstein condensates, some twenty years ago. We here review some analytical and numerical results concerning the influence of fluctuations , either arising from fluctuations of the confining parameters, or due to temperature effects, in the models describing the dynamics of such condensates. (10.1051/proc/201861055)
  DOI : 10.1051/proc/201861055
• Hilbert and Thompson geometries isometric to infinite-dimensional Banach spaces
  
  • Walsh Cormac
  Annales de l'Institut Fourier, Association des Annales de l'Institut Fourier, 2018, 68 (5), pp.1831-1877. We study the horofunction boundaries of Hilbert and Thompson geome-tries, and of Banach spaces, in arbitrary dimension. By comparing the boundaries of these spaces, we show that the only Hilbert and Thompson geometries that are isometric to Banach spaces are the ones defined on the cone of positive continuous functions on a compact space.
• Dynamic programming approach to principal-agent problems
  
  • Cvitanić Jakša
  • Possamaï Dylan
  • Touzi Nizar
  Finance and Stochastics, Springer Verlag (Germany), 2018, 22, pp.1-37. We consider a general formulation of the Principal-Agent problem with a lump-sum payment on a finite horizon, providing a systematic method for solving such problems. Our approach is the following: we first find the contract that is optimal among those for which the agent's value process allows a dynamic programming representation, for which the agent's optimal effort is straightforward to find. We then show that the optimization over the restricted family of contracts represents no loss of generality. As a consequence, we have reduced this non-zero sum stochastic differential game to a stochastic control problem which may be addressed by the standard tools of control theory. Our proofs rely on the backward stochastic differential equations approach to non-Markovian stochastic control, and more specifically, on the recent extensions to the second order case. (10.1007/s00780-017-0344-4)
  DOI : 10.1007/s00780-017-0344-4
• Modal basis approaches in shape and topology optimization of frequency response problems
  
  • Allaire Grégoire
  • Michailidis Georgios
  International Journal for Numerical Methods in Engineering, Wiley, 2018, 113 (8), pp.1258-1299. The optimal design of mechanical structures subject to periodic excitations within a large frequency interval is quite challenging. In order to avoid bad performances for non-discretized frequencies, it is necessary to finely discretize the frequency interval, leading to a very large number of state equations. Then, if a standard adjoint-based approach is used for optimization, the computational cost (both in terms of CPU and memory storage) may be prohibitive for large problems, especially in three space dimensions. The goal of the present work is to introduce two new non-adjoint approaches for dealing with frequency response problems in shape and topology optimization. In both cases, we rely on a classical modal basis approach to compute the states, solutions of the direct problems. In the first method, we do not use any adjoint but rather directly compute the shape derivatives of the eigenmodes in the modal basis. In the second method, we compute the adjoints of the standard approach by using again the modal basis. The numerical cost of these two new strategies are much smaller than the usual ones if the number of modes in the modal basis is much smaller than the number of discretized excitation frequencies. We present numerical examples for the minimization of the dynamic compliance in two and three space dimensions. (10.1002/nme.5504)
  DOI : 10.1002/nme.5504
• Relaxation Limit and Initial-Layers for a Class of Hyperbolic-Parabolic Systems
  
  • Giovangigli Vincent
  • Yang Zai-Bao
  • Yong Wen-An
  SIAM Journal on Mathematical Analysis, Society for Industrial and Applied Mathematics, 2018. We consider a class of hyperbolic-parabolic systems with small diffusion terms and stiff sources. Existence of solutions to the Cauchy problem with ill prepared initial data is established by using composite expansions including initial-layer correctors and a convergence-stability lemma. New multitime expansions are introduced and lead to second-order error estimates between the composite expansions and the solution. Reduced equilibrium systems of second-order accuracy are also investigated as well as initial-layers of Chapman-Enskog expansions. (10.1137/18M1170091)
  DOI : 10.1137/18M1170091
• Solving generic nonarchimedean semidefinite programs using stochastic game algorithms
  
  • Allamigeon Xavier
  • Gaubert Stephane
  • Skomra Mateusz
  Journal of Symbolic Computation, Elsevier, 2018, 85, pp.25-54. A general issue in computational optimization is to develop combinatorial algorithms for semidefinite programming. We address this issue when the base field is nonarchimedean. We provide a solution for a class of semidefinite feasibility problems given by generic matrices. Our approach is based on tropical geometry. It relies on tropical spectrahedra, which are defined as the images by the valuation of nonarchimedean spectrahedra. We establish a correspondence between generic tropical spectrahedra and zero-sum stochastic games with perfect information. The latter have been well studied in algorithmic game theory. This allows us to solve nonarchimedean semidefinite feasibility problems using algorithms for stochastic games. These algorithms are of a combinatorial nature and work for large instances. (10.1016/j.jsc.2017.07.002)
  DOI : 10.1016/j.jsc.2017.07.002
• Variational methods for tomographic reconstruction with few views
  
  • Bergounioux Maïtine
  • Abraham Isabelle
  • Abraham Romain
  • Carlier Guillaume
  • Le Pennec Erwan
  • Trélat Emmanuel
  Milan Journal of Mathematics, Springer Verlag, 2018, 86 (2), pp.157--200. We deal with a severe ill posed problem, namely the reconstruction process of an image during tomography acquisition with (very) few views. We present different methods that we have been investigated during the past decade. They are based on variational analysis. This is a survey paper and we refer to the quoted papers for more details. Mathematics Subject Classification (2010). 49K40, 45Q05,65M32.