防止标准水平集方法中的质量损失:变分分析的新见解

IF 3.8 2区 物理与天体物理 Q2 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS
Kaustubh Khedkar , Amirreza Charchi Mamaghani , Pieter Ghysels , Neelesh A. Patankar , Amneet Pal Singh Bhalla
{"title":"防止标准水平集方法中的质量损失:变分分析的新见解","authors":"Kaustubh Khedkar ,&nbsp;Amirreza Charchi Mamaghani ,&nbsp;Pieter Ghysels ,&nbsp;Neelesh A. Patankar ,&nbsp;Amneet Pal Singh Bhalla","doi":"10.1016/j.jcp.2024.113495","DOIUrl":null,"url":null,"abstract":"<div><div>For decades, the computational multiphase flow community has grappled with mass loss in the level set method. Numerous solutions have been proposed, from fixing the reinitialization step to combining the level set method with other conservative schemes. However, our work reveals a more fundamental culprit: the smooth Heaviside and delta functions inherent to the standard formulation. Even if reinitialization is done exactly, i.e., the zero contour interface remains stationary, the use of smooth functions lead to violation of mass conservation. We propose a novel approach using variational analysis to incorporate a mass conservation constraint. This introduces a Lagrange multiplier that enforces overall mass balance. Notably, as the delta function sharpens, i.e., approaches the Dirac delta limit, the Lagrange multiplier approaches zero. However, the exact Lagrange multiplier method disrupts the signed distance property of the level set function. This motivates us to develop an approximate version of the Lagrange multiplier that preserves both overall mass and signed distance property of the level set function. Our framework even recovers existing mass-conserving level set methods, revealing some inconsistencies in prior analyses. We extend this approach to three-phase flows for fluid-structure interaction (FSI) simulations. We present variational equations in both immersed and non-immersed forms, demonstrating the convergence of the former formulation to the latter when the body delta function sharpens. Rigorous test problems confirm that the FSI dynamics produced by our simple, easy-to-implement immersed formulation with the approximate Lagrange multiplier method are accurate and match state-of-the-art solvers.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"520 ","pages":"Article 113495"},"PeriodicalIF":3.8000,"publicationDate":"2024-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Preventing mass loss in the standard level set method: New insights from variational analyses\",\"authors\":\"Kaustubh Khedkar ,&nbsp;Amirreza Charchi Mamaghani ,&nbsp;Pieter Ghysels ,&nbsp;Neelesh A. Patankar ,&nbsp;Amneet Pal Singh Bhalla\",\"doi\":\"10.1016/j.jcp.2024.113495\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><div>For decades, the computational multiphase flow community has grappled with mass loss in the level set method. Numerous solutions have been proposed, from fixing the reinitialization step to combining the level set method with other conservative schemes. However, our work reveals a more fundamental culprit: the smooth Heaviside and delta functions inherent to the standard formulation. Even if reinitialization is done exactly, i.e., the zero contour interface remains stationary, the use of smooth functions lead to violation of mass conservation. We propose a novel approach using variational analysis to incorporate a mass conservation constraint. This introduces a Lagrange multiplier that enforces overall mass balance. Notably, as the delta function sharpens, i.e., approaches the Dirac delta limit, the Lagrange multiplier approaches zero. However, the exact Lagrange multiplier method disrupts the signed distance property of the level set function. This motivates us to develop an approximate version of the Lagrange multiplier that preserves both overall mass and signed distance property of the level set function. Our framework even recovers existing mass-conserving level set methods, revealing some inconsistencies in prior analyses. We extend this approach to three-phase flows for fluid-structure interaction (FSI) simulations. We present variational equations in both immersed and non-immersed forms, demonstrating the convergence of the former formulation to the latter when the body delta function sharpens. Rigorous test problems confirm that the FSI dynamics produced by our simple, easy-to-implement immersed formulation with the approximate Lagrange multiplier method are accurate and match state-of-the-art solvers.</div></div>\",\"PeriodicalId\":352,\"journal\":{\"name\":\"Journal of Computational Physics\",\"volume\":\"520 \",\"pages\":\"Article 113495\"},\"PeriodicalIF\":3.8000,\"publicationDate\":\"2024-10-11\",\"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/S0021999124007435\",\"RegionNum\":2,\"RegionCategory\":\"物理与天体物理\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q2\",\"JCRName\":\"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational Physics","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0021999124007435","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0

摘要

几十年来,计算多相流界一直在努力解决水平集方法中的质量损失问题。从固定重新初始化步骤到将水平集方法与其他保守方案相结合,已经提出了许多解决方案。然而,我们的工作揭示了一个更根本的罪魁祸首:标准公式固有的平滑海维塞函数和德尔塔函数。即使重新初始化精确完成,即零轮廓界面保持静止,光滑函数的使用也会导致违反质量守恒。我们提出了一种使用变分分析的新方法,将质量守恒约束纳入其中。这就引入了一个拉格朗日乘数,强制实现整体质量平衡。值得注意的是,随着三角函数的尖锐化,即接近狄拉克三角极限,拉格朗日乘数趋近于零。然而,精确的拉格朗日乘法破坏了水平集函数的符号距离特性。这促使我们开发一种近似版本的拉格朗日乘法器,它既能保持总体质量,又能保持水平集函数的符号距离特性。我们的框架甚至恢复了现有的质量保证水平集方法,揭示了之前分析中的一些不一致之处。我们将这种方法扩展到流固耦合(FSI)模拟的三相流。我们提出了沉浸式和非沉浸式两种形式的变分方程,证明了当体三角函数变得尖锐时,前者的公式会向后者收敛。严格的测试问题证实,我们简单易用的浸没式公式与近似拉格朗日乘法器方法所产生的 FSI 动力学是精确的,与最先进的求解器不相上下。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
Preventing mass loss in the standard level set method: New insights from variational analyses
For decades, the computational multiphase flow community has grappled with mass loss in the level set method. Numerous solutions have been proposed, from fixing the reinitialization step to combining the level set method with other conservative schemes. However, our work reveals a more fundamental culprit: the smooth Heaviside and delta functions inherent to the standard formulation. Even if reinitialization is done exactly, i.e., the zero contour interface remains stationary, the use of smooth functions lead to violation of mass conservation. We propose a novel approach using variational analysis to incorporate a mass conservation constraint. This introduces a Lagrange multiplier that enforces overall mass balance. Notably, as the delta function sharpens, i.e., approaches the Dirac delta limit, the Lagrange multiplier approaches zero. However, the exact Lagrange multiplier method disrupts the signed distance property of the level set function. This motivates us to develop an approximate version of the Lagrange multiplier that preserves both overall mass and signed distance property of the level set function. Our framework even recovers existing mass-conserving level set methods, revealing some inconsistencies in prior analyses. We extend this approach to three-phase flows for fluid-structure interaction (FSI) simulations. We present variational equations in both immersed and non-immersed forms, demonstrating the convergence of the former formulation to the latter when the body delta function sharpens. Rigorous test problems confirm that the FSI dynamics produced by our simple, easy-to-implement immersed formulation with the approximate Lagrange multiplier method are accurate and match state-of-the-art solvers.
求助全文
通过发布文献求助,成功后即可免费获取论文全文。 去求助
来源期刊
Journal of Computational Physics
Journal of Computational Physics 物理-计算机:跨学科应用
CiteScore
7.60
自引率
14.60%
发文量
763
审稿时长
5.8 months
期刊介绍: 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.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
copy
已复制链接
快去分享给好友吧!
我知道了
右上角分享
点击右上角分享
0
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术官方微信