We consider the numerical approximation of Gaussian random fields on closed surfaces defined as the solution to a fractional stochastic partial differential equation (SPDE) with additive white noise. The SPDE involves two parameters controlling the smoothness and the correlation length of the Gaussian random field. The proposed numerical method relies on the Balakrishnan integral representation of the solution and does not require the approximation of eigenpairs. Rather, it consists of a sinc quadrature coupled with a standard surface finite element method. We provide a complete error analysis of the method and illustrate its performances in several numerical experiments.
{"title":"Numerical Approximation of Gaussian Random Fields on Closed Surfaces","authors":"Andrea Bonito, Diane Guignard, Wenyu Lei","doi":"10.1515/cmam-2022-0237","DOIUrl":"https://doi.org/10.1515/cmam-2022-0237","url":null,"abstract":"We consider the numerical approximation of Gaussian random fields on closed surfaces defined as the solution to a fractional stochastic partial differential equation (SPDE) with additive white noise. The SPDE involves two parameters controlling the smoothness and the correlation length of the Gaussian random field. The proposed numerical method relies on the Balakrishnan integral representation of the solution and does not require the approximation of eigenpairs. Rather, it consists of a sinc quadrature coupled with a standard surface finite element method. We provide a complete error analysis of the method and illustrate its performances in several numerical experiments.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2024-01-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139104590","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}
The problem of optimal recovering high-order mixed derivatives of bivariate functions with finite smoothness is studied. Based on the truncation method, an algorithm for numerical differentiation is constructed, which is order-optimal both in the sense of accuracy and in terms of the amount of involved Galerkin information. Numerical examples are provided to illustrate the fact that our approach can be implemented successfully.
{"title":"An Optimal Method for High-Order Mixed Derivatives of Bivariate Functions","authors":"Evgeniya V. Semenova, Sergiy G. Solodky","doi":"10.1515/cmam-2023-0137","DOIUrl":"https://doi.org/10.1515/cmam-2023-0137","url":null,"abstract":"The problem of optimal recovering high-order mixed derivatives of bivariate functions with finite smoothness is studied. Based on the truncation method, an algorithm for numerical differentiation is constructed, which is order-optimal both in the sense of accuracy and in terms of the amount of involved Galerkin information. Numerical examples are provided to illustrate the fact that our approach can be implemented successfully.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2024-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139079091","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}
This work is concerned with the classical wave equation with a high-contrast coefficient in the spatial derivative operator. We first treat the periodic case, where we derive a new limit in the one-dimensional case. The behavior is illustrated numerically and contrasted to the higher-dimensional case. For general unstructured high-contrast coefficients, we present the Localized Orthogonal Decomposition and show a priori error estimates in suitably weighted norms. Numerical experiments illustrate the convergence rates in various settings.
{"title":"Wave Propagation in High-Contrast Media: Periodic and Beyond","authors":"Élise Fressart, Barbara Verfürth","doi":"10.1515/cmam-2023-0066","DOIUrl":"https://doi.org/10.1515/cmam-2023-0066","url":null,"abstract":"This work is concerned with the classical wave equation with a high-contrast coefficient in the spatial derivative operator. We first treat the periodic case, where we derive a new limit in the one-dimensional case. The behavior is illustrated numerically and contrasted to the higher-dimensional case. For general unstructured high-contrast coefficients, we present the Localized Orthogonal Decomposition and show a priori error estimates in suitably weighted norms. Numerical experiments illustrate the convergence rates in various settings.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2024-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139079450","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}
The JPEG algorithm is a defacto standard for image compression. We investigate whether adaptive mesh refinement can be used to optimize the compression ratio and propose a new adaptive image compression algorithm. We prove that it produces a quasi-optimal subdivision grid for a given error norm with high probability. This subdivision can be stored with very little overhead and thus leads to an efficient compression algorithm. We demonstrate experimentally, that the new algorithm can achieve better compression ratios than standard JPEG compression with no visible loss of quality on many images. The mathematical core of this work shows that Binev’s optimal tree approximation algorithm is applicable to image compression with high probability, when we assume small additive Gaussian noise on the pixels of the image.
{"title":"Adaptive Image Compression via Optimal Mesh Refinement","authors":"Michael Feischl, Hubert Hackl","doi":"10.1515/cmam-2023-0097","DOIUrl":"https://doi.org/10.1515/cmam-2023-0097","url":null,"abstract":"The JPEG algorithm is a defacto standard for image compression. We investigate whether adaptive mesh refinement can be used to optimize the compression ratio and propose a new adaptive image compression algorithm. We prove that it produces a quasi-optimal subdivision grid for a given error norm with high probability. This subdivision can be stored with very little overhead and thus leads to an efficient compression algorithm. We demonstrate experimentally, that the new algorithm can achieve better compression ratios than standard JPEG compression with no visible loss of quality on many images. The mathematical core of this work shows that Binev’s optimal tree approximation algorithm is applicable to image compression with high probability, when we assume small additive Gaussian noise on the pixels of the image.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2024-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139079452","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}
This work presents the multiharmonic analysis and derivation of functional type a posteriori estimates of a distributed eddy current optimal control problem and its state equation in a time-periodic setting. The existence and uniqueness of the solution of a weak space-time variational formulation for the optimality system and the forward problem are proved by deriving inf-sup and sup-sup conditions. Using the inf-sup and sup-sup conditions, we derive guaranteed, sharp and fully computable bounds of the approximation error for the optimal control problem and the forward problem in the functional type a posteriori estimation framework. We present here the first computational results on the derived estimates.
{"title":"A Posteriori Error Estimation for the Optimal Control of Time-Periodic Eddy Current Problems","authors":"Monika Wolfmayr","doi":"10.1515/cmam-2023-0119","DOIUrl":"https://doi.org/10.1515/cmam-2023-0119","url":null,"abstract":"This work presents the multiharmonic analysis and derivation of functional type a posteriori estimates of a distributed eddy current optimal control problem and its state equation in a time-periodic setting. The existence and uniqueness of the solution of a weak space-time variational formulation for the optimality system and the forward problem are proved by deriving inf-sup and sup-sup conditions. Using the inf-sup and sup-sup conditions, we derive guaranteed, sharp and fully computable bounds of the approximation error for the optimal control problem and the forward problem in the functional type a posteriori estimation framework. We present here the first computational results on the derived estimates.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2023-11-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138505899","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}
In this paper, we apply a two-grid scheme to the DG formulation of the 2D transient Navier–Stokes model. The two-grid algorithm consists of the following steps: Step 1 involves solving the nonlinear system on a coarse mesh with mesh size 𝐻, and Step 2 involves linearizing the nonlinear system by using the coarse grid solution on a fine mesh of mesh size ℎ and solving the resulting system to produce an approximate solution with desired accuracy. We establish optimal error estimates of the two-grid DG approximations for the velocity and pressure in energy and L2L^{2}-norms, respectively, for an appropriate choice of coarse and fine mesh parameters. We further discretize the two-grid DG model in time, using the backward Euler method, and derive the fully discrete error estimates. Finally, numerical results are presented to confirm the efficiency of the proposed scheme.
{"title":"Discontinuous Galerkin Two-Grid Method for the Transient Navier–Stokes Equations","authors":"Kallol Ray, Deepjyoti Goswami, Saumya Bajpai","doi":"10.1515/cmam-2023-0035","DOIUrl":"https://doi.org/10.1515/cmam-2023-0035","url":null,"abstract":"In this paper, we apply a two-grid scheme to the DG formulation of the 2D transient Navier–Stokes model. The two-grid algorithm consists of the following steps: Step 1 involves solving the nonlinear system on a coarse mesh with mesh size 𝐻, and Step 2 involves linearizing the nonlinear system by using the coarse grid solution on a fine mesh of mesh size ℎ and solving the resulting system to produce an approximate solution with desired accuracy. We establish optimal error estimates of the two-grid DG approximations for the velocity and pressure in energy and <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:msup> <m:mi>L</m:mi> <m:mn>2</m:mn> </m:msup> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2023-0035_ineq_0001.png\" /> <jats:tex-math>L^{2}</jats:tex-math> </jats:alternatives> </jats:inline-formula>-norms, respectively, for an appropriate choice of coarse and fine mesh parameters. We further discretize the two-grid DG model in time, using the backward Euler method, and derive the fully discrete error estimates. Finally, numerical results are presented to confirm the efficiency of the proposed scheme.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2023-11-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138505924","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}
We construct a C1C^{1}-P7P_{7} Bell finite element by restricting its normal derivative from a P6P_{6} polynomial to a P5P_{5} polynomial, and its second normal derivative from a P5P_{5} polynomial to a P4P_{4} polynomial, on the three edges of every triangle. On one triangle, the finite element space contains the P6P_{6} polynomial space. We show the method converges at order 7 in L2
{"title":"A 𝐶1-𝑃7 Bell Finite Element on Triangle","authors":"Xuejun Xu, Shangyou Zhang","doi":"10.1515/cmam-2023-0068","DOIUrl":"https://doi.org/10.1515/cmam-2023-0068","url":null,"abstract":"We construct a <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:msup> <m:mi>C</m:mi> <m:mn>1</m:mn> </m:msup> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2023-0068_ineq_0001.png\" /> <jats:tex-math>C^{1}</jats:tex-math> </jats:alternatives> </jats:inline-formula>-<jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:msub> <m:mi>P</m:mi> <m:mn>7</m:mn> </m:msub> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2023-0068_ineq_0002.png\" /> <jats:tex-math>P_{7}</jats:tex-math> </jats:alternatives> </jats:inline-formula> Bell finite element by restricting its normal derivative from a <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:msub> <m:mi>P</m:mi> <m:mn>6</m:mn> </m:msub> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2023-0068_ineq_0003.png\" /> <jats:tex-math>P_{6}</jats:tex-math> </jats:alternatives> </jats:inline-formula> polynomial to a <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:msub> <m:mi>P</m:mi> <m:mn>5</m:mn> </m:msub> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2023-0068_ineq_0004.png\" /> <jats:tex-math>P_{5}</jats:tex-math> </jats:alternatives> </jats:inline-formula> polynomial, and its second normal derivative from a <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:msub> <m:mi>P</m:mi> <m:mn>5</m:mn> </m:msub> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2023-0068_ineq_0004.png\" /> <jats:tex-math>P_{5}</jats:tex-math> </jats:alternatives> </jats:inline-formula> polynomial to a <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:msub> <m:mi>P</m:mi> <m:mn>4</m:mn> </m:msub> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2023-0068_ineq_0006.png\" /> <jats:tex-math>P_{4}</jats:tex-math> </jats:alternatives> </jats:inline-formula> polynomial, on the three edges of every triangle. On one triangle, the finite element space contains the <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:msub> <m:mi>P</m:mi> <m:mn>6</m:mn> </m:msub> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2023-0068_ineq_0003.png\" /> <jats:tex-math>P_{6}</jats:tex-math> </jats:alternatives> </jats:inline-formula> polynomial space. We show the method converges at order 7 in <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:msup> <m:mi>L</m:mi> <m:mn>2</m:mn> </m:msup> </m:math> <jats:inline-graphic xmlns:xlin","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2023-11-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138505895","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}
Reconstructing the pressure from given flow velocities is a task arising in various applications, and the standard approach uses the Navier–Stokes equations to derive a Poisson problem for the pressure p. That method, however, artificially increases the regularity requirements on both solution and data. In this context, we propose and analyze two alternative techniques to determine p∈L2(Ω){pin L^{2}(Omega)}. The first is an ultra-weak variational formulation applying integration by parts to shift all derivatives to the test functions. We present conforming finite element discretizations and prove optimal convergence of the resulting Galerkin–Petrov method. The second approach is a least-squares method for the original gradient equation, reformulated and solved as an artificial Stokes system. To simplify the incorporation of the given velocity within the right-hand side, we assume in the derivations that the velocity field is solenoidal. Yet this assumption is not restrictive, as we can use non-divergence-free approximations and even compressible velocities. Numerical experiments confirm the optimal a priori error estimates for both methods considered.
从给定的流速重建压力是一项在各种应用中出现的任务,标准方法使用Navier-Stokes方程来推导压力p的泊松问题。然而,这种方法人为地增加了对解和数据的规则性要求。在这种情况下,我们提出并分析了两种替代技术来确定p∈l2¹(Ω) {p In L^{2}(Omega)}。第一个是一个超弱变分公式,应用分部积分法将所有导数平移到测试函数。给出了合适的有限元离散化方法,并证明了Galerkin-Petrov方法的最佳收敛性。第二种方法是对原始梯度方程采用最小二乘法,将其重新表述并求解为人工Stokes系统。为了简化将给定的速度并入右边,我们在推导中假定速度场是螺线形的。然而,这个假设不是限制性的,因为我们可以使用非散度近似,甚至是可压缩速度。数值实验证实了两种方法的最优先验误差估计。
{"title":"Optimal Pressure Recovery Using an Ultra-Weak Finite Element Method for the Pressure Poisson Equation and a Least-Squares Approach for the Gradient Equation","authors":"Douglas R. Q. Pacheco, Olaf Steinbach","doi":"10.1515/cmam-2021-0242","DOIUrl":"https://doi.org/10.1515/cmam-2021-0242","url":null,"abstract":"Reconstructing the pressure from given flow velocities is a task arising in various applications, and the standard approach uses the Navier–Stokes equations to derive a Poisson problem for the pressure <jats:italic>p</jats:italic>. That method, however, artificially increases the regularity requirements on both solution and data. In this context, we propose and analyze two alternative techniques to determine <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:mrow> <m:mi>p</m:mi> <m:mo>∈</m:mo> <m:mrow> <m:msup> <m:mi>L</m:mi> <m:mn>2</m:mn> </m:msup> <m:mo></m:mo> <m:mrow> <m:mo stretchy=\"false\">(</m:mo> <m:mi mathvariant=\"normal\">Ω</m:mi> <m:mo stretchy=\"false\">)</m:mo> </m:mrow> </m:mrow> </m:mrow> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2021-0242_eq_0228.png\" /> <jats:tex-math>{pin L^{2}(Omega)}</jats:tex-math> </jats:alternatives> </jats:inline-formula>. The first is an ultra-weak variational formulation applying integration by parts to shift all derivatives to the test functions. We present conforming finite element discretizations and prove optimal convergence of the resulting Galerkin–Petrov method. The second approach is a least-squares method for the original gradient equation, reformulated and solved as an artificial Stokes system. To simplify the incorporation of the given velocity within the right-hand side, we assume in the derivations that the velocity field is solenoidal. Yet this assumption is not restrictive, as we can use non-divergence-free approximations and even compressible velocities. Numerical experiments confirm the optimal a priori error estimates for both methods considered.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2023-11-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138505856","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 an adaptive absorbing boundary layer technique for the nonlinear Schrödinger equation that is used in combination with the Time-splitting Fourier spectral method (TSSP) as the discretization for the NLS equations. We propose a new complex absorbing potential (CAP) function based on high order polynomials, with the major improvement that an explicit formula for the coefficients in the potential function is employed for adaptive parameter selection. This formula is obtained by an extension of the analysis in [R. Kosloff and D. Kosloff, Absorbing boundaries for wave propagation problems, J. Comput. Phys. 63 1986, 2, 363–376]. We also show that our imaginary potential function is more efficient than what is used in the literature. Numerical examples show that our ansatz is significantly better than existing approaches. We show that our approach can very accurately compute the solutions of the NLS equations in one dimension, including in the case of multi-dominant wave number solutions.
{"title":"Adaptive Absorbing Boundary Layer for the Nonlinear Schrödinger Equation","authors":"Hans Peter Stimming, Xin Wen, Norbert J. Mauser","doi":"10.1515/cmam-2023-0096","DOIUrl":"https://doi.org/10.1515/cmam-2023-0096","url":null,"abstract":"Abstract We present an adaptive absorbing boundary layer technique for the nonlinear Schrödinger equation that is used in combination with the Time-splitting Fourier spectral method (TSSP) as the discretization for the NLS equations. We propose a new complex absorbing potential (CAP) function based on high order polynomials, with the major improvement that an explicit formula for the coefficients in the potential function is employed for adaptive parameter selection. This formula is obtained by an extension of the analysis in [R. Kosloff and D. Kosloff, Absorbing boundaries for wave propagation problems, J. Comput. Phys. 63 1986, 2, 363–376]. We also show that our imaginary potential function is more efficient than what is used in the literature. Numerical examples show that our ansatz is significantly better than existing approaches. We show that our approach can very accurately compute the solutions of the NLS equations in one dimension, including in the case of multi-dominant wave number solutions.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135161159","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 This article develops and analyses a conforming virtual element scheme for the spatial discretization of parabolic integro-differential equations combined with backward Euler’s scheme for temporal discretization. With the help of Ritz–Voltera and L2 L^{2} projection operators, optimal a priori error estimates are established. Moreover, several numerical experiments are presented to confirm the computational efficiency of the proposed scheme and validate the theoretical findings.
{"title":"A Conforming Virtual Element Method for Parabolic Integro-Differential Equations","authors":"Sangita Yadav, Meghana Suthar, Sarvesh Kumar","doi":"10.1515/cmam-2023-0061","DOIUrl":"https://doi.org/10.1515/cmam-2023-0061","url":null,"abstract":"Abstract This article develops and analyses a conforming virtual element scheme for the spatial discretization of parabolic integro-differential equations combined with backward Euler’s scheme for temporal discretization. With the help of Ritz–Voltera and <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:msup> <m:mi>L</m:mi> <m:mn>2</m:mn> </m:msup> </m:math> L^{2} projection operators, optimal a priori error estimates are established. Moreover, several numerical experiments are presented to confirm the computational efficiency of the proposed scheme and validate the theoretical findings.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136059360","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}