二维Euler方程的自适应非结构直角网格算法

本文用四叉树型非结构直角网格和有限体积法求解了二维Ｅｕｌｅｒ方程,在基本网格生成过程中,采用了以几何外形为基础的网格加细方法,包括物面边界网格加细和曲率加细;并提出了一个用自适应的方法在激波处及其附近局部加细基本网格过程,提高了激波分辨率,在计算中将结构网格经常采用的以中心差分为基础的Ｊａｍｅｓｏｎ的有限体积法做了推广,减少了计算工作量,节省了机时,为今后推广到三维创造条件,本文对单个翼型,多段翼     

二维Euler方程的自适应非结构网格生成方法

发展了一种用于数值模拟二维Euler方程的自适应非结构网格的生成方法。这一方法根据背景网格提供的信息直接确定内部结点的位置，然后依照Delaunay准则，将新结点加入到已有网格划分中。用于平观流及管道流的实际计算中，结果表明，该方法不受流场形状及流动特征的限制，具有较好的通用性。     

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

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

二维Euler方程的非结构网格再生成自适应方法

给出了提高二维Ｅｕｌｅｒ方程定常解质量的非结构网格自适应方法和模拟结果。流场算法为结点型有限体积法和Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式，网格用阵面推进法生成，生成初始网格时使用均匀的直角坐标背景网格。     

二维DRAGON网格自动生成技术的研究

对重叠网格技术中新出现的Cut-Paste方法进行了研究和改进,结合网格洞边界标识的缝合与重组织技术,首次将其应用在DRAGON网格生成中,提出了二维DRAGON网格的自动生成方法,最后编制出通用程序,给出了多体网格的复杂示例.实际计算表明,本文方法具有自动化程度高、可靠性好、DRAGON网格质量高等特点.     

基于非结构重叠网格的二维N-S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

跨声速机翼操纵面非定常气动力Euler方程计算

采用代数法与椭圆型方程优化相结合的结构生成方法为带操纵面F5机翼生成面对接块结构网格;弹性变形技术生成动态网格;有限体积法、双时间推进法求解非定常Euler方程.给出一种适于面对接块结构网格计算的守恒量插值方法.数值计算了F5机翼操纵面作简谐运动时的非定常气动力,计算值与试验值吻合较好.     

二维可压缩Euler方程组的SUPG有限元数值解

本文在文献[1]的基础上,构造了二维非定常可压缩Euler方程组的SUPG变分方程组。文中对由文献[1]提出的预估——校正算法作了分析,并用该算法对变分方程组进行解算,所做算例分别为:无粘性激波在固壁上的反射、零攻角10％抛物翼型亚、跨声速绕流、带攻角NACA 0012翼型跨声速绕流,取得了较好的数值结果。     

一种计算非定常二维流动的无网格算法

主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。     

基于非结构重叠网格的二维外挂物投放模拟

本文发展了基于非结构重叠网格的动态网格技术,并提出了新的相邻单元搜索算法。非结构重叠网格集非结构网格和重叠网格的优点于一身,适合处理复杂外形和多体相对运动问题。新的相邻单元搜索算法在碰到物面时能继续搜索下去,搜索效率大大提高。非定常Euler方程的求解,采用有限体积法,时间上采用双时间步长推进。为验证算法的有效性,先对绕振荡NACA0012翼形的流动进行了模拟;然后通过耦合求解刚体三自由度运动方程,模拟了二维外挂物投放过程,结果表明非结构重叠网格在处理多体相对问题时具有优越性。     

预处理法求解定常/非定常混合网格的全速流场

运用时间导数预处理法在混合网格上求解三维定常/非定常N-S方程。预处理后的N-S方程,用二阶迎风格式的有限体积法离散,隐式、LU-SGS迭代求解,双时间步推进求解非定常流。本文以一维预处理欧拉方程为基础进行分析,提出了一种特征变量的边界条件以适应预处理后的特征系统。从极低速到跨声速范围的众多算例表明,本文的边界条件处理稳定有效,预处理法在各种速度下都具备同样快速的收敛性,在较宽的使用范围内均能得到理想的计算结果。     

二维高亚声速空腔流激振荡的数值模拟研究

采用基于P .L .Roe的近似Riemann解的修正Osher&Chakravarthy (MOC)三阶TVD有限差分格式 ,数值求解二维雷诺平均全Navier Stokes方程 ,并用Cebeci Smith代数湍流模型 (对腔内区域作修正 )来模拟湍流效应 ,时间方向的积分采用四阶Runge Kutta方法 ,通过对GAMM超音速前台阶绕流的计算 ,验证了格式及程序的有效性 ,对高亚声速来流下的空腔流动作了数值模拟研究 ,取得了较好的效果     

旋翼/机身干扰非定常流场数值模拟(英文)

基于非结构动态嵌套网格方法对直升机旋翼/机身干扰非定常流场进行了三维Euler方程数值模拟。嵌套网格系统由包含旋翼的运动区网格和包含机身的静止区网格组成,两区网格均来自于同一套初始的整体非结构网格,人工边界以及初始的插值贡献单元对应关系在网格分区时同时确定。在非定常计算中两区网格始终存在一定的重叠区域,无需重新建立人工边界。高效的贡献单元搜索算法可以根据初始的插值贡献单元对应关系很快更新每一个真实时间步上的贡献单元对应关系。为了准确计入旋翼桨叶的挥舞和变距变形,在运动区网格中引入线性弹簧方法。采用非定常Euler方程和双时间推进方法模拟了美国佐治亚理工学院旋翼/机身干扰实验外形非定常流场,计算结果与实验值以及文献的计算值相吻合。     

地效飞行器近波浪地面大迎角飞行分离流动数值研究

数值模拟了地效飞行器近波浪地面大迎角飞行的分离流场,研究飞行高度和地面波浪对气动性能的影响.使用Fluent软件求解非定常可压缩流动的质量加权平均NS方程和标准k-ε湍流模型,使用滑移网格处理地效飞行器和波浪地面间的相对运动.飞行高度降低,机翼上翼面的附着流动转化为分离流动.波浪地面飞行气动力平均值和水平地面飞行气动力随飞行高度变化的规律完全一致.     

螺旋桨非定常滑流数值模拟

运用非结构动态嵌套网格技术,通过求解三维非定常Euler方程,对螺旋桨飞机的螺旋桨滑流进行了真实数值模拟.动态嵌套网格方法通过建立运动区子域网格来模拟螺旋桨的旋转运动,静止的机身包含在静止区网格中,有效解决了螺旋桨滑流模拟中的相对运动问题.采用双时间方法求解控制方程,在每一个物理时间步上运用经典四步Runge-Kutt...     

Computations of unsteady transonic flow over airfoil at low Reynolds number

Transonic flow over a thin airfoil at low Reynolds number was studied numerically by directly solving two-dimensional full Navier-Stokes equations through 5th order weighted essentially non-oscillatory(WENO) scheme without using any turbulence model.A series of distinguished unsteady phenomena for a thin 2-D transonic airfoil flow were presented.Due to continuous adverse pressure gradient in the subsonic flow downstream of the sonic line, the unsteady separated boundary layer with main vortex and secondary vortex was developed at the rear of the airfoil.At the trailing edge,the vortex-shedding was characterized by periodical connection of the main vortex and secondary vortex on the other side of the airfoil.The unsteady separation and vortex-shedding occurred with the same period.On the airfoil surface,the average pulse pressure related to the unsteady supersonic region was obviously smaller than that related to the vortex-shedding at the trailing edge.With the attack angle increasing from 0° to 2°, the frequency of vortex-shedding decreases about 4.2%.At last, the turbulence intensity and many second-order statistics in the wake region were investigated.      

透平机械空穴漏气非定常流动的数值模拟

为了解空穴漏气在透平中非定常流动,对某型带空穴漏气单级透平进行非定常数值模拟,计算格式采用包括LU-SGS-GE隐式格式和改良型高精度、高分辨率的MUSCL TVD格式的时间推进算法,非定常计算采用Jameson双时间步隐式迭代法.计算结果表明空穴漏气随着二次流的影响沿着动叶吸力面一侧向下游输运,输运过程中会不断沿径向向外迁移.由于空穴漏气的射入,非定常效应在动叶叶根处要比中间截面处明显,而且随着漏气向下游的迁移,非定常效应也逐渐削弱.      

多段翼型局部主动变形流动控制的非定常数值模拟

对30P30N三段翼型失速攻角附近的分离流动进行了数值模拟研究。为了抑制大攻角时背风区的流动分离,在主翼段上表面引入了行波壁变形模型和抛物型局部主动振动模型,利用作者以往发展的动态混合网格技术和相应的非定常计算方法,对变形过程中的非定常分离流动进行了数值模拟,分析了各种变形参数对流动分离的影响。计算结果表明,在适当的条件下,局部主动变形能够抑制翼型背风区的分离,由此可以起到增升减阻的作用,改善翼型的气动性能。     

低亚声速和跨声速矩形柱绕流的大涡模拟研究

采用基于非结构网格的有限体/有限元混合格式和大涡模拟的方法求解可压缩Navier-Stokes方程,研究了不同长宽比矩形柱低亚声速和跨声速绕流的流动特性。在雷诺数为22000时,对来流马赫数等于0.1和0.75,截面长宽比分别为1∶1、2∶1、3∶1和4∶1的矩形柱绕流进行了大涡模拟,以研究长宽比和压缩性对矩形柱绕流流场的影响。马赫数为0.1时,Strouhal数随着长宽比的增大先降低再增大然后再降低;长宽比为3∶1和4∶1时会有流动的再附产生;柱体上表面的三维特性在长宽比大时更明显。马赫数为0.75时,Strouhal数随着长宽比的增大逐渐减小;湍流脉动和涡脱落受到抑制;方柱的近尾迹区域,有两种形成机制不同的局部超声速区。     

振荡叶栅三维非定常流动数值分析

以全三维、非线性、非定常欧拉方程计算作用于振荡叶栅的非定常气动载荷。欧拉方程的数值离散采用了具有时、空二阶精度的MacCormack显式分裂的两步格式,用时间推进方法得到其数值解。应用动网格技术捕获叶片的运动,保证叶片振动过程中网格边界与叶片表面保持一致。以相移周期条件处理具有不同相位的叶片振动问题,同时发展了伪边界形状修正方法以简化周期边界条件处理。