用浮阻力模型研究Richtmyer-Meshkov不稳定性诱导混合

采用浮阻力模型对激波管低压缩和激光加载高压缩情况下的Richtmyer-Meshkov不稳定性诱导混合现象进行了研究。通过与实验和理论分析结果进行比较发现:为了达到好的吻合, Richtmyer-Meshkov不稳定性情况下阻力系数的取值范围(2.0~5.36)比Rayleigh-Taylor不稳定性情况下的值(3.3~4.0)宽得多; 而在Richtmyer-Meshkov不稳定性情况下, 高压缩时阻力系数的不确定度(约为3.36)明显高于低压缩时的值(约为1.46), 模型的进一步完善还有待于更精确实验的验证。研究显示:指数律经验公式中指数随工况的不同而显著变化, 目前工程设计中采用指数律经验公式是粗糙的。     

欧拉流体力学数值方法中的界面计算

本文所提算法适用于二维和三维多介质流体力学两步欧拉数值方法中输运计算的混合网格(包括自由面网格)界面处理。在一个混合网格中，界面被近似地看作直线(二维)或平面(三维)。整个方法分为三步：(1)第一步，用混合网格周围的八个网格的介质面积份额(二维)或二十六个网格的介质体积份额(三维)确定界面的法线方向；第二步，用混合网格的本身的介质面积份额(二维)或体积份额(三维)确定界面的方程(位置)；第三步，用此直线方程求出通过网格边界的流以及下一时刻网格的面积份额(二维)或体积份额(三维)。最后给出了用此方法所做的一些数值计算及与SLIC算法的比较。     

圆筒试验JWL状态方程参数的贝叶斯标定

在研究数值模拟的输入参数引入的不确定性时，通常需要人为给定每个输入参数的概率分布，且输入参数概率分布的选择可能会对分析结果产生直接影响。利用贝叶斯方法标定了圆筒试验JWL状态方程参数，得到了标定参数的估计值和后验分布，并研究了不同统计模型假设对标定参数的估计值和后验分布的影响。贝叶斯后验分布融合了基准试验的试验数据的信息，因此将其作为不确定度量化分析时输入参数的初始概率分布，可以尽量减少分布选择引入的认知不确定性。     

基于多项式混沌的全局敏感度分析

回顾基于多项式混沌和方差分解的全局敏感度分析方法,针对高维数随机空间和高阶多项式混沌展开面临的"维数灾难"问题,采用回归法、稀疏网格积分及基于l1优化的稀疏重构技术(即压缩感知技术)来减少非嵌入式多项式混沌方法所需的样本配置点数目.针对几个典型响应面模型(包括Ishigami函数、Sobol函数、Corner peak函数和Morris函数)进行Sobol全局敏感度指标计算,展示多项式混沌方法在基于方差分解的全局敏感度分析中的有效性.     

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

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

拉氏数值模拟中的人工粘性方法

较好的人工粘性需要满足较小的计算开销、不能去除真实具有的涡运动等条件.提出一种应用于拉氏数值模拟中基于Lew人工粘性,同时增加了限制器的人工牯性方法.可以有效减少数值模拟结果对网格的依赖;采用特征值限制器控制施加的人工粘性大小,通过限制器能够区分激波压缩和等熵压缩;方便应用在二维、三维,结构网格或者非结构网格上.     

含强间断导热系数的三维扩散方程数值计算

针对惯性约束聚变领域需要求解的含强间断导热系数的扩散方程,借鉴各向异性导热的思想,给出了一种导热系数定义在网格节点的支撑算子改进格式。改进后的方法保留了支撑算子格式的算子相容性优点。数值算例表明,新方法可以很好地处理含强间断和强非线性的热扩散问题,计算结果保持较高的精度,适合于辐射流体力学问题的模拟。     

ALE中的重映算法

介绍了二维非结构网格上的守恒重映算法，重点是基于SFB／DC思想的通量重映算法。用统一的公式表示不同的单元量重映算法，包括原始的贡献网格法、Barth—Jespersen方法、最小二乘法，不同算法间的区别体现为梯度求法的差异。对于交错网格上速度的重映，介绍了SALE和HIS算法。此外，为保证重映算法的有界性，引入了修补方法。     

SALE速度重映算法的改进

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

缺陷形状及角度对于铝材料微喷射的影响

目前，一般认为引起微喷射的主要物理机制有：(1)试件的表面缺陷(如凹坑、划痕等)形成的微射流机制，很可能是微喷射的主要来源。(2)在冲击波作用下，晶界处产生局部能量沉积，而形成“热点”，致使物质局部熔化甚至汽化。(3)材料表面的杂质可能产生微层裂或剥落而形成微喷射。在材料发生微喷射的过程中，这些机理经常同时存在，对于不同的材料，在不同的条件下，可能是不同的机制起主导作用。在实验中发现，铝试件在冲击压为40GPa时，其表面温度大约为500K，这大大低于931K的熔化温度。因此，可以认为铝材料在冲击压力低于40GPa时，没有熔化现象发生。故而对于纯铝试件，当冲击压力达到约30GPa时，铝材料表面缺陷形成的微射流是其微喷射的主要机制。