首页 > 最新文献

Computational Methods in Applied Mathematics最新文献

英文 中文
A Domain Decomposition Scheme for Couplings between Local and Nonlocal Equations 局部与非局部方程耦合的一种区域分解格式
IF 1.3 4区 数学 Q3 MATHEMATICS, APPLIED Pub Date : 2022-12-12 DOI: 10.48550/arXiv.2212.06093
Gabriel Acosta, Francisco M. Bersetche, J. Rossi
Abstract We study a natural alternating method of Schwarz type (domain decomposition) for a certain class of couplings between local and nonlocal operators. We show that our method fits into Lions’s framework and prove, as a consequence, convergence in both the continuous and the discrete settings.
摘要研究了一类局部算子与非局部算子耦合的Schwarz型自然交替方法(域分解)。我们证明了我们的方法符合Lions的框架,并因此证明了在连续和离散环境中的收敛性。
{"title":"A Domain Decomposition Scheme for Couplings between Local and Nonlocal Equations","authors":"Gabriel Acosta, Francisco M. Bersetche, J. Rossi","doi":"10.48550/arXiv.2212.06093","DOIUrl":"https://doi.org/10.48550/arXiv.2212.06093","url":null,"abstract":"Abstract We study a natural alternating method of Schwarz type (domain decomposition) for a certain class of couplings between local and nonlocal operators. We show that our method fits into Lions’s framework and prove, as a consequence, convergence in both the continuous and the discrete settings.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.3,"publicationDate":"2022-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45441883","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}
引用次数: 1
Positivity-Preserving Numerical Method for a Stochastic Multi-Group SIR Epidemic Model 随机多群SIR流行病模型的保正数值方法
IF 1.3 4区 数学 Q3 MATHEMATICS, APPLIED Pub Date : 2022-12-07 DOI: 10.1515/cmam-2022-0143
Han Ma, Qimin Zhang, X. Xu
Abstract The stochastic multi-group susceptible–infected–recovered (SIR) epidemic model is nonlinear, and so analytical solutions are generally difficult to obtain. Hence, it is often necessary to find numerical solutions, but most existing numerical methods fail to preserve the nonnegativity or positivity of solutions. Therefore, an appropriate numerical method for studying the dynamic behavior of epidemic diseases through SIR models is urgently required. In this paper, based on the Euler–Maruyama scheme and a logarithmic transformation, we propose a novel explicit positivity-preserving numerical scheme for a stochastic multi-group SIR epidemic model whose coefficients violate the global monotonicity condition. This scheme not only results in numerical solutions that preserve the domain of the stochastic multi-group SIR epidemic model, but also achieves the “ order - 1 2 {mathrm{order}-frac{1}{2}} ” strong convergence rate. Taking a two-group SIR epidemic model as an example, some numerical simulations are performed to illustrate the performance of the proposed scheme.
随机多群体易感感染恢复(SIR)流行病模型是非线性的,通常难以得到解析解。因此,通常需要寻找数值解,但大多数现有的数值方法都不能保持解的非负性或正性。因此,迫切需要一种合适的通过SIR模型研究传染病动力学行为的数值方法。本文基于Euler-Maruyama格式和对数变换,对系数违反全局单调性条件的随机多群SIR流行病模型,提出了一种新的显式保正数值格式。该方案不仅得到了保持随机多群SIR流行病模型定域的数值解,而且实现了“order - 1 2 { mathm {order}-frac{1}{2}}”的强收敛速率。以两组SIR流行病模型为例,进行了数值模拟,验证了该方法的有效性。
{"title":"Positivity-Preserving Numerical Method for a Stochastic Multi-Group SIR Epidemic Model","authors":"Han Ma, Qimin Zhang, X. Xu","doi":"10.1515/cmam-2022-0143","DOIUrl":"https://doi.org/10.1515/cmam-2022-0143","url":null,"abstract":"Abstract The stochastic multi-group susceptible–infected–recovered (SIR) epidemic model is nonlinear, and so analytical solutions are generally difficult to obtain. Hence, it is often necessary to find numerical solutions, but most existing numerical methods fail to preserve the nonnegativity or positivity of solutions. Therefore, an appropriate numerical method for studying the dynamic behavior of epidemic diseases through SIR models is urgently required. In this paper, based on the Euler–Maruyama scheme and a logarithmic transformation, we propose a novel explicit positivity-preserving numerical scheme for a stochastic multi-group SIR epidemic model whose coefficients violate the global monotonicity condition. This scheme not only results in numerical solutions that preserve the domain of the stochastic multi-group SIR epidemic model, but also achieves the “ order - 1 2 {mathrm{order}-frac{1}{2}} ” strong convergence rate. Taking a two-group SIR epidemic model as an example, some numerical simulations are performed to illustrate the performance of the proposed scheme.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"23 1","pages":"671 - 694"},"PeriodicalIF":1.3,"publicationDate":"2022-12-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42594497","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}
引用次数: 0
Local Discontinuous Galerkin Method for a Third-Order Singularly Perturbed Problem of Convection-Diffusion Type 一类三阶对流扩散型奇摄动问题的局部不连续伽辽金方法
IF 1.3 4区 数学 Q3 MATHEMATICS, APPLIED Pub Date : 2022-12-06 DOI: 10.1515/cmam-2022-0176
Li Yan, Zhoufeng Wang, Yao Cheng
The local discontinuous Galerkin (LDG) method is studied for a third-order singularly perturbed problem of convection-diffusion type. Based on a regularity assumption for the exact solution, we prove almost O ( N - ( k + 1 2 ) ) {O(N^{-(k+frac{1}{2})})} (up to a logarithmic factor) energy-norm convergence uniformly in the perturbation parameter. Here, k 0 {kgeq 0} is the maximum degree of piecewise polynomials used in discrete space, and N is the number of mesh elements. The results are valid for the three types of layer-adapted meshes: Shishkin-type, Bakhvalov–Shishkin-type, and Bakhvalov-type. Numerical experiments are conducted to test the theoretical results.
研究了一类三阶对流扩散型奇异摄动问题的局部不连续伽辽金方法。基于精确解的正则性假设,我们证明了几乎O(N -(k+ 1 2)) {O(N^{-(k+ frac{1}{2}))}(直到一个对数因子)能量范数在扰动参数上一致收敛。其中,k≥0 }k{geq 0为}离散空间中使用分段多项式的最大程度,N为网格单元个数。结果适用于三种类型的层适应网格:shishkin型、bakhvalov - shishkin型和bakhvalov型。数值实验对理论结果进行了验证。
{"title":"Local Discontinuous Galerkin Method for a Third-Order Singularly Perturbed Problem of Convection-Diffusion Type","authors":"Li Yan, Zhoufeng Wang, Yao Cheng","doi":"10.1515/cmam-2022-0176","DOIUrl":"https://doi.org/10.1515/cmam-2022-0176","url":null,"abstract":"The local discontinuous Galerkin (LDG) method is studied for a third-order singularly perturbed problem of convection-diffusion type. Based on a regularity assumption for the exact solution, we prove almost <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:mrow> <m:mi>O</m:mi> <m:mo>⁢</m:mo> <m:mrow> <m:mo stretchy=\"false\">(</m:mo> <m:msup> <m:mi>N</m:mi> <m:mrow> <m:mo>-</m:mo> <m:mrow> <m:mo stretchy=\"false\">(</m:mo> <m:mrow> <m:mi>k</m:mi> <m:mo>+</m:mo> <m:mfrac> <m:mn>1</m:mn> <m:mn>2</m:mn> </m:mfrac> </m:mrow> <m:mo stretchy=\"false\">)</m:mo> </m:mrow> </m:mrow> </m:msup> <m:mo stretchy=\"false\">)</m:mo> </m:mrow> </m:mrow> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2022-0176_eq_0280.png\" /> <jats:tex-math>{O(N^{-(k+frac{1}{2})})}</jats:tex-math> </jats:alternatives> </jats:inline-formula> (up to a logarithmic factor) energy-norm convergence uniformly in the perturbation parameter. Here, <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:mrow> <m:mi>k</m:mi> <m:mo>≥</m:mo> <m:mn>0</m:mn> </m:mrow> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2022-0176_eq_0414.png\" /> <jats:tex-math>{kgeq 0}</jats:tex-math> </jats:alternatives> </jats:inline-formula> is the maximum degree of piecewise polynomials used in discrete space, and <jats:italic>N</jats:italic> is the number of mesh elements. The results are valid for the three types of layer-adapted meshes: Shishkin-type, Bakhvalov–Shishkin-type, and Bakhvalov-type. Numerical experiments are conducted to test the theoretical results.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"11 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2022-12-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138537980","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}
引用次数: 0
A Posteriori Error Estimator for Weak Galerkin Finite Element Method for Stokes Problem Using Diagonalization Techniques 基于对角化技术的Stokes问题弱Galerkin有限元后验误差估计
IF 1.3 4区 数学 Q3 MATHEMATICS, APPLIED Pub Date : 2022-12-06 DOI: 10.1515/cmam-2022-0087
Jiachuan Zhang, Ran Zhang, Jingzhi Li
Abstract Based on a hierarchical basis a posteriori error estimator, an adaptive weak Galerkin finite element method (WGFEM) is proposed for the Stokes problem in two and three dimensions. In this paper, we propose two novel diagonalization techniques for velocity and pressure, respectively. Using diagonalization techniques, we need only to solve two diagonal linear algebraic systems corresponding to the degree of freedom to get the error estimator. The upper bound and lower bound of the error estimator are also shown to address the reliability of the adaptive method. Numerical simulations are provided to demonstrate the effectiveness and robustness of our algorithm.
摘要基于层次基后验误差估计器,提出了一种求解二维和三维Stokes问题的自适应弱Galerkin有限元方法。在本文中,我们分别提出了两种新的速度和压力对角化技术。利用对角化技术,我们只需要求解两个与自由度相对应的对角线性代数系统,就可以得到误差估计器。误差估计器的上界和下界也被示出,以解决自适应方法的可靠性问题。通过数值模拟验证了算法的有效性和鲁棒性。
{"title":"A Posteriori Error Estimator for Weak Galerkin Finite Element Method for Stokes Problem Using Diagonalization Techniques","authors":"Jiachuan Zhang, Ran Zhang, Jingzhi Li","doi":"10.1515/cmam-2022-0087","DOIUrl":"https://doi.org/10.1515/cmam-2022-0087","url":null,"abstract":"Abstract Based on a hierarchical basis a posteriori error estimator, an adaptive weak Galerkin finite element method (WGFEM) is proposed for the Stokes problem in two and three dimensions. In this paper, we propose two novel diagonalization techniques for velocity and pressure, respectively. Using diagonalization techniques, we need only to solve two diagonal linear algebraic systems corresponding to the degree of freedom to get the error estimator. The upper bound and lower bound of the error estimator are also shown to address the reliability of the adaptive method. Numerical simulations are provided to demonstrate the effectiveness and robustness of our algorithm.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"23 1","pages":"783 - 811"},"PeriodicalIF":1.3,"publicationDate":"2022-12-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45425406","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}
引用次数: 0
A Time-Adaptive Space-Time FMM for the Heat Equation 热方程的时间自适应时空FMM
IF 1.3 4区 数学 Q3 MATHEMATICS, APPLIED Pub Date : 2022-12-06 DOI: 10.1515/cmam-2022-0117
R. Watschinger, G. Of
Abstract We present a new time-adaptive FMM for a space-time boundary element method for the heat equation. The method extends the existing parabolic FMM by adding new operations that allow for an efficient treatment of tensor product meshes which are adaptive in time. We analyze the efficiency of the new operations and the approximation quality of the related kernel expansions and present numerical experiments that reveal the benefits of the new method.
摘要我们为热方程的时空边界元方法提出了一种新的时间自适应FMM。该方法通过添加新的运算来扩展现有的抛物型FMM,这些运算允许对时间自适应的张量积网格进行有效处理。我们分析了新运算的效率和相关核展开的近似质量,并通过数值实验揭示了新方法的优点。
{"title":"A Time-Adaptive Space-Time FMM for the Heat Equation","authors":"R. Watschinger, G. Of","doi":"10.1515/cmam-2022-0117","DOIUrl":"https://doi.org/10.1515/cmam-2022-0117","url":null,"abstract":"Abstract We present a new time-adaptive FMM for a space-time boundary element method for the heat equation. The method extends the existing parabolic FMM by adding new operations that allow for an efficient treatment of tensor product meshes which are adaptive in time. We analyze the efficiency of the new operations and the approximation quality of the related kernel expansions and present numerical experiments that reveal the benefits of the new method.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"23 1","pages":"445 - 471"},"PeriodicalIF":1.3,"publicationDate":"2022-12-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48073717","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}
引用次数: 0
A Cost-Efficient Space-Time Adaptive Algorithm for Coupled Flow and Transport 一种具有成本效益的流与输运耦合时空自适应算法
IF 1.3 4区 数学 Q3 MATHEMATICS, APPLIED Pub Date : 2022-12-06 DOI: 10.48550/arXiv.2212.02954
Marius Paul Bruchhäuser, M. Bause
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.
摘要本文针对流动和对流主导传输的耦合模型问题,提出并研究了一种基于双加权残差(DWR)方法的高效时空自适应算法。关键因素是一种适用于子问题随时间变化的动力学的多速率方法,分别用于传输和流动问题的加权和非加权误差指标,以及用于数据结构的基于张量积空间的时空板的概念。在数值例子中,研究了底层算法在基准问题和实际应用中的性能。此外,还研究了强对流主导输运的稳定性和面向目标的自适应性的相互作用。
{"title":"A Cost-Efficient Space-Time Adaptive Algorithm for Coupled Flow and Transport","authors":"Marius Paul Bruchhäuser, M. Bause","doi":"10.48550/arXiv.2212.02954","DOIUrl":"https://doi.org/10.48550/arXiv.2212.02954","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":" ","pages":""},"PeriodicalIF":1.3,"publicationDate":"2022-12-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42984459","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}
引用次数: 1
Space-Time Approximation of Local Strong Solutions to the 3D Stochastic Navier–Stokes Equations 三维随机Navier-Stokes方程局部强解的时空逼近
IF 1.3 4区 数学 Q3 MATHEMATICS, APPLIED Pub Date : 2022-11-30 DOI: 10.48550/arXiv.2211.17011
D. Breit, Alan Dodgson
Abstract We consider the 3D stochastic Navier–Stokes equation on the torus. Our main result concerns the temporal and spatio-temporal discretisation of a local strong pathwise solution. We prove optimal convergence rates for the energy error with respect to convergence in probability, that is convergence of order (up to) 1 in space and of order (up to) 1/2 in time. The result holds up to the possible blow-up of the (time-discrete) solution. Our approach is based on discrete stopping times for the (time-discrete) solution.
摘要我们考虑环面上的三维随机Navier-Stokes方程。我们的主要结果涉及局部强路径解的时间和时空离散化。我们证明了能量误差相对于概率收敛的最优收敛速度,即在空间上(高达)1阶的收敛和在时间上(高至)1/2阶的收敛。这一结果证明了(时间离散的)解可能会爆炸。我们的方法基于(时间离散的)解的离散停止时间。
{"title":"Space-Time Approximation of Local Strong Solutions to the 3D Stochastic Navier–Stokes Equations","authors":"D. Breit, Alan Dodgson","doi":"10.48550/arXiv.2211.17011","DOIUrl":"https://doi.org/10.48550/arXiv.2211.17011","url":null,"abstract":"Abstract We consider the 3D stochastic Navier–Stokes equation on the torus. Our main result concerns the temporal and spatio-temporal discretisation of a local strong pathwise solution. We prove optimal convergence rates for the energy error with respect to convergence in probability, that is convergence of order (up to) 1 in space and of order (up to) 1/2 in time. The result holds up to the possible blow-up of the (time-discrete) solution. Our approach is based on discrete stopping times for the (time-discrete) solution.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.3,"publicationDate":"2022-11-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"47859877","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}
引用次数: 0
Bernoulli Free Boundary Problems Under Uncertainty: The Convex Case 不确定条件下的Bernoulli自由边界问题:凸情形
IF 1.3 4区 数学 Q3 MATHEMATICS, APPLIED Pub Date : 2022-11-24 DOI: 10.1515/cmam-2022-0038
M. Dambrine, H. Harbrecht, B. Puig
Abstract The present article is concerned with solving Bernoulli’s exterior free boundary problem in the case of an interior boundary that is random. We provide a new regularity result on the map that sends a parametrization of the inner boundary to a parametrization of the outer boundary. Moreover, assuming that the interior boundary is convex, also the exterior boundary is convex, which enables to identify the boundaries with support functions and to determine their expectations. We in particular construct a confidence region for the outer boundary based on Aumann’s expectation and provide a numerical method to compute it.
摘要本文研究了随机内边界情况下的伯努利外自由边界问题的求解。我们在映射上提供了一个新的正则性结果,将内边界的参数化发送到外边界的参数化。此外,假设内部边界是凸的,那么外部边界也是凸的,从而可以识别具有支持函数的边界,并确定其期望。特别地,我们基于Aumann期望构造了外边界的置信区域,并给出了计算置信区域的数值方法。
{"title":"Bernoulli Free Boundary Problems Under Uncertainty: The Convex Case","authors":"M. Dambrine, H. Harbrecht, B. Puig","doi":"10.1515/cmam-2022-0038","DOIUrl":"https://doi.org/10.1515/cmam-2022-0038","url":null,"abstract":"Abstract The present article is concerned with solving Bernoulli’s exterior free boundary problem in the case of an interior boundary that is random. We provide a new regularity result on the map that sends a parametrization of the inner boundary to a parametrization of the outer boundary. Moreover, assuming that the interior boundary is convex, also the exterior boundary is convex, which enables to identify the boundaries with support functions and to determine their expectations. We in particular construct a confidence region for the outer boundary based on Aumann’s expectation and provide a numerical method to compute it.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"23 1","pages":"333 - 352"},"PeriodicalIF":1.3,"publicationDate":"2022-11-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41622520","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}
引用次数: 0
Computational Multiscale Methods for Nondivergence-Form Elliptic Partial Differential Equations 非发散型椭圆型偏微分方程的计算多尺度方法
IF 1.3 4区 数学 Q3 MATHEMATICS, APPLIED Pub Date : 2022-11-24 DOI: 10.48550/arXiv.2211.13731
P. Freese, D. Gallistl, D. Peterseim, Timo Sprekeler
Abstract This paper proposes novel computational multiscale methods for linear second-order elliptic partial differential equations in nondivergence form with heterogeneous coefficients satisfying a Cordes condition. The construction follows the methodology of localized orthogonal decomposition (LOD) and provides operator-adapted coarse spaces by solving localized cell problems on a fine scale in the spirit of numerical homogenization. The degrees of freedom of the coarse spaces are related to nonconforming and mixed finite element methods for homogeneous problems. The rigorous error analysis of one exemplary approach shows that the favorable properties of the LOD methodology known from divergence-form PDEs, i.e., its applicability and accuracy beyond scale separation and periodicity, remain valid for problems in nondivergence form.
摘要本文提出了一种新的计算多尺度方法,用于求解具有满足Cordes条件的异质系数的非微分形式的线性二阶椭圆型偏微分方程。该构造遵循局部正交分解(LOD)的方法,并通过在数值均匀化的精神下在精细尺度上解决局部单元问题来提供算子自适应的粗略空间。粗空间的自由度与齐次问题的非协调和混合有限元方法有关。一种示例性方法的严格误差分析表明,从发散形式的偏微分方程中已知的LOD方法的有利性质,即其超出尺度分离和周期性的适用性和准确性,对于非发散形式的问题仍然有效。
{"title":"Computational Multiscale Methods for Nondivergence-Form Elliptic Partial Differential Equations","authors":"P. Freese, D. Gallistl, D. Peterseim, Timo Sprekeler","doi":"10.48550/arXiv.2211.13731","DOIUrl":"https://doi.org/10.48550/arXiv.2211.13731","url":null,"abstract":"Abstract This paper proposes novel computational multiscale methods for linear second-order elliptic partial differential equations in nondivergence form with heterogeneous coefficients satisfying a Cordes condition. The construction follows the methodology of localized orthogonal decomposition (LOD) and provides operator-adapted coarse spaces by solving localized cell problems on a fine scale in the spirit of numerical homogenization. The degrees of freedom of the coarse spaces are related to nonconforming and mixed finite element methods for homogeneous problems. The rigorous error analysis of one exemplary approach shows that the favorable properties of the LOD methodology known from divergence-form PDEs, i.e., its applicability and accuracy beyond scale separation and periodicity, remain valid for problems in nondivergence form.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.3,"publicationDate":"2022-11-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43778213","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}
引用次数: 1
Finite Element Approximations for PDEs with Irregular Dirichlet Boundary Data on Boundary Concentrated Meshes 边界集中网格上具有不规则Dirichlet边界数据的偏微分方程的有限元逼近
IF 1.3 4区 数学 Q3 MATHEMATICS, APPLIED Pub Date : 2022-11-11 DOI: 10.1515/cmam-2022-0129
J. Pfefferer, M. Winkler
Abstract This paper is concerned with finite element error estimates for second order elliptic PDEs with inhomogeneous Dirichlet boundary data in convex polygonal domains. The Dirichlet boundary data are assumed to be irregular such that the solution of the PDE does not belong to H 2 ⁢ ( Ω ) {H^{2}(Omega)} but only to H r ⁢ ( Ω ) {H^{r}(Omega)} for some r ∈ ( 1 , 2 ) {rin(1,2)} . As a consequence, a discretization of the PDE with linear finite elements exhibits a reduced convergence rate in L 2 ⁢ ( Ω ) {L^{2}(Omega)} and H 1 ⁢ ( Ω ) {H^{1}(Omega)} . In order to restore the best possible convergence rate we propose and analyze in detail the usage of boundary concentrated meshes. These meshes are gradually refined towards the whole boundary. The corresponding grading parameter does not only depend on the regularity of the Dirichlet boundary data and their discrete implementation but also on the norm, which is used to measure the error. In numerical experiments we confirm our theoretical results.
摘要本文研究凸多边形域中具有非齐次Dirichlet边界数据的二阶椭圆偏微分方程的有限元误差估计。假设Dirichlet边界数据是不规则的,使得PDE的解不属于H2(Ω){H^{2}(Omega)},而只属于某些r∈(1,2){r(1,2)}的HR(Ω)}{H^{r}( Omega)}。因此,具有线性有限元的PDE的离散化在L2(Ω){L^{2}(Omega)}和H1(Ω){H^{1}(Omega)}中表现出降低的收敛速度。为了恢复最佳的收敛速度,我们提出并详细分析了边界集中网格的使用。这些网格会朝着整个边界逐渐细化。相应的分级参数不仅取决于狄利克雷边界数据的正则性及其离散实现,还取决于用于测量误差的范数。在数值实验中,我们证实了我们的理论结果。
{"title":"Finite Element Approximations for PDEs with Irregular Dirichlet Boundary Data on Boundary Concentrated Meshes","authors":"J. Pfefferer, M. Winkler","doi":"10.1515/cmam-2022-0129","DOIUrl":"https://doi.org/10.1515/cmam-2022-0129","url":null,"abstract":"Abstract This paper is concerned with finite element error estimates for second order elliptic PDEs with inhomogeneous Dirichlet boundary data in convex polygonal domains. The Dirichlet boundary data are assumed to be irregular such that the solution of the PDE does not belong to H 2 ⁢ ( Ω ) {H^{2}(Omega)} but only to H r ⁢ ( Ω ) {H^{r}(Omega)} for some r ∈ ( 1 , 2 ) {rin(1,2)} . As a consequence, a discretization of the PDE with linear finite elements exhibits a reduced convergence rate in L 2 ⁢ ( Ω ) {L^{2}(Omega)} and H 1 ⁢ ( Ω ) {H^{1}(Omega)} . In order to restore the best possible convergence rate we propose and analyze in detail the usage of boundary concentrated meshes. These meshes are gradually refined towards the whole boundary. The corresponding grading parameter does not only depend on the regularity of the Dirichlet boundary data and their discrete implementation but also on the norm, which is used to measure the error. In numerical experiments we confirm our theoretical results.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.3,"publicationDate":"2022-11-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42712860","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}
引用次数: 0
期刊
Computational Methods in Applied Mathematics
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1