Computers & Mathematics with Applications最新文献_第4页

ITPINN : Integral-trainable physics informed neural network for solving high-dimensional evolution non-local partial differential equations 求解高维演化非局部偏微分方程的积分可训练物理通知神经网络

IF 2.5 2区数学 Q1 MATHEMATICS, APPLIED

Computers & Mathematics with Applications

Pub Date : 2026-01-21 DOI: 10.1016/j.camwa.2026.01.019

Wenkai Liu , Fanhai Zeng , Hong Li , Yang Liu

In this paper, we develop a new physics informed neural network (PINN) method, named integral-trainable PINN (ITPINN), to solve high-dimensional non-local partial differential equations (PDEs), involving PDEs with fractional derivatives (Caputo derivative and Riemann-Liouville derivative) or multiple integral. In the ITPINN framework, we perform integration by parts on the original integral to obtain a new constraint condition, which forms a coupled system with the original equation. We consider the integral terms as unknown functions in the coupled system and construct a neural network with three output terms, one for predicting the exact solution, one for predicting the original integral term, and one for approximating the new integral obtained by integration by parts. The network is used as a surrogate model for fractional derivatives or multiple integral, which allows approximation of the fractional derivatives or multiple integral to be achieved by training the network. The proposed method omits the process of discretizing the integral term using traditional numerical methods, such as finite difference method or interpolation approximation. Moreover, the physical information obtained from integration by parts is used to construct a new supervised learning task to further constrain the surrogate model for the integral terms. Several experiments are used to illustrate the performance of the ITPINN. The numerical results confirm that our proposed method can effectively solve high-dimensional evolution non-local PDEs, such as 50D problems. Compared to fractional PINN (fPINN) and auxiliary PINN (A-PINN), the ITPINN can achieve higher prediction accuracy and save more training time. In particular, we also test the robustness of the ITPINN under interference with noise intensities ranging from 0.01% to 50% and further discuss its scalability in 100D and 1000D problems.

本文提出了一种新的物理通知神经网络（PINN）方法，称为积分可训练神经网络（ITPINN），用于求解含有分数阶导数（Caputo导数和Riemann-Liouville导数）或多重积分的高维非局部偏微分方程（PDEs）。在ITPINN框架中，我们对原积分进行分部积分，得到一个新的约束条件，与原方程形成一个耦合系统。我们将积分项视为耦合系统中的未知函数，并构造了一个具有三个输出项的神经网络，一个用于预测精确解，一个用于预测原始积分项，一个用于逼近由分部积分得到的新积分。该网络被用作分数阶导数或多重积分的代理模型，它允许通过训练网络来实现分数阶导数或多重积分的近似。该方法省去了用有限差分法或插值逼近等传统数值方法对积分项进行离散化的过程。此外，利用分部积分获得的物理信息构建新的监督学习任务，进一步约束积分项的代理模型。用几个实验说明了ITPINN的性能。数值结果表明，该方法可以有效地求解高维演化非局部偏微分方程，如50D问题。与分数PINN （fPINN）和辅助PINN （A-PINN）相比，ITPINN可以达到更高的预测精度，节省更多的训练时间。特别是，我们还测试了ITPINN在0.01%至50%噪声强度干扰下的鲁棒性，并进一步讨论了其在100D和1000D问题中的可扩展性。

{"title":"ITPINN : Integral-trainable physics informed neural network for solving high-dimensional evolution non-local partial differential equations","authors":"Wenkai Liu , Fanhai Zeng , Hong Li , Yang Liu","doi":"10.1016/j.camwa.2026.01.019","DOIUrl":"10.1016/j.camwa.2026.01.019","url":null,"abstract":"<div><div>In this paper, we develop a new physics informed neural network (PINN) method, named integral-trainable PINN (ITPINN), to solve high-dimensional non-local partial differential equations (PDEs), involving PDEs with fractional derivatives (Caputo derivative and Riemann-Liouville derivative) or multiple integral. In the ITPINN framework, we perform integration by parts on the original integral to obtain a new constraint condition, which forms a coupled system with the original equation. We consider the integral terms as unknown functions in the coupled system and construct a neural network with three output terms, one for predicting the exact solution, one for predicting the original integral term, and one for approximating the new integral obtained by integration by parts. The network is used as a surrogate model for fractional derivatives or multiple integral, which allows approximation of the fractional derivatives or multiple integral to be achieved by training the network. The proposed method omits the process of discretizing the integral term using traditional numerical methods, such as finite difference method or interpolation approximation. Moreover, the physical information obtained from integration by parts is used to construct a new supervised learning task to further constrain the surrogate model for the integral terms. Several experiments are used to illustrate the performance of the ITPINN. The numerical results confirm that our proposed method can effectively solve high-dimensional evolution non-local PDEs, such as 50D problems. Compared to fractional PINN (fPINN) and auxiliary PINN (A-PINN), the ITPINN can achieve higher prediction accuracy and save more training time. In particular, we also test the robustness of the ITPINN under interference with noise intensities ranging from 0.01% to 50% and further discuss its scalability in 100D and 1000D problems.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"205 ","pages":"Pages 282-317"},"PeriodicalIF":2.5,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146014491","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 parameter-free approach for 3D staggered Lagrangian hydrodynamics 三维交错拉格朗日流体力学的无参数方法

IF 2.5 2区数学 Q1 MATHEMATICS, APPLIED

Computers & Mathematics with Applications

Pub Date : 2026-01-20 DOI: 10.1016/j.camwa.2026.01.002

Xihua Xu , Xuan Zhou

In engineering applications, excessive artificial parameters can consume significant computational resources, particularly in 3D. This study presents a novel parameter-free approach for 3D staggered Lagrangian hydrodynamics that eliminates the need for artificial parameters and offers advantages over traditional methods. We construct a new form of artificial viscosity and give a new explanation for preventing the hourglass phenomenon. The proposed scheme ensures conservation of total mass, total momentum, and total energy through its unique formulation. Extensive testing has demonstrated the high robustness of this scheme, making it well-suited for multi-physics problems and various engineering applications.

在工程应用中，过多的人工参数会消耗大量的计算资源，特别是在3D中。该研究提出了一种新的三维交错拉格朗日流体动力学无参数方法，消除了对人工参数的需要，具有传统方法无法比拟的优点。我们构造了一种新的人工黏度形式，并对防止沙漏现象给出了新的解释。该方案通过其独特的公式保证了总质量、总动量和总能量的守恒。大量的测试证明了该方案的高鲁棒性，使其非常适合于多物理场问题和各种工程应用。

引用次数: 0

MS-PINN: A physics-informed neural network for multi-field coupled evolution modeling in metal solidification MS-PINN：用于金属凝固多场耦合演化建模的物理信息神经网络

IF 2.5 2区数学 Q1 MATHEMATICS, APPLIED

Computers & Mathematics with Applications

Pub Date : 2026-01-20 DOI: 10.1016/j.camwa.2026.01.015

Chen Bai , Yunhu Zhang , Hongxing Zheng , Quan Qian

Simulating the metal solidification process is crucial for improving product quality, optimizing manufacturing processes, and developing new materials. Traditional numerical methods like the Finite Element Method (FEM) and Physics-Informed Neural Networks (PINNs) face significant challenges when applied to metal solidification simulations due to inefficiencies and inaccuracies in dealing with multiphysics coupling, nonlinearity, and spatio-temporal complexity. Despite their potential, PINNs require further optimization to accurately capture complex physical phenomena in practical simulations. In this study, we propose a novel method based on PINNs, termed MS-PINN, which integrates Fourier Feature Embedding (FFE), Residual-based Adaptive Resampling (RAD), and Self-adaptive Loss Balanced methods (SAL) to significantly enhance simulation accuracy and efficiency. FFE improves the model’s ability to capture high-frequency features, RAD increases learning efficiency in high-gradient regions, and SAL dynamically adjusts loss function weights to optimize the training process. Experimental results show that MS-PINN outperforms traditional PINNs and other advanced approaches, achieving average error reductions of approximately 81.00% compared to Conv-LSTM, 77.11% compared to TCN, and 61.56% compared to PINN in reconstruction experiments. In predictive experiments, MS-PINN reduces errors by 53.23%, 68.81%, and 72.54% compared to PINN, TCN, and CONV-LSTM methods, respectively. Additionally, we developed a general PDE-solving software, NeuroPDE, based on this method. NeuroPDE has demonstrated success not only in the solidification process of Cu-1wt.%Ag alloy but also in solving Burgers, diffusion, and Navier-Stokes (NS) equations, including turbulent datasets characterized by high Reynolds numbers, and their inverse problems. This highlights NeuroPDE’s versatility and broad applicability in solving complex forward and inverse problems in fluid dynamics and other fields.

模拟金属凝固过程对于提高产品质量、优化制造工艺和开发新材料至关重要。传统的数值方法，如有限元法（FEM）和物理信息神经网络（pinn），由于在处理多物理场耦合、非线性和时空复杂性方面的低效率和不准确性，在应用于金属凝固模拟时面临着重大挑战。尽管具有潜力，但在实际模拟中需要进一步优化才能准确捕获复杂的物理现象。在这项研究中，我们提出了一种基于pinn的新方法，称为MS-PINN，它集成了傅里叶特征嵌入（FFE），基于残差的自适应重采样（RAD）和自适应损失平衡方法（SAL），以显着提高仿真精度和效率。FFE提高了模型捕获高频特征的能力，RAD提高了高梯度区域的学习效率，SAL动态调整损失函数权值以优化训练过程。实验结果表明，MS-PINN优于传统的PINN和其他先进的方法，在重建实验中，与convl - lstm相比，MS-PINN的平均误差降低约81.00%，与TCN相比，MS-PINN的平均误差降低约77.11%，与PINN相比，MS-PINN的平均误差降低约61.56%。在预测实验中，MS-PINN比PINN、TCN和convl - lstm方法分别降低了53.23%、68.81%和72.54%的误差。此外，我们还基于该方法开发了通用的pde求解软件NeuroPDE。NeuroPDE不仅在Cu-1wt的凝固过程中取得了成功。%Ag合金，但也在解决汉堡，扩散，和Navier-Stokes （NS）方程，包括湍流数据集特征的高雷诺数，和他们的反问题。这突出了NeuroPDE在解决流体动力学和其他领域复杂的正反问题方面的多功能性和广泛适用性。

{"title":"MS-PINN: A physics-informed neural network for multi-field coupled evolution modeling in metal solidification","authors":"Chen Bai , Yunhu Zhang , Hongxing Zheng , Quan Qian","doi":"10.1016/j.camwa.2026.01.015","DOIUrl":"10.1016/j.camwa.2026.01.015","url":null,"abstract":"<div><div>Simulating the metal solidification process is crucial for improving product quality, optimizing manufacturing processes, and developing new materials. Traditional numerical methods like the Finite Element Method (FEM) and Physics-Informed Neural Networks (PINNs) face significant challenges when applied to metal solidification simulations due to inefficiencies and inaccuracies in dealing with multiphysics coupling, nonlinearity, and spatio-temporal complexity. Despite their potential, PINNs require further optimization to accurately capture complex physical phenomena in practical simulations. In this study, we propose a novel method based on PINNs, termed MS-PINN, which integrates Fourier Feature Embedding (FFE), Residual-based Adaptive Resampling (RAD), and Self-adaptive Loss Balanced methods (SAL) to significantly enhance simulation accuracy and efficiency. FFE improves the model’s ability to capture high-frequency features, RAD increases learning efficiency in high-gradient regions, and SAL dynamically adjusts loss function weights to optimize the training process. Experimental results show that MS-PINN outperforms traditional PINNs and other advanced approaches, achieving average error reductions of approximately 81.00% compared to Conv-LSTM, 77.11% compared to TCN, and 61.56% compared to PINN in reconstruction experiments. In predictive experiments, MS-PINN reduces errors by 53.23%, 68.81%, and 72.54% compared to PINN, TCN, and CONV-LSTM methods, respectively. Additionally, we developed a general PDE-solving software, NeuroPDE, based on this method. NeuroPDE has demonstrated success not only in the solidification process of Cu-1wt.%Ag alloy but also in solving Burgers, diffusion, and Navier-Stokes (NS) equations, including turbulent datasets characterized by high Reynolds numbers, and their inverse problems. This highlights NeuroPDE’s versatility and broad applicability in solving complex forward and inverse problems in fluid dynamics and other fields.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"207 ","pages":"Pages 60-78"},"PeriodicalIF":2.5,"publicationDate":"2026-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145996321","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

Mean-square contractivity and convergence rate of stochastic theta Milstein schemes for non-autonomous SDEs with non-globally Lipschitz diffusion coefficients 具有非全局Lipschitz扩散系数的非自治SDEs随机theta Milstein格式的均方收缩性和收敛速度

IF 2.5 2区数学 Q1 MATHEMATICS, APPLIED

Computers & Mathematics with Applications

Pub Date : 2026-01-20 DOI: 10.1016/j.camwa.2026.01.007

Jinran Yao, Zhengwei Yin, Luozhong Gong

This study is dedicated to investigating the mean-square contractivity and convergence rate of stochastic theta Milstein (STM) schemes for non-autonomous stochastic differential equations (SDEs) in the non-globally Lipschitz setting. More precisely, we begin by analyzing the mean-square contractivity of the numerical discretization for the STM schemes with θ ∈ [1/2, 1] under the coupled monotonicity condition along with the polynomial growth condition with respect to the state variable. Additionally, we employ an iterative strategy to deduce the strong convergence order of the STM schemes under non-globally Lipschitz conditions. It is demonstrated that the numerical approximations of the STM schemes with θ ∈ [1/2, 1] converges strongly with order one in the mean-square sense to the exact solution of the considered SDEs, under Assumptions 2.1 and 4.1 with respect to both the spatial and time variables. In the present study, the convergence analysis specifically addresses the case where

L \in R

in the coupled monotonicity condition, rather than L ≥ 0, as commonly presented in the literature. This analysis does not count upon a priori high-order moment estimates of the considered numerical schemes in establishing the convergence theorem. Finally, numerical experiments are performed to empirically substantiate the theoretical findings.

本文研究了非全局Lipschitz环境下非自治随机微分方程（SDEs）的随机theta Milstein （STM）格式的均方收缩性和收敛率。更准确地说，我们首先分析了关于状态变量的耦合单调性条件下，具有θ ∈ [1/2,1]的STM格式的数值离散化的均方收缩性。此外，我们采用迭代策略推导了非全局Lipschitz条件下STM格式的强收敛阶。证明了在假设2.1和4.1条件下，具有θ ∈ [1/2,1]的STM格式的数值逼近在均方意义上强收敛于所考虑的SDEs的空间和时间变量的精确解，且收敛速度为1阶。在本研究中，收敛性分析专门针对耦合单调性条件下L∈R的情况，而不是文献中常见的L ≥ 0。这种分析不依赖于在建立收敛定理时所考虑的数值格式的先验高阶矩估计。最后，通过数值实验对理论结果进行了实证验证。

{"title":"Mean-square contractivity and convergence rate of stochastic theta Milstein schemes for non-autonomous SDEs with non-globally Lipschitz diffusion coefficients","authors":"Jinran Yao, Zhengwei Yin, Luozhong Gong","doi":"10.1016/j.camwa.2026.01.007","DOIUrl":"10.1016/j.camwa.2026.01.007","url":null,"abstract":"<div><div>This study is dedicated to investigating the mean-square contractivity and convergence rate of stochastic theta Milstein (STM) schemes for non-autonomous stochastic differential equations (SDEs) in the non-globally Lipschitz setting. More precisely, we begin by analyzing the mean-square contractivity of the numerical discretization for the STM schemes with <em>θ</em> ∈ [1/2, 1] under the coupled monotonicity condition along with the polynomial growth condition with respect to the state variable. Additionally, we employ an iterative strategy to deduce the strong convergence order of the STM schemes under non-globally Lipschitz conditions. It is demonstrated that the numerical approximations of the STM schemes with <em>θ</em> ∈ [1/2, 1] converges strongly with order one in the mean-square sense to the exact solution of the considered SDEs, under Assumptions 2.1 and 4.1 with respect to both the spatial and time variables. In the present study, the convergence analysis specifically addresses the case where <span><math><mrow><mi>L</mi><mo>∈</mo><mi>R</mi></mrow></math></span> in the coupled monotonicity condition, rather than <em>L</em> ≥ 0, as commonly presented in the literature. This analysis does not count upon a priori high-order moment estimates of the considered numerical schemes in establishing the convergence theorem. Finally, numerical experiments are performed to empirically substantiate the theoretical findings.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"206 ","pages":"Pages 228-256"},"PeriodicalIF":2.5,"publicationDate":"2026-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146014831","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 uniform superconvergent error estimates of anisotropic nonconforming energy-dissipative SAV-CN FEM for fourth-order singular perturbation Bi-flux diffusion model 四阶奇异摄动双通量扩散模型各向异性非协调能量耗散SAV-CN有限元的无条件一致超收敛误差估计

IF 2.5 2区数学 Q1 MATHEMATICS, APPLIED

Computers & Mathematics with Applications

Pub Date : 2026-01-20 DOI: 10.1016/j.camwa.2026.01.012

Dongyang Shi , Sihui Zhang

This investigation aims to delve into the unconditional uniform superconvergence behavior of the scalar auxiliary variable (SAV) Crank-Nicolson (CN) fully discrete scheme for the fourth-order singular perturbation Bi-flux diffusion equation with anisotropic constrained nonconforming rotated

Q_{1}

(

C N Q_{1}^{r o t}

) finite element method (FEM). Innovative high-accuracy estimates tailored specifically for the

C N Q_{1}^{r o t}

element on anisotropic meshes are established as the cornerstone for achieving superconvergent outcomes. Following this, a novel SAV-CN scheme is formulated, with the energy dissipation theoretically given to guarantee that the numerical solutions remain bounded in the broken H¹-norm. Then, by leveraging the aforementioned novel estimates and sidestepping the reliance on inverse inequalities, the unconditional uniform superclose result is rigorously derived, which transcends the constraints of the negative power of the perturbation parameter

A = \sqrt{β (1 - β) R}

and the limitations imposed by the ratio between the temporal step τ and spatial partition parameter h. Furthermore, an application of the economic interpolation post-processing technique enables the procurement of superconvergence estimates for the devised numerical schemes. Ultimately, the theoretical analysis is strengthened through the conduct of numerical experiments.

研究了四阶各向异性约束非协调旋转Q1 （CNQ1rot）有限元法四阶奇异摄动双通量扩散方程的标量辅助变量（SAV） Crank-Nicolson （CN）完全离散格式的无条件一致超收敛行为。在各向异性网格上为CNQ1rot元素量身定制的创新高精度估计是实现超收敛结果的基石。在此基础上，提出了一种新的SAV-CN格式，并在理论上给出了能量耗散，以保证数值解在h1 -范数破碎时保持有界。然后，利用上述新估计，避开对逆不等式的依赖，严格推导了无条件均匀超接近结果，该结果超越了扰动参数A=β(1−β)R的负幂约束和时间步长τ与空间分割参数h之比的限制。经济插值后处理技术的应用使所设计的数值格式获得了超收敛估计。最后，通过数值实验加强理论分析。

{"title":"Unconditional uniform superconvergent error estimates of anisotropic nonconforming energy-dissipative SAV-CN FEM for fourth-order singular perturbation Bi-flux diffusion model","authors":"Dongyang Shi , Sihui Zhang","doi":"10.1016/j.camwa.2026.01.012","DOIUrl":"10.1016/j.camwa.2026.01.012","url":null,"abstract":"<div><div>This investigation aims to delve into the unconditional uniform superconvergence behavior of the scalar auxiliary variable (SAV) Crank-Nicolson (CN) fully discrete scheme for the fourth-order singular perturbation Bi-flux diffusion equation with anisotropic constrained nonconforming rotated <span><math><msub><mi>Q</mi><mn>1</mn></msub></math></span> (<span><math><mrow><mi>C</mi><mi>N</mi><msubsup><mi>Q</mi><mrow><mn>1</mn></mrow><mrow><mi>r</mi><mi>o</mi><mi>t</mi></mrow></msubsup></mrow></math></span>) finite element method (FEM). Innovative high-accuracy estimates tailored specifically for the <span><math><mrow><mi>C</mi><mi>N</mi><msubsup><mi>Q</mi><mrow><mn>1</mn></mrow><mrow><mi>r</mi><mi>o</mi><mi>t</mi></mrow></msubsup></mrow></math></span> element on anisotropic meshes are established as the cornerstone for achieving superconvergent outcomes. Following this, a novel SAV-CN scheme is formulated, with the energy dissipation theoretically given to guarantee that the numerical solutions remain bounded in the broken <em>H</em><sup>1</sup>-norm. Then, by leveraging the aforementioned novel estimates and sidestepping the reliance on inverse inequalities, the unconditional uniform superclose result is rigorously derived, which transcends the constraints of the negative power of the perturbation parameter <span><math><mrow><mi>A</mi><mo>=</mo><msqrt><mrow><mi>β</mi><mo>(</mo><mn>1</mn><mo>−</mo><mi>β</mi><mo>)</mo><mi>R</mi></mrow></msqrt></mrow></math></span> and the limitations imposed by the ratio between the temporal step <em>τ</em> and spatial partition parameter <em>h</em>. Furthermore, an application of the economic interpolation post-processing technique enables the procurement of superconvergence estimates for the devised numerical schemes. Ultimately, the theoretical analysis is strengthened through the conduct of numerical experiments.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"207 ","pages":"Pages 41-59"},"PeriodicalIF":2.5,"publicationDate":"2026-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145996272","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

Unconditionally convergent fourth-order compact difference methods for Quasi-linear convection-diffusion equations 拟线性对流扩散方程的无条件收敛四阶紧致差分方法

IF 2.5 2区数学 Q1 MATHEMATICS, APPLIED

Computers & Mathematics with Applications

Pub Date : 2026-01-19 DOI: 10.1016/j.camwa.2026.01.018

Mengling Wu , Hongling Hu , Kejia Pan , Wei Wu

In this paper, we propose an unconditionally convergent implicit compact difference method for the quasi-linear convective diffusion equation. First, we use Taylor expansion combined with the truncation error residual correction method to discretize the first-order derivative in the temporal direction, while the first- and second-order derivatives in the spatial direction are approximated by the consistent fourth-order boundary schemes and the fourth-order difference method, respectively. We develop a two-level scheme with

O (τ^{2} + h^{4})

. Secondly, the convergence of the scheme in the L^∞ norm has been proven. To further improve computational efficiency, we apply Richardson extrapolation to improve the temporal accuracy of the scheme to fourth-order, i.e.,

O (τ^{4} + h^{4})

. Finally, we present numerical examples to verify the theoretical analysis and the efficiency of our proposed method.

本文给出了拟线性对流扩散方程的一种无条件收敛隐式紧差分法。首先，采用Taylor展开式结合截断误差残差校正方法对时间方向上的一阶导数进行离散化，而空间方向上的一阶导数和二阶导数分别采用一致的四阶边界格式和四阶差分法进行逼近。我们开发了一个具有O（τ2+h4）的两能级格式。其次，证明了该方案在L∞范数上的收敛性。为了进一步提高计算效率，我们采用Richardson外推将格式的时间精度提高到四阶，即O（τ4+h4）。最后，通过数值算例验证了理论分析和所提方法的有效性。

{"title":"Unconditionally convergent fourth-order compact difference methods for Quasi-linear convection-diffusion equations","authors":"Mengling Wu , Hongling Hu , Kejia Pan , Wei Wu","doi":"10.1016/j.camwa.2026.01.018","DOIUrl":"10.1016/j.camwa.2026.01.018","url":null,"abstract":"<div><div>In this paper, we propose an unconditionally convergent implicit compact difference method for the quasi-linear convective diffusion equation. First, we use Taylor expansion combined with the truncation error residual correction method to discretize the first-order derivative in the temporal direction, while the first- and second-order derivatives in the spatial direction are approximated by the consistent fourth-order boundary schemes and the fourth-order difference method, respectively. We develop a two-level scheme with <span><math><mrow><mi>O</mi><mrow><mo>(</mo><mrow><msup><mi>τ</mi><mn>2</mn></msup><mo>+</mo><msup><mi>h</mi><mn>4</mn></msup></mrow><mo>)</mo></mrow></mrow></math></span>. Secondly, the convergence of the scheme in the <em>L</em><sup>∞</sup> norm has been proven. To further improve computational efficiency, we apply Richardson extrapolation to improve the temporal accuracy of the scheme to fourth-order, i.e., <span><math><mrow><mi>O</mi><mrow><mo>(</mo><mrow><msup><mi>τ</mi><mn>4</mn></msup><mo>+</mo><msup><mi>h</mi><mn>4</mn></msup></mrow><mo>)</mo></mrow></mrow></math></span>. Finally, we present numerical examples to verify the theoretical analysis and the efficiency of our proposed method.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"205 ","pages":"Pages 268-281"},"PeriodicalIF":2.5,"publicationDate":"2026-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146001484","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 corrected Alikhanov scheme for a subdiffusion equation 亚扩散方程的修正Alikhanov格式

IF 2.5 2区数学 Q1 MATHEMATICS, APPLIED

Computers & Mathematics with Applications

Pub Date : 2026-01-19 DOI: 10.1016/j.camwa.2026.01.014

Chaobao Huang , Yujie Yu , Na An , Hu Chen

This paper considers the subdiffusion equation with a weakly singular solution. To achieve the optimal accuracy with a smaller grading parameter r than that of the standard Alikhanov scheme, it’s essential to correct the Alikhanov scheme and the approximation

v^{n - θ} \approx θ v^{n} + (1 - θ) v^{n - 1}

(0 < θ < 1) by investigating the corrected terms

β_{n}^{σ} (v^{1} - v^{0})

and

μ_{n}^{σ} (v^{1} - v^{0})

, respectively. After that, the truncation error of the corrected Alikhanov scheme and the corrected approximation of

v^{n - θ}

are given. By adopting the corrected Alikhanov scheme for the Caputo derivative and the corrected approximation for

v^{n - θ}

, we construct a fully discrete scheme for the subdiffusion equation, employing a standard finite element method in space. Furthermore, the stability analysis and the optimal convergent analysis for the proposed scheme are investigated. Finally, numerical experiments are conducted to verify the theoretical results.

研究一类具有弱奇异解的次扩散方程。为了获得比标准Alikhanov格式更小的分级参数r的最优精度，必须对Alikhanov格式和近似vn−θ≈θvn+(1−θ)vn−1（0 <； θ <； 1）分别通过研究校正项βnσ(v1−v0)和μnσ(v1−v0)进行校正。然后给出了修正后的Alikhanov格式的截断误差和修正后的vn−θ近似值。采用Caputo导数的修正Alikhanov格式和vn−θ的修正近似，利用空间标准有限元方法构造了亚扩散方程的完全离散格式。此外，还研究了该方案的稳定性分析和最优收敛性分析。最后通过数值实验对理论结果进行了验证。

{"title":"A corrected Alikhanov scheme for a subdiffusion equation","authors":"Chaobao Huang , Yujie Yu , Na An , Hu Chen","doi":"10.1016/j.camwa.2026.01.014","DOIUrl":"10.1016/j.camwa.2026.01.014","url":null,"abstract":"<div><div>This paper considers the subdiffusion equation with a weakly singular solution. To achieve the optimal accuracy with a smaller grading parameter <em>r</em> than that of the standard Alikhanov scheme, it’s essential to correct the Alikhanov scheme and the approximation <span><math><mrow><msup><mi>v</mi><mrow><mi>n</mi><mo>−</mo><mi>θ</mi></mrow></msup><mo>≈</mo><mi>θ</mi><msup><mi>v</mi><mi>n</mi></msup><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>−</mo><mi>θ</mi><mo>)</mo></mrow><msup><mi>v</mi><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow></msup></mrow></math></span> (0 < <em>θ</em> < 1) by investigating the corrected terms <span><math><mrow><msubsup><mi>β</mi><mrow><mi>n</mi></mrow><mi>σ</mi></msubsup><mrow><mo>(</mo><msup><mi>v</mi><mn>1</mn></msup><mo>−</mo><msup><mi>v</mi><mn>0</mn></msup><mo>)</mo></mrow></mrow></math></span> and <span><math><mrow><msubsup><mi>μ</mi><mi>n</mi><mi>σ</mi></msubsup><mrow><mo>(</mo><msup><mi>v</mi><mn>1</mn></msup><mo>−</mo><msup><mi>v</mi><mn>0</mn></msup><mo>)</mo></mrow></mrow></math></span>, respectively. After that, the truncation error of the corrected Alikhanov scheme and the corrected approximation of <span><math><msup><mi>v</mi><mrow><mi>n</mi><mo>−</mo><mi>θ</mi></mrow></msup></math></span> are given. By adopting the corrected Alikhanov scheme for the Caputo derivative and the corrected approximation for <span><math><msup><mi>v</mi><mrow><mi>n</mi><mo>−</mo><mi>θ</mi></mrow></msup></math></span>, we construct a fully discrete scheme for the subdiffusion equation, employing a standard finite element method in space. Furthermore, the stability analysis and the optimal convergent analysis for the proposed scheme are investigated. Finally, numerical experiments are conducted to verify the theoretical results.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"207 ","pages":"Pages 15-26"},"PeriodicalIF":2.5,"publicationDate":"2026-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145996271","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

Mixed finite element method for a system of hemivariational inequalities of Navier-Stokes equations coupled with the heat equation Navier-Stokes方程与热方程耦合的半变分不等式系统的混合有限元方法

IF 2.5 2区数学 Q1 MATHEMATICS, APPLIED

Computers & Mathematics with Applications

Pub Date : 2026-01-19 DOI: 10.1016/j.camwa.2026.01.013

Hailing Xuan , Wensi Wang , Xiaoliang Cheng

This paper explores a mixed finite element method for a system of hemivariational inequalities arising from the stationary Navier-Stokes equations coupled with the heat equation. The hemivariational inequalities system models the motion of a viscous incompressible fluid with thermal effects, where both the boundary conditions for the velocity and temperature fields involve the generalized Clarke gradient. We discuss the equivalence between different variational formulations and establish an existence theorem based on an alternative iterative approach and some results on hemivariational inequalities. The Navier-Stokes hemivariational inequalities system is solved using the mixed finite element method, and corresponding error estimates are derived. Numerical outcomes are presented based on the Uzawa method, which illustrates the optimal convergence rate as indicated by the error analysis.

本文探讨了求解由稳态Navier-Stokes方程与热方程耦合引起的半变分不等式系统的混合有限元方法。半变分不等式系统模拟了具有热效应的粘性不可压缩流体的运动，其中速度场和温度场的边界条件都涉及广义Clarke梯度。本文讨论了不同变分形式之间的等价性，并基于一种可选迭代方法和关于半变分不等式的一些结果，建立了一个存在性定理。采用混合有限元法求解了Navier-Stokes半变不等式系统，并给出了相应的误差估计。给出了基于Uzawa方法的数值结果，并通过误差分析说明了该方法的最优收敛速度。

引用次数: 0

A hybrid MR-WENO scheme with a simplified troubled-cell indicator for hyperbolic conservation laws 双曲型守恒律的带有简化故障单元指示器的MR-WENO混合方案

IF 2.5 2区数学 Q1 MATHEMATICS, APPLIED

Computers & Mathematics with Applications

Pub Date : 2026-01-19 DOI: 10.1016/j.camwa.2026.01.017

Jinming Zhang , Zhanjing Tao , Jun Zhu , Jianxian Qiu

In this paper, we propose a finite difference hybrid weighted essentially non-oscillatory (WENO) scheme to solve the hyperbolic conservation laws. A simplified troubled-cell indicator is designed for the hybrid scheme. The multi-resolution WENO (MR-WENO) reconstruction is applied in the troubled-cells, while a simple linear reconstruction is used in the remaining regions. The new hybrid scheme inherits the excellent characteristics of the original MR-WENO scheme [1], and can reduce the expensive computational cost of the WENO reconstruction. Compared to the previous hybrid WENO scheme [2] which used the high-degree polynomial in the troubled-cell indicator, the new scheme can reduce the percentage of the troubled-cells, leading to higher computational efficiency. Moreover, our scheme can effectively identify the troubled-cells, and has better resolution for certain problems than the previous scheme. Extensive numerical examples demonstrate the accuracy, efficiency and high resolution of the proposed method.

本文提出了一种有限差分混合加权本质非振荡（WENO）格式来求解双曲守恒律。针对混合方案，设计了一种简化的故障单元指示器。故障单元采用多分辨率WENO （MR-WENO）重建，其余区域采用简单的线性重建。新的混合方案继承了原MR-WENO方案[1]的优良特性，降低了WENO重构的计算成本。与以往在故障单元指标中使用高次多项式的混合WENO方案[2]相比，新方案可以减少故障单元的百分比，从而提高计算效率。此外，该方案能有效地识别出故障单元，对某些问题的解决效果优于原有方案。大量的数值算例证明了该方法的精度、效率和高分辨率。

引用次数: 0

Symmetric direct discontinuous Galerkin method for the biharmonic equation with non-homogeneous boundary condition 非齐次边界条件双调和方程的对称直接不连续Galerkin方法

IF 2.5 2区数学 Q1 MATHEMATICS, APPLIED

Computers & Mathematics with Applications

Pub Date : 2026-01-19 DOI: 10.1016/j.camwa.2026.01.010

Hongying Huang , Huanhuan Wang , Lin Zhang

The direct discontinuous Galerkin method is proposed and analyzed for the biharmonic equation with non-homogeneous boundary conditions. Numerical fluxes on all edges of arbitrary triangles are defined, which are concerned with four parameters (β₁, β₂, β₃, β₄), and then the weak variational form of the original problem is derived. The corresponding bilinear form is proved to be bounded and coercive for sufficiently large values of the parameters β₁ and β₂. The well-posedness of the variational problem in finite element space

V_{h}^{k}

and optimal error estimates of the approximate solution under L²-norm and energy norm are obtained. The accuracy and reliability of the proposed method were verified through numerical experiments. Although the theoretical analysis requires the solution u ∈ H⁴, numerical results indicate that the method remains effective for solutions with low regularity.

提出并分析了具有非齐次边界条件的双调和方程的直接不连续伽辽金方法。定义了任意三角形所有边上与四个参数（β1、β2、β3、β4）有关的数值通量，导出了原问题的弱变分形式。当参数β1和β2的值足够大时，证明了相应的双线性形式是有界的和强制的。得到了有限元空间Vhk中变分问题的适定性，以及在l2范数和能量范数下近似解的最优误差估计。通过数值实验验证了该方法的准确性和可靠性。虽然理论分析需要解u ∈ H4，但数值结果表明，该方法对于低正则性解仍然有效。

{"title":"Symmetric direct discontinuous Galerkin method for the biharmonic equation with non-homogeneous boundary condition","authors":"Hongying Huang , Huanhuan Wang , Lin Zhang","doi":"10.1016/j.camwa.2026.01.010","DOIUrl":"10.1016/j.camwa.2026.01.010","url":null,"abstract":"<div><div>The direct discontinuous Galerkin method is proposed and analyzed for the biharmonic equation with non-homogeneous boundary conditions. Numerical fluxes on all edges of arbitrary triangles are defined, which are concerned with four parameters (<em>β</em><sub>1</sub>, <em>β</em><sub>2</sub>, <em>β</em><sub>3</sub>, <em>β</em><sub>4</sub>), and then the weak variational form of the original problem is derived. The corresponding bilinear form is proved to be bounded and coercive for sufficiently large values of the parameters <em>β</em><sub>1</sub> and <em>β</em><sub>2</sub>. The well-posedness of the variational problem in finite element space <span><math><msubsup><mi>V</mi><mi>h</mi><mi>k</mi></msubsup></math></span> and optimal error estimates of the approximate solution under <em>L</em><sup>2</sup>-norm and energy norm are obtained. The accuracy and reliability of the proposed method were verified through numerical experiments. Although the theoretical analysis requires the solution <em>u</em> ∈ <em>H</em><sup>4</sup>, numerical results indicate that the method remains effective for solutions with low regularity.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"207 ","pages":"Pages 27-40"},"PeriodicalIF":2.5,"publicationDate":"2026-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145996273","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