Double-decoupled inverse design of natural laminar flow nacelle under transonic conditions

The inverse design based on the pressure distribution is an essential approach to realize the improvement of Natural Laminar Flow(NLF) performance for nacelles. However, the direct definition of target pressure distribution at design point is challenging for the dilemma to consider the constraints of shock wave and laminar flow at the same time. In addition, the universality of method will be limited when the inverse design is strongly coupled with the solver. Thus, a double-decoupled methodolog...     

Conjugate heat transfer investigations of turbine vane based on transition models

The accurate simulation of boundary layer transition process plays a very important role in the prediction of turbine blade temperature field. Based on the Abu-Ghannam and Shaw (AGS) and c-Re h transition models, a 3D conjugate heat transfer solver is developed, where the fluid domain is discretized by multi-block structured grids, and the solid domain is discretized by unstructured grids. At the unmatched fluid/solid interface, the shape function interpolation method is adopted to ensure the conservation of the interfacial heat flux. Then the shear stress transport (SST) model, SST & AGS model and SST & c-Re h model are used to investigate the flow and heat transfer characteristics of Mark II turbine vane. The results indicate that compared with the full turbulence model (SST model), the transition models could improve the prediction accuracy of temperature and heat transfer coefficient at the laminar zone near the blade leading edge. Compared with the AGS transition model, the c-Re h model could predict the transition onset location induced by shock/boundary layer interaction more accurately, and the prediction accuracy of temperature field could be greatly improved.     

Transition prediction and sensitivity analysis for a natural laminar flow wing glove flight experiment

Natural laminar flow technology can significantly reduce aircraft aerodynamic drag and has excellent technical appeal for transport aircraft development with high aerodynamic efficiency. Accurately and efficiently predicting the laminar-to-turbulent transition and revealing the maintenance mechanism of laminar flow in a transport aircraft’s flight environment are significant for developing natural laminar flow wings. In this research, we carry out natural laminar flow flight experiments with different Reynolds numbers and angles of attack. The critical N-factor is calibrated as 9.0 using flight experimental data and linear stability theory from a statistical perspective, which makes sure that the relative error of transition location is within 5%. We then implement a simplified e^N transition prediction method with a similar accuracy compared with linear stability theory. We compute the sensitivity information for the simplified e^N method with an adjoint-based method, using the automatic differentiation technique (ADjoint). The impact of Reynolds numbers and pressure distributions on TS waves is analyzed using the sensitivity information. Through the sensitivity analysis, we find that: favorable pressure gradients not only suppress the development of TS waves but also decrease their sensitivity to Reynolds numbers; there exist three special regions which are very sensitive to the pressure distribution, and the sensitivity decreases as the local favorable pressure gradient increases. The proposed sensitivity analysis method enables robust natural laminar flow wings design.     

模型的转捩流动计算分析简

 为了在黏性流动数值模拟中实现边界层转捩的自动预测,将γ-Re_θt转捩模型引入到三维非结构混合网格的雷诺平均Navier-Stokes方程求解程序（HUNS3D）。该转捩模型由两个依赖当地变量定义的关于间歇因子和当地化转捩起始动量厚度雷诺数的输运方程组成,其数值求解算法与流场求解程序中湍流模型的求解方法相同。为了考察和验证HUNS3D程序中γ-Re_θt转捩模型对航空工程中的常见附面层自由转捩问题的预测精度,对低速平板流动、Aerospatial-A翼型、NLR 7301超临界翼型和NASA Trap wing 高升力构型等典型外形的自由转捩流动进行了计算,并将计算结果与相关试验结果进行了对比分析。算例结果表明：γ-Re_θt转捩模型对于转捩位置具有很好的敏感性,能比较准确地预测自然转捩和分离转捩,可以有效提高HUNS3D程序对实际流动的模拟能力和预测精度。     

Robust design of NLF airfoils

 A robust optimization design approach of natural laminar airfoils is developed in this paper. First, the non-uniform rational B-splines (NURBS) free form deformation method based on NURBS basis function is introduced to the airfoil parameterization. Second, aerodynamic characteristics are evaluated by solving Navier-Stokes equations, and the γ-Re_θt transition model coupling with shear-stress transport (SST) turbulent model is introduced to simulate boundary layer transition. A numerical simulation of transition flow around NLF0416 airfoil is conducted to test the code. The comparison between numerical simulation results and wind tunnel test data approves the validity and applicability of the present transition model. Third, the optimization system is set up, which uses the separated particle swarm optimization (SPSO) as search algorithm and combines the Kriging models as surrogate model during optimization. The system is applied to carry out robust design about the uncertainty of lift coefficient and Mach number for NASA NLF-0115 airfoil. The data of optimized airfoil aerodynamic characteristics indicates that the optimized airfoil can maintain laminar flow stably in an uncertain range and has a wider range of low drag.     

γ-Reθ转捩模型在二维低速问题中的应用

介绍了基于当地关联的γ-Reθ转捩模型.通过Spalart建议的环境源项法对自由流湍流度衰减的有效控制,实现了该转捩模型的实用化.利用中国空气动力研究与发展中心的TRIP软件平台上实现的该转捩模型,对S809翼型、NLR7301两段翼型和30P30N三段翼型3个低速翼型绕流问题进行了计算,将不同迎角下的计算结果与试验结...     

Effect of scavenge port angles on flow distribution and performance of swirl-loop scavenging in 2-stroke aircraft diesel engine

Swirl-Loop Scavenging (SLS) improves the performance of 2-stroke aircraft diesel engine because the involved swirl may not only benefit the scavenging process, but also facilitate the fuel atomization and combustion. The arrangement of scavenge port angles greatly influences in-cylinder flow distribution and swirl intensity, as well as the performance of the SLS engine. However, the mechanism of the effect and visualization experiment are rarely mentioned in the literature. To further investigate the SLS, Particle Image Velocimetry (PIV) experiment and Computational Fluid Dynamics (CFD) simulation are adopted to obtain its swirl distribution characteristics, and the effect of port angles on scavenging performance is discussed based on engine fired cycle simulation. The results illustrate that Reynolds Stress Turbulence model is accurate enough for in-cylinder flow simulation. Tangential and axial velocity distribution of the flow, as well as the scavenging performance, are mainly determined by geometric scavenge port angles α_geom and β_geom. For reinforcement of scavenging on cross-sections and meridian planes, α_geom value of 27° and β_geom value of 60° are preferred, under which the scavenging efficiency reaches as high as 73.7%. Excessive swirl intensity has a negative effect on SLS performance, which should be controlled to a proper extent.     

Assessment and improvement of k-ω-γ model for separation-induced transition prediction

The purpose of this work is to improve the k-ω-γ transition model for separationinduced transition prediction. The fundamental cause of the excessively small separation bubble predicted by k-ω-γ model is scrutinized from the perspective of model construction. On the basis,three rectifications are conducted to improve the k-ω-γ model for separation-induced transition.Firstly, a damping function is established via comparing the molecular diffusion timescale with the rapid pressure-strain timescale...     

RBCC引射/亚燃模态过渡工作过程数值模拟

针对模态过渡过程中的瞬变特性,开展了引射/亚燃模态过渡工作过程中的RBCC进气道/燃烧室一体化数值模拟,计算比较了四种不同的模态过渡方案,分析了模态过渡工作过程中流动、燃烧模式和进气道状态的瞬变过程,研究结果表明:(1)推力产生较大波动的主要原因是燃烧室内燃烧/流动参数匹配性差和在流动方向上燃烧放热间断引起的;(2)主火箭保留的燃气在模态过渡过程中起到了火焰稳定和自持燃烧的作用,保留合理流量的燃气不仅可以缩短模态过渡时间,而且可以提高发动机的比冲;(3)提出了通过调节燃料喷注策略和主火箭节流方式实现模态平稳过渡的方案,并对该方案进行了数值验证,可望为进一步的实验研究提供了燃烧控制方法。     

转捩对平流层飞艇外形优化的影响

采用梯度算法优化了平流层飞艇外形以减少气动阻力,目标函数飞艇总阻力系数通过求解Navier-Stokes方程获得.优化过程中采用了SSTκ-ω转捩湍流模型和Realizable κ-ε全湍流模型求解绕流流场,采用大涡模拟方法详细分析了获得的两种飞艇外形的气动性能和特点,并且比较了两种湍流模拟方法的流场求解对梯度法优化的目标函数以及飞艇形状特点的影响.结果表明:是否考虑转捩对优化的结果有明显的影响,考虑转捩优化方案的气动性能明显优于全湍流流动求解优化方案.采用转捩模型的优化由于层流区域的存在改变了飞艇尾部的形状,从而获得了尾部分离区较小的飞艇外形轮廓.      

A numerical and experimental study of a transitional separation bubble

A combined numerical and experimental study of a two-dimensional transitional separation bubble due to an adverse pressure gradient is reported. The experiments have been performed in the MTL wind tunnel with a contoured wall imposing an adverse pressure gradient on the flow over a flat plate. The separated shear-layer is highly unstable and transition to turbulence occurs in the flow. The experimental separation bubble flow is modelled numerically using two-dimensional direct numerical simulations (DNS). Prescribing free stream boundary conditions in the wall normal velocity the experimental bubble is reproduced. The development of artificially forced two-dimensional instability waves is investigated and good agreement is found between experiments, simulations and linear stability theory (LST). The performance of several engineering transition prediction methods applied on the present separation bubble is presented and compared. Methods based on simplifications of the eⁿ-method yield predictions in accordance with experiments and DNS.     

Physics-based model for boundary layer transition prediction in a wide speed range

A new physics-based model employing three transport equations is developed for the simulation of boundary layer transitions in a wide speed range. The laminar kinetic energy is used to represent pretransitional streamwise velocity fluctuations, taking account of different instability modes. The fluctuation velocity components normal to the streamwise direction are modeled by another transport equation. Transition is triggered automatically with the development of the pretransitional velocity fluctuations. In the fully turbulent region, the model reverts to the k-ω turbulence model. Different test cases, including subsonic, supersonic and hypersonic flows around flat plates, airfoils and straight cones, are numerically simulated to validate the performance of the model. The results demonstrate the excellent predictive capabilities of the model in different paths of transition. The model can serve as a basis for the extension of additional transition mechanisms, such as rotation and curvature effects, roughness-induced transition and crossflow-induced transition.     

某型燃气轮机涡轮过渡段流场的数值模拟

对某型燃气轮机涡轮过渡流道流场进行数值模拟,得到了与实验结果比较吻合的计算结果。结果表明:过渡段内流动为亚声速流动,沿流道有较大的逆压梯度,整流支板的作用类似于扩压叶栅。该方法与结果为此类过渡流道的设计提供了参考,并为以后的进一步改进和优化打下了坚实的基础。     

Quantification of parametric uncertainty in γ-Reθ model for typical flat plate and airfoil transitional flows

Owing to the lack of physical knowledge of boundary layer transition, the γ-Re_θ transition model introduces closure parameters, which increase the uncertainty of transition prediction. The objective of this work is to quantify the uncertainties of closure parameters in the quantities of interests and identify the key parameters. The six closure parameters in the uncertainty intervals are used as input variables, and the uncertainties of the output results are propagated by a stochastic expansion based on the point-collocation nonintrusive polynomial chaos method. The relative contribution of each parameter to uncertainty is evaluated by the Sobol index. The computational cases include natural and bypass transitional flows on zero-pressure-gradient flat plates, and subsonic and transonic flows around airfoils. For most cases, c_e2, c_a2, and c_a1 dominate the uncertainty, and the influence of σ_θt is also significant when the history effects of flow are evident. The contribution of parameters in airfoils is more complex than that in flat plates. The transonic airfoil case shows that flow separation dramatically changes the distribution of Sobol indices, which poses a challenge to the accurate prediction of transition. Generally, c_e2 and c_a2 are the key parameters of the γ-Re_θ model.     

Extended Pile Driving Model to Predict the Penetration of the Insight/HP<Superscript>3</Superscript> Mole into the Martian Soil

The NASA InSight mission will provide an opportunity for soil investigations using the penetration data of the heat flow probe built by the German Aerospace Center DLR. The Heat flow and Physical Properties Probe (HP³) will penetrate 3 to 5 meter into the Martian subsurface to investigate the planetary heat flow. The measurement of the penetration rate during the insertion of the HP³ will be used to determine the physical properties of the soil at the landing site. For this purpose, numerical simulations of the penetration process were performed to get a better understanding of the soil properties influencing the penetration performance of HP³. A pile driving model has been developed considering all masses of the hammering mechanism of HP³. By cumulative application of individual stroke cycles it is now able to describe the penetration of the Mole into the Martian soil as a function of time, assuming that the soil parameters of the material through which it penetrates are known. We are using calibrated materials similar to those expected to be encountered by the InSight/HP³ Mole when it will be operated on the surface of Mars after the landing of the InSight spacecraft. We consider various possible scenarios, among them a more or less homogeneous material down to a depth of 3–5 m as well as a layered ground, consisting of layers with different soil parameters. Finally we describe some experimental tests performed with the latest prototype of the InSight Mole at DLR Bremen and compare the measured penetration performance in sand with our modeling results. Furthermore, results from a 3D DEM simulation are presented to get a better understanding of the soil response.     

壁面温度对斜爆震发动机进气道流动特性影响的数值研究

为研究壁面温度条件对层流、转捩、湍流状态下斜爆震发动机进气道流场结构、流场参数的影响,选取Ma10级、具有曲面压缩段的斜爆震发动机进气道为研究对象进行数值模拟,对进气道壁面附近激波诱导分离区、热边界层的变化进行了深入探讨。数值模拟结果表明,进气道肩部圆弧过渡段出现的再层流化现象,壁面冷却对其起抑制作用,绝热壁面条件下再层流化程度最为严重。壁面温度的增加有利于延缓流动转捩,同时也导致了分离区尺寸的增加以及转捩、湍流状态下分离区主体位置逐渐前移,进气道内通道的转捩为分离诱导转捩,转捩位置主要受到分离点位置的影响,整体表现为壁面温度增加转捩位置前移。进气道出口顶板侧热边界层厚度随着壁面温度的增加逐渐变厚,转捩状态下热边界层厚度变化可达5%,温度峰值也随着壁面温度的增加逐渐增加,且峰值位置逐渐靠近壁面。壁面温度条件相同时,层流状态下热流、热边界层厚度均较小。转捩、湍流状态下进气道出口顶板侧热边界层较厚,约为层流状态3倍,同时转捩、湍流状态下热边界层厚度相差可达2%。     

锯齿形转捩片触发高超声速进气道边界层转捩的大涡模拟

为促进锯齿形转捩片在高超声速进气道中的应用,以地面风洞条件下的二元进气道为研究对象,采用高精度大涡模拟方法对锯齿形转捩片在三级压缩楔面上触发的边界层转捩现象开展了研究。数值方法基于隐式亚格子模型,空间离散采用高精度通量限制型紧致格式,时间推进采用显式Runge-Kutta方法。数值模拟清晰捕捉到了边界层转捩的空间发展演化过程,并获得了统计平均流场以及流场脉动特征。数值模拟结果表明转捩片能够有效触发进气道压缩面边界层转捩;通过与等熵压缩面及单楔面数值模拟结果的对比分析,获得了转捩片触发边界层转捩的内在机理,为后续研究工作奠定了基础。     

基于放大因子与Spalart-Allmaras湍流模型的转捩预测

为了验证放大因子输运方程与Spalart-Allmaras(S-A)湍流模型耦合对转捩现象的模拟精度,选取Schubauer and Klebanoff(S-K)平板、S809低速翼型、30p30n多段翼型以及复杂的三维HiLiftPW-1构型进行自由转捩计算,并将计算结果与实验进行比较分析,其中针对S809算例,还与Langtry-Menter(L-M)转捩模型进行了比较.算例结果表明:放大因子输运方程与S-A湍流模型的耦合能够较好的捕捉转捩位置以及转捩发展过程,对分离泡诱导的转捩的模拟相比L-M转捩模型更精确,转捩位置的捕捉精度提升了10%;对比实验,多段翼转捩位置的捕捉误差最大为6.5%;针对三维高升力增升构型,以实验作为参考,全湍流计算与考虑边界层转捩的对比显示考虑边界层转捩能够更加精确的模拟气动力系数,升力和表面摩擦阻力系数的模拟精度精度提升1%.      

进气畸变对跨声速轴流压气机气动影响的机理研究

利用插板模拟跨声速轴流压气机进口流场畸变.通过实验和数值模拟两种手段详细分析了进口畸变下跨声速压气机的流动特点、激波沿周向的分布以及压气机的失速机理.实验和数值模拟的结果表明:插板使压气机前的不同周向位置产生正、负预旋流动;转子前速度的周向分布沿轴向是不断变化的,在距离转子较近时,畸变区的轴向速度反而较大;气流在过渡区A（转子进入低压区）时为正预旋,为非激波区;在过渡区B（转子退出低压区）和畸变区时为负预旋,为强激波区;过渡区的泄漏涡最强,是不稳定流动的诱发区域.      

超临界RP-3航空煤油喷嘴内部流动与相变特性研究

未来航空燃气轮机采用燃油作为冷却剂,对发动机冷却空气及热端部件进行冷却,使得燃油在进入燃烧室之前超过其热力学临界点,燃油的热力学性质发生巨大变化,将对燃油在喷嘴内的流动和喷嘴下游的喷射掺混过程产生重要的影响,因此有必要对超临界RP-3在喷嘴内的流动与相变特性展开研究。本文采用自主设计的喷嘴内部收缩通道模拟试验件,对超临界RP-3在喷嘴内的沿程压力分布进行测量,并采用基于一维等熵假设的计算方法进行数值模拟。试验首次获得了超临界航空煤油RP-3在喷嘴收缩通道内的沿程静压分布随喷射压力、温度变化的规律。通过与模型对比,发现未发生相变时,计算结果与试验值拟合精度较高,可以较好地预测RP-3在喷嘴内的流动参数,喷嘴内的静压分布也趋于一致;发生冷凝相变时,计算结果与试验值产生误差较大,喷射参数对于静压分布存在着较大的影响。本文试验获得的超临界RP-3航空煤油在喷嘴内的流动与相变特性,为航空发动机超临界RP-3喷射的喷嘴设计提供了重要的设计依据。