首页 > 最新文献

Journal of Computational and Applied Mathematics最新文献

英文 中文
Third order two-step Runge–Kutta–Chebyshev methods 三阶两步 Runge-Kutta-Chebyshev 方法
IF 2.1 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2024-09-26 DOI: 10.1016/j.cam.2024.116291
The well-known high order stabilized codes (such as DUMKA and ROCK) have several drawbacks: numerically obtained stability polynomials (which do not have a closed analytic form), poor internal stability and convergence. RKC-type methods have much better computational properties. However, these types of methods currently have a second order maximum. In this paper, a family of third order stabilized methods with an explicit analytical solution of stability polynomials is presented. This was made possible by usage of two-step Runge–Kutta methods. A new code TSRKC3 is proposed, illustrated by several examples, and compared to existing programs.
众所周知的高阶稳定代码(如 DUMKA 和 ROCK)有几个缺点:数值获得的稳定多项式(没有封闭的解析形式)、内部稳定性和收敛性差。RKC 类方法的计算性能要好得多。然而,这类方法目前只有二阶最大值。本文提出了一系列三阶稳定方法,这些方法具有明确的稳定多项式解析解。这是通过使用两步 Runge-Kutta 方法实现的。本文提出了一种新代码 TSRKC3,并通过几个例子进行了说明,还与现有程序进行了比较。
{"title":"Third order two-step Runge–Kutta–Chebyshev methods","authors":"","doi":"10.1016/j.cam.2024.116291","DOIUrl":"10.1016/j.cam.2024.116291","url":null,"abstract":"<div><div>The well-known high order stabilized codes (such as DUMKA and ROCK) have several drawbacks: numerically obtained stability polynomials (which do not have a closed analytic form), poor internal stability and convergence. RKC-type methods have much better computational properties. However, these types of methods currently have a second order maximum. In this paper, a family of third order stabilized methods with an explicit analytical solution of stability polynomials is presented. This was made possible by usage of two-step Runge–Kutta methods. A new code TSRKC3 is proposed, illustrated by several examples, and compared to existing programs.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-09-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142326350","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
Finite difference methods for stochastic Helmholtz equation driven by white noise 白噪声驱动的随机亥姆霍兹方程的有限差分法
IF 2.1 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2024-09-20 DOI: 10.1016/j.cam.2024.116286
In this paper, we propose two numerical methods for the stochastic Helmholtz equation driven by white noise. We obtain the approximate stochastic problem by approximating the white noise with piecewise constant process, provide some regularity of its solution and the truncation error between the approximate stochastic problem and the original problem. The limitation on the wave number k of the finite difference method (FDM) is analyzed and a stochastic finite difference (SFD) scheme is presented. The error analysis shows that the stochastic finite difference method is efficient with a certain convergence rate. Numerical experiments are provided to examine our theoretical results.
本文针对白噪声驱动的随机亥姆霍兹方程提出了两种数值方法。我们通过用片断常数过程逼近白噪声得到近似随机问题,给出了其解的一些规律性以及近似随机问题与原问题之间的截断误差。分析了有限差分法(FDM)对波数 k 的限制,并提出了一种随机有限差分法(SFD)方案。误差分析表明,随机有限差分法具有一定的收敛效率。我们还提供了数值实验来检验我们的理论结果。
{"title":"Finite difference methods for stochastic Helmholtz equation driven by white noise","authors":"","doi":"10.1016/j.cam.2024.116286","DOIUrl":"10.1016/j.cam.2024.116286","url":null,"abstract":"<div><div>In this paper, we propose two numerical methods for the stochastic Helmholtz equation driven by white noise. We obtain the approximate stochastic problem by approximating the white noise with piecewise constant process, provide some regularity of its solution and the truncation error between the approximate stochastic problem and the original problem. The limitation on the wave number <span><math><mi>k</mi></math></span> of the finite difference method (FDM) is analyzed and a stochastic finite difference (SFD) scheme is presented. The error analysis shows that the stochastic finite difference method is efficient with a certain convergence rate. Numerical experiments are provided to examine our theoretical results.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-09-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142320199","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
Poisson noise removal based on non-convex hybrid regularizers 基于非凸混合正则的泊松噪声消除
IF 2.1 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2024-09-19 DOI: 10.1016/j.cam.2024.116289
The presence of TV regularizer always induces an unsatisfactory staircase effect. To overcome the staircase while better sustaining edge information, this work proposes a novel model for Poisson noise removal. The model is based on non-convex mixed regularizers, which involves introducing a non-convex penalty into a composition of the total variation and the higher-order total variation. The iterative reweighted l1 algorithm was used to convert the non-convex model into a convex one. The classic alternating direction method of multipliers was then employed to obtain approximate solutions of the model. When applying this model to degraded images contaminated by Poisson noise of medium to high intensity, its performance in noise suppression was tested. The regularizer was compared with others in the terms of the visual effect of the picture, time cost, and several commonly accepted quantitative indicators for evaluation, such as peak signal-to-noise ratio, feature similarity index and structural similarity index. Numerical experiments showed that the present model not only eliminates block artifacts but also retains sharp edges.
电视正则的存在总是会引起令人不满意的阶梯效应。为了克服阶梯效应,同时更好地保留边缘信息,本研究提出了一种新的泊松噪声去除模型。该模型基于非凸混合正则,即在总变化和高阶总变化的组合中引入非凸惩罚。迭代加权 l1 算法用于将非凸模型转换为凸模型。然后采用经典的交替方向乘法来获得模型的近似解。在将该模型应用于受中高强度泊松噪声污染的退化图像时,对其噪声抑制性能进行了测试。该正则化器与其他正则化器在图像的视觉效果、时间成本以及几个公认的量化评估指标(如峰值信噪比、特征相似性指数和结构相似性指数)方面进行了比较。数值实验表明,本模型不仅能消除块状伪影,还能保留锐利的边缘。
{"title":"Poisson noise removal based on non-convex hybrid regularizers","authors":"","doi":"10.1016/j.cam.2024.116289","DOIUrl":"10.1016/j.cam.2024.116289","url":null,"abstract":"<div><div>The presence of TV regularizer always induces an unsatisfactory staircase effect. To overcome the staircase while better sustaining edge information, this work proposes a novel model for Poisson noise removal. The model is based on non-convex mixed regularizers, which involves introducing a non-convex penalty into a composition of the total variation and the higher-order total variation. The iterative reweighted <span><math><msub><mrow><mi>l</mi></mrow><mrow><mn>1</mn></mrow></msub></math></span> algorithm was used to convert the non-convex model into a convex one. The classic alternating direction method of multipliers was then employed to obtain approximate solutions of the model. When applying this model to degraded images contaminated by Poisson noise of medium to high intensity, its performance in noise suppression was tested. The regularizer was compared with others in the terms of the visual effect of the picture, time cost, and several commonly accepted quantitative indicators for evaluation, such as peak signal-to-noise ratio, feature similarity index and structural similarity index. Numerical experiments showed that the present model not only eliminates block artifacts but also retains sharp edges.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142326348","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 H∞ control for LFC of discrete T–S fuzzy MAPS with DFIG and time-varying delays 带双馈变流器和时变延迟的离散 T-S 模糊 MAPS LFC 的鲁棒 H∞ 控制
IF 2.1 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2024-09-18 DOI: 10.1016/j.cam.2024.116271
The H LFC control problem for a class of nonlinear power systems with time-varying delays is under study. Considering uncertainties arising from nonlinear issues such as the generation rate constraint (GRC) and governor dead band (GDB), as well as the high variability of renewable energy sources like wind power, the model is transformed into a discrete-time Takagi–Sugeno (T–S) fuzzy model with parameter uncertainty. By constructing a Lyapunov–Krasovskii functional, and employing difference inequalities and generalized cross-convex matrix inequalities, sufficient conditions for the asymptotic stability of power systems are provided. Based on the obtained conditions, a controller is designed to ensure the asymptotic stability of Multi-Area Power Systems (MAPS), with the performance index being H. Finally, simulation results demonstrate the correctness and effectiveness of the theorem.
我们正在研究一类具有时变延迟的非线性电力系统的 H∞ LFC 控制问题。考虑到发电率约束(GRC)和调速器死区(GDB)等非线性问题所产生的不确定性,以及风电等可再生能源的高变化性,该模型被转化为具有参数不确定性的离散时间高木-菅野(T-S)模糊模型。通过构建 Lyapunov-Krasovskii 函数,并利用差分不等式和广义交叉凸矩阵不等式,提供了电力系统渐近稳定性的充分条件。根据所获得的条件,设计了一种控制器来确保多区域电力系统(MAPS)的渐近稳定性,其性能指标为 H∞。最后,仿真结果证明了定理的正确性和有效性。
{"title":"Robust H∞ control for LFC of discrete T–S fuzzy MAPS with DFIG and time-varying delays","authors":"","doi":"10.1016/j.cam.2024.116271","DOIUrl":"10.1016/j.cam.2024.116271","url":null,"abstract":"<div><div>The <span><math><msub><mrow><mi>H</mi></mrow><mrow><mi>∞</mi></mrow></msub></math></span> LFC control problem for a class of nonlinear power systems with time-varying delays is under study. Considering uncertainties arising from nonlinear issues such as the generation rate constraint (GRC) and governor dead band (GDB), as well as the high variability of renewable energy sources like wind power, the model is transformed into a discrete-time Takagi–Sugeno (T–S) fuzzy model with parameter uncertainty. By constructing a Lyapunov–Krasovskii functional, and employing difference inequalities and generalized cross-convex matrix inequalities, sufficient conditions for the asymptotic stability of power systems are provided. Based on the obtained conditions, a controller is designed to ensure the asymptotic stability of Multi-Area Power Systems (MAPS), with the performance index being <span><math><msub><mrow><mi>H</mi></mrow><mrow><mi>∞</mi></mrow></msub></math></span>. Finally, simulation results demonstrate the correctness and effectiveness of the theorem.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142312507","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
Fading regularization method for an inverse boundary value problem associated with the biharmonic equation 与双谐方程相关的反边界值问题的消隐正则化方法
IF 2.1 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2024-09-18 DOI: 10.1016/j.cam.2024.116285
In this paper, we propose a numerical algorithm that combines the fading regularization method with the method of fundamental solutions (MFS) to solve a Cauchy problem associated with the biharmonic equation. We introduce a new stopping criterion for the iterative process and compare its performance with previous criteria. Numerical simulations using MFS validate the accuracy of this stopping criterion for both compatible and noisy data and demonstrate the convergence, stability, and efficiency of the proposed algorithm, as well as its ability to deblur noisy data.
在本文中,我们提出了一种结合消隐正则化方法和基本解法(MFS)的数值算法,用于解决与双谐波方程相关的柯西问题。我们为迭代过程引入了一个新的停止准则,并将其性能与之前的准则进行了比较。使用 MFS 进行的数值模拟验证了这一停止准则对于兼容数据和噪声数据的准确性,并证明了所提算法的收敛性、稳定性和效率,以及其消除噪声数据的能力。
{"title":"Fading regularization method for an inverse boundary value problem associated with the biharmonic equation","authors":"","doi":"10.1016/j.cam.2024.116285","DOIUrl":"10.1016/j.cam.2024.116285","url":null,"abstract":"<div><div>In this paper, we propose a numerical algorithm that combines the fading regularization method with the method of fundamental solutions (MFS) to solve a Cauchy problem associated with the biharmonic equation. We introduce a new stopping criterion for the iterative process and compare its performance with previous criteria. Numerical simulations using MFS validate the accuracy of this stopping criterion for both compatible and noisy data and demonstrate the convergence, stability, and efficiency of the proposed algorithm, as well as its ability to deblur noisy data.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142322079","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Error analysis of the explicit-invariant energy quadratization (EIEQ) numerical scheme for solving the Allen–Cahn equation 用于求解艾伦-卡恩方程的显式不变能量四分法(EIEQ)数值方案的误差分析
IF 2.1 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2024-09-16 DOI: 10.1016/j.cam.2024.116224
This paper focuses on the error analysis of a first-order, time-discrete scheme for solving the nonlinear Allen–Cahn equation. The discretization of the nonlinear potential is achieved through the EIEQ method, which employs an auxiliary variable to linearize the nonlinear double-well potential effectively. The energy stability of the scheme is demonstrated, along with its decoupled type implementation. Under a set of reasonable assumptions related to boundedness and continuity, an extensive error analysis is performed. This analysis results in the establishment of L2 and H1 error bounds for the numerical solution. Furthermore, a variety of numerical examples are conducted to illustrate the accuracy of the EIEQ scheme, highlighting its effectiveness in addressing complex dynamical systems governed by the Allen–Cahn equation.
本文重点分析了求解非线性 Allen-Cahn 方程的一阶时间离散方案的误差。非线性势的离散化是通过 EIEQ 方法实现的,该方法采用了一个辅助变量来有效地线性化非线性双阱势。演示了该方案的能量稳定性及其解耦类型的实现。在一系列与有界性和连续性相关的合理假设下,进行了广泛的误差分析。通过分析,建立了数值解的 L2 和 H1 误差边界。此外,还通过各种数值示例说明了 EIEQ 方案的准确性,突出了它在处理受 Allen-Cahn 方程控制的复杂动力系统时的有效性。
{"title":"Error analysis of the explicit-invariant energy quadratization (EIEQ) numerical scheme for solving the Allen–Cahn equation","authors":"","doi":"10.1016/j.cam.2024.116224","DOIUrl":"10.1016/j.cam.2024.116224","url":null,"abstract":"<div><div>This paper focuses on the error analysis of a first-order, time-discrete scheme for solving the nonlinear Allen–Cahn equation. The discretization of the nonlinear potential is achieved through the EIEQ method, which employs an auxiliary variable to linearize the nonlinear double-well potential effectively. The energy stability of the scheme is demonstrated, along with its decoupled type implementation. Under a set of reasonable assumptions related to boundedness and continuity, an extensive error analysis is performed. This analysis results in the establishment of <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> and <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> error bounds for the numerical solution. Furthermore, a variety of numerical examples are conducted to illustrate the accuracy of the EIEQ scheme, highlighting its effectiveness in addressing complex dynamical systems governed by the Allen–Cahn equation.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-09-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142312508","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
Generalized Multiscale Finite Element Method for discrete network (graph) models 离散网络(图)模型的广义多尺度有限元法
IF 2.1 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2024-09-13 DOI: 10.1016/j.cam.2024.116275

In this paper, we consider a time-dependent discrete network model with highly varying connectivity. The approximation by time is performed using an implicit scheme. We propose the coarse scale approximation construction of network models based on the Generalized Multiscale Finite Element Method. An accurate coarse-scale approximation is generated by solving local spectral problems in sub-networks. Convergence analysis of the proposed method is presented for semi-discrete and discrete network models. We establish the stability of the multiscale discrete network. Numerical results are presented for structured and random heterogeneous networks.

在本文中,我们考虑了一个具有高度变化连接性的随时间变化的离散网络模型。时间近似采用隐式方案。我们提出了基于广义多尺度有限元法的网络模型粗尺度近似构造。通过求解子网络中的局部谱问题,生成精确的粗尺度近似值。对半离散和离散网络模型进行了收敛分析。我们建立了多尺度离散网络的稳定性。还给出了结构化和随机异构网络的数值结果。
{"title":"Generalized Multiscale Finite Element Method for discrete network (graph) models","authors":"","doi":"10.1016/j.cam.2024.116275","DOIUrl":"10.1016/j.cam.2024.116275","url":null,"abstract":"<div><p>In this paper, we consider a time-dependent discrete network model with highly varying connectivity. The approximation by time is performed using an implicit scheme. We propose the coarse scale approximation construction of network models based on the Generalized Multiscale Finite Element Method. An accurate coarse-scale approximation is generated by solving local spectral problems in sub-networks. Convergence analysis of the proposed method is presented for semi-discrete and discrete network models. We establish the stability of the multiscale discrete network. Numerical results are presented for structured and random heterogeneous networks.</p></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142241449","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
Unconditional error analysis of the linearized transformed L1 virtual element method for nonlinear coupled time-fractional Schrödinger equations 非线性耦合时分数薛定谔方程线性化变换 L1 虚拟元素法的无条件误差分析
IF 2.1 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2024-09-13 DOI: 10.1016/j.cam.2024.116283
This paper constructs a linearized transformed L1 virtual element method for the generalized nonlinear coupled time-fractional Schrödinger equations. The solutions to such problems typically exhibit singular behavior at the beginning. To avoid this pitfall, we introduce an identical s-fractional differential system derived from a smoothing transformation of variables t=s1/α, 0<α<1. By utilizing the discrete complementary convolution kernels, we prove the boundedness and error estimates of the solution of time-discrete system. Moreover, the unconditionally optimal error bounds of the proposed fully discrete scheme are derived in L2-norm without restriction on the grid ratio. Finally, numerical tests on a set of polygonal meshes are presented to verify the theoretical results.
本文为广义非线性耦合时间分数薛定谔方程构建了线性化变换 L1 虚拟元素方法。此类问题的解通常在开始时表现出奇异行为。为了避免这一缺陷,我们引入了一个由变量 t=s1/α, 0<α<1 的平滑变换导出的相同 s 分式微分方程系统。此外,我们还在不限制网格比的情况下,以 L2 规范推导出了所提出的完全离散方案的无条件最优误差边界。最后,对一组多边形网格进行了数值测试,以验证理论结果。
{"title":"Unconditional error analysis of the linearized transformed L1 virtual element method for nonlinear coupled time-fractional Schrödinger equations","authors":"","doi":"10.1016/j.cam.2024.116283","DOIUrl":"10.1016/j.cam.2024.116283","url":null,"abstract":"<div><div>This paper constructs a linearized transformed <span><math><mrow><mi>L</mi><mn>1</mn></mrow></math></span> virtual element method for the generalized nonlinear coupled time-fractional Schrödinger equations. The solutions to such problems typically exhibit singular behavior at the beginning. To avoid this pitfall, we introduce an identical <span><math><mi>s</mi></math></span>-fractional differential system derived from a smoothing transformation of variables <span><math><mrow><mi>t</mi><mo>=</mo><msup><mrow><mi>s</mi></mrow><mrow><mn>1</mn><mo>/</mo><mi>α</mi></mrow></msup></mrow></math></span>, <span><math><mrow><mn>0</mn><mo>&lt;</mo><mi>α</mi><mo>&lt;</mo><mn>1</mn></mrow></math></span>. By utilizing the discrete complementary convolution kernels, we prove the boundedness and error estimates of the solution of time-discrete system. Moreover, the unconditionally optimal error bounds of the proposed fully discrete scheme are derived in <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span>-norm without restriction on the grid ratio. Finally, numerical tests on a set of polygonal meshes are presented to verify the theoretical results.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142312509","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 new space transformed finite element method for elliptic interface problems in Rn Rn 中椭圆界面问题的新空间变换有限元法
IF 2.1 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2024-09-12 DOI: 10.1016/j.cam.2024.116277

Interface problems, where distinct materials or physical domains meet, pose significant challenges in numerical simulations due to the discontinuities and sharp gradients across interfaces. Traditional finite element methods struggle to capture such behavior accurately. A new space transformed finite element method (ST-FEM) is developed for solving elliptic interface problems in Rn. A homeomorphic stretching transformation is introduced to obtain an equivalent problem in the transformed domain which can be solved easily, and the solution can be projected back to original domain by the inverse transformation. Compared with the existing methods, this new scheme has capability of handling discontinuities across the interface. The proposed approach has advantages in circumventing interface approximation properties and reducing the degree of freedom. We initially develop ST-FEM for elliptic problems and subsequently expand upon this concept to address elliptic interface problems. We prove optimal a priori error estimates in the H1 and L2 norms, and quasi-optimal error estimate for the maximum norm. Finally, numerical experiments demonstrate the superior accuracy and convergence properties of the ST-FEM when compared to the standard finite element method. The interface is assumed to be a (n1)-sphere, nevertheless, our analysis can cover symmetric domains such as an ellipsoid or a cylinder.

界面问题,即不同材料或物理域相遇的地方,由于界面上的不连续性和急剧梯度,给数值模拟带来了巨大挑战。传统的有限元方法难以准确捕捉这种行为。本文开发了一种新的空间变换有限元方法(ST-FEM),用于求解 Rn 中的椭圆界面问题。该方法引入了同构拉伸变换,从而在变换域中得到一个等效问题,该等效问题可以轻松求解,并且求解结果可以通过反变换投影回原始域。与现有方法相比,这一新方案具有处理跨界面不连续性的能力。所提出的方法在规避界面逼近特性和降低自由度方面具有优势。我们最初为椭圆问题开发了 ST-FEM,随后将这一概念扩展到椭圆界面问题。我们证明了 H1 和 L2 规范的最优先验误差估计,以及最大规范的准最优误差估计。最后,数值实验证明,与标准有限元方法相比,ST-FEM 具有更高的精度和收敛性。我们假设界面是一个 (n-1)- 球体,然而,我们的分析可以涵盖对称域,如椭圆体或圆柱体。
{"title":"A new space transformed finite element method for elliptic interface problems in Rn","authors":"","doi":"10.1016/j.cam.2024.116277","DOIUrl":"10.1016/j.cam.2024.116277","url":null,"abstract":"<div><p>Interface problems, where distinct materials or physical domains meet, pose significant challenges in numerical simulations due to the discontinuities and sharp gradients across interfaces. Traditional finite element methods struggle to capture such behavior accurately. A new space transformed finite element method (ST-FEM) is developed for solving elliptic interface problems in <span><math><msup><mrow><mi>R</mi></mrow><mrow><mi>n</mi></mrow></msup></math></span>. A homeomorphic stretching transformation is introduced to obtain an equivalent problem in the transformed domain which can be solved easily, and the solution can be projected back to original domain by the inverse transformation. Compared with the existing methods, this new scheme has capability of handling discontinuities across the interface. The proposed approach has advantages in circumventing interface approximation properties and reducing the degree of freedom. We initially develop ST-FEM for elliptic problems and subsequently expand upon this concept to address elliptic interface problems. We prove optimal a priori error estimates in the <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> and <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> norms, and quasi-optimal error estimate for the maximum norm. Finally, numerical experiments demonstrate the superior accuracy and convergence properties of the ST-FEM when compared to the standard finite element method. The interface is assumed to be a <span><math><mrow><mo>(</mo><mi>n</mi><mo>−</mo><mn>1</mn><mo>)</mo></mrow></math></span>-sphere, nevertheless, our analysis can cover symmetric domains such as an ellipsoid or a cylinder.</p></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-09-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142241453","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
Representation computation for the hypergeometric function of a Hermitian matrix argument 赫米矩阵参数的超几何函数的表示计算
IF 2.1 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2024-09-12 DOI: 10.1016/j.cam.2024.116258

We establish the exact expressions for the hypergeometric function of a Hermitian matrix argument. This result allows for the eigenvalues of the matrix argument to occur with arbitrary multiplicities and can be used for numerical computation. These exact expressions are particularly important since they provide the key ingredient which allows many results which involve this function to be useful from a practical engineering perspective.

我们建立了赫米特矩阵参数的超几何函数的精确表达式。这一结果允许矩阵参数的特征值以任意倍数出现,并可用于数值计算。这些精确表达式尤为重要,因为它们提供了关键要素,使许多涉及该函数的结果在实际工程中发挥作用。
{"title":"Representation computation for the hypergeometric function of a Hermitian matrix argument","authors":"","doi":"10.1016/j.cam.2024.116258","DOIUrl":"10.1016/j.cam.2024.116258","url":null,"abstract":"<div><p>We establish the exact expressions for the hypergeometric function of a Hermitian matrix argument. This result allows for the eigenvalues of the matrix argument to occur with arbitrary multiplicities and can be used for numerical computation. These exact expressions are particularly important since they provide the key ingredient which allows many results which involve this function to be useful from a practical engineering perspective.</p></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-09-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142241451","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
期刊
Journal of Computational and Applied Mathematics
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1