{"title":"A high-order accurate unconditionally stable bound-preserving numerical scheme for the Cahn-Hilliard-Navier-Stokes equations","authors":"Yali Gao , Daozhi Han , Sayantan Sarkar","doi":"10.1016/j.apnum.2025.06.004","DOIUrl":null,"url":null,"abstract":"<div><div>A high order numerical method is developed for solving the Cahn-Hilliard-Navier-Stokes equations with the Flory-Huggins potential. The scheme is based on the <span><math><msub><mrow><mi>Q</mi></mrow><mrow><mi>k</mi></mrow></msub></math></span> finite element with mass lumping on rectangular grids, the second-order convex splitting method and the pressure correction method. The unique solvability, unconditional stability, and bound-preserving properties are rigorously established. The key for bound-preservation is the discrete <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> estimate of the singular potential. Ample numerical experiments are performed to validate the desired properties of the proposed numerical scheme.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"217 ","pages":"Pages 96-111"},"PeriodicalIF":2.4000,"publicationDate":"2025-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Applied Numerical Mathematics","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0168927425001217","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 0
Abstract
A high order numerical method is developed for solving the Cahn-Hilliard-Navier-Stokes equations with the Flory-Huggins potential. The scheme is based on the finite element with mass lumping on rectangular grids, the second-order convex splitting method and the pressure correction method. The unique solvability, unconditional stability, and bound-preserving properties are rigorously established. The key for bound-preservation is the discrete estimate of the singular potential. Ample numerical experiments are performed to validate the desired properties of the proposed numerical scheme.
期刊介绍:
The purpose of the journal is to provide a forum for the publication of high quality research and tutorial papers in computational mathematics. In addition to the traditional issues and problems in numerical analysis, the journal also publishes papers describing relevant applications in such fields as physics, fluid dynamics, engineering and other branches of applied science with a computational mathematics component. The journal strives to be flexible in the type of papers it publishes and their format. Equally desirable are:
(i) Full papers, which should be complete and relatively self-contained original contributions with an introduction that can be understood by the broad computational mathematics community. Both rigorous and heuristic styles are acceptable. Of particular interest are papers about new areas of research, in which other than strictly mathematical arguments may be important in establishing a basis for further developments.
(ii) Tutorial review papers, covering some of the important issues in Numerical Mathematics, Scientific Computing and their Applications. The journal will occasionally publish contributions which are larger than the usual format for regular papers.
(iii) Short notes, which present specific new results and techniques in a brief communication.