{"title":"Regularized lattice Boltzmann method based maximum principle and energy stability preserving finite-difference scheme for the Allen-Cahn equation","authors":"Ying Chen , Xi Liu , Zhenhua Chai , Baochang Shi","doi":"10.1016/j.jcp.2025.113831","DOIUrl":null,"url":null,"abstract":"<div><div>The Allen-Cahn equation (ACE) inherently possesses two crucial properties: the maximum principle and the energy dissipation law. How to preserve these two properties at the discrete level is of significant importance in the numerical methods for the ACE. In this paper, unlike the traditional macroscopic numerical schemes that directly discretize the ACE, we first propose a novel mesoscopic regularized lattice Boltzmann method based macroscopic numerical scheme for the <em>d</em> (=1, 2, 3)-dimensional ACE, where the D<em>d</em>Q<span><math><mo>(</mo><mn>2</mn><mi>d</mi><mo>+</mo><mn>1</mn><mo>)</mo></math></span> [(<span><math><mn>2</mn><mi>d</mi><mo>+</mo><mn>1</mn></math></span>) discrete velocities in <em>d</em>-dimensional space] lattice structure is employed. In particular, the proposed numerical scheme has a second-order accuracy in space, and can also be regarded as an implicit-explicit finite-difference scheme for the ACE. In this scheme, the nonlinear term is discretized semi-implicitly, the temporal derivative term and the dissipation term are discretized via the explicit Euler and second-order central difference methods, respectively. Compared to the implicit schemes for the ACE, the present scheme proves to be more efficient as it is actually explicit and avoids the use of iterative methods when dealing with the nonlinear term. In addition, we also demonstrate that under some certain conditions, the proposed scheme can preserve the maximum bound principle and the original energy dissipation law at the discrete level. Finally, we conduct numerical simulations of several benchmark problems, and find that the numerical results are not only in agreement with our theoretical analysis, but also show that in comparison with the mesoscopic lattice Boltzmann method, the proposed macroscopic scheme has a great advantage in reducing the memory usage.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"528 ","pages":"Article 113831"},"PeriodicalIF":3.8000,"publicationDate":"2025-02-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational Physics","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0021999125001147","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
The Allen-Cahn equation (ACE) inherently possesses two crucial properties: the maximum principle and the energy dissipation law. How to preserve these two properties at the discrete level is of significant importance in the numerical methods for the ACE. In this paper, unlike the traditional macroscopic numerical schemes that directly discretize the ACE, we first propose a novel mesoscopic regularized lattice Boltzmann method based macroscopic numerical scheme for the d (=1, 2, 3)-dimensional ACE, where the DdQ [() discrete velocities in d-dimensional space] lattice structure is employed. In particular, the proposed numerical scheme has a second-order accuracy in space, and can also be regarded as an implicit-explicit finite-difference scheme for the ACE. In this scheme, the nonlinear term is discretized semi-implicitly, the temporal derivative term and the dissipation term are discretized via the explicit Euler and second-order central difference methods, respectively. Compared to the implicit schemes for the ACE, the present scheme proves to be more efficient as it is actually explicit and avoids the use of iterative methods when dealing with the nonlinear term. In addition, we also demonstrate that under some certain conditions, the proposed scheme can preserve the maximum bound principle and the original energy dissipation law at the discrete level. Finally, we conduct numerical simulations of several benchmark problems, and find that the numerical results are not only in agreement with our theoretical analysis, but also show that in comparison with the mesoscopic lattice Boltzmann method, the proposed macroscopic scheme has a great advantage in reducing the memory usage.
期刊介绍:
Journal of Computational Physics thoroughly treats the computational aspects of physical problems, presenting techniques for the numerical solution of mathematical equations arising in all areas of physics. The journal seeks to emphasize methods that cross disciplinary boundaries.
The Journal of Computational Physics also publishes short notes of 4 pages or less (including figures, tables, and references but excluding title pages). Letters to the Editor commenting on articles already published in this Journal will also be considered. Neither notes nor letters should have an abstract.