DLR-F6翼身组合体跨音速流场

针对跨音速运输机经典算例DLR-F6翼身组合体模型,采用CFD方法对其气动特性进行了粘性流动数值模拟,流动模型为雷诺平均NS(RANS)方程。首先采用“超立方体”概念生成绕DLR-F6翼身组合体的高质量多块结构拼接网格,研究网格拓扑结构对气动特性的影响,在此基础上通过网格细分和粗分考查了网格密度对计算结果的影响,最后进行了湍流模型的影响研究。通过与实验数据对比分析,得出了适宜DLR-F6翼身组合体跨音速粘性流动的计算网格,并总结出了能较好模拟其跨音速流场特性的湍流模型。结果表明：网格拓扑结构的合理设计会对计算结果产生一定的影响;网格密度对机翼表面压力分布没有明显影响,但对阻力系数影响显著;湍流模型对机翼表面压力系数分布的影响主要体现在激波位置上,对翼根处的分离也有一定的影响;SST 模型计算的气动力系数比SA模型接近实验值。     

DLR-F6翼身组合体跨声速阻力计算

采用美国航空航天学会阻力测试小组提供的多块对接网格,结合Spalart-Allmaras、Wilcox’sk-ω和Menter’s k-ω SST三种湍流模型,通过求解雷诺平均Navier-Stokes方程、数值模拟DLR-F6翼身组合体的流场来研究阻力计算精度,考查网格和湍流模型对翼身组合体构型气动特性的影响.结果表明:三种湍流模型得到的机翼表面压力系数分布与实验数据吻合良好,气动力随攻角的变化趋势与实验结果一致;Spalart-Allmaras模型得到了网格收敛结果,所得阻力优于其他软件的结果;网格密度对阻力有影响,对机翼表面压力系数分布无明显影响;湍流模型对机翼表面压力系数分布的影响主要体现在激波位置上,对升力影响较小,对阻力(尤其是摩擦阻力)影响显著,对翼根处的流动分离有一定影响;在跨声速流动中,Menter’s k-ω SST模型的结果最接近实验数据.     

翼身组合体定常亚音速升力特性的数值研究

用有限基本解方法研制成翼身组合体升力和力矩的计算程序。以当量回转体代替机身,机翼可以有后掠、厚度及弯度。厚度问题用线源模拟,先予解决。然后求解升力问题,这时考虑了厚度的影响。此法简单且有一定的精度。本文还计算了一些实例,并与实验作了比较。对机身影响、机翼的后掠和根梢比等作了简要的讨论。     

民用飞机翼身整流粘性绕流数值模拟方法研究

利用CFD方法,针对民用飞机开展翼身整流黏性绕流数值模拟方法研究.采用Roe三阶迎风偏置通量差分分裂方法(FDS)和隐式近似因子分解方法(AF)对三维非定常雷诺平均N-S方程进行数值离散求解,采用多重网格技术加速流场求解收敛过程,并比较了不同湍流对翼身结合处流动分离的模拟能力.计算结果与风洞实验的对比表明,所采用的数值方法可以有效地模拟翼身结合处黏性分离流动,具有一定的工程实用价值.     

弧形翼-身组合体绕流流场数值模拟及实验研究

为了研究标准技术协作计划(TTCP)弧形翼-身组合体的空气动力学特性,对其在超声速下的三维绕流流场进行了数值模拟.应用雷诺应力湍流模型和平流迎风劈裂法离散格式,对弧形翼-身组合体在Ma∈[2.0,5.5]范围内的超声速绕流流场结构进行了分析.计算分析了弧形翼-身组合体的各气动力系数随马赫数、攻角的变化规律.阐述了弧形翼在零攻角状态下产生自诱导滚转力矩的机理.利用超声速风洞进行了标准TTCP弧形翼-身组合体在Ma=3.0条件下的吹风实验.实验结果与数值模拟结果的计算误差小于10％,验证了所采用数值模拟方法的可信度.     

翼-身组合体的网格生成与N-S方程数值模拟

以DPW4的CRM模型为例,给出了一种复杂外形飞行器网格生成的分块对接网格拓扑方案以及实现过程,在无限插值法生成的初始网格基础上,用椭圆微分方程生成空间网格。基于该网格求解N-S方程,获得了较为合理的全机气动力系数与压力分布。     

多翼空空导弹复杂流场气动力计算

为了高精度地确定导弹的气动力参数,为导弹合理的气动布局提供理论依据,对绕某外形多翼空-空导弹全弹体三维绕流流场进行了雷诺平均Navier-Stokes方程数值模拟,为了使方程封闭采用了Spalart-Allmaras一方程湍流模型,计算方法为Jameson的有限体积多步Runge-Kutta时间推进法和空间中心差分耦合人工粘性项,使用分区技术,通过求解椭圆型偏微分方程方法生成了连续拼接式多块结构化网格,在超音速带有攻角和有滚转角或无滚转角的条件下,获得了导弹弹体、舵翼和副翼上的气动力和力矩,计算结果与实验数据以及其它程序的计算结果吻合较好.结果表明:文中算法可用于真实复杂外形导弹、复杂流态下的气动力计算.     

翼身融合飞机结构研究

主要研究了翼身融合飞机结构受力,讨论了BWB结构与传统飞机结构最主要的区别,提出了一种新的BWB非圆形截面结构.通过计算分析,比较了圆形截面、矩形截面以及混合型截面三种不同结构受力.计算结果表明,在满足重量限制的条件下,混合型结构是一种高效的非圆柱BWB结构.该结构减小了结构所承受的弯曲应力,提高了BWB飞机结构的安全性.     

航空CFD高性能计算应用研究

应用万核级并行流场数值模拟软件CCFD,在深腾7000高性能计算机上对飞机高速巡航构型阻力和低速增升装置性能预测进行了计算研究,比较了不同规模网格、湍流模型对计算结果的影响,并且将高性能计算与常规计算的工作效率进行了比较.飞机高速巡航构型阻力和低速增升装置性能预测的计算结果与试验结果较为接近,可见:CCFD软件的常规气动计算精度较高;高性能计算技术可有效提高数值计算的可靠性及工作效率;单纯增加网格规模,而不进行算法升级,对复杂问题计算精度的提高有限.     

大冲角翼栅绕流流场及压力场的边界元方法

采用粘性流场边界元方法，对大冲角的翼栅绕流问题进行了研究。得到了速度场、内点的压力分布及翼栅边界面上的压力分布。为了提高计算精度，在对积分方程离散化时，加密了刻分网格，并且用混合单元方法对角点进行处理。计算方法简洁结果合理，适应性较强，显示了本文方法的优越性。     

“钻石背”弹翼颤振的数值计算

研究"钻石背"弹翼的颤振特性.采用Newmark算法简化颤振方程,耦合结构模型和非定常Navier-Stokes模型在时域内求解颤振方程.使用FLUENT 6.2软件求解非定常Navier-Stokes方程,得到非定常气动力;使用NASTRAN软件弹翼进行结构分析.用AGARD 445.6翼算例验证了计算方法的可行性,进行了"钻石背"弹翼的颤振特性计算,得到了颤振速度和颤振频率.     

MPC中分块耦合与边界正交的贴体网格生成

在计算ＭＰＣ模块多分支联接的流场时，网格划分的分块耦合技术可使网格线与边界成任意角度，这样既能实现网格线与边界正交，也能保证在分界处网格光滑连续，还可实现对内场网格的疏密控制．应用贴体网格生成方法，分别生成了引射喷管和排气总管的二维网格，证实了该方法的可行性     

用改进非结构动网格方法模拟跨音速非定常绕流

对用于非结构动网格生成的弹簧近似方法进行了研究.通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入防扭转、防挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进.改进后的方法可以大大提高网格变形能力和网格质量.应用本文发展的非结构动网格生成方法并通过耦合求解基于ALE(Arbitrary Lagrangian-Eulerian)描述的Euler方程,模拟了做谐和振动的矩形刚性机翼及大展弦比后掠机翼弯扭耦合振动跨音速非定常绕流,计算结果与参考文献提供的结果及实验结果吻合良好.     

MPC模块三维分块贴体网格的耦合生成

采用分块耦合技术生成了ＭＰＣ模块多分支联接的三维非正交贴体网格，为进行流场计算作了必要的准备．所介绍的网格生成方法保证了在分块生成网格时穿过分界面的网格线是光滑连续的，并且可通过边界面上网格疏密或者规定边界面上网格线之间的夹角来控制内部区域的网格疏密．     

动力系统流场计算动态网格生成模型研究

运用网格分块生成与整体集成方法,提出了复杂动力机械系统内流场三维流动计算网格自动生成模型,并以K100摩托车发动机扫描气系统为对象,实现了换气过程瞬态数值模拟动态网格的自动生成,数值计算结果,该方法可用于复杂机械系统内流场三维瞬态流动的数值计算。     

全隐式E-CUSP迎风格式大攻角分离流数值模拟

采用全隐式、低耗散E-CUSP格式,通过求解雷诺平均Navier-Stokes方程耦合Spalart-Allmaras(SA)湍流模型,模拟了细长旋成体在超声速、大攻角下的流场,分析了背风面涡的发展过程.结果表明:E-CUSP格式耦合SA湍流模型能够准确地模拟背风面的流动分离和精细的二次涡,对横向分离具有较高的模拟精度;预测的物面压力系数分布和激波位置与实验数据吻合良好,力和力矩的相对误差在1.98%之内;E-CUSP格式可用于模拟复杂的分离流动,具有高的计算精度和效率.     

构型配分函数的简化计算与讨论

统计物理教材中分子构型配分函数的计算,一般是将双矢量积分化简为单矢量积分。文章对该双矢量积分进行了详细推导,并采用Lennard—Jones势模型进行数值计算,分析讨论了构型配分函数简化计算的适用范围。     

同位网格有限体积法在变密度浆液流场计算中的应用

为研究自膨胀浆液的扩散机理,基于现代计算流体动力学理论,建立了一种变密度浆液二维流场计算方法.该方法采用结构化同位网格剖分计算区域,利用有限体积法离散浆液流动控制方程,借助动量插值技术构建压力修正方程;采用压力耦合方程组的半隐式算法(SIMPLE)迭代求解动量离散方程和压力、速度修正方程.通过变密度浆液在二维矩形狭槽中自由膨胀算例对该算法进行检验.结果表明,该方法求解正确,并具有较高的计算精度,为建立膨胀性高聚物注浆材料在二维裂隙中流动扩散仿真分析方法奠定了基础.