旋转爆轰流场的数值模拟

基于带化学反应的二维Euler方程,采用氢气-空气的9组分19步基元反应简化模型,对充有当量比的氢气-空气预混气和空气的环形旋转爆轰流场,从点火燃烧到发展成旋转爆轰的过程进行了数值模拟。根据数值结果分析了波后流场中爆轰产物受激波、高温和离心力等作用而挤向外壁,形成有利于充入燃料,实现持续稳定旋转爆轰的流场特征。还讨论流场中爆轰波、激波与间断面和内外壁面反射或折射,从而形成多个激波相交的波系特征。为认识和理解旋转爆轰流场,开展旋转爆轰的实验研究等具有现实的指导意义。     

气流出口型面对凹面腔内激波聚焦起爆的影响

为提高两级脉冲爆震发动机工作的可靠性,探寻强化凹面腔内激波聚焦起爆爆震波的方法,通过数值模拟,采用氢气作为燃料、空气作为氧化剂,探讨了3种凹面腔气流出口型面及出口面积对激波聚焦起爆的影响。结果表明:采用垂直出口壁面有助于提高凹面腔内的激波聚焦起爆压力和温度,使起爆时刻提前,随着出口面积的减小,激波聚焦峰值压力有较大的提高;但效果不如垂直壁面,随着倾斜壁面出口面积减小,起爆点的压力和温度只有小幅度增加,起爆时刻基本无变化;采用弧形出口壁面,不能对向凹面腔底部传播的激波产生反射强化作用,不能有效起爆爆震波。对比3种出口壁面型面,垂直壁面能更有利于凹面腔内的能量聚集,强化激波聚焦起爆爆震波。      

异形腔体压电声衬声学性能

为提高压电声衬对低频噪声的抑制范围,对声衬腔体进行结构优化。利用平面波理论构建了两种曲线管道的声学物理模型,并分别建立了两种模型的传递矩阵,以此作为异形腔体亥姆霍兹共振器传递损失计算的理论依据,并通过仿真验证其正确性。结合压电振子的形变对声衬进行有限元仿真分析,结果表明:在压电振子施加500V驱动电压时,两种声衬频率偏移量分别为115Hz和120Hz。与圆柱形腔体声衬进行对比结果表明:在相同腔体厚度范围内,由曲率越大的曲线所生成的腔体,在相同驱动电压条件下,频率变化率越高,这为今后对声衬腔体结构优化提供一种有效的依据。      

前后缘型线同时可控的乘波体设计

提出了一种前后缘型线同时可控的乘波体设计方法,在马赫数可控的外锥形曲面基准流场中,结合流线追踪技术和混合函数,实现了椭圆前缘转椭圆后缘的乘波体设计,并在设计点(Ma=6.0)和接力点(Ma=4.0)对其进行数值仿真研究。前后缘同时可控的乘波体在型面剧烈过渡处产生了较弱的激波,出口两侧存在高温高压区,后部对称面附近的激波形状由圆弧变为平直线且出口处流场基本均匀,非常有利于与进气道匹配设计。另外,该乘波体具有较高的容积率和预压缩效率,附面层修正后的容积率为0.24,设计点时乘波特性较好,接力点时前部完全乘波,具有较高的升阻比,有黏条件下设计点和接力点的升阻比分别为2.54和2.41。此外,与给定前缘的乘波体相比,其升力、阻力、俯仰力矩和出口增压比都有明显增加,但是升阻比和出口总压恢复系数有所降低,在设计点无黏升阻比由3.56降为3.00。以上研究表明,本文的设计方法可行且更加灵活,拓宽了乘波体的选择范围。     

一种部分森林覆盖环境电波传播特性预测方法

针对部分森林覆盖环境中传播媒介的非均匀性,提出一种基于动态网格的抛物方程(PE)预测方法,该方法采用宽角PE法进行大区域电波传播特性预测。为平衡PE求解的准确性与时效性,引入动态网格技术。通过仿真算例,探讨了基于动态网格的PE法在部分森林覆盖环境电波传播特性预测中的有效性。结果表明:在计算精度相似的情况下,相对于均匀细网格,动态网格会大大减少PE所需的计算时间;而相比均匀粗网格,动态网格具有更高的计算精度。     

高速气流冲击柱状泡沫多孔体的传热特性分析

采用流体/多孔区域一体化单区域算法,数值研究了高速绕流条件下前置于圆柱体前缘表面的柱状泡沫多孔体内部的传热特性。基于蒙特卡罗法考虑多孔域内的辐射热效应,分析了变化多孔区域长度和多孔阻力特性对模型激波阻力和前缘多孔区域气动热的影响。结果表明:在圆柱体前缘安置一定长度及带有适当阻力特性的泡沫多孔材料,可同时减小整体激波阻力并降低前缘表面的气动热效应。在模拟工况下,无量纲长度1.0、黏性阻力系数0.2×107m-2及惯性阻力系数200m-1的前缘泡沫多孔可减小激波阻力13.5%,降低约75%的前缘表面的平均气动热流密度。保持无量纲长度不变,减小泡沫多孔区域惯性阻力系数会降低激波阻力,但会略微增加前缘壁面气动热流密度。      

流向强迫作用下的液体初始雾化机制及动力学特征

使用Volume of fluid (VOF)方法和基于树形数据结构的自适应算法来研究射流雾化的破碎过程以及扰动对射流破碎机理产生的影响。在无扰动情况下,液体射流的头部、液丝和液滴随着射流时间的发展不断发生演变。射流头部先呈现蘑菇状外形,随后液丝生成,并慢慢转变成网兜状,直至断裂形成小液滴。在周期性流向强迫的作用下,射流液柱的表面会形成周期波,其液丝破裂形成液滴的时机与稳定射流情形相比会有所提前,射流形成的头部更趋于扁平,最终生成的液滴数量更多。低中频阶段,随着扰动频率的增大,射流未扰动液柱长度(L)逐渐缩短,液滴直径的概率密度分布(PDF)趋于尖锐,液滴平均直径(SMD)增大。在高频阶段,随着扰动频率的增大,L会随之增大,液滴直径的PDF分布变得平缓,SMD会减小。     

基于模拟退火粒子群算法的波浪发电系统最大功率跟踪控制

波浪发电系统最大功率点跟踪控制中,传统粒子群算法存在早熟收敛和局部搜索能力不足问题,为此提出基于模拟退火算法的粒子群优化方案。该算法每次更新粒子的速度和位置时,通过比较当前温度下各个粒子的适配值与随机数的大小,从所有粒子中确定全局最优解的替代值,从而使粒子群算法在发生早熟收敛时能够跳出局部最优并快速找到全局最优解。仿真结果表明,与传统粒子群优化算法相比,模拟退火粒子群算法可有效避免波浪发电系统陷入局部最大功率点,并快速实现全局最大功率跟踪,提高了波浪能捕获率。     

Numerical study of compression corner flowfield using Gao-Yong turbulence model

A numerical simulation of shock wave turbulent boundary layer interaction induced by a 24° compression corner based on Gao-Yong compressible turbulence model was presented.The convection terms and the diffusion terms were calculated using the second-order AUSM (advection upstream splitting method) scheme and the second-order central difference scheme,respectively.The Runge-Kutta time marching method was employed to solve the governing equations for steady state solutions.Significant flow separation-region which indicates highly non-isotropic turbulence structure has been found in the present work due to intensity interaction under the 24° compression corner.Comparisons between the calculated results and experimental data have been carried out,including surface pressure distribution,boundary-layer static pressure profiles and mean velocity profiles.The numerical results agree well with the experimental values,which indicate Gao-Yong compressible turbulence model is suitable for the prediction of shock wave turbulent boundary layer interaction in two-dimensional compression corner flows.      

流体控制矢量喷管启动过程的瞬态响应

采用三维非定常数值模拟方法对流体矢量喷管启动的瞬态响应过程进行了研究,给出了流场的演化特性及性能参数的响应时间.根据不同时刻的流场和壁面压力分布,分析后得到了注气附近射流扰动波系及漩涡的演化过程,并与上壁面压力分布相互印证.由性能参数(包括流量,推力系数,矢量角)的时间响应曲线发现:流量和推力系数在响应初始阶段存在剧烈波动,而矢量角是平稳变化,三者的响应时间在毫秒量级.