旋翼非定常平行桨-涡干扰流场的数值模拟

建立了一个基于Euler方程的旋翼平行桨-涡干扰流场的数值模拟方法。该方法在时间离散上使用隐式的双时间法,以避免因物理时间步长过小而引起的数值不稳定问题;为减小由于空间离散格式精度和网格密度引起的涡数值耗散,采用了预定涡方法。以Kitaplioglu的平行桨-涡干扰试验模型为算例,计算了表面瞬时压强分布,并与可得到的试验结果相比较,验证了方法的有效性,且描述了干扰流场的气动力变化特性以及桨叶表面扰动压力波的产生和传播过程,然后分析了马赫数、涡强度、干扰距离等参数对桨-涡干扰流场特性的影响。在此基础上,得出了一些有意义的结论。     

湍流模型对喷流计算精度的影响研究(英文)

采用三维RANS代码MBNS3D对包含了内流、外流及后体干扰的三种简化尾喷流,进行了数值模拟研究。通过壁面压力分布及喷流截面速度型与实验数据的对比,研究了SA和SST两种湍流模型及其修正(SACC,SADES,SSTCC)对喷流计算精度的影响。本文计算表明,SSTCC模拟效果最为精确;利用SA模型进行脱体涡模拟,在同样的网格下优势不明显。湍流模型可压缩修正改善了模拟喷流的膨胀趋势,但对附面层分离位置的预测,在本文的工作中,还不如原始模型。两类模型相比,SST模型的结果好于SA模型。     

基于响应面法的低速翼型气动优化设计

响应面方法相较于其它直接优化方法有其高效、实用的优势,此前的研究更多地将响应面方法用于超音速和跨音速翼型的减阻优化设计中。本文将此方法应用于低速翼型优化设计中,进行了基于RANS(Reynolds-Aver-aged Navier-Stokes)方程和自由转捩预测耦合求解的低速翼型气动优化设计。通过计算附面层方程得到附面层参数并用en方法计算转捩位置,并考虑了T-S波和层流分离造成的转捩。RANS方程计算中,使用了转捩过渡区模型,以保证附面层外边界压力分布的精度。RANS方程和转捩预测迭代进行至转捩位置收敛。在响应面模型计算中,使用不含二阶交叉项的二阶多项式模型,减少了构造模型所需的计算量;合理的选择设计空间保证了构造的响应面模型具有较高的精度。使用三个设计点的多目标优化设计,保证了设计的合理性。通过对NACA64(1)-112翼型优化计算结果表明,本文的方法可以有效地进行低速翼型的气动优化,各设计点上转捩位置也得到了改善,有较好的工程实用前景。     

含扰流板结构的固液火箭发动机数值模拟

为进一步提高固液火箭发动机的燃烧效率,在FLUENT软件平台上,利用二维轴对称的N-S方程和组分方程,采用有限速率化学反应模型和S-A单方程湍流模型,对药柱和后燃室中添加不同数量和位置的扰流板对燃烧室和喷管温度和效率的影响进行了数值模拟研究.在数值模拟中,对流场进行假设,假设流动为纯气相流动,燃烧室中气体为理想气体.数值模拟结果表明,固液火箭发动机在纯气相的反应条件下,在固体药柱中添加扰流板可以提高燃烧效率,但提高的程度有限,且不能改变喷管入口处温度分布不均匀的问题;在后燃室中添加扰流板可以明显地提高喷管入口处的平均温度,而且温度分布基本均匀.由于固液火箭发动机燃烧的特点是反应发生在燃料表面上的边界层中的火焰层中,在固体药柱中添加扰流板可以改变火焰层的位置,在后燃室中添加扰流板后,由于火焰层位置相对固定,所以反应开始时和反应进行中,喷管入口处的温度分布没有发生太大的变化,可以使固液火箭发动机维持一个相对稳定的燃烧情况.     

电脉冲除冰系统电磁脉冲力仿真分析

为了获得更精确的电脉冲除冰(EIDI)系统中电磁力在飞机蒙皮上的分布及其随时间变化的情况,针对系统的工作原理与工作过程进行了仿真研究。计算分析了EIDI系统电路中电流随时间变化的关系;通过麦克斯韦方程组的求解完成了EIDI系统电磁场的仿真;计算分析了电磁力随时间变化的规律,获得了蒙皮上的磁感应强度分布、涡电流密度分布和瞬态电磁力密度分布。考虑线圈工作时蒙皮对其电感的影响,建立了更接近真实电路的激励加载方法,得到了准确的电磁力变化趋势;考虑趋肤效应对蒙皮厚度方向上涡电流密度、电磁力密度分布的影响,用有限元分析方法使蒙皮上电磁力分布的计算更为全面合理。      

基于非结构网格的跨声速静气动弹性计算

基于非结构网格开发了静气动弹性计算程序,通过气动力与弹性变形的交替迭代计算实现流固耦合。气动力求解采用RANS方程,结构变形计算采用柔度系数法,通过面样条(IPS)方法进行柔度系数插值,使用改进弹簧近似法实现网格变形。通过对AGARD 445.6机翼的静气弹计算验证了方法的正确性,同时对机翼的静气动弹性特性进行了分析,计算表明,弹性变形使其升力线斜率下降、焦点前移。     

某飞翼布局隐身飞行器的翼型优化

对某飞翼布局无人机的翼型进行考虑隐身迎角的气动、隐身综合优化研究。采用基于非均匀有理B样条曲线的自由变形方法实现翼型的参数化表达;无限插值动网格方法生成气动求解的计算网格;分别采用基于雷诺平均N—S方程求解的计算流体力学方法及物理光学法计算翼型的气动、隐身性能;选择多目标遗传算法进行综合寻优设计研究。优化翼型兼顾了气动与隐身的要求,其综合性能优于原翼型。结果表明,方法具有较好的实用性。     

单相机迎角测量中振动影响研究

从描述像空间坐标与物理空间坐标关系的共线方程出发,推导出单相机迎角测量中模型迎角关于相机外方位元素的计算公式。在单相机迎角测量中模型表面标记点展向坐标Y不变的假设下,研究了2m量级风洞中,相机位于5个不同位置时,3个角外方位元素(φ,ω,κ)和模型侧滑角β在-2°～2°之间振动对迎角测量结果的影响。模拟结果表明:φ和κ的振动对迎角测量结果影响较大,且该影响对相机位置和模型迎角的变化不敏感;在相机位于法向投影位置(相机光轴与Y轴平行)时,ω与β的振动对迎角测量结果影响非常小,可以忽略,但影响量随相机位置偏离法向投影位置迅速增加。     

Computational aerodynamic analysis on perimeter reinforced （PR）-compliant wing

Abstract Implementing the morphing technique on a micro air vehicle （MAV） wing is a very chal- lenging task, due to the MAWs wing size limitation and the complex morphing mechanism. As a result, understanding aerodynamic characteristics and flow configurations, subject to wing structure deformation of a morphing wing MAV has remained obstructed. Thus, this paper presents the investigation of structural deformation, aerodynamics performance and flow formation on a pro- posed twist morphing MAV wing design named perimeter reinforced （PR）-compliant wing. The numerical simulation of two-way fluid structure interaction （FSI） investigation consist of a quasi- static aeroelastic structural analysis coupled with 3D incompressible Reynolds-averaged Navier- Stokes and shear-stress-transport （RANS-SST） solver utilized throughout this study. Verification of numerical method on a rigid rectangular wing achieves a good correlation with available exper- imental results. A comparative aeroelastic study between PR-compliant to PR and rigid wing per- formance is organized to elucidate the morphing wing performances. Structural deformation results show that PR-compliant wing is able to alter the wing＇s geometric twist characteristic, which has directly influenced both the overall aerodynamic performance and flow structure behavior. Despite the superior lift performance result, PR-compliant wing also suffers from massive drag penalty, which has consequently affected the wing efficiency in general. Based on vortices investigation, the results reveal the connection between these aerodynamic performances with vortices formation on PR-comoliant wing.     

高铝Ni3Al基合金的热变形行为

采用Gleeble-3800型热模拟试验机研究了高铝Ni3Al基合金在变形温度为1200～1240℃,应变速率为0.01～1s-1条件下的热压缩变形,结果表明:在应变速率为0.01s-1时,高铝Ni3Al基合金对应的热变形本构方程为σ=28.57(lnε+6.72×105/RT-44.08),而当应变速率为0.1s-1和1s-1时,热变形本构方程为σ=28.57(lnε+1.28×106/RT-92.76)。变形过程中只有γ’相发生不同程度上回溶,但未发生动态再结晶。合金的最佳变形区间位于变形温度为1200～1215℃,应变速率为0.01s-1范围内;而当提高速率至1s-1附近,γ’相中塞积的位错容易造成单相γ’区中β/γ’界面的开裂,对应变形过程中的"失稳区"。