基于CFD技术的螺旋桨气动特性研究

通过求解绝对坐标系下的雷诺平均Navier-Stokes方程数值模拟螺旋桨侧流黏性绕流,计算螺旋桨非定常气动特性.为了保证时间精度和提高计算效率,采用含牛顿型LU-SGS(lower-upper symmetric-Gauss-Seidel)子迭代的全隐式双时间法并引入多重网格技术.为了描述桨叶的相对运动,采用运动嵌套网格技术.对于轴流状态,计算结果与实验值吻合良好;对于侧流状态,与工程算法1P(once-per-revolution)载荷法进行了比较,验证了方法的有效性.应用该方法计算螺旋桨在不同来流偏角下的周期性气动特性,可有效地用于螺旋桨的设计.      

平板自由下落的数值模拟

采用动网格技术,通过耦合求解N-S方程和运动方程的方法进行数值模拟,研究二维平板的自由下落运动。选取了不同质量分布的二维平板作为研究对象,研究初始角度、转动惯量和质心位置等因素对其自由下落运动轨迹的影响,以及其中的流体力学和动力学机理。数值模拟结果显示了二维平板的自由下落运动,除了摆动运动、翻滚运动和定常下落等经典运动轨迹外,随着质心位置的偏移还会衍生出非对称摆动、摆翻运动等特殊运动轨迹。研究结果表明:无量纲转动惯量和质心位置是决定二维平板自由下落运动方式的重要参数,但两者对下落运动影响的物理机制存在差异。这些结果对相关物理问题,特别是低雷诺数流体力学和动力学研究,具有一定的指导意义。     

用粒子仿真法计算含缺陷固体发动机内流场

利用“N-S方程 粒子仿真”方法来计算含缺陷固体发动机的内流场，计算结果表明：该模型可以很好模拟裂纹对发动机内流场的影响，同时指出环形裂纹内的流动对主流的影响大于纵向裂纹，环形裂纹内上游位置流动对主流的影响较下游位置的大，并可以保持到喷管喉部，翼槽附近的环形裂纹内的流动，对主流的影响相对较弱。     

先进直升机旋翼悬停状态气动性能计算

 为更好地模拟先进直升机旋翼流场,同时准确计算其悬停状态的气动性能,建立了一套基于Navier-Stokes/Euler方程的混合CFD方法。该方法的求解域由两部分组成：一是围绕旋翼桨叶周围的黏性区域,采用可压缩Navier-Stokes方程来模拟旋翼附近的黏性流动和近场尾涡的捕捉;二是离桨叶较远、黏性可以忽略的远场区域,用Euler方程来描述其流动。在该方法中,将三阶逆风格式（MUSCL）与通量差分分裂方法相结合,无需添加人工黏性,因而可有效地减少旋翼尾迹数值耗散。?阌诹鞒》智蠼庖约爸芷谛员呓缣跫氖凳?采用了嵌套网格方法,并给出旋翼网格与背景网格交界面的信息传递方式。应用所建立的计算方法,首先对二维翼型、三维M6机翼的流场进行了数值模拟,以验证计算方法;然后,着重计算了有实验结果可供对比的具有先进气动外形的HELISHAPE 7A模型旋翼和UH-60A直升机旋翼,通过计算得到旋翼表面压力分布、桨叶展向拉力系数分布、桨叶表面细节流动以及气动性能等,进一步验证了该方法的有效性。     

大迎角下翼型振动的气动性能研究

众所周知,风洞试验中的飞机模型,尤其是带有大展弦比机翼的模型有时会出现翼梢振动现象,振动模式主要表现为翼梢沉浮和俯仰形式,以致影响实验结果的精度和可靠性.选取相对厚度较小的NACA0008翼型,在求解非定常Navier-Stokes方程的基础上,采用改进的无限插值理论和绕翼型的C型网格,模拟风洞实验中模型振动条件下的流场,研究振动模式及其不同耦合对流场、尤其是大迎角流场的影响,并考虑了模型弹性轴不同位置对结果的影响.研究结果表明:在临近传统定常失速迎角的大迎角条件下,翼型的振动可以引起翼型大尺度的分离,导致翼型失速的提前发生;振动在不同的相位滞后条件下,对翼型流场的分离程度不尽相同.     

机体/动力装置一体化分析中的动力影响效应数值模拟

采用合适的进、排气边界条件, 对带动力的单独涡扇发动机模型和带动力的DLR-F6翼身架舱复杂组合体(WBPN)模型的流场进行了数值模拟和分析.运用Roe三阶迎风偏置通量差分裂方法和隐式近似因子分解方法求解雷诺平均Navier-Stokes方程, 基于点搭接多块网格技术生成高质量计算网格, 通过多重网格技术来加速收敛.着重对动力发动机进、排气边界条件的给定进行了推导, 以准确地模拟发动机动力效应.计算结果与实验值吻合良好, 验证了所采用计算方法、网格分块策略和边界条件处理方法的正确性.      

超燃冲压发动机总体化性能分析

针对由后掠侧压进气道、矩形燃烧室和后掠侧板尾喷管构成的超燃冲压发动机模型,为了快速准确得到超燃冲压发动机的整机性能,建立了一套分析发动机总体性能的耦合流场计算方法。其中进气道和尾喷管流场分别求解三维单一气体和多组分抛物化N-S方程,燃烧室流场求解一维N-S方程,各部件在交接面处传递几何和流动参数。运用总体性能计算方法,研究了发动机构型的气动特性,给出了飞行高度、马赫数和攻角改变对发动机性能的影响规律。结果表明:该总体性能计算方法可用于超燃冲压发动机总体性能的快速预估;飞行高度对推力、升力和俯仰力矩的影响均较大且程度相近,三者随飞行高度增加均减小;飞行马赫数对发动机推力、升力的影响相近,对俯仰力矩的影响最大,三者随飞行马赫数增加均增加;飞行攻角对升力的影响最大,俯仰力矩次之,对推力的影响最小,随飞行攻角增加,发动机升力近似线性增加,推力和俯仰力矩分别呈增加和减小趋势。     

Robust Navier-Stokes method for predicting unsteady flowfield and aerodynamic characteristics of helicopter rotor

A robust unsteady rotor flowfield solver CLORNS code is established to predict the complex unsteady aerodynamic characteristics of rotor flowfield. In order to handle the difficult problem about grid generation around rotor with complex aerodynamic shape in this CFD code,a parameterized grid generated method is established, and the moving-embedded grids are constructed by several proposed universal methods. In this work, the unsteady Reynolds-Averaged Navier-Stokes(RANS) equations with Spalart-Allmaras are selected as the governing equations to predict the unsteady flowfield of helicopter rotor. The discretization of convective fluxes is accomplished by employing the second-order central difference scheme, third-order MUSCL-Roe scheme, and fifth-order WENO-Roe scheme. Aimed at simulating the unsteady aerodynamic characteristics of helicopter rotor, the dual-time scheme with implicit LU-SGS scheme is employed to accomplish the temporal discretization. In order to improve the computational efficiency of holecells and donor elements searching of the moving-embedded grid technology, the ‘‘disturbance diffraction method and ‘‘minimum distance scheme of donor elements method are established in this work. To improve the computational efficiency, Message Passing Interface(MPI) parallel method based on subdivision of grid, local preconditioning method and Full Approximation Storage(FAS) multi-grid method are combined in this code. By comparison of the numerical results simulated by CLORNS code with test data, it is illustrated that the present code could simulate the aerodynamic loads and aerodynamic noise characteristics of helicopter rotor accurately.     

基于非惯性系的悬停状态旋翼CFD/CSD耦合气动分析

旨在提高先进旋翼气动特性的分析精度,在旋翼高精度CFD分析中耦合气动弹性效应,取代传统方法中的刚性桨叶假设,并考虑悬停状态旋翼流场准定常的特性,在非惯性坐标系下建立了一套适合于悬停状态旋翼气动特性计算的CFD/CSD耦合分析方法。旋翼气动载荷通过求解三维Navier-Stokes方程求得,空间离散及通量计算采用Jameson中心格式,时间方向则选用五步Runge-Kutta迭代求解,湍流模型采用B-L模型;基于Hamilton原理建立了描述旋翼弹性运动的非线性微分方程,针对旋翼悬停状态的工作特点,采用Raphson迭代方法求解获得旋翼桨叶的弹性变形量。在CFD/CSD耦合计算中,旋翼桨叶交接面载荷及变形信息通过CFD与CSD模块进行传递,同时为提高桨叶弹性变形后贴体网格生成的效率和质量,采用基于网格点坐标转换的网格变形方法。在CFD和CSD程序分别验证基础上,采用建立的旋翼CFD/CSD耦合分析方法计算了先进的UH-60A直升机旋翼的表面压强及气动载荷。计算结果表明,与刚性旋翼CFD模拟结果比较,本文建立的CFD/CSD耦合分析模型可以更准确地预估旋翼气动载荷和性能。     

高超声速内收缩进气道分步优化设计方法

提出了基准流场与唇口平面形状分步优化的高超声速内收缩进气道设计方法。基准流场以反射激波不均匀性最小和总压恢复最大进行多目标优化设计,使用结合Tayler-Maccoll方程的有旋特征线方法(MOC)进行流场计算,获得双拐点母线内收缩锥基准流场。进气道唇口形状以沿流线积分(Streamline Integral Method, SIM)获得的进气道无黏阻力最小为目标进行优化设计,获得类椭圆形唇口平面形状。针对优化设计结果进行数值模拟,与传统直母线基准流场相比,双拐点母线基准流场反射激波后流动不均匀性下降40%左右,总压损失减少35%左右,总体性能提升明显。类椭圆唇口进气道在设计点的单位质量流量无黏阻力相较于圆形唇口降低6%,具有良好的压缩特性和气动效率,能够减弱进气系统对飞行器气动性能的不利影响。研究结果表明该方法是一种高效且实用的高超声速内收缩进气道设计方法。