应用等离子体进行翼型前缘分离控制的数值模拟研究

将等离子体对中性气体的作用模型化为彻体力矢量,求解带源项的Navier-Stokes方程,数值模拟了等离子体激励器在NACA0015翼型大迎角下的分离控制效果,彻体力为净电荷在外加电场作用下的电场力.解拉普拉斯方程得到外加电场分布,等离子体中的净电荷分布由泊松方程给出.为了较好地模拟分离涡的发展,采用了雷诺平均与大涡模拟相结合的脱体涡模拟(Detached Eddy Simulation)方法.通过与实验结果对比,发现该模型能较好地模拟等离子体激励器的控制作用.     

单级跨声速压气机进气畸变数值模拟

基于彻体力模型的基本思想,发展了一套能够反映进气畸变对风扇/压气机气动性能影响的三维数值计算程序CSAC.首先详细研究了彻体力模型的具体建模方法,然后发展了相应的三维数值计算程序.利用该程序针对某单级跨声速压气机在均匀进气下的流场进行了计算验证,计算结果与雷诺平均Navier-Stokes(RANS)计算吻合得很好,流场主要参数分布与相关实验数据也有很好的一致性.利用该程序分别计算分析了进口存在稳态周向总压畸变及周向总温畸变的全环三维畸变流场,结果表明:压气机进口畸变流场存在大尺度、强三维特征,将导致压气机总压比及稳定裕度下降.所有计算结果均显示该模型能够在降低对计算资源和工程经验依赖的同时,有效分析进气畸变对压气机特性的影响.      

跨声速轴流压气机特性全三维黏性流动分析

应用Steger-Warming通量分裂技术,将守恒型方程中的流通向量分裂成两部分,应用NND(non-oscillatory,containing no free parameters and dissipative)格式对无黏项进行数值离散,黏性项采用中心差分格式,采用LU-SGS(lower-upper symmetric Gauss-Seidel)隐式推进迭代建立了跨声速轴流压气机性能及内流场三维Navier-Stokes方程的高分辨率和高效率的数值方法.数值模拟给出NASA rotor 37和NASArotor 67性能及内流场流动结构.并进一步基于混合平面法建立了多级轴流压气机性能计算模型,详细分析了某一级跨声速轴流压气机性能和内流特性,计算与实验结果比较表明了本模型与方法是可行的和可靠的.      

应用等离子体进行翼型升力控制的数值模拟研究

将等离子体对中性气体的作用模型化为彻体力矢量,求解带源项的Navier-Stokes方程,数值模拟了在NLF(1)-0213翼型上表面60%弦长处安装等离子体激励器对升力的控制效果.彻体力为净电荷在外加电场作用下的电场力.解拉普拉斯方程得到外加电场分布,等离子体中的净电荷分布由泊松方程给出.升力线计算结果与实验值吻合,激励器工作时,升力线向上平移,控制效果与襟翼类似.     

基于IPDG方法的超声速混合层流动数值模拟

为了满足超声速混合层高精度模拟需求,实现了基于内罚方法的间断伽辽金(IPDG)方法数值模拟.通过将黏性通量作为辅助变量使得Navier-Stokes方程降阶,并利用间断伽辽金方法进行空间离散,最后采用Newton-Krylov隐式方法对空间半离散方程进行时间推进.相对于有限体积法数值精度提高到了 3阶.将该方法应用于对...     

基于混合网格的三维Navier-Stokes方程并行算法

提出了一种基于混合网格的三维Navier-Stokes方程的并行计算方法。Navier-Stokes的求解采用了基于面的有限体积方法,该方法适用于任何网格类型。采用一方程Spalart-Allmaras模型来计算紊流黏性。并行计算采用区域分裂的方法,利用METIS网格分区系统实现了各节点的加载平衡。节点间的数据交换通过调用MPI库函数来实现,采用非阻断通讯的方式来减少数据交换时间。充分利用FORTRAN90的动态存储特性来减少对内存的需求。最后,通过对绕DLR-F6外形(翼身组合体+挂架+发动机短舱)黏性流动的数值模拟,验证了该并行程序的准确性,高性能并行计算以及处理复杂几何外形的能力。     

一种叶轮机三维叶型黏性反问题方法

提出了一种应用于叶轮机三维叶片气动设计的黏性反问题方法,该方法将叶片表面边界视为是运动的壁面,通过建立并求解壁面法向动量方程,从而可根据给定的压力分布或压力面与吸力面之间的压力差分布直接求解叶型几何.结果表明:该方法可与定常雷诺平均Navier-Stokes(RANS)模拟配合进行求解,鲁棒性强、计算耗时短,且符合真实黏性流动环境,无需对叶型进行黏性修正.算例验证结果显示,该方法在单排和多排叶片环境下都能获得反问题解,计算稳定性和收敛性良好.      

一种改进的类DES湍流模拟方法

构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格子湍流涡黏性输运方程,并采用混合函数将其与S-A模型方程进行混合,从而改进了DES(detached eddy simulation)模型的亚格子行为,同时克服了其依靠网格控制模型转换的缺点.模拟了超声速的带斜坡凹腔流动,并与相同网格下的LES及DES结果进行了比较,结果表明该混合RANS/LES方法在远离壁面的自由剪切流区域与LES行为一致,而在附着边界层区域表现优于LES和DES方法.      

三维动态非结构重叠网格Navier-Stokes方程并行算法

 提出了一种三维动态非结构重叠网格Navier-Stokes(N-S)方程的并行计算方法。N-S方程的空间离散采用格点有限体积方法,时间离散采用隐式的双时间步长方法。应用一方程Spalart-Allmaras(S-A)模型来计算湍流黏性。并行计算采用动态的区域分裂方法,在每一物理时间步利用METIS网格分区系统对网格进行分区。为了实现各CPU之间的负载均衡,每块网格都按CPU个数进行分区并对活动节点和非活动节点进行了加权处理。最后,通过对外挂物投放无黏流动的数值模拟和内埋武器弹舱开启黏性流动的数值模拟,验证了该并行程序的准确性、高性能并行计算以及处理复杂几何外形的能力。     

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

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

风扇进气畸变三维非定常数值模拟技术研究

采用三维欧拉方程加体积力模型的方法, 发展了一套三维非定常进气畸变数值模拟程序, 以实现风扇/压气机在进气畸变条件下的全环面三维非定常数值模拟.利用该程序模拟计算了某航空发动机风扇进口稳态总压畸变条件下特性和流场结构, 并与试验结果进行了对比, 结果表明该程序用较小的计算机资源即可以实现对风扇全环面三维非定场畸变流场的模拟, 能够以较高的精度预估气流畸变对压气机/风扇特性的影响.      

基于CFD方法对微型扑翼翼型设计的研究

基于CFD方法从扑翼设计的角度出发对扑翼的外形进行深入的研究.使用基于广义无限插值理论的代数网格生成方法开发适合扑翼运动的三维结构N-S动网格,并应用有限体积法结合双时间推进技术求解预处理后的低速三维非定常N-S方程,完成了对低速扑翼绕流的数值模拟.通过与文献算例的比较验证了本文方法的可靠性和正确性.在此基础上,对文献[8]中介绍的一种混合翼型的相关参数变化对扑翼气动特性的影响开展了深入的研究.     

某发动机风扇进气畸变数值模拟

采用周向平均体积力的方法将叶片对气流的作用简化为体积力源项研究进气畸变对某航空发动机风扇的影响.分别对该风扇在均匀进气条件下和进口存在稳态总压畸变条件下的流场进行了模拟计算,得到了和雷诺平均Navier-Stockes(RANS)计算几乎完全一致的特性线以及流场主要结构特征,在进行进气畸变条件下流场模拟计算时成功地捕捉到了进口畸变扰动在流场中的传播情况,并能够正确反映进气畸变对压气机的影响.结果表明:压气机进气畸变流场有着大尺度、强耦合、信息传播三维性明显等特征,采用周向平均体积力研究方法可以利用非常少的计算资源加深对进气畸变流动的认识.      

跨声速轴流风扇转子流场三维粘性数值分析

采用一种快速求解三维粘性流场的计算方法求解跨声速压气机风扇转子内部流场。该方法以ＬＵ型隐式格式和ＭＵＳＣＬ ＴＶＤ迎风格式为基础,结合壁面函数方法和简单的混合长度湍流模型,使用多重网格迭代加速收敛技术,对三维可压缩雷诺平均Ｎ－Ｓ方程进行求解。计算得到了ＮＡＳＡ Ｌｅｗｉｓ ６７跨声速、低展弦比轴流风扇转子的性能曲线,并重点分析了最高效率点工况和近失速工况下的内部流场。计算与实验结果的对比表明,此方     

双级风扇三维粘性非定常流动的数值解

针对一台双级跨音速轴流风扇动静叶片排相互干扰的三维非定常流动,进行了N-S方程的数值求解。首先计算风扇定常流场,以定常流收敛结果为起始场计算非定常流动。计算结果表明,非定常流动对各叶片排进口气流角有较大的影响,其中静叶的变化大于动叶,与定常流结果相比,非定常的瞬态攻角增大的量大于减少的量,最多时增大近10°左右。静子叶片受到的非定常气动力变化幅度也大于转子叶片,叶片排所处的轴向位置不同,非定常流动的影响会有较大的差别。      

基于混合笛卡儿网格方法的可压流动数值模拟

以可压缩黏性流动的数值模拟为研究背景,发展了一套自适应混合笛卡儿网格(AHCG)方法以及基于有限体积方法的雷诺数平均Navier-Stokes(RANS)的数值求解方法.为更好地模拟边界层的黏性流动在近壁面处采用贴体结构网格,剩余计算区域自动生成与之相重叠的笛卡儿网格,并同时发展了基于流场特征的笛卡儿网格自适应技术.结合ADT(alternating digital tree)算法显著减少了网格生成中“挖洞”和“贡献单元”搜索的消耗机时,50万左右的网格数目下,搜索耗时为0.062s,仅为普通遍历方法的1/1847.通过二维圆柱与两段翼型绕流的数值算例显示,定常AHCG方法能够准确地预测物面压力分布与升阻力系数并且具备处理复杂外形的能力;同时通过二维非定常圆柱绕流问题与NACA0015矩形机翼翼尖尾涡的捕捉算例显示,结合了动态自适应网格加密的非定常AHCG方法尤其适用于旋涡主导流动.      

用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

绕大攻角翼型低速非定常流动的SUPG有限元数值模拟

采用既有良好稳定性又有较高精度的ＳＵＰＧ有限元方法,通过数值求解涡量一流函数形式的不可压缩流ＮＳ方程组,模拟大攻角翼型在高Ｒｅ数下几种典型情况的非定常绕流流场,同时验证了ＳＵＰＧ有限元方法在求解对流占优的二维非定常对流的－扩散方程时的有铲性。     

用粘性体积力方法计算高速离心压气机叶轮内部流场

高速离心压气机是用途非常广泛的一种动力机械,充分了解其内部复杂的三维流场对设计高压比、高效率离心压气机具有很重大的意义。本文采用一种简化的Navier-Stokes方程的计算方法-粘性体积力方法,计算了Krain叶轮内部流场。使用H型网格离散微分方程,并采用多重网格方法和当地时间步长方法加速计算收敛。粘性应力通过Baldwin-Lomax模型计算。为了证明计算结果的有效性,文中也给出了相应的激光测量结果。      

低速压气机静叶非定常Clocking效应数值模拟

采用三维粘性非定常数值模拟方法研究了重复级低速压气机内部非定常流动,其中包括静叶Clocking(时序)效应以及尾迹-边界层、尾迹-尾迹和尾迹-泄漏流的相互作用.对静叶时序效应的结果分析显示,当两级静叶相对周向位置改变时最大效率变化只有0.1%,表明低速压气机中静叶时序效应不明显.另外,非定常结果流动图画显示,在非定常环境下上游周期性尾迹的通过对下游叶片边界层、尾迹以及叶尖泄漏流的发展都有一定影响.