R-Adaptive Reconnection-Based Arbitrary Lagrangian Eulerian Method-R-ReALE

We present a new R-adaptive Arbitrary Lagrangian Eulerian (ALE) method, based on the reconnection-based ALE - ReALE methodology [5, 41, 42]. The main elements in a standard ReALE method are: an explicit Lagrangian phase on an arbitrary polygonal (in 2D) mesh, followed by a rezoning phase in which a new grid is defined, and a remapping phase in which the Lagrangian solution is transferred onto the new grid. The rezoned mesh is smoothed by using one or several steps toward centroidal Voronoi tessellation, but it is not adapted to the solution in any way. We present a new R-adaptive ReALE method (R-ReALE, where R stands for Relocation). The new method is based on the following design principles. First, a monitor function (or error indicator) based on Hessian of some flow parameter(s), is utilized. Second, the new algorithm uses the equidistribution principle with respect to the monitor function as criterion for defining an adaptive mesh. Third, centroidal Voronoi tessellation is used for the construction of the adaptive mesh. Fourth, we modify the raw monitor function (scale it to avoid extremely small and large cells and smooth it to create a smooth mesh), in order to utilize theoretical results related to centroidal Voronoi tessellation. In the R-ReALE method, the number of mesh cells is chosen at the beginning of the calculation and does not change with time, but the mesh is adapted according to the modified monitor function during the rezone stage at each time step. We present all details required for implementation of the new adaptive R-ReALE method and demonstrate its performance relative to standard ReALE method on a series of numerical examples.     

A high order ENO conservative Lagrangian type scheme for the compressible Euler equations

We develop a class of Lagrangian type schemes for solving the Euler equations of compressible gas dynamics both in the Cartesian and in the cylindrical coordinates. The schemes are based on high order essentially non-oscillatory (ENO) reconstruction. They are conservative for the density, momentum and total energy, can maintain formal high order accuracy both in space and time and can achieve at least uniformly second-order accuracy with moving and distorted Lagrangian meshes, are essentially non-oscillatory, and have no parameters to be tuned for individual test cases. One and two-dimensional numerical examples in the Cartesian and cylindrical coordinates are presented to demonstrate the performance of the schemes in terms of accuracy, resolution for discontinuities, and non-oscillatory properties.     

Simulation of the spread of a viscous fluid using a bidimensional shallow water model

In this paper we propose a numerical method to solve the Cauchy problem based on the viscous shallow water equations in an horizontally moving domain. More precisely, we are interested in a flooding and drying model, used to modelize the overflow of a river or the intrusion of a tsunami on ground. We use a nonconservative form of the two-dimensional shallow water equations, in eight velocity formulation and we build a numerical approximation, based on the Arbitrary lagrangian eulerian formulation, in order to compute the solution in the moving domain.     

多介质五方程简化模型及界面捕捉的人工压缩算法

根据两介质五方程简化模型的基本假设,发展了适用于任意多种介质的体积分数方程。为了捕捉多介质界面,将HLLC-HLLCM混合型数值通量的计算格式推广应用于二维平面和柱几何的多介质复杂流动问题,在高阶精度的数据重构过程中采用斜率修正型人工压缩方法ACM。通过一维、二维多介质黎曼问题算例测试,结果表明：发展的计算格式能够较好地分辨接触间断和激波,间断附近物理量无振荡;对于添加了初始扰动的激波问题,能够有效抑制激波数值不稳定性;使用二维柱球SOD问题和接触间断型黎曼问题检验计算格式对多介质复杂流动问题的适应性。     

金属圆管内爆轰波相互作用效应的数值模拟研究

采用动力学有限元程序,针对大变形带来的网格扭曲造成计算的锁死现象进行了ALE技术处理,结合高速摄影和X光照相实验结果,对金属圆管内爆轰波的相互作用效应进行数值模拟,给出了爆轰波对碰作用下金属圆管运动计算结果,复现了金属圆管对碰区鼓包现象,计算结果与实验吻合较好。可以预见,基于该方法能够初步分析影响金属圆管运动规律的诸多因素。     

From Precursors to Thin Films - Thermoanalytical Techniques in the Thin Film Technology

Processing thin films for advanced applications, for instance in electronics and optoelectronics, involves several steps starting from precursor synthesis and ending up with the devices. Especially when optimizing the first steps of this chain of processes, thermoanalytical techniques play an important role. The review will focus on the main chemical deposition methods (CVD, ALE, spray pyrolysis, sol-gel) giving selected examples of problem-solving by thermal analysis. The techniques discussed are TG, DTA/DSC, EGA and their combinations. High-temperature X-ray diffraction (HTXRD) is also a powerful tool for in situ studies of thin films. The examples are taken from solar cell, superconductor and flat panel electroluminescent display technologies. This revised version was published online in July 2006 with corrections to the Cover Date.     

两种人工机械心瓣启闭过程的ALE有限元分析*

本文采用任意拉格朗日-欧拉(ALE)有限单元法,将血液视为不可压缩粘性流体,同时将人工机械瓣简化为定轴转动刚体,建立了机械心瓣/血液耦合运动的二维计算模型。在此基础上。计算了顺向人工机械瓣(DDM)的开闭过程,并与St.Jude瓣(SJM)进行了对比。由本文研究可得以下主要结论:1.DDM瓣与SJM相比,开启较为迅速,关闭则较为柔和。2.DDM瓣的峰值回流量较之SJM瓣小。本文研究表明DDM瓣由于更好利用了天然心瓣的关闭机理,而具有较高的耐久性潜力。     

Unitive analysis schemes on problems of multiple moving boundaries with 3-D liquid-solid multiple nonlinear coupling for uplift of anchored liquid storage tanks

I.Intr0ducti0nAnchored1iquidstoragetanksareextensivelyusedinbusiness,nuclearpowerplants,spaceflightandpetrochemicalindustry.Theystoragepetrochemicalproductions0fover9Opercentintheworld.Theanchoredliquidstoragetanksarecontinuelydamagedduringearthquakestoha…     

A CELL-CENTERED ALE METHOD WITH HLLC-2D RIEMANN SOLVER IN 2D CYLINDRICAL GEOMETRY

This paper presents a second-order direct arbitrary Lagrangian Eulerian (ALE) method for compressible flow in two-dimensional cylindrical geometry.This algorithm has half-face fluxes and a nodal velocity solver,which can ensure the compatibility between edge fluxes and the nodal flow intrinsically.In two-dimensional cylindrical geometry,the control vol-ume scheme and the area-weighted scheme are used respectively,which are distinguished by the discretizations for the source term in the momentum equation.The two-dimensional second-order extensions of these schemes are constructed by employing the monotone up-wind scheme of conservation law (MUSCL) on unstructured meshes.Numerical results are provided to assess the robustness and accuracy of these new schemes.