Numerical Methods for Partial Differential Equations最新文献

筛选
英文 中文
Entropy stable discontinuous Galerkin methods for the shallow water equations with subcell positivity preservation 浅水方程的熵稳定非连续伽勒金方法与子单元实在性保持
IF 3.9 3区 数学
Numerical Methods for Partial Differential Equations Pub Date : 2024-07-25 DOI: 10.1002/num.23129
Xinhui Wu, Nathaniel Trask, Jesse Chan
{"title":"Entropy stable discontinuous Galerkin methods for the shallow water equations with subcell positivity preservation","authors":"Xinhui Wu, Nathaniel Trask, Jesse Chan","doi":"10.1002/num.23129","DOIUrl":"https://doi.org/10.1002/num.23129","url":null,"abstract":"High order schemes are known to be unstable in the presence of shock discontinuities or under‐resolved solution features, and have traditionally required additional filtering, limiting, or artificial viscosity to avoid solution blow up. Entropy stable schemes address this instability by ensuring that physically relevant solutions satisfy a semi‐discrete entropy inequality independently of discretization parameters. However, additional measures must be taken to ensure that solutions satisfy physical constraints such as positivity. In this work, we present a high order entropy stable discontinuous Galerkin (ESDG) method for the nonlinear shallow water equations (SWE) on two‐dimensional (2D) triangular meshes which preserves the positivity of the water heights. The scheme combines a low order positivity preserving method with a high order entropy stable method using convex limiting. This method is entropy stable and well‐balanced for fitted meshes with continuous bathymetry profiles.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":"67 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-07-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141777545","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}
引用次数: 0
Convergence and stability analysis of energy stable and bound‐preserving numerical schemes for binary fluid‐surfactant phase‐field equations 二元流体-表面活性相场方程的能量稳定和保界数值方案的收敛性和稳定性分析
IF 3.9 3区 数学
Numerical Methods for Partial Differential Equations Pub Date : 2024-07-01 DOI: 10.1002/num.23125
Jiayi Duan, Xiao Li, Zhonghua Qiao
{"title":"Convergence and stability analysis of energy stable and bound‐preserving numerical schemes for binary fluid‐surfactant phase‐field equations","authors":"Jiayi Duan, Xiao Li, Zhonghua Qiao","doi":"10.1002/num.23125","DOIUrl":"https://doi.org/10.1002/num.23125","url":null,"abstract":"In this article, we develop stable and efficient numerical schemes for a binary fluid‐surfactant phase‐field model which consists of two Cahn–Hilliard type equations with respect to the free energy containing a Ginzburg–Landau double‐well potential, a logarithmic Flory–Huggins potential and a nonlinear coupling entropy. The numerical schemes, which are decoupled and linear, are established by the central difference spatial approximation in combination with the first‐ and second‐order exponential time differencing methods based on the convex splitting of the free energy. For the sake of the linearity of the schemes, the nonlinear terms, especially the logarithmic term, are approximated explicitly, which requires the bound preservation of the numerical solution to make the algorithm robust. We conduct the convergence analysis and prove the bound‐preserving property in details for both first‐ and second‐order schemes, where the high‐order consistency analysis is applied to the first‐order case. In addition, the energy stability is also obtained by the nature of the convex splitting. Numerical experiments are performed to verify the accuracy and stability of the schemes and simulate the dynamics of phase separation and surfactant adsorption.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":"67 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141551400","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}
引用次数: 0
Semi‐implicit method of high‐index saddle dynamics and application to construct solution landscape 高指数鞍动力学的半隐式方法及其在构建解景观中的应用
IF 3.9 3区 数学
Numerical Methods for Partial Differential Equations Pub Date : 2024-06-29 DOI: 10.1002/num.23123
Yue Luo, Lei Zhang, Pingwen Zhang, Zhiyi Zhang, Xiangcheng Zheng
{"title":"Semi‐implicit method of high‐index saddle dynamics and application to construct solution landscape","authors":"Yue Luo, Lei Zhang, Pingwen Zhang, Zhiyi Zhang, Xiangcheng Zheng","doi":"10.1002/num.23123","DOIUrl":"https://doi.org/10.1002/num.23123","url":null,"abstract":"We analyze the semi‐implicit scheme of high‐index saddle dynamics, which provides a powerful numerical method for finding the any‐index saddle points and constructing the solution landscape. Compared with the explicit schemes of saddle dynamics, the semi‐implicit discretization relaxes the step size and accelerates the convergence, but the corresponding numerical analysis encounters new difficulties compared to the explicit scheme. Specifically, the orthonormal property of the eigenvectors at each time step could not be fully employed due to the semi‐implicit treatment, and computations of the eigenvectors are coupled with the orthonormalization procedure, which further complicates the numerical analysis. We address these issues to prove error estimates of the semi‐implicit scheme via, for example, technical splittings and multi‐variable circulating induction procedure. We further analyze the convergence rate of the generalized minimum residual solver for solving the semi‐implicit system. Extensive numerical experiments are carried out to substantiate the efficiency and accuracy of the semi‐implicit scheme in constructing solution landscapes of complex systems.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":"14 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-06-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141551401","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}
引用次数: 0
Analysis of a class of spectral volume methods for linear scalar hyperbolic conservation laws 线性标量双曲守恒定律的一类谱量法分析
IF 3.9 3区 数学
Numerical Methods for Partial Differential Equations Pub Date : 2024-06-28 DOI: 10.1002/num.23126
Jianfang Lu, Yan Jiang, Chi‐Wang Shu, Mengping Zhang
{"title":"Analysis of a class of spectral volume methods for linear scalar hyperbolic conservation laws","authors":"Jianfang Lu, Yan Jiang, Chi‐Wang Shu, Mengping Zhang","doi":"10.1002/num.23126","DOIUrl":"https://doi.org/10.1002/num.23126","url":null,"abstract":"In this article, we study the spectral volume (SV) methods for scalar hyperbolic conservation laws with a class of subdivision points under the Petrov–Galerkin framework. Due to the strong connection between the DG method and the SV method with the appropriate choice of the subdivision points, it is natural to analyze the SV method in the Galerkin form and derive the analogous theoretical results as in the DG method. This article considers a class of SV methods, whose subdivision points are the zeros of a specific polynomial with a parameter in it. Properties of the piecewise constant functions under this subdivision, including the orthogonality between the trial solution space and test function space, are provided. With the aid of these properties, we are able to derive the energy stability, optimal a priori error estimates of SV methods with arbitrary high order accuracy. We also study the superconvergence of the numerical solution with the correction function technique, and show the order of superconvergence would be different with different choices of the subdivision points. In the numerical experiments, by choosing different parameters in the SV method, the theoretical findings are confirmed by the numerical results.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":"44 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141525514","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}
引用次数: 0
Minimizers for the de Gennes–Cahn–Hilliard energy under strong anchoring conditions 强锚定条件下的 de Gennes-Cahn-Hilliard 能量最小值
IF 3.9 3区 数学
Numerical Methods for Partial Differential Equations Pub Date : 2024-06-27 DOI: 10.1002/num.23127
Shibin Dai, Abba Ramadan
{"title":"Minimizers for the de Gennes–Cahn–Hilliard energy under strong anchoring conditions","authors":"Shibin Dai, Abba Ramadan","doi":"10.1002/num.23127","DOIUrl":"https://doi.org/10.1002/num.23127","url":null,"abstract":"In this article, we use the Nehari manifold and the eigenvalue problem for the negative Laplacian with Dirichlet boundary condition to analytically study the minimizers for the de Gennes–Cahn–Hilliard energy with quartic double‐well potential and Dirichlet boundary condition on the bounded domain. Our analysis reveals a bifurcation phenomenon determined by the boundary value and a bifurcation parameter that describes the thickness of the transition layer that segregates the binary mixture's two phases. Specifically, when the boundary value aligns precisely with the average of the pure phases, and the bifurcation parameter surpasses or equals a critical threshold, the minimizer assumes a unique form, representing the homogeneous state. Conversely, when the bifurcation parameter falls below this critical value, two symmetric minimizers emerge. Should the boundary value be larger or smaller from the average of the pure phases, symmetry breaks, resulting in a unique minimizer. Furthermore, we derive bounds of these minimizers, incorporating boundary conditions and features of the de Gennes–Cahn–Hilliard energy.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":"14 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-06-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141525463","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}
引用次数: 0
A linear finite difference scheme with error analysis designed to preserve the structure of the 2D Boussinesq paradigm equation 带误差分析的线性有限差分方案,旨在保留二维布森斯克范式方程的结构
IF 3.9 3区 数学
Numerical Methods for Partial Differential Equations Pub Date : 2024-06-27 DOI: 10.1002/num.23119
K. Poochinapan, P. Manorot, T. Mouktonglang, B. Wongsaijai
{"title":"A linear finite difference scheme with error analysis designed to preserve the structure of the 2D Boussinesq paradigm equation","authors":"K. Poochinapan, P. Manorot, T. Mouktonglang, B. Wongsaijai","doi":"10.1002/num.23119","DOIUrl":"https://doi.org/10.1002/num.23119","url":null,"abstract":"Use of the finite difference method has produced successful solutions to the general partial differential equations due to its efficiency and effectiveness with wide applications. For example, the 2D Boussinesq paradigm equation can be numerically studied using a linear‐implicit finite difference scheme based on the Crank‐Nicolson/Adams‐Bashforth technique. First, conservative quantities are derived and preserved through numerical scheme. Then, the convergence and stability analysis is then provided to simulate a numerical solution whose existence and uniqueness are proved based on the boundedness of the numerical solution. Analysis of spatial accuracy is found to be second order on a uniform grid. Numerical results from simulations indicate that these proposed scheme provide satisfactory second‐order accuracy both in time and space with an ‐norm, and also preserve discrete invariants. Additionally, previous scientific literature review has provided little evidence of studied terms with dispersive effect in 2D Boussinesq paradigm equation. The current study explores solution behavior by applying the proposed scheme to numerically analyze initial Gaussian condition.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":"107 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-06-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141531971","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}
引用次数: 0
A flux‐based HDG method 基于通量的 HDG 方法
IF 3.9 3区 数学
Numerical Methods for Partial Differential Equations Pub Date : 2024-05-31 DOI: 10.1002/num.23117
Issei Oikawa
{"title":"A flux‐based HDG method","authors":"Issei Oikawa","doi":"10.1002/num.23117","DOIUrl":"https://doi.org/10.1002/num.23117","url":null,"abstract":"In this article, we present a flux‐based formulation of the hybridizable discontinuous Galerkin (HDG) method for steady‐state diffusion problems and propose a new method derived by letting a stabilization parameter tend to infinity. Assuming an inf‐sup condition, we prove its well‐posedness and error estimates of optimal order. We show that the inf‐sup condition is satisfied by some triangular elements. Numerical results are also provided to support our theoretical results.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":"1 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141191692","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}
引用次数: 0
Extensions and investigations of space‐time generalized Riemann problems numerical schemes for linear systems of conservation laws with source terms 带源项线性守恒定律系统的时空广义黎曼问题数值方案的扩展与研究
IF 3.9 3区 数学
Numerical Methods for Partial Differential Equations Pub Date : 2024-05-31 DOI: 10.1002/num.23118
Rodolphe Turpault
{"title":"Extensions and investigations of space‐time generalized Riemann problems numerical schemes for linear systems of conservation laws with source terms","authors":"Rodolphe Turpault","doi":"10.1002/num.23118","DOIUrl":"https://doi.org/10.1002/num.23118","url":null,"abstract":"The space‐time generalized Riemann problems method allows to obtain numerical schemes of arbitrary high order that can be used with very large time steps for systems of linear hyperbolic conservation laws with source term. They have been introduced in Berthon et al. (J. Sci. Comput. 55 (2013), 268–308.) in 1D and on 2D unstructured meshes made of triangles. The objective of this article is to complement them in order to answer some important questions arising when they are involved. The formulation is described in detail on quadrangle meshes, the choice of approximation basis is discussed and Legendre polynomials are used in practical cases. The addition of a limiter to preserve certain properties without compromising accuracy is also considered. Finally, the asymptotic behavior of the scheme in the diffusion regime is studied.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":"43 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141191684","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}
引用次数: 0
A parareal exponential integrator finite element method for semilinear parabolic equations 半线性抛物方程的准指数积分有限元法
IF 3.9 3区 数学
Numerical Methods for Partial Differential Equations Pub Date : 2024-05-29 DOI: 10.1002/num.23116
Jianguo Huang, Lili Ju, Yuejin Xu
{"title":"A parareal exponential integrator finite element method for semilinear parabolic equations","authors":"Jianguo Huang, Lili Ju, Yuejin Xu","doi":"10.1002/num.23116","DOIUrl":"https://doi.org/10.1002/num.23116","url":null,"abstract":"In this article, we present a parareal exponential finite element method, with the help of variational formulation and parareal framework, for solving semilinear parabolic equations in rectangular domains. The model equation is first discretized in space using the finite element method with continuous piecewise multilinear rectangular basis functions, producing the semi‐discrete system. We then discretize the temporal direction using the explicit exponential Runge–Kutta approach accompanied by the parareal framework, resulting in the fully‐discrete numerical scheme. To further improve computational speed, we design a fast solver for our method based on tensor product spectral decomposition and fast Fourier transform. Under certain regularity assumption, we successfully derive optimal error estimates for the proposed parallel‐based method with respect to ‐norm. Extensive numerical experiments in two and three dimensions are also carried out to validate the theoretical results and demonstrate the performance of our method.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":"55 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-05-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141191824","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}
引用次数: 0
Discontinuous Galerkin finite element method for dynamic viscoelasticity models of power‐law type 幂律型动态粘弹性模型的非连续 Galerkin 有限元方法
IF 3.9 3区 数学
Numerical Methods for Partial Differential Equations Pub Date : 2024-05-06 DOI: 10.1002/num.23107
Yongseok Jang, Simon Shaw
{"title":"Discontinuous Galerkin finite element method for dynamic viscoelasticity models of power‐law type","authors":"Yongseok Jang, Simon Shaw","doi":"10.1002/num.23107","DOIUrl":"https://doi.org/10.1002/num.23107","url":null,"abstract":"Linear viscoelasticity can be characterized by a stress relaxation function. We consider a power‐law type stress relaxation to yield a fractional order viscoelasticity model. The governing equation is a Volterra integral problem of the second kind with a weakly singular kernel. We employ spatially discontinuous Galerkin methods, <jats:italic>symmetric interior penalty Galerkin method</jats:italic> (SIPG) for spatial discretization, and the implicit finite difference schemes in time, <jats:italic>Crank–Nicolson method</jats:italic>. Further, in order to manage the weak singularity in the Volterra kernel, we use a linear interpolation technique. We present a priori stability and error analyses without relying on Grönwall's inequality, and so provide high quality bounds that do not increase exponentially in time. This indicates that our numerical scheme is well‐suited for long‐time simulations. Despite the limited regularity in time, we establish suboptimal fractional order accuracy in time as well as optimal convergence of SIPG. We carry out numerical experiments with varying regularity of exact solutions to validate our error estimates. Finally, we present numerical simulations based on real material data.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":"11 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-05-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140881626","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}
引用次数: 0
0
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
相关产品
×
本文献相关产品
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术官方微信