We propose and analyze a new mixed finite element method for the nonlinear problem given by the coupling of the steady Brinkman--Forchheimer and double-diffusion equations. Besides the velocity, temperature, and concentration, our approach introduces the velocity gradient, the pseudostress tensor, and a pair of vectors involving the temperature/concentration, its gradient and the velocity, as further unknowns. As a consequence, we obtain a fully mixed variational formulation presenting a Banach spaces framework in each set of equations. In this way, and differently from the techniques previously developed for this and related coupled problems, no augmentation procedure needs to be incorporated now into the formulation nor into the solvability analysis. The resulting non-augmented scheme is then written equivalently as a fixed-point equation, so that the well-known Banach theorem, combined with classical results on nonlinear monotone operators and the Babusska-Brezzi theory in Banach spaces, are applied to prove the unique solvability of the continuous and discrete systems. Appropriate finite element subspaces satisfying the required discrete inf-sup conditions are specified, and optimal a priori error estimates are derived. Several numerical examples confirm the theoretical rates of convergence and illustrate the performance and flexibility of the method.
{"title":"A fully-mixed formulation in Banach spaces for the coupling of the steady Brinkman-Forchheimer and double-diffusion equations","authors":"Sergio Caucao, G. Gatica, J. Ortega","doi":"10.1051/m2an/2021072","DOIUrl":"https://doi.org/10.1051/m2an/2021072","url":null,"abstract":"We propose and analyze a new mixed finite element method for the nonlinear problem given by the coupling of the steady Brinkman--Forchheimer and double-diffusion equations. Besides the velocity, temperature, and concentration, our approach introduces the velocity gradient, the pseudostress tensor, and a pair of vectors involving the temperature/concentration, its gradient and the velocity, as further unknowns. As a consequence, we obtain a fully mixed variational formulation presenting a Banach spaces framework in each set of equations. In this way, and differently from the techniques previously developed for this and related coupled problems, no augmentation procedure needs to be incorporated now into the formulation nor into the solvability analysis. The resulting non-augmented scheme is then written equivalently as a fixed-point equation, so that the well-known Banach theorem, combined with classical results on nonlinear monotone operators and the Babusska-Brezzi theory in Banach spaces, are applied to prove the unique solvability of the continuous and discrete systems. Appropriate finite element subspaces satisfying the required discrete inf-sup conditions are specified, and optimal a priori error estimates are derived. Several numerical examples confirm the theoretical rates of convergence and illustrate the performance and flexibility of the method.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2021-10-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"88715966","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}
R. Lipton, Robert Viator, Silvia Jiménez Bolanos, A. Adili
Analytic representation formulas and power series are developed describing the band structure inside non-magnetic periodic photonic three-dimensional crystals made from high dielectric contrast inclusions. Central to this approach is the identifcation and utilization of a resonance spectrum for quasiperiodic source-free modes. These modes are used to represent solution operators associated with electromagnetic and acoustic waves inside periodic high contrast media. A convergent power series for the Bloch wave spectrum is recovered from the representation formulas. Explicit conditions on the contrast are found that provide lower bounds on the convergence radius. These conditions are sufficient for the separation of spectral branches of the dispersion relation for any fixed quasi-momentum.
{"title":"Bloch waves in high contrast electromagnetic crystals","authors":"R. Lipton, Robert Viator, Silvia Jiménez Bolanos, A. Adili","doi":"10.1051/m2an/2022045","DOIUrl":"https://doi.org/10.1051/m2an/2022045","url":null,"abstract":"Analytic representation formulas and power series are developed describing the band structure inside non-magnetic periodic photonic three-dimensional crystals made from high dielectric contrast inclusions. Central to this approach is the identifcation and utilization of a resonance spectrum for quasiperiodic source-free modes. These modes are used to represent solution operators associated with electromagnetic and acoustic waves inside periodic high contrast media. A convergent power series for the Bloch wave spectrum is recovered from the representation formulas. Explicit conditions on the contrast are found that provide lower bounds on the convergence radius. These conditions are sufficient for the separation of spectral branches of the dispersion relation for any fixed quasi-momentum.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2021-10-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"83711153","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}
M. Badsi, B. Després, M. Campos-Pinto, Ludovic Godard-Cadillac
We construct a stationary gyrokinetic variational model for sheaths close to the metallic wall of a magnetized plasma, following a physical extremalization principle for the natural energy. By considering a reduced set of parameters we show that our model has a unique minimal solution, and that the resulting electric potential has an infinite number of oscillations as it propagates towards the core of the plasma. We prove this result for the non linear problem and also provide a simpler analysis for a linearized problem, based on the construction of exact solutions. Some numerical illustrations show the well-posedness of the model after numerical discretization. They also exhibit the oscillating behavior.
{"title":"A variational sheath model for gyrokinetic Vlasov-Poisson equations","authors":"M. Badsi, B. Després, M. Campos-Pinto, Ludovic Godard-Cadillac","doi":"10.1051/m2an/2021067","DOIUrl":"https://doi.org/10.1051/m2an/2021067","url":null,"abstract":"We construct a stationary gyrokinetic variational model for sheaths close to the metallic wall of a magnetized plasma, following a physical extremalization principle for the natural energy. By considering a reduced set of parameters we show that our model has a unique minimal solution, and that the resulting electric potential has an infinite number of oscillations as it propagates towards the core of the plasma. We prove this result for the non linear problem and also provide a simpler analysis for a linearized problem, based on the construction of exact solutions. Some numerical illustrations show the well-posedness of the model after numerical discretization. They also exhibit the oscillating behavior.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2021-10-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"79554765","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}
We study unsaturated poroelasticity, i.e., coupled hydro-mechanical processes in variably saturated porous media, here modeled by a non-linear extension of Biot's well-known quasi-static consolidation model. The coupled elliptic-parabolic system of partial differential equations is a simplified version of the general model for multi-phase flow in deformable porous media, obtained under similar assumptions as usually considered for Richards' equation. In this work, existence of weak solutions is established in several steps involving a numerical approximation of the problem using a physically-motivated regularization and a finite element/finite volume discretization. Eventually, solvability of the original problem is proved by a combination of the Rothe and Galerkin methods, and further compactness arguments. This approach in particular provides the convergence of the numerical discretization to a regularized model for unsaturated poroelasticity. The final existence result holds under non-degeneracy conditions and natural continuity properties for the constitutive relations. The assumptions are demonstrated to be reasonable in view of geotechnical applications.
{"title":"Global existence of weak solutions to unsaturated poroelasticity","authors":"J. Both, Iuliu Sorin Pop, I. Yotov","doi":"10.1051/m2an/2021063","DOIUrl":"https://doi.org/10.1051/m2an/2021063","url":null,"abstract":"We study unsaturated poroelasticity, i.e., coupled hydro-mechanical processes in variably saturated porous media, here modeled by a non-linear extension of Biot's well-known quasi-static consolidation model. The coupled elliptic-parabolic system of partial differential equations is a simplified version of the general model for multi-phase flow in deformable porous media, obtained under similar assumptions as usually considered for Richards' equation. In this work, existence of weak solutions is established in several steps involving a numerical approximation of the problem using a physically-motivated regularization and a finite element/finite volume discretization. Eventually, solvability of the original problem is proved by a combination of the Rothe and Galerkin methods, and further compactness arguments. This approach in particular provides the convergence of the numerical discretization to a regularized model for unsaturated poroelasticity. The final existence result holds under non-degeneracy conditions and natural continuity properties for the constitutive relations. The assumptions are demonstrated to be reasonable in view of geotechnical applications.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2021-10-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"80711520","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}
In this article, we consider the convection-diffusion-reaction problem coupled the Darcy-Forchheimer problem by a non-linear external force depending on the concentration. We establish existence of a solution by using a Galerkin method and we prove uniqueness. We introduce and analyse a numerical scheme based on the finite element method. An optimal a priori error estimate is then derived for each numerical scheme. Numerical investigation are performed to confirm the theoretical accuracy of the discretization.
{"title":"Finite element methods for the Darcy-Forchheimer problem coupled with the convection-diffusion-reaction problem","authors":"Toni Sayah, G. Semaan, Faouzi Triki","doi":"10.1051/m2an/2021066","DOIUrl":"https://doi.org/10.1051/m2an/2021066","url":null,"abstract":"In this article, we consider the convection-diffusion-reaction problem coupled the Darcy-Forchheimer problem by a non-linear external force depending on the concentration. We establish existence of a solution by using a Galerkin method and we prove uniqueness. We introduce and analyse a numerical scheme based on the finite element method. An optimal a priori error estimate is then derived for each numerical scheme. Numerical investigation are performed to confirm the theoretical accuracy of the discretization.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2021-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90833245","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}
In this work, a new Virtual Element Method (VEM) of arbitrary order $k geq 2$ for the time dependent Navier-Stokes equations in stream-function form is proposed and analyzed. Using suitable projection operators, the bilinear and trilinear terms are discretized by only using the proposed degrees of freedom associated with the virtual space. Under certain assumptions on the computational domain, error estimations are derived and shown that the method is optimally convergent in both space and time variables. Finally, to justify the theoretical analysis, four benchmark examples are examined numerically.
{"title":"A virtual element discretization for the time dependent Navier-Stokes equations in stream-function formulation","authors":"D. Adak, D. Mora, S. Natarajan, A. Silgado","doi":"10.1051/m2an/2021058","DOIUrl":"https://doi.org/10.1051/m2an/2021058","url":null,"abstract":"In this work, a new Virtual Element Method (VEM) of arbitrary order $k geq 2$ for the time dependent Navier-Stokes equations in stream-function form is proposed and analyzed. Using suitable projection operators, the bilinear and trilinear terms are discretized by only using the proposed degrees of freedom associated with the virtual space. Under certain assumptions on the computational domain, error estimations are derived and shown that the method is optimally convergent in both space and time variables. Finally, to justify the theoretical analysis, four benchmark examples are examined numerically.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2021-10-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78844140","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}
A recent paper introduced the network element method (NEM) where the usual mesh was replaced by a discretization network. Using the associated network geometric coefficients and following the virtual element framework, a consistent and stable numerical scheme was proposed. The aim of the present paper is to derive a convergence theory for the NEM under mild assumptions on the exact problem. We also derive basic error estimates, which are sub-optimal in the sense that we have to assume more regularity than usual.
{"title":"Basic convergence theory for the network element method","authors":"J. Coatléven","doi":"10.1051/m2an/2021062","DOIUrl":"https://doi.org/10.1051/m2an/2021062","url":null,"abstract":"A recent paper introduced the network element method (NEM) where the usual mesh was replaced by a discretization network. Using the associated network geometric coefficients and following the virtual element framework, a consistent and stable numerical scheme was proposed. The aim of the present paper is to derive a convergence theory for the NEM under mild assumptions on the exact problem. We also derive basic error estimates, which are sub-optimal in the sense that we have to assume more regularity than usual.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2021-09-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90678348","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}
We present a structured coagulation-fragmentation model which describes the population dynamics of oceanic phytoplankton. This model is formulated on the space of Radon measures equipped with the bounded Lipschitz norm and unifies the study of the discrete and continuous coagulation-fragmentation models. We prove that the model is well-posed and show it can reduce down to the classic discrete and continuous coagulation-fragmentation models. To understand the interplay between the physical processes of coagulation and fragmentation and the biological processes of growth, reproduction, and death, we establish a regularity result for the solutions and use it to show that stationary solutions are absolutely continuous under some conditions on model parameters. We develop a semi-discrete approximation scheme which conserves mass and prove its convergence to the unique weak solution. We then use the scheme to perform numerical simulations for the model.
{"title":"Structured coagulation-fragmentation equation in the space of radon measures: Unifying discrete and continuous models","authors":"Azmy S. Ackleh, Rainey Lyons, N. Saintier","doi":"10.1051/m2an/2021061","DOIUrl":"https://doi.org/10.1051/m2an/2021061","url":null,"abstract":"We present a structured coagulation-fragmentation model which describes the population dynamics of oceanic phytoplankton. This model is formulated on the space of Radon measures equipped with the bounded Lipschitz norm and unifies the study of the discrete and continuous coagulation-fragmentation models. We prove that the model is well-posed and show it can reduce down to the classic discrete and continuous coagulation-fragmentation models. To understand the interplay between the physical processes of coagulation and fragmentation and the biological processes of growth, reproduction, and death, we establish a regularity result for the solutions and use it to show that stationary solutions are absolutely continuous under some conditions on model parameters. We develop a semi-discrete approximation scheme which conserves mass and prove its convergence to the unique weak solution. We then use the scheme to perform numerical simulations for the model.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2021-09-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78745240","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}
Electrical Impedance Tomography (EIT) is a well-known imaging technique for detecting the electrical properties of an object in order to detect anomalies, such as conductive or resistive targets. More specifically, EIT has many applications in medical imaging for the detection and location of bodily tumors since it is an affordable and non-invasive method, which aims to recover the internal conductivity of a body using voltage measurements resulting from applying low frequency current at electrodes placed at its surface. Mathematically, the reconstruction of the internal conductivity is a severely ill-posed inverse problem and yields a poor quality image reconstruction. To remedy this difficulty, at least in part, we regularize and solve the nonlinear minimization problem by the aid of a Krylov subspace-type method for the linear sub problem during each iteration. In EIT, a tumor or general anomaly can be modeled as a piecewise constant perturbation of a smooth background, hence, we solve the regularized problem on a subspace of relatively small dimension by the Flexible Golub-Kahan process that provides solutions that have sparse representation. For comparison, we use a well-known modified Gauss-Newton algorithm as a benchmark. Using simulations, we demonstrate the effectiveness of the proposed method. The obtained reconstructions indicate that the Krylov subspace method is better adapted to solve the ill-posed EIT problem and results in higher resolution images and faster convergence compared to reconstructions using the modified Gauss-Newton algorithm.
{"title":"A Krylov subspace type method for Electrical Impedance Tomography","authors":"M. Pasha, Shyla Kupis, S. Ahmad, T. Khan","doi":"10.1051/m2an/2021057","DOIUrl":"https://doi.org/10.1051/m2an/2021057","url":null,"abstract":"Electrical Impedance Tomography (EIT) is a well-known imaging technique for detecting the electrical properties of an object in order to detect anomalies, such as conductive or resistive targets. More specifically, EIT has many applications in medical imaging for the detection and location of bodily tumors since it is an affordable and non-invasive method, which aims to recover the internal conductivity of a body using voltage measurements resulting from applying low frequency current at electrodes placed at its surface.\u0000\u0000Mathematically, the reconstruction of the internal conductivity is a severely ill-posed inverse problem and yields a poor quality image reconstruction. To remedy this difficulty, at least in part, we regularize and solve the nonlinear minimization problem by the aid of a Krylov subspace-type method for the linear sub problem during each iteration. In EIT, a tumor or general anomaly can be modeled as a piecewise constant perturbation of a smooth background, hence, we solve the regularized problem on a subspace of relatively small dimension by the Flexible Golub-Kahan process that provides solutions that have sparse representation. For comparison, we use a well-known modified Gauss-Newton algorithm as a benchmark. Using simulations, we demonstrate the effectiveness of the proposed method. The obtained reconstructions indicate that the Krylov subspace method is better adapted to solve the ill-posed EIT problem and results in higher resolution images and faster convergence compared to reconstructions using the modified Gauss-Newton algorithm.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2021-09-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"83167036","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}
In this paper, the a prior error estimates of an embedded discontinuous Galerkin method for the Oseen equations are presented. It is proved that the velocity error in the L 2 (Ω) norm, has an optimal error bound with convergence order k + 1, where the constants are dependent on the Reynolds number (or ν − 1 ), in the diffusion-dominated regime, and in the convection-dominated regime, it has a Reynolds-robust error bound with quasi-optimal convergence order k +1 / 2. Here, k is the polynomial order of the velocity space. In addition, we also prove an optimal error estimate for the pressure. Finally, we carry out some numerical experiments to corroborate our analytical results.
{"title":"An embedded discontinuous Galerkin method for the Oseen equations","authors":"Yongbin Han, Yanren Hou","doi":"10.1051/m2an/2021059","DOIUrl":"https://doi.org/10.1051/m2an/2021059","url":null,"abstract":"In this paper, the a prior error estimates of an embedded discontinuous Galerkin method for the Oseen equations are presented. It is proved that the velocity error in the L 2 (Ω) norm, has an optimal error bound with convergence order k + 1, where the constants are dependent on the Reynolds number (or ν − 1 ), in the diffusion-dominated regime, and in the convection-dominated regime, it has a Reynolds-robust error bound with quasi-optimal convergence order k +1 / 2. Here, k is the polynomial order of the velocity space. In addition, we also prove an optimal error estimate for the pressure. Finally, we carry out some numerical experiments to corroborate our analytical results.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2021-09-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"85053316","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}
京公网安备 11010802042870号
京ICP备2023020795号-1
扫码关注我们
求助内容：
应助结果提醒方式：