Applied Numerical Mathematics最新文献

英文中文

Numerical method and analysis for fluid-structure model on unbounded domains 无界域流固模型的数值方法与分析

IF 2.2 2区数学 Q1 MATHEMATICS, APPLIED

Applied Numerical Mathematics

Pub Date : 2025-06-25 DOI: 10.1016/j.apnum.2025.06.009

Hongwei Li, Xinyue Chen

Numerically solving the fluid-structure model on unbounded domains poses a challenge, due to the unbounded nature of the physical domain. To overcome this challenge, the artificial boundary method is specifically applied to numerically solve the fluid-structure model on unbounded domains, which can be used to analyze fluid-structure interactions in various scientific and engineering fields. Drawing inspiration from the artificial boundary method, we employ artificial boundaries to truncate the unbounded domain, subsequently designing the high order local artificial boundary conditions thereon based on the Padé approximation. Then, the initial value problem on the unbounded domain is reduced into an initial boundary value problem on the computational domain, which can be efficiently solved by adopting the finite difference method. Furthermore, a series of auxiliary variables is introduced specifically to address the issue of mixed derivatives arising in the artificial boundary conditions, and the stability, convergence and solvability of the reduced problem are rigorously analyzed. Numerical experiments are reported to demonstrate the effectiveness of artificial boundary conditions and theoretical analysis.

由于物理域的无界性，在无界域上对流固模型进行数值求解是一项挑战。为了克服这一挑战，专门采用人工边界法对无界域上的流固耦合模型进行数值求解，可用于各种科学和工程领域的流固耦合分析。受人工边界法的启发，我们采用人工边界截断无界域，然后基于pad近似在无界域上设计高阶局部人工边界条件。然后，将无界域上的初值问题转化为计算域上的初边值问题，采用有限差分法进行有效求解。此外，针对人工边界条件下的混合导数问题，引入了一系列辅助变量，并对简化后的问题的稳定性、收敛性和可解性进行了严格的分析。数值实验证明了人工边界条件和理论分析的有效性。

{"title":"Numerical method and analysis for fluid-structure model on unbounded domains","authors":"Hongwei Li, Xinyue Chen","doi":"10.1016/j.apnum.2025.06.009","DOIUrl":"10.1016/j.apnum.2025.06.009","url":null,"abstract":"<div><div>Numerically solving the fluid-structure model on unbounded domains poses a challenge, due to the unbounded nature of the physical domain. To overcome this challenge, the artificial boundary method is specifically applied to numerically solve the fluid-structure model on unbounded domains, which can be used to analyze fluid-structure interactions in various scientific and engineering fields. Drawing inspiration from the artificial boundary method, we employ artificial boundaries to truncate the unbounded domain, subsequently designing the high order local artificial boundary conditions thereon based on the Padé approximation. Then, the initial value problem on the unbounded domain is reduced into an initial boundary value problem on the computational domain, which can be efficiently solved by adopting the finite difference method. Furthermore, a series of auxiliary variables is introduced specifically to address the issue of mixed derivatives arising in the artificial boundary conditions, and the stability, convergence and solvability of the reduced problem are rigorously analyzed. Numerical experiments are reported to demonstrate the effectiveness of artificial boundary conditions and theoretical analysis.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"217 ","pages":"Pages 255-277"},"PeriodicalIF":2.2,"publicationDate":"2025-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144489531","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Numerical scheme for the solution of the “bad” Boussinesq equation “坏”Boussinesq方程解的数值格式

IF 2.2 2区数学 Q1 MATHEMATICS, APPLIED

Applied Numerical Mathematics

Pub Date : 2025-06-23 DOI: 10.1016/j.apnum.2025.06.011

Christophe Charlier , Daniel Eriksson , Jonatan Lenells

We present a numerical scheme for the solution of the initial-value problem for the “bad” Boussinesq equation. The accuracy of the scheme is tested by comparison with exact soliton solutions as well as with recently obtained asymptotic formulas for the solution.

给出了一类“坏”Boussinesq方程初值问题的数值解。通过与精确孤子解以及最近得到的解的渐近公式的比较，验证了该格式的准确性。

引用次数: 0

Numerical simulation of a dynamic human capital model with demographic delays via the local discrete Galerkin method 基于局部离散伽辽金方法的人口统计学时滞动态人力资本模型数值模拟

IF 2.2 2区数学 Q1 MATHEMATICS, APPLIED

Applied Numerical Mathematics

Pub Date : 2025-06-20 DOI: 10.1016/j.apnum.2025.06.007

Yadollah Ordokhani , Alireza Hosseinian , Pouria Assari

A strong and dynamic economy depends on various factors, with human capital playing a crucial role in fostering resilience and adaptability in an ever-changing world. Human capital, which depends on the past behavior of the system, requires strategic investments in education, health, and skill development. This study presents a numerical approach for solving the human capital model with age-structured delays, formulated as integro-differential equations with double delays and difference kernels. The proposed method employs a local meshless discrete Galerkin approach based on the moving least squares (MLS) technique, which can work with irregular or non-uniform data. The localized nature of the MLS scheme enhances computational efficiency by focusing on small neighborhoods. Moreover, the stabilized MLS framework, achieved by using shifted and scaled polynomial basis functions, enhances numerical stability and reduces sensitivity to the distribution of nodes, thereby transferring these advantageous properties to the method. The simplicity of the proposed algorithm makes it easy to implement on standard personal computers and extend to a wider class of delay integro-differential equations. To assess its reliability, we analyzed its error and determined the convergence order of the presented method. We have applied it to solve several numerical examples, and the obtained results confirm the method's accuracy, stability, and alignment with theoretical findings.

一个强大而充满活力的经济取决于各种因素，而人力资本在培养适应不断变化的世界的韧性和适应性方面发挥着至关重要的作用。人力资本取决于系统过去的行为，需要在教育、卫生和技能发展方面进行战略投资。本文提出了一种求解具有年龄结构延迟的人力资本模型的数值方法，该模型被表述为具有双延迟和差分核的积分-微分方程。该方法采用基于移动最小二乘（MLS）技术的局部无网格离散伽辽金方法，可以处理不规则或非均匀数据。MLS方案的局部特性通过关注小邻域来提高计算效率。此外，通过使用移位和缩放的多项式基函数实现的稳定MLS框架提高了数值稳定性，降低了对节点分布的敏感性，从而将这些优点转移到该方法中。该算法的简单性使其易于在标准个人计算机上实现，并扩展到更广泛的延迟积分微分方程。为了评估其可靠性，我们分析了该方法的误差，并确定了该方法的收敛顺序。应用该方法对若干数值算例进行了求解，得到的结果证实了该方法的准确性、稳定性以及与理论结果的一致性。

{"title":"Numerical simulation of a dynamic human capital model with demographic delays via the local discrete Galerkin method","authors":"Yadollah Ordokhani , Alireza Hosseinian , Pouria Assari","doi":"10.1016/j.apnum.2025.06.007","DOIUrl":"10.1016/j.apnum.2025.06.007","url":null,"abstract":"<div><div>A strong and dynamic economy depends on various factors, with human capital playing a crucial role in fostering resilience and adaptability in an ever-changing world. Human capital, which depends on the past behavior of the system, requires strategic investments in education, health, and skill development. This study presents a numerical approach for solving the human capital model with age-structured delays, formulated as integro-differential equations with double delays and difference kernels. The proposed method employs a local meshless discrete Galerkin approach based on the moving least squares (MLS) technique, which can work with irregular or non-uniform data. The localized nature of the MLS scheme enhances computational efficiency by focusing on small neighborhoods. Moreover, the stabilized MLS framework, achieved by using shifted and scaled polynomial basis functions, enhances numerical stability and reduces sensitivity to the distribution of nodes, thereby transferring these advantageous properties to the method. The simplicity of the proposed algorithm makes it easy to implement on standard personal computers and extend to a wider class of delay integro-differential equations. To assess its reliability, we analyzed its error and determined the convergence order of the presented method. We have applied it to solve several numerical examples, and the obtained results confirm the method's accuracy, stability, and alignment with theoretical findings.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"217 ","pages":"Pages 234-254"},"PeriodicalIF":2.2,"publicationDate":"2025-06-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144489530","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

Note on an inertial projection-based approach for solving extended general quasi-variational inequalities and its convergence analysis 基于惯性投影的广义拟变分不等式求解方法及其收敛性分析

IF 2.2 2区数学 Q1 MATHEMATICS, APPLIED

Applied Numerical Mathematics

Pub Date : 2025-06-20 DOI: 10.1016/j.apnum.2025.06.012

Saudia Jabeen , Siegfried Macías , Saleem Ullah , Muhammad Aslam Noor , Jorge E. Macías-Díaz

In this work, an inertial projection-based method is proposed to find approximate solutions to a new class of quasi-variational inequalities. The approach utilizes projection techniques to develop an iterative scheme, and its convergence properties are rigorously examined under appropriate conditions. Various specific cases are derived from the general framework, highlighting the adaptability of the proposed method. The comparative performance of this approach with existing techniques remains an open area for exploration. It is expected that the theoretical framework and techniques proposed in this study will stimulate further research in this domain.

在这项工作中，提出了一种基于惯性投影的方法来求一类新的拟变分不等式的近似解。该方法利用投影技术开发迭代格式，并在适当的条件下严格检验了其收敛性。从总体框架中衍生出各种具体案例，突出了所提方法的适应性。这种方法与现有技术的性能比较仍然是一个有待探索的开放领域。期望本研究提出的理论框架和技术将促进该领域的进一步研究。

引用次数: 0

Analysis of Tseng algorithm with inertial extrapolation step for stochastic variational inequality problem 随机变分不等式问题的惯性外推Tseng算法分析

IF 2.2 2区数学 Q1 MATHEMATICS, APPLIED

Applied Numerical Mathematics

Pub Date : 2025-06-20 DOI: 10.1016/j.apnum.2025.06.008

Olawale K. Oyewole , Seithuti P. Moshokoa , Sani Salisu , Yekini Shehu

In this paper, we design an inertial version of the Tseng extragradient algorithm (also called the Forward-Backward-Forward Algorithm) with self-adaptive step sizes to solve the stochastic variational inequality problem. We prove that the sequence of iterates generated by our proposed algorithm converges to a solution of the stochastic variational inequality problem under mild conditions. Furthermore, we obtain some convergence rates and numerical simulations of our proposed algorithm with comparisons with related algorithms to show the superiority of our algorithm.

本文设计了一种具有自适应步长的惯性版本的Tseng extraggradient算法（也称为Forward-Backward-Forward算法）来解决随机变分不等式问题。证明了该算法生成的迭代序列在温和条件下收敛于随机变分不等式问题的一个解。最后给出了算法的一些收敛速度和数值模拟，并与相关算法进行了比较，证明了算法的优越性。

引用次数: 0

An efficient multi projection methods for systems of Fredholm integral equations with mixed weakly singular kernels: A superconvergence approach 混合弱奇异核Fredholm积分方程组的一种有效多投影方法：一种超收敛方法

IF 2.2 2区数学 Q1 MATHEMATICS, APPLIED

Applied Numerical Mathematics

Pub Date : 2025-06-19 DOI: 10.1016/j.apnum.2025.06.006

Krishna Murari Malav , Kapil Kant , Joydip Dhar , Samiran Chakraborty

In this article, we develop the multi-Galerkin and iterated multi-Galerkin methods to solve systems of second-kind linear Fredholm integral equations (FIEs) with smooth and mixed weakly singular kernels. First, we develop the mathematical formulation of the multi-Galerkin and iterated multi-Galerkin methods using piecewise polynomial approximations to solve such systems and obtain superconvergence results. These methods transform the linear system of FIEs into corresponding matrix equations. We derive error estimates and obtain the convergence analysis. We prove that the convergence rates for the multi-Galerkin method are

O (h^{3 r})

for smooth kernels and

O (h^{1 + r - α} \log h)

for mixed weakly singular kernels, where r denote the degree of the piecewise polynomials, h is the norm of partitions and

α = \max_{i, j} α_{i j}

. Moreover, we establish that the iterated multi-Galerkin method achieves improved convergence rates of

O (h^{4 r})

for smooth kernels and

O (h^{r + 2 (1 - α)} {(\log h)}^{2})

for mixed weakly singular kernels. Hence, the results show that the iterated multi-Galerkin method improves the multi-Galerkin method. Finally, the theoretical results are validated through the numerical examples.

利用多重galerkin和迭代多重galerkin方法，求解了一类具有光滑和混合弱奇异核的第二类线性Fredholm积分方程组。首先，我们利用分段多项式逼近建立了多重galerkin和迭代多重galerkin方法的数学表达式，并得到了超收敛结果。这些方法将FIEs的线性方程组转化为相应的矩阵方程。我们推导了误差估计，并得到了收敛性分析。证明了多重伽辽金方法对于光滑核的收敛速度为O(h3r)，对于混合弱奇异核的收敛速度为O(h1+r−αlog (h))，其中r表示分段多项式的阶数，h表示分区范数，α=maxi,j∈αij。此外，我们还证明了迭代多重伽辽金方法对于光滑核的收敛速度为O(h4r)，对于混合弱奇异核的收敛速度为O(hr+2（1−α)(log (h)2)）。结果表明，迭代多重伽辽金方法是对多重伽辽金方法的改进。最后通过数值算例对理论结果进行了验证。

{"title":"An efficient multi projection methods for systems of Fredholm integral equations with mixed weakly singular kernels: A superconvergence approach","authors":"Krishna Murari Malav , Kapil Kant , Joydip Dhar , Samiran Chakraborty","doi":"10.1016/j.apnum.2025.06.006","DOIUrl":"10.1016/j.apnum.2025.06.006","url":null,"abstract":"<div><div>In this article, we develop the multi-Galerkin and iterated multi-Galerkin methods to solve systems of second-kind linear Fredholm integral equations (FIEs) with smooth and mixed weakly singular kernels. First, we develop the mathematical formulation of the multi-Galerkin and iterated multi-Galerkin methods using piecewise polynomial approximations to solve such systems and obtain superconvergence results. These methods transform the linear system of FIEs into corresponding matrix equations. We derive error estimates and obtain the convergence analysis. We prove that the convergence rates for the multi-Galerkin method are <span><math><mi>O</mi><mo>(</mo><msup><mrow><mi>h</mi></mrow><mrow><mn>3</mn><mi>r</mi></mrow></msup><mo>)</mo></math></span> for smooth kernels and <span><math><mi>O</mi><mo>(</mo><msup><mrow><mi>h</mi></mrow><mrow><mn>1</mn><mo>+</mo><mi>r</mi><mo>−</mo><mi>α</mi></mrow></msup><mi>log</mi><mo>⁡</mo><mi>h</mi><mo>)</mo></math></span> for mixed weakly singular kernels, where <em>r</em> denote the degree of the piecewise polynomials, <em>h</em> is the norm of partitions and <span><math><mi>α</mi><mo>=</mo><munder><mi>max</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></munder><mo>⁡</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></math></span>. Moreover, we establish that the iterated multi-Galerkin method achieves improved convergence rates of <span><math><mi>O</mi><mo>(</mo><msup><mrow><mi>h</mi></mrow><mrow><mn>4</mn><mi>r</mi></mrow></msup><mo>)</mo></math></span> for smooth kernels and <span><math><mi>O</mi><mo>(</mo><msup><mrow><mi>h</mi></mrow><mrow><mi>r</mi><mo>+</mo><mn>2</mn><mo>(</mo><mn>1</mn><mo>−</mo><mi>α</mi><mo>)</mo></mrow></msup><msup><mrow><mo>(</mo><mi>log</mi><mo>⁡</mo><mi>h</mi><mo>)</mo></mrow><mrow><mn>2</mn></mrow></msup><mo>)</mo></math></span> for mixed weakly singular kernels. Hence, the results show that the iterated multi-Galerkin method improves the multi-Galerkin method. Finally, the theoretical results are validated through the numerical examples.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"217 ","pages":"Pages 172-189"},"PeriodicalIF":2.2,"publicationDate":"2025-06-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144322646","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

Localization of tumor through a non-conventional numerical shape optimization technique 利用非常规数值形状优化技术对肿瘤进行定位

IF 2.2 2区数学 Q1 MATHEMATICS, APPLIED

Applied Numerical Mathematics

Pub Date : 2025-06-16 DOI: 10.1016/j.apnum.2025.06.005

Julius Fergy Tiongson Rabago

This paper presents a method for estimating the shape and location of an embedded tumor using shape optimization techniques, specifically through the coupled complex boundary method. The inverse problem—characterized by a measured temperature profile and corresponding heat flux (e.g., from infrared thermography)—is reformulated as a complex boundary value problem with a complex Robin boundary condition, thereby simplifying its over-specified nature. The geometry of the tumor is identified by optimizing an objective functional that depends on the imaginary part of the solution throughout the domain. The shape derivative of the functional is derived through shape sensitivity analysis. An iterative algorithm is developed to numerically recover the tumor shape, based on the Riesz representative of the gradient and implemented using the finite element method. In addition, the mesh sensitivity of the finite element solution to the state problem is analyzed, and bounds are established for its variation with respect to mesh deformation and its gradient. Numerical examples are presented to validate the theoretical results and to demonstrate the accuracy and effectiveness of the proposed method.

本文提出了一种利用形状优化技术，特别是通过耦合复杂边界方法来估计嵌入肿瘤的形状和位置的方法。反问题——以测量的温度剖面和相应的热通量（例如，来自红外热成像）为特征——被重新表述为具有复杂Robin边界条件的复杂边值问题，从而简化了其过度指定的性质。通过优化目标函数来识别肿瘤的几何形状，目标函数依赖于整个域内解的虚部。通过形状灵敏度分析，导出了泛函的形状导数。提出了一种基于Riesz梯度表示的肿瘤形状数值恢复迭代算法，并采用有限元法实现。此外，分析了状态问题有限元解的网格灵敏度，并建立了其随网格变形及其梯度变化的边界。数值算例验证了理论结果，验证了所提方法的准确性和有效性。

{"title":"Localization of tumor through a non-conventional numerical shape optimization technique","authors":"Julius Fergy Tiongson Rabago","doi":"10.1016/j.apnum.2025.06.005","DOIUrl":"10.1016/j.apnum.2025.06.005","url":null,"abstract":"<div><div>This paper presents a method for estimating the shape and location of an embedded tumor using shape optimization techniques, specifically through the coupled complex boundary method. The inverse problem—characterized by a measured temperature profile and corresponding heat flux (e.g., from infrared thermography)—is reformulated as a complex boundary value problem with a complex Robin boundary condition, thereby simplifying its over-specified nature. The geometry of the tumor is identified by optimizing an objective functional that depends on the imaginary part of the solution throughout the domain. The shape derivative of the functional is derived through shape sensitivity analysis. An iterative algorithm is developed to numerically recover the tumor shape, based on the Riesz representative of the gradient and implemented using the finite element method. In addition, the mesh sensitivity of the finite element solution to the state problem is analyzed, and bounds are established for its variation with respect to mesh deformation and its gradient. Numerical examples are presented to validate the theoretical results and to demonstrate the accuracy and effectiveness of the proposed method.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"217 ","pages":"Pages 135-171"},"PeriodicalIF":2.2,"publicationDate":"2025-06-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144307834","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 high-order accurate unconditionally stable bound-preserving numerical scheme for the Cahn-Hilliard-Navier-Stokes equations Cahn-Hilliard-Navier-Stokes方程的高阶精确无条件稳定保界数值格式

IF 2.2 2区数学 Q1 MATHEMATICS, APPLIED

Applied Numerical Mathematics

Pub Date : 2025-06-10 DOI: 10.1016/j.apnum.2025.06.004

Yali Gao , Daozhi Han , Sayantan Sarkar

A high order numerical method is developed for solving the Cahn-Hilliard-Navier-Stokes equations with the Flory-Huggins potential. The scheme is based on the

Q_{k}

finite element with mass lumping on rectangular grids, the second-order convex splitting method and the pressure correction method. The unique solvability, unconditional stability, and bound-preserving properties are rigorously established. The key for bound-preservation is the discrete

L^{1}

estimate of the singular potential. Ample numerical experiments are performed to validate the desired properties of the proposed numerical scheme.

提出了一种求解具有Flory-Huggins势的Cahn-Hilliard-Navier-Stokes方程的高阶数值方法。该方案基于矩形网格质量集总的Qk有限元、二阶凸分裂法和压力修正法。严格地建立了唯一可解性、无条件稳定性和保界性。保持边界的关键是奇异势的离散L1估计。进行了大量的数值实验来验证所提出的数值格式的预期性能。

引用次数: 0

Analysis of a divergence-free element-free Galerkin method for the Navier-Stokes equations Navier-Stokes方程的无散度无单元Galerkin方法分析

IF 2.2 2区数学 Q1 MATHEMATICS, APPLIED

Applied Numerical Mathematics

Pub Date : 2025-06-10 DOI: 10.1016/j.apnum.2025.06.002

Xiaolin Li , Haiyun Dong

In this paper, an efficient divergence-free element-free Galerkin (DFEFG) method is proposed for the numerical analysis of the incompressible Navier-Stokes equations. In this method, a divergence-free moving least squares (DFMLS) approximation is used to obtain the meshless approximation of the divergence-free velocity field. The properties, stability and error of the DFMLS approximation are analyzed firstly, and then the stability and error estimation of the DFEFG method are derived theoretically. Finally, numerical results demonstrate the efficiency of the proposed methods and verify the theoretical analysis.

本文提出了一种有效的无散度无单元伽辽金（DFEFG）方法，用于不可压缩Navier-Stokes方程的数值分析。该方法采用无散度移动最小二乘（DFMLS）近似，得到无散度速度场的无网格近似。首先分析了DFMLS近似的性质、稳定性和误差，然后从理论上推导了DFEFG方法的稳定性和误差估计。最后，数值结果验证了所提方法的有效性和理论分析的正确性。

引用次数: 0

Generalized Clenshaw-Curtis quadrature method for systems of linear ODEs with constant coefficients 常系数线性微分方程系统的广义Clenshaw-Curtis正交法

IF 2.2 2区数学 Q1 MATHEMATICS, APPLIED

Applied Numerical Mathematics

Pub Date : 2025-06-09 DOI: 10.1016/j.apnum.2025.06.003

Fu-Rong Lin, Xi Yang, Gui-Rong Zhang

In this paper, we consider high precision numerical methods for the initial problem of systems of linear ordinary differential equations (ODEs) with constant coefficients. It is well-known that the analytic solution of such a system of linear ODEs involves a matrix exponential function and an integral whose integrand is the product of a matrix exponential and a vector-valued function. We mainly consider numerical quadrature methods for the integral term in the analytic solution and propose a generalized Clenshaw-Curtis (GCC) quadrature method. The proposed method is then applied to the initial-boundary value problem for a heat conduction equation and a Riesz space fractional diffusion equation, respectively. Numerical results are presented to demonstrate the effectiveness of the proposed method.

本文研究了常系数线性常微分方程系统初始问题的高精度数值解法。众所周知，这种线性微分方程系统的解析解涉及一个矩阵指数函数和一个积分，其被积是一个矩阵指数函数和一个向量值函数的乘积。本文主要研究了解析解中积分项的数值求积分方法，提出了一种广义的clclenshaw - curtis求积分方法。然后将该方法分别应用于热传导方程和Riesz空间分数扩散方程的初边值问题。数值结果验证了该方法的有效性。

引用次数: 0

首页上一页

下一页尾页

类型

全部化学•材料生命科学医学物理工程技术环境•农林材料科学地球科学法学管理学化学环境科学与生态学计算机科学教育学经济学农林科学人文科学生物学数学物理与天体物理心理学综合性期刊其他工业工程理学历史学农学文学信息工程

数据库

全部 ACS Publications Elsevier ieeexplore Springer The Royal Society of Chemistry Wiley

期刊

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.

﹀