Numerical Methods for Partial Differential Equations最新文献

{"title":"Unconditional H2$$ {H}^2 $$‐stability of the Euler implicit/explicit SAV‐based scheme for the 2D Navier–Stokes equations with smooth or nonsmooth initial data","authors":"Teng‐Yuan Chang, Ming‐Cheng Shiue","doi":"10.1002/num.23099","DOIUrl":"https://doi.org/10.1002/num.23099","url":null,"abstract":"In this article, we propose ‐unconditional stable schemes for solving time‐dependent incompressible Navier–Stokes equations with smooth or nonsmooth initial data, , . The ‐stability analysis is established by leveraging the scalar auxiliary variable (SAV) approach. When dealing with nonsmooth initial data, we utilize a limited number of iteration of the semi‐implicit scheme followed by the SAV scheme. The overall efficiency is greatly enhanced due to the minimal computational cost of the semi‐implicit scheme and the explicit treatment of the nonlinear term within the SAV approach. The proposed schemes investigate two types of scalar auxiliary variables: the energy‐based variable and the exponential‐based variable. Rigorous proofs of the ‐unconditional stability of both schemes have been provided. Notice that both proposed numerical schemes enjoy unconditional long time stability for smooth and nonsmooth initial data when . Numerical experiments have been conducted to demonstrate the theoretical results.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140313437","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A posteriori error estimate of the weak Galerkin finite element method solving the Stokes problems on polytopal meshes","authors":"Shipeng Xu","doi":"10.1002/num.23102","DOIUrl":"https://doi.org/10.1002/num.23102","url":null,"abstract":"In this paper, we propose an a posteriori error estimate of the weak Galerkin finite element method (WG-FEM) solving the Stokes problems with variable coefficients. Its error estimator, based on the property of Stokes' law conservation, Helmholtz decomposition and bubble functions, yields global upper bound and local lower bound for the approximation error of the WG-FEM. Error analysis is proved to be valid under the mesh assumptions of the WG-FEM and the way can be extended to other FEMs with the property of Stokes' law conservation, for example, discontinuous Galerkin (DG) FEMs. Finally, we verify the performance of error estimator by performing a few numerical examples.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140325209","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"An a posteriori error analysis for an augmented discontinuous Galerkin method applied to Stokes problem","authors":"Tomás P. Barrios, Rommel Bustinza","doi":"10.1002/num.23100","DOIUrl":"https://doi.org/10.1002/num.23100","url":null,"abstract":"This paper deals with the a posteriori error analysis for an augmented mixed discontinuous formulation for the stationary Stokes problem. By considering an appropriate auxiliary problem, we derive an a posteriori error estimator. We prove that this estimator is reliable and locally efficient, and consists of just five residual terms. Numerical experiments confirm the theoretical properties of the augmented discontinuous scheme as well as of the estimator. They also show the capability of the corresponding adaptive algorithm to localize the singularities and the large stress regions of the solution.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-03-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140196270","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A locally calculable P3‐pressure in a decoupled method for incompressible Stokes equations","authors":"Chunjae Park","doi":"10.1002/num.23101","DOIUrl":"https://doi.org/10.1002/num.23101","url":null,"abstract":"This article will suggest a new finite element method to find a ‐velocity and a ‐pressure solving incompressible Stokes equations at low cost. The method solves first the decoupled equation for a ‐velocity. Then, using the calculated velocity, a locally calculable ‐pressure will be defined component‐wisely. The resulting ‐pressure is analyzed to have the optimal order of convergence. Since the pressure is calculated by local computation only, the chief time cost of the new method is on solving the decoupled equation for the ‐velocity. Besides, the method overcomes the problem of singular vertices or corners.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-03-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140196276","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A new optimal error analysis of a mixed finite element method for advection–diffusion–reaction Brinkman flow","authors":"Huadong Gao, Wen Xie","doi":"10.1002/num.23097","DOIUrl":"https://doi.org/10.1002/num.23097","url":null,"abstract":"This article deals with the error analysis of a Galerkin‐mixed finite element methods for the advection–reaction–diffusion Brinkman flow in porous media. Numerical methods for incompressible miscible flow in porous media have been studied extensively in the last several decades. In practical applications, the lowest‐order Galerkin‐mixed method is the most popular one, where the linear Lagrange element is used for the concentration and the lowest‐order Raviart–Thomas element, the lowest‐order Nédélec edge element and piece‐wise constant discontinuous Galerkin element are used for the velocity, vorticity and pressure, respectively. The existing error estimate of this lowest‐order finite element method is only for all variables in spatial direction, which is not optimal for the concentration variable. This paper focuses on a new and optimal error estimate of a linearized backward Euler Galerkin‐mixed FEMs, where the second‐order accuracy for the concentration in spatial directions is established unconditionally. The key to our optimal error analysis is a new negative norm estimate for Nédélec edge element. Moreover, based on the computed numerical concentration, we propose a simple one‐step recovery technique to obtain a new numerical velocity, vorticity and pressure with second‐order accuracy. Numerical experiments are provided to confirm our theoretical analysis.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-03-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140153443","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Unconditional optimal first‐order error estimates of a full pressure segregation scheme for the magnetohydrodynamics equations","authors":"Yun‐Bo Yang, Yao‐Lin Jiang","doi":"10.1002/num.23098","DOIUrl":"https://doi.org/10.1002/num.23098","url":null,"abstract":"In this article, a first‐order linear fully discrete pressure segregation scheme is studied for the time‐dependent incompressible magnetohydrodynamics (MHD) equations in three‐dimensional bounded domain. Based on an incremental pressure projection method, this scheme allows us to decouple the MHD system into two sub‐problems at each time step, one is the velocity‐magnetic field system, the other is the pressure system. Firstly, a coupled linear elliptic system is solved for the velocity and the magnetic field. Next, a Poisson‐Neumann problem is treated for the pressure. We analyze the temporal error and the spatial error, respectively, and derive the temporal‐spatial error estimates of for the velocity and the magnetic field in the discrete space and for the pressure in the discrete space without imposing constraints on the mesh width and the time step size . Finally, some numerical results are presented to confirm the theoretical predictions and demonstrate the efficiency of the method.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-03-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140153447","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Adaptive finite element methods for scalar double‐well problem","authors":"Bingzhen Li, Dongjie Liu","doi":"10.1002/num.23096","DOIUrl":"https://doi.org/10.1002/num.23096","url":null,"abstract":"Some scalar double‐well problems eventually lead to a degenerate convex minimization problem with unique stress. We consider the adaptive conforming and nonconforming finite element methods for the scalar double‐well problem with the reliable a posteriori error analysis. A number of experiments confirm the effective decay rates of the methods.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-03-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140070170","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"On strong convergence of a fully discrete scheme for solving stochastic strongly damped wave equations","authors":"Chengqiang Xu, Yibo Wang, Wanrong Cao","doi":"10.1002/num.23094","DOIUrl":"https://doi.org/10.1002/num.23094","url":null,"abstract":"This article develops an efficient fully discrete scheme for a stochastic strongly damped wave equation (SSDWE) driven by an additive noise and presents its error estimates in the strong sense. We use the truncated spectral expansion of the noise to get an approximate equation and prove its regularity. Then we establish a spatio-temporal discretization of the approximate equation by a finite element method in space and an exponential trapezoidal scheme in time. We prove that the combination can derive higher strong convergence order in time than the use of the piecewise approximation of the noise and the exponential Euler scheme or the implicit Euler scheme in time. Particularly, the temporal strong convergence order of the fully discrete scheme reaches <mjx-container aria-label=\"5 divided by 4 minus epsilon\" ctxtmenu_counter=\"0\" ctxtmenu_oldtabindex=\"1\" jax=\"CHTML\" role=\"application\" sre-explorer- style=\"font-size: 103%; position: relative;\" tabindex=\"0\"><mjx-math aria-hidden=\"true\"><mjx-semantics><mjx-mrow data-semantic-children=\"5,4\" data-semantic-content=\"3\" data-semantic- data-semantic-role=\"subtraction\" data-semantic-speech=\"5 divided by 4 minus epsilon\" data-semantic-type=\"infixop\"><mjx-mrow data-semantic-children=\"0,2\" data-semantic-content=\"1\" data-semantic- data-semantic-parent=\"6\" data-semantic-role=\"division\" data-semantic-type=\"infixop\"><mjx-mn data-semantic-annotation=\"clearspeak:simple\" data-semantic-font=\"normal\" data-semantic- data-semantic-parent=\"5\" data-semantic-role=\"integer\" data-semantic-type=\"number\"><mjx-c></mjx-c></mjx-mn><mjx-mo data-semantic- data-semantic-operator=\"infixop,/\" data-semantic-parent=\"5\" data-semantic-role=\"division\" data-semantic-type=\"operator\" rspace=\"1\" space=\"1\"><mjx-c></mjx-c></mjx-mo><mjx-mn data-semantic-annotation=\"clearspeak:simple\" data-semantic-font=\"normal\" data-semantic- data-semantic-parent=\"5\" data-semantic-role=\"integer\" data-semantic-type=\"number\"><mjx-c></mjx-c></mjx-mn></mjx-mrow><mjx-mo data-semantic- data-semantic-operator=\"infixop,−\" data-semantic-parent=\"6\" data-semantic-role=\"subtraction\" data-semantic-type=\"operator\" rspace=\"1\" style=\"margin-left: 0.056em;\"><mjx-c></mjx-c></mjx-mo><mjx-mi data-semantic-annotation=\"clearspeak:simple\" data-semantic-font=\"italic\" data-semantic- data-semantic-parent=\"6\" data-semantic-role=\"greekletter\" data-semantic-type=\"identifier\"><mjx-c></mjx-c></mjx-mi></mjx-mrow></mjx-semantics></mjx-math><mjx-assistive-mml aria-hidden=\"true\" display=\"inline\" unselectable=\"on\"><math altimg=\"/cms/asset/246ca0b6-0dc3-49d5-80f8-c979b4835343/num23094-math-0001.png\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><semantics><mrow data-semantic-=\"\" data-semantic-children=\"5,4\" data-semantic-content=\"3\" data-semantic-role=\"subtraction\" data-semantic-speech=\"5 divided by 4 minus epsilon\" data-semantic-type=\"infixop\"><mrow data-semantic-=\"\" data-semantic-children=\"0,2\" data-semantic-content=\"1\" data-semantic-parent=\"6\" data-semantic-role=\"division\" data-semantic-type=\"infixo","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-02-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139948964","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Discrete null field equation methods for solving Laplace's equation: Boundary layer computations","authors":"Li-Ping Zhang, Zi‐Cai Li, Ming-Gong Lee, Hung‐Tsai Huang","doi":"10.1002/num.23092","DOIUrl":"https://doi.org/10.1002/num.23092","url":null,"abstract":"Consider Dirichlet problems of Laplace's equation in a bounded simply‐connected domain , and use the null field equation (NFE) of Green's representation formulation, where the source nodes are located on a pseudo‐boundary outside but not close to its boundary . Simple algorithms are proposed in this article by using the central rule for the NFE, and the normal derivatives of the solutions on the boundary can be easily obtained. These algorithms are called the discrete null field equation method (DNFEM) because the collocation equations are, indeed, the direct discrete form of the NFE. The bounds of the condition number are like those by the method of fundamental solutions (MFS) yielding the exponential growth as the number of unknowns increases. One trouble of the DNFEM is the near singularity of integrations for the solutions in boundary layers in Green's representation formulation. The traditional BEM also suffers from the same trouble. To deal with the near singularity, quadrature by expansions and the sinh transformation are often used. To handle this trouble, however, we develop two kinds of new techniques: (I) the interpolation techniques by Taylor's formulas with piecewise ‐degree polynomials and the Fourier series, and (II) the mini‐rules of integrals, such as the mini‐Simpson's and the mini‐Gaussian rules. Error analysis is made for technique I to achieve optimal convergence rates. Numerical experiments are carried out for disk domains to support the theoretical analysis made. The numerical performance of the DNFEM is excellent for disk domains to compete with the MFS. The errors with can be obtained by combination algorithms, which are satisfactory for most engineering problems. In summary, the new simple DNFEM is based on the NFE, which is different from the boundary element method (BEM). The theoretical basis in error and stability has been established in this article. One trouble in seeking the numerical solutions in boundary layers has been handled well; this is also an important contribution to the BEM. Besides, numerical experiments are encouraging. Hence the DNFEM is promising, and it may become a new boundary method for scientific/engineering computing.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-02-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139836928","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"New quadratic/serendipity finite volume element solutions on arbitrary triangular/quadrilateral meshes","authors":"Yanhui Zhou","doi":"10.1002/num.23093","DOIUrl":"https://doi.org/10.1002/num.23093","url":null,"abstract":"By postprocessing quadratic and eight‐node serendipity finite element solutions on arbitrary triangular and quadrilateral meshes, we obtain new quadratic/serendipity finite volume element solutions for solving anisotropic diffusion equations. The postprocessing procedure is implemented in each element independently, and we only need to solve a 6‐by‐6 (resp. 8‐by‐8) local linear algebraic system for triangular (resp. quadrilateral) element. The novelty of this paper is that, by designing some new quadratic dual meshes, and adding six/eight special constructed element‐wise bubble functions to quadratic/serendipity finite element solutions, we prove that the postprocessed solutions satisfy local conservation property on the new dual meshes. In particular, for any full anisotropic diffusion tensor, arbitrary triangular and quadrilateral meshes, we present a general framework to prove the existence and uniqueness of new quadratic/serendipity finite volume element solutions, which is better than some existing ones. That is, the existing theoretical results are improved, especially we extend the traditional rectangular assumption to arbitrary convex quadrilateral mesh. As a byproduct, we also prove that the new solutions converge to exact solution with optimal convergence rates under and norms on primal arbitrary triangular/quasi–parallelogram meshes. Finally, several numerical examples are carried out to validate the theoretical findings.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-02-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139777251","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}