螺旋桨流场Euler/NS数值解对比分析

采用多块点对点网格生成技术,生成螺旋桨流场空间计算网格。针对不同前进比,分别求解Euler和NS方程对螺旋桨流场进行数值模拟,并与实验结果进行了对比分析。研究发现无粘和有粘的计算结果均接近于实验结果,因此基于Euler方程的流场数值模拟更适用于工程设计。     

采用近似边界求解壁板颤振的CBS有限元算法

为了简化Euler方程求解壁板颤振问题,将近似边界条件成功地推广到CBS(characteristic-based split)有限元方法求解壁板颤振领域,把运动壁面的影响转化为Euler方程的边界条件,避免了对动网格的处理.通过对近似边界方法、动网格方法以及3阶活塞理论的数值模拟结果进行对比发现:首先,加载近似边界的CBS有限元方法操作非常简便,且模拟结果与动网格方法完全吻合;其次,Euler方程可以准确捕捉到壁板前后两端不光滑尖点所引起的压力突跳,而3阶活塞理论则无法做到;最后,采用近似边界的CBS有限元方法在求解无黏可压缩流动中细长体小变形的流-固耦合问题上具有潜力和优势.      

基于Euler方程的二维高速海面效应数值分析

以近海面飞行器的海面效应为研究背景,通过数值求解可压缩Euler方程模拟了波形壁对翼型高亚声速气动力特性的影响.计算网格布局采用了重叠网格方法,避免了翼型与波形壁之间网格生成的困难.详细分析了波形壁边界条件,给出了正确的通量计算公式.时间推进方法采用了双时间方法,子迭代过程由LU-SGS方法完成.通过数值分析,给出了飞行高度、波长和波幅等因素对翼型气动力特性的影响.     

基于控制理论的Euler方程翼型减阻优化设计

在Jameson的控制论气动优化设计思想下，应用Euler方程研究了翼型的反设计和减阻问题。设计过程中的梯度是通过求解一个伴随偏微分方程而得到的，每个设计迭代只需求解一次Euler方程和一次伴随方,一与设计变量数无关。在具体数值执行中：①利用分部积分将目标函数变分的计算转化成了类似于计算Euler方程残值的形式，节省了机器时间；②根据控制理论的要求，将减阻问题转化成了相应地修改伴随方程的物面边界条件和目标函数的变分表达，使问题得到了简化；③利用特征不变量分析法，处理了伴承方程的物面和远场边界条件。设计算例证明了本文方法可靠性好、收敛快、大大节约设计时间。     

Euler方程的自适应多重非结构网格计算

介绍一种基于四面体网格的多重网格技术,用于求解三维Euler方程定常解,具有十分明显的加速收敛效果。求解Euler方程的基本算法采用有限体积中心差分格点格式。多重网格技术采用一组几何上相互独立的粗细不同的网格。流动变量、余量和修正量在相邻粗细网格间的相互传递通过线性插值实现。采用自适应技术对网格进行加密,可以获得较优的网格,并作为一套细网格,进行多重网格计算。数值实验表明本文自适应多重网格法十分有效。     

Euler及N-S方程的二维四边形网格自适应算法

描述了一种基于贴体网格的Euler和N-S方程自适应算法。采用中心有限体积法,对Euler及N-S方程进行数值求解,N-S方程的计算中加入了基于SA湍流模型的DES模型。在流场中,运用了网格自适应算法,提高了数值计算对激波、流动分离、涡等特性的捕捉和分辨能力。采用上述方法,数值分析了单段翼型的绕流问题。文末给出的2个二维问题的对比算例表明本文方法有效。     

隐式无网格算法及其应用研究

本文的主要目的在于研究求解Euler方程的隐式无网格算法，并应用于复杂的流场计算。采用无网格算法，计算区域用点云离散代替通常的网格划分；计算点上的空间导数，用当地点云上引入的二次极小曲面逼近。求解的Euler方程用隐式时间后差离散，结合用Roe的近似Riemann解确定通量，并用LU－SGS算法分步计算，数值求解了单翼型或双翼型模拟的复杂绕流。     

多体分离数值模拟及方法验证

利用自主研发的混合网格Euler／RANs方程解算软件WoF90以及六自由度动力学与运动学方程解算程序,采用脚本语言技术实现网格的自动生成,建立了一套可用于多体分离问题的准定常Euler流动非结构网格方法。通过AEDC算例对非结构准定常Euler方程方法解决多体分离问题进行了计算研究,验证了方法的可行性。     

Euler方程与边界层积分方程耦合求解跨音速翼型绕流

 应用Euler方程求解跨音速翼型特性时考虑了粘性影响,粘性影响是通过边界层动量和能量积分方程求解的,即粘流/无粘流迭代方法。其中Euler方程采用LU-ADI方法求解;边界层方程均由正解法过渡到反解法,以解决强激波干扰区出现小分离泡的计算问题。计算中使用了贴体C网格,通过一定变换使其保持基本正交。计算结果表明,压力分布、摩阻系数分布与实验结果符合较好。     

三维Euler方程的自适应非结构网络计算

介绍一种三维自适应非结构网格的Euler方程求解技术。采用有限体积中心差分的格点格式对Euler方程进行空间离散,多步龙格-库塔时间推进,并采用当地时间步长、焓阻尼修正等加速收敛。采用Delaunay三角剖分技术生成三维流场的四面体网格。网格自适应技术中的网格局部加密也采用Delaunay三角剖分技术。最后用ONERAM6机翼的跨音速绕流计算说明本文方法的有效性。     

叶栅绕流计算的高分辨率隐式解法

以一种性能优越的高分辨率格式为基础, 发展了其隐式解法, 提出并实现了与这种格式相适应的隐式边界处理过程, 从而显着提高了高分辨率格式的收敛速度。用不同来流条件下叶栅绕流数组计算结果与实验做了对比, 以说明解的准确程度和本文算法的效果。      

旋成体轴对称跨声速全位势流的高效差分算法

本文用守恒型全位势方程,贴体坐标网格,对旋成体轴对称跨声速绕流的差分数值计算方法进行了研究;根据最佳收敛准则,提出了轴对称情形的AF2迭代算法,并将此算法应用于半球头柱体、弹体等各种外形的旋成体。与一般方法仅适用于亚声速自由流不同,本计算可从亚声速、跨声速自由流一直到低超声速自由流。计算结果表明,本文方法收敛快,与实验及其它方法的结果符合较好。     

三维机翼可压缩层湍流边界层计算

 从三维时间相依可压缩边界层动量积分方程和平均流动能积分方程出发,将二维可压缩层、湍流边界层积分方程算法推广到求解有限翼展后掠机翼的边界层流动。利用四步Runge-Kutta时间步进方案数值求解积分方程组,并利用当地时间步长加速迭代收敛。文中分析了数值方法的稳定性与收敛性,并考查了横向耗散项对计算结果的影响。算例表明,能获得令人满意的三维机翼定常可压缩层、湍流边界层的计算结果。     

多重网格法求二元跨声速流在一种半自适应网格上的欧拉解

本文采用MacCormack两步显式格式,用有限体积法求解了二元跨声速欧拉流。推导了物面边界条件,采用了特征远场边界条件及远场环量修正。利用保角变换方法生成O型贴体网格,并修改得到了一种在激波处局部加密的半自适应贴体网格。采用多重网格及焓阻尼加速收敛技术计算了NACA0012翼型的跨声速气动特性,得到了十分满意的收敛过程和计算结果。     

一维可压缩流的有限元状态空间模型

通过引入位置交错的两种有限控制体积单元——状态单元和速度单元,对其分别应用质量、能量守恒方程和动量守恒方程, 建立了一维可压缩流的有限元状态空间模型, 模型单元边界条件的处理方式简单有效, 适用于处理任意的物理边界条件。与分析解的比较表明该模型是精确的     

亚跨超音速喷管内流场的整体矢通量分裂算法

从N－S方程出发，采用先进的矢通量分裂算法，对亚跨超音速喷管内流场进行了整体模拟计算。文中对几种不同计算格式的计算效率进行了比较，给出了超音速占优混合流动的直接分割求解方法。计算采用代数湍流模型，跨音速段的计算结果同其它实验和计算结果进行了比较。计算中发现对出口边界条件进行亚音速、超音速和有回流分别处理很有必要；边界层内的M数沿壁面法向变化非常快；不同壁温条件和不同的流动Re数均会改变边界层的厚度。     

涡轮叶栅叶型损失的数值模拟方法

给出了一个计算亚、跨音涡轮叶栅叶型损失的数值计算方法。主流采用时间推进有限体积法求解积分型欧拉方程,并采用了局部网格修正技术;附面层采用全隐格式求解有限差分方程;在叶栅出口与远后方均匀流之间进行了叶片尾迹与主流的掺混损失计算。算例表明本文的数值方法可准确地预测涡轮叶栅的叶型损失。     

固体火箭喷管两相粘性跨音速流场计算

本文进行了固体火箭喷管两相粘性湍流跨音速流场计算.粘性的气相控制方程用隐式近似因子分解法求解,粒子方程采用跟踪粒子轨迹的特征线法求解,粘性湍流选用代数模型,气相和凝相充分地偶合.CFL数可以取500左右,收敛速度快,使粘性两相跨音速喷管流场计算耗费的机时达到工程计算可以接受的程度.通过计算,获得了在粘性和粒子同时作用下的流场参数分布.这对固体火箭发动机流场需要同时考虑粒子和粘性作用的专题研究大有裨益.     

用多层网格法解二维无粘跨音叶栅流动

本文给出一种较有效的多层网格构成形式,并提出“预测—修正”壁面边界处理法。通过采用当地时间步和引入余差光滑处理法,进一步提高了格式的收敛性和稳定性。作为算例,本文对二维无粘跨音速叶栅流动求解,计算结果与试验数据对比表明,本算法收敛快和计算精度好。     

多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。