微波瓣燃烧器中流向涡混合数值研究

为了检验波瓣喷管下游流向涡能否促进燃烧, 对尺度为毫米量级的微波瓣喷管混合燃烧器中流向涡的行为进行了数值模拟研究.研究发现, 微波瓣喷管燃烧器中远没有微波瓣喷管混合器中非常强烈的流向涡形成和发展, 造成差别的根本原因在于本该形成流向涡的主次气流接触界面变成了燃烧火焰的锋面, 燃烧反应气流体积的急剧膨胀使流向涡环量迅速耗散衰减.另外, 流场计算结果显示在4波瓣和8波瓣微喷管燃烧器中主次流接触初期有较明显的流向涡生成、发展和很快消失, 而在12波瓣微喷管燃烧器中基本没有明显的流向涡形成和发展.但是, 对比模拟表明, 在同样边界条件下12波瓣微喷管混合器中却有非常明显的流向涡的形成和发展.      

多压电膜式零质量射流激励器的设计与性能研究

设计了一种多压电膜式零质量射流激励器,并利用热线风速仪,得到了对应于射流中心线上固定位置处速度最大时的最佳频率.在该频率下,测量了零质量射流在不同孔径时沿中心线的速度分布;测量、分析和比较了不同激励膜数目对射流出口速度分布的影响.表明:在其他条件完全相同的情况下,一片振动膜时的激励器出口中心线上的平均速度和峰值速度,都较5片膜的情况小一半左右,而3片和5片压电膜时的激励器性能差别却不大.也就是说,虽然压电膜的数目增加能使激励器出口射流的能量增加,但这两者之间并非呈线性关系.此外还用PIV测量方法,直观地显示了射流出口的涡结构.      

高超声速钝楔边界层转捩大涡模拟

 以五阶迎风和八阶对称格式混合差分格式求解三维可压缩滤波Navier-Stokes方程,对来流Mach数为6.0、半锥角5°的高超声速空间发展钝楔边界层转捩至完全湍流进行了大涡模拟。时间推进采用紧致存储三阶Runge-Kutta方法,亚格子尺度模型为Driest因子修正的Smagorinsky涡黏性模型。通过定常流场入口边界附近吹／吸引入不稳定扰动斜波的方法数值模拟得到了层流失稳转捩直至完全湍流的空间发展全过程。对扰动的线性、非线性增长以及湍流斑的形成和发展进行了分析,给出了转捩及完全湍流下的速度相关量统计并与实验、DNS结果进行了对比分析,计算结果与理论及实验吻合。     

头部鼓包对不同截面机身侧力影响研究

 飞机或导弹在大迎角下飞行时会形成非对称涡,从而产生很大的侧滑力和偏航力矩。为了寻求消除侧滑力和偏航力矩或利用流动的非对称来提供飞机飞行所需的航向控制能力,研究了机头上放置小鼓包对3种不同截面的机身前体模型的影响。研究过程中采用了风洞测力和水洞流动显示的实验方法。结果表明,在20°~70°的迎角范围内,鼓包对圆锥的影响最大,对椭圆截面模型有一定的影响,对带棱截面模型的影响最小。鼓包越靠近前体头部,对侧力的影响越大。对圆锥柱体模型鼓包只能改变侧力最大值出现的迎角,而对减小最大侧力值作用不明显。对椭圆截面模型,可使最大侧力系数从1.98降到0.53,很大程度上降低了侧力。带棱截面模型对鼓包的大小和位置均不敏感。     

Bush有限单元原理及其在航天器结构建模中的应用

利用大量的数值计算确定了MSC.NASTRAN软件中广义弹簧单元（Bush）刚度阵的解析表达式,通过对比Bush元与剪切梁元的刚度阵,推导了Bush元参数与剪切梁刚度参数的对应关系,并利用展开状态的太阳翼结构和大型桁架结构的静力和动力分析验证了推导结果。该研究揭示了Bush有限单元的原理和所有输入参数的物理含义,有助于Bush元在结构有限元分析中的准确使用。     

雾化过程的一种Euler-Lagrangian耦合算法（英文）

对雾化过程的直接数值模拟需要巨大的计算资源和时间,而工程中的简化模型则经常会给出错误的结果。因此,可以折中地采用一种混合方法,即在不同尺度上采用不同的模型。提出了一种雾化过程的欧拉——拉格朗日耦合算法。较大的液团采用VOF法直接求解,与网格尺度相当或更小的液滴则采用双向耦合的拉格朗日粒子法进行追踪。而该方法要求粒子的体积小于网格体积的10%,为此又提出了一种虚网格粒子追踪法。由于湍流结构对雾化过程的影响很大,故湍流采用了大涡模拟模型。采用多个算例对开发的算法进行了验证,并对部分关键参数的影响进行了深入研究。采用新算法对两股撞击射流的雾化过程进行了研究,瞬态和统计结果均表明新算法能够给出良好的预测。     

基于MSFLE与Liutex的涡合并实验研究

应用具有拉格朗日性质的运动单帧长曝光(MSFLE)图像测量方法,对矩形管内湍流边界层涡结构合并现象开展了实验研究。研究了雷诺数Reθ在97到194之间管道湍流边界层流向-法向平面内涡合并现象。在实验中,测量系统以与涡运动速度相近的速度匀速移动,采用连续长曝光记录示踪粒子运动轨迹,以捕捉涡随时空的演变过程,骨架提取图像处理算法获得示踪粒子速度,并应用Liutex物理量表征涡的旋转强度。研究表明:MSFLE方法测量装置简单,对实验条件要求低,可从拉格朗日视角直观测量湍流边界层涡结构及周围流场的时空演变过程,MSFLE图像测量方法与Liutex涡识别算法结合可以很好地应用于湍流边界层涡结构的可视化与量化。矩形管内湍流边界层涡合并的条件是两个涡相邻、强度和尺寸基本相同且为同向旋转涡,合并中两个涡的强度呈反向变化,合并生成的新涡的强度和尺寸基本为初始合并时两个涡之和,且旋转方向与两个涡同向。     

桥梁若干流致振动与卡门涡街

首先介绍了与桥梁流致振动相关的卡门涡街、涡致振动和颤振的概念。在此基础上,讲述历史上塔科马大桥风毁事件的研究争论,并再次澄清塔科马大桥风毁的直接原因是颤振而不是涡激共振。对虎门大桥振动问题,针对风工程领域中涡致振动问题的性质进行了讨论,指出限幅涡振问题,并伴随的锁频现象,本质上仍是结构自激振动问题,卡门涡街起到的是诱发而不是直接的推动作用。最后,对形形色色流致振动问题的差异和共性进行了讨论。     

同轴旋转倒置圆台环隙间流动特性分析

针对同轴旋转倒置圆台环隙间流体复杂流动问题, 对其环隙间流动特性进行了实验研究。重点进行染色液流动显示实验和PIV流场测速实验, 对实验结果做定性及定量分析, 研究内筒转速和环隙宽度对环隙间流动特性的影响。染色液流动显示实验和PIV流场测速实验分别定性和定量地展现了环隙间螺旋涡的产生及变化过程。对不同内筒转速和环隙宽度下的螺旋涡涡心运动周期进行分析, 结果表明, 内筒转速升高, 周期减小;环隙宽度增大, 周期增大。运用瞬时流动和时均流场解析了环隙间螺旋涡运动产生机制, 探究内筒转速和环隙宽度对3种雷诺应力大小的影响与分布情况。内筒转速变化, 雷诺径向正应力始终最大;环隙宽度变化, 雷诺切应力始终最小。      

基于释热率分析的钝体贫油熄火过程分析

针对锥形钝体稳定的甲烷-空气预混湍流火焰复杂的熄火过程,采用大涡模拟(LES)与输运方程概率密度函数(TPDF)湍流燃烧模型相结合的模拟方法,研究远离熄火、近熄火及熄火点的火焰和释热率(HRR)数值变化情况,定量分析熄火判据。结果表明:冷态速度模拟结果和实验的相对均方根误差在10%以内,热态误差在20%以内;释热率是否出现在OH和CH₂O重叠的区域,是判断熄火的一个重要参数;远离熄火时,释热率高的区域主要沿内侧剪切层出现;近熄火工况下,释热率在流向轴上闭合,回流区下游也出现较大的释热率;熄火点工况下,较大释热率的区域在回流区下游和上游均有出现;模拟预测的熄火情况和实验PLIF结果一致;平均释热率可作为判断熄火的定量依据,即当钝体后方0.2d处内侧剪切层平均释热率与回流区平均释热率的比值小于4时,发生熄火。