激光微弯曲成形机理的数值研究

利用自行研制的含热传导、冲击动力学大、变形有限元程序，模拟了小尺寸梁在脉冲激光加热条件下的变形过程。在此基础上，利用商用程序模拟了冷却及残余应力的产生，研究了激光参数（强度及分布）等对于微弯曲的影响。数值模拟结果与文献中的实验观察相吻合。

快速凝固/粉末冶金高阻尼铝合金挤压成形过程的数值模拟

采用有限元分析软件ANSYS,分析了快速凝固/粉末冶金(Rapidly Solidified/Powder Metallurgy,RS/PM)高阻尼铝合金FMS0714/10(Zn-30Al)的挤压成形过程,研究了模具与坯料间的摩擦条件和挤压比对阻尼铝合金成形过程的影响,探讨了挤压过程材料表面产生裂纹的机理,进行了成形工艺优化.研究结果表明:挤压过程中FMS0714/10(Zn-30Al)阻尼铝合金材料的流动与变形不均匀,导致了应力应变分布的不均匀:材料表面应力应变较大,芯部应力应变较小;表面过大的应变和轴向拉应力是表面裂纹的诱因;减小模具和坯料间的摩擦,增大挤压比,可以减小材料表面应力,使应变分布趋于均匀,从而减少材料表面损伤,优化材料表面质量,提高成品合格率.数值模拟的研究结果将为FMS0714/10(Zn-30Al)阻尼铝合金材料挤压工艺的制订以及新材料的设计和研制提供有益的参考.

气相爆轰波数值模拟中化学反应模型研究

基于改进的时-空守恒元解元算法对气相爆轰波数值模拟中3种常用化学反应模型(二步模型,基元反应模型和Sichel的二步模型)进行了考察。对平面爆轰波和具有胞格结构的爆轰波进行了数值模拟,并对数值结果进行了比较和讨论。结果表明:3种化学反应模型得到的爆轰参数准确性有所差异,但得到的胞格结构均能和实验结果较好吻合。3种化学反应模型在爆轰波数值模拟中各有优缺点,应视具体问题决定使用哪种化学反应模型。

氢氧混合气体爆轰波的真实化学反应模型数值模拟

采用高精度的ENO格式和基于基元化学反应的真实化学反应模型求解氢氧混合气体一维爆轰波的精细结构。采用直接起爆方法得到稳定传播的爆轰波，计算的爆轰波阵面参数和实验相当符合。对爆轰波反应区化学反应的研究表明，参与反应的不同组分具有不同类型的变化特征。网格尺寸影响的研究表明，计算结果的精度随着网格尺寸的增加而增加，并能保持较好的收敛性。移动网格研究结果表明，网格运动速度和爆轰速度接近时，两者的相互作用对计算结果产生一定影响。

水中悬浮隧道在冲击载荷作用下的计算模型与数值模拟

建立了水中悬浮隧道在冲击载荷作用下的简化计算模型。用等效质量法将圆柱壳分布质量折算成冲击点处的集中质量,模型中考虑流体附加质量和系统阻尼的影响。根据碰撞过程中的动量守恒、变形过程中的能量守恒以及结构的位移与内力关系,得到问题的解析解。为验证解析解,在ANSYS/LS-DYNA中建立了动态冲击有限元分析模型。通过算例分别考察了在忽略和考虑流体附加质量两种情况下,冲击点位置和冲击速度对冲击点处最大径向位移的影响,将解析解与数值解进行对比,结果吻合较好。然后采用数值模拟方法得到了系统阻尼对计算结果的影响规律。数值模拟过程中还可以得到冲击点处的最大Mises应力。

突扩圆管内液-固两相流固体颗粒运动特性的DPM数值模拟

采用考虑颗粒碰撞的欧拉-拉格朗日数值模拟方法(DPM),对水平突扩圆管中液固两相流固体颗粒的碰撞过程进行了数值计算。在模型中,对液相采用欧拉法建立控制方程,对离散颗粒采用拉格朗日方法模拟。采用硬球模型描述颗粒间的碰撞作用。计算结果表明,该模型可以真实地模拟液固两相流中固体颗粒运动的动态变化过程以及颗粒的非均匀分布特征,从单颗粒层次上提供颗粒的运动信息,这有助于深入研究液固两相流中固体颗粒的运动规律

外部流场对激光加热运动目标影响的数值模拟

考虑激光与运动目标相互作用的基础上。利用有限元方法分析了亚声速条件下运动目标在激光辐照全过程的温度场和热应力场的分布与演化规律。结果表明：高速流场的存在，导致了明显的冷却效应；加热过程中目标材料出现了屈服，导致激光熄灭后结构内出现残余应力和变形；激光辐照区边缘产生很高的温度梯度和应力梯度，并且由于气流影响，受辐照区域前后两端应力分布不对称。

微重力环境下蒸发液层热毛细对流的数值模拟

提出了一种新模型来研究由单一物质构成的液层在其纯蒸气中的蒸发．液层置于微重力环境中并且受到水平方向温度梯度的作用，液层的热毛细对流和蒸发耦合在一起，使得气液界面的传热传质规律更加复杂．用理论分析的方法求解了不考虑热毛细效应的纯蒸发模型，得出温度场分布和界面质量流量的解析表达式．对于热毛细对流和蒸发耦合情况，采用有限差分的投影算法同时求解Navier-Stokes方程和能量方程，得到了不同蒸发Blot数和Marangoni数下流场和温度场的稳态数值解．论述了蒸发Biot数和Marangoni数对界面传热传质的影响，提出并解释了蒸发和热毛细对流耦合的三种模式．

无缓冲气COIL扩压器流场数值模拟

以无缓冲气化学氧碘激光器（COIL）实验器件的数据计算得到的混合喷管出口参数平均值作为光腔入口条件，对6种不同构型的扩压器从光腔入口至扩压器出口的流场进行了数值模拟，得出了各流场参数分布；对不同构型扩压器的流场特点、总压恢复性能进行了分析；研究了扩压器出口背压对流场参数的影响。结果表明：对于主流无缓冲气的COIL，等截面扩压器具有较好的压力恢复性能；增大扩压器出口背压可以使扩压器的压力恢复性能提高，然而，较高的背压使激波串向光腔方向移动，从而使光腔流场受到干扰，影响光腔的光束质量。

直径任意分布球填充的数值模拟

提出球填充数值算法的新分类方法.改进原有的松弛算法,使其能够模拟直径任意分布的球填充问题,采用可变循环周期使不同球数情形下的填充率基本保持不变.算例数据表明,该算法的填充率和配位数均高于原算法.由于采用背景网格搜索和双向链表组数据结构,使得邻接球搜索效率有相当大的提高,算法的时间复杂度为O（N）（N为球数）.在一台AMD Athlon 3200＋PC上,对于10000个等径球的随机密排列,只需217s,填充率即可达到0.64.

注蒸汽井套管应力的数值计算与工业应用

在注蒸汽吞吐井的井筒应力模型中引入井筒温度场，利用ANSYS软件对不同约束条件下井简的应力场进行了计算．结果表明，最大热应力发生在套管内壁，超过了N80套管的弹性屈服极限，最大热膨胀发生在温度过渡区．当套管周围掏空时，热应变远大于材料的弹性极限应变，是诱发热采井套管变形损坏的主要原因．为此提出了相应的套管损坏防治措施，在油田现场采用厚壁合金钢TP120TH套管完井，试验结果验证了防治措施的有效性．

气相爆轰波传播特性的数值模拟及实验对照

本文应用基元反应模型和频散可控耗散格式(DCD)对氢氧爆轰波进行了二维数值模拟.氢氧混合物的化学反应模型考虑了8种组分20个反应方程式.在处理化学反应引起的刚性问题时采用了时间算子分裂的方法.本文首先对爆轰波数值结果和实验结果进行了对照验证,然后对爆轰波在楔面反射由马赫反射向规则反射转变的过程进行了数值分析,得到了反射转变临界角,并和实验结果及理论分析进行了比较,结果是令人满意的;本文还对爆轰波的多波结构进行了初步的数值分析.

磨料流在离心研磨过程中运动特征的数值模拟

用液固两相流体动力学的方法对离心研磨过程中磨料流的运动特征进行数值模拟，得出磨料流运动的动态结果，并与高速摄影实验结果进行对比，指出计算机数值模拟结果对研磨机的设计和研磨工艺参数的选择具有实际的指导意义。该方法经济简捷，可以代替大量的工艺实验实现对离心研磨技术参数的优化，从而提高生产效率。

微重力蜡烛火焰特征数值模拟

建立了微重力蜡烛火焰的数学模型。计算与分析表明，火焰的形状由空气动力学特征决定，火焰是温度取决于化学反应动力学特征和火焰的热损失。在静止微重力环境中，自然对流的消失使火焰为半球形。辐射热损失对烛火焰温度（颜色）特征的形成有重要贡献，在静止重力环境下，化学反应放热速率受氧气扩散速率控制，辐射热损失的冷却使火焰温度低于正常重力温度值。但当环境气体的流动速度加大时，辐射热损失的影响逐渐减小，蜡烛火焰的温度逐渐接近正常重力蜡烛火焰的温度。当氧浓度较小时，火焰峰值温度小于烟黑形成的阈值温度（1300 K）；当氧浓度较大时，火焰温度大于烟黑形成的阈值温度。