二维等离子体电弧随时间演化的数值模拟

采用交错均匀网格、全隐欧拉后差格式,对二维轴对称等离子体电弧进行随时间演化的数值模拟.得到了平衡时的稳态温度场、电流场和随时间演化的气流速度场.发现其中有周期性传播的磁流体激波,激波速度各点不等,约为700～1000m/s,这与实验观测到的等离子体电弧不稳定性现象是一致的 关键词：     

沙漏计时原理二维数值模拟

提出沙漏计时二维数值模型,并用元胞自动机的方法对该模型进行计算机模拟.结果表明由于开口处沙粒的流动而引起的沙堆崩塌行为是一个准周期性振荡过程,开口大小与崩塌总时间近似成线性关系.特殊的开口位置会引起崩塌"拍",且随沙堆高度呈周期性变化.     

爆轰产物驱动飞片运动数值模拟

化爆加载装置是研究材料动态响应特性的重要冲击加载装置之一。如果炸药与飞片直接贴附，将难以获得较低的飞片速度，而且飞片中载荷的上升时间较短。为研究材料在冲击压缩低压状态和准等熵压缩条件下的动态响应特性，对利用爆轰产物经过空气隙驱动飞片运动的加载装置进行了二维数值模拟，获得了一些规律性认识。     

二维二相不规则网格的油藏数值模拟

本文采用不规则网格的差分方法,对平面油水模型进行了数值模拟,并用IMPES方法计算了两个实例。例1的计算结果与精确解吻合良好,比Pedrosa的局部加密网格法使用的网格数少而且计算更简便。例2为矩形网格难以计算的问题,应用本文方法也得出了较好的结果。这表明本文的不规则网格差分方法是可靠而简便的,它可以有效地应用于复杂外边界和复杂地层结构的油藏模拟问题。     

靶球压缩二维数值模拟比对研究

利用轴对称二维数值模拟程序Ｌａｒｅｄ－Ⅰ,对两端单环激光入射产生的辐射场驱动下的靶球进行了数值模拟,给出了靶球二维变形的计算结果。随着入射激光焦斑位置的改变,辐射场在靶球表面的温度分布随之相应改变,驱动靶球内爆压缩氘气区分别成为铁饼形、香肠形或球形。将数值模拟结果与美国ＬＬＮＬ实验室的实验结果及其Ｌａｓｎｅｘ程序的数值模拟结果进行了比对,得到了定性上完全相同的结论和接近的数值结果。从数值模拟结果可以看出,在单环入射的情况下,尚无法消除辐射场的高阶不对称性,无法满足球形对称均匀压缩的要求。     

低温离心泵的二维流场数值模拟及分析

低温离心泵是使用于低温工况下的离心泵,它具有尺寸小、转速高、升压高等特点。文中使用fluent软件,对自主设计的低温泵内流体的流动情况进行了数值模拟。通过划分计算区域、选择合适的计算网格、采用合理的湍流模型和边界条件等方面的处理,得到了一种模拟值与设计值比较吻合的数值模拟方法。利用数值模拟方法真实反映了低温泵内部的复杂流动,验证了泵的传统理论设计方法对于低温工况的适用性,并为低温泵内流道的设计、改进提供了理论依据。     

低速轴流压气机旋转失速的二维数值模拟

本文通过求解二维不可压N-S方程,对某三级低速轴流压气机的第一级进行数值模拟。首先用定常计算得到了该级的稳态性能曲线,然后在级出口加上节流阀进行非定常计算,模拟压气机进入失速的整个过程,重点是先兆的发展和内部流场的分析。计算结果表明,当阀门关到某个位置,无外加扰动,像数值误差这样的小扰动就能使压气机失速。本文还讨论了不同轴向计算域、关阀门速率等对模拟结果的影响。     

二维光子晶体中的无序数值模拟研究

应用on-shell方法研究弱无序情况下的二维光子晶体，无序的引入使得原来不透电磁波的禁带中，出现了共振模，即某一波长的电磁波可以透过，计算结果表明共振模的出现位置，条数多少及透过率随着无序程度的不同而变化。     

一维爆轰波不稳定性的数值模拟

 对当量比为1的乙烷与空气的混合气体的一维爆轰不稳定问题进行数值模拟，得到了不同大小的网格对爆轰不稳定问题数值模拟结果的影响。网格大小Δx由从ZND模型分析得到的导引长度L_in确定，网格大小从0.2L_in变到0.002 5L_in。随着网格变细，没有得到振幅趋于一致的解，每一种网格尺寸得到的解的振幅都互不相同。当网格大小为Δx=0.01L_in、0.005L_in时，得到有规则的爆轰激波阵面压力的振荡，振荡的模式是峰值一大一小的振荡。网格更细时，爆轰波的振荡在计算范围内由一些有规则的振荡和一些较不规则的振荡组成。但爆轰激波阵面压力振荡的波长最后趋于一致，为91~93 mm，与实验得到的胞格长度88 mm很接近。     

双层多孔介质内预混燃烧的实验和二维数值模拟

对双层多孔介质燃烧器内丙烷/空气预混燃烧进行实验和数值研究.实验对多孔介质燃烧器内固体温度场分布进行测量分析;数值计算利用商业软件FLUENT6.2,通过添加用户自定义标量方程和用户自定义函数,对有壁面散热的双层多孔介质内预混燃烧进行了二维模拟,并与实验结果进行了比较.结果表明,双层多孔介质燃烧器具有良好的稳定燃烧范围和较低的污染物排放;壁面散热对多孔介质燃烧的影响不可忽视.     

JO-9159炸药的初始密度对爆轰波反应区宽度影响的实验研究

 利用电磁法对两种具有不同初始密度的JO-9159炸药的爆轰反应区宽度进行了测量。结果表明,对于JO-9159这种爆轰反应区非常窄的炸药,既使爆轰达到了定常状态,化学反应仍然与炸药的初始物理状态有关。     

滑移爆轰作用下内爆柱形钢壳层裂的数值模拟

 Иванов等发表了滑移爆轰作用下内爆柱形钢壳层裂实验的回收试件测量结果，也给出了相应的数值模拟结果，但两者的差别较大。本文利用损伤度函数方法，对上述实验做了数值研究，所得到的回收试件尺寸数值结果与Иванов实验结果十分相符。这说明，损伤度函数方法对研究动态断裂问题有一定的普适性。本文还对Иванов实验的柱形钢壳内爆运动过程进行了数值模拟，讨论了临界损伤度α_c取值对断裂面位置、钢层运动轨迹及钢壳运动速度历史的影响，指出了在分析界面不稳定性问题时要注意α_c取值的影响。     

激光与靶非平衡耦合的数值模拟

对激光与高Z靶耦合主要物理过程,其中包括激光吸收,电子热传导,等离子体运动,非平衡原子物理,X光转换和辐射输运等现象的数值模拟作了全面的概括性描述。提出了总体微分方程组,建立了差分计算格式和求解方法,给出了程序数值计算的结果。     

非平衡Al等离子体三温自由程的数值模拟

在非平衡系统中,采用三温(光子温度t_γ,电子温度t_e和离子温度t_i)模型来计算光子能流的自由程;通过解速率方程得到电子的占据几率P_i,讨论了它与用Femi-Dirac分布函数算得的n_f之间的差别;给出了Al等离子体三温自由程的数值模拟结果。     

喷动流化床的数值模拟研究

本文借鉴 Ｃｕｎｄａｌｌ＆ Ｓｔｒａｃｋ提出的碰撞模型,将欧拉方法和拉格朗日法结合起来,建立了喷动流化床气固两相流动的数学模型,在计算机上对试验台进行了模拟,探讨了床内颗粒的运动特征。按１：１的比例建立了喷动流化床的模型试验台并进行了试验,对模拟结果和试验结果进行了比较与分析。颗粒在喷动床内的运动特征,模拟结果都能够再现出来,模拟结果与试验结果取得了很好的一致性。     

电子束辐照效应的数值模拟

论述了电子束辐照材料时的能量沉积Monte-Carlo计算,在此基础上,与弹塑性流体力学方程组相结合,给出电子束激发的热一力学效应的数值模拟方法.对一系列电子束打铝靶实验进行了实际计算,得到与实验比较一致的结果.     

单个球形生物质颗粒热解过程的数值模拟

本文对单个球形生物质颗粒的热解过程进行数值模拟．利用双倒易边界元法和四阶龙格-库塔方法分别对非线性的导热方程和化学反应动力学方程组进行求解．讨论了生物质颗粒的大小与环境温度对热解时间和热解产物中相关组分质量分数的影响．     

平面叶栅中的湿蒸汽两相结流动数值模拟

对存在自发凝结的湿蒸汽两相流动建立了完全欧拉坐标系下的数值模型.考虑水蒸汽与理想气体的偏差,引入维里型状态方程对模型作了进一步完善.对平面叶栅中的两相凝结流动进行了数值模拟,计算结果与实验值比较表明本文计算模型正确,可以扩展到三维复杂两相凝结流场的计算.     

具有隔板的平行通道内空气混合对流换热数值模拟

本文以锅炉干排渣装置为背景,对抽象的理论模型具有隔板的平行通道内空气混合对流换热进行了数值模拟.数值计算表明,在Re＞1000时应采用非稳态数学模型进行数值模拟;在Re＞500时,自然对流机制对流动和换热的影响基本可以忽略.数值计算给出了不同Re时的进出口无量纲压差、局部的Nux和平均Nu以及流线图.这些结果可为深入研究干排渣装置中流动和换热特性提供参考.     

横向电动效应作用下微混合的数值模拟

微流体机械系统中电渗流的雷诺数一般都比较低,微通道内的混合主要通过扩散来完成,因此需要较长的距离和时间才能混合均匀.本文在一个二维微通道内施加横向电场,并讨论了不同横向电场强度、微电极对布置情况、壁面zeta电势等因素对微混合的影响.结果表明此方法可以有效地提高混合效率,增加横向电场强度和电极对密度可以明显提高混合效率,较小的zeta电势可以获得更好的混合效果.