基于修正面元法的机翼焦点位置变化研究

针对静气动弹性问题中机翼焦点位置求解精度与效率的矛盾,提出采用分段斜率修正面元法来计算气动力,该方法以刚体机翼在不同迎角下的CFD数据作为修正基础,对弹性变形后的机翼下洗采用分段修正,克服了传统面元法无法考虑非线性因素的缺点。通过M6机翼气动焦点位置的计算表明:相对于CFD方法,该方法保持了较高的准确度并且计算效率较高;相对于传统面元法,该方法能够准确地计算焦点位置和弹性机翼焦点位置变化量,适用于工程设计。     

叶轮机计算流体动力学技术现状与发展趋势

 通过对 863高负荷风扇、2 0 0MW汽轮机组改进、组合压气机改进和轴流大小叶片压气机的研制,说明了全三维定常CFD技术在我国叶轮机研究中的作用,并指出了当前的主要问题;通过对风扇 -吊架干扰,单级压气机和处理机匣的研究,说明了非定常CFD技术对叶轮机气动力学研究的作用和价值。结合上述全三维定常和非定常CFD应用的实例,论述了我国叶轮机CFD技术面临的挑战,及其发展趋势。     

改进的综合CFD与畸变合成法预测动态畸变(英文)

对综合CFD与畸变合成方法进行了改进。通过CFD求解,获得AIP上各稳态气动参数。采用BP人工神经网络直接求解高度非线性的湍流相关方程,该方程将CFD求解的稳态气动参数与试验测得湍流度进行关联。经过训练的网络可根据稳态CFD结果预测湍流度。之后根据CFD计算的稳态总压和神经网络预测的湍流度合成动态压力。最后根据合成的动态压力寻找最大瞬时畸变。文中使用了六个工况的飞行数据验证该方法的准确度并且与原方法的计算结果进行了比较。结果表明该方法预测的畸变值与试验畸变值吻合较好,且预测结果较原方法与试验值更接近,该方法能够在试验前用来预测动态畸变,具有很强的工程价值。     

基于降阶模型和梯度优化的流场加速收敛方法

CFD是解决流体相关复杂工程问题的重要手段，随着流体力学研究和工程应用问题的日益复杂，更高效、精确的CFD算法研究至关重要。提出了一种基于降阶模型（ROM）和梯度优化用于加速稳态流场的收敛方法。该方法首先收集CFD迭代求解过程中的流场快照构造基于正交分解（POD）的降阶模型，然后使用梯度优化方法求解降阶模型得到残差更低的流场，最后以该流场为初值继续进行CFD迭代求解从而加速CFD的收敛过程。将该方法用于基于有限体积法的CFD求解器中，针对无黏、湍流算例验证了所提加速收敛方法的有效性。结果表明，加速收敛方法相比于初始CFD方法迭代步数能显著减少2～3倍，是一种高效且鲁棒的加速收敛方法，且在复杂工程问题中具有广泛的应用前景。     

计算流体力学在航空叶轮机设计中的应用

计算流体力学(CFD)在航空发动机叶轮机设计和分析中 的重要性日益明显。本文介绍了在叶轮机设计中CFD的意义、叶轮机 CFD发展概况及CFD发展中的几个问题。     

基于Richardson插值法的CFD验证和确认方法的研究

针对典型高超声速流动问题(T2-97模型),利用Richardson插值法,研究了多套连续变化网格下数值解的空间离散误差、收敛性。通过网格收敛性研究,完成了CFD的验证过程。利用不确定度分析方法,结合实验数据,开展了该问题的确认。研究表明:针对所选问题,在现有实验数据及计算条件下,该问题的CFD得到了验证与确认。     

CFD计算网格误差分析的一个算例

为了给一般工程问题的CFD仿真的网格划分提供参考,研究了光滑直圆管中粘性稳态等温流动在不同尺寸、不同纵横比网格情况下的CFD计算数值解的误差.采用FLUENT软件进行数值计算,应用CFD仿真验证和确认的基本方法,与精确解比较,计算了各种网格的数值计算结果的误差,得到了网格独立解的网格尺寸,给出了误差随网格尺寸变化的规律.     

CFD在我国城市风热环境中的研究与应用

城市风热环境研究对于城市可持续发展、居民健康和气候变化适应性都具有重要意义。随着高性能计算能力的不断提升和数值模拟方法的改进,计算流体力学（Computational Fluid Dynamics, CFD）已经广泛应用于城市风热环境预测。本文概括性地总结了城市风热环境实验研究和CFD研究的进展,包括典型的实验研究方法、CFD数值模型、模型验证等,并重点从建筑形态布局、植被、水体以及冷性材料4个方面对CFD应用研究成果进行归纳总结,最后从多季节、多尺度、精细化模拟以及多因素耦合等角度对未来城市风热环境CFD模拟研究进行了展望,并提出了参考建议。     

基于CFD/CSD的非线性气动弹性分析方法

提出了一种基于计算流体力学/计算结构动力学(CFD/CSD)的非线性气动弹性分析方法,并应用于切尖三角翼的非线性颤振和极限环振荡(LCO)研究。该方法将非线性有限元(FEM)和CFD计算相结合,并辅以高精度的界面插值,能够分析结构和气动非线性共存的气动弹性问题。结构部分以四边形平板壳元为基础,采用更新的拉格朗日(UL)方法分析结构大变形引起的几何非线性问题。气动部分以Navier-Stokes(N-S)方程作为控制方程,采用CFD方法计算跨声速气动力。机翼的非线性颤振计算表明了方法的有效性。最后应用该方法研究了切尖三角翼的LCO现象,其计算精度明显优于已有结果。     

一种CFD/CSD耦合计算方法

针对柔性大展弦比机翼气动弹性分析和主动弹性机翼(AAW)设计发展了一种计算流体动力学(CFD)和计算结构动力学(CSD)的耦合计算方法。其主要思想是采用在同一物理时间弱耦合求解CFD/CSD技术。气动力采用非定常N-S方程的双时间有限体积求解技术,结构响应则采用有限元数值求解技术。CFD和CSD耦合计算的边界信息(气动力和网格)由所设计的界面程序传输。网格信息传输采用守恒体积转换(CVT)方法将CSD计算结构响应位移插值到CFD网格点上。变形已有的CFD网格技术用以确定CFD的变形网格。以位移或载荷的迭代误差为判断耦合计算的收敛标准。最后得到了机翼在Ma=0.8395,α=5.06°时CFD/CSD耦合计算的收敛值。针对计算结果分析了机翼受静气动弹性过程中结构响应和气动特性随时间变化的效应。初步研究结果表明:这种弱耦合方法求解非线性气动弹性问题是可行的。     

基于广义Richardson外插方法的阻力预测精度分析

利用计算流体力学(CFD)方法预测阻力是飞行器气动设计中的关键环节。采用广义Richardson外插方法分别对数值方法预测二维简单构型的压差阻力、摩擦阻力和三维复杂构型的总阻力的精度进行了分析。在NACA 0012翼型无黏绕流和平板湍流边界层两个算例中验证了NSAWET程序和广义Richardson外插方法,分别得到数值算法预测压差阻力和摩擦阻力能达到的名义精度。进而模拟三维通用研究模型(CRM)翼身组合体绕流,得到的阻力名义精确值在DPW 5的统计误差带范围之内;综合DPW 5的计算结果来看,不同CFD解算器的结果之间存在一定差别,阻力预测精度总体上不符合二阶。可见,标准Richardson方法采用的二阶精度假设难以普遍适用,有必要采用广义Richardson外插方法得到名义精度。针对不合理的名义精度,采用Roache建议的方法加以限制。广义Richardson外插方法有助于提高误差分析的合理性,可以进一步降低网格对阻力预测的影响。     

基于非定常气动力辨识技术的气动弹性数值模拟

选择离散型输入输出差分模型,运用最小二乘方法进行非定常气动力建模,并将辨识得到的降阶模型用于气动弹性的数值模拟。1个马赫数下的颤振临界点的计算仅需调用一次非定常流场求解器。计算精度保持与非定常欧拉方程计算方法相当的同时计算效率提高了1~2个量级。计算了跨声速具有S型颤振边界的气动弹性标准算例－Isogaiwing和三维气动弹性标模算例AGARD445.6,辨识模型计算边界与非定常Euler方程计算结果吻合。证明非定常气动力辨识技术可以提供高效的高精度的气动弹性分析。     

AUSM+方法在 粘性内流模拟中的应用

本文阐述了三维RANS方程求解程序中使用的AUSM+离散方法以及Spalart-Allmaras湍流模型，然后对两个算例—Sajben收缩扩散通道内跨声速流动和NASA TP-3252文献中给出的透平静叶排内部流动，进行了数值模拟，并把数值结果与实验结果进行了比较，说明三维RANS方程求解程序可以准确模拟复杂的内部流场。     

基于CFD/CSD耦合方法的旋翼气动弹性载荷计算分析

为提高直升机前飞状态下旋翼非定常气动弹性载荷的预估精度,在旋翼气动弹性综合分析方法中引入旋翼CFD模块,建立了一套基于CFD/CSD松耦合分析的计算方法和程序。为高效解决流固耦合方法中由于桨叶挥舞、扭转等弹性变形带来的旋翼贴体网格变形问题,采用基于代数变换方法的网格变形技术,桨叶运动变形量和旋翼气动力信息通过流固交接面传递。旋翼流场分析方法中,主控方程采用耦合S-A湍流模型的Navier-Stokes方程,围绕旋翼流场的网格采用结构嵌套网格方法生成,无黏通量计算采用Roe格式,时间推进采用双时间法。旋翼结构分析中,考虑旋翼配平,基于Hamilton变分原理和20自由度Timoshenko梁模型求解弹性旋翼非线性运动方程。分别对CSD和CFD方法进行验证,在此基础上,计算了SA349/2旋翼桨叶在前飞状态下的非定常气动力、挥舞弯矩、摆振弯矩和扭转力矩,并与飞行测试数据进行了对比。计算表明:CFD/CSD耦合方法可以显著提高旋翼非定常气动弹性载荷的分析精度,精确捕捉桨叶表面压强峰值、激波位置等,表明本文发展的旋翼CFD/CSD耦合方法可以有效地运用到旋翼气动弹性载荷的预测分析中。     

基于CFD的圆形静压油垫数值仿真技术研究

应用CFD技术,采用二维轴对称流动模型对圆形静压油垫内流场特性及其承载能力和油膜刚度进行了定量分析与计算,并与静压设计公式得到的计算结果进行了比对,二者计算结果吻合.研究表明,采用CFD技术能够实现对静压油垫的设计计算.提高了设计效率和计算准确性.     

CFD中统计误差的数值精度分析

在计算流体力学（CFD）的算法研究中经常会对离散误差进行数值精度分析，通常以统计误差的各范数为研究对象，最常用的统计误差范数为L₁范数、L₂范数和L_∞范数，一般认为各范数在数值精度上具有等价性。实际上，由于流场局部存在间断、网格局部不连续或者是在极值点附近采用非线性加权插值等可能使数值方法存在局部降阶问题，导致统计误差各范数所表达的数值精度并不一致。通过详细的理论分析，揭示了统计误差各范数所表达的数值精度之间的关系，并通过相应的数值试验予以验证。研究结果不仅能够指导CFD算法的数值精度验证工作，而且也可为更为复杂流动模拟的数值精度判定提供理论依据。     

壁判据用于计算流体力学（CFD）可信度评估

本文把作者提出的近壁干扰剪切流动(ISF)全域理论与流体运动方程组及流体在壁面上无滑移条件相结合导出一组壁面判据.壁判据为计算流体力学(CFD)仿真的可信度评估提供了基于流体理论的一条直接验证途径.对不可压缩流动的十一个熟知的NS方程组精确解,包括二维驻点和斜入射再附点流,二维分离点和背风驻点流,轴对称驻点和背风驻点流,旋转圆盘附近的三维Von Karman 流、收缩和扩张渠道流和非定常斜入射三维驻点流;以及经典边界层及其无粘外流和相似性边界层及其粘外流的NS方程组解,提出可用于验证近壁流动计算的几个壁相关函数.证实它们准确满足所有壁判据;说明壁判据可用来检验NS方程组数值解近壁计算结果的计算精度并验证其可信度.     

基于高精度模态气动力的跨声速静弹高效分析方法

跨声速静弹分析一直是工程设计中的难点问题,以模态坐标系下的线性静弹方程为基础,提出了基于高精度模态气动力的跨声速静弹高效分析方法,该方法仍需求解线性静弹方程,但对于其中关键的模态变形引起的弹性气动力增量,采用由结构变形到气动力的单向计算流体动力学(Computational Fluid Dynamics,CFD)/计算结构动力学(Computational Structural Dynamics,CSD)耦合方法获得,实现了高效线性方法与高精度CFD/CSD耦合方法的有效融合。以某小展弦比机翼基本状态、舵偏状态以及某型战斗机跨声速副翼效率的静弹分析为例,对比分析了本文方法、经典线性方法、CFD/CSD耦合方法的计算结果以及某型机的试飞辨识结果。分析结果表明,所提出的方法在计算效率、精度和鲁棒性方面具备综合优势,具有较高的工程应用价值。     

Prediction of high lift: review of present CFD capability

A survey is conducted of CFD methods applied to the computation of high-lift multi-element configurations over the last 10–15 years. Both 2-D and 3-D configurations are covered. The review is organized by configuration, in an effort to glean useful insights with respect to particular successes or failings of CFD methods as a whole. In general, for both 2-D and 3-D flows, if certain guidelines regarding grid, transition, and turbulence model are followed, then surface pressures, skin friction, lift, and drag can be predicted with reasonably good accuracy at angles of attack below stall. Velocity profiles can generally be predicted in 2-D flow fields, with the exception of the slat wake, which tends to be predicted too deep by most CFD codes for a range of different configurations. CFD codes can usually predict trends due to Reynolds number in 2-D, but they are inconsistent in the prediction of trends due to configuration changes. On the whole, 2-D CFD is unreliable for predicting stall (maximum lift and the angle of attack at which it occurs); in most cases, maximum lift is overpredicted, but for some configurations the opposite occurs. However, there is some evidence that stall misprediction of nominally 2-D experiments may be caused by 3-D effects, which are obviously not modeled by 2-D CFD. In general, 3-D computations are also inconsistent with respect to computing stall, but there have been fewer of these applications to date. The paper concludes with a list of challenges that confront CFD at the start of the next decade, which should witness a dramatic increase in the number of CFD applications for 3-D high-lift configurations.