首页 > 最新文献

IMA Journal of Numerical Analysis最新文献

英文 中文
Discrete Gagliardo–Nirenberg inequality and application to the finite volume approximation of a convection–diffusion equation with a Joule effect term 离散gagliado - nirenberg不等式及其在具有焦耳效应项的对流扩散方程有限体积近似中的应用
IF 2.1 2区 数学 Q1 Mathematics Pub Date : 2023-08-30 DOI: 10.1093/imanum/drad063
C. Calgaro, C. Cancès, E. Creusé
A discrete order-two Gagliardo–Nirenberg inequality is established for piecewise constant functions defined on a two-dimensional structured mesh composed of rectangular cells. As in the continuous framework, this discrete Gagliardo–Nirenberg inequality allows to control in particular the $L^4$ norm of the discrete gradient of the numerical solution by the $L^2$ norm of its discrete Hessian times its $L^infty $ norm. This result is crucial for the convergence analysis of a finite volume method for the approximation of a convection–diffusion equation involving a Joule effect term on a uniform mesh in each direction. The convergence proof relies on compactness arguments and on a priori estimates under a smallness assumption on the data, which is essential also in the continuous framework.
对于定义在矩形单元组成的二维结构网格上的分段常数函数,建立了离散二阶Gagliardo-Nirenberg不等式。就像在连续框架中一样,这个离散的伽利亚多-尼伦伯格不等式允许通过其离散的海森的$L^2$范数乘以其$L^infty $范数来控制数值解的离散梯度的$L^4$范数。这一结果对于用有限体积法逼近每个方向均匀网格上包含焦耳效应项的对流扩散方程的收敛性分析是至关重要的。收敛性证明依赖于紧性参数和数据在小假设下的先验估计,这在连续框架中也是必不可少的。
{"title":"Discrete Gagliardo–Nirenberg inequality and application to the finite volume approximation of a convection–diffusion equation with a Joule effect term","authors":"C. Calgaro, C. Cancès, E. Creusé","doi":"10.1093/imanum/drad063","DOIUrl":"https://doi.org/10.1093/imanum/drad063","url":null,"abstract":"A discrete order-two Gagliardo–Nirenberg inequality is established for piecewise constant functions defined on a two-dimensional structured mesh composed of rectangular cells. As in the continuous framework, this discrete Gagliardo–Nirenberg inequality allows to control in particular the $L^4$ norm of the discrete gradient of the numerical solution by the $L^2$ norm of its discrete Hessian times its $L^infty $ norm. This result is crucial for the convergence analysis of a finite volume method for the approximation of a convection–diffusion equation involving a Joule effect term on a uniform mesh in each direction. The convergence proof relies on compactness arguments and on a priori estimates under a smallness assumption on the data, which is essential also in the continuous framework.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2023-08-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48381868","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
An arbitrary-order discrete rot-rot complex on polygonal meshes with application to a quad-rot problem 多边形网格上任意阶离散rot-rot复形及其在四次rot问题中的应用
2区 数学 Q1 Mathematics Pub Date : 2023-08-28 DOI: 10.1093/imanum/drad045
Daniele Antonio Di Pietro
Abstract In this work, following the discrete de Rham approach, we develop a discrete counterpart of a two-dimensional de Rham complex with enhanced regularity. The proposed construction supports general polygonal meshes and arbitrary approximation orders. We establish exactness on a contractible domain for both the versions of the complex with and without boundary conditions and, for the former, prove a complete set of Poincaré-type inequalities. The discrete complex is then used to derive a novel discretization method for a quad-rot problem, which, unlike other schemes in the literature, does not require the forcing term to be prepared. We carry out complete stability and convergence analyses for the proposed scheme and provide numerical validation of the results.
在这项工作中,遵循离散de Rham方法,我们开发了具有增强规律性的二维de Rham复合体的离散对应物。该结构支持一般多边形网格和任意近似顺序。对于有边界条件和无边界条件的复合体,我们在可缩域上建立了精确性,对于有边界条件的复合体,我们证明了一组完备的poincar型不等式。然后利用离散复形导出了一种新的四次问题的离散化方法,与文献中的其他方案不同,该方法不需要准备强迫项。我们对所提出的方案进行了完整的稳定性和收敛性分析,并对结果进行了数值验证。
{"title":"An arbitrary-order discrete rot-rot complex on polygonal meshes with application to a quad-rot problem","authors":"Daniele Antonio Di Pietro","doi":"10.1093/imanum/drad045","DOIUrl":"https://doi.org/10.1093/imanum/drad045","url":null,"abstract":"Abstract In this work, following the discrete de Rham approach, we develop a discrete counterpart of a two-dimensional de Rham complex with enhanced regularity. The proposed construction supports general polygonal meshes and arbitrary approximation orders. We establish exactness on a contractible domain for both the versions of the complex with and without boundary conditions and, for the former, prove a complete set of Poincaré-type inequalities. The discrete complex is then used to derive a novel discretization method for a quad-rot problem, which, unlike other schemes in the literature, does not require the forcing term to be prepared. We carry out complete stability and convergence analyses for the proposed scheme and provide numerical validation of the results.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-08-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136348995","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Correction to: On the stability of totally upwind schemes for the hyperbolic initial boundary value problem 修正:关于双曲型初边值问题的全逆风格式的稳定性
IF 2.1 2区 数学 Q1 Mathematics Pub Date : 2023-08-26 DOI: 10.1093/imanum/drad068
{"title":"Correction to: On the stability of totally upwind schemes for the hyperbolic initial boundary value problem","authors":"","doi":"10.1093/imanum/drad068","DOIUrl":"https://doi.org/10.1093/imanum/drad068","url":null,"abstract":"","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2023-08-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48242424","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A nodally bound-preserving finite element method 节点保界有限元方法
2区 数学 Q1 Mathematics Pub Date : 2023-08-26 DOI: 10.1093/imanum/drad055
Gabriel R Barrenechea, Emmanuil H Georgoulis, Tristan Pryer, Andreas Veeser
Abstract This work proposes a nonlinear finite element method whose nodal values preserve bounds known for the exact solution. The discrete problem involves a nonlinear projection operator mapping arbitrary nodal values into bound-preserving ones and seeks the numerical solution in the range of this projection. As the projection is not injective, a stabilisation based upon the complementary projection is added in order to restore well-posedness. Within the framework of elliptic problems, the discrete problem may be viewed as a reformulation of a discrete obstacle problem, incorporating the inequality constraints through Lipschitz projections. The derivation of the proposed method is exemplified for linear and nonlinear reaction-diffusion problems. Near-best approximation results in suitable norms are established. In particular, we prove that, in the linear case, the numerical solution is the best approximation in the energy norm among all nodally bound-preserving finite element functions. A series of numerical experiments for such problems showcase the good behaviour of the proposed bound-preserving finite element method.
摘要本文提出了一种非线性有限元方法,其节点值保持精确解的已知界。离散问题涉及一个非线性投影算子将任意节点值映射为保界节点值,并在该投影范围内寻求数值解。由于投影不是内射,为了恢复适位性,在互补投影的基础上增加了一个稳定化。在椭圆问题的框架内,离散问题可以看作是一个离散障碍问题的重新表述,通过Lipschitz投影纳入不等式约束。对线性和非线性反应扩散问题给出了该方法的推导。在合适的规范下建立了近似的最佳近似结果。特别地,我们证明了在线性情况下,数值解是所有节点保界有限元函数在能量范数上的最佳逼近。对这类问题的一系列数值实验表明了所提出的保界有限元方法的良好性能。
{"title":"A nodally bound-preserving finite element method","authors":"Gabriel R Barrenechea, Emmanuil H Georgoulis, Tristan Pryer, Andreas Veeser","doi":"10.1093/imanum/drad055","DOIUrl":"https://doi.org/10.1093/imanum/drad055","url":null,"abstract":"Abstract This work proposes a nonlinear finite element method whose nodal values preserve bounds known for the exact solution. The discrete problem involves a nonlinear projection operator mapping arbitrary nodal values into bound-preserving ones and seeks the numerical solution in the range of this projection. As the projection is not injective, a stabilisation based upon the complementary projection is added in order to restore well-posedness. Within the framework of elliptic problems, the discrete problem may be viewed as a reformulation of a discrete obstacle problem, incorporating the inequality constraints through Lipschitz projections. The derivation of the proposed method is exemplified for linear and nonlinear reaction-diffusion problems. Near-best approximation results in suitable norms are established. In particular, we prove that, in the linear case, the numerical solution is the best approximation in the energy norm among all nodally bound-preserving finite element functions. A series of numerical experiments for such problems showcase the good behaviour of the proposed bound-preserving finite element method.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-08-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135236854","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
High order approximations of the Cox–Ingersoll–Ross process semigroup using random grids Cox-Ingersoll-Ross过程半群的随机网格高阶逼近
2区 数学 Q1 Mathematics Pub Date : 2023-08-24 DOI: 10.1093/imanum/drad059
Aurélien Alfonsi, Edoardo Lombardo
Abstract We present new high order approximations schemes for the Cox–Ingersoll–Ross (CIR) process that are obtained by using a recent technique developed by Alfonsi and Bally (2021, A generic construction for high order approximation schemes of semigroups using random grids. Numer. Math., 148, 743–793) for the approximation of semigroups. The idea consists in using a suitable combination of discretization schemes calculated on different random grids to increase the order of convergence. This technique coupled with the second order scheme proposed by Alfonsi (2010, High order discretization schemes for the CIR process: application to affine term structure and Heston models. Math. Comp., 79, 209–237) for the CIR leads to weak approximations of order $2k$, for all $kin{{mathbb{N}}}^{ast }$. Despite the singularity of the square-root volatility coefficient, we show rigorously this order of convergence under some restrictions on the volatility parameters. We illustrate numerically the convergence of these approximations for the CIR process and for the Heston stochastic volatility model and show the computational time gain they give.
本文提出了一种新的Cox-Ingersoll-Ross (CIR)过程的高阶近似格式,该格式是利用Alfonsi和Bally(2021)最近开发的技术获得的,这是一种使用随机网格的半群高阶近似格式的一般构造。号码。数学。半群的近似。数学学报,14,743-793)。其思想在于使用在不同随机网格上计算的离散化方案的适当组合来提高收敛阶。该技术与Alfonsi(2010)提出的二阶方案相结合,CIR过程的高阶离散化方案:应用于仿射期限结构和Heston模型。数学。对于{{mathbb{N}}}^{ast}$中的所有$k, Comp., 79, 209-237)对于CIR的弱近似为$2k$。尽管平方根波动系数具有奇异性,但在波动参数的某些限制下,我们严格地证明了这种收敛顺序。我们用数值说明了CIR过程和赫斯顿随机波动模型的这些近似的收敛性,并显示了它们给出的计算时间增益。
{"title":"High order approximations of the Cox–Ingersoll–Ross process semigroup using random grids","authors":"Aurélien Alfonsi, Edoardo Lombardo","doi":"10.1093/imanum/drad059","DOIUrl":"https://doi.org/10.1093/imanum/drad059","url":null,"abstract":"Abstract We present new high order approximations schemes for the Cox–Ingersoll–Ross (CIR) process that are obtained by using a recent technique developed by Alfonsi and Bally (2021, A generic construction for high order approximation schemes of semigroups using random grids. Numer. Math., 148, 743–793) for the approximation of semigroups. The idea consists in using a suitable combination of discretization schemes calculated on different random grids to increase the order of convergence. This technique coupled with the second order scheme proposed by Alfonsi (2010, High order discretization schemes for the CIR process: application to affine term structure and Heston models. Math. Comp., 79, 209–237) for the CIR leads to weak approximations of order $2k$, for all $kin{{mathbb{N}}}^{ast }$. Despite the singularity of the square-root volatility coefficient, we show rigorously this order of convergence under some restrictions on the volatility parameters. We illustrate numerically the convergence of these approximations for the CIR process and for the Heston stochastic volatility model and show the computational time gain they give.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-08-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"134983787","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Numerical approximation of singular-degenerate parabolic stochastic partial differential equations 奇异退化抛物型随机偏微分方程的数值逼近
IF 2.1 2区 数学 Q1 Mathematics Pub Date : 2023-08-24 DOI: 10.1093/imanum/drad061
L. Baňas, B. Gess, C. Vieth
We study a general class of singular degenerate parabolic stochastic partial differential equations (SPDEs) that include, in particular, the stochastic porous medium equations and the stochastic fast diffusion equation. We propose a fully discrete numerical approximation of the considered SPDEs based on the very weak formulation. By exploiting the monotonicity properties of the proposed formulation we prove the convergence of the numerical approximation towards the unique solution. Furthermore, we construct an implementable finite element scheme for the spatial discretization of the very weak formulation and provide numerical simulations to demonstrate the practicability of the proposed discretization.
研究了一类广义的退化抛物型随机偏微分方程(SPDEs),主要包括随机多孔介质方程和随机快速扩散方程。我们提出了基于非常弱公式的考虑的spde的完全离散数值近似。利用所提公式的单调性,证明了数值逼近对唯一解的收敛性。此外,我们构建了一个可实现的非常弱公式空间离散化的有限元方案,并提供数值模拟来证明所提出的离散化的实用性。
{"title":"Numerical approximation of singular-degenerate parabolic stochastic partial differential equations","authors":"L. Baňas, B. Gess, C. Vieth","doi":"10.1093/imanum/drad061","DOIUrl":"https://doi.org/10.1093/imanum/drad061","url":null,"abstract":"\u0000 We study a general class of singular degenerate parabolic stochastic partial differential equations (SPDEs) that include, in particular, the stochastic porous medium equations and the stochastic fast diffusion equation. We propose a fully discrete numerical approximation of the considered SPDEs based on the very weak formulation. By exploiting the monotonicity properties of the proposed formulation we prove the convergence of the numerical approximation towards the unique solution. Furthermore, we construct an implementable finite element scheme for the spatial discretization of the very weak formulation and provide numerical simulations to demonstrate the practicability of the proposed discretization.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2023-08-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44342261","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A second-order bulk–surface splitting for parabolic problems with dynamic boundary conditions 具有动态边界条件的抛物型问题的二阶体面分裂
2区 数学 Q1 Mathematics Pub Date : 2023-08-12 DOI: 10.1093/imanum/drad062
Robert Altmann, Christoph Zimmer
Abstract This paper introduces a novel approach for the construction of bulk–surface splitting schemes for semilinear parabolic partial differential equations with dynamic boundary conditions. The proposed construction is based on a reformulation of the system as a partial differential–algebraic equation and the inclusion of certain delay terms for the decoupling. To obtain a fully discrete scheme, the splitting approach is combined with finite elements in space and a backward differentiation formula in time. Within this paper, we focus on the second-order case, resulting in a $3$-step scheme. We prove second-order convergence under the assumption of a weak CFL-type condition and confirm the theoretical findings by numerical experiments. Moreover, we illustrate the potential for higher-order splitting schemes numerically.
摘要本文介绍了一种构造具有动态边界条件的半线性抛物型偏微分方程体面分裂格式的新方法。所提出的构造是基于将系统重新表述为偏微分代数方程,并包含解耦的某些延迟项。为了得到完全离散格式,将空间上的有限元和时间上的后向微分公式相结合。在本文中,我们关注二阶情况,得到一个$3$步方案。在弱cfl型条件下证明了二阶收敛性,并通过数值实验验证了理论结果。此外,我们在数值上说明了高阶分裂方案的潜力。
{"title":"A second-order bulk–surface splitting for parabolic problems with dynamic boundary conditions","authors":"Robert Altmann, Christoph Zimmer","doi":"10.1093/imanum/drad062","DOIUrl":"https://doi.org/10.1093/imanum/drad062","url":null,"abstract":"Abstract This paper introduces a novel approach for the construction of bulk–surface splitting schemes for semilinear parabolic partial differential equations with dynamic boundary conditions. The proposed construction is based on a reformulation of the system as a partial differential–algebraic equation and the inclusion of certain delay terms for the decoupling. To obtain a fully discrete scheme, the splitting approach is combined with finite elements in space and a backward differentiation formula in time. Within this paper, we focus on the second-order case, resulting in a $3$-step scheme. We prove second-order convergence under the assumption of a weak CFL-type condition and confirm the theoretical findings by numerical experiments. Moreover, we illustrate the potential for higher-order splitting schemes numerically.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-08-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"134977702","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Optimized Schwarz methods for the time-dependent Stokes–Darcy coupling 优化了时变Stokes-Darcy耦合的Schwarz方法
2区 数学 Q1 Mathematics Pub Date : 2023-08-05 DOI: 10.1093/imanum/drad057
Marco Discacciati, Tommaso Vanzan
Abstract This paper derives optimized coefficients for optimized Schwarz iterations for the time-dependent Stokes–Darcy problem using an innovative strategy to solve a nonstandard min-max problem. The coefficients take into account both physical and discretization parameters that characterize the coupled problem, and they guarantee the robustness of the associated domain decomposition method. Numerical results validate the proposed approach in several test cases with physically relevant parameters.
摘要本文利用一种创新的策略求解非标准最小-最大问题,导出了时变Stokes-Darcy问题的最优Schwarz迭代的最优系数。该系数考虑了表征耦合问题的物理参数和离散参数,保证了相关区域分解方法的鲁棒性。数值结果验证了几种具有物理相关参数的测试用例的有效性。
{"title":"Optimized Schwarz methods for the time-dependent Stokes–Darcy coupling","authors":"Marco Discacciati, Tommaso Vanzan","doi":"10.1093/imanum/drad057","DOIUrl":"https://doi.org/10.1093/imanum/drad057","url":null,"abstract":"Abstract This paper derives optimized coefficients for optimized Schwarz iterations for the time-dependent Stokes–Darcy problem using an innovative strategy to solve a nonstandard min-max problem. The coefficients take into account both physical and discretization parameters that characterize the coupled problem, and they guarantee the robustness of the associated domain decomposition method. Numerical results validate the proposed approach in several test cases with physically relevant parameters.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-08-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136083409","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A posteriori error analysis of space-time discontinuous Galerkin methods for the ε-stochastic Allen–Cahn equation ε-随机Allen–Cahn方程时空间断Galerkin方法的后验误差分析
IF 2.1 2区 数学 Q1 Mathematics Pub Date : 2023-08-02 DOI: 10.1093/imanum/drad052
D. Antonopoulou, Bernard A. Egwu, Yubin Yan
In this work, we apply an a posteriori error analysis for the space-time, discontinuous in time, Galerkin scheme, which has been proposed in Antonopoulou (2020, Space-time discontinuous Galerkin methods for the $varepsilon $-dependent stochastic Allen–Cahn equation with mild noise. IMA J. Num. Analysis, 40, 2076–2105) for the $varepsilon $-dependent stochastic Allen–Cahn equation with mild noise $dot{W}^varepsilon $ tending to rough as $varepsilon rightarrow 0$. Our results are derived under low regularity since the noise even smooth in space is assumed only one-time continuously differentiable in time, according to the minimum regularity properties of Funaki (1999, Singular limit for stochastic reaction–diffusion equation and generation of random interfaces. Acta Math. Sinica, 15, 407–438). We prove a posteriori error estimates for the $m$-dimensional problem, $mleq 4$ for a general class of space-time finite element spaces. The a posteriori bound is growing only polynomially in $varepsilon ^{-1}$ if the step length $h$ is bounded by a positive power of $varepsilon $. This agrees with the restriction posed so far in the a priori error analysis of continuous finite element schemes for the $varepsilon $-dependent deterministic Allen–Cahn or deterministic and stochastic Cahn–Hilliard equation. As an application, we examine tensorial elements where the discrete solution is approximated by polynomial functions of separated space and time variables; the a posteriori estimates there involve dimensions, and the space, time discretization parameters. We then consider the special case of the mild noise $dot{W}^varepsilon $ as defined in Weber (2010, On the short time asymptotic of the stochastic Allen–Cahn equation. Ann. Inst. Henri Poincare Probab. Stat., 46, 965–975) through the convolution of a Gaussian process with a proper mollifying kernel, which is then numerically constructed. Finally, we provide some useful insights for the numerical algorithm, and present for the first time some numerical experiments of the scheme for both one- and two-dimensional problems in various cases of interest, and compare with the deterministic ones.
在这项工作中,我们对Antonopoulou(2020)中提出的具有轻微噪声的$varepsilon $依赖随机Allen-Cahn方程的时空不连续Galerkin格式进行了后测误差分析。李建军,李建军,李建军,等。具有噪声的$varepsilon $依赖随机Allen-Cahn方程$dot{W}^varepsilon $趋近于$varepsilon rightarrow 0$。根据Funaki(1999)的最小正则性,随机反应-扩散方程的奇异极限和随机界面的生成,我们的结果是在低正则性下得到的,因为假设噪声在空间上均匀光滑,只是在时间上一次连续可微。数学学报。中国科学,15,407-438)。我们证明了一个后验误差估计的$m$维问题,$mleq 4$对于一类一般的时空有限元空间。如果步长$h$以$varepsilon $的正幂为界,则后验界仅在$varepsilon ^{-1}$中多项式地增长。这与迄今为止对$varepsilon $依赖的确定性Allen-Cahn或确定性和随机Cahn-Hilliard方程的连续有限元格式的先验误差分析所提出的限制一致。作为一个应用,我们研究了离散解由分离空间变量和时间变量的多项式函数近似的张量元素;其中的后验估计涉及维度,以及空间、时间离散参数。然后,我们考虑Weber(2010)在随机Allen-Cahn方程的短时间渐近中定义的轻微噪声$dot{W}^varepsilon $的特殊情况。安。亨利·庞加莱研究所。Stat., 46, 965-975),通过一个高斯过程的卷积和一个适当的缓和核,然后数值构造。最后,我们为数值算法提供了一些有用的见解,并首次给出了该方案在各种感兴趣的情况下的一维和二维问题的数值实验,并与确定性问题进行了比较。
{"title":"A posteriori error analysis of space-time discontinuous Galerkin methods for the ε-stochastic Allen–Cahn equation","authors":"D. Antonopoulou, Bernard A. Egwu, Yubin Yan","doi":"10.1093/imanum/drad052","DOIUrl":"https://doi.org/10.1093/imanum/drad052","url":null,"abstract":"\u0000 In this work, we apply an a posteriori error analysis for the space-time, discontinuous in time, Galerkin scheme, which has been proposed in Antonopoulou (2020, Space-time discontinuous Galerkin methods for the $varepsilon $-dependent stochastic Allen–Cahn equation with mild noise. IMA J. Num. Analysis, 40, 2076–2105) for the $varepsilon $-dependent stochastic Allen–Cahn equation with mild noise $dot{W}^varepsilon $ tending to rough as $varepsilon rightarrow 0$. Our results are derived under low regularity since the noise even smooth in space is assumed only one-time continuously differentiable in time, according to the minimum regularity properties of Funaki (1999, Singular limit for stochastic reaction–diffusion equation and generation of random interfaces. Acta Math. Sinica, 15, 407–438). We prove a posteriori error estimates for the $m$-dimensional problem, $mleq 4$ for a general class of space-time finite element spaces. The a posteriori bound is growing only polynomially in $varepsilon ^{-1}$ if the step length $h$ is bounded by a positive power of $varepsilon $. This agrees with the restriction posed so far in the a priori error analysis of continuous finite element schemes for the $varepsilon $-dependent deterministic Allen–Cahn or deterministic and stochastic Cahn–Hilliard equation. As an application, we examine tensorial elements where the discrete solution is approximated by polynomial functions of separated space and time variables; the a posteriori estimates there involve dimensions, and the space, time discretization parameters. We then consider the special case of the mild noise $dot{W}^varepsilon $ as defined in Weber (2010, On the short time asymptotic of the stochastic Allen–Cahn equation. Ann. Inst. Henri Poincare Probab. Stat., 46, 965–975) through the convolution of a Gaussian process with a proper mollifying kernel, which is then numerically constructed. Finally, we provide some useful insights for the numerical algorithm, and present for the first time some numerical experiments of the scheme for both one- and two-dimensional problems in various cases of interest, and compare with the deterministic ones.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2023-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48726746","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Optimal numerical integration and approximation of functions on ℝd equipped with Gaussian measure 具有高斯测度的函数的最优数值积分与逼近
2区 数学 Q1 Mathematics Pub Date : 2023-08-02 DOI: 10.1093/imanum/drad051
Dinh Dũng, Van Kien Nguyen
Abstract We investigate the numerical approximation of integrals over $mathbb{R}^{d}$ equipped with the standard Gaussian measure $gamma $ for integrands belonging to the Gaussian-weighted Sobolev spaces $W^{alpha }_{p}(mathbb{R}^{d}, gamma )$ of mixed smoothness $alpha in mathbb{N}$ for $1 < p < infty $. We prove the asymptotic order of the convergence of optimal quadratures based on $n$ integration nodes and propose a novel method for constructing asymptotically optimal quadratures. As for related problems, we establish by a similar technique the asymptotic order of the linear, Kolmogorov and sampling $n$-widths in the Gaussian-weighted space $L_{q}(mathbb{R}^{d}, gamma )$ of the unit ball of $W^{alpha }_{p}(mathbb{R}^{d}, gamma )$ for $1 leq q < p < infty $ and $q=p=2$.
摘要本文研究了$mathbb{R}^{d}$上具有标准高斯测度$gamma $的混合光滑高斯加权Sobolev空间$W^{alpha }_{p}(mathbb{R}^{d}, gamma )$(对于$1 < p < infty $) $alpha in mathbb{N}$的积分的数值逼近。基于$n$积分节点证明了最优正交的渐近阶收敛性,提出了一种构造渐近最优正交的新方法。对于相关问题,我们用类似的方法建立了$1 leq q < p < infty $和$q=p=2$在单位球$W^{alpha }_{p}(mathbb{R}^{d}, gamma )$的高斯加权空间$L_{q}(mathbb{R}^{d}, gamma )$中线性宽度、Kolmogorov宽度和采样$n$ -宽度的渐近阶。
{"title":"Optimal numerical integration and approximation of functions on ℝ<i>d</i> equipped with Gaussian measure","authors":"Dinh Dũng, Van Kien Nguyen","doi":"10.1093/imanum/drad051","DOIUrl":"https://doi.org/10.1093/imanum/drad051","url":null,"abstract":"Abstract We investigate the numerical approximation of integrals over $mathbb{R}^{d}$ equipped with the standard Gaussian measure $gamma $ for integrands belonging to the Gaussian-weighted Sobolev spaces $W^{alpha }_{p}(mathbb{R}^{d}, gamma )$ of mixed smoothness $alpha in mathbb{N}$ for $1 &amp;lt; p &amp;lt; infty $. We prove the asymptotic order of the convergence of optimal quadratures based on $n$ integration nodes and propose a novel method for constructing asymptotically optimal quadratures. As for related problems, we establish by a similar technique the asymptotic order of the linear, Kolmogorov and sampling $n$-widths in the Gaussian-weighted space $L_{q}(mathbb{R}^{d}, gamma )$ of the unit ball of $W^{alpha }_{p}(mathbb{R}^{d}, gamma )$ for $1 leq q &amp;lt; p &amp;lt; infty $ and $q=p=2$.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135064454","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
期刊
IMA Journal of Numerical Analysis
全部 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