弹塑性流体力学拉氏自适应网格和网格重分守恒重映数值模拟方法

拉氏自适应重分弹塑性流体力学有限元程序实现了网格完全自适应，具有良好、灵活的非结构自适应网格数据结构，实现了滑移界面两边(接触间断)网格动态调整，网格的细分和合并处理灵活，网格重分和网格自适应模块兼容、守恒重映，网格重分中采用多种方法控制新网格的质量，爆轰计算可采用Lee-Tarver的化学反应率模式。初步数值计算结果表明，弹塑性流体力学拉氏自适应重分数值模拟方法合理，计算结果正确，基本反映了流场的物理结构。     

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

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

散乱物理量逼近的插值重映算法

随着流体问题的复杂和数值模拟高精度的要求,常常需要不同方法的耦合或不同网格体系的转化、叠加.只要改动网格或不同体系耦合(网格重构、细化与粗化、结构网格与非结构网格的互换)就涉及两种体系网格物理量的耦合或重映.将曲面拟合或插值的方法引入到散乱物理量重映计算中,给出几种物理量重映的算法.这些算法不受网格体系的限制,简单易用.     

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

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

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

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

散乱物理量逼近的重映算法

在计算流体力学(CFD)领域，几乎所有的方法都离不开网格，网格是各种数值方法的基础。网格质量的好坏直接影响数值结果的精度，甚至影响到数值计算的成败。为此CFD工作者发展了许多方法。如迭合网格、贴体网格和非结构网格，为了更好地数值模拟大变形问题，又进一步发展了结构/非结构混合网格的技术，尤其是发展了网格跟随流场智能化调整的网格自适应技术。这些网格技术的发展，几乎都涉及网格的变动，只要改动网格就涉及物理量的重映。重映方法一般被分为两种类型：插植重映和积分重映。所谓插植重映方法就是在计算区域D上，用已知网格上的物理量分布，通过插值理论把它插值到新网格或任意定义的规则网格上的一个过程，通过这个过程给出新网格的物理量的分布。所谓积分重映方法就是用积分的形式把旧网格上的守恒量重新映射到新网格上。如对某种体积密度分布q，简单的积分形式为     

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

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

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

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

二维拉氏自适应流体动力学（LAD2D）软件研制

二维拉氏自适应流体动力学软件LAD2D，是采用建立在拉氏自适应结构和非结构网格上的有限体积格式，可以计算平面二维和柱对称二维多物质大变形弹塑性流体动力学问题。LAD2D软件系统主要由5部分组成：主控程序、数据模块、前处理模块、主体计算模块、网格模块和后处理模块。其中主体计算采用了结构网格与非结构网格联合使用的拉氏网格体系，计算格式采用了有限体积格式。网格模块包括网格生成、自适应网格加密（AMR技术）和网格重分技术，以及网格改变后物理量守恒重映技术。LAD2D软件系统由主体程序、二维网格生成程序（GRID2D）、二维自适应网格加密程序（AMR2D）、二维自适应程序（ADAPT2D）f[I-维物理量重映程序（REMAP2D）组成。     

一种高精度保界的守恒重映方法

Lagrange方法中，当流场发生大变形时，跟踪流体运动的Lagrange网格发生扭曲，使计算无法进行下去，此时必须重分网格，把网格修复成较好的形状。另外，网格自适应技术中的重构、合并与加密，以及同一问题不同程序相继计算的连接，并行计算中相邻块边界区域的数据传递等，这些情况都需要利用旧网格上的物理量来确定新网格上的物理量，是一个物理量重映过程。质点重映方法是基于物理上守恒规律的一种离散的物理量守恒映射方法，既可实现分片常数分布的一阶精度重映计算，又可实现分片线性分布的二阶精度重映计算。这种方法可严格保证守恒量的守恒性，且可以实现任意多边形网格以及节点上物理量的守恒重映。但是，基于分片线性分布的二阶精度重映方法，如果新网格的守恒量没有进行保界调整，那么相应的强度量有可能在其局部的限制范围之外，破坏了原网格物理量的单调性。因而，对二阶精度的质点重映方法进行了进一步研究。在分片线性分布的基础上，将基于结构网格的保界算法扩展到非结构网格上，给出了二阶保界的质点守恒重映方法。     

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

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

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

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

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

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

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.     

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

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