首页 > 最新文献

Applied Numerical Mathematics最新文献

英文 中文
CSCMO: Relative permittivity-based complex shifted operator preconditioning method for solving time-harmonic Maxwell equations 基于相对介电常数的复移算子预处理方法求解时谐Maxwell方程
IF 2.4 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-09-02 DOI: 10.1016/j.apnum.2025.08.009
Zikang Qin , Xiaoyu Duan , Hengbin An , Haoxuan Zhang , Shaoliang Hu
In electromagnetic field modeling and simulation, a critical challenge lies in solving discretized time-harmonic Maxwell equations. The inherent complexity of these systems – including matrix indefiniteness and solution oscillations – poses significant difficulties for efficient numerical solutions. Preconditioned Krylov subspace methods have emerged as a standard approach for solving the large scale discretized time-harmonic Maxwell equations, where the construction of effective preconditioners remains pivotal. The shifted operator method represents a prominent preconditioning technique for Maxwell equations. Typically, a purely imaginary shift to the original differential operator is used, since this capitalizes on the observed phenomenon where solution oscillatory behavior diminishes with increasing modulus of the imaginary component of relative permittivity. In this paper, we propose a novel preconditioning technique by shifting both the real and imaginary parts of the relative permittivity. The motivation for proposing this kind of shifted operator preconditioning is based on theoretical analysis, which shows that decreasing the real part of the relative permittivity will improve the positive definiteness of the discretized matrix. The resulting preconditioned system admits efficient solution via multigrid methods. Some analysis shows that the spectral distribution of the preconditioned matrix is more clustered than the purely imaginary shift preconditioned matrix. Also, theoretical analysis for a special model shows that the condition number of the preconditioned system is less than its purely imaginary-shifted counterpart. Numerical results demonstrate that the proposed preconditioning is more effective than the other two state-of-the-art shift operator preconditioning methods.
在电磁场建模与仿真中,求解离散时谐麦克斯韦方程是一个关键问题。这些系统固有的复杂性——包括矩阵的不确定性和解的振荡——给有效的数值解带来了很大的困难。预条件Krylov子空间方法已成为求解大规模离散时谐Maxwell方程的标准方法,其中有效预条件的构造仍然是关键。移位算子法是求解麦克斯韦方程组的一种重要的预处理方法。通常,使用纯虚移到原始微分算子,因为这利用了观察到的现象,即溶液振荡行为随着相对介电常数虚分量的模量的增加而减少。在本文中,我们提出了一种新的预处理技术,通过移动相对介电常数的实部和虚部。提出这种移位算子预处理的动机是基于理论分析,理论分析表明减小相对介电常数的实部将提高离散矩阵的正确定性。所得到的预处理系统可以通过多网格方法进行有效的求解。分析表明,预条件矩阵的谱分布比纯虚移预条件矩阵更具有聚类性。对一个特殊模型的理论分析表明,预条件系统的条件数小于纯虚移系统的条件数。数值结果表明,该预处理方法比其他两种最先进的移位算子预处理方法更有效。
{"title":"CSCMO: Relative permittivity-based complex shifted operator preconditioning method for solving time-harmonic Maxwell equations","authors":"Zikang Qin ,&nbsp;Xiaoyu Duan ,&nbsp;Hengbin An ,&nbsp;Haoxuan Zhang ,&nbsp;Shaoliang Hu","doi":"10.1016/j.apnum.2025.08.009","DOIUrl":"10.1016/j.apnum.2025.08.009","url":null,"abstract":"<div><div>In electromagnetic field modeling and simulation, a critical challenge lies in solving discretized time-harmonic Maxwell equations. The inherent complexity of these systems – including matrix indefiniteness and solution oscillations – poses significant difficulties for efficient numerical solutions. Preconditioned Krylov subspace methods have emerged as a standard approach for solving the large scale discretized time-harmonic Maxwell equations, where the construction of effective preconditioners remains pivotal. The shifted operator method represents a prominent preconditioning technique for Maxwell equations. Typically, a purely imaginary shift to the original differential operator is used, since this capitalizes on the observed phenomenon where solution oscillatory behavior diminishes with increasing modulus of the imaginary component of relative permittivity. In this paper, we propose a novel preconditioning technique by shifting both the real and imaginary parts of the relative permittivity. The motivation for proposing this kind of shifted operator preconditioning is based on theoretical analysis, which shows that decreasing the real part of the relative permittivity will improve the positive definiteness of the discretized matrix. The resulting preconditioned system admits efficient solution via multigrid methods. Some analysis shows that the spectral distribution of the preconditioned matrix is more clustered than the purely imaginary shift preconditioned matrix. Also, theoretical analysis for a special model shows that the condition number of the preconditioned system is less than its purely imaginary-shifted counterpart. Numerical results demonstrate that the proposed preconditioning is more effective than the other two state-of-the-art shift operator preconditioning methods.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"219 ","pages":"Pages 145-169"},"PeriodicalIF":2.4,"publicationDate":"2025-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145106356","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
Robust globally divergence-free HDG finite element method for steady thermally coupled incompressible MHD flow 稳定热耦合不可压缩MHD流动的鲁棒全局无散度HDG有限元方法
IF 2.4 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-08-29 DOI: 10.1016/j.apnum.2025.08.007
Min Zhang , Zimo Zhu , Qijia Zhai , Xiaoping Xie
This paper develops a hybridizable discontinuous Galerkin (HDG) finite element method of arbitrary order for the steady thermally coupled incompressible Magnetohydrodynamics (MHD) flow. The HDG scheme uses piecewise polynomials of degrees k(k1),k,k1,k1 and k respectively for the approximations of the velocity, the magnetic field, the pressure, the magnetic pseudo-pressure, and the temperature in the interior of elements, and uses piecewise polynomials of degree k for their numerical traces on the interfaces of elements. The method is shown to yield globally divergence-free approximations of the velocity and magnetic fields. Existence and uniqueness results for the discrete scheme are given and O(hk)- optimal error estimates are derived for all the variables. Numerical experiments are provided to verify the obtained theoretical results.
本文针对稳态热耦合不可压缩磁流体动力学(MHD)流动问题,提出了一种任意阶的杂化不连续Galerkin (HDG)有限元方法。HDG格式分别采用k(k≥1)、k、k−1、k−1、k−1和k次分段多项式来逼近元件内部的速度、磁场、压力、磁伪压力和温度,并采用k次分段多项式来逼近元件界面上的数值迹线。该方法可以得到速度和磁场的全局无发散近似。给出了离散格式的存在唯一性结果,并得到了所有变量的0 (hk)-最优误差估计。数值实验验证了所得理论结果。
{"title":"Robust globally divergence-free HDG finite element method for steady thermally coupled incompressible MHD flow","authors":"Min Zhang ,&nbsp;Zimo Zhu ,&nbsp;Qijia Zhai ,&nbsp;Xiaoping Xie","doi":"10.1016/j.apnum.2025.08.007","DOIUrl":"10.1016/j.apnum.2025.08.007","url":null,"abstract":"<div><div>This paper develops a hybridizable discontinuous Galerkin (HDG) finite element method of arbitrary order for the steady thermally coupled incompressible Magnetohydrodynamics (MHD) flow. The HDG scheme uses piecewise polynomials of degrees <span><math><mrow><mi>k</mi><mo>(</mo><mi>k</mi><mo>≥</mo><mn>1</mn><mo>)</mo><mo>,</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>−</mo><mn>1</mn><mo>,</mo><mi>k</mi><mo>−</mo><mn>1</mn></mrow></math></span> and <span><math><mi>k</mi></math></span> respectively for the approximations of the velocity, the magnetic field, the pressure, the magnetic pseudo-pressure, and the temperature in the interior of elements, and uses piecewise polynomials of degree <span><math><mi>k</mi></math></span> for their numerical traces on the interfaces of elements. The method is shown to yield globally divergence-free approximations of the velocity and magnetic fields. Existence and uniqueness results for the discrete scheme are given and <span><math><mrow><mi>O</mi><mo>(</mo><msup><mi>h</mi><mi>k</mi></msup><mo>)</mo></mrow></math></span>- optimal error estimates are derived for all the variables. Numerical experiments are provided to verify the obtained theoretical results.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"219 ","pages":"Pages 19-40"},"PeriodicalIF":2.4,"publicationDate":"2025-08-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144989796","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
Compact difference method for Euler-Bernoulli beams and plates with nonlinear nonlocal damping 具有非线性非局部阻尼的Euler-Bernoulli梁和板的紧致差分法
IF 2.4 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-08-26 DOI: 10.1016/j.apnum.2025.08.008
Tao Guo , Yiqun Li , Wenlin Qiu
We investigate the numerical approximation to the Euler-Bernoulli (E-B) beams and plates with nonlinear nonlocal damping, which describes the damped mechanical behavior of beams and plates in real applications. We discretize the damping term by the composite Simpson’s rule and the six-point Simpson’s formula in the beam and plate problems, respectively, and then construct the fully discrete compact difference scheme for these problems. To account for the nonlinear-nonlocal term, we design several novel discrete norms to facilitate the error estimates of the damping term and the numerical scheme. The stability, convergence, and energy dissipation properties of the proposed scheme are proved, and numerical experiments are carried out to substantiate the theoretical findings.
研究了具有非线性非局部阻尼的欧拉-伯努利(E-B)梁和板的数值近似,该近似描述了梁和板在实际应用中的阻尼力学行为。分别采用复合辛普森规则和六点辛普森公式对梁和板问题中的阻尼项进行离散化,并构造了这些问题的完全离散紧致差分格式。为了考虑非线性非局部项,我们设计了几个新的离散范数来方便阻尼项的误差估计和数值格式。证明了该格式的稳定性、收敛性和能量耗散性,并进行了数值实验来验证理论结果。
{"title":"Compact difference method for Euler-Bernoulli beams and plates with nonlinear nonlocal damping","authors":"Tao Guo ,&nbsp;Yiqun Li ,&nbsp;Wenlin Qiu","doi":"10.1016/j.apnum.2025.08.008","DOIUrl":"10.1016/j.apnum.2025.08.008","url":null,"abstract":"<div><div>We investigate the numerical approximation to the Euler-Bernoulli (E-B) beams and plates with nonlinear nonlocal damping, which describes the damped mechanical behavior of beams and plates in real applications. We discretize the damping term by the composite Simpson’s rule and the six-point Simpson’s formula in the beam and plate problems, respectively, and then construct the fully discrete compact difference scheme for these problems. To account for the nonlinear-nonlocal term, we design several novel discrete norms to facilitate the error estimates of the damping term and the numerical scheme. The stability, convergence, and energy dissipation properties of the proposed scheme are proved, and numerical experiments are carried out to substantiate the theoretical findings.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"219 ","pages":"Pages 1-18"},"PeriodicalIF":2.4,"publicationDate":"2025-08-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144926147","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 biharmonic solver based on Fourier extension with oversampling technique for arbitrary domain 基于傅里叶扩展和过采样技术的任意域双调和求解器
IF 2.4 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-08-16 DOI: 10.1016/j.apnum.2025.08.005
Wenbin Li, Tinggang Zhao, Zhenyu Zhao
The biharmonic equation is commonly encountered in various fields such as elasticity theory, fluid dynamics, and image processing. Solving it on irregular domain presents a significant challenge. In this paper, Fourier extension method is used to solve the biharmonic equation on arbitrary domain. The method involves the oversampling collocation technique with the truncated singular value decomposition regularization, which comes out a spectral convergence rate for the smooth solution. This method only uses the function values on equidistant nodes and has the characteristics of less computation, strong universality and better accuracy. The effectiveness of the proposed method is demonstrated by a variety of numerical experiments.
双调和方程在弹性理论、流体动力学和图像处理等各个领域都经常遇到。在不规则域上求解这一问题是一个重大的挑战。本文采用傅里叶扩展法求解任意域上的双调和方程。该方法将过采样配置技术与截断奇异值分解正则化相结合,得到光滑解的谱收敛速率。该方法只使用等距节点上的函数值,具有计算量少、通用性强、精度高等特点。各种数值实验证明了该方法的有效性。
{"title":"A biharmonic solver based on Fourier extension with oversampling technique for arbitrary domain","authors":"Wenbin Li,&nbsp;Tinggang Zhao,&nbsp;Zhenyu Zhao","doi":"10.1016/j.apnum.2025.08.005","DOIUrl":"10.1016/j.apnum.2025.08.005","url":null,"abstract":"<div><div>The biharmonic equation is commonly encountered in various fields such as elasticity theory, fluid dynamics, and image processing. Solving it on irregular domain presents a significant challenge. In this paper, Fourier extension method is used to solve the biharmonic equation on arbitrary domain. The method involves the oversampling collocation technique with the truncated singular value decomposition regularization, which comes out a spectral convergence rate for the smooth solution. This method only uses the function values on equidistant nodes and has the characteristics of less computation, strong universality and better accuracy. The effectiveness of the proposed method is demonstrated by a variety of numerical experiments.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"218 ","pages":"Pages 261-274"},"PeriodicalIF":2.4,"publicationDate":"2025-08-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144888764","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
Fractional convection-diffusion systems in complex 2D and 3D geometries: A Bernoulli polynomial-based kernel method 复杂二维和三维几何中的分数对流扩散系统:基于伯努利多项式的核方法
IF 2.4 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-08-15 DOI: 10.1016/j.apnum.2025.08.004
Mojtaba Fardi , Mahmoud A. Zaky , Babak Azarnavid
This study presents an accurate meshless method for the efficient solution of nonlinear time-fractional convection-diffusion systems in complex two- and three-dimensional geometries. The proposed approach combines spatial discretization using a Bernoulli polynomial kernel function with temporal discretization via the backward differentiation formula. By employing positive definite kernels, the method achieves high spatial accuracy, while the use of the backward differentiation formula ensures high-order temporal accuracy. Convergence conditions and error bounds are rigorously analyzed using the Mittag-Leffler function. Error estimates are derived based on the spectral properties of the associated matrices, and inequalities describing error propagation over time are established. The method is tested on a variety of benchmark problems, including the Brusselator model and nonlinear coupled convection-diffusion systems, across both 2D and 3D domains. Extensive numerical experiments are carried out on various geometries-such as rectangular, circular, and spherical shapes-demonstrating the method’s robustness and accuracy in handling both regular and irregular computational domains.
本文提出了一种精确的无网格方法,用于求解复杂二维和三维几何结构的非线性时分式对流扩散系统。该方法结合了使用伯努利多项式核函数的空间离散化和通过后向微分公式的时间离散化。该方法利用正定核实现了较高的空间精度,同时利用后向微分公式保证了高阶时间精度。利用Mittag-Leffler函数严格分析了收敛条件和误差界。误差估计是根据相关矩阵的谱特性推导出来的,并且建立了描述误差随时间传播的不等式。该方法在各种基准问题上进行了测试,包括布鲁塞尔模型和非线性耦合对流扩散系统,跨越二维和三维领域。广泛的数值实验进行了各种几何形状-如矩形,圆形和球形-证明了该方法的鲁棒性和准确性在处理规则和不规则计算域。
{"title":"Fractional convection-diffusion systems in complex 2D and 3D geometries: A Bernoulli polynomial-based kernel method","authors":"Mojtaba Fardi ,&nbsp;Mahmoud A. Zaky ,&nbsp;Babak Azarnavid","doi":"10.1016/j.apnum.2025.08.004","DOIUrl":"10.1016/j.apnum.2025.08.004","url":null,"abstract":"<div><div>This study presents an accurate meshless method for the efficient solution of nonlinear time-fractional convection-diffusion systems in complex two- and three-dimensional geometries. The proposed approach combines spatial discretization using a Bernoulli polynomial kernel function with temporal discretization via the backward differentiation formula. By employing positive definite kernels, the method achieves high spatial accuracy, while the use of the backward differentiation formula ensures high-order temporal accuracy. Convergence conditions and error bounds are rigorously analyzed using the Mittag-Leffler function. Error estimates are derived based on the spectral properties of the associated matrices, and inequalities describing error propagation over time are established. The method is tested on a variety of benchmark problems, including the Brusselator model and nonlinear coupled convection-diffusion systems, across both 2D and 3D domains. Extensive numerical experiments are carried out on various geometries-such as rectangular, circular, and spherical shapes-demonstrating the method’s robustness and accuracy in handling both regular and irregular computational domains.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"218 ","pages":"Pages 275-297"},"PeriodicalIF":2.4,"publicationDate":"2025-08-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144896067","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
Adaptive SIPG method for approximations of parabolic boundary control problems with bilateral box constraints on Neumann boundary Neumann边界上双侧框约束抛物型边界控制问题逼近的自适应SIPG方法
IF 2.4 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-08-13 DOI: 10.1016/j.apnum.2025.08.002
Ram Manohar , B․ V․ Rathish Kumar , Kedarnath Buda , Rajen Kumar Sinha
This study presents an a posteriori error analysis of adaptive finite element approximations of parabolic boundary control problems with bilateral box constraints that act on a Neumann boundary. The control problem is discretized using the symmetric interior penalty Galerkin (SIPG) technique. We derive both reliable and efficient type residual-based error estimators coupling with the data oscillations. The implementation of these error estimators serves as a guide for the adaptive mesh refinement process, indicating whether or not more refinement is required. Although the control error estimator accurately captured control approximation errors, it had limitations in terms of guiding refinement localization in critical circumstances. To overcome this, an alternative control indicator was used in numerical tests. The results demonstrated the clear superiority of adaptive refinements over uniform refinements, confirming the proposed approach’s effectiveness in achieving accurate solutions while optimizing computational efficiency. Numerical experiments showcase the effectiveness of the derived error estimators.
本研究提出了具有双边框约束作用于诺伊曼边界的抛物边界控制问题的自适应有限元近似的后检验误差分析。采用对称内罚伽辽金(SIPG)技术对控制问题进行离散化。我们得到了与数据振荡耦合的可靠和有效的基于残差的误差估计器。这些误差估计的实现可以作为自适应网格细化过程的指南,表明是否需要更多的细化。虽然控制误差估计器能准确捕获控制逼近误差,但在指导关键情况下的精化定位方面存在局限性。为了克服这一点,在数值试验中使用了一种替代控制指示器。结果表明,自适应细化明显优于均匀细化,证实了所提方法在获得精确解的同时优化计算效率的有效性。数值实验证明了该误差估计方法的有效性。
{"title":"Adaptive SIPG method for approximations of parabolic boundary control problems with bilateral box constraints on Neumann boundary","authors":"Ram Manohar ,&nbsp;B․ V․ Rathish Kumar ,&nbsp;Kedarnath Buda ,&nbsp;Rajen Kumar Sinha","doi":"10.1016/j.apnum.2025.08.002","DOIUrl":"10.1016/j.apnum.2025.08.002","url":null,"abstract":"<div><div>This study presents an a posteriori error analysis of adaptive finite element approximations of parabolic boundary control problems with bilateral box constraints that act on a Neumann boundary. The control problem is discretized using the symmetric interior penalty Galerkin (SIPG) technique. We derive both reliable and efficient type residual-based error estimators coupling with the data oscillations. The implementation of these error estimators serves as a guide for the adaptive mesh refinement process, indicating whether or not more refinement is required. Although the control error estimator accurately captured control approximation errors, it had limitations in terms of guiding refinement localization in critical circumstances. To overcome this, an alternative control indicator was used in numerical tests. The results demonstrated the clear superiority of adaptive refinements over uniform refinements, confirming the proposed approach’s effectiveness in achieving accurate solutions while optimizing computational efficiency. Numerical experiments showcase the effectiveness of the derived error estimators.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"219 ","pages":"Pages 170-201"},"PeriodicalIF":2.4,"publicationDate":"2025-08-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145106357","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 estimates for the finite element approximation of the convection–diffusion–reaction equation based on the variational multiscale concept 基于变分多尺度概念的对流扩散反应方程有限元近似的后验误差估计
IF 2.4 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-08-11 DOI: 10.1016/j.apnum.2025.08.003
Ramon Codina , Hauke Gravenkamp , Sheraz Ahmed Khan
In this study, we employ the variational multiscale (VMS) concept to develop a posteriori error estimates for the stationary convection-diffusion-reaction equation. The variational multiscale method is based on splitting the continuous part of the problem into a resolved scale (coarse scale) and an unresolved scale (fine scale). The unresolved scale (also known as the sub-grid scale) is modeled by choosing it proportional to the component of the residual orthogonal to the finite element space, leading to the orthogonal sub-grid scale (OSGS) method. The idea is then to use the modeled sub-grid scale as an error estimator, considering its contribution in the element interiors and on the edges. We present the results of the a priori analysis and two different strategies for the a posteriori error analysis for the OSGS method. Our proposal is to use a scaled norm of the sub-grid scales as an a posteriori error estimate in the so-called stabilized norm of the problem. This norm has control over the convective term, which is necessary for convection-dominated problems. Numerical examples show the reliable performance of the proposed error estimator compared to other error estimators belonging to the variational multiscale family.
在这项研究中,我们采用变分多尺度(VMS)的概念来建立稳态对流-扩散-反应方程的后验误差估计。变分多尺度方法是将问题的连续部分分解为一个已解尺度(粗尺度)和一个未解尺度(细尺度)。通过选择与有限元空间正交的残差分量成比例的未解析尺度(也称为子网格尺度)来建模,从而产生正交子网格尺度(OSGS)方法。然后,我们的想法是使用建模的子网格尺度作为误差估计器,考虑到它在元素内部和边缘的贡献。我们给出了先验分析的结果和两种不同的策略,用于OSGS方法的后验误差分析。我们的建议是使用子网格尺度的缩放范数作为问题的所谓稳定范数的后验误差估计。该范数可以控制对流项,这对于对流主导的问题是必要的。数值算例表明,与其他的变分多尺度误差估计器相比,所提出的误差估计器具有可靠的性能。
{"title":"A posteriori error estimates for the finite element approximation of the convection–diffusion–reaction equation based on the variational multiscale concept","authors":"Ramon Codina ,&nbsp;Hauke Gravenkamp ,&nbsp;Sheraz Ahmed Khan","doi":"10.1016/j.apnum.2025.08.003","DOIUrl":"10.1016/j.apnum.2025.08.003","url":null,"abstract":"<div><div>In this study, we employ the variational multiscale (VMS) concept to develop a posteriori error estimates for the stationary convection-diffusion-reaction equation. The variational multiscale method is based on splitting the continuous part of the problem into a resolved scale (coarse scale) and an unresolved scale (fine scale). The unresolved scale (also known as the sub-grid scale) is modeled by choosing it proportional to the component of the residual orthogonal to the finite element space, leading to the orthogonal sub-grid scale (OSGS) method. The idea is then to use the modeled sub-grid scale as an error estimator, considering its contribution in the element interiors and on the edges. We present the results of the a priori analysis and two different strategies for the a posteriori error analysis for the OSGS method. Our proposal is to use a scaled norm of the sub-grid scales as an a posteriori error estimate in the so-called stabilized norm of the problem. This norm has control over the convective term, which is necessary for convection-dominated problems. Numerical examples show the reliable performance of the proposed error estimator compared to other error estimators belonging to the variational multiscale family.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"218 ","pages":"Pages 238-260"},"PeriodicalIF":2.4,"publicationDate":"2025-08-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144886791","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
2D and 3D reconstructions in acousto-electric tomography via two-point gradient Kaczmarz-type algorithm 基于两点梯度kaczmarz型算法的声电层析成像二维和三维重建
IF 2.4 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-08-11 DOI: 10.1016/j.apnum.2025.07.015
Kai Zhu, Min Zhong
This paper presents a Kaczmarz type two-point gradient algorithm with the general convex penalty functional Θ (KTPG-Θ), for efficient reconstruction of conductivity in acousto-electric tomography (AET). The algorithm optimizes a convex functional with flexible non-smooth regularization terms, such as L1-like and total variation-like, to handle sparse and discontinuous conductivity distributions. By cyclically processing the measurement equations and incorporating an acceleration strategy, the proposed method achieves high computational efficiency while ensuring convergence. Numerical experiments on both synthetic and realistic phantoms demonstrate the method’s superior accuracy, strong noise robustness, and ability to resolve fine details. Beyond AET, the KTPG-Θ framework can be applied to a wide range of nonlinear inverse problems involving systems of equations, showcasing its potential for broader applications in science and engineering.
本文提出了一种具有一般凸罚函数Θ (KTPG-Θ)的Kaczmarz型两点梯度算法,用于声电断层扫描(AET)中电导率的有效重建。该算法优化了一个具有柔性非光滑正则化项的凸泛函,如类l1和类总变差,以处理稀疏和不连续的电导率分布。该方法通过对测量方程进行循环处理,并结合加速策略,在保证收敛性的同时获得了较高的计算效率。仿真实验结果表明,该方法具有较好的精度、较强的噪声鲁棒性和较好的细节处理能力。除了AET之外,KTPG-Θ框架还可以应用于涉及方程系统的广泛非线性逆问题,展示了其在科学和工程领域更广泛应用的潜力。
{"title":"2D and 3D reconstructions in acousto-electric tomography via two-point gradient Kaczmarz-type algorithm","authors":"Kai Zhu,&nbsp;Min Zhong","doi":"10.1016/j.apnum.2025.07.015","DOIUrl":"10.1016/j.apnum.2025.07.015","url":null,"abstract":"<div><div>This paper presents a Kaczmarz type two-point gradient algorithm with the general convex penalty functional <span><math><mstyle><mi>Θ</mi></mstyle></math></span> (KTPG-<span><math><mstyle><mi>Θ</mi></mstyle></math></span>), for efficient reconstruction of conductivity in acousto-electric tomography (AET). The algorithm optimizes a convex functional with flexible non-smooth regularization terms, such as <span><math><msup><mi>L</mi><mn>1</mn></msup></math></span>-like and total variation-like, to handle sparse and discontinuous conductivity distributions. By cyclically processing the measurement equations and incorporating an acceleration strategy, the proposed method achieves high computational efficiency while ensuring convergence. Numerical experiments on both synthetic and realistic phantoms demonstrate the method’s superior accuracy, strong noise robustness, and ability to resolve fine details. Beyond AET, the KTPG-<span><math><mstyle><mi>Θ</mi></mstyle></math></span> framework can be applied to a wide range of nonlinear inverse problems involving systems of equations, showcasing its potential for broader applications in science and engineering.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"218 ","pages":"Pages 220-237"},"PeriodicalIF":2.4,"publicationDate":"2025-08-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144886790","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 boundary-corrected weak Galerkin mixed finite method for elliptic interface problems with curved interfaces 具有弯曲界面的椭圆界面问题的边界修正弱Galerkin混合有限方法
IF 2.4 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-08-06 DOI: 10.1016/j.apnum.2025.08.001
Yongli Hou , Yi Liu , Yanqiu Wang
We propose a boundary-corrected weak Galerkin mixed finite element method for solving elliptic interface problems in 2D domains with curved interfaces. The method is formulated on body-fitted polygonal meshes, where interface edges are straight and may not align exactly with the curved physical interface. To address this discrepancy, a boundary value correction technique is employed to transfer the interface conditions from the physical interface to the approximate interface using a Taylor expansion approach. The Neumann interface condition is then weakly imposed in the variational formulation. This approach eliminates the need for numerical integration on curved elements, thereby reducing implementation complexity. We establish optimal-order convergence in the energy norm for arbitrary-order discretizations. Numerical results are provided to support the theoretical findings.
提出了一种边界修正的弱Galerkin混合有限元方法,用于求解具有曲面界面的二维区域中的椭圆界面问题。该方法是在贴体多边形网格上制定的,其中界面边缘是直的,可能与弯曲的物理界面不完全对齐。为了解决这种差异,采用边界值校正技术,使用泰勒展开方法将界面条件从物理界面转移到近似界面。然后在变分公式中弱地施加诺伊曼界面条件。这种方法消除了对曲面元素进行数值积分的需要,从而降低了实现的复杂性。建立了任意阶离散化的能量范数的最优阶收敛性。数值结果支持了理论结果。
{"title":"A boundary-corrected weak Galerkin mixed finite method for elliptic interface problems with curved interfaces","authors":"Yongli Hou ,&nbsp;Yi Liu ,&nbsp;Yanqiu Wang","doi":"10.1016/j.apnum.2025.08.001","DOIUrl":"10.1016/j.apnum.2025.08.001","url":null,"abstract":"<div><div>We propose a boundary-corrected weak Galerkin mixed finite element method for solving elliptic interface problems in 2D domains with curved interfaces. The method is formulated on body-fitted polygonal meshes, where interface edges are straight and may not align exactly with the curved physical interface. To address this discrepancy, a boundary value correction technique is employed to transfer the interface conditions from the physical interface to the approximate interface using a Taylor expansion approach. The Neumann interface condition is then weakly imposed in the variational formulation. This approach eliminates the need for numerical integration on curved elements, thereby reducing implementation complexity. We establish optimal-order convergence in the energy norm for arbitrary-order discretizations. Numerical results are provided to support the theoretical findings.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"218 ","pages":"Pages 201-219"},"PeriodicalIF":2.4,"publicationDate":"2025-08-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144842654","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
Semi-implicit fully exactly well-balanced schemes for the two-layer shallow water system 双层浅水系统的半隐式完全精确平衡方案
IF 2.4 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-07-29 DOI: 10.1016/j.apnum.2025.07.014
C. Caballero-Cárdenas , M.J. Castro , C. Chalons , T. Morales de Luna , M.L. Muñoz-Ruiz
This work addresses the design of semi-implicit numerical schemes that are fully exactly well-balanced for the two-layer shallow water system, meaning that they are capable of preserving every possible steady state, and not only the lake-at-rest ones. The proposed approach exhibits better performance compared to standard explicit methods in low-Froude number regimes, where wave propagation speeds significantly exceed flow velocities, thereby reducing the computational cost associated with long-time simulations. The methodology relies on a combination of splitting strategies and relaxation techniques to construct first- and second-order semi-implicit schemes that satisfy the fully exactly well-balanced property.
这项工作解决了半隐式数值方案的设计,这些方案对于两层浅水系统来说是完全平衡的,这意味着它们能够保持每一种可能的稳态,而不仅仅是静止的湖泊。与标准显式方法相比,该方法在低傅鲁德数条件下表现出更好的性能,在低傅鲁德数条件下,波的传播速度明显超过流的速度,从而减少了与长时间模拟相关的计算成本。该方法采用分裂策略和松弛技术相结合的方法来构造满足完全正平衡性质的一阶和二阶半隐式格式。
{"title":"Semi-implicit fully exactly well-balanced schemes for the two-layer shallow water system","authors":"C. Caballero-Cárdenas ,&nbsp;M.J. Castro ,&nbsp;C. Chalons ,&nbsp;T. Morales de Luna ,&nbsp;M.L. Muñoz-Ruiz","doi":"10.1016/j.apnum.2025.07.014","DOIUrl":"10.1016/j.apnum.2025.07.014","url":null,"abstract":"<div><div>This work addresses the design of semi-implicit numerical schemes that are fully exactly well-balanced for the two-layer shallow water system, meaning that they are capable of preserving every possible steady state, and not only the lake-at-rest ones. The proposed approach exhibits better performance compared to standard explicit methods in low-Froude number regimes, where wave propagation speeds significantly exceed flow velocities, thereby reducing the computational cost associated with long-time simulations. The methodology relies on a combination of splitting strategies and relaxation techniques to construct first- and second-order semi-implicit schemes that satisfy the fully exactly well-balanced property.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"218 ","pages":"Pages 128-147"},"PeriodicalIF":2.4,"publicationDate":"2025-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144720945","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
期刊
Applied Numerical 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学术文献互助群
群 号:604180095
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1