一种CFD/CSD耦合计算方法

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

CFD商用软件在民用支线飞机高速计算中的应用

本文主要是采用用CATIA V5.10版本制作的民用支线飞机计算数模.用ICEM-CFD网格软件生成各种相应的结构化全六面体计算网格.数值计算的求解器采用FLUENT6.1.22、CFX5.6两种商用软件.对民用支线飞机的纵、横向高速飞行状况(包括飞机各操纵面的不同偏转角)进行了较为全面的计算.计算结果良好.     

采用TRIP2.0软件计算DLR-F6构型的阻力

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了DLR-F6构型,主要目的是通过计算DLR-F6构型的安装阻力考察TRIP2.0软件的数值模拟精度,并为运输机构型的气动特性计算积累经验.本文数值模拟采用的多块对接网格,测压和测力的试验结果均来自AIAA CFD Drag Prediction Workshop II(DPWII),对比计算结果采用了CFL3D的结果.本文详细研究了网格密度、湍流模型对DLR-F6翼身组合体和翼/身/架/舱复杂组合体两种构型的的总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致.本文采用SST两方程模型计算两种构型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布影响较小.     

基于粘性笛卡尔网格的N-S方程计算

本文采用基于叉树数据结构的粘性笛卡尔网格求解N-S方程。为了使网格的分布更为合理,采用基于物体几何外形的方法对网格进行自适应加密。在物面处采用光顺、投影、分割加密的方法保证了物面处网格的正交性与贴体性,从而生成高质量粘性笛卡尔网格。将有限体积法应用于粘性笛卡尔网格求解N-S方程,采用了SA紊流模型。算例计算结果与实验数据符合良好。     

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

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

基于结构动网格的无人机地面效应研究

为了更好地分析地面效应对无人机气动特性的影响,采用一种新的结构动网格生成技术,运用CFD手段对某无人机的地面效应进行了研究。在一套静态网格CFD计算基础上,运用结构动网格生成技术,获得非定常的网格移动,得到一系列的不同攻角、地面效应高度的CFD网格及初始流场。在此基础上进行定常流场求解,计算得到地面效应的定常CFD结果。最后对某双尾撑布局无人机的地面效应进行了CFD计算研究,对相应的全机气动特性的影响进行了分析。计算结果表明：地面效应对平尾的影响大于对机翼的影响并使平尾的失速迎角提前;该方法是可行的,具有一定的工程应用价值。     

一种适合于旋翼涡流场计算的非结构自适应嵌套网格方法

针对旋翼涡流场的特点,提出了一种新的非结构自适应运动嵌套网格方法。在该运动嵌套网格方法中,贴体网格采用非结构网格,背景网格使用可自适应的笛卡尔网格,以便发挥非结构网格的形状描述能力和笛卡尔网格的空间自适应优势,方便网格生成过程并提高计算精度,同时克服在物体表面附近自适应笛卡尔网格生成工作量大和过程烦琐的不足。在此基础上,使用求解N-S主控方程的方法,分别对旋翼的悬停和前飞流场进行了数值模拟,计算结果表明本文的非结构自适应运动嵌套网格方法对提高旋翼涡流场的细节捕捉能力和计算精度是很有效的。     

多面体网格在静气动弹性计算中的应用

介绍了多面体网格技术在静气动弹性问题中的应用.基于M6算例和径向基函数插值法,探究了多面体网格的CFD计算效率和变形能力,使用AGARD445.6机翼进行了静气动弹性数值计算.结果表明,多面体网格具有计算效率高、变形能力强、计算结果准确等优点,可适用于静气动弹性问题的数值计算.     

基于混合网格的DLR-F6模型跨声速阻力预测粘流计算

该文详细介绍了针对第二届AIAA阻力预测会议公布的标准模型及其网格进行的基于混合网格的纳维尔-斯托克斯粘流计算考题,我们将自主开发的WoF90程序的计算结果与DLR的TAU程序以及FOI的EDGE程序的计算结果进行了对比分析,发现三种程序的升力系数普遍比实验结果偏高,但偏高的程度大致相同,阻力计算的精度范围大体一致。结果确认了300万网格点的计算结果尚不足以得到网格收敛的效果,阻力预测需要更大规模的网格。     

一个CFD计算作业调度系统的设计与实现

针对目前CFD行业中计算作业的特点、用户计算资源状况和使用特点,以及对现有相关作业调度系统的分析研究,设计出采用共享目录方式进行远程计算的作业调度系统.该系统不仅具有对局域网内计算资源优化调度的基本功能,而且具有对计算作业的快速、安全、实时监控运行以及不在远程计算机上产生垃圾文件等特点,更好地解决了CFD行业局域网中作业优化调度的问题.本文阐述了该系统的总体架构,系统功能,以及关键技术,并对系统的可用性进行了分析.     

旋转坐标系下分区计算的LU隐式方法

将隐式时间迭代方法应用于并行计算是CFD研究中的热点问题,研究发现LU-SGS格式在旋转坐标系下对分区边界比较敏感,对旋翼悬停状态进行分区并行计算时碰到了计算发散的问题。针对此问题,对基于LU分解的3种隐式时间迭代格式（LU-SGS、DP-LUR和HLU-SGS）进行了对比研究,设计了静止流场下的旋转网格算例对计算方法进行测试。结果表明,LU-SGS格式在网格边界采用简化处理方法,当边界的逆变速度增大时会引起数值误差放大,误差的不断积累导致计算发散。DP-LUR和HLU-SGS格式通过在边界单元采用雅克比迭代算法,能有效消除分区边界影响,使计算格式在大CFL数条件下保持稳定,其中HLU-SGS继承了LU-SGS迭代效率高的特点。在此基础上,采用并行化的LU隐式方法对Caradonna-Tung旋翼进行了并行计算,针对两套稀密程度不同的背景网格,在产生复杂分区边界的计算网格条件下,均获得了气动力和旋翼尾迹都充分收敛的流场。通过计算,分析了背景网格密度对预测气动力性能和捕捉尾迹流场的影响,计算结果与试验值吻合良好,验证了当前隐式计算方法适用于旋转坐标系并行计算,适合于推广至其他大规模并行分区的流场计算。     

亚跨超CFD软件平台中的网格生成技术

在亚跨超ＣＦＤ软件平台中,提供了一个网格生成平台。利用它,可以快速生成常规兵器、常规导弹和正常布局的常规战斗机的计算网格。有多种拓扑结构和多种方法可供选择。结果表明,生成的网格实用,质量令人满意。气动特性的计算结果也是可靠的。     

基于CFD/CSD耦合的叶轮机叶片失速颤振计算

为了研究叶轮机叶片的失速颤振特性,发展了一种计算流体力学与计算结构力学(CFD/CSD)时域耦合方法。该方法通过每一物理时刻CFD和CSD的循环迭代实现了耦合计算。在CFD分析中,采用鲁棒性较好的空间离散格式AUSM+-UP,并基于延迟脱体涡模型(DDES)模拟了带分离流动。在结构分析中,通过模态法构建了旋转叶片动力学方程并运用杂交多步方法进行求解。以孤立转子Rotor37为例,计算了不同工况下流场总体与细节参数,与实验结果的对比验证了CFD算法的精度。对某转子叶片进行了颤振特性研究,计算所得的广义位移时间响应曲线表明该叶片在近失速工况下会发生失速颤振,其表现形式为一阶弯曲模态发散且各阶模态之间不耦合。分析表明,流场不稳定和非定常效应是引起失速颤振的关键因素,同时折合频率的降低也会导致原本气动弹性稳定的叶片发生失速颤振。     

人员离机轨迹数值模拟及快速预测研究

离机轨迹研究在分析外挂物分离、导弹投放以及人员空中脱离过程等方面有着重要的参考意义。运用动态嵌套网格技术,通过CFD方法来模拟人员空中离机后的运动轨迹。提出一套基于本征正交分解(POD)方法的人员离机轨迹快速预测算法。将CFD数值模拟得出的离机轨迹作为样本,以飞行马赫数和攻角双参数变化为例,详细介绍POD方法预测离机轨迹的实现步骤。结果表明:POD方法相对于CFD方法可以大大减少计算时间,并且POD预测结果与CFD方法得到的结果吻合较好,为离机轨迹预测在工程应用上提供了一种快速简便的方法。     

三维机翼的静气动弹性特性数值模拟研究

发展了一种计算流体动力学（CFD）和计算结构动力学（CSD）的耦合计算方法,对三维机翼的静气动弹性进行了数值模拟研究。采用三维欧拉方程为控制方程基于直角网格计算气动力,并耦合结构静平衡方程进行静气动弹性数值模拟,设计了CFD/CSD耦合计算的数据交换的方式。以M6机翼作为算例,进行了机翼静气动弹性的数值模拟,计算结果表明所发展的三维机翼静气动弹性数值模拟方法是合理可行的,并可作为机翼结构优化设计和考虑结构弹性变形影响的气动外形优化设计的基础。     

开放式气动力数值模拟系统研究

气动力数值模拟系统是ＣＦＤ流场解算技术、网格生成技术、数据可视化技术和网络技术相结合的产物。描述了一种运行于计算机网格环境下多任务、多用户、交互式的数值模拟软件环境－ 气动力数值模拟系统,该环境能计算模拟亚、跨、超音速绕流的气动力,也能方便地集成模拟不同复杂气动布局飞机流场的解算程序。介绍了其网格生成、数据可视化、系统接口等核心部分的结构和功能。     

多块网格网络并行计算中的负载分配研究

针对CFD中多块网格计算的特点，并使用MPI网络并行系统，对某战斗机绕流进行了基于三维Euler方程的并行计算，主要研究了多块网格并行计算中负载的分配方法，发展了负载自动分配和网格自动重分区程序，计算结果表明：并行计算结果和实验结果完全吻合，8个节点机的并行效率达到了89％。     

一种基于CFD的叶轮机非定常气动力组合建模方法

为了获得一个准确高效的非定常空气动力学模型并将其应用于叶轮机叶片颤振特性分析中去,论文发展了一种基于CFD方法的叶轮机非定常气动力组合建模方法,可以快速计算叶轮机叶片在等相角差振动时的气动阻尼系数。运用小扰动流场的叠加原理,通过不同通道数模型的非定常流场求解(通常需要两次或三次),针对流场的周期性边界条件,组合分析得到一系列更多通道数情况下的非定常气动力低阶模型。基于这种降阶模型计算的气动阻尼系数与直接的CFD方法计算结果吻合很好,计算效率提高10倍以上。