激波与堆积粉尘相互作用

实验中利用带示踪粒子的X光脉冲摄影和阴影摄影技术,成功拍摄了激波与堆积粉尘作用后粉尘内部的波系图像和波后卷扬粉尘云的形状,进而测得入射激波倾斜角、透射激波倾斜角和波后物质界面折射角。在此基础上,利用相关理论对该流场进行分析,确定出该堆积粉尘状态方程中中的待定系数。计算结果与测试结果符合较好。     

激波与火焰相互作用过程的实验研究与数值模拟

在长4.8m、截面为60mm×60mm的卧式激波管中,进行激波与火焰相互作用过程的实验研究。利用YA16高速摄影仪拍摄了该过程的照片,显示了激波作用下火焰的失稳、变形和破碎。基于带化学反应的NavierStokes方程和有关热力学和反应动力学数据,利用改进的VLS格式,对该过程进行了数值模拟。根据实验和数值计算结果,讨论了火焰稳定性和漩涡生成机理,且理论计算与实验结果较为一致。     

可压Euler方程的格子Boltzmann方法:二维情形

根据微观和宏观之间的质量、动量和能量守恒原理，构建了一个可用于模拟可压流体的D2V19格子Boltzmann模型，利用Chapman—Enskog分析可从该模型推导出宏观流体力学方程-Euler方程．利用该模型对三种典型激波管进行数值模拟并与Riemann解析解进行比较和分析，发现该模型能较好的捕捉激波，精度令人满意，证明了该模型对可压流体的适用性．该模型适用于气体多方指数为任意值的系统．     

甲烷-空气预混气体泄爆过程的数值模拟与实验验证

通过求解带化学反应的Euler方程，对甲烷-空气预混气体泄爆过程进行了数值模拟。采用分裂格式处理方程，其中采用二阶TVD格式求解流场，采用基于Gear算法的微分方程解法器求解化学反应过程。计算结果显示在外流场中出现了破膜激波和第二道激波。破膜激波逐渐衰退，第二道激波经历了由增强到衰退的过程。计算结果与实验结果在定性上是一致的。     

改进的强隐式格式及其在三维Euler与N-S方程中的应用

将 Ｓｔｏｎｅ强隐式格式由单个方程推广到方程组，并用于三维高速可压缩流的 Ｅｕｌｅｒ和Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组的数值计算。为较好地捕捉激波、提高激波的分辨率，对离散后方程右端项的数值通量进行了改进，采用了ＮＮＤ格式的数值通量。典型三维算例表明，本格式具有高效率、高分辨率的特点，流场与实验数据较接近。     

二阶非线性奇摄动方程的激波解

考虑了一个二阶非线性奇摄动方程的Dirichlet问题.利用微分不等式理论证明了该问题解的存在性,并指出了解的上下界,还利用匹配法求得了该问题的激波解,并给出了该问题的激波解、激波位置与边值的关系.     

双波长NO-PLIF法测量运动激波前后温度场的实验研究

利用双波长NO-PLIF法,研究激波管内运动激波前后温度场分布.实验中以NO作为温度测量的示踪剂,分别选取波长为225.072 nm和225.282 nm的染料激光作为激发光,获得了较好的荧光图像.对入射激波流场的1区和2区温度分布进行了理论估算和实验测量,结果表明NO的荧光强度对温度的敏感程度与激发波长有关,当波长为225.072 nm时,荧光强度随温度变化不明显,当波长为225.282 nm时,荧光强度随温度变化而剧烈变化.实验测得的温度分布相对理论估算值的误差约为10％左右,这表明该方法具有较高的精度,具有一定的工程应用价值.     

一类二阶非线性奇摄动方程的激波解

利用边界层理论和匹配原理，讨论了一类二阶非线性奇摄动方程的激波问题。分别对区域的两个端点和内部进行了分析，得出了与边界条件相对应的激波位置解的不同表达式。     

变截面激波管中激波运动规律的初步数值研究

利用有限差分方法求解变截面激波管中存在激波的二维流体流动问题。着重研究初始流场为静止和超音速两种情形。分析讨论所得计算结果,并和已有实验结果进行了比较。     

一种新的自由曲线的造型方法

研究了型值点与曲线的位置，利用曲线方程的已有属性，构造出一种样条曲线，使曲线通过型值点又不必作反求计算。并根据曲线方程的特点，对参数的选取进行了分析。通过采用适当的参数化方法，曲线的光顺性得到了进一步提高。理论分析和编程实验证明以上方案是可行的。该算法的运算速度快，生成的曲线有较好的光顺性，有一定的应用价值。     

激光支持等离子体爆轰波流场研究

采用高速纹影摄影法记录激光支持等离子体爆轰波流场的演化,得到流场冲击波阵面形状和尺寸随时间变化规律。建立激光支持等离子体爆轰波流场演化的物理模型,采用计算流体力学方法模拟计算激光击穿空气后,激光作用时间内及激光作用结束后等离子体流场演化过程。结果表明,纹影摄影法能有效观测到激光支持等离子体爆轰波的传播情况。等离子体爆轰波在激光作用初期为椭球形,随着时间的增加逐渐转变为球形。计算的等离子体爆轰波流场密度分布图与实验结果基本吻合。     

三维强旋流场中的激波捕捉

为了研究激波对气体离心机取料支臂附近流场的影响,采用数值求解三维N av ier-S tokes方程组的方法来模拟三维强旋流场。使用有限体积法中的二阶迎风格式离散方程,隐式迭代求解。计算捕捉到了障碍物前方的脱体激波,随着径向位置的增大,激波的强度逐渐变大,激波面与支臂表面的距离越小;支臂与转子侧壁之间形成了斜激波,向下游延伸,同时激波面向内弯曲;激波对转子侧壁有影响;激波对上边界的影响很小。计算结果表明该文的方法对强旋流场的计算是有效的。     

Stability in channel flow with fiber suspensions

The constitutive equation of fiber suspensions is established on the basis of fiber orientation tensors. The modified Orr-Sommerfeld stability equation is obtained further and numerically solved by aid of spectral method and finite difference method. The computational results of channel flow without fibers agree well with the experimental data with a higher degree of accuracy than previous numerical results. The results of the channel flow with fiber suspensions indicate that the presence of fibers attenuates the instability of flow, increases the critical Reynolds number, reduces the growth rate of perturbations and narrows the range of unstable waves. The extent of the effect of fibers on the flow stability is in direct proportion to the volume fraction and aspect-ratio of the fibers.     

炸药颗粒床中致密波状态方程研究

对一维两相流模型进行简化,引入颗粒守恒假设,建立了炸药颗粒床的定常致密波计算模型.分别使用Tait状态方程和Hayes状态方程,研究了致密波区内的压力、固相体积分数、固相密度、颗粒半径和颗粒数密度的分布.通过比较发现,使用Tait状态方程,固相密度、颗粒半径和颗粒数密度的分布以及它们之间的对应关系与活塞驱动颗粒床产生稳态致密波的过程不相符.结果还表明,大于颗粒材料音速的致密波(超音速致密波),波阵面存在间断;而小于颗粒材料音速的致密波(亚音速致密波),从波前到波后,所有物理参量光滑过度.分析得到了活塞速度、初始体积分数与各物理参量之间的关系.     

颗粒流本构关系的实验研究

基于"拟流体"的思想,给出了测量颗粒流本构关系的实验方法.在颗粒斜槽流实验中,通过颗粒的抛物线运动计算流层速度分布,根据非牛顿流体理论求得颗粒流黏性的本构关系;建立了颗粒斜槽流的数学模型,流层的速度为指数分布,流量为斜槽倾角和流层厚度的函数.以小麦颗粒为例,实验结果与用"拟流体"方法所得的预测值进行比较,相对误差在13%以内.     

水平管入口段内水煤浆流动特性数值模拟

以Eulerian多相流方法为基础建立了水煤浆流动的液固两相流体动力学模型,并对水平管入口段内水煤浆流动进行了模拟.模型采用颗粒动理学理论求解固相本构方程,采用RNGk-ε湍流混合模型描述颗粒间具有强烈作用的两相湍流流动.针对水煤浆中煤粉颗粒的双峰分布特性,将煤粉看作2种大小不同的固相,同时考虑固相与液相、固相与固相之间的动量交换.模型有效性通过Kaushal等试验结果和水煤浆的压降试验验证.通过模拟考察了入口段内速度及浓度分布过程和入口段长度的变化规律.结果表明:重力和颗粒间的相互作用对浓度和速度分布过程具有重要影响;固相体积分数在30%~49.5%范围内,入口段长度随浓度的增加而减小,随平均流速(0.2~5.0m/s)先增加后减小.     

广义二阶流体的脉冲泊肃叶流动分析

讨论了广义二阶流体的脉冲泊肃叶流动,引入黎曼-刘维尔分数阶微分建立本构方程,结合不可压缩流体时间分数阶动量方程,得到控制方程。利用傅里叶正弦变换和分数阶拉普拉斯变换,获得流体速度解析解。利用Stehfest算法对结果进行数值模拟,通过图像讨论了分数阶参数以及延迟时间对流动的影响。结果表明速度过冲现象主要取决于动量方程的时间分数阶参数。     

基于改进桥联模型的CFRP层合结构抗冲击载荷动态力学性能研究

为研究碳纤维增强树脂基复合材料（CFRP）在冲击载荷下的力学响应,本文基于桥联模型,采用含应变率效应的动态本构方程以及屈服准则,对CFRP层合板进行了优化设计和动态压缩强度分析.利用相关文献中的层合板准静态及动态试验结果,对改进的桥联模型程序进行验证,结果吻合良好.在此基础上,对比分析了几种CFRP层合板铺层设计方案,优选并制备其中一铺层的层合板,实施了中等应变率的压缩试验,前期理论计算结果与试验结果有较好一致性.本文工作可用于抗高冲击载荷CFRP层合结构的优化设计及性能预测.