ENO守恒插值(重映)方法及其在流体计算中的应用

在ENO(Essentially Non-oscillatory)守恒插值方法的基础上,分析和研究现今流体力学计算中涉及的几类网格技术:重叠网格技术、自适应加密技术和运动网格技术.基于ENO插值多项式构造的重映方法具有良好的守恒性,可以有效保证数据传递中物理量的总体守恒.提出该类守恒插值方法在以上几种网格技术中的一些应用前景,并给出一些数值算例.     

Two-step hybrid conservative remapping for multimaterial arbitrary Lagrangian–Eulerian methods

We present a new hybrid conservative remapping algorithm for multimaterial Arbitrary Lagrangian–Eulerian (ALE) methods. The hybrid remapping is performed in two steps. In the first step, only nodes of the grid that lie inside subdomains occupied by single materials are moved. At this stage, computationally cheap swept-region remapping is used. In the second step, nodes that are vertices of mixed cells (cells containing several materials) and vertices of some cells in a buffer zone around mixed cells are moved. At this stage, intersection-based remapping is used. The hybrid algorithm results in computational expense that lies between swept-region and intersection-based remapping We demonstrate the performance of our new method for both structured and unstructured polygonal grids in two dimensions, as well as for cell-centered and staggered discretizations.     

Numerical simulations of instabilities in the implosion process of inertial confined fusion in 2D cylindrical coordinates

In this paper, we introduce a multi-material arbitrary Lagrangian and Eulerian method for the hydrodynamic radiative multi-group diffusion model in 2D cylindrical coordinates. The basic idea in the construction of the method is the following: In the Lagrangian step, a closure model of radiation-hydrodynamics is used to give the states of equations for materials in mixed cells. In the mesh rezoning step, we couple the rezoning principle with the Lagrangian interface tracking method and an Eulerian interface capturing scheme to compute interfaces sharply according to their deformation and to keep cells in good geometric quality. In the interface reconstruction step, a dual-material Moment-of-Fluid method is introduced to obtain the unique interface in mixed cells. In the remapping step, a conservative remapping algorithm of conserved quantities is presented. A number of numerical tests are carried out and the numerical results show that the new method can simulate instabilities in complex fluid field under large deformation,and are accurate and robust.     

一类基于RBF插值的守恒重映算法

为保证重映过程的高守恒精度和单调性,并且在间断处具有极高的分辨率,基于径向基函数(RBF)插值方法构造了一类适用于任意网格的RBF守恒重映算法,通过计算守恒误差测试重映算法的守恒精度。将该方法用于光滑函数和含有间断的函数,并与其它守恒重映方法比较,表明该方法数值结果较好。     

Metric-based mesh adaptation for 2D Lagrangian compressible flows

In this paper, we present a method to compute compressible flows in 2D. It uses two steps: a Lagrangian step and a metric-based triangular mesh adaptation step. Computational mesh is locally adapted according to some metric field that depends on physical or geometrical data. This mesh adaptation step embeds a conservative remapping procedure to satisfy consistency with Euler equations. The whole method is no more Lagrangian.After describing mesh adaptation patterns, we recall the metric formalism. Then, we detail an appropriate remapping procedure which is first-order and relies on exact intersections.We give some hints about the parallel implementation. Finally, we present various numerical experiments which demonstrate the good properties of the algorithm.     

A high order cell-centered semi-Lagrangian scheme for multi-dimensional kinetic simulations of neutral gas flows

The term ‘Convected Scheme’ (CS) refers to a family of algorithms, most usually applied to the solution of Boltzmann’s equation, which uses a method of characteristics in an integral form to project an initial cell forward to a group of final cells. As such the CS is a ‘forward-trajectory’ semi-Lagrangian scheme. For multi-dimensional simulations of neutral gas flows, the cell-centered version of this semi-Lagrangian (CCSL) scheme has advantages over other options due to its implementation simplicity, low memory requirements, and easier treatment of boundary conditions. The main drawback of the CCSL-CS to date has been its high numerical diffusion in physical space, because of the 2nd order remapping that takes place at the end of each time step. By means of a modified equation analysis, it is shown that a high order estimate of the remapping error can be obtained a priori, and a small correction to the final position of the cells can be applied upon remapping, in order to achieve full compensation of this error. The resulting scheme is 4th order accurate in space while retaining the desirable properties of the CS: it is conservative and positivity-preserving, and the overall algorithm complexity is not appreciably increased. Two monotone (i.e. non-oscillating) versions of the fourth order CCSL-CS are also presented: one uses a common flux-limiter approach; the other uses a non-polynomial reconstruction to evaluate the derivatives of the density function. The method is illustrated in simple one- and two-dimensional examples, and a fully 3D solution of the Boltzmann equation describing expansion of a gas into vacuum through a cylindrical tube.     

多边形交错网格的守恒重映算法

提出基于细分和数值积分思想的一种离散的守恒重映方法——质点重映方法.密度分布可采用一阶精度的分片常数分布,或二阶精度的分片线性分布.分片线性密度分布函数采用面平均方法构造.重映过程中,借助四边形辅助网格,实现了交错网格节点量的重映.质点重映方法既适用于结构网格,也适用于非结构网格,且不要求新旧网格之间一一对应.数值结果表明,一阶精度重映算法健壮性好,但会产生较大的扩散效应;二阶精度重映算法可较好地保持密度分布的特性,但存在单调性问题.为改善二阶精度重映方法单调性,将结构网格质量守恒调整算法推广到非结构网格上,以限制新网格的质量密度.给出了一些重映的例子,并进行了误差分析.     

一类基于ENO插值的守恒重映算法

在大变形流体力学问题的数值模拟中,常常会涉及到计算网格的重分.基于不同网格的物理量传递便是所谓的重映技术.基于ENO插值的思想,发展了一类适用于任意网格的ENO守恒重映算法,并给出了数值结果.     

一种基于MOF界面重构的二维中心型MMALE方法

发展了一种基于MOF（Moment of Fluid）界面重构的二维中心型MMALE（Multi-Material Arbitrary Lagrangian-Eulerian）方法.其中,流体力学方程组采用中心型拉氏方法进行离散求解.混合网格的热力学封闭采用Tipton压力松弛模型.混合网格内的界面重构采用MOF方法,并对MOF方法作了简化和改进.重映步采用一种基于多边形剪裁算法的精确积分守恒重映方法.计算了若干数值例子,包括二维漩涡发展问题、Sedov问题、激波与氦气泡相互作用问题、水中强激波与空气泡相互作用问题、二维RT不稳定性问题等.数值算例表明,该方法具有二阶精度,能够计算界面两侧密度比和压力比很大的问题,并且其健壮性优于交错型MMALE方法,适合计算多介质复杂流体动力学问题.     

分片抛物线对流守恒重映方法

引进一种守恒的分片抛物线对流重映方法,通过交替扫描平均法提高对流重映方法的对称性,使用分片抛物线分布函数提高对流重映方法的精度.给出一维算例和二维算例检验分片抛物线对流重映方法的精度和对称性.     

非结构多面体二阶局部保界全局重映算法

提出一种三维非结构多面体二阶保界全局重映算法.在旧网格上选取模板利用最小二乘构造插值多项式,采用凸包算法计算多面体相交部分,最后使用局部保界修正技术修补重映后的越界量.多项数值实验表明这种格式同时具有高精度、高分辨率和高效率的特点.     

二维Lagrange网格的积分守恒重映方法

研究了解决二维Lagrange网格大变形的一种有效的网格重分方法——积分守恒重映方法.详细地介绍了算法,并给出了数值实验结果.     

多尺度模拟中网格守恒重映算法

针对耦合微观分子动力学(MD)和宏观有限元方法(FE)的多尺度模拟,提出一类新的基于贡献单元法的网格守恒重映算法.由于物理量是由有限元节点以及相应区域的原子信息通过积分重构得到的,对结构和非结构网格都能适用.对于未知量定义在顶点的情形,引入辅助网格.数值例子验证了算法的准确性和有效性.     

激光靶耦合问题的二阶精度重映算法研究

 任意拉格朗日欧拉方法被认为是解决2维流体力学大变形数值模拟的有效途径之一。以激光黑腔靶耦合问题为背景,结合质量流的物理意义,提出了一种二阶守恒物理量重映的算法,并利用1维激光靶耦合程序对方法及其可行性进行了考察,结合物理问题调整网格,用较少的网格得到了与较多网格拉氏计算同样的结果。为算法在2维LARED-H程序中的应用奠定了基础。     

One-step hybrid remapping algorithm for multi-material arbitrary Lagrangian–Eulerian methods

In this paper, a new flux-based one-step hybrid remapping method for multi-material arbitrary Lagrangian–Eulerian (ALE) approach is introduced. In the vicinity of material interfaces, the swept region is intersected with pure material polygons in the Lagrangian mesh to construct the material fluxes. Far from interfaces, the fluxes are constructed in a standard swept-region manner without intersections. This method is conservative, second-order accurate and linearity-preserving (in case of straight material interfaces), and faster than method based on intersections, as shown on selected numerical examples.     

基于MOF界面重构的多物质ALE方法

提出一种基于MOF(Moment-of-Fluid)界面重构的多物质ALE(Arbitrary Lagrangian-Eulerian)方法.流体力学方程组采用相容有限元方法进行空间离散.提出一种新的二维子网格力学模型,用来计算混合网格中的物理量经过一个拉氏步后发生的变化,混合网格内的界面重构采用MOF方法.提出一种精确积分守恒重映方法.给出数值算例,如空气和水的Riemann问题,Dukowicz问题,水中强激波与空气泡相互作用问题等.结果表明,方法具有较高的精度,能够处理物质界面和网格的大变形问题.     

An Eulerian–Lagrangian WENO finite volume scheme for advection problems

We develop a locally conservative Eulerian–Lagrangian finite volume scheme with the weighted essentially non-oscillatory property (EL–WENO) in one-space dimension. This method has the advantages of both WENO and Eulerian–Lagrangian schemes. It is formally high-order accurate in space (we present the fifth order version) and essentially non-oscillatory. Moreover, it is free of a CFL time step stability restriction and has small time truncation error. The scheme requires a new integral-based WENO reconstruction to handle trace-back integration. A Strang splitting algorithm is presented for higher-dimensional problems, using both the new integral-based and pointwise-based WENO reconstructions. We show formally that it maintains the fifth order accuracy. It is also locally mass conservative. Numerical results are provided to illustrate the performance of the scheme and verify its formal accuracy.     

SALE速度重映算法的改进

在一维交错网格上基于SALE算法的速度重映策略,提出了3种动量通量的改进算法:①采用迎风斜率加修补的SUR目的;②采用minmod斜率限制器的SM目的和③采用一种新的斜率限制器加修补的SLR目的.3种目的具有二阶守恒保界的性质,同时继承了SALE算法简单高效的特点,可以直接推广到二维情况.     

任意网格重映的样条逼近算法

在大变形流体力学问题的数值模拟中,任何方法都必须考虑网格重分或网格自适应,只要改动网格就涉及重分,或自适应后从旧的、扭曲的网格到新网格的守恒量重映,包括质量、动量和能量.在研究样条函数逼近的基础上,给出一种物理量重映的对结构网格和非结构网格均适应的算法,并给出了数值结果.     

Optimization-based synchronized flux-corrected conservative interpolation (remapping) of mass and momentum for arbitrary Lagrangian–Eulerian methods

A new optimization-based synchronized flux-corrected conservative interpolation (remapping) of mass and momentum for arbitrary Lagrangian–Eulerian hydro methods is described. Fluxes of conserved variables – mass and momentum – are limited in a synchronous way to preserve local bounds of primitive variables – density and velocity.