基于气体射流的气液两相流动减阻特性

为了明确射流表面减阻特性,基于气液两相流动原理,建立具有射流功能的气液两相流动减阻计算模型,采用数值计算方法,通过RNG k-ε湍流模型对射流表面不同射流方向角、以及不同气体射流速度情况下的流场进行数值模拟及实验验证,分析气体射流对壁面黏性阻力、减阻率及流场的影响,研究气体射流的气液两相的减阻特性。结果表明:气体射流表面能够减小壁面黏性阻力,具有较好的减阻效果;四种射流方向角下,射流表面减阻率均在气体射流速度为5m/s时达到极大值,并且当射流方向角在xoz平面为-30°时,减阻率最大,为7.18%;数值模拟的减阻率效果略高于实验结果;气体射流的气液两相流动对壁面边界层进行了有效控制,减小了壁面的黏性剪应力和雷诺应力,导致壁面的黏性阻力减小,使得射流表面具有减阻性能。     

飞行器大迎角下俯仰静、动导数的数值计算

 采用一种适用于非定常流动计算的逐次超松弛（SOR）时间推进方法耦合求解非定常Navier-Stokes方程和强迫运动方程,数值模拟了NACA0012翼型、弹道外形(HBS)和有翼导弹等飞行器在不同初始迎角下做强迫俯仰振荡的黏性动态流场,给出动态气动力及力矩系数的时间历程;采用积分法获得飞行器的静导数和动导数,将计算结果与实验结果及半经验理论预测方法进行比较并分析飞行器在大迎角下动态特性的非线性特征;计算结果与实验值相吻合,表明本文采用的非定常数值模拟方法在大迎角下能够较准确地模拟飞行器的动态特性,具有工程应用价值。     

弹体穿越平面激波三维非定常流场的数值模拟

用数值方法求解三维非定常可压缩Navier-Stokes方程，计算了复杂弹体在相同马赫数、不同迎角下的流场和气动力特征，计算结果与实验结果吻合，证明本文采用的方法能够准确计算此类复杂流场；然后采用该方法模拟了弹体穿越正激波的非定常过程，得到了弹体周围流场及气动力特征随穿越过程的非定常变化，阻力特征穿越后急剧下降。     

复杂外形气动力计算的阻力修正方法

基于气动力计算的位流线化理论,发展了一种亚音速附面层阻力修正方法。采用联合流场求解方法获得复杂外形的表面压力分布,在此基础上,将绕复杂外形物体沿流向面元划分为片条,在片条上,把三维附面层转换为二维附面层,对二维附面层进行数值求解,得到附面层粘性阻力修正。算例表明,该方法计算量小,与实验和N-S数值解结果吻合较好。为复杂气动外形飞行器初期设计、方案选型及评估提供了较为理想的工程设计手段。     

基于直接力控制的阵风响应及阵风减缓研究

考察某无人机模型的阵风气动响应特性和基于直接力操纵的阵风减缓效果.在求解非定常Euler方程时引入“网格速度”方法模拟阵风边界条件,通过动态嵌套网格技术实现舵面运动.首先对NACA0006翼型迎角阶跃型阵风的气动力响应进行了验证计算,计算结果与理论结果、文献计算结果吻合良好.在此基础上对某无人机模型在迎角阶跃型、One-minus-cosine型阵风作用下的气动力响应过程进行了数值模拟,并分析评估了不同舵面运动方式对飞机气动性能的影响.最终研究比较了不同舵面运动方式下阵风气动响应的减缓效果.结果表明:通过合理设计舵面运动,可以有效达到抑制阵风引起的非定常气动干扰的目的.     

高超声速黏性干扰相关性参数的修正

为了更精确地将黏性干扰相关性理论应用到实际工程上,基于高超声速平板黏性干扰相关性理论和斜激波理论,针对当前广泛采用的黏性干扰相关性参数(∨)∞进行了修正,建立了气动力关联性更好的参数(∨)∞.采用计算流体力学(CFD)数值方法模拟了航天飞机轨道器外形OV-102在典型状态下的黏性干扰效应,考察了修正后黏性干扰相关性参数(∨)∞对于气动力数据的关联特性.验证结果表明:修正后的黏性干扰相关性参数(∨)∞对于黏性诱导的各气动力数据均有良好的关联特性,并且相比未修正的黏性干扰相关性参数(∨)∞更加有效.     

飞机大迎角机动非线性气动力模型的辨识

给出一种基于付里叶分析的方法，用以分析飞机模型大迎角、大幅值强迫简谐振动获得的气动力和力矩数据。由一组不同频率飞机模型简谐振动的气动力响应建立升力、阻力、俯仰力矩系数的非线性模型。这些模型的最终表达形式包含阶跃类型函数的时间积分。利用飞机模型的试验数据验证本方法的计算结果。实际计算结果表明，本文给出的非线性气动力模型，能够较为精确地计算飞机模型大幅值迎角简谐振动和斜坡运动产生的气动力响应。     

基于CFD的机翼结冰过程分析

采用一种适用于分析机翼结冰过程的数值模拟方法,该方法在飞机结冰气象条件下采用计算流体力学（Computational Fluid Dynamics,CFD）对机翼流场进行计算,获得流场中二维机翼翼型的水滴撞击特性,并遵循质量和能量守恒方程,计算出机翼表面的结冰量,预测出机翼的积冰形状。所得冰形与文献中的试验数据吻合良好,表明数值模拟方法有效。另外还分析了结冰对机翼升力、阻力和压力系数的影响,其结果可以指导飞机机翼防／除冰系统的设计研究     

捆绑式运载火箭气动力计算

 应用细长体摄动理论和边界元数值解法发展了一种适用于计算捆绑式运载火箭气动力的工程数值计算方法。本法可利用已有芯级火箭的气动力数据,加上安装助推器后净增值,得到捆绑状态火箭的气动力,例如法向力、侧力以及力矩等。本法还可分别计算出捆绑状态下芯级和各助推器所受的气动力。对我国C火箭和日本N火箭的计算结果与实验数据比较,说明该方法能满足方案设计阶段精度要求。是方案设计过程中快速而经济的计算方法。     

临近空间高速飞行器微量气动力试验及计算

针对临近空间大气环境微量气动力风洞试验准确测量需求，研制了微量天平测力系统，实现了微量天平的结构设计和静态校准。采用钝锥简单外形进行试验验证，并对三角翼升力体复杂外形做了探索，结果表明：钝锥外形验证试验中，3次试验时天平各载荷单元的气动数据重复性精度均优于4.8%；探索试验中，三角翼升力体外形采用该风洞目前能达到最低密度状态，试验结果表明：该微量测力天平在极限状态下表现较好。在此基础上，利用数值计算方法对上述外形进行模拟，并将模拟结果与试验结果进行对比，表明数值计算得到的钝锥和三角翼升力体的气动力均与微量天平测力结果吻合较好；对于简单钝锥外形，在其试验条件下，钝锥表面压阻远高于摩阻；对于三角翼升力体外形，其试验条件下环境大气更加稀薄，三角翼表面摩阻占比与压阻相当。     

简化微小液滴在超声速气流横向喷射中的运动探索

从研究超声速气流中简化液滴(刚性小球)的气动力着手,比较分析使用CFD方法与使用两相流理论中已有的颗粒阻力系数模型得到的简化液滴气动力结果,得出Charles B.Henderson给出的阻力系数经验关系式适用于计算简化液滴在超声速气流中的运动。进一步,对不同直径简化液滴的运动开展工程方法的计算。在来流Ma=2.7的二维平板超声速流场中选取一个截面,作为气相流场,结果显示:(1)简化液滴与主气流存在相对超声速运动,当简化液滴直径dk≤0.12mm时,纵向相对超声速运动区域约为0.15m~0.4m,当dk0.12mm时,作用区域明显增大;(2)以10m/s喷射出的简化液滴,其横向穿透深度与纵向运动距离比约为0.004m/1m~0.021m/1m;(3)以100m/s的速度喷射出的30~120μm直径简化液滴,其横向穿透深度与纵向运动距离比约为0.02m/1m~0.055m/1m,实际过程中,小尺寸简化液滴的汽化很快,其穿透深度很小。     

降雨条件下翼型的气动特性

采用计算流体力学方法对降雨条件下翼型的气动特性进行了研究.通过求解定常可压的Navier-Stokes方程来计算流场,采用拉格朗日法对流场中的雨滴轨迹进行了模拟跟踪,得到了翼型表面的雨滴收集率;建立了水膜层数学模型,并假定水膜层沿着翼型表面的法向方向增长,预测了降雨引起的翼型外形的变化,并得到了水膜层表面的粗糙度高度;采用k-ωSST(shear-stress transport)两方程湍流模型结合增强的壁面函数,研究了降雨对翼型气动特性的影响.结果表明,在低迎角范围内,降雨对翼型的升力系数和阻力系数影响很小;当达到失速迎角后,降雨会引起上翼面边界层气流的提前分离,造成翼型气动性能的严重损失.      

低雷诺数翼型局部振动非定常气动特性

针对低雷诺数翼型特殊的气动特性,采用基于动网格的非定常数值模拟方法,研究翼型表面不同弦向位置的局部蒙皮以不同频率及振幅振动时对低雷诺数翼型气动特性及流场结构的影响,揭示蒙皮振动增升减阻的机理。研究表明,在低雷诺数条件下局部蒙皮振动可有效提高翼型气动特性,与刚性翼型相比蒙皮局部振动可使翼型升力系数提高,阻力系数降低,升阻比提高。振动位置对翼型气动特性及流场结构有显著的影响,振动表面位于翼型前缘附近或位于层流分离泡中心时可有效控制翼型层流分离,从而提高翼型气动特性。振动频率对翼型表面层流分离及转捩位置均有显著的影响,随着振动频率增加,翼型气动特性出现最优值。与刚性翼型相比,表面振动使翼型转捩位置略向上游移动,摩擦阻力增加,但振动使等效翼型相对厚度减小,压差阻力明显减小。在小幅振动范围内,随着振幅增加,流场非定常特性更加显著,翼型升阻比增加。     

A Structured Mesh Euler and Interactive Boundary Layer Method for Wing/Body Configurations

To compute transonic flows over a complex 3D aircraft configuration, a viscous/inviscid interaction method is developed by coupling an integral boundary-layer solver with an Eluer solver in a ＂semi-inverse＂ manner. For the turbulent boundary-layer, an integral method using Green＇s lag equation is coupled with the outer inviscid flow. A blowing velocity approach is used to simulate the displacement effects of the boundary layer. To predict the aerodynamic drag, it is developed a numerical technique called far-field method that is based on the momentum theorem, in which the total drag is divided into three component drags, i.e. viscous, induced and wave-formed. Consequently, it can provide more physical insight into the drag sources than the often-used surface integral technique. The drag decomposition can be achieved with help of the second law of thermodynamics, which implies that entropy increases and total pressure decreases only across shock wave along a streamline of an inviscid non-isentropic flow. This method has been applied to the DLR-F4 wing/body configuration showing results in good agreement with the wind tunnel data.     

水滴动力学行为对其撞击特性影响

设计了一种压电式等径离散水滴发生器,利用高速摄影的手段对水滴变形现象进行了实验研究。建立了水滴变形后的非球形阻力模型,同时引入了水滴破碎和飞溅模型,采用欧拉法对不同条件下水滴撞击NACA23012翼型的过程进行了数值模拟,并与传统模型下水滴撞击特性计算进行了对比。研究表明:水滴变形会导致其阻力系数增大,水滴阻力系数曲线在雷诺数大于500时逐渐脱离球形阻力曲线,不断增大接近圆盘阻力系数;水滴阻力特性的改变会影响水滴运动轨迹,增强水滴的气流跟随性;水滴飞溅会造成水滴质量损失,减小翼型表面局部水收集系数,这与实验数据更加吻合。      

A new non-linear vortex lattice method:Applications to wing aerodynamic optimizations

This paper presents a new non-linear formulation of the classical Vortex Lattice Method (VLM) approach for calculating the aerodynamic properties of lifting surfaces. The method accounts for the effects of viscosity, and due to its low computational cost, it represents a very good tool to perform rapid and accurate wing design and optimization procedures. The mathematical model is constructed by using two-dimensional viscous analyses of the wing span-wise sections, according to strip theory, and then coupling the strip viscous forces with the forces generated by the vortex rings distributed on the wing camber surface, calculated with a fully three-dimensional vortex lifting law. The numerical results obtained with the proposed method are validated with experimental data and show good agreement in predicting both the lift and pitching moment, as well as in predicting the wing drag. The method is applied to modifying the wing of an Unmanned Aerial System to increase its aerodynamic efficiency and to calculate the drag reductions obtained by an upper surface morphing technique for an adaptable regional aircraft wing.     

基于阻力气动特性计算的飞艇艇身外形研究

采用SIMPLE算法、S-A湍流模型求解了三维不可压雷诺平均N-S方程,数值模拟了绕飞艇的低速不可压粘性流动,计算了不同艇身外形的阻力特性.本文方法计算结果与实验结果或其他文献的计算结果符合良好,可以用于飞艇的气动特性计算.通过对6个不同外形的相同体积飞艇阻力特性计算,得到了最佳艇身外形,表明本文方法可以用于艇身外形的选型设计.     

俯仰振荡翼型推阻力转变滞后机制数值研究

扑翼产生的反卡门涡街被认为是一种推力型尾迹,但已有研究指出,随着斯特劳哈尔数(St)增大,低雷诺数下俯仰振荡翼型的净推力产生明显滞后于反卡门涡街的出现。为探究该现象背后的物理机制,对NACA0012翼型在雷诺数1000条件下作简谐俯仰运动的流场进行了数值模拟。采用翼型表面积分方法和基于有限控制体的气动力估计方法分别研究了翼面分布力特性和尾迹流场特性变化对阻力-推力转变临界点的影响。翼面分布力积分结果表明,当翼型振荡参数进入对应反卡门涡街尾迹形态区域时,翼面压力分布产生的推力分量无法克服剪切层的黏滞阻力,是造成俯仰翼型推力产生滞后于反卡门涡街出现的主要原因。对尾迹流场及相应的推阻力特性变化的分析表明,尽管反卡门涡街会产生“喷流效应”,但在St较小时,其产生的动量诱导推力无法克服反卡门涡自身的涡致阻力和尾迹流场压力损失引入的附加阻力,因此即使存在反卡门涡街也不能产生净推力,进而从流场分析的角度进一步解释了这一滞后现象的发生机制。