基于欧拉方法的2 维翼型冰晶结冰数值计算

针对现有航空发动机冰晶结冰试验方法存在数据采集困难,且难以揭示冰晶结冰机理问题,对某楔形翼型进行冰晶结 冰数值计算。采用S-A 湍流模型获得流场结果,应用欧拉方法获得冰晶和液滴轨迹结果,基于Messinger 模型获得冰形,并以 NASA-NRC 的第139 号试验结果验证了数值方法的可行性。结果表明：在其他工况相同时,压力越低（飞行高度越高）,冰晶/ 液滴蒸 发越强烈,结冰情况越严重;分析不同液态水含量与总水含量比例（CLW/CTW）对结冰的影响可知,液态水过少或过多都不利于冰层形成。     

壁面催化效应对高超声速气动热影响研究

高超声速飞行器壁面催化效应会导致激波层中原子在壁面处复合释热,加剧周围气动热环境。针对高超声速流动壁面催化特性,选择不同飞行马赫数及高度条件,采用完全催化和非催化两种条件对球锥模型壁面热流率进行数值模拟计算,研究壁面催化效应对气动热的影响规律。结果表明,固定飞行高度时,壁面催化效应对气动热的影响随马赫数增加而加强,Ma=25条件下驻点处完全催化与非催化热流比值高达1.92;固定飞行马赫数时,在50km高空以上壁面催化效应对气动热的影响随高度增加而减弱;壁面催化效应不仅会影响壁面附近的流场特性及组分分布状态,而且对整个激波层都有一定的影响作用。     

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

空心气冷低压涡轮动叶气热耦合数值模拟

以空心气冷低压涡轮动叶为研究对象,采用高质量的流体域和固体域网格控制技术,带转捩模型的双方程SST湍流模型,开展了基于CFD方法的叶片气热耦合问题研究。获得了不同冷气流量比(分别为1.0%,1.38%,1.8%和2.2%)、温比(分别为2.1,2.25,2.3,2.4和2.5)和压比(分别为1.4,1.6和1.8)对叶片换热特性的影响规律,设计状态中截面按弧长平均的叶片壁面金属温度计算值较试验值偏小0.3%,气热耦合计算的叶片壁面温度分布与试验结果吻合良好,验证了气热强耦合计算方法的精度,为涡轮叶片温度场分析提供了一种有效的方法。     

基于并行CFD和优化技术的返回舱外形多目标优化设计北大核心CSCD

在给定的质心设计范围内,围绕球冠倒锥返回舱外形的高超声速气动单点静稳定性、配平升阻特性、质心横偏量的综合设计问题,提出了多点多目标优化设计数学模型。通过多目标优化设计方法结合并行数值模拟技术,对该多点多目标气动外形优化设计问题进行研究,为了加快多点数值计算进度,采用了嵌套并行方法,通过有效利用硬件资源来提高多个状态气动数值求解效率。根据以上方法给出的最优设计边界指出了返回舱单点静稳定性与配平升阻比和质心横偏量的矛盾关系,改善单点静稳定性会导致配平升阻比下降,使质心横偏量增加;反之,配平升阻比增加,质心横偏量减少都会使单点静稳定性变差。     

微型凸起物高超声速气动热特性研究

利用N-S数值解方法对高超声速气流中微型凸起物的气动热环境进行了计算研究,并对前斜坡峰值压力和热流跃升进行了分析,结果表明对于高雷诺数流动,由于其边界层很薄,气流对凸起物的冲击作用比低雷诺数流动大很多,给凸起物前缘和上表面带来较严重的气动热问题。研究还表明利用斜面侧扩张角可以降低凸起物后斜坡的气动热环境。最后对球锥与平板物型的差异给微型凸起物热环境特性带来的影响进行了研究,结果表明在本文计算状态下,风洞平板实验模型得到的凸起物前缘峰值压力和热流跃升比实际球锥外形要高,因此地面实验结果外推到天上实际情况还是比较保守的。     

纳米-微米纤维复合形式对惯性粒子过滤性能影响的数值研究

考虑纳米纤维表面流体滑移效应,在二维简化模型下采用数值方法研究了纳米-微米复合纤维对惯性粒子的捕集行为.分析了纳米-微米纤维复合形式对粒子捕集效率、效率增长比例因子以及综合过滤性能的影响.结果表明:在微米纤维迎风面45°和90°处放置纳米纤维,可显著增大复合纤维对粒子的捕集效率,并对弱惯性或中等惯性粒子捕集表现出较高的...     

基于广义预测控制的网络控制系统设计

网络控制系统中,网络诱导时延和数据包丢失会导致系统控制性能的下降甚至不稳定。针对上述问题,提出了一种基于广义预测控制的控制方法。该方法采用时间戳机制和模型的在线辨识,能够根据测量系统输出及时对预测值进行在线修正。在执行器端设计了网络延迟补偿器,可根据当前时刻的最新控制信号对网络系统实现有效控制。仿真结果验证了该方法的有效性。     

轴对称双喉道流体控制矢量喷管三维数值模拟

对轴对称双喉道流体控制矢量喷管的流场进行了数值模拟,研究了喷管气动参数(主、次流落压比)对其内特性(流量系数、推力系数、推力矢量角和推力矢量效率)的影响,并与试验数据进行了对比.研究结果表明,在计算工况下,当喷管次流流量比恒定时,喷管推力矢量角随主流落压比的升高而降低,喷管流量系数和推力系数先是随主流落压比的升高而升高,在某-主流落压比时达到最大值,而后喷管推力系数随主流落压比的升高逐渐降低,流量系数则基本维持不变.当喷管主流落压比固定时,随着次流落压比的升高,推力矢量角也随之增加,推力系数无明显变化,而流量系数则呈下降趋势.