Computational modelling of multi-material energetic materials and systems

ABSTRACT

In this paper, we present a systematic roadmap for developing a robust and parallel multi-material reactive hydrodynamic solver that integrates historically stable algorithms with new and current modern methods to solve explosive system design problems. The Ghost Fluid Method and Riemann solvers were used to enforce appropriate interface boundary conditions. Improved performance in terms of computational work and convergence properties was achieved by modifying a local node sorting strategy that decouples ghost nodes, allowing us to set material boundary conditions via an explicit procedure, removing the need to solve a coupled system of equations numerically. The locality and explicit nature of the node sorting concept allows for greater levels of parallelism and lower computational cost when populating ghost nodes. Non-linear numerical issues endemic to the use of real Equations of State in hydro-codes were resolved by using more thermodynamically consistent forms allowing us to accurately resolve large density gradients associated with high energy detonation problems at material interfaces. Pre-computed volume tables were implemented adding to the robustness of the solver base.     

各向异性复合材料尖劈和接头的奇性应力指数研究

提出了一个新的、基于位移的、求解三维尖劈端部奇性应力指数问题的非协调元特征分析法。该方法假定尖劈端部邻域内的位移场没有采用奇异变换技术，导出虚功方程的出发点不同于过去原有求解裂纹尖端近似场的有限元特征分析法，在有限元离散时采用的单元形式为非协调元。文中运用该方法给出了若干求解各向异性复合材料尖劈／接头端部奇性应力指数的算例。所有的计算结果表明，本文方法能够求解复杂尖劈／接头的全部奇性应力指数，使用的单元少而且精度高。     

轴对称几何下耦合MOF界面重构的多介质ALE方法

推导了轴对称几何下的MOF(Moment of Fluid)界面重构,将其与多介质ALE方法相耦合,形成MOFMMALE方法,并应用于多介质大变形流动问题的数值模拟研究。数值算例表明,耦合MOF界面重构的多介质ALE方法是求解多介质大变形流动问题的有效手段,并且具有很好的界面精度和分辨率。     

基于level set的Eulerian-Lagrangian耦合方法及其应用

提出了一种基于水平集的Eulerian-Lagrangian耦合方法,其中Lagrangian方法采用相容显式有限元拉氏方法,Eulerian方法采用基于近似Riemann解的有限体积Eulerian方法,多介质界面处理采用新的水平集和Ghost方法计算. 给出了若干数值算例,包括激波管问题以及金属和气体的运动界面及其大变形问题,并分别与精确解和相容显式有限元拉氏方法的计算结果进行了对比. 数值结果表明,该方法计算结果正确,精度较高,能够准确捕捉物质界面,适用于处理大变形问题.      

聚能射流三维数值模拟

研究三维多介质界面处理及数值模拟问题。采用拉格朗日法、欧拉法相结合的方法在矩形网格上离散差分基本方程组;在欧拉步中引入模糊方法处理界面,计算各输运量;编写了数值模拟程序,并对线型装药金属罩聚能射流模型进行模拟。证明模糊界面描述和模糊输运计算有效、可行。     

The Plane Wedge Problem Loaded at Its Apex – the Self-similarity Property and the Characteristic Vector

In the present paper the elastostatic problem of a generally anisotropic and angularly inhomogeneous plane wedge loaded at its apex by a concentrated force, is studied in linear elasticity. At first the self-similarity property is formulated and the stress field of the inhomogeneous anisotropic self-similar wedge problem, is deduced. The wedge is radially separated and the plane wedge problem is reformulated by the introduction of a characteristic vector. Furthermore, the angular distribution of the load is determined. The multi-material wedge problem in terms of a formulation based on the isotropic angularly inhomogeneous wedge, is confronted, and necessary conditions that ensure the self-similarity property, are found. Finally, the similar elastostatic wedge problems and the involution between stresses, are studied. Mathematics Subject Classifications (2000) 74B05, 74K30, 34B05, 51N15.     

On an extension of the d’Alembert solution to initial–boundary value problems in multi-layered,multi-material domains

This work introduces a method for the exact solution of initial–boundary value problems for linear, one-dimensional conservation laws in multi-layered, multi-material domains. The method is based on the geometry of the solutions of such conservation laws and represents an extension of the d’Alembert solution to initial–boundary value problems in multi-layered, multi-material domains.     

多组份计算中封闭性模型研究

在使用单速度多组份方法计算多介质混合网格中的物理量时,需要给出一种封闭模型来使控制方程封闭。本文分析了压力增量相等封闭性模型存在基础遭到破坏的原因,之后在保证流场守恒性的前提下构造可以保证假设基础的可解方程组对压力增量相等封闭性模型进行修正,给出了一种改进方法。使用改进的模型进行了爆轰驱动飞片问题的计算,计算结果表明,改进的模型能很好的追踪物质界面和冲击波位置,同时较好的抑制界面处压力、内能震荡。     

基于CE/SE方法的二维Euler型多物质流体弹塑性问题计算

将CE/SE方法推广到二维固体流体弹塑性问题的数值计算,同时结合杂交粒子水平集方法追踪物质界面和合适的边界条件,提出一套完整的二维Euler型流体弹塑性计算方案.通过长钨杆侵彻装甲钢实验的数值模拟,对方法的精度和有效性进行验证.     

A 2D Staggered Multi-Material ALE Code Using MOF Interface Reconstruction

Hydrocodes are necessary numerical tools in the fields of implosion and high-velocity impact, which often involve large deformations with changing-topology interfaces. It is very difficult for Lagrangian or Simplified Arbitrary Lagrangian-Eulerian (SALE) codes to tackle these kinds of large-deformation problems, so a staggered Multi-Material ALE (MMALE) code is developed in this paper, which is the explicit time-marching Lagrange plus remap type. We use the Moment Of Fluid (MOF) method to reconstruct the interfaces of multi-material cells and present an adaptive bisection method to search for the global minimum value of the nonlinear objective function. To keep the Lagrangian computations as long as possible, we develop a robust rezoning method named as Combined Rezoning Method (CRM) to generate the convex, smooth grids for the large-deformation domain. Regarding the staggered remap phase, we use two methods to remap the variables of Lagrangian mesh to the rezoned one. One is the first-order intersection-based remapping method that doesn't limit the distances between the rezoned and Lagrangian meshes, so it can be used in the applications of wide scope. The other one is the conservative second-order flux-based remapping method developed by Kucharika and Shashkov [22] that requires the rezoned element to locate in its adjacent old elements. Numerical results of triple point problem show that the result of first-order remapping method using ALE computations is gradually convergent to that of second-order remapping method using Eulerian computations with the decrease of rezoning, thereby telling us that MMALE computations should be performed as few as possible to reduce the errors of the interface reconstruction and the remapping. Numerical results provide a clear evidence of the robustness and the accuracy of this MMALE scheme, and that our MMALE code is powerful for the large-deformation problems.