基于数值模拟的轴对称矢量喷管性能预测数学模型

基于三维数值模拟,对不同轴对称矢量喷管在多种工作状态下的内外流场进行了研究,分析了扩张角及扩张段长度对喷管有效矢量角的影响.基于最小二乘曲面拟合理论,建立了自变量包括扩张角、扩张段长度/喉道直径、落压比/设计落压比、几何偏转矢量角的多变量轴对称矢量喷管性能预测数学模型,并根据已有实验数据,对该模型进行了验证结果表明:推力系数误差最大为0.41%,流量系数误差最大为1.58%,矢量角误差最大为1.76°.建立的数学模型通用性较强,实现了用统一的模型对不同喷管性能参数进行预测和分析,具有一定的工程意义.      

2099铝锂合金型材热压下陷模具结构优选

2099铝锂合金型材热压下陷成形在飞机零件制造中应用广泛,下陷成形时容易出现温度分布不均、型材上表面回弹不一致以及裂纹等缺陷,模具结构对其下陷零件的质量具有较大影响.通过热力耦合有限元数值模拟预测2099铝锂合金型材热压下陷成形过程,分析上述缺陷产生的原因,对模具加热结构、测压块和固定下模进行优选.有限元数值模拟和试验结果表明,优选后的模具能够较好地消除或缓解上述缺陷,提高了下陷零件的质量.     

小型垂直轴风力机叶片结冰风洞试验与数值计算

利用风洞试验和数值模拟相结合的手段,研究了小型垂直轴风力机叶片在旋转状态下的结冰特性以及结冰后翼型与风力机的气动特性变化,以期为建立较复杂的大型水平轴风力机叶片旋转试验系统和研究其结冰机理、防除冰技术提供参考。试验在东北农业大学自行设计的利用自然低温的结冰风洞中进行,获得了采用NACA0018翼型的小型2叶片垂直轴风力机风轮在5种尖速比下的结冰分布：风力机叶片结冰遍布叶片整个表面,随着尖速比的增大,结冰形状出现不对称性。同时,数值模拟结果表明：叶片结冰后,随着尖速比的增加和结冰量的增多,升力系数降低阻力系数增大的趋势明显,风力机的功率系数也随之下降。分析发现,叶片结冰导致不同旋转角下叶片翼型周围的压力场和速度场发生了不同程度的变化,从而气动特性发生变化,影响了风力机性能。     

基于涡限制法的旋翼气动噪声计算

在直升机旋翼流动的数值模拟中,流动求解精度对旋翼气动噪声的预测有着重要影响.采用中心格式有限体积法求解可压缩Euler方程数值模拟旋翼绕流,并采用基于Ffowcs Williams- Hawkings( FW-H)方程的声类比方法计算旋翼气动噪声.为了抑制中心格式的数值耗散,保持机翼翼尖以及旋翼桨尖的尾涡结构,引入了涡...     

三种典型翼型边界层稳定性对比分析

为了便于工程上翼型的选取及优化,本文针对三个典型翼型边界层稳定性进行对比分析。首先采用eN方法对翼型SD8020进行转捩预测,并采用萘升华实验法检验数值计算的准确性。再利用该数值方法对比了三种典型翼型上表面的转捩位置、不同弦长雷诺数下的压力系数分布、扰动增长率以及最不稳定扰动波的频率。分析结果表明,在相同弦长雷诺数下,NACA0012最先发生转捩,而NACA64-204和RAE2822都保持着较长的层流区;压力梯度对扰动增长有很大影响;同一翼型随弦长雷诺数增长,转捩点雷诺数变大,而不是工程上常采用的固定值。     

旋翼翼型多目标多约束气动优化设计

为了克服传统旋翼翼型优化设计方法的不足,发展了一种基于Kriging模型与遗传算法的旋翼翼型多目标多约束气动优化设计方法。采用基于雷诺平均Navier-Stokes方程的数值模拟获得样本翼型气动性能,并建立目标函数和状态函数的Kriging模型,采用遗传算法搜索Kriging模型最小值和相应的EI(Expected Improvement)函数最大值,更新Kriging 模型直至找到满足约束的最优翼型。运用加权目标函数法进行了旋翼翼型的多设计点优化设计。优化结果表明,优化后旋翼翼型在满足约束的同时,与基准旋翼翼型OA209相比:在悬停状态下,阻力系数下降了2.1%;在机动状态下,最大升力系数提高了4.2%;在前飞状态下,阻力系数在不同马赫数下均有所减小。     

辅助进气激励器强化喷流混合的数值模拟

对采用高动量的辅助进气激励器和“常规”激励器强化的高速、大尺度喷流混合进行了数值模拟研究。对比了辅助进气激励器和“常规”激励器吹气射流对流向涡、喷流拍动、阻塞面积的影响,分析了合成射流强化喷流混合的流动滞后现象。结果表明:合成射流强化喷流混合中的辅助进气激励器在吹气时间、吹气动量、峰值动量比方面要明显优于“常规”激励器。合成射流强化喷流混合由激励器迟滞和流场迟滞两部分组成。“常规”激励器和辅助进气激励器射流峰值动量都滞后于活塞峰值速度时刻19.5%,流场滞后峰值动量时刻5%左右。相比于“常规”激励器,采用辅助进气激励器的“小突片”堵塞面积较大,强化喷流混合的效果更好。      

渗透壁面竖直平板Blasius流熵产特性

对渗透壁面竖直平板Blasius流层流边界层流动熵产及传热熵产与总熵产的比值(Bejan数)进行了深入研究,将耦合的动量与能量偏微分方程组通过相似变换转换为非线性常微分方程组,用Runge-Kutta法求解获得了无量纲流函数与无量纲温度的相似解,进而研究了无量纲参数对熵产及Bejan数的影响.结果表明:吸入时最大熵产位于壁面处,吸入速度越大最大熵产越大,熵产随相似变量的增加单调下降;喷注速度越大壁面处熵产越小;近壁面处熵产随变黏度参数增大而增大,远离壁面处随喷注/吸入参数或变黏度参数的增大而减小;远离壁面处Bejan数随喷注/吸入参数或变黏度参数的增大而增大,近壁面处与此相反;熵产和Bejan数随着辐射参数和Reynolds数增大而减小,随着Biot数增大而增大;Biot数在0~2范围内对传热熵产的影响很剧烈,Biot数较小时总熵产由流动熵产控制,Biot数较大时总熵产由传热熵产控制;熵产随滑移参数的减小而增大,Bejan数随滑移参数的增大而增大;Grashof数对熵产的影响较大,在较大的Grashof数条件下熵产随相似变量的变化较剧烈;滑移参数越大,流动熵产越小,Bejan数越大.      

二次空气系统盘腔一维非定常计算方法

通过对盘腔结构的一维算法展开研究,提出适用于具有非标准结构、包含多个进出口的复杂盘腔的一维建模规则,并采用特征线法,实现盘腔结构的一维非定常流动与换热的计算,通过将广义的变截面、有摩擦、考虑对外热交换的盘腔一维可压缩非定常流动与换热的偏微分方程组转换为常微分方程组,进而获得盘腔内部沿流动方向压力、流量、温度等系统参数的特征线数值解,采用商业计算流体软件CFX对其进行验证.结果表明:所提出的盘腔建模与计算规则能够直观且准确地描述盘腔一维非定常流动及换热特征,仿真结果与CFX计算结果最大偏差约为9.4%,可应用于具有非标准结构的多接口复杂盘腔的模拟,可以作为二次空气系统盘腔结构一维非定常计算研究的有效手段.      

柔性扑翼全机气动结构耦合特性及稳定性研究

由于微型扑翼飞行器具有体积小、重量轻等特点,在飞行过程中会产生明显的柔性变形,因此在对扑翼气动特性进行数值模拟时,有必要考虑柔性变形的影响。基于结构动力学理论,发展了一种适用于扑翼全机气动结构耦合特性的数值模拟方法与一种适用于扑翼的结构动力学求解方法。方法利用哈密顿原理,对动能和应变能进行变分,进而得到扑翼动力学方程,采用结构有限元方法对运动方程进行离散并求解;采用基于嵌套网格的CFD求解器对微型扑翼全机非定常绕流进行数值模拟。采用CFD/CSD耦合求解器对微型扑翼飞行器全机气动结构耦合特性进行数值模拟,分析了结构变形对柔性扑翼气动特性的影响,并分别对扑动过程中刚性扑翼和柔性扑翼的压心变化范围进行研究,分析了结构变形对柔性扑翼稳定性的影响。