首页 > 最新文献

Applied Numerical Mathematics最新文献

英文 中文
Error estimates of implicit-explicit compact BDF2 schemes for the pseudo parabolic equations with logarithmic nonlinearity
IF 2.2 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-01-23 DOI: 10.1016/j.apnum.2025.01.010
Qifeng Zhang , Haiyan Cao , Hongyu Qin
In this paper, two classes of linearized difference schemes are presented for the one and two-dimensional pseudo parabolic equations with logarithmic nonlinearity. These schemes are derived based on the implicit-explicit second-order backward differential formula (BDF2) for the temporal discretization and fourth-order compact/second-order difference schemes for spatial discretization. With the help of the truncation function method and regularization technique, error estimates for the logarithmic nonlinear term are handled rigorously. As a result, the convergence of the fully-discrete schemes is obtained based on the energy argument. Extensive numerical examples are presented to confirm the theoretical results.
{"title":"Error estimates of implicit-explicit compact BDF2 schemes for the pseudo parabolic equations with logarithmic nonlinearity","authors":"Qifeng Zhang ,&nbsp;Haiyan Cao ,&nbsp;Hongyu Qin","doi":"10.1016/j.apnum.2025.01.010","DOIUrl":"10.1016/j.apnum.2025.01.010","url":null,"abstract":"<div><div>In this paper, two classes of linearized difference schemes are presented for the one and two-dimensional pseudo parabolic equations with logarithmic nonlinearity. These schemes are derived based on the implicit-explicit second-order backward differential formula (BDF2) for the temporal discretization and fourth-order compact/second-order difference schemes for spatial discretization. With the help of the truncation function method and regularization technique, error estimates for the logarithmic nonlinear term are handled rigorously. As a result, the convergence of the fully-discrete schemes is obtained based on the energy argument. Extensive numerical examples are presented to confirm the theoretical results.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"212 ","pages":"Pages 135-154"},"PeriodicalIF":2.2,"publicationDate":"2025-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143155012","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
Energy stability and error estimate of the RKMK2e scheme for the extended Fisher–Kolmogorov equation
IF 2.2 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-01-23 DOI: 10.1016/j.apnum.2025.01.014
Haifeng Wang, Yan Wang, Hong Zhang, Songhe Song
In this paper, we develop a second-order accurate scheme for the extended Fisher–Kolmogorov (EFK) equation and investigate its global-in-time energy stability and convergence. The proposed scheme uses the Fourier spectral collocation method in space and the stabilization Runge–Kutta–Munthe–Kaas-2e (RKMK2e) method for temporal approximation. To demonstrate the global-in-time energy stability of the proposed scheme, we first verify that, under the assumption that all numerical solutions are uniformly bounded, the scheme is energy stable when using a sufficiently large stabilization parameter. Then, to establish the uniform-in-time boundedness of the numerical solutions, we fully utilize the nonlinear operator estimates and discrete Sobolev embedding in each stage of the scheme. Moreover, we conduct an optimal rate convergence analysis with a sufficient regularity assumption for the exact solution. Several numerical examples are presented to validate the accuracy, computational efficiency, and energy stability of the proposed scheme.
{"title":"Energy stability and error estimate of the RKMK2e scheme for the extended Fisher–Kolmogorov equation","authors":"Haifeng Wang,&nbsp;Yan Wang,&nbsp;Hong Zhang,&nbsp;Songhe Song","doi":"10.1016/j.apnum.2025.01.014","DOIUrl":"10.1016/j.apnum.2025.01.014","url":null,"abstract":"<div><div>In this paper, we develop a second-order accurate scheme for the extended Fisher–Kolmogorov (EFK) equation and investigate its global-in-time energy stability and convergence. The proposed scheme uses the Fourier spectral collocation method in space and the stabilization Runge–Kutta–Munthe–Kaas-2e (RKMK2e) method for temporal approximation. To demonstrate the global-in-time energy stability of the proposed scheme, we first verify that, under the assumption that all numerical solutions are uniformly bounded, the scheme is energy stable when using a sufficiently large stabilization parameter. Then, to establish the uniform-in-time boundedness of the numerical solutions, we fully utilize the nonlinear operator estimates and discrete Sobolev embedding in each stage of the scheme. Moreover, we conduct an optimal rate convergence analysis with a sufficient regularity assumption for the exact solution. Several numerical examples are presented to validate the accuracy, computational efficiency, and energy stability of the proposed scheme.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"212 ","pages":"Pages 60-76"},"PeriodicalIF":2.2,"publicationDate":"2025-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143155010","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
Refined variational iteration method through optimal linear operator and Lagrange multiplier: Local/global convergence
IF 2.2 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-01-23 DOI: 10.1016/j.apnum.2025.01.013
Tapas Roy, Dilip Kumar Maiti
Our focus is on refining the Variational Iteration Method (VIM), which will henceforth be called R-VIM. We begin by introducing a general linear operator and consequently the Lagrange's multiplier; then, minimizing the residual error, we subsequently find the optimal linear operator as well as the best fitting Lagrange's multiplier for a given differential equation. By means of Banach Fixed Point Theorem, the necessary conditions for the convergence of the solutions of VIM, Optimal VIM, and R-VIM are established involving the Lagrange's multiplier; it has been demonstrated (and subsequently numerically certified) that solutions of VIM and Optimal VIM locally converge. While R-VIM solutions, with the aid of necessary and sufficient criteria, have been shown to converge globally for appropriate values of the unknown parameters involved in the general linear operator. We make sure that the convergence of the VIM solution depends critically on the optimal selection of the linear operator and, consequently, the best suited Lagrange's multiplier.
Our bird's eye view would be to differentiate between the contributions of the homotopy equation in homotopy-based approaches and Lagrange's multiplier-based correction functional of the VIM for the fastest convergence. Furthermore, in order to achieve fast convergence, we additionally examine if the optimal VIM's regulating parameter offers any further benefit above proposed R-VIM. We aim to apply R-VIM to problems with multiple solutions to ensure that all solutions are accurately identified. Our proposed method has been proven superior and more widely applicable through both theoretical and numerical illustrations. However, it also has its limits. Hence, we suggest the implementation of the multistep R-VIM, referred to as MR-VIM, as a means to solve chaotic dynamical systems.
{"title":"Refined variational iteration method through optimal linear operator and Lagrange multiplier: Local/global convergence","authors":"Tapas Roy,&nbsp;Dilip Kumar Maiti","doi":"10.1016/j.apnum.2025.01.013","DOIUrl":"10.1016/j.apnum.2025.01.013","url":null,"abstract":"<div><div>Our focus is on refining the Variational Iteration Method (VIM), which will henceforth be called R-VIM. We begin by introducing a general linear operator and consequently the Lagrange's multiplier; then, minimizing the residual error, we subsequently find the optimal linear operator as well as the best fitting Lagrange's multiplier for a given differential equation. By means of Banach Fixed Point Theorem, the necessary conditions for the convergence of the solutions of VIM, Optimal VIM, and R-VIM are established involving the Lagrange's multiplier; it has been demonstrated (and subsequently numerically certified) that solutions of VIM and Optimal VIM locally converge. While R-VIM solutions, with the aid of necessary and sufficient criteria, have been shown to converge globally for appropriate values of the unknown parameters involved in the general linear operator. We make sure that the convergence of the VIM solution depends critically on the optimal selection of the linear operator and, consequently, the best suited Lagrange's multiplier.</div><div>Our bird's eye view would be to differentiate between the contributions of the homotopy equation in homotopy-based approaches and Lagrange's multiplier-based correction functional of the VIM for the fastest convergence. Furthermore, in order to achieve fast convergence, we additionally examine if the optimal VIM's regulating parameter offers any further benefit above proposed R-VIM. We aim to apply R-VIM to problems with multiple solutions to ensure that all solutions are accurately identified. Our proposed method has been proven superior and more widely applicable through both theoretical and numerical illustrations. However, it also has its limits. Hence, we suggest the implementation of the multistep R-VIM, referred to as MR-VIM, as a means to solve chaotic dynamical systems.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"212 ","pages":"Pages 1-28"},"PeriodicalIF":2.2,"publicationDate":"2025-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143154831","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
An unconditionally stable second-order scheme for Maxwell's equations in the Cole–Cole dispersive medium
IF 2.2 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-01-22 DOI: 10.1016/j.apnum.2025.01.009
Jingjing Xiao , Desong Kong
Coupling the averaged L1 scheme and the Crank–Nicolson scheme for temporal derivatives, we study Maxwell's equations in the Cole–Cole dispersive medium. A rigorous analysis is carried out to show that the proposed scheme is unconditionally stable and has a second-order convergence in time for sufficiently smooth solutions. We avoid using the mathematical induction method, which simplifies the analysis than the existing schemes. A fully discrete scheme with a finite difference method at Yee's grid is proposed. Numerical experiments are carefully designed to illustrate our theoretical analysis.
{"title":"An unconditionally stable second-order scheme for Maxwell's equations in the Cole–Cole dispersive medium","authors":"Jingjing Xiao ,&nbsp;Desong Kong","doi":"10.1016/j.apnum.2025.01.009","DOIUrl":"10.1016/j.apnum.2025.01.009","url":null,"abstract":"<div><div>Coupling the averaged L1 scheme and the Crank–Nicolson scheme for temporal derivatives, we study Maxwell's equations in the Cole–Cole dispersive medium. A rigorous analysis is carried out to show that the proposed scheme is unconditionally stable and has a second-order convergence in time for sufficiently smooth solutions. We avoid using the mathematical induction method, which simplifies the analysis than the existing schemes. A fully discrete scheme with a finite difference method at Yee's grid is proposed. Numerical experiments are carefully designed to illustrate our theoretical analysis.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"211 ","pages":"Pages 211-227"},"PeriodicalIF":2.2,"publicationDate":"2025-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143141931","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
Optimal L2 error estimates of two structure-preserving finite element methods for Schrödinger-Boussinesq equations
IF 2.2 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-01-22 DOI: 10.1016/j.apnum.2025.01.012
Houchao Zhang, Junjun Wang, Xueran Gong
This paper is concerned with the construction and analysis of two structure-preserving finite element approximation schemes for solving one and two dimensional coupled Schrödinger-Boussinesq (SBq) equations. Firstly, two finite element approximation schemes are developed and the total mass and energy preserving properties of the proposed schemes are demonstrated in the discrete sense. Secondly, by use of the innovative cut-off function method, an auxiliary fully-discrete system is established, which leads to the unique solvability with the Brouwer's fixed-pointed theorem and the optimal L2 error estimates for the proposed nonlinear scheme. Finally, linearized iterative algorithms are showed rigorously and extensive numerical results are presented to check the theoretical analysis of the structure-preserving numerical methods.
{"title":"Optimal L2 error estimates of two structure-preserving finite element methods for Schrödinger-Boussinesq equations","authors":"Houchao Zhang,&nbsp;Junjun Wang,&nbsp;Xueran Gong","doi":"10.1016/j.apnum.2025.01.012","DOIUrl":"10.1016/j.apnum.2025.01.012","url":null,"abstract":"<div><div>This paper is concerned with the construction and analysis of two structure-preserving finite element approximation schemes for solving one and two dimensional coupled Schrödinger-Boussinesq (SBq) equations. Firstly, two finite element approximation schemes are developed and the total mass and energy preserving properties of the proposed schemes are demonstrated in the discrete sense. Secondly, by use of the innovative cut-off function method, an auxiliary fully-discrete system is established, which leads to the unique solvability with the Brouwer's fixed-pointed theorem and the optimal <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> error estimates for the proposed nonlinear scheme. Finally, linearized iterative algorithms are showed rigorously and extensive numerical results are presented to check the theoretical analysis of the structure-preserving numerical methods.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"211 ","pages":"Pages 193-210"},"PeriodicalIF":2.2,"publicationDate":"2025-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143141887","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
Gradient descent-based parameter-free methods for solving coupled matrix equations and studying an application in dynamical systems
IF 2.2 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-01-22 DOI: 10.1016/j.apnum.2025.01.011
Akbar Shirilord, Mehdi Dehghan
This paper explores advanced gradient descent-based parameter-free methods for solving coupled matrix equations and examines their applications in dynamical systems. We focus on the coupled matrix equations{AX+YB=C,DX+YE=F, where A,DRm×m,B,ERn×n,C,FRm×n are given matrices, and X,YRm×n are the unknown matrices to be determined. We propose a novel gradient descent-based approach with parameter-free, enhancing convergence through an accelerated technique related to momentum methods. A comprehensive analysis of the convergence and characteristics of these methods is provided. Our convergence analysis demonstrates that if the spectrum of a block matrix constructed from the matrices A, B, D, and E is confined within a horizontal ellipse in the complex plane, centered at (0,0) with a major axis length of 3 and a minor axis length of 1, then the accelerated momentum method will converge to the exact solution of the discussed model. The numerical results indicate that proposed methods significantly improve efficiency, showing faster convergence and reduced computational time compared to traditional approaches. Additionally, we apply these methods to linear dynamic systems, demonstrating their effectiveness in real-world scenarios.
{"title":"Gradient descent-based parameter-free methods for solving coupled matrix equations and studying an application in dynamical systems","authors":"Akbar Shirilord,&nbsp;Mehdi Dehghan","doi":"10.1016/j.apnum.2025.01.011","DOIUrl":"10.1016/j.apnum.2025.01.011","url":null,"abstract":"<div><div>This paper explores advanced gradient descent-based parameter-free methods for solving coupled matrix equations and examines their applications in dynamical systems. We focus on the coupled matrix equations<span><span><span><math><mrow><mo>{</mo><mtable><mtr><mtd><mi>A</mi><mi>X</mi><mo>+</mo><mi>Y</mi><mi>B</mi><mo>=</mo><mi>C</mi><mo>,</mo></mtd></mtr><mtr><mtd><mi>D</mi><mi>X</mi><mo>+</mo><mi>Y</mi><mi>E</mi><mo>=</mo><mi>F</mi><mo>,</mo></mtd></mtr></mtable></mrow></math></span></span></span> where <span><math><mi>A</mi><mo>,</mo><mi>D</mi><mo>∈</mo><msup><mrow><mi>R</mi></mrow><mrow><mi>m</mi><mo>×</mo><mi>m</mi></mrow></msup><mo>,</mo><mi>B</mi><mo>,</mo><mi>E</mi><mo>∈</mo><msup><mrow><mi>R</mi></mrow><mrow><mi>n</mi><mo>×</mo><mi>n</mi></mrow></msup><mo>,</mo><mi>C</mi><mo>,</mo><mi>F</mi><mo>∈</mo><msup><mrow><mi>R</mi></mrow><mrow><mi>m</mi><mo>×</mo><mi>n</mi></mrow></msup></math></span> are given matrices, and <span><math><mi>X</mi><mo>,</mo><mi>Y</mi><mo>∈</mo><msup><mrow><mi>R</mi></mrow><mrow><mi>m</mi><mo>×</mo><mi>n</mi></mrow></msup></math></span> are the unknown matrices to be determined. We propose a novel gradient descent-based approach with parameter-free, enhancing convergence through an accelerated technique related to momentum methods. A comprehensive analysis of the convergence and characteristics of these methods is provided. Our convergence analysis demonstrates that if the spectrum of a block matrix constructed from the matrices <em>A</em>, <em>B</em>, <em>D</em>, and <em>E</em> is confined within a horizontal ellipse in the complex plane, centered at <span><math><mo>(</mo><mn>0</mn><mo>,</mo><mn>0</mn><mo>)</mo></math></span> with a major axis length of 3 and a minor axis length of 1, then the accelerated momentum method will converge to the exact solution of the discussed model. The numerical results indicate that proposed methods significantly improve efficiency, showing faster convergence and reduced computational time compared to traditional approaches. Additionally, we apply these methods to linear dynamic systems, demonstrating their effectiveness in real-world scenarios.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"212 ","pages":"Pages 29-59"},"PeriodicalIF":2.2,"publicationDate":"2025-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143155008","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A posteriori error analysis of a mixed finite element method for the stationary convective Brinkman–Forchheimer problem
IF 2.2 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-01-17 DOI: 10.1016/j.apnum.2025.01.007
Sergio Caucao , Gabriel N. Gatica , Luis F. Gatica
We consider a Banach spaces-based mixed variational formulation that has been recently proposed for the nonlinear problem given by the stationary convective Brinkman–Forchheimer equations, and develop a reliable and efficient residual-based a posteriori error estimator for the 2D and 3D versions of the associated mixed finite element scheme. For the reliability analysis, we utilize the global inf-sup condition of the problem, combined with appropriate small data assumptions, a stable Helmholtz decomposition in nonstandard Banach spaces, and the local approximation properties of the Raviart–Thomas and Clément interpolants. In turn, inverse inequalities, the localization technique based on bubble functions in local Lp-spaces, and known results from previous works, are the main tools yielding the efficiency estimate. Finally, several numerical results confirming the theoretical properties of the estimator and illustrating the performance of the associated adaptive algorithm are reported. In particular, the case of flow through a 2D porous medium with fracture networks is considered.
{"title":"A posteriori error analysis of a mixed finite element method for the stationary convective Brinkman–Forchheimer problem","authors":"Sergio Caucao ,&nbsp;Gabriel N. Gatica ,&nbsp;Luis F. Gatica","doi":"10.1016/j.apnum.2025.01.007","DOIUrl":"10.1016/j.apnum.2025.01.007","url":null,"abstract":"<div><div>We consider a Banach spaces-based mixed variational formulation that has been recently proposed for the nonlinear problem given by the stationary convective Brinkman–Forchheimer equations, and develop a reliable and efficient residual-based <em>a posteriori</em> error estimator for the 2D and 3D versions of the associated mixed finite element scheme. For the reliability analysis, we utilize the global inf-sup condition of the problem, combined with appropriate small data assumptions, a stable Helmholtz decomposition in nonstandard Banach spaces, and the local approximation properties of the Raviart–Thomas and Clément interpolants. In turn, inverse inequalities, the localization technique based on bubble functions in local <span><math><msup><mrow><mi>L</mi></mrow><mrow><mi>p</mi></mrow></msup></math></span>-spaces, and known results from previous works, are the main tools yielding the efficiency estimate. Finally, several numerical results confirming the theoretical properties of the estimator and illustrating the performance of the associated adaptive algorithm are reported. In particular, the case of flow through a 2D porous medium with fracture networks is considered.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"211 ","pages":"Pages 158-178"},"PeriodicalIF":2.2,"publicationDate":"2025-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143141885","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 robust mapping spectral method for elastic equations in curved fan-shaped domains
IF 2.2 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-01-16 DOI: 10.1016/j.apnum.2025.01.008
Yuling Guo , Zhongqing Wang , Chao Zhang
In this paper, we introduce a robust mapping Legendre spectral-Galerkin method for addressing elastic problems in simply-connected, fan-shaped domains with curved boundaries. By employing a polar coordinate transformation, we map the fan-shaped domain onto a rectangle, which transforms the original elastic equation into a variable coefficient equation. We then develop a Legendre spectral-Galerkin scheme for this variable coefficient equation. Additionally, we demonstrate the optimal convergence of the numerical solution in the H1-norm as the Lamé coefficient λ remains bounded. Numerical examples illustrate the high accuracy and robustness of our method, even as λ approaches infinity.
{"title":"A robust mapping spectral method for elastic equations in curved fan-shaped domains","authors":"Yuling Guo ,&nbsp;Zhongqing Wang ,&nbsp;Chao Zhang","doi":"10.1016/j.apnum.2025.01.008","DOIUrl":"10.1016/j.apnum.2025.01.008","url":null,"abstract":"<div><div>In this paper, we introduce a robust mapping Legendre spectral-Galerkin method for addressing elastic problems in simply-connected, fan-shaped domains with curved boundaries. By employing a polar coordinate transformation, we map the fan-shaped domain onto a rectangle, which transforms the original elastic equation into a variable coefficient equation. We then develop a Legendre spectral-Galerkin scheme for this variable coefficient equation. Additionally, we demonstrate the optimal convergence of the numerical solution in the <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span>-norm as the Lamé coefficient <em>λ</em> remains bounded. Numerical examples illustrate the high accuracy and robustness of our method, even as <em>λ</em> approaches infinity.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"211 ","pages":"Pages 144-157"},"PeriodicalIF":2.2,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143141888","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
Variable-time-step weighted IMEX FEMs for nonlinear evolution equations
IF 2.2 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-01-15 DOI: 10.1016/j.apnum.2025.01.005
Meng Li , Dan Wang , Junjun Wang , Xiaolong Zhao
In this paper, the variable-time-step weighted implicit-explicit (IMEX) finite element methods (FEMs) are developed for some types of nonlinear real- or complex-valued evolution equations. Extensive research is conducted on the discrete orthogonal convolution (DOC) kernels and the discrete complementary convolution (DCC) kernels of the variable-time-step weighted IMEX scheme, elucidating their crucial properties in both real- and complex-valued scenarios. We prove that the scheme exhibits optimal convergence without any restrictions on the time-space step ratio. At last, several numerical examples are provided to demonstrate our theoretical results. With the weighted parameter θ=1, the scheme in this work can degenerate into a special case: variable-time-step two-step backward differentiation formula (BDF2) scheme, and the convergence analysis in this special case was introduced in Liao et al. (2020) [15] and Liao et al. (2021) [29].
{"title":"Variable-time-step weighted IMEX FEMs for nonlinear evolution equations","authors":"Meng Li ,&nbsp;Dan Wang ,&nbsp;Junjun Wang ,&nbsp;Xiaolong Zhao","doi":"10.1016/j.apnum.2025.01.005","DOIUrl":"10.1016/j.apnum.2025.01.005","url":null,"abstract":"<div><div>In this paper, the variable-time-step weighted implicit-explicit (IMEX) finite element methods (FEMs) are developed for some types of nonlinear real- or complex-valued evolution equations. Extensive research is conducted on the discrete orthogonal convolution (DOC) kernels and the discrete complementary convolution (DCC) kernels of the variable-time-step weighted IMEX scheme, elucidating their crucial properties in both real- and complex-valued scenarios. We prove that the scheme exhibits optimal convergence without any restrictions on the time-space step ratio. At last, several numerical examples are provided to demonstrate our theoretical results. With the weighted parameter <span><math><mi>θ</mi><mo>=</mo><mn>1</mn></math></span>, the scheme in this work can degenerate into a special case: variable-time-step two-step backward differentiation formula (BDF2) scheme, and the convergence analysis in this special case was introduced in Liao et al. (2020) <span><span>[15]</span></span> and Liao et al. (2021) <span><span>[29]</span></span>.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"211 ","pages":"Pages 123-143"},"PeriodicalIF":2.2,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143141889","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
Analysis of average bound preserving time-implicit discretizations for convection-diffusion-reaction equation
IF 2.2 2区 数学 Q1 MATHEMATICS, APPLIED Pub Date : 2025-01-15 DOI: 10.1016/j.apnum.2025.01.006
Fengna Yan , Yinhua Xia
We propose a high-order average bound-preserving limiter for implicit backward differentiation formula (BDF) and local discontinuous Galerkin (LDG) discretizations applied to convection-diffusion-reaction equations. Our approach first imposes cell average bounds of the numerical solution using the Karush-Kuhn-Tucker (KKT) limiter and then enforces pointwise bounds with an explicit bound-preserving limiter. This method reduces the number of constraints compared to using only the KKT system to directly ensure pointwise bounds, resulting in a relatively small system of nonlinear equations to solve at each time step. We prove the unique solvability of the proposed average bound-preserving BDF-LDG discretizations. Furthermore, we establish the stability and optimal error estimates for the second-order average bound-preserving BDF2-LDG discretization. The unique solvability and stability are derived by transforming the KKT-limited cell average bounds-preserving LDG discretizations into a variational inequality. The error estimates are derived using the cell average bounds-preserving inequality constraints. Numerical results are presented to validate the accuracy and effectiveness of the proposed method in preserving the bounds.
{"title":"Analysis of average bound preserving time-implicit discretizations for convection-diffusion-reaction equation","authors":"Fengna Yan ,&nbsp;Yinhua Xia","doi":"10.1016/j.apnum.2025.01.006","DOIUrl":"10.1016/j.apnum.2025.01.006","url":null,"abstract":"<div><div>We propose a high-order average bound-preserving limiter for implicit backward differentiation formula (BDF) and local discontinuous Galerkin (LDG) discretizations applied to convection-diffusion-reaction equations. Our approach first imposes cell average bounds of the numerical solution using the Karush-Kuhn-Tucker (KKT) limiter and then enforces pointwise bounds with an explicit bound-preserving limiter. This method reduces the number of constraints compared to using only the KKT system to directly ensure pointwise bounds, resulting in a relatively small system of nonlinear equations to solve at each time step. We prove the unique solvability of the proposed average bound-preserving BDF-LDG discretizations. Furthermore, we establish the stability and optimal error estimates for the second-order average bound-preserving BDF2-LDG discretization. The unique solvability and stability are derived by transforming the KKT-limited cell average bounds-preserving LDG discretizations into a variational inequality. The error estimates are derived using the cell average bounds-preserving inequality constraints. Numerical results are presented to validate the accuracy and effectiveness of the proposed method in preserving the bounds.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"211 ","pages":"Pages 103-122"},"PeriodicalIF":2.2,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143141891","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
期刊
Applied Numerical Mathematics
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1