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

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

NS方程中耦合转捩判断的风力机翼型气动力数值模拟

在Reynolds-Averaged Navier-Stokes(RANS)方程计算中耦合Michel经验判据和Chen-Thy-son转捩模型对风力机翼型S809进行计算。由RANS方程求得翼型表面压力分布作为层流边界层方程求解的输入参数,然后使用Michel经验判据分析判断层流边界层的解得到转捩点的位置,使用Chen-Thyson转捩模型得到从层流到湍流的转捩过渡区,这样随着流场的迭代求解,求解器自动更新判断出转捩点位置。对S809翼型进行全湍流和耦合转捩判断的RANS方程计算,可以看到考虑转捩判断后得到的升阻力系数与实验值吻合较好,验证了方法的可行性。     

耦合Navier-Stokes方程的可压缩边界层数值求解

为了在Reynolds-averaged Navier-Stokes(RANS)方程计算中耦合eN方法进行转捩判断,在RANS方程求解过程中耦合求解了可压缩层流边界层方程为判断转捩提供了精确、可靠的边界层信息.利用等熵关系由RANS方程求出的物面压力分布确定边界层外边界的速度分布,进一步确定出边界层外边界.边界层方程的求解使用Keller提出的二阶BOX方法.为了验证方法求解边界层方程的正确性,在低速流动状态下将计算结果和XFOIL的计算得到的边界层解进行了比较,二者吻合较好.     

NS方程计算中耦合转捩自动判断的阻力精确计算方法初探

在Reynolds-Averaged Navier-Stokes(RANS)方程计算中耦合了流动转捩的自动判断以提高现有求解器预测翼型阻力的准确性.由RANS方程求得翼型表面压力分布作为层流边界层方程求解的输入参数,然后使用简化的eN-数据库转捩判断方法分析层流边界层的解得到转捩点的位置,这样随着流场的迭代求解求解器自动判断转捩点的位置.在对NLF0416翼型的气动性能计算中考虑流动转捩的因素后得到的翼型升阻力特性和实验吻合较好,验证了本文方法的正确性.     

一种基于动模态分解的翼型流动转捩预测新方法

考虑自由转捩的定常/非定常流动Navier-Stokes方程数值求解,对于翼型流动细节的精确模拟和气动力的精确预测均具有十分重要的意义。采用动模态分解(DMD)方法进行流动稳定性分析,再结合e N方法,提出了一套适用于翼型绕流的转捩预测新方法,称为DMD/e N方法。相比于传统的线性稳定性分析方法,DMD方法不需要求解附面层方程和线性稳定性方程,也没有引入平行流假设,具有更好的理论适用性和算法鲁棒性。开展了NLF0416、S809和SD7003等翼型的转捩预测数值验证研究,通过与实验结果以及与传统的基于线性稳定性分析的e N方法的比较,验证了本文所发展的转捩预测新方法在预测翼型的定常流动和非定常流动转捩方面的正确性,也表明了该方法具有解决含层流分离泡的翼型绕流转捩预测的能力。     

抛物化稳定性方程在曲面边界层中的应用

为了研究抛物化稳定性方程在曲面边界层中的应用,选取了等曲率圆环管道、等曲率板和NACA0012翼型3种典型的曲面边界层.通过计算小幅值扰动波的演化,抛物化稳定性方程的计算结果和线性稳定性理论、数值模拟扰动方程结果相符,说明可以使用抛物化稳定性方程研究曲面边界层的小幅值扰动波的演化和稳定性分析;通过计算有限幅值扰动波的演化,线性阶段抛物化稳定性方程计算结果和扰动方程相符,在非线性很强的阶段,抛物化稳定性方程计算发散,发散的位置可作为转捩位置.      

耦合转捩自动判断的机翼粘性绕流计算与优化设计

具有边界层转捩自动判断功能的三维雷诺平均Navier-Stokes方程(RANS)方程求解器是自然层流机翼设计的必要手段。将一种基于线性稳定性理论的NTS/NCF转捩判断方法与RANS求解器耦合,以提高求解器的计算精度。通过求解三维稳定性方程,然后使用包络线方法计算对应Tollmien-Schlichting不稳定性的NTS,使用固定β方法计算对应横流不稳定性的NCF,然后根据转捩阀值(NTS)tr和(NCF)tr判断转捩位置。在上述工作的基础上,开展了自然层流机翼平面形状优化的气动减阻设计研究。     

湍流附面层方程的有限差分解

本文用三点隐式有限差分法数值求解了轴对称流动的层流、转捩和湍流可压缩附面层方程,湍流附面层用两层概念处理,而每一层用适当的涡粘性模型代替Reynolds应力项,用常值湍流Prandtl数将湍流热流项与Reynolds应力关联起来。支配方程用迭代法求解。计算结果表明,采用本文改进的涡粘性模型和涡粘性分布的平滑办法,对文中所计算的几种情况,都给出了满意的结果。     

压缩性修正对γ-Re_θ转捩模型的影响研究北大核心CSCD

将基于流场当地变量的γ-Reθ转捩模型加入可压缩RANS方程计算程序,对平均流场控制方程、湍流和转捩控制方程进行基于LU-SGS的隐式紧耦合求解。结合三种压缩性修正方法对高超声速平板、双楔和圆锥流动进行了数值模拟,给出了壁面斯坦顿数、热流值和实验值的比较,以及湍动能和间歇因子分布云图。数值计算结果表明,相同来流湍流度下不同压缩性修正方法得到的转捩位置差别较大。相比于未修正的模型,基于当地马赫数压缩性修正后的转捩位置最靠后,其对带有分离的双楔流动预测较为准确;基于来流马赫数的压缩性修正仅使得转捩位置稍有延迟;对湍动能源项的压缩性修正可提高模型在转捩后湍流段的热流预测精度。     

基于eN-数据库方法复杂构型飞机转捩预测

为探索边界层转捩对大型运输机在起降条件和有较大层流区的巡航条件下的气动力精确计算问题，通过在三维RANS求解器中引入eN-数据库方法来预测飞行器表面的转捩位置，并探索转捩对气动力的影响规律。方法与目前流行的基于间歇因子控制方程的转捩预测方法相比，具有计算效率高、易于工程应用、且考虑TS不稳定性转捩因素的特点。在此基础上，通过计算NASA梯形翼来分析起降构型条件下气动力受转捩影响的规律，并通过计算DLR－F6翼身组合体来探索三维构型在巡航条件下的气动力精度。使用eN-数据库转捩判断方法的计算结果与实验值吻合较好，验证了所构建的基于RANS求解器的eN-数据库转捩预测方法的有效性，并为大型运输机气动力精确计算提供了分析工具。     

大厚度钝后缘翼型气动性能计算研究

大厚度钝后缘翼型由于其气动和结构的优点,近来被作为大型风力机叶片设计时的内侧翼型成为风力机翼型设计的热点之一EULL由于此类翼型厚度和风洞阻塞度的限制,大雷诺数的风洞实验数据很少。耦合RANS方程和基于线性稳定性分析EULL捩预测方法,进行了大厚度钝后缘翼型气动性能计算研究。使用Gridgen程序生成计算网格,湍流模型为S-A模型。对常用的钝后缘修形方式如直接截断、对称增加厚度等并且有风洞实验结果的几种翼型进行了CFD计算并和实验结果及文献计算结果进行了比较。结果表明升力系数在线性段,计算结果和实验结果吻合很好,但是对失速迎角的捕捉能力差;阻力系数计算结果和实验值吻合较好。不同修形方式的计算结果为钝后缘修形设计提供了参考。     

低雷诺数翼型和跨声速翼型边界层计算

本文利用边界层动量积分方程和平均流动能积分方程兼容计算了翼面的层流和湍流边界层流动。文中采用一种e~9型转捩判别公式预测翼面存在层流分离气泡情况的转捩位置。并引入对剪应力系数C_r的滞后方程以体现湍流中雷诺应力的历程效应。为避免在翼面流动分离时边界层方程的奇异,采用反解法计算。算例表明,计算与测量结果吻合良好。     

高超声速边界层的转捩问题

取40km高空处的气体参数,对来流马赫数为8、10和12等三种不同情况的平板边界层和楔角为20°的楔型体边界层,做了转捩发生位置的预测。所用方法为改进的eN方法。结果发现,壁面条件对是第一还是第二模态波决定转捩位置有很大关系。如对于传统的eN方法,如果取N为10作为转捩判断标准,对马赫数为12的楔体边界层,等温壁条件下由第二模态波决定的转捩位置距前缘13m左右,而平板边界层中两模态波所决定的转捩位置距前缘均超过50m。而采用改进的eN方法,平板边界层中两模态波所决定的转捩位置距前缘均超过10m,而对马赫数为12的楔体边界层,两模态波所决定的转捩位置距前缘均出现小于5m的情况。     

高超声速边界层感受性问题高精度数值模拟

通过五阶WENO格式和六阶对称紧致格式以及三阶TVD R-K结合的方法,对存在强激波和小扰动相互干扰的高超声速边界层感受性问题进行研究.结果表明:此方法能够模拟边界层内不稳定波的产生和发展以及小扰动和激波相互干扰等现象,因此能广泛应用于含激波的感受性等问题的可压缩湍流直接数值模拟以及具有激波和边界层干扰等复杂流场的计算.     

Application of focused laser differential interferometer to hypersonic boundary-layer instability study

Accurate prediction of hypersonic boundary-layer transition plays an important role in thermal protection system design of hypersonic vehicles. Restricted by the capability of spatial diagnostics for hypersonic boundary-layer study, quite a lot of problems of hypersonic boundary-layer transition, such as nonlinearity and receptivity, remain outstanding. This work reports the application of focused laser differential interferometer to instability wave development across hypersonic boundary-layer on a flared cone model. To begin with, the focused laser differential interferometer is designed and set up in a Mach number 6 hypersonic quiet wind tunnel with the focal point in the laminar boundary-layer of a 5 degree half-angle flared cone model. Afterwards, instability experiments are carried out by traversing the focal point throughout the hypersonic boundary-layer and the density fluctuation along the boundary-layer profile is measured and analyzed. The results show that three types of instability waves ranging from 10 kHz to over 1 MHz are co-existing in the hypersonic boundary-layer, indicating the powerful capability of focused laser differential interferometer in dynamic response resolution for instability wave study in hypersonic flow regime; furthermore, quantitative analyses including spectra and bicoherence analysis of instability waves throughout the hypersonic boundary-layer for both cold and heated cone models are performed.     

Instability analysis and drag coefficient prediction on a swept RAE2822 wing with constant lift coefficient

Swept wing is widely used in civil aircraft,whose airfoil is chosen,designed and optimized to increase the cruise speed and decrease the drag coefficient.The parameters of swept wing,such as sweep angle and angle of attack,are determined according to the cruise lift coefficient requirement,and the drag coefficient is expected to be predicted accurately,which involves the instability characteristics and transition position of the flow.The pressure coefficient of the RAE2822 wing with given constant lift coefficient is obtained by solving the three-dimensional Navier-Stokes equation numerically,and then the mean flow is calculated by solving the boundary layer(BL) equation with spectral method.The cross-flow instability characteristic of boundary layer of swept wing in the windward and leeward is analyzed by linear stability theory(LST),and the transition position is predicted by eNmethod.The drag coefficient is numerically predicted by introducing a laminar/turbulent indicator.A simple approach to calculate the lift coefficient of swept wing is proposed.It is found that there is a quantitative relationship between the angle of attack and sweep angle when the lift coefficient keeps constant;when the angle of attack is small,the flow on the leeward of the wing is stable.when the angle of attack is larger than 3°,the flow becomes unstable quickly;with the increase of sweep angle or angle of attack the disturbance on the windward becomes more unstable,leading to the moving forward of the transition position to the leading edge of the wing;the drag coefficient has two significant jumping growth due to the successive occurrence of transition in the windward and the leeward;the optimal range of sweep angle for civil aircraft is suggested.     

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

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

基于γ-Reθt转捩模型的翼型气动特性计算精度改进研究

将结合SSTk-ω湍流模型的γ-Reθt转捩模型与二维雷诺平均NS（RANS）方程数值求解相耦合,发展了考虑转捩影响的翼型气动特性计算方法。得到的计算结果与没有考虑转捩影响的计算结果的比较表明,与实验更加吻合,计算精度得到了明显提高。