DLR-F6翼身组合体跨声速绕流的CFD计算

合集下载

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

C2&O" 能求解二 P 三维雷诺平均纳维尔 B 斯托克斯方程E 可进行基于区域分裂策略的高性能并行计算。它调用了基于消息传递机
该文为“!""# 东西方高速流场会议”论文。— ——编者注
航空科学技术·!""# 年第 $ 期
学术交流
航空科学技术
制的 %&’ 并行包来实现并行化，并调用了开放源码的 %()*+ 库函数对混合网格进行负载平衡的分区。
航空科学技术·!""# 年第 $ 期
五、计算结果所有工况所采用的马赫数都
是 "IJ#，雷诺数为 D(=F &K"$4（基于平均气动弦长 KLKI!MME。每个计算网格都是半模，对称面上采用了对称边界条件。其他格点格式求解器 0CD:N@O 、8058 的结果都
来自公开的文献。对于 ,; 模型的中尺度网格，
一、标准模型数据介绍 ’@AB&$+ 模型是一种双发动
机的宽体飞机，多年来已是多种风洞实验和数值计算的研究对象。围绕它已经针对发动机的多种尺度和安装位置进行了物理实验。阻力预测会议要求的实验马赫数为 "HI#，目标升力系数是 "H#""。这里采用了第二次阻力预测会议组委会公布的几何数据以及相应工况的实验数据。
三、湍流模型
计算采用了混合显式代数应力
模型的 !: ! 两方程湍流模型。这
是一个非线性涡粘性模型<4 它通过
求解两个关于 ! 和 ! 的输运方程
来得到涡粘性 " "。这个模型是直接由完整的雷诺应力方程通过简化

DLR—F6外形计算网格及湍流模型影响因素研究

DLR—F6外形计算网格及湍流模型影响因素研究作者：孙悦李仁府熊有德周玲王亮来源：《航空兵器》2017年第05期摘要：本文选用DLR-F6翼身组合体模型，分析了不同类型网格及湍流模型对机翼表面压力分布和翼根分离区的预测精度影响。

分析结果表明，六面体、四面体和多面体网格预测得到的机翼表面压力分布和翼根分离区大小基本一致，在保证相同计算结果精度的前提下，多面体网格使用的网格单元数量最少，计算效率最高，且网格生成十分方便。

此外， S-A，SST和RSM湍流模型均能较准确地预测出DLR-F6机翼表面压力分布，但S-A和SST湍流模型预测得到的翼根处分离区较实验结果明显偏大，而RSM湍流模型预测结果与实验结果更加吻合。

从湍流模型构造上分析发现， S-A和SST模型基于湍流各向同性假设，忽略了角区分离流动处的雷诺正应力之差，而RSM湍流模型由于反映了雷诺应力的各向异性，因此预测得到的翼根处分离区与实验结果更加接近。

关键词： DLR-F6；计算网格；湍流模型；多面体网格；分离区中图分类号： V211.3； O354.4文献标识码： A文章编号： 1673-5048（2017）05-0060-080引言为了更好地评估计算流体力学（Computational Fluid Dynamics， CFD）技术在亚音速飞机阻力预测方面的应用情况，从2001年到2016年，美国航空航天学会（American Institute of Aeronautics and Astronautics， AIAA）应用空气动力学技术委员会（Applied Aerodynamics Technical Committee）共组织召开了六次阻力预测会议（Drag Prediction Workshop， DPW）[1]。

会议邀请全世界的高校、研究院所、工业部门参加，前三届会议统一选用德国宇航公司（Deutsches Zentrum für Luftund Raumfahrt， DLR）的翼身组合模型，从第四届会议开始，改用波音公司和NASA共同设计的CRM（Common Research Model）模型。

三种湍流模型在跨声速绕流中的计算精度分析

三种湍流模型在跨声速绕流中的计算精度分析郑秋亚;左大海;刘三阳【摘要】通过求解雷诺平均Navier-Stokes方程考察Spalart-Allmaras、Wilcox's k-ω和Menter's SST三种湍流模型在跨声速流动中的计算精度.结果表明:Menter's SST模型预测的力和力矩最接近实验数据;Spalart-Allmaras模型的压力分布和激波位置与Menter,s SST模型的一致,Wilcox's k-ω模型的激波位置相对偏后,且预测的升力和摩擦阻力偏大.【期刊名称】《弹箭与制导学报》【年(卷),期】2011(031)005【总页数】3页(P152-154)【关键词】计算精度;Navier Stokes方程;跨声速流动;湍流模型【作者】郑秋亚;左大海;刘三阳【作者单位】西安电子科技大学理学院,西安710071;长安大学理学院,西安710064;长安大学理学院,西安710064;西安电子科技大学理学院,西安710071【正文语种】中文【中图分类】V211.30 引言在雷诺平均 Navier－Stokes（RANS）方程计算中，计算结果的精度不仅与所用格式的精度有关，而且与湍流模型对流动的模拟精度有很大的关系。

目前工程中应用最为广泛的湍流模型主要包括：Spalart－Allmaras（SA）一方程模型[1]，Wilcox’s k－ω[2]（Wilcox）两方程模型和Menter’s SST[3]（SST）两方程模型。

长期以来，人们在工程应用中发现不同湍流模型表现出了不同的特性，因此研究湍流模型对流动的模拟精度，对CFD模型的选取和CFD计算精度的提高有着重要的意义。

文中以ONERA M6机翼和DLR－F6翼身组合体[4]为例，采用具有高精度和稳定性良好的Roe格式[5]，通过求解RANS方程分析SA、Wilcox和SST湍流模型对跨声速流动的模拟精度，为更高准确度湍流模型方程的建立和CFD模型的选取提供参考。

DLR-F4翼身组合体跨声速绕流CFD计算

DLR-F4翼⾝组合体跨声速绕流CFD计算☆作业⼆☆DLR-F4翼⾝组合体跨声速绕流的CFD计算姓名：涂飞学号：SZ1101020⽬录⼀．本作业概述 (1)1. 所⽤模型 (1)2 来流条件 (1)3 ⽹格划分 (2)⼆. 求解结果 (3)1 流场云图 (3)2动⼒学参数 (3)3截⾯压强系数分布 (5)三．总结 (7)参考⽂献 (8)图表1 DLR-F4模型（1） (1)图表2模型参考数据 (1)图表3 来流参数 (1)图表4 ⽹格划分 (2)图表5 物⾯压强系数云图 (3)图表6 动⼒学参数 (3)图表7 升⼒系数计算结果和实验数据对⽐图 (4)图表8 阻⼒系数计算结果和实验数据对⽐图 (4)图表9俯仰⼒矩系数计算结果和实验数据对⽐图 (4)图表10 升阻⽐计算结果和实验数据对⽐图 (4)图表11 极曲线计算结果和实验数据对⽐图 (4)图表12 残差收敛历史 (4)图表13 Cl，Cd，Cm收敛历史 (5)图表14 机翼7个不同位置截⾯图 (5)图表15 机翼各个截⾯表⾯压强系数分布与实验数据对⽐图 (7)⼀．本作业概述1.所⽤模型DLR-F4翼⾝组合体，如图：图表 1 DLR-F4模型（1）参考⾯积0.1454 m2参考长度C ref141.2 mm参考半展长b ref /2 585.647 mm参考点（CAD 坐标）x = 157.9 mm, z = -33.92 mm 图表2模型参考数据2 来流条件⽓体属性理想，粘性来流温度T∞255.6 K来流压强P∞101325 pa来流雷诺数Re 3×106来流马赫数Ma 0.75攻⾓（dge）-1°， -0.5°， 0°， 0.5°， 1°图表 3 来流参数3 ⽹格划分⽹格要求：粘性⽹格，y＋≈30图表 4 ⽹格划分本作业我采⽤块结构⽹格进⾏划分，⽹格单元600万，⽹格划分⽤Pointwise软件完成。

_DLR-F6复杂组合体跨声速阻力计算研究

， F6 WB 和 WBNP 对 DLR-
组合体流场进行了全湍流和固定转捩 2 种边界层流动处理方式的数值模拟，以分析阻力对网格尺寸的敏感性，评估挂架和吊舱引起的阻力增量，考察转捩对阻力和阻力增量的影响．
1
控制方程
控制方程为与时间相关的三维守恒型 RANS 方
程，在一般曲线坐标系下其形式为
－ 2 /3 分别给出了 WB 和 WBNP 的总阻力系数随着 N
a 为当地声速，式中： Re 为雷诺数， μ T 为湍流粘性系数．计算采用格心式有限体积方法，无粘项使用高分辨率的 Roe 格式耦合保单调性的 MUSCL 插值方 SGS 时间推进方法，法，时间采用隐式 LU湍流模型采用鲁棒性良好的 SA 模型，其模型方程详见文献［ 10］．
．
本文以 DPW Ⅱ 为背景，通过求解雷诺平均 NavierStokes （ RANS ）方程耦合 SpalartAllmaras （ SA）一方程湍流模型
［10 ］
处，上表面 10% 弦长处；计算网格选用 DPWⅡ提供的 ICEM 1to1 结构网格， F6 WBNP 组合体表图 1 是 DLRF6 WB 的表面网格．面网格，取掉挂架和吊舱为 DLR-
^ （ F － F v ）（ G － G v ） Q （ E － E v ） + + + = 0． t ξ η ζ ^ ^ ^ ；式中： t 为时间； Q 守恒变量； E、 Ev 、 Fv 、 G v 为粘性矢通量， Ev 、 Fv 、 G v 中包含的应力张量分量和热通量分别由下式给出： τij = （ μl + μT ） M∞ Re
［69 ］
C L = 0． 5 粗、算例 1 ：网格细分研究 M ∞ = 0 ． 75 ， WB ， WBNP，中、细网格，全湍流； CL = 算例 2 ：转捩对阻力的影响 M ∞ = 0 ． 75 ， 0 ． 5 中网格， WB ， WBNP，固定转捩．

CFD课程设计说明书翼身组合体流场分析

CFD课程设计翼身组合体流场分析院系：航空航天工程学部专业：飞行器设计与工程班级： 24030301学号： 2012040303023姓名：摘要此次课程设计是利用ANSYS软件中的ICEM和Fluent求解器计算不同迎角下，翼身组合体的升力系数，阻力系数，力矩系数以及各个状态下的流场分布情况，机身为方截面机身，机翼为三角上单翼，翼型选择NACA4412，计算结束后，利用tecplot软件绘制Cy-α，Cy-Cx，Mz-Cy曲线，得出Cy0，最大升阻比等气动力特征参数。

关键词ICEM Fluent 翼身组合体tecplot目录第一章绪论 (1)1.1 ANSYS软件介绍 (1)1.2主要内容 (1)第二章模型的建立 (2)2.1 CATIA建立模型及导出 (8)第三章ANSYS.ICEM处理 (4)3.1 导入模型 (4)3.2 网格划分 (4)3.3 导出网格 (8)第四章Fluent计算 (9)4.1 设置参数计算 (9)4.2 计算结果 (12)第五章数据处理分析 (18)4.1气动参数曲线 (18)参考文献 (21)第一章绪论1.1 ANSYS软件介绍ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件，是一个多用途的有限元法计算机设计程序，可以用来求解结构、流体、电力、电磁场及碰撞等问题。

由世界上最大的有限元分析软件公司之一的美国ANSYS开发，它能与多数CAD软件接口，实现数据的共享和交换，如Pro/Engineer，NASTRAN，Alogor，I－DEAS， AutoCAD等，是现代产品设计中的高级CAD工具之一。

在此次的课题中，主要用到其中的ICEM及Fluent部分。

1.2 主要内容本次课程设计的主要内容就是通过CATIA建立机身和机翼的组合体模型，通过fluent解算器进行有限元分析，从而得到该组合体的一些相关的气动数据。

此次课程设计的重点在于模型的建立，通过CATIA建立基础的模型，然后导入到ANSYS.ICEM中进行模型的处理以及网格包括壳网格、体网格及附面层网格的划分。

cfd-17

文章编号: (2009)-17大客设计中的结构网格划分规范研究刘杰平1张淼2张玉东1张卫民11中国航天空气动力技术研究院，北京，1000742上海飞机设计研究所，上海）摘要：本文选取AIAA阻力预测工作组（DPW-II）所提供的翼身组合体DLR-F6为研究对象，它同时也是大型客机设计中最基本和最典型的外形。

为满足工程设计对结果精度和计算周期的双重要求，在计算机硬件条件和湍流模型、数值求解器以及具体参数设置都相同的前提下，网格技术就显得尤为关键。

本文根据实践中摸索出的原则所独立完成的网格，与DPW-II提供的三套密度由粗到细的标准网格，分别对DPW-II所要求的相同状态进行了认真的计算和分析。

结果表明：本文所使用的网格技术更适于大型客机工程设计的需求，同时也为进一步研究更复杂和更真实的大型客机奠定了基础。

关键词：大型客机；CFD；网格；DPW；DLR-F60 引言大型客机是国家“十一五”计划的重点支持项目。

在大型客机的方案选型和初样设计阶段，计算流体力学（CFD）发挥着其独有的重要作用。

在气动优化设计、雷诺数修正以及局部流场精细刻画等方面，CFD较风洞试验具有成本低、效率高的天然优势；但如何使CFD更好的满足大型客机工程设计对计算精度和周期的双重要求，这是CFD工作者必须面对和解决的问题。

计算精度主要体现在阻力的数值预测上，为研究阻力的计算精度问题，AIAA阻力预测工作组（DP W）于2001年和2003年先后召开了两次工作会（DPW-I和DPW-II），分别针对两组外形（DLR-F4翼身组合体和DLR-F6翼身组合体及加吊挂/发动机舱的组合体），各提供三套密度由粗到细的网格展开盲比[1,2]。

相关试验是20世纪90年代在法国ONERA S2MA 1.77m×1.75m跨声速风洞中完成的。

两次会议都有近20家单位使用近20种软件给出了计算结果，尽管与试验结果的比较情况千差万别，但基本反映出了CFD 对阻力计算精度问题的研究现状。

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

ｔｏｏｋｈｐｉｅｌｙｄｉｈｏｕａｉｎｉｎＷｒｓｏｓｍｐｏｅｎｔｅｃｍｐｔｔｏ．ＴｈｏｕａｉｎｌｅｕｔｈｗｈｔｔｅｐｅｓｒｅｃｍｐｔｔｏａｓｌｓｓｏｔａｈｒｓｕｅｒｄｓｒｂｔｎｎｔｅｗｉｇｓｒａｅａｒｅｗｅｌｔｈｘｅｉｅｔＴｈａｙｎｒｎｓｏｈＯＴ— ｉｔｉｕｉｓｏｈｎｕｆｃｇｅｌｗｉｈｔｅｅｐｒｏｍｎ．ｅｖｒｉｇｔｅｄｆｔｅＣＩＩｐｔｄａｒｄｎｍｉｏｃｓｗｉｈｔｅａｇｅｏｔａｋａｅｃｎｉｔｎｔｈｘｅｉｅｔｕｅｅｏｙａｃｆｒｅｔｈｎｌｆｔｃｒｏｓｓｅｔｗｉｈｔｅｅｐｒｍｎ．Ｇｒｄｒｆｅａｉｅｉ — ｎ
（．１西安电子科技大学理学院，１０１７０７，西安；．２长安大学理学院，１０４西安；７０６，３中国航空计算技术研究所航空气动力数值模拟重点实验室，７０６，西安）．１０８
摘要：采用美国航空航天学会阻力测试小组提供的多块对接网格，合Ｓａａｔｌｒｓ结ｐｌｒＡｌａ、～ｍａ
ｖｎｙｏｄ－ｖｒｇｄＮａｉｒＳｏｅ（ｉｇＲｅｎｌｓａｅａｅｖｅ－ｔｋｓＲＡＮＳｅｕｔｏｓｃｕｌｄｗｉｈｔｒｅｋｎｓｏｕｂｌｎｅ）ｑａｉｎｏｐｅｔｈｅｉｄｆｔｒｕｅｃｍｏｅｓｔｅＳａａｔＡｌａａＳ）ｏｅｅｕｔｏｄｌｈｉｏＳ ∞ ｍｏｅ，ａｄｔｅＭｅｔｄｌ－ｈ ห้องสมุดไป่ตู้ ｐｌｒ— ｌｒｓ（Ａｍｎｑａｉｎｍｏｅ，ｔｅＷｌｘ’ ｃｄｌｎｈｎ —

DLR_F6翼身组合体阻力计算

[文章编号] 1001-246X(2008)02-0145-06[收稿日期]2006-11-14;[修回日期]2007-03-17[作者简介]王运涛(1967-),男,黑龙江密山,研究员,博士,从事计算空气动力学方面的研究.DLR -F 6翼身组合体阻力计算王运涛, 王光学, 张玉伦(中国空气动力研究与发展中心计算空气动力学研究所,四川绵阳 621000)[摘要] 采用/亚跨超CFD 软件平台0(TRIP210)数值模拟DLR -F6翼身组合体构型,采用的多块对接网格、测压和测力的试验结果均来自美国AIAA 阻力计算小组,对比计算结果采用CFL3D 的结果.详细研究网格密度、湍流模型对DLR -F6翼身组合体构型的总体气动特性和压力分布的影响,计算结果与相应的试验结果较一致.采用SST 两方程模型得到网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布影响较小.[关键词] TRIP210;DLR -F6;DPWII;阻力;湍流模型[中图分类号] V21117[文献标识码] A0 引言随着计算机软硬件技术和计算流体力学(CFD)的飞速发展,CFD 在航空航天飞行器方案选型和初样设计阶段发挥的作用愈来愈重要,在雷诺数修正、优化设计、细致流场刻画等方面,CFD 较风洞试验更具有优势,已经与风洞试验一样,成为飞行器气动设计最重要的研究手段[1].随着大型网格前置处理软件、C FD 软件系统和后置处理软件的发展,CFD 几乎可以模拟所有高度复杂飞行器外形的绕流流场[2].飞行器设计者使用CFD 工具面临的首要问题是软件的精度和效率问题.CFD 软件的计算效率问题依靠多重网格技术和并行计算技术已得到较好解决,但计算精度问题(Verification &Valida tion)依然是当前研究的热点[3].为此,国际上先后组织了许多专题研讨会,如欧洲的计算空气动力学研究项目ECARP[4](European C omputational Aerodynamics Research Project)和AI AA 的DPW(Dra g Prediction Workshop).EC ARP 项目的目标是对欧洲的软件系统进行统一的确认研究,并对各种湍流模型进行评估.为研究C FD 的阻力计算精度问题,AI AA 阻力计算工作小组在2001年6月召开了第一次阻力计算的工作会议(DPW I),该次会议选择DLR -F4翼身组合体作为标准算例.会上18家单位提供了采用14种软件的计算结果[5].2003年6月召开了第二次阻力计算的工作会议(DPW II),该次会议选择DLR -F6翼P 身P 挂P 舱组合体作为算例,会议的重点是阻力计算精度,共有22家研究机构提供了20种CFD 软件的计算结果[6].试验结果是90年代在法国ONERA S2MA 1177m @1175m 跨声速风洞中完成的.近年来,国内也先后组织了全国范围内的CFD 软件可信度确认工作,如2004年8月中国空气动力研究与发展中心与中国航空工业六三一所联合举办的/CFD 统一算例研讨0活动和2005年6月中国空气动力研究与发展中心举办的CT-1大攻角气动特性研讨会.本文研究目的是考核中国空气动力研究与发展中心(C ARDC)自行研发的CFD 软件TRIP210(TRIsonic Platform Version210)的数值模拟精度.软件的确认算例选择了DLR -F6翼身组合体外形,利用DPWII 提供的3套粗细不等的结构网格,首先详细研究了网格密度、湍流模型对翼身组合体构型典型气动特性的影响,在此基础上,数值模拟了固定马赫数下的极曲线,对比的试验结果包括DPWII 提供的测力和测压的试验结果以及Rumsey 等人采用CFL3D 得到的计算结果[7].1 TRIP210软件简介TRIP210是中国空气动力研究与发展中心自行研发的C FD 软件,该软件采用结构网格技术和有限体积第25卷第2期2008年3月计算物理C HI NESE JOURNAL OF C OMP UTATI ONAL PHYSICSVol.25,No.2M ar.,2008方法,通过数值求解三维任意坐标系下的Reynolds-Averaged Navier-Stokes(RANS)方程,主要面向以下几类问题的数值计算:全机中等攻角以下纵、横向气动特性计算;飞机全机进排气系统一体化匹配计算;喷管、进气道内部流动的模拟及特性计算;导弹全弹纵、横向气动特性计算;航天器亚跨超气动特性数值模拟.TRIP210软件已经具备的功能如下:1)多块对接P拼接网格结构流场计算能力,初步具备对接P重叠结构的计算能力;2)多窗口技术和多重网格加速收敛技术;3)LU-SGS分解、SSOR求解和四步Runge-Kutta法,多种五点差分格式和七点差分格式接口;4)多种湍流模型,Spalar-t Allmaras一方程和Menter.s k-Omega SST等多种湍流模型;5)网格文件标准的Gridgen接口和ICE M接口,后置处理标准的Tecplot接口;6)开发了TRIP210软件的工具库和算例库,以及简洁实用的用户界面.2DLR-F6翼身组合体计算网格和计算状态本文采用的计算网格均来自DPWII,该结构网格由ICE M软件生成,网格结构为多块对接网格(1-to-1),分为粗网格、中等网格和细网格3种,网格详细信息如表1所示,DLR-F6翼身组合体的计算构型和表面网格分布(中等)见图1.由表1可以看出,在该网格序列中,各套网格之间并没有2倍数关系,这也许会给网格收敛性研究带来一定的影响.表1标准网格统计列表Table1Grid statistics for standard grids网格类型网格节点网格单元计算分区Coarse3579589337484821 Medium6038635571596827Fine10421976996659232图1DLR-F6翼身组合体的表面网格(中等)Fig11Surface grids for DLR-F6wing P body(medi um)本文计算状态如下[6]:算例1:网格收敛性研究.M=0175,Re=310@106(基于M.A.C=011412m),C L=01500?01001,全湍流计算(fully turbulent).算例2:极曲线.M=0175,Re=310@106(基于M.A.C=011412m),A=-3b,-2b,-115b,-1b,0b,1b, 115b,全湍流计算(fully turbulent).其中M.A.C为平均气动弦长.3计算结果与分析采用TRIP210软件,选择二阶精度的通量差分(FDS)类型的MUSCL差分格式,Spalart-Allmaras(SA)一方程和Menter.s k-X SST(SST)两方程湍流模型,模拟上述两种工况的流场,采用3~4重网格加速收敛.测力与测压试验是1990年8月在法国ONERA S2MA风洞中完成的,对比计算结果采用Chris Rumsey等人采用CFL3D软件得到的结果[7].选择这一计算结果,首先是由于该计算采用的网格与本文相同,可以避免网格因素的影响,其次是因为CFL3D软件的基本功能与本文采用的软件类似.需要说明的是,CFL3D在算例2计算中采用了DPW II建议的转捩位置,而本文采用的是全湍流计算.311算例1计算结果及分析采用DPWII提供的3套(Coarse,Medium,Fine)多块对接网格和SA与SST两种湍流模型,数值模拟了算例1的工况,主要目的是研究网格收敛性问题和湍流模型的影响两个方面.图2给出采用粗网格、中等网格和密网格3套网格以及两种湍流模型得到的网格收敛性研究结果,同时146计算物理第25卷还给出CFL3D 采用SA 模型的全湍流计算结果,气动横坐标为网格节点的(-2P 3)幂,纵坐标为阻力系数.由图2可以看出,本文采用SST 两方程模型得到了网格收敛性结果,而采用SA 模型时,中等网格计算得到的阻力系数偏大.SST 模型从粗网格到密网格阻力变化只有3Counts(1Counts=10-4),而SA 模型的变化则大很多,说明SST 模型的计算结果对网格的依赖性较弱.网格收敛性研究常用的策略是通过在物面法向加密使得网格序列之间有2倍数的关系,但本文采用的网格序列之间不具有这种关系(表1),如何设计网格序列,研究网格收敛性是一个值得研究的问题.图2 DLR-F6翼身组合体的网格收敛性研究Fi g 12 Grid convergence for DLR -F6wing Pbody表2给出在固定升力系数下,采用三套网格、两种湍流模型的详细计算结果,包括相应的攻角、总阻力、压差阻力、摩擦阻力和力矩系数等,同时给出的还有C FL3D 采用SA 模型的计算结果和相应的试验结果.由表中可以看出,在固定升力系数条件下,计算攻角普遍小于试验攻角;相同湍流模型条件下,网格密度对摩擦阻力和压差阻力均有一定影响,相比较而言,对压差阻力的影响大一些;在相同网格密度条件下,湍流模型对阻力的影响主要体现在摩擦阻力上;在固定升力系数条件下,不同的网格密度对力矩系数影响较小,不同的湍流模型对力矩系数影响较大,计算得到的低头力矩普遍大于试验结果,在以下的研究中将看到,这主要是由于网格密度对压力分布影响较小,不同湍流模型对压力分布有一定影响(主要是上表面激波位置)所致.表2 算例1计算结果与试验结果对比Table 2 Numecrical and experimental results in Case 1软件P 湍流模型攻角升力阻力压差阻力摩擦阻力俯仰力矩COAR SE (3374848)TRIP210(SA)0116701500010291010163010129-01146TRIP210(SST)0129401500010274010163010111-01138CFL3D(SA)0113901500010294010160010133-01138MEDIUM (5715968)TRIP210(SA)0116901501010298010161010137-01145TRIP210(SST)0129401500010273010159010114-01140CFL3D(SA)0111601500010288010156010132-01141FINE (9966592)TRIP210(SA)0115201500010280010155010124-01149TRIP210(SST)0129001500010271010156010114-01139CFL3D(SA)0110801500010285010154010132-01141TEST-0152001500010295---01121312 算例2计算结果及分析图3给出采用SA 模型得到的升力曲线与极曲线,同时给出相应的试验结果及采用C FL3D 软件SA 模型得到的计算结果.可以看到,不同的网格密度对升力特性基本没有影响,对阻力特性影响比较显著,最大影响量在17Counts 左右,采用SA 模型时,粗网格的极曲线与试验吻合最好,我们认为这是一个巧合.研究表明[5],模拟转捩位置与采用全湍流计算阻力相差13C ounts,考虑到C FL3D 模拟转捩位置而本文采用全湍流计算,因此,虽然CFL3D 中等网格的结果与本文密网格的结果接近,但实际上本文中等网格的计算结果与CFL3D 中等网格的计算结果是吻合的.图4给出采用SST 模型得到的升力曲线和极曲线,图中同时给出相应的试验结果和采用CFL3D 软件SST 模型的计算结果.可以看到,采用SST 模型时,不同的网格密度对升力特性和阻力特性基本没有影响,在计算的攻角范围内一直保持网格收敛性,网格密度对阻力的影响量在4Counts 以下.相同升力系数条件下,计算得到的阻力普遍小于试验的阻力系数,考虑到模型转捩位置的影响,与CFL3D 的计算结果相比较,本文采用中等网格的阻力系数偏大了13Counts 左右.图5给出采用中等网格,SA 和SST 两种湍流模型得到的压差阻力和摩擦阻力随攻角的变化曲线.由图中可以看出,在本文的计算范围内,摩擦阻力系数基本不随攻角变化;不同的湍流模型对压差阻力影响很小,147第2期王运涛等:DLR -F6翼身组合体阻力计算对摩擦阻力的影响较大,采用SA 模型的摩擦阻力比采用SST 模型得到的摩擦阻力大24Counts 左右.图3 升力和阻力极曲线(SA)Fig 13 Lift curve and drag polar(SA)图4 升力和阻力极曲线(SS T )Fig 14 L i ft curve and drag polar(SS T )图5 压差阻力和摩擦阻力随攻角的变化Fi g 15 Fraction of drag vs.angle ofattack313 压力系数与测压试验的比较法国ONERA 完成的测压试验是在升力系数C L =014984的条件下完成的,相应的攻角为0149b ,为了与相应的测压数据相比较,本节研究究竟是在相同升力下比较测压结果,还是在相同攻角下比较测压结果,及网格密度和湍流模型对压力分布的影响.图6给出在4个典型站位(2319%,3717%,5114%,8417%)计算得到的压力分布与相应试验结果的比较.计算采用中等网格和SST 两方程湍流模型,分别给出C L =015003,A =01294b 和C L =015224,A =0149b 条件下的压力分布.可以看到,在KINK 以内的站位(2319%),两种来流条件下的压力分布都与试验结果吻合较好;在KINK 以外的站位上,与试验相同攻角下的压力分布显然更接近试验结果.从计算结果看,在相同攻角下比较计算与试验的压力分布更合理一些.图6 压力分布与试验的比较Fi g 16 C p distri bu tion图7给出在4个典型站位(2319%,3717%,5114%,8417%)计算得到的压力分布与相应试验结果的比较.采用粗网格、密网格和SA 一方程湍流模型,计算来流攻角均为0149b ,相应的升力系数分别为015364和015385.由此可以看出,网格分布对压力分布的影响非常小,对阻力系数的影响比较大.上述来流状态下,粗网格下的阻力系数为010308,密网格为010297,两者相差11Counts.图8给出在4个典型站位(2319%,3717%,5114%,8417%)计算得到的压力分布与相应试验结果的比较.采用中等网格和SA 一方程、SST 两方程湍流模型,计算来流攻角均为0149b ,相应的升力系数分别为015362和015224.可以看出,湍流模型对压力分布的影响不明显,只是在翼梢的站位上,SST 模型得到激波位置略微靠前,必须指出,此时对阻力系数的影响比较大.上述来流状态下,SA 模型的阻力系数为010314,SST 模型为010283,两者相差31Counts!148计算物理第25卷图7 网格对压力分布的影响Fig 17 C p distribution with differentgrids图8 湍流模型对压力分布的影响Fig 18 C p distribution in various turbulent models图9 DLR -F6翼身组合体表面流线Fig 19 DLR -F6wing -body surfacestreamline314 局部流动细节分析图9给出采用中等网格、SA 一方程模型和SST 两方程模型得到的翼身接合部的流动细节,来流条件M =0175,A =0149b .可以看到,在上述来流条件下,翼身结合处上表面产生了小的分离区,比较两种湍流模型的计算结果,从定性方面来讲,分离区的位置、大小等方面没有明显的区别,详细研究湍流模型阻力方面的差异,应该从边界层内的速度型等方面入手.4 结论本文采用TRIP210软件,利用DPWII 提供的多块对接网格,数值模拟DLR -F6翼身组合体的外形,主要研究了网格密度、湍流模型对数值模拟结果的影响,通过与相应的试验结果和CFL3D 的计算结果相对比,得到一些基本结论:1)采用SST 湍流模型得到了网格收敛性的结果;2)相对于SA 一方程模型,SST 两方程湍流模型对网格依赖性较弱;在相同网格密度条件下,不同湍流模型对阻力的影响主要体现在摩擦阻力上;3)比较压力分布时,计算与试验建议在相同攻角下比较;网格密度对压力分布影响较小,湍流模型对压力分布有一定的影响.软件的确认工作具有相当大的难度,首先是难以得到精确的试验结果,其次是如何保证计算模型与试验模型的一致性,再次是影响计算结果的因素众多.本工作只是TRIP210软件确认工作的一个开端,在今后的工作中我们还将给出DLR -F6翼P 身P 架P 舱组合体的数值模拟结果.149第2期王运涛等:DLR -F6翼身组合体阻力计算150计算物理第25卷[参考文献][1]Johnson Forres ter T,Tinoco Ed ward T,Yu N Jong.Thirty years of developmen t and application of CFD at Boeing commercial airplaneSeattle[R].AIAA2003-3439,2003.[2]Tinoco E N,Bogue D R.Progress toward CFD for full flight envelope[J].Aeronautical Journal,2005,109(1100):451-460.[3]Oberkampf William L,Trucano Timothy G.Verification and validation i n computational fluid d ynamics[J].Progress in AerospaceSciences,2002,38:209-272.[4]Haase W,Chaput E,Leschziner M A.ECARP:European computational aerodynamics rearch project:Validation of C FD codes andassessment of turbulence models[R].Vieweg&Sohn Verlagsgesellschaft mb H,Germany,1997.[5]Levy David W,Zickuhr Tom,et al.Summary of data from the firs t AIAA CFD drag p rediction workshop[R].AIAA2002-0841,2002.[6]Laflin K R,Klausmeyer S M,Zickuhr T.Summary of data from the second AIAA C FD drag prediction workshop[R].AIAA2004-0555,2004..[7]Ru msey C,Rivers S,Morrison J.Study of CFD variation on transport configurations for the second drag-prediction workshop[R].AIAA2004-0393,2004.Drag of DLR-F6Wing-Body Con figurationWANG Yuntao,WANG Guangxue,ZHANG Yulun(Computational Aerodynamics Institute,China Aerodynamics Research and Development Center,Mian y an g621000,China)Abstract:Navier-Stokes calculation on mult-i block is performed to calculate drag of the DLR-F6wing-body configuration with CFD software TRIP210.Structured grids(1-to-1)and test results used are from Drag prediction workshopÒ(DPWII).Referenced numerical resul ts are obtained with CFL3D.Effects of mesh densi ty and turbulent models on aerodynamic characteristics and pressure distribution are carefully studied.Results are verified by experimen tal data and CFL3D results.Grid refinement leads to convergent results in SS T turbulent models.It is demonstrated that turbulent models have li ttle influence on pressure drag,bu t obvious influence on friction d rag.The turbulent model and grid densi ty have li ttle influence on pressure distribution.Key words:TRIP210;DLR-F6;d rag prediction workshopÒ;drag;turbulent modelsReceived date:2006-11-14;Revised date:2007-03-17。

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

用ＳＴ两方程模型得到网格收敛结果；不同的湍流模型对压差阻力影响较小，对摩擦阻力影响较大；不同的网Ｓ
格密度和湍流模型对压力分布影响较小．
［键词］关ＴＩ２０Ｌ —６ＤＷＩ；阻力；流模型ＲＰ．；ＤＲＦ；ＰＩ湍［献标识码］文Ａ
翼身组合体作为标准算例．上１会８家单位提供了采用１４种软件的计算结果．０３年６月召开了第二次阻力２０
计算的工作会议（ＰＩ）该次会议选择ＤＲＦ翼／挂／ＤＷ，ＩＬ．６身／舱组合体作为算例，议的重点是阻力计算精度，会共
势，已经与风洞试验一样，为飞行器气动设计最重要的研究手段ｎ．成随着大型网格前置处理软件、ＦＣＤ软件系统和后置处理软件的发展，ＦＣＤ几乎可以模拟所有高度复杂飞行器外形的绕流流场．飞行器设计者使用ＣＤ工具面临的首要问题是软件的精度和效率问题．ＦＦＣＤ软件的计算
【章编号］１０ —４Ｘ（０８０—１５０文０１２６２０）２０４－６
ＤＲ—６翼身组合体阻力计算ＬＦ
王运涛，王光学，张玉伦
（国空气动力研究与发展中心计算空气动力学研究所，四川绵阳６１０）中２００
［图分类号］Ｖ１．中２１７
Ｏ引言
随着计算机软硬件技术和计算流体力学（Ｆ）飞速发展，ＦＣＤ的ＣＤ在航空航天飞行器方案选型和初样设计阶段发挥的作用愈来愈重要，雷诺数修正、化设计、致流场刻画等方面，Ｆ在优细ＣＤ较风洞试验更具有优

DLR—F6翼身组合体跨音速流场CFD应用计算研究

问题的ＡＡ第二届工作会议（Ｐ Ⅱ）的ＩＡＤＷ中
ＤＲ＿６翼身组合体模型。首先利用 “ 立方体 ” Ｌ－Ｆ超
概念构建了绕ＤＲ６身组合体的高质量多块Ｌ —Ｆ翼
结构拼接网格，究了网格拓扑结构对ＤＲＦ研Ｌ — ６翼
动力学研究项目ＥＡＰ和致力于研究阻力精度ＣＲ…
法一直无法精确计算跨音速流场问题，一方面这是
由于跨音速流场的速度范围具有特殊性，行器表飞面跨音速区的气流是既有亚音速又有超音速的“ 混合流动 ” ，且跨音速问题对算法的耗散特性十区而分敏感，常会出现计算收敛难，靠性差等问题。通可
响具有十分重要的意义。
Hale Waihona Puke ＤＲＦＬ — ６翼身组合体是一个跨音速运输机的经典算例，国外多个风洞中进行过吹风试验，在其试验数据可靠，到了国际流体力学工作者的一致得认可。本文采用ＣＤ数值方法，用欧洲计算空气Ｆ选
摘
要
针对跨音速运输机经典算例ＤＲ６翼身组合体模型，Ｌ —Ｆ采用ＣＤ方法对其气动特性进行了黏性流动数值模拟，Ｆ流
动模型为雷诺平均ＮＳＲＮ）．（ＡＳ方程。首先采用 “ 超立方体” 概念生成绕ＤＲＦＬ — ６翼身组合体的高质量多块结构拼接网格，研究网格拓扑结构对气动特性的影响。在此基础上通过网格细分和粗分考查了网格密度对计算结果的影响，后进行了湍流最模型的影响研究。通过与实验数据对比分析，出了适宜ＤＲＦ得Ｌ — ６翼身组合体跨音速黏性流动的计算网格，并总结出了能较好模拟其跨音速流场特性的湍流模型。结果表明：网格拓扑结构的合理设计会对计算结果产生一定的影响。网格密度对机

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

ＴｒａｎｓｏｎｉｃＤｒａＣｏｍｕｔａｔｉｏｎＡｒｏｕｎｄＤＬＲＦ６ＷｉｎＢｏｄＣｏｎｆｉｕｒａｔｉｏｎｓ－ｇｐｇｙｇ
１，２１３ＺＨＥＮＧＱｉｕａ，ＬＩＵＳａｎａｎＨＯＵＴｉａｎｘｉａｏｙｙｇ，Ｚ
：ＴＡｂｓｔｒａｃｔｈｅａｃｃｕｒａｃｏｆｃｏｍｕｔｅｄｄｒａａｒｏｕｎｄｔｈｅＤＬＲ－Ｆ６ｗｉｎｂｏｄｃｏｎｆｉｕｒａｔｉｏｎａｎｄｔｈｅｙｐｇｇｙｇｅｆｆｅｃｔｓｏｆｒｉｄａｎｄｔｕｒｂｕｌｅｎｃｅｍｏｄｅｌｓｏｎｔｈｅａｅｒｏｄｎａｍｉｃｃｈａｒａｃｔｅｒｉｓｔｉｃｓａｒｅｉｎｖｅｓｔｉａｔｅｄｂｓｏｌ－ｇｙｇｙ）ｖｉｎＲｅｎｏｌｄｓ－ａｖｅｒａｅｄＮａｖｉｅｒ－Ｓｔｏｋｅｓ（ＲＡＮＳｅｕａｔｉｏｎｓｃｏｕｌｅｄｗｉｔｈｔｈｒｅｅｋｉｎｄｓｏｆｔｕｒｂｕｌｅｎｃｅｇｙｇｑｐ， ’ ，ｋｍｏｄｅｌｓ－ｔｈｅＳａｌａｒｔ－Ａｌｌｍａｒａｓ（ＳＡ）ｏｎｅｅｕａｔｉｏｎｍｏｄｅｌｔｈｅＷｉｌｃｏｘｓ－ ωｍｏｄｅｌａｎｄｔｈｅＭｅｎｔ－ｐｑ ’ ｋｅｒｓ－ ωＳＳＴｍｏｄｅｌ．Ａｍｕｌｔｉ－ｂｌｏｃｋｓｔｒｕｃｔｕｒｅｄａｔｃｈｅｄｒｉｄｒｏｖｉｄｅｄｂｔｈｅＡＩＡＡＤｒａＰｒｅｄｉｃ－ｐｇｐｙｇｔｉｏｎＷｏｒｋｓｈｏｉｓｅｍｌｏｅｄｉｎｔｈｅｃｏｍｕｔａｔｉｏｎ．Ｔｈｅｃｏｍｕｔａｔｉｏｎａｌｒｅｓｕｌｔｓｓｈｏｗｔｈａｔｔｈｅｒｅｓｓｕｒｅｐｐｙｐｐｐｄｉｓｔｒｉｂｕｔｉｏｎｓｏｎｔｈｅｗｉｎｓｕｒｆａｃｅａｒｅｅｗｅｌｌｗｉｔｈｔｈｅｅｘｅｒｉｍｅｎｔ．Ｔｈｅｖａｒｉｎｔｒｅｎｄｓｏｆｔｈｅｃｏｍ－ｇｇｐｙｇａｅｒｏｄｎａｍｉｃｆｏｒｃｅｓｗｉｔｈｔｈｅａｎｌｅｏｆａｔｔａｃｋａｒｅｃｏｎｓｉｓｔｅｎｔｗｉｔｈｔｈｅｅｘｅｒｉｍｅｎｔ．Ｇｒｉｄｒｅｆｉｎｅ－ｕｔｅｄｙｇｐｐ，ｍｅｎｔｌｅａｄｓｔｏｃｏｎｖｅｒｅｎｔｒｅｓｕｌｔｓｆｏｒＳａｌａｒｔ－Ａｌｌｍａｒａｓｍｏｄｅｌａｎｄｔｈｅｄｒａｂｅｔｔｅｒｔｈａｎｔｈｏｓｅｅｔｓｇｐｇｇｆｒｏｍｏｔｈｅｒｓｏｆｔｗａｒｅｓ．Ｔｈｅｄｉｆｆｅｒｅｎｔｔｕｒｂｕｌｅｎｃｅｍｏｄｅｌｓｅｘｅｒｔｃｅｒｔａｉｎｅｆｆｅｃｔｓｏｎｔｈｅｒｅｓｓｕｒｅｄｉｓｔｒｉ－ｐ，；；ｂｕｔｉｏｎｓｅｓｅｃｉａｌｌｏｎｔｈｅｏｓｉｔｉｏｎｓｏｆｓｈｏｃｋｗａｖｅａｌｉｔｔｌｅｅｆｆｅｃｔｓｏｎｌｉｆｔａｎｄｏｂｖｉｏｕｓｅｆｆｅｃｔｓｏｎｐｙｐ，ｄｒａｅｓｅｃｉａｌｌｆｒｉｃｔｉｏｎｄｒａ．Ｔｈｅｔｕｒｂｕｌｅｎｃｅｍｏｄｅｌｓａｌｓｏｅｘｅｒｔｅｆｆｅｃｔｓｏｎｔｈｅｆｌｏｗｓｅａｒａｔｉｏｎｇｐｙｇｐ， ’ ｋ－ｎｅａｒｔｈｅｒｏｏｔｏｆｗｉｎ．ＩｎｔｈｅｔｒａｎｓｏｎｉｃｆｌｏｗｃｏｍｕｔａｔｉｏｎｓｔｈｅＭｅｎｔｅｒｓ ωＳＳＴｍｏｄｅｌａｃｈｉｅｖｅｓｇｐｔｈｅｈｉｈｅｓｔａｃｃｕｒａｃ．ｇｙ：ｄ；Ｎ；；Ｋｅｗｏｒｄｓｒａａｖｉｅｒ－Ｓｔｏｋｅｓｅｕａｔｉｏｎｓｔｕｒｂｕｌｅｎｃｅｍｏｄｅｌｗｉｎｂｏｄｃｏｎｆｉｕｒａｔｉｏｎｓｙｇｑｇｙｇ

考虑短舱安装参数的全机气动外形优化设计

考虑短舱安装参数的全机气动外形优化设计许丹阳;白俊强;雷锐午【摘要】针对跨声速宽体客机机翼、机身、短舱、挂架、平尾、垂尾全机构型中各部件之间存在的气动干扰和力矩耦合问题,采用多块拼接自由变形(FFD)方法进行全机多部件的参数化,利用雷诺平均Navier-Stokes方程(RANS)进行流场求解,结合离散伴随技术和序列二次规划优化算法进行考虑短舱安装参数变化和平尾配平的气动外形优化设计.经过优化之后,短舱附近机翼的压力分布形态变光顺,设计点阻力减小3.82％,并且获得了无激波的压力分布形态.【期刊名称】《西安交通大学学报》【年(卷),期】2019(053)004【总页数】8页(P150-157)【关键词】气动外形优化;平尾配平;短舱安装参数;自由变形;离散伴随【作者】许丹阳;白俊强;雷锐午【作者单位】西北工业大学航空学院,710072,西安;西北工业大学航空学院,710072,西安;西北工业大学航空学院,710072,西安【正文语种】中文【中图分类】V211近年来,高性能计算机硬件平台发展突飞猛进,基于雷诺平均Navier-Stokes方程的气动优化设计方法在飞行器外形设计领域得到越来越广泛的应用[1-3]。

宽体客机作为当前民航市场上的最新一代客机,巡航马赫数达到0.85,处于典型的跨声速非线性区间。

客机机翼外形的微小改变往往会产生较大的气动力变化;同时全机多部件之间的互相干扰,也使得气动外形设计工程师无法进行精确定量的外形设计。

为了能够合理地对设计空间进行探索,以得到气动特性较好的构型,跨声速客机气动外形优化设计问题往往需要大规模设计变量[4];同时也需要考虑多部件之间的耦合来获得综合性能最优的气动外形。

国外大量的研究成果表明,基于伴随的气动优化设计方法在复杂大规模气动外形优化设计问题中有较好的应用前景。

Alonso等基于斯坦福大学非结构网格求解器(SU2)实现了旋翼以及全机构型的气动优化设计[5];Zingg等采用离散伴随方程结合Newton-Krylov隐式求解方法,开展了气动外形优化设计研究取得了显著的成果,并将该方法拓展到了自然层流翼型设计中,实现了层流翼型的优化减阻[6];Martins 等利用高可信度多学科优化设计平台(MACH)进行了机翼、翼身组合体以及风力机的优化设计研究[7]。

TRIP2.0软件的确认:DPWⅡ复杂组合体的数值模拟

TRIP2.0软件的确认：DPWⅡ复杂组合体的数值模拟王运涛;王光学;张玉伦【期刊名称】《航空学报》【年(卷),期】2008(029)001【摘要】采用"亚跨超声速计算流体力学软件平台"(TRIP2.0)数值模拟了阻力预测小组(AIAA CFD Drag Prediction Workshop Ⅱ,DPWⅡ)翼/身/架/舱复杂组合体运输机构型,数值模拟采用的多块对接网格、测压和测力的试验结果均来自DPWⅡ,对比计算采用了CFL3D的结果.重点针对DLR-F6翼/身/架/舱复杂组合体构型,详细研究了网格密度和湍流模型对总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致.采用SA一方程和SST两方程模型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布有一定的影响.【总页数】7页(P34-40)【作者】王运涛;王光学;张玉伦【作者单位】中国空气动力研究与发展中心,计算空气动力学研究所,四川,绵阳,621000;中国空气动力研究与发展中心,计算空气动力学研究所,四川,绵阳,621000;中国空气动力研究与发展中心,计算空气动力学研究所,四川,绵阳,621000【正文语种】中文【中图分类】V211.7【相关文献】1.DPW4翼/身/平尾组合体的数值模拟 [J], 王运涛;张书俊;孟德虹2.DPWⅢ机翼和翼身组合体构型数值模拟 [J], 王运涛;王光学;张玉伦3.采用TRIP2.0软件计算DLR-F6构型的阻力 [J], 王运涛;王光学;张玉伦4.TRIP2.0软件湍流模块的开发与确认 [J], 王光学;王运涛5.TRIP2.0软件非定常模块的开发与确认 [J], 孟德虹;王运涛;张玉伦;王光学因版权原因，仅展示原文概要，查看原文内容请购买。

湍流模型在民机跨声速绕流计算的应用

ＥＡＲＳＭｋ．【１）ｍｏｄｅ１．ＳｐａｌａｒｔＡｌｌｍａｒａｓｍｏｄｅｌ。ＳＳＴｍｏｄｅｌ，Ｓｔａｎｄａｒｄｋ一 ∞ ｍｏｄｅｌａｎｄＬｏｗ．ＲｅｙｎｏｌｄｓＮｕｍ—
ｂｅｒｋ一 ‘ １）ｍｏｄｅ１．ＴｈｅｒｅｓｕｌｔｓａｒｅｐｒｅｓｅｎｔｅｄｉｎｔｈｉｓｐａｐｅｒｂｙｃｏｍｐａｒｉｎｇｔｈｅｓｅｐａｒａｔｉｏｎｏｆｔｈｅｗｉｔｉｏｎｏｎＤｉｆｅｒｅｎｔＴｕｒｂｕｌｅｎｃｅＭｏｄｅｌｓｉｎＴｒａｎｓｏｎｉｃＦｌｏｗｏｆＣｉｖｉｌＡｉｒｃｒａｆｔ
ＨＡＮＴａｏ
ｃｏｍｂｉｎａｔｉｏｎ．ｏｎｃｒｕｉｓｅ，ａｎｄｉｔｉｓｈｅｌｐｆｕｌｆｏｒａｐｐｌｉｃａｔｉｏｎｏｆｔｕｒｂｕｌｅｎｃｅｍｏｄｅｌｓ．Ｋｅｙｗｏｒｄｓ：ｓｉｍｕｌａｔｉｏｎ；ＤＬＲ－Ｆ６ｍｏｄｅｌ；ｔｕｒｂｕｌｅｎｃｅｍｏｄｅｌ；ｌｏｆｗｐａｔｔｅｒｎ
韩涛
（海军驻沈阳地区航空军事代表室，辽宁沈阳１１００３４）
摘要：针对民机构型ＤＬＲ．Ｆ６翼身组合体，研究了不同湍流模型对气动特性和流场特征的影响。控制方程采用

超临界翼型跨声速抖振CFD计算和POD分析

超临界翼型跨声速抖振ＣＦＤ计算和ＰＯＤ分析
董圣华，史爱明，叶正寅，田海涛
（西北工业大学航空学院，陕西西安７１００７２）
摘
要：使用二阶迎风Ｒｏｅ格式、隐式时间推进和孓Ａ（Ｓｐａｌａｒｔ－Ａｌｌｍａｒａｓ）一方程湍流模型，通过求解基于格心格式
ｓｏｌｕｔｉｏｎｏｆｔｈｅｃｅｌｌ－ｃｅｎｔｅｒｅｄｆｉｎｉｔｅ－ｖｏｌｕｍｅｍｅｔｈｏｄ（ＦＶＭ）ｂａｓｅｄＲＡＮＳ（Ｒｅｙｎｏｌｄｓ－ａｖｅｒａｇｅｄＮａｖｉｅｒ — Ｓｔｏｋｅｓ）ｅｑｕａｔｉｏｎｓｗｉｔｈｔｈｅＳ — Ａ（Ｓｐａｌａｒｔ — Ａｌｌｍａｒａｓ）ｏｎｅ — ｅｑｕａｔｉｏｎｔｕｒｂｕｌｅｎｃｅｍｏｄｅｌ，ａｓｗｅｌｌａｓｔｈｅｉｍｐｌｉｃｉｔｔｉｍｅ — ｓｔｅｐｐｉｎｇｓｃｈｅｍｅａｎｄｔｈｅＲｏｅｓｃｈｅｍｅ．Ｏｎｔｈｅｂａｓｉｓｏｆｓｕｃｃｅｓｓｆｕｌｓｉｍｕｌａ —
ｅｖｏｌｕｔｉｏｎｏｆｆｌｏｗｓｅｐａｒａｔｉｏｎａｎｄｉｔｒｅｖｅａｌｓａｈｙｓｔｅｒｅｓｉｓｐｈｅｎｏｍｅｎｏｎｔｈａｔｔｈｅｐｈａｓｅｏｆｓｈｏｃｋｍｏｖｅ —

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

基于混合网格的DLR-F6模型跨声速阻力预测粘流计算周磊;杨振虎;邱征;白文
【期刊名称】《航空科学技术》
【年(卷),期】2005(000)006
【摘要】该文详细介绍了针对第二届AIAA阻力预测会议公布的标准模型及其网格进行的基于混合网格的纳维尔-斯托克斯粘流计算考题,我们将自主开发的WoF90程序的计算结果与DLR的TAU程序以及FOI的EDGE程序的计算结果进行了对比分析,发现三种程序的升力系数普遍比实验结果偏高,但偏高的程度大致相同,阻力计算的精度范围大体一致.结果确认了300万网格点的计算结果尚不足以得到网格收敛的效果,阻力预测需要更大规模的网格.
【总页数】4页(P37-40)
【作者】周磊;杨振虎;邱征;白文
【作者单位】中国航空计算技术研究所;中国航空计算技术研究所;中国航空计算技术研究所;中国航空计算技术研究所
【正文语种】中文
【相关文献】
1.基于非结构平台的DLR－F6标模阻力预测 [J], 章锦威;戚姝妮;郭承鹏;董军
2.基于BP神经网络的模型铲刀推土阻力预测模型 [J], 黄葵;李因武;李建桥;彭剑波
3.基于流通阻力的DPF碳载量预测模型研究 [J], 石秀勇;蒋得刚;梁云芳;梁鹏飞
4.基于巷道摩擦阻力系数BP神经网络预测模型的矿井风网风量预测研究 [J], 刘
彦青
5.基于神经网络的杨树苗木注射针头阻力预测模型研究 [J], 孙炜;陆一枝;商庆清因版权原因，仅展示原文概要，查看原文内容请购买。

机翼—机身—吊架—短舱的跨声速流计算

机翼—机身—吊架—短舱的跨声速流计算
孙刚;乔志德
【期刊名称】《空气动力学学报》
【年(卷),期】1995(013)004
【摘要】本文通过使用多块网格的嵌合全技术，计算了机翼－机身－吊架－短舱
的跨声速流场。

将机翼－机身－吊架－短舱流场分为机翼机－机身和吊架－短舱两个子区域，根据嵌套网格结构建立机翼－机身－吊架－短舱的组合网格，分别对这两个子区域根据相应的边界条件独立求解，然后采用嵌合体技术传递两个子区域的干扰信息，通过两个流场的耦合迭代得到机翼－机身－吊架－短舱流场的计算结果，数值模拟结果与实验结果进行了比较，两者吻合较好。

【总页数】8页(P374-381)
【作者】孙刚;乔志德
【作者单位】不详;不详
【正文语种】中文
【中图分类】V211.3
【相关文献】
1.结合部网格系统在机身-边条-机翼组合体绕流计算中的应用 [J], 李孝伟;范绪箕
2.机翼机身跨声速绕流的计算 [J], 孙刚;乔志德
3.绕侧滑锥体,三维机翼与机身的超音速非线性全位势流计算 [J], 陈红全;黄明恪
4.大展弦比直机翼双机身飞机中机身与机翼的最佳布置问题 [J], 邓扬晨;吴敬涛;贺
集乐
5.机翼-机身-挂架-短舱复杂流场数值模拟 [J], 高正红;黄江涛
因版权原因，仅展示原文概要，查看原文内容请购买。