凝聚介质中斜激波的反射

结合数值模拟和理论分析,研究凝聚介质中斜激波反射.采用龙格-库塔控制体积间断有限元方法,数值求解"刚性气体"状态方程形式下的欧拉方程组;理论分析凝聚介质中斜激波反射模式;运用激波极曲线理论,给出典型激波强度下正规反射向非正规反射过渡的临界角及波后状态.比较数值模拟结果和激波极曲线理论分析结果,得到典型弱、强斜激波的反射图像.     

非定常欠膨胀射流的正格式数值模拟

 采用二阶正格式方法对非定常欠膨胀射流进行了数值模拟。将二维守恒方程的正格式方法推广到轴对称Euler方程组的求解,并对不同马赫数下的燃气射流进行了数值计算。计算结果表明,该方法能够较好地捕捉到包含膨胀波、入射激波、反射激波、马赫盘、射流边界以及三波点等复杂射流流场的波系结构,与实验照片反映的流动特征以及已有的数值结果相吻合。表明该方法对间断解具有较强的捕捉能力,在激波阵面上不会出现数值振荡。     

应用多介质PPM方法计算斜激波与物质交界面的相互作用

利用多介质PPM方法研究斜激波与物质交界面的相互作用.采用与体积分数耦合的Euler方程组作为计算模型,用双波近似来求解一般刚性气体状态方程Riemann问题.通过体积分数的计算来获得界面的位置,在整个流场采用统一的高阶PPM格式进行计算.文中对斜激波与不同物质界面相互作用进行了数值模拟,并给出了交界面上由于斜压效应产生的涡列的演化过程,特别是强斜激波与不同物质界面的相互作用的情况.     

对称TVD差分格式的构造及其应用

本文从ＧＣＩＲ差分格式出发，通过附加中心加权形式的反扩散流量，构造出一类具有二阶精度的对称型ＴＶＤ格式以及相应的一般形式的隐式格式，并针对双曲型方程组问题进行了讨论。同时文中还采用这类格式求解Ｅｕｌｅｒ方程，成功地数值计算了具有激波的流动问题。     

DG方法求解可压缩气固两相流动

给出求解双向耦合可压缩气固两相流的间断有限元方法,对所得到的气固两相流方程组不需要采用分裂的方法离散,对气相、颗粒相方程及其对流部分和源项可以统一处理,两相都采用基于近似Riemann解的数值通量.数值模拟低压含尘激波管内的两相非平衡流动,并与平衡流、冻结流的结果进行比较.分析颗粒相的存在对气体运动的影响,及激波后松弛区域内两相间相互作用规律.发现颗粒质量比决定两相平衡后的最终状态,而颗粒直径决定两相流从非平衡到平衡的过渡过程,即不同尺寸颗粒对应的驰豫时间、松弛距离不同.结果表明：本文提出的计算方法对求解可压缩气固两相流是可行的,为研究复杂的气固两相流动问题奠定了基础.     

刚体定点转动的欧拉定理

刚体定点转动是一个著名的力学问题．在Ｈ．Ｇｏｌｄｓｔｅｉｎ所著的《经典力学》一书中，有一段内容利用矩阵论和代数方程论论述了这个问题的欧拉定理［１］．本文在同一数学框架中讨论这个问题，而比文献［１］的论述显得简要；另外，还严格证明了描述刚体定点转动的转动（变换）矩阵的复特征值之幅角就是欧拉定理中的旋转角     

非定常空气流场中污染气体的对流与粘性扩散

基于离散涡方法求得的非定常、不稳定空气流场，数值求解了污染气体在空气中的对流湍流扩散方程．两者在每一个时间步耦合，以求非定常对流项．考虑速度脉动，改进了扩散模式．用SF6作为示踪气体，进行了现场测量实验．测得了各测点的SF6浓度随时间的变化．实验结果与数值仿真进行了对比，两者吻合．     

N2分子X^1Σg^＋,A^3Σu^＋和B^3Пg的势能函数

用文献［１，２，８］介绍的方法推导了Ｎ２分子的基态（Ｘ＾１Σｇ＾＋）的激发态（Ａ＾３Σｕ＾＋和Ｂ＾３Пｇ）的合理离解极限。计算并比较了在６－３１１Ｇ＾＊基集合，ＵＨＦ、ＣＩＤ、ＵＣＩＳＤ和ＱＣＩＳＤ水平下Ｎ２分子上述三个电子态的平衡结构和谐振频率；并用ＱＣＩＳＤ／６－３１１Ｇ＾＊计算了各态的系列单点势能值，由正规方程组拟合Ｍｕｒｒｅｌｌ－Ｓｏｒｂｉｅ函数得到了相应各态的完整势能函数，结果与实验数计     

数值模拟二维槽道内Oldroyd—B流体的流动

提出了一种求解Ｏｌｄｒｏｙｄ－Ｂ流体流动的数值方法。该方法采用有限体积离散，在求解连续方程时引进了Ｃｈｏｒｉｎ的人工可压缩性；求解本构方程时考虑了迎风处理的技巧；离散后的方程采用Ｇａｕｓｓ－Ｓｅｉｄｅｌ线松弛迭代求解全流场。     

基于双重介质模型的正仲氢转化器流动规律研究

为了研究正仲氢转化器中流体的流动规律，建立了催化剂堆积床的颗粒孔隙-堆床空隙的双重介质模型。基于双重介质模型分别建立颗粒孔隙和堆床空隙中正氢和仲氢的质量守恒方程，两种介质之间通过流量交换源项相互耦合。采用控制容法离散控制方程，并数值求解离散方程组。数值算例表明，由于颗粒孔隙与堆床空隙的差异，两种介质中的流动规律存在明显差异。     

N_2分子X~1Σ_g~ ,A~3Σ_u~ 和B~3Π_g的势能函数

用文献［１，２，８］介绍的方法推导了Ｎ２分子的基态（Ｘ１Σ＋ｇ）和激发态（Ａ３Σ＋ｕ和Ｂ３Πｇ）的合理离解极限。计算并比较了在６－３１１Ｇ基集合，ＵＨＦ、ＣＩＤ、ＵＣＩＳＤ和ＱＣＩＳＤ水平下Ｎ２分子上述三个电子态的平衡结构和谐振频率；并用ＱＣＩＳＤ／６－３１１Ｇ计算了各态的系列单点势能值，由正规方程组拟合Ｍｕｒｒｅｌ－Ｓｏｒｂｉｅ函数得到了相应各态的完整势能函数，结果与实验数计算值符合得比较好。     

超声速等离子体射流的数值模拟

基于可压缩的全Naiver-Stokes方程，利用PHOENICS程序对由会聚 辐射阳极形状等离子体炬产生的超声速等离子体射流进行了数值模拟.考虑了等离子体的黏性、可压缩性以及变物性对等离子体射流特性影响.研究了超声速等离子体射流的流场结构特性以及不同环境压力对等离子体射流产生激波结构的影响.结果表明，超声速等离子体射流在喷口附近形成的周期性激波结构是其和环境气体相互作用的结果. 关键词： 等离子体炬 超声速等离子体射流 PHOENICS     

振动叶栅中三维振荡粘性流场的压力密度消去法

1前言三维非定常流场的求解是目前国内外的一个热点研究课题山。文献［2］完成了三维可压非定常欧拉流场的求解,这一方法在求解三维非定常欧拉流场时,应用了四阶fringe－Kutta方法对控制方程进行积分,用中心差分进行空间离散,采用了四阶人工粘性项来保证计算格式稳定,计算稳定性要求严格,时间步长不能大,计算时间长。本文从非定常三维粘性N－S方程组出发,通过合理的数学方法,消去压力及密度项,得到只包含振荡速度矢量项对空间的偏微分方程组,在已知定常速度场后,这一方程组很容易求解。2基本方程在以角速度为n作旋转的相对坐…     

薄膜干涉条纹定域的确定

薄膜干涉条纹定域的确定，文献［１］～［４］等已有讨论．其中［４］较简单直观。本文在直角坐标系中推出定域中心曲面方程，作适当近似与文献［１］、［４］的结果相同，并在此基础上对条纹定域深度进行了讨论     

考虑转动能的一维/二维Boltzmann-Rykov模型方程数值算法

研究考虑转动能的Boltzmann-Rykov模型方程,基于转动自由度对气体分子速度分布函数矩积分,引入约化速度分布函数,应用离散速度坐标法与数值积分技术,将气体运动论模型方程化为在离散速度坐标点处关于三个约化速度分布函数的联立方程组.应用拓展计算流体力学有限差分方法,数值计算考虑转动自由度的双原子气体一维、二维Boltzmann模型方程,得到高、低Knudsen数一维激波管内流动和二维竖直平板绕流问题的流场,分析验证考虑转动能的Boltzmann-Rykov模型方程全流域统一算法求解一维/二维气体流动问题的可靠性.结果表明,气体稀薄程度与分子内自由度对流场具有较大影响,且Knudsen数较高的稀薄气体流动呈现严重的非平衡流动特点.     

等离子体电弧炉内流动与传热数值模拟

本文采用磁流体动力学(MHD)模型对直流等离子体自由燃烧电弧和电弧炉内部的流动与传热进行了数值模拟研究。通过对基于磁矢量势描述的电磁场方程组和流体力学方程组的耦合迭代计算,求解得到了流体的温度场和速度场等,计算结果清晰地反映出等离子体电弧的高温阴极射流现象,并与同行的实验和数值结果进行了对比。本模拟方法和结果对于电弧炉的工业应用和优化设计有重要的指导意义。     

自相似Euler方程的数值方法

Euler方程某些问题的解具有自相似特点，可以使用更准确的方法求解.提出了两种数值方法，分别称为自相似和准自相似方法，新方法可以使用现有守恒律方程的数值格式，无须设计特殊方法.对一维激波管问题、二维Riemann问题、激波反射以及激波折射问题进行了数值计算.对自相似Euler方程，一维计算结果显示数值解基本等同于精确解，二维结果也比现有文献计算的结果有更高的分辨率.对准自相似Euler方程，新方法可以求解不具有自相似性但接近自相似的问题，并在计算时间足够长时可以取得自相似Euler方程的效果.数值求解自相似Euler方程对自相似问题的研究，高分辨率、高精度格式的设计乃至Euler方程的精确解都有重要启示.      

关于铱(Ir)的晶胞常数a0值的更正

关于铱（Ｉｒ）的晶胞常数ａ０值的更正兰田徐飞岳罗敬良夏都灵Ｉｒ的ａ０值与文献值稍有差别（如３．８５，３．８４和３．８３９等），文献［１］表４和文献［２］表１误用为３．５８，虽然各ｄｎ值未错，但其驰豫量（％）也应修正于后。园括号中为错误数值。Ｉ...     

压缩因子Z〉1的定容高压气体的压力与温度的关系

改良了文献［１］的近似处理方法，从求解新的非线性常微分方程中，导出了Ｚ＞１的定容高压气体在（Ｔ０，ｐ０）附近的压力与温度的关系。对照Ｒｅｄｌｉｃｈ－Ｋｗｏｎｇ方程的计算值表明，在｜ΔＴ｜＜１５Ｋ时，其压力的平均相对误差小于０．９％。     

非平衡化学反应流场NS方程计算研究

介绍一种新的与化学问题相耦合的NS方程数值方法,用来计算有限速率化学反应再入流场.采用迎风TVD差分格式的激波捕捉法,求解控制方程。本项研究针对轴对称流动。然而,可以直接推广到三维流动和其它的化学反应系统。所采用的整个泰氟隆-空气化学反应模型是由23个化学组元和35个化学反应所组成。对钝头球锥体全流场(包括前体和近尾迹)进行了数值模拟,获得了充分收敛的数值解。对清洁空气流场计算的峰值电子数密度与实验结果符合很好。