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

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

应用傅里叶时间谱方法求解二维跨音速流动问题解的精度研究

为了从时间离散精度的角度评估傅里叶时间谱方法的实际计算效率,对用该方法求解二维跨音速流动问题的解进行了精度研究。通过求解Euler方程模拟俯仰振动NACA0012翼型的周期性非定常流动。结果表明:傅里叶时间谱方法可以使用很大的时间间隔较为准确地模拟有激波存在的周期性跨音速流动;对于此类流动问题,相比二阶向后差分公式法,傅里叶时间谱方法在预测翼型表面压力系数积分量特别是升力系数的时间变化规律方面具有明显的计算效率优势;用于模拟翼型表面压力系数自身的时间变化规律时,傅里叶时间谱方法具有不低于二阶向后差分公式法的计算效率。     

球形点火腔对脉冲爆震DDT过程影响的数值研究

为了减小脉冲爆震发动机PDE（pulse detonation engine）单次循环时间,缩短爆燃向爆震转变DDT（deflagration to detonation transition）距离是关键.为此,提出了一种在头部添加球形点火腔的新型爆震管道结构设计方案,并采用丙烷和空气为可爆混合物,通过对新型爆震管DDT过程的二维数值模拟,研究其对DDT距离以及DDT时间的影响.数值模拟结果表明,头部采用球形点火腔后,点火腔中的压缩波经过多次反射后,能够在爆震管中更快地促使爆震波形成,DDT距离大大缩短;当球形点火腔直径为1.5倍等直爆震管直径时,相对于常规等直管爆震管结构,其DDT距离和过程时间分别减少了14%和16.26%.      

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

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

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

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

软式空中加油受油机头波数值仿真分析

针对软式空中加油对接过程中受油机头波加剧软管锥套飘摆、降低对接成功率的问题,建立了受油机三维模型,利用Delaunay方法自下向上生成非结构网格.通过计算ONERA M6三维翼型表面压力分布和压力系数,验证了三维Navier-Stokes方程和SST k-ω湍流模型的有效性.在此基础上对不同飞行条件下受油机头波进行静态数值计算,分析了头波随马赫数和迎角的变化规律,并通过仿真解释了头波影响下软管锥套的飘摆原因,预测了对接时锥套的运动趋势.     

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

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

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

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

侧板边界层对正激波/边界层干扰特性影响的数值模拟研究

为了研究外压式进气道处于临界工况时结尾正激波同时与压缩面和侧板上边界层的相互干扰,专门设计了三组具有不同前伸长度侧板的简化构型作为研究对象,利用数值模拟手段评估了侧板边界层厚度对正激波/边界层干扰特性的影响。仿真结果表明,当没有侧板边界层参与干扰时流动呈现较好的准二维特性,但当侧板边界层参与干扰后将形成较强的角区干扰结构,该角区三维干扰结构与对称面上的准二维干扰结构存在此消彼长的关系。此外,波系的空间形态也将由“准二维λ波”结构变为“双λ波+角区压缩波”结构,波系形态的改变则进一步导致壁面回流区分布以及摩擦力线拓扑结构变得更加复杂。     

捕获加速新生离子对Alfvén波的激发

根据准线性方程,求出了在离子持续产生并被太阳风捕获加速的情况下离子的渐近态分布函数.采用此分布函数分析了斜向传播的Alfven波的稳定性.