Abstract In this paper, using Sinc-Galerkin method and TSVD regularization, an approximation of the quasi-solution to an inverse source problem is obtained. To do so, the solution of direct problem is obtained by the Sinc-Galerkin method, and this solution is applied in a least squares cost functional. Then, to obtain an approximation of the quasi-solution, we minimize the cost functional by TSVD regularization. Error analysis and convergence of the proposed method are investigated. In addition, at the end, four numerical examples are given in details to show the efficiency and accuracy of the proposed method.
{"title":"Identification of an Inverse Source Problem in a Fractional Partial Differential Equation Based on Sinc-Galerkin Method and TSVD Regularization","authors":"A. Safaie, A. H. Salehi Shayegan, M. Shahriari","doi":"10.1515/cmam-2022-0178","DOIUrl":"https://doi.org/10.1515/cmam-2022-0178","url":null,"abstract":"Abstract In this paper, using Sinc-Galerkin method and TSVD regularization, an approximation of the quasi-solution to an inverse source problem is obtained. To do so, the solution of direct problem is obtained by the Sinc-Galerkin method, and this solution is applied in a least squares cost functional. Then, to obtain an approximation of the quasi-solution, we minimize the cost functional by TSVD regularization. Error analysis and convergence of the proposed method are investigated. In addition, at the end, four numerical examples are given in details to show the efficiency and accuracy of the proposed method.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.3,"publicationDate":"2023-06-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41993212","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract The Miura ori is a very classical origami pattern used in numerous applications in engineering. A study of the shapes that surfaces using this pattern can assume is still lacking. A constrained nonlinear partial differential equation (PDE) that models the possible shapes that a periodic Miura tessellation can take in the homogenization limit has been established recently and solved only in specific cases. In this paper, the existence and uniqueness of a solution to the unconstrained PDE is proved for general Dirichlet boundary conditions. Then an H 2 H^{2} -conforming discretization is introduced to approximate the solution of the PDE coupled to a Newton method to solve the associated discrete problem. A convergence proof for the method is given as well as a convergence rate. Finally, numerical experiments show the robustness of the method and that nontrivial shapes can be achieved using periodic Miura tessellations.
{"title":"𝐻2-Conformal Approximation of Miura Surfaces","authors":"F. Marazzato","doi":"10.1515/cmam-2022-0259","DOIUrl":"https://doi.org/10.1515/cmam-2022-0259","url":null,"abstract":"Abstract The Miura ori is a very classical origami pattern used in numerous applications in engineering. A study of the shapes that surfaces using this pattern can assume is still lacking. A constrained nonlinear partial differential equation (PDE) that models the possible shapes that a periodic Miura tessellation can take in the homogenization limit has been established recently and solved only in specific cases. In this paper, the existence and uniqueness of a solution to the unconstrained PDE is proved for general Dirichlet boundary conditions. Then an H 2 H^{2} -conforming discretization is introduced to approximate the solution of the PDE coupled to a Newton method to solve the associated discrete problem. A convergence proof for the method is given as well as a convergence rate. Finally, numerical experiments show the robustness of the method and that nontrivial shapes can be achieved using periodic Miura tessellations.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.3,"publicationDate":"2023-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45319008","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Julian Roth, Jan Philipp Thiele, Uwe Köcher, Thomas Wick
Abstract In this work, the dual-weighted residual method is applied to a space-time formulation of nonstationary Stokes and Navier–Stokes flow. Tensor-product space-time finite elements are being used to discretize the variational formulation with discontinuous Galerkin finite elements in time and inf-sup stable Taylor–Hood finite element pairs in space. To estimate the error in a quantity of interest and drive adaptive refinement in time and space, we demonstrate how the dual-weighted residual method for incompressible flow can be extended to a partition-of-unity based error localization. We substantiate our methodology on 2D benchmark problems from computational fluid mechanics.
{"title":"Tensor-Product Space-Time Goal-Oriented Error Control and Adaptivity With Partition-of-Unity Dual-Weighted Residuals for Nonstationary Flow Problems","authors":"Julian Roth, Jan Philipp Thiele, Uwe Köcher, Thomas Wick","doi":"10.1515/cmam-2022-0200","DOIUrl":"https://doi.org/10.1515/cmam-2022-0200","url":null,"abstract":"Abstract In this work, the dual-weighted residual method is applied to a space-time formulation of nonstationary Stokes and Navier–Stokes flow. Tensor-product space-time finite elements are being used to discretize the variational formulation with discontinuous Galerkin finite elements in time and inf-sup stable Taylor–Hood finite element pairs in space. To estimate the error in a quantity of interest and drive adaptive refinement in time and space, we demonstrate how the dual-weighted residual method for incompressible flow can be extended to a partition-of-unity based error localization. We substantiate our methodology on 2D benchmark problems from computational fluid mechanics.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"12 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135394952","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract In the current work, we build a difference analog of the Caputo fractional derivative with generalized memory kernel (𝜇L2-1𝜎 formula). The fundamental features of this difference operator are studied, and on its ground, some difference schemes generating approximations of the second order in time for the generalized time-fractional diffusion equation with variable coefficients are worked out. We have proved stability and convergence of the given schemes in the grid L 2 L_{2} -norm with the rate equal to the order of the approximation error. The achieved results are supported by the numerical computations performed for some test problems.
{"title":"A Second-Order Difference Scheme for Generalized Time-Fractional Diffusion Equation with Smooth Solutions","authors":"A. Khibiev, A. Alikhanov, Chengming Huang","doi":"10.1515/cmam-2022-0089","DOIUrl":"https://doi.org/10.1515/cmam-2022-0089","url":null,"abstract":"Abstract In the current work, we build a difference analog of the Caputo fractional derivative with generalized memory kernel (𝜇L2-1𝜎 formula). The fundamental features of this difference operator are studied, and on its ground, some difference schemes generating approximations of the second order in time for the generalized time-fractional diffusion equation with variable coefficients are worked out. We have proved stability and convergence of the given schemes in the grid L 2 L_{2} -norm with the rate equal to the order of the approximation error. The achieved results are supported by the numerical computations performed for some test problems.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.3,"publicationDate":"2023-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43628392","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract We study variants of the mixed finite element method (mixed FEM) and the first-order system least-squares finite element (FOSLS) for the Poisson problem where we replace the load by a suitable regularization which permits to use H − 1 H^{-1} loads. We prove that any bounded H − 1 H^{-1} projector onto piecewise constants can be used to define the regularization and yields quasi-optimality of the lowest-order mixed FEM resp. FOSLS in weaker norms. Examples for the construction of such projectors are given. One is based on the adjoint of a weighted Clément quasi-interpolator. We prove that this Clément operator has second-order approximation properties. For the modified mixed method, we show optimal convergence rates of a postprocessed solution under minimal regularity assumptions—a result not valid for the lowest-order mixed FEM without regularization. Numerical examples conclude this work.
{"title":"On a Mixed FEM and a FOSLS with 𝐻−1 Loads","authors":"T. Führer","doi":"10.1515/cmam-2022-0215","DOIUrl":"https://doi.org/10.1515/cmam-2022-0215","url":null,"abstract":"Abstract We study variants of the mixed finite element method (mixed FEM) and the first-order system least-squares finite element (FOSLS) for the Poisson problem where we replace the load by a suitable regularization which permits to use H − 1 H^{-1} loads. We prove that any bounded H − 1 H^{-1} projector onto piecewise constants can be used to define the regularization and yields quasi-optimality of the lowest-order mixed FEM resp. FOSLS in weaker norms. Examples for the construction of such projectors are given. One is based on the adjoint of a weighted Clément quasi-interpolator. We prove that this Clément operator has second-order approximation properties. For the modified mixed method, we show optimal convergence rates of a postprocessed solution under minimal regularity assumptions—a result not valid for the lowest-order mixed FEM without regularization. Numerical examples conclude this work.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"0 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2023-05-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"66871017","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract We study frequency domain acoustic scattering at a bounded, penetrable, and inhomogeneous obstacle Ω − ⊂ R d Omega^{-}subsetmathbb{R}^{d} , d = 2 , 3 d=2,3 . By defining constant reference coefficients, a representation formula for the pressure field is derived. It contains a volume integral operator, related to the one in the Lippmann–Schwinger equation. Besides, it features integral operators defined on ∂ Ω − partialOmega^{-} and closely related to boundary integral equations of single-trace formulations (STF) for transmission problems with piecewise constant coefficients. We show well-posedness of the continuous variational formulation and asymptotic convergence of Galerkin discretizations. Numerical experiments in 2D validate our expected convergence rates.
{"title":"Volume Integral Equations and Single-Trace Formulations for Acoustic Wave Scattering in an Inhomogeneous Medium","authors":"Ignacio Labarca, R. Hiptmair","doi":"10.1515/cmam-2022-0119","DOIUrl":"https://doi.org/10.1515/cmam-2022-0119","url":null,"abstract":"Abstract We study frequency domain acoustic scattering at a bounded, penetrable, and inhomogeneous obstacle Ω − ⊂ R d Omega^{-}subsetmathbb{R}^{d} , d = 2 , 3 d=2,3 . By defining constant reference coefficients, a representation formula for the pressure field is derived. It contains a volume integral operator, related to the one in the Lippmann–Schwinger equation. Besides, it features integral operators defined on ∂ Ω − partialOmega^{-} and closely related to boundary integral equations of single-trace formulations (STF) for transmission problems with piecewise constant coefficients. We show well-posedness of the continuous variational formulation and asymptotic convergence of Galerkin discretizations. Numerical experiments in 2D validate our expected convergence rates.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.3,"publicationDate":"2023-05-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46463057","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract In this work, a cost-efficient space-time adaptive algorithm based on the Dual Weighted Residual (DWR) method is developed and studied for a coupled model problem of flow and convection-dominated transport. Key ingredients are a multirate approach adapted to varying dynamics in time of the subproblems, weighted and non-weighted error indicators for the transport and flow problem, respectively, and the concept of space-time slabs based on tensor product spaces for the data structure. In numerical examples, the performance of the underlying algorithm is studied for benchmark problems and applications of practical interest. Moreover, the interaction of stabilization and goal-oriented adaptivity is investigated for strongly convection-dominated transport.
{"title":"A Cost-Efficient Space-Time Adaptive Algorithm for Coupled Flow and Transport","authors":"Marius Paul Bruchhäuser, Markus Bause","doi":"10.1515/cmam-2022-0245","DOIUrl":"https://doi.org/10.1515/cmam-2022-0245","url":null,"abstract":"Abstract In this work, a cost-efficient space-time adaptive algorithm based on the Dual Weighted Residual (DWR) method is developed and studied for a coupled model problem of flow and convection-dominated transport. Key ingredients are a multirate approach adapted to varying dynamics in time of the subproblems, weighted and non-weighted error indicators for the transport and flow problem, respectively, and the concept of space-time slabs based on tensor product spaces for the data structure. In numerical examples, the performance of the underlying algorithm is studied for benchmark problems and applications of practical interest. Moreover, the interaction of stabilization and goal-oriented adaptivity is investigated for strongly convection-dominated transport.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"235 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-05-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"134922897","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract In this paper, we develop a new immersed finite element method (IFEM) for two-phase incompressible Stokes flows. We allow the interface to cut the finite elements. On the noninterface element, the standard Crouzeix–Raviart element and the P 0 {P_{0}} element pair is used. On the interface element, the basis functions developed for scalar interface problems (Kwak et al., An analysis of a broken P 1 {P_{1}} -nonconforming finite element method for interface problems, SIAM J. Numer. Anal. (2010)) are modified in such a way that the coupling between the velocity and pressure variable is different. There are two kinds of basis functions. The first kind of basis satisfies the Laplace–Young condition under the assumption of the continuity of the pressure variable. In the second kind, the velocity is of bubble type and is coupled with the discontinuous pressure, still satisfying the Laplace–Young condition. We remark that in the second kind the pressure variable has two degrees of freedom on each interface element. Therefore, our methods can handle the discontinuous pressure case. Numerical results including the case of the discontinuous pressure variable are provided. We see optimal convergence orders for all examples.
摘要本文提出了一种新的求解两相不可压缩Stokes流的浸入式有限元方法。我们允许界面切割有限元。在非接口元素上,使用标准的Crouzeix–Raviart元素和P0{P_{0}}元素对。在界面单元上,对为标量界面问题开发的基函数(Kwak et al.,An analysis of a breaked P1{P_{1}}-conformant finite element method for interface problems,SIAM J.Numer.Anal.(2010))进行了修改,使得速度和压力变量之间的耦合不同。基函数有两种。在压力变量连续性的假设下,第一类基满足拉普拉斯-杨条件。在第二种情况下,速度是气泡型的,并与不连续压力耦合,仍然满足拉普拉斯-杨条件。我们注意到,在第二类中,压力变量在每个界面元件上有两个自由度。因此,我们的方法可以处理不连续压力的情况。给出了包括不连续压力变量情况下的数值结果。我们看到所有例子的最优收敛阶。
{"title":"A New Immersed Finite Element Method for Two-Phase Stokes Problems Having Discontinuous Pressure","authors":"Gwanghyun Jo, D. Kwak","doi":"10.1515/cmam-2022-0122","DOIUrl":"https://doi.org/10.1515/cmam-2022-0122","url":null,"abstract":"Abstract In this paper, we develop a new immersed finite element method (IFEM) for two-phase incompressible Stokes flows. We allow the interface to cut the finite elements. On the noninterface element, the standard Crouzeix–Raviart element and the P 0 {P_{0}} element pair is used. On the interface element, the basis functions developed for scalar interface problems (Kwak et al., An analysis of a broken P 1 {P_{1}} -nonconforming finite element method for interface problems, SIAM J. Numer. Anal. (2010)) are modified in such a way that the coupling between the velocity and pressure variable is different. There are two kinds of basis functions. The first kind of basis satisfies the Laplace–Young condition under the assumption of the continuity of the pressure variable. In the second kind, the velocity is of bubble type and is coupled with the discontinuous pressure, still satisfying the Laplace–Young condition. We remark that in the second kind the pressure variable has two degrees of freedom on each interface element. Therefore, our methods can handle the discontinuous pressure case. Numerical results including the case of the discontinuous pressure variable are provided. We see optimal convergence orders for all examples.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.3,"publicationDate":"2023-04-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"47608032","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract We present the self-consistent Pauli equation, a semi-relativistic model for charged spin- 1 / 2 1/2 particles with self-interaction with the electromagnetic field. The Pauli equation arises as the O ( 1 / c ) O(1/c) approximation of the relativistic Dirac equation. The fully relativistic self-consistent model is the Dirac–Maxwell equation where the description of spin and the magnetic field arises naturally. In the non-relativistic setting, the correct self-consistent equation is the Schrödinger–Poisson equation which does not describe spin and the magnetic field and where the self-interaction is with the electric field only. The Schrödinger–Poisson equation also arises as the mean field limit of the 𝑁-body Schrödinger equation with Coulomb interaction. We propose that the Pauli–Poisson equation arises as the mean field limit N → ∞ Ntoinfty of the linear 𝑁-body Pauli equation with Coulomb interaction where one has to pay extra attention to the fermionic nature of the Pauli equation. We present the semiclassical limit of the Pauli–Poisson equation by the Wigner method to the Vlasov equation with Lorentz force coupled to the Poisson equation which is also consistent with the hierarchy in 1 / c 1/c of the self-consistent Vlasov equation. This is a non-trivial extension of the groundbreaking works by Lions & Paul and Markowich & Mauser, where we need methods like magnetic Lieb–Thirring estimates.
{"title":"Nonlinear PDE Models in Semi-relativistic Quantum Physics","authors":"Jakob Möller, N. Mauser","doi":"10.1515/cmam-2023-0101","DOIUrl":"https://doi.org/10.1515/cmam-2023-0101","url":null,"abstract":"Abstract We present the self-consistent Pauli equation, a semi-relativistic model for charged spin- 1 / 2 1/2 particles with self-interaction with the electromagnetic field. The Pauli equation arises as the O ( 1 / c ) O(1/c) approximation of the relativistic Dirac equation. The fully relativistic self-consistent model is the Dirac–Maxwell equation where the description of spin and the magnetic field arises naturally. In the non-relativistic setting, the correct self-consistent equation is the Schrödinger–Poisson equation which does not describe spin and the magnetic field and where the self-interaction is with the electric field only. The Schrödinger–Poisson equation also arises as the mean field limit of the 𝑁-body Schrödinger equation with Coulomb interaction. We propose that the Pauli–Poisson equation arises as the mean field limit N → ∞ Ntoinfty of the linear 𝑁-body Pauli equation with Coulomb interaction where one has to pay extra attention to the fermionic nature of the Pauli equation. We present the semiclassical limit of the Pauli–Poisson equation by the Wigner method to the Vlasov equation with Lorentz force coupled to the Poisson equation which is also consistent with the hierarchy in 1 / c 1/c of the self-consistent Vlasov equation. This is a non-trivial extension of the groundbreaking works by Lions & Paul and Markowich & Mauser, where we need methods like magnetic Lieb–Thirring estimates.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.3,"publicationDate":"2023-04-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"47920744","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}