FAE街道爆炸的三维数值模拟与面向对象编程

采用气液两相方程,对燃料空气炸药（FAE）气液两相爆轰的爆轰波发展过程进行数值模拟,以此结果作为计算的初值,用三维有限体积TVD方法和面向对象编程技术,模拟FAE爆炸产生的冲击波在街道中的传播规律以及对建筑物的毁伤效应,计算结果表明,三维有限体积TVD方法描述了FAE爆炸冲击波的反射,绕射,传播等物理现象,在计算机模型的构造和管理中,C＋＋和面向对象编程技术体现出极高的优越性,对连续边界的处理,采用了垂直扰动法,得到了清晰的物理图像。     

基于非结构化网格有限体积的LBM

提出了一种基于非结构化网格有限体积的LBM.采用Cell-vertex有限体积法离散控制方程.该方法在时间上采用伪、实二时间步,其中伪时间步采用向前差分,实时间步采用二阶向后差分方法;空间上采用edge-based通量计算方法,采用高阶TVD格式计算控制体边界通量.离散后的控制方程采用隐式迭代,控制变量采用五层二阶Runge-Kutta方法求解.二维同心圆环内圆柱间Couette流与顶盖驱动方腔流的数值结果显示该方法为一种有效求解不规则边界流体动力特性的实用工具.     

无结构三角形网格下运动界面追踪的WENO有限体积法

在处理运动界面追踪问题的流体体积函数(VOF)法的基础上,给出了一种无结构三角形网格下的高分辨率的运动界面捕捉方法.该方法采用高精度的加权本质无振荡(WENO)有限体积格式离散VOF函数的空间导数,采用三阶TVD Runge-kutta方法离散时间导数,采用Lax-Friedrichs通量作为数值流通量.数值试验结果表明,用该方法来进行旋转速度场和剪切速度场的运动界面追踪,可以得到与理论解非常一致的追踪结果.     

空心弹丸流场数值模拟与阻力特性

采用隐式有限体积TVD(总变差减小)数值格式,在超声速条件下数值模拟了某空心弹丸绕流流场。系统研究了流场波系结构与空气阻力特性,证实空心弹丸阻力系数比普通弹丸小得多,同时发现发射角度对于空心弹丸的速度降影响不大。     

基于四边形剖分的多尺度有限体积方法

在一类多尺度方法（heterogeneous multi-scale method,HMM）的基础上,结合有限体积方法（finitevolume method,FVM）,提出多尺度有限体积方法（HMM-FVM）,用于数值计算一类多尺度抛物型偏微分方程,同时给出相关误差估计.计算结果表明,多尺度有限体积方法比传统的有限元法、有限体积法有效,既节省计算量,又有较高的精度.     

轴对称二维氢氧爆轰驱动激波管的模拟

对由前向爆轰产生驱动气体并具有局部锥形收缩截面的轴对称激波管利用两相常比热完全气体的欧拉方程组和有限体积TVD格式在局部正交的网格上进行了二维数值模拟.将消除两种介质界面处数值振荡的新方法成功地推广到二维情况.分析了该激波管的力学特性.     

轴对称地维氢氧爆轰驱动激波管的模拟

对由前向爆轰产生驱动气体并具有局部锥形收缩截面的轴对称激波管利用两相常比热完全气体的欧位方程组和有限体积TVD格式在局部下在交的网格上进行了二维数值模拟。将消除两种介质界面处数值振荡的新方法成功地推广到二维情况。分析了该激波管力学特性。     

三维Euler方程的非结构网格非振荡有限体积法

考虑三维Euler方程,对非结构四面体网格给出了一类时空分离的非震荡有限体积方法。空间方面采用Roe通量进行离散,同时基于物理量守恒运用最小二乘和极值原理在每个四面体单元构造非振荡线性重构多项式,时间离散采用两步Runge-Kutta TVD方法。文末对三维球形爆炸问题进行数值模拟,得到了满意的结果。     

火箭喷流对平板冲击的数值模拟

该文针对欠膨胀喷流对平板冲击产生复杂激波结构的特点 ,采用高精度的TVD有限体积格式对其进行数值模拟 ,获得清晰的冲击激波结构 ,如 :相交激波、马赫盘、曲面激波及马赫反射波等 ,同时给出不同时刻平板上的压力分布规律。为了加速收敛 ,文中应用了时间推进多步方法     

HTR－10控制棒及第二停堆系统反应性当量计算

介绍了ＨＴＲ－１０高温气冷实验堆位于反射层中的控制棒和作为第二停堆系统的硼球孔道的反应性当量计算方法。用ＷＩＭＳ（积分输运碰撞几率法）程序、ＧＡＭ、ＴＨＥＲＭＯＳ程序方法产生控制棒区的截面。用一维离散坐标法程序ＯＮＥＳＮ１产生泄漏谱，在控制棒区的Ｂ４Ｃ从１０群截面归并为４群截面；用二维离散坐标法程序ＳＮ２ＤＢ，在控制棒区（包括Ｂ４Ｃ、空隙和部分石墨反射层）考虑空间效应做均匀化４群截面。有控制棒和无控制棒情况下的Ｋｅｆｆ本征值，是由三维扩散有限差分程序ＣＩＴＡＴＩＯＮ在（θ，ｒ，ｚ）坐标系条件下做出的，并给出控制棒的反应性当量。     

低耗散TVD格式的构造

目的　为进一步提高双曲型方程数值解的精度,减少数值耗散,同时保证稳定性,处理有强间断的流动问题; 方法　从总变差减小（ＴＶＤ）定义出发,分析保证数值解ＴＶＤ的条件,找出满足研究目的的可能性和途径; 结果　首次提出和证明双曲型方程数值解ＴＶＤ的充分必要条件,并据此简化导出适用于时间显示格式的ＴＶＤ充分条件,并据此构造了一种低耗散ＴＶＤ 格式; 结论　进行的数值实验证明,高分辨格式可以进一步减小数值耗散,提高数值解精度;     

双曲型线性方程三阶和四阶TVD格式的新构造

利用Taylor级数理论和总变差减小(TVD)格式的充分条件构造了时间二阶、空间五点三阶和四阶新TVD格式．给出了新TVD格式与传统TVD格式及近期建立的二阶新TVD格式用于线性双曲型方程的计算结果，表明本文新格式特别是四阶TVD格式具有比二阶新TVD格式和传统TVD格式峰值衰减更慢、间断更陡，而计算工作量具有与传统二阶TVD格式相当的良好数值性能．     

新型高分辨率格式及其在CFD的应用

 利用小波奇异分析构造了一种新型的高分辨率高效计算格式,借助小波奇异分析所找到的流场中奇异点区域和光滑区域,合理选择高分辨率、高精度WENO格式和高精度中心差分格式进行流场的离散求解;利用该格式对超声速二维前台阶流动问题、跨声速RAE 2822典型翼型的二维外流场、跨声速VKI-LS59涡轮叶栅的二维流场以及NASA Rotor 37跨声速压气机转子和NASA Rotor 67跨声速风扇转子的三维流场进行数值模拟,计算结果表明,采用这种新型格式可以用于内外流复杂流场计算。数值计算还发现,在采用相同网格的情况下,本文格式不仅具有通常WENO格式高分辨率的特点,而且计算效率较WENO格式高出3～5倍。     

TVD格式在多组分爆炸波计算中的应用

该文将原始的高精度TVD格式推广到含多组分的爆炸波的数值计算，建立了多组分混合气体的状态方程，并推出了特征值和特征向量矩阵。对一维球对称爆炸场（爆裂球问题）的数值计算结果表明：采用文中推广的多组分TVD格式获得的数值结果比在Lagrangian坐标系下用von Nuemann-Richtmyer格式得到的结果更为合理，克服了由人工粘性格式计算引起的爆炸波远场超压偏低的不足，同时较好地解决了在Euler坐标系下多组分的变比热比γ的问题。     

双曲方程的一种二阶TVD差分格式构造方法

本文采用一阶迎风差分格式作Ｔａｙｌｏｒ展开,消去低阶项,给出求解守恒型双曲方程初值问题的一种二阶ＴＶＤ 差分格式的构造方法,并导出可能形式上用于非守恒型双曲方程的差分公式     

用动力学的流矢量分裂法计算溃坝波问题

采用动力学的流矢量分裂法模拟浅水长波方程. 利用气体动力学方程与浅水长波方程 之间的比拟关系, 并结合TVD差分格式给出浅水长波方程的数值计算方法. 计算了溃坝波问 题. 数值结果与精确解的比较表明, 该方法能很好地模拟浅水长波方程中的间断解问题.     

溃坝洪水演进的数值模拟

开发了溃坝洪水演进数值模型,对模型控制方程采用有限体积法进行空间离散,采用具有总变差减小特性的MacCormack预测-校正格式进行时间离散.采用非规则局部不连续网格划分计算区域,增强了模型对复杂地形的适应能力;具有总变差减小特性的数值格式使模型具备了强大的间断波捕捉能力,并消除了模型稳定性对人工黏性的依赖;基于单元属性与基于界面属性的动边界处理技术使模型能合理地模拟非恒定性强、干湿变动剧烈的水陆边界.实际算例表明：模型在强非恒定流状态下具有稳定性及守恒性,可作为溃坝洪水预报及其灾害防范的有力工具.     

全三维粘性流场计算的快速收敛方法

以新的高收敛率的 L U型隐式格式和高精度、高分辨率的 MU SCL TVD迎风格式为基础 ,提出了一种快速求解三维粘性流场的计算方法。为有效减少计算工作量 ,近固壁区采用了壁面函数方法。湍流模型采用简单的混合长度模型。为提高迭代收敛速度 ,从多重网格理论出发 ,结合改良型L U隐式和迎风格式 ,对 Reynolds平均三维可压缩 Navier-Stokes方程进行求解。用此方法数值模拟了 NASA L ewis37跨音速压气机转子流场。计算得到的流场与实验结果进行了对比。结果表明 ,此方法可快速得到三维粘性流场的流动特性且计算结果精度较高 ,可望在叶轮机械多叶片排全三维粘性流场校核计算中发挥作用。     

一维溃坝问题的间断Galerkin方法

将间断Galerkin方法应用于一维溃坝问题中,采用了两种不同数值流函数,对其进行比较。使用TVDM,TVBM型的限制器,同时利用一种改进的限制器来消除振荡。与一种高分辨率TVD性质的差分方法做了比较。结果表明间断Galerkin方法得到了更加逼真的图像。     

一种二阶TVD差分格式构造方法

采用一阶迎风差分格式，作Taylor展开．消去低阶项。给出了求解一维和二维等熵气动力学方程组的一种双参数二阶精度TVD差分格式的构造方法．