{"title":"A LWR model with constraints at moving interfaces","authors":"A. Sylla","doi":"10.1051/m2an/2022030","DOIUrl":"https://doi.org/10.1051/m2an/2022030","url":null,"abstract":"We propose a mathematical framework to the study of scalar conservation laws with moving interfaces. This framework is developed on a LWR model with constraint on the flux along these moving interfaces. Existence is proved by means of a finite volume scheme. The originality lies in the local modification of the mesh and in the treatment of the crossing points of the trajectories.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"77855698","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Analysis of compressible bubbly flows. Part II: Derivation of a macroscopic model.","authors":"M. Hillairet, H. Mathis, N. Seguin","doi":"10.1051/m2an/2023046","DOIUrl":"https://doi.org/10.1051/m2an/2023046","url":null,"abstract":"This paper is the second of the series of two papers, which focuses on the derivation of an averaged 1D model for compressible bubbly flows. For this, we start from a microscopic description of the interactions between a large but finite number of small bubbles with a surrounding compressible fluid. This microscopic model has been derived and analysed in the first paper. In the present one, provided physical parameters scale according to the number of bubbles, we prove that solutions to the microscopic model exist on a timespan independent of the number of bubbles. Considering then that we have a large number of bubbles, we propose a construction of the macroscopic variables and derive the averaged system satisfied by these quantities. Our method is based on a compactness approach in a strong-solution setting. In the last section, we propose the derivation of the Williams-Boltzmann equation corresponding to our setting.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"74868093","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Monotone discretization of the Monge-Ampère equation of optimal transport","authors":"G. Bonnet, J. Mirebeau","doi":"10.1051/m2an/2022029","DOIUrl":"https://doi.org/10.1051/m2an/2022029","url":null,"abstract":"We design a monotone finite difference discretization of the second boundary value problem for the Monge-Ampère equation, whose main application is optimal transport. We prove the existence of solutions to a class of monotone numerical schemes for degenerate elliptic equations whose sets of solutions are stable by addition of a constant, and we show that the scheme that we introduce for the Monge-Ampère equation belongs to this class. We prove the convergence of this scheme, although only in the setting of quadratic optimal transport. The scheme is based on a reformulation of the Monge-Ampère operator as a maximum of semilinear operators. In dimension two, we recommend to use Selling's formula, a tool originating from low-dimensional lattice geometry, in order to choose the parameters of the discretization. We show that this approach yields a closed-form formula for the maximum that appears in the discretized operator, which allows the scheme to be solved particularly efficiently. We present some numerical results that we obtained by applying the scheme to quadratic optimal transport problems as well as to the far field refractor problem in nonimaging optics.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-03-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"79898689","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Convergence analysis of a fully discrete finite element method for thermally coupled incompressible MHD problems with temperature-dependent coefficients","authors":"Qianqian Ding, Xiaonian Long, S. Mao","doi":"10.1051/m2an/2022028","DOIUrl":"https://doi.org/10.1051/m2an/2022028","url":null,"abstract":"In this paper, we study a fully discrete finite element scheme of thermally coupled incompressible magnetohydrodynamic with temperature-dependent coefficients in Lipschitz domain. The variable coefficients in the MHD system and possible nonconvex domain may cause nonsmooth solutions. We propose a fully discrete Euler semi-implicit scheme with the magnetic equation approximated by Ne´de´lec edge elements to capture the physical solutions. The fully discrete scheme only needs to solve one linear system at each time step and is unconditionally stable. Utilizing the stability of the numerical scheme and the compactness method, the existence of weak solution to the thermally coupled MHD model in three dimensions is established. Furthermore, the uniqueness of weak solution and the convergence of the proposed numerical method are also rigorously derived. Under the hypothesis of a low regularity for the exact solution, we rigorously establish the error estimates for the velocity, temperature and magnetic induction unconditionally in the sense that the time step is independent of the spacial mesh size.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-03-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"88308365","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Formulation and Convergence Analysis of New Methods for Reduced Fragmentation Model: Illustrative Application to Depolymerization","authors":"Mehakpreet Singh, Gavin Walker, V. Ranade","doi":"10.1051/m2an/2022023","DOIUrl":"https://doi.org/10.1051/m2an/2022023","url":null,"abstract":"In this work, two discrete formulations based on the finite volume approach for a reduced fragmentation model are developed. The important features such as mass conservation and accurate prediction of the zeroth order moments are accomplished by the modification of the selection function. The new schemes can compute the second order moment, which plays a significant role in predicting the area of the particles in real life applications, with high accuracy without taking any specific measures. A thorough convergence analysis of both schemes including Lipschitz condition and consistency is presented and exhibit second order convergence. The accuracy and efficiency of both schemes is demonstrated by applying them to the depolymerization problem which commonly arises in polymer sciences and chemical engineering. It is demonstrated that the new schemes are easy to implement, computationally efficient and able to compute the numerical results with higher precision even on a coarser grid.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-03-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"89339398","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An algorithm for the grade-two rheological model","authors":"Sara N. Pollock, L. Ridgway Scott","doi":"10.1051/m2an/2022024","DOIUrl":"https://doi.org/10.1051/m2an/2022024","url":null,"abstract":"We develop an algorithm for solving the general grade-two model of non-Newtonian fluids which for the first time includes inflow boundary conditions. The algorithm also allows for both of the rheological parameters to be chosen independently. The proposed algorithm couples a Stokes equation for the velocity with a transport equation for an auxiliary vector-valued function. We prove that this model is well posed using the algorithm that we show converges geometrically in suitable Sobolev spaces for sufficiently small data. We demonstrate computationally that this algorithm can be successfully discretized and that it can converge to solutions for the model parameters of order one. We include in the appendix a description of appropriate boundary conditions for the auxiliary variable in standard geometries.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"79318506","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Linear programming fictitious play algorithm for mean field games with optimal stopping and absorption","authors":"Roxana Dumitrescu, Marcos Leutscher, P. Tankov","doi":"10.1051/m2an/2023019","DOIUrl":"https://doi.org/10.1051/m2an/2023019","url":null,"abstract":"We develop the fictitious play algorithm in the context of the linear programming approach for mean field games of optimal stopping and mean field games with regular control and absorption. This algorithm allows to approximate the mean field game population dynamics without computing the value function by solving linear programming problems associated with the distributions of the players still in the game and their stopping times/controls. We show the convergence of the algorithm using the topology of convergence in measure in the space of subprobability measures, which is needed to deal with the lack of continuity of the flows of measures. Numerical examples are provided to illustrate the convergence of the algorithm.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-02-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81268226","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Continuous Lambertian shape from shading: a primal-dual algorithm","authors":"Hamza Ennaji, N. Igbida, Van Thanh Nguyen","doi":"10.1051/m2an/2022014","DOIUrl":"https://doi.org/10.1051/m2an/2022014","url":null,"abstract":"The continuous Lambertian shape from shading is studied using a PDE approach\u0000\u0000in terms of Hamilton–Jacobi equations. The latter will then be characterized by a maximization\u0000\u0000problem. In this paper we show the convergence of discretization and propose to use the wellknown\u0000\u0000Chambolle–Pock primal-dual algorithm to solve numerically the shape from shading\u0000\u0000problem. The saddle-point structure of the problem makes the Chambolle–Pock algorithm\u0000\u0000suitable to approximate solutions of the discretized problems.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-02-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"86103206","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Model adaptation for non-linear elliptic equations in mixed form: existence of solutions and numerical strategies","authors":"A. Fumagalli, F. Patacchini","doi":"10.1051/m2an/2022016","DOIUrl":"https://doi.org/10.1051/m2an/2022016","url":null,"abstract":"Depending on the physical and geometrical properties of a given porous medium, fluid flow can behave differently, going from a slow Darcian regime to more complicated Brinkman or even Forchheimer regimes for high velocity. The main problem is to determine where in the medium one regime is more adequate than others. In order to determine the low-speed and high-speed regions, this work proposes an adaptive strategy which is based on selecting the appropriate constitutive law linking velocity and pressure according to a threshold criterion on the magnitude of the fluid velocity itself. Both theoretical and numerical aspects are considered and investigated, showing the potentiality of the proposed approach. From the analytical viewpoint, we show existence of weak solutions to such model under reasonable hypotheses on the constitutive laws. To this end, we use a variational approach identifying solutions with minimizers of an underlying energy functional. From the numerical viewpoint, we propose a one-dimensional algorithm which tracks the transition zone between the low- and high-speed regions. By running numerical experiments using this algorithm, we illustrate some interesting behaviors of our adaptive model on academic cases and on small networks of intersecting fractures.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"80681323","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Mathematical analysis of Goldstein's model for time harmonic acoustics in flow","authors":"A. Bensalah, P. Joly, J. Mercier","doi":"10.1051/m2an/2022007","DOIUrl":"https://doi.org/10.1051/m2an/2022007","url":null,"abstract":"Goldstein’s equations have been introduced in 1978 as an alternative model to linearized Euler equations to model acoustic waves in moving fluids. This new model is particularly attractive since it appears as a perturbation a simple scalar model: the potential model. In this work we propose a mathematical analysis of boundary value problems associated with Goldstein’s equations in the time-harmonic regime.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-01-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"73365632","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}