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.     

ReALE: A reconnection-based arbitrary-Lagrangian–Eulerian method

We present a new reconnection-based arbitrary-Lagrangian–Eulerian (ALE) method. The main elements in a standard ALE simulation are an explicit Lagrangian phase in which the solution and grid are updated, a rezoning phase in which a new grid is defined, and a remapping phase in which the Lagrangian solution is transferred (conservatively interpolated) onto the new grid. In standard ALE methods the new mesh from the rezone phase is obtained by moving grid nodes without changing connectivity of the mesh. Such rezone strategy has its limitation due to the fixed topology of the mesh. In our new method we allow connectivity of the mesh to change in rezone phase, which leads to general polygonal mesh and allows to follow Lagrangian features of the mesh much better than for standard ALE methods. Rezone strategy with reconnection is based on using Voronoi tessellation. We demonstrate performance of our new method on series of numerical examples and show it superiority in comparison with standard ALE methods without reconnection.     

Remapping-free ALE-type kinetic method for flow computations

Based on the integral form of the fluid dynamic equations, a finite volume kinetic scheme with arbitrary control volume and mesh velocity is developed. Different from the earlier unified moving mesh gas-kinetic method [C.Q. Jin, K. Xu, An unified moving grid gas-kinetic method in Eulerian space for viscous flow computation, J. Comput. Phys. 222 (2007) 155–175], the coupling of the fluid equations and geometrical conservation laws has been removed in order to make the scheme applicable for any quadrilateral or unstructured mesh rather than parallelogram in 2D case. Since a purely Lagrangian method is always associated with mesh entangling, in order to avoid computational collapsing in multidimensional flow simulation, the mesh velocity is constructed by considering both fluid velocity (Lagrangian methodology) and diffusive velocity (Regenerating Eulerian mesh function). Therefore, we obtain a generalized Arbitrary-Lagrangian–Eulerian (ALE) method by properly designing a mesh velocity instead of re-generating a new mesh after distortion. As a result, the remapping step to interpolate flow variables from old mesh to new mesh is avoided. The current method provides a general framework, which can be considered as a remapping-free ALE-type method. Since there is great freedom in choosing mesh velocity, in order to improve the accuracy and robustness of the method, the adaptive moving mesh method [H.Z. Tang, T. Tang, Adaptive mesh methods for one-and two-dimensional hyperbolic conservation laws, SIAM J. Numer. Anal. 41 (2003) 487–515] can be also used to construct a mesh velocity to concentrate mesh to regions with high flow gradients.     

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

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

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

任意拉格朗日欧拉方法被认为是解决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.     

基于任意拉格朗日—欧拉型移动网格的反应流广义黎曼问题方法

提出一种在自由重映移动网格下的广义黎曼问题方法模拟反应流.该方法基于显式的自由重映移动网格广义黎曼问题的解.为保证在时间和空间上的高精度，应用广义黎曼问题方法构造数值通量.为保证反应区的高分辨率，采用变分法生成自适应移动网格.该方法不仅能够保证网格质量，而且能有效地避免任意拉格朗日—欧拉方法中由于显式重映过程而带来的数值误差.包括CJ爆轰及不稳定爆轰的数值实验说明该格式的精确性和鲁棒性，证明这种移动网格下的二阶广义黎曼问题方法可以较好地捕捉反应流的间断与光滑结构.     

Semi-implicit surface tension formulation with a Lagrangian surface mesh on an Eulerian simulation grid

We present a method for applying semi-implicit forces on a Lagrangian mesh to an Eulerian discretization of the Navier Stokes equations in a way that produces a sparse symmetric positive definite system. The resulting method has semi-implicit and fully-coupled viscosity, pressure, and Lagrangian forces. We apply our new framework for forces on a Lagrangian mesh to the case of a surface tension force, which when treated explicitly leads to a tight time step restriction. By applying surface tension as a semi-implicit Lagrangian force, the resulting method benefits from improved stability and the ability to take larger time steps. The resulting discretization is also able to maintain parasitic currents at low levels.     

A new anisotropic mesh adaptation method based upon hierarchical a posteriori error estimates

A new anisotropic mesh adaptation strategy for finite element solution of elliptic differential equations is presented. It generates anisotropic adaptive meshes as quasi-uniform ones in some metric space, with the metric tensor being computed based on hierarchical a posteriori error estimates. A global hierarchical error estimate is employed in this study to obtain reliable directional information of the solution. Instead of solving the global error problem exactly, which is costly in general, we solve it iteratively using the symmetric Gauß–Seidel method. Numerical results show that a few GS iterations are sufficient for obtaining a reasonably good approximation to the error for use in anisotropic mesh adaptation. The new method is compared with several strategies using local error estimators or recovered Hessians. Numerical results are presented for a selection of test examples and a mathematical model for heat conduction in a thermal battery with large orthotropic jumps in the material coefficients.     

Second-order sign-preserving conservative interpolation (remapping) on general grids

An accurate conservative interpolation (remapping) algorithm is an essential component of most arbitrary Lagrangian–Eulerian (ALE) methods. In this paper we describe a local remapping algorithm for a positive scalar function. This algorithm is second-order accurate, conservative, and sign preserving. The algorithm is based on estimating the mass exchanged between cells at their common interface, and so is equally applicable to structured and unstructured grids. We construct the algorithm in a series of steps, clearly delineating the assumptions and errors made at each step. We validate our theory with a suite of numerical examples, analyzing the results from the viewpoint of accuracy and order of convergence.     

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

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

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

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

一种基于预估-修正策略的非定常流网格自适应加密/稀疏方法

给出一种非定常流动数值模拟的网格自适应处理方法.在"求解流动方程-自适应调整网格"的流程中,引入预估-修正步.根据自适应周期内每个时间步上的流场预估解,计算单元上的事后误差估算值.建立考虑解演变的网格自适应指示器,并进行多层次单元加密-稀疏的动态网格自适应处理.在自适应网格上重新计算流场.每个自适应周期中,流动演变区域的网格获得加密;而前一个周期中的特征现象已离开区域的网格被稀疏.应用边界非协调的当地DFD(Domain-Free Discretization)方法求解流动方程.为验证网格自适应处理方法,针对静止圆柱和自推进游鱼的流动进行了数值实验.     

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.     

Coupling Parameters and the Form of the Potential via Noether Symmetry

We explore the conditions for the existence of Noether symmetries in the dynamics of FRW metric, non minimally coupled with a scalar field, in the most general situation, and with nonzero spatial curvature. When such symmetries are present we find a general exact solution for the Einstein equations. We also show that non Noether symmetries can be found. Finally, we present an extension of the procedure to the Kantowski-Sachs metric which is particularly interesting in the case of degenerate Lagrangian.     

Metric construction by length distribution tensor and edge based error for anisotropic adaptive meshing

Metric tensors play a key role to control the generation of unstructured anisotropic meshes. In practice, the most well established error analysis enables to calculate a metric tensor on an element basis. In this paper, we propose to build a metric field directly at the nodes of the mesh for a direct use in the meshing tools. First, the unit mesh metric is defined and well justified on a node basis, by using the statistical concept of length distribution tensors. Then, the interpolation error analysis is performed on the projected approximate scalar field along the edges. The error estimate is established on each edge whatever the dimension is. It enables to calculate a stretching factor providing a new edge length distribution, its associated tensor and the corresponding metric. The optimal stretching factor field is obtained by solving an optimization problem under the constraint of a fixed number of edges in the mesh. Several examples of interpolation error are proposed as well as preliminary results of anisotropic adaptation for interface and free surface problem using a level set method.     

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

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

A two-dimensional unstructured cell-centered multi-material ALE scheme using VOF interface reconstruction

We present a new cell-centered multi-material arbitrary Lagrangian–Eulerian (ALE) scheme to solve the compressible gas dynamics equations on two-dimensional unstructured grid. Our ALE method is of the explicit time-marching Lagrange plus remap type. Namely, it involves the following three phases: a Lagrangian phase wherein the flow is advanced using a cell-centered scheme; a rezone phase in which the nodes of the computational grid are moved to more optimal positions; a cell-centered remap phase which consists of interpolating conservatively the Lagrangian solution onto the rezoned grid. The multi-material modeling utilizes either concentration equations for miscible fluids or the Volume Of Fluid (VOF) capability with interface reconstruction for immiscible fluids. The main original feature of this ALE scheme lies in the introduction of a new mesh relaxation procedure which keeps the rezoned grid as close as possible to the Lagrangian one. In this formalism, the rezoned grid is defined as a convex combination between the Lagrangian grid and the grid resulting from condition number smoothing. This convex combination is constructed through the use of a scalar parameter which is a scalar function of the invariants of the Cauchy–Green tensor over the Lagrangian phase. Regarding the cell-centered remap phase, we employ two classical methods based on a partition of the rezoned cell in terms of its overlap with the Lagrangian cells. The first one is a simplified swept face-based method whereas the second one is a cell-intersection-based method. Our multi-material ALE methodology is assessed through several demanding two-dimensional tests. The corresponding numerical results provide a clear evidence of the robustness and the accuracy of this new scheme.     

Simulation of Vortex Convection in a Compressible Viscous Flow with Dynamic Mesh Adaptation

In this work, vortex convection is simulated using a dynamic mesh adaptation procedure. In each adaptation period, the mesh is refined in the regions where the phenomena evolve and is coarsened in the regions where the phenomena deviate since the last adaptation. A simple indicator of mesh adaptation that accounts for the solution progression is defined. The generation of dynamic adaptive meshes is based on multilevel refinement/coarsening. The efficiency and accuracy of the present procedure are validated by simulating vortex convection in a uniform flow. Two unsteady compressible turbulent flows involving blade-vortex interactions are investigated to demonstrate further the applicability of the procedure. Computed results agree well with the published experimental data or numerical results.     

A Lagrangian integration point finite element method for large deformation modeling of viscoelastic geomaterials

We review the methods available for large deformation simulations of geomaterials before presenting a Lagrangian integration point finite element method designed specifically to tackle this problem. In our Ellipsis code, the problem domain is represented by an Eulerian mesh and an embedded set of Lagrangian integration points or particles. Unknown variables are computed at the mesh nodes and the Lagrangian particles carry history variables during the deformation process. This method is ideally suited to model fluid-like behavior of continuum solids which are frequently encountered in geological contexts. We present benchmark examples taken from the geomechanics area.