{"title":"基于移动网格求解器的高阶保正性非连续伽勒金重绘方法,用于可压缩流体流动的 ALE 仿真","authors":"Xiaolu Gu,Juan Cheng, Chiwang Shu","doi":"10.4208/cicp.oa-2023-0083","DOIUrl":null,"url":null,"abstract":"The arbitrary Lagrangian-Eulerian (ALE) method is widely used in the field\nof compressible multi-material and multi-phase flow problems. In order to implement the indirect ALE approach for the simulation of compressible flow in the context\nof high order discontinuous Galerkin (DG) discretizations, we present a high order\npositivity-preserving DG remapping method based on a moving mesh solver in this\npaper. This remapping method is based on the ALE-DG method developed by Klingenberg et al. [17, 18] to solve the trivial equation $\\frac{∂u}{∂t} = 0$ on a moving mesh, which\nis the old mesh before remapping at $t = 0$ and is the new mesh after remapping at $t = T.$ An appropriate selection of the final pseudo-time $T$ can always satisfy the relatively mild smoothness requirement (Lipschitz continuity) on the mesh movement\nvelocity, which guarantees the high order accuracy of the remapping procedure. We\nuse a multi-resolution weighted essentially non-oscillatory (WENO) limiter which can\nkeep the essentially non-oscillatory property near strong discontinuities while maintaining high order accuracy in smooth regions. We further employ an effective linear\nscaling limiter to preserve the positivity of the relevant physical variables without sacrificing conservation and the original high order accuracy. Numerical experiments are\nprovided to illustrate the high order accuracy, essentially non-oscillatory performance\nand positivity-preserving of our remapping algorithm. In addition, the performance\nof the ALE simulation based on the DG framework with our remapping algorithm is\nexamined in one- and two-dimensional Euler equations.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"25 1","pages":""},"PeriodicalIF":2.6000,"publicationDate":"2023-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"A High Order Positivity-Preserving Discontinuous Galerkin Remapping Method Based on a Moving Mesh Solver for ALE Simulation of the Compressible Fluid Flow\",\"authors\":\"Xiaolu Gu,Juan Cheng, Chiwang Shu\",\"doi\":\"10.4208/cicp.oa-2023-0083\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"The arbitrary Lagrangian-Eulerian (ALE) method is widely used in the field\\nof compressible multi-material and multi-phase flow problems. In order to implement the indirect ALE approach for the simulation of compressible flow in the context\\nof high order discontinuous Galerkin (DG) discretizations, we present a high order\\npositivity-preserving DG remapping method based on a moving mesh solver in this\\npaper. This remapping method is based on the ALE-DG method developed by Klingenberg et al. [17, 18] to solve the trivial equation $\\\\frac{∂u}{∂t} = 0$ on a moving mesh, which\\nis the old mesh before remapping at $t = 0$ and is the new mesh after remapping at $t = T.$ An appropriate selection of the final pseudo-time $T$ can always satisfy the relatively mild smoothness requirement (Lipschitz continuity) on the mesh movement\\nvelocity, which guarantees the high order accuracy of the remapping procedure. We\\nuse a multi-resolution weighted essentially non-oscillatory (WENO) limiter which can\\nkeep the essentially non-oscillatory property near strong discontinuities while maintaining high order accuracy in smooth regions. We further employ an effective linear\\nscaling limiter to preserve the positivity of the relevant physical variables without sacrificing conservation and the original high order accuracy. Numerical experiments are\\nprovided to illustrate the high order accuracy, essentially non-oscillatory performance\\nand positivity-preserving of our remapping algorithm. In addition, the performance\\nof the ALE simulation based on the DG framework with our remapping algorithm is\\nexamined in one- and two-dimensional Euler equations.\",\"PeriodicalId\":50661,\"journal\":{\"name\":\"Communications in Computational Physics\",\"volume\":\"25 1\",\"pages\":\"\"},\"PeriodicalIF\":2.6000,\"publicationDate\":\"2023-12-01\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Communications in Computational Physics\",\"FirstCategoryId\":\"101\",\"ListUrlMain\":\"https://doi.org/10.4208/cicp.oa-2023-0083\",\"RegionNum\":3,\"RegionCategory\":\"物理与天体物理\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q1\",\"JCRName\":\"PHYSICS, MATHEMATICAL\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Communications in Computational Physics","FirstCategoryId":"101","ListUrlMain":"https://doi.org/10.4208/cicp.oa-2023-0083","RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"PHYSICS, MATHEMATICAL","Score":null,"Total":0}
引用次数: 0
摘要
任意拉格朗日-欧勒(ALE)方法被广泛应用于可压缩多材料和多相流问题领域。为了在高阶非连续伽勒金(DG)离散化背景下实现可压缩流动模拟的间接 ALE 方法,我们在本文中提出了一种基于移动网格求解器的高阶正向保留 DG 重映射方法。这种重映射方法基于 Klingenberg 等人开发的 ALE-DG 方法[17, 18],在移动网格上求解三元方程 $\frac{∂u}{∂t} = 0$,移动网格是在 $t = 0$ 时重映射前的旧网格和在 $t = T 时重映射后的新网格。最终伪时间 $T$ 的适当选择总能满足对网格移动速度相对温和的平滑性要求(Lipschitz 连续性),从而保证重映射过程的高阶精度。我们使用了多分辨率加权本质非振荡(WENO)限制器,它可以在强不连续性附近保持本质非振荡特性,同时在平滑区域保持高阶精度。我们进一步采用了有效的线性缩放限制器,在不牺牲守恒性和原有高阶精度的情况下,保持相关物理变量的正向性。我们提供了数值实验,以说明我们的重映射算法具有高阶精度、基本无振荡性能和保留正性的特点。此外,还在一维和二维欧拉方程中检验了基于 DG 框架的 ALE 仿真与我们的重映射算法的性能。
A High Order Positivity-Preserving Discontinuous Galerkin Remapping Method Based on a Moving Mesh Solver for ALE Simulation of the Compressible Fluid Flow
The arbitrary Lagrangian-Eulerian (ALE) method is widely used in the field
of compressible multi-material and multi-phase flow problems. In order to implement the indirect ALE approach for the simulation of compressible flow in the context
of high order discontinuous Galerkin (DG) discretizations, we present a high order
positivity-preserving DG remapping method based on a moving mesh solver in this
paper. This remapping method is based on the ALE-DG method developed by Klingenberg et al. [17, 18] to solve the trivial equation $\frac{∂u}{∂t} = 0$ on a moving mesh, which
is the old mesh before remapping at $t = 0$ and is the new mesh after remapping at $t = T.$ An appropriate selection of the final pseudo-time $T$ can always satisfy the relatively mild smoothness requirement (Lipschitz continuity) on the mesh movement
velocity, which guarantees the high order accuracy of the remapping procedure. We
use a multi-resolution weighted essentially non-oscillatory (WENO) limiter which can
keep the essentially non-oscillatory property near strong discontinuities while maintaining high order accuracy in smooth regions. We further employ an effective linear
scaling limiter to preserve the positivity of the relevant physical variables without sacrificing conservation and the original high order accuracy. Numerical experiments are
provided to illustrate the high order accuracy, essentially non-oscillatory performance
and positivity-preserving of our remapping algorithm. In addition, the performance
of the ALE simulation based on the DG framework with our remapping algorithm is
examined in one- and two-dimensional Euler equations.
期刊介绍:
Communications in Computational Physics (CiCP) publishes original research and survey papers of high scientific value in computational modeling of physical problems. Results in multi-physics and multi-scale innovative computational methods and modeling in all physical sciences will be featured.