螺旋槽干气密封微尺度流动场的近似计算及其参数优化

应用PH线性化方法、迭代法,近似求解了螺旋槽内稳态微尺度流动场的非线性雷诺方程,求得了气体动压和速度分布的解析解。继而利用多目标优化方法构建了气膜刚度与泄漏量之比的协调函数,并对该目标函数进行了近似求解,获得了最佳的螺旋槽几何参数值。     

非对称射流形成的几何理论(二)

将平板非对称正碰撞形成射流满足的定常方程组封闭条件推广到非正碰情况，提出了平面二维条件下非对称射流形成的一般理论，给出了射流形成的临界条件和出流射流的理论公式。对斜碰撞情形下非对称射流的出流方向、厚度随碰撞角和来流初始位形的变化规律进行了理论预测和数值模拟，理论解与二维欧拉程序的仿真结果符合较好。     

液态金属熔池中气泡-液体两相湍流的数值模拟研究

应用欧拉-欧拉双流体模型,以及分别描述气泡和液体湍流的两相湍流模型,描述液态金属中气体射流发泡过程的两相流动及气泡分布。与实验结果的对比表明,本模型的预测能力优于前人未考虑气相湍流的模型。考查了液体物性和气泡尺寸对流场特性的影响,结果发现,液体密度越低,气泡尺寸越小,气泡分布就越均匀。     

人工冻结黏土承载力模型试验及数值模拟分析

在天然冻土地区进行工程建设需研究冻土地基承载力。在人工冻结黏土蠕变试验的基础上建立广义开尔文模型,通过改造人工冻土大型冻结试验平台,进行人工冻结黏土地基承载力试验,发现载荷–沉降曲线可以分为3个阶段,基础呈局部剪切破坏。在ANSYS中将考虑相变的温度场转化为温度载荷施加在三维实体模型上,并将广义开尔文模型应用于软件中,采用增量加载方式分级加载。对比数值模拟与模型试验结果,发现地基载荷–沉降曲线规律一致,极限承载力大小基本相等,基础破坏模式相同。结果表明ANSYS软件可以分析复杂边界、初始条件下冻土材料的性质,且能得到不同载荷下冻土材料的应力、位移、塑性区发展状况。     

Bingham流体在旋转圆盘上流动的数值解

导出了描述Ｂｉｎｇｈａｍ流体在旋转圆盘上流动的基本方程，用差分方法数值解薄膜厚度分布方程，得到二种类型的厚度分布。数值解分别和计算机磁盘的厚度分布，Ｊｅｎｅｋｈｅ等的实验结果定性一致。     

炸药驱动强电流短路开关的研

由于电磁轨道炮和其它强电流装置中短路运行的需要,研制了一种条形炸药驱动的大电流闭合开关。利用爆炸成型技术使金属飞板负电极在炸药部件作用下,大面积地同时切断绝缘薄膜,嵌入另一电极模槽中构成牢固的接触。开关的绝缘强度大于30kV,电感(包括传输线)小于50nH,最大峰值电流达到400kA,闭合接触电阻小于0.1m,动作时间约为14us。对于开关的爆轰,撞击及变形过程进行了一维和二维数值模拟计算,并提出了一个描述爆炸开关行为的经验模型。利用这个模型进行电磁轨道炮整个系统的数值模拟,同实验结果符合较好。     

复杂结构优化选型的计算方法

本文为复杂结构优化选型提供了一个比较实用的计算方法:根据原有结构型式和设计经验来建立优化选型计算模型;对目标函数和约束条件作数学上近似处理后,利用Kuhn-Tucker条件和Lagrange乘子,把选型问题转化为求解典型二次规划问题;引进加权因子和利用改进后的变尺度法使计算过程简化、可靠,而且精度较高.     

爆炸容器内冲击波系演化及壳体响应的数值研究

对中心装药爆炸后冲击波的产生、传播和壳体动态响应全过程进行了数值研究。认为RDX瞬时爆炸 ,爆炸近场采用自相似解 ;冲击波传播和波系演化采用PPM (the Piecewise Parabolic Method)格式求解Eu ler方程 ;壳体响应采用有限元方法求解拉氏坐标系下由虚功原理得到的动力学方程。壳体内壁面边界条件分别采用强耦合和弱耦合方法处理。结果表明 :(1)当装药量相同时 ,薄壁壳体振型比厚壁壳体复杂得多 ,振幅也大 ;(2 )当装药量不同 ,壳体厚度相同时 ,爆炸场冲击波的演化过程不同 ;(3)对少量装药 ,产生的冲击波强度低 ,壳体变形小 ,是否考虑内边界运动 ,对计算结果的影响不大 ;(4 )在本文条件下 ,爆炸容器封头顶点所受的载荷最大 ,是最易发生破坏的地方 ,侧壁与爆点所在横截面的交线 ,也易破坏。     

后台阶流动的数值模拟

访述了大涡模拟的基本思想，指出大涡模拟的效率主要取决于四个因素，即流动中须有大尺度涡存在、合理的计算格式、合适的滤波器和亚格子应力模型。在深入考虑粘性不可压缩流Navier—Stokes方程各个子项作用的基础上，提出二阶全展开Euler—Taylor—Galerkin有限元方法作为大涡模拟的离散格式，并采用Gauss滤波器，对典型算例——后台阶处的流动进行大涡模拟，计算结果与相关文献符合的很好。从计算结果还可以看出大涡模拟与二阶全展开ETG有限元方法的结合在捕捉涡系及反映涡动时变过程方面具有明显的优势，说明大涡模拟适合于边界几何形状复杂区域流动的模拟。同时应用二阶全展开ETG有限元方法对低雷诺数粘性不可压缩后台阶流动进行了计算，得到与相关文献符合良好的计算结果，即该方法也可独立用于对低雷诺数粘性不可压缩流动的计算。     

湍流加速火焰的三维数值模拟

火焰在设有障碍物的管内传播时会自身加速，并可能导致爆炸。本文基于湍流κ-ε模型和改进的EBU—Arrhenius反应模型，对该现象进行了三维空间的数值模拟。计算结果反映了障碍物、湍流和火焰之间相互作用的正反馈机理，描绘了火焰在管内加速传播的三维图像。