热化学非平衡流场辐射特性数值计算

设计了耦合辐射的三维热化学非平衡流场计算方法,以适用于任意形状的网格。采用Jameson有限体积法求解耦合辐射源项的三维N-S方程。化学模型包含11个组元,20个化学反应。辐射源项通过直接求解辐射输运方程RTE获得。在空间和方向上分别离散后,利用有限体积法求解不同方向上的辐射输运方程,精确计算了波长为100-1500nm所有主要组元的吸收系数,逐条计算了目前已观测到的空气组元的原子谱线,而对分子的带状谱采用了带模型计算,并最终得到了驻点处辐射密度随波长变化的分布情况。     

某支线客机总体方案中增升装置的设计与优化

针对某支线客机,从工程设计的角度,完成了该客机增升装置的气动优化设计。首先根据飞机的相关参数确定增升装置所要达到的气动目标,设计增升装置在机翼的展向和弦向的类型和布置方式,从而确定增升装置的总体设计方案。在二维增升装置设计中,使用NURBS曲线参数化拟合多段翼型缝道部分的外形,通过数值模拟方法和优化算法对三段翼型的缝道外形和缝道参数进行优化设计,最终确定了三段翼型起飞、着陆构型,提出一种二维三段翼型(前缘缝翼、主翼、后缘襟翼)的参数化优化设计方法。在三维增升装置设计中,根据增升装置在机翼展向和弦向的布置方式,取若干关键截面,按照三段翼型的优化结果生成三维增升装置。对三维增升起飞构型和着陆构型的气动特性分析评估表明,所设计的三维增升装置满足设定的气动目标。     

用于弹性旋翼流场模拟的网格变形方法研究

建立了一个基于改进"Ball-Vertex"的具有拓扑无关性的网格变形方法,用来精确模拟桨叶发生弹性变形后贴体网格的运动以及旋翼近场气动力环境,为更好地分析弹性变形对桨叶贴体网格的影响提供了一套解决方案.该网格变形方法从基本的弹簧模拟思路出发,用一网状弹簧系统来模拟桨叶贴体网格的变形和运动.在网格结点周围引入"Ball-Vertex"概念,通过添加相关的冗余约束,该方法保证网格结点在变形过程中不会越过与该点所对应的对角面,从而避免了畸形网格单元的产生,可有效地用于较大变形情况下的流场数值计算.使用所建立的方法,首先对绕二维翼型的结构网格和非结构网格进行了网格变形模拟,以表明方法的有效性;然后,将该方法应用于振荡翼型的流场数值求解,进一步验证了所建立的方法;最后,结合该网格变形方法和旋翼流场求解方法对弹性旋翼流场进行了数值计算.在此基础上,得出了一些有意义的结论.     

基于上游渐扩管安装条件的内锥流量计性能预测

针对内锥流量计使用灵活性要求,利用计算流体动力学数值仿真和实流实验相结合的方法,研究上游渐扩管安装条件对内锥流量计性能的影响,以获取所需的最短直管段长度.研究对象是lOOmm口径、β值分别为0.45,0.65,0.85三种结构类型的样机.开展了基线和渐扩管两种类型的实验,仿真和实验的介质均为常温水,雷诺数范围分别为0.2488×105~2.488×105和0.3843×105~2.479×105,仿真结果和实验结论一致.利用附加不确定度和流出系数相对误差作为主要的评价标准,给出了上游渐扩管安装条件内锥流量计所需的直管段长度.     

弹尾超声速轴对称喷流的正格式数值模拟(英文)

采用二阶正格式方法对超声速轴对称喷流流动进行数值模拟。将二维守恒方程的正格式方法发展到轴对称Euler方程组的求解,并对导弹尾部超声速伴随射流进行了数值计算。数值结果与实验照片所反映的流动特征吻合较好,与三阶精度间断有限元方法计算结果吻合,相比于二阶高分辨率TVD格式的计算结果,正格式方法的计算结果在间断点处具有较大的梯度变化。这表明该方法对激波具有较强的捕捉能力,在间断处不会出现数值振荡,对弹尾声音速伴随射流的数值模拟真实、有效。     

螺旋桨安装效应对无人机气动特性影响

为解决某型飞翼布局无人机(UAV)带动力构型风洞试验最大升阻比相对无动力状态大幅下降的问题,采用计算流体动力学(CFD)方法对无人机无动力与带动力构型进行了数值模拟,数值模拟结果分别与无动力以及带动力风洞试验数据吻合良好,在此基础上深入研究了螺旋桨安装效应对无人机气动特性的影响。结果表明:推力螺旋桨与机身之间气动干扰产生的低压区致使阻力增加,从而导致飞机最大升阻比相比无动力状态下降了30.7%。针对无人机在推力螺旋桨影响下出现的最大升阻比下降问题,采用增大螺旋桨与机身之间距离的方法可以有效地消除机身后部出现的低压区,减小了阻力,提升了无人机最大升阻比。桨毂拉长方案在8°和9°迎角下最大升阻比分别提升了17.3%和15.4%。     

计算域大小对大涡模拟计算结果的影响

为了节省计算资源,提升计算效率,在保证计算结果准确的前提下,尽可能选取较小的计算域显得尤为重要。针对该问题,分别沿槽道的展向和流向选取不同的计算尺度(053πδ×2δ×0177πδ、πδ×2δ×05πδ、2πδ×2δ×πδ、3πδ×2δ×15πδ,δ为槽高的一半),在相同的计算条件且设置均满足大涡模拟计算要求的前提下,对摩擦雷诺数为540和930两种摩擦雷诺数下计算域大小对大涡模拟数值计算结果的影响进行了研究。结果发现:计算尺度不小于2πδ×2δ×πδ的计算域能很好的呈现湍流计算结果。此外,在摩擦雷诺数为930时,相同计算域下所能检测的范围均有所提升,得到的计算结果更加准确。     

任意多面体网格上的欧拉方程数值算法

发展并验证了一种新的支持多面体网格的欧拉方程离散算法,采用Fortran 95编写了支持任意网格拓扑及三维复杂外形的通用求解器.在空间离散上,基于径向基函数理论提出了一种新的梯度计算方法,并采用动能通量分裂格式来得到准确并且稳定的二阶精度重构.该方法不依赖于网格几何形状,因此对网格质量不敏感.由于在时间方向采用了点隐松...     

Experimental and Numerical Investigation of the Effects of Passive Vortex Generators on Aludra UAV Performance

A study of the effects of passive vortex generators (VGs) on Aludra unmanned aerial vehicle (UAV) aerodynamic characteristics is presented. Both experimental and numerical works are carried out where an array of VGs is attached on Aludra UAV’s wing. The flow measurements are made at various angles of attack by using 3-axis component balance system. In the numerical investigation, the Reynolds-averaged Navier-Stokes (RANS) code FLUENT 6.3TM is used in the simulations with fully structured mesh with Spalart-Allmaras (S-A) turbulence model and standard wall function. The comparison between the experimental and numerical results reveals a satisfactory agreement. The parametric study shows that higher maximum lift coefficient is achieved when the VGs are placed nearer to the separation point. In addition to this, shorter spanwise distance between the VGs also increases the maximum lift coefficient, rectangular and curve-edge VG performs better than triangular VG.     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...