等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

高精度ENO格式在喷管流动模拟中的应用

本文分析了ＥＮＯ格式的特点,并应用于Ｅｕｌｅｒ方程和全ＮＳ方程的迁移项和压力项,模拟了各种喷管流动。首先计算了ＪＰＬ喷管流动,取得了和实验一致的结果,分析了喷管内产生弱激波的原因;其次计算了非定常二维喷管流动,给出了相应的计算结果,并分析了壁面附近流场的拓扑结构;最后计算了摆动喷管的流动及推力矢量计算,取得了推力随摆角变化的规律,和实验结果表明基本一致。     

具有无源控制空腔时正激波/湍流附面层干扰的数值模拟

采用雷诺平均N-S方程和B/L代数湍流模型计算了具有无源控制空腔时正激波/湍流附面层干扰流场。计算与实验结果的比较表明,本文方法可较准确地预测激波结构、激波与附面层干扰区流动基本特征及波后流动分离状态、激波位置、波前马赫数等参数。      

超音速底部流动的大涡模拟研究

超音速底部流动是CFD的难点问题之一。采用基于可压缩流动的大涡模拟方法求解N-S方程,对超音速底部流动进行了模拟研究,分析了其流动结构。计算结果表明：得到的速度型、底面压力分布与实验值吻合,优于雷诺平均方法的结果。     

管内正激波/湍流边界层干扰的数值模拟

本文采用时均Ｎ－Ｓ方程和Ｂａｌｄｗｉｎ／Ｌｏｍａｘ代数湍流模型计算了典型拉伐尔喷管内正激波与湍流边界层干扰流场。计算与实验结果的比较表明，方法可准确地预测激波结构，激波与边界层干扰区流动特征，激波位置、激波前马赫数和壁面压力分布等。     

基于DES的二维和三维空腔流动特性研究

空腔流动广泛存在于航空航天工程中,对其流动特性的研究具有十分重要的工程意义。采用数值方法对比模拟了二维和三维空腔流动,控制方程采用N—S方程,空间离散采用有限体积方法,对流通量计算采用Roe格式,非定常时间离散采用双时间步长方法,湍流粘性计算采用基于SA模型的DES方法。数值计算所得的三维空腔底面压强分布与实验结果一致,所得的二维和三维空腔内的流动结构与相关文献中的分析相吻合。将二维和三维的计算结果进行比较,发现空腔底部压强分布、空腔内的压强振荡、空腔内部流线图及声压级等均有所不同,得出在空腔宽度有限时,横向流动简化了流动结构．削弱了振荡幅值．但使振荡过程复杂化．其三维效应不能忽略．     

点火管流量的确定

本文就固体火箭发动机中采用的点火管提出了一维两相流动模型,并用差分方法求解模型方程,得点火管流量。为了验证计算结果进行了实验。计算和实验结果对比表明趋势一致。本文还给出了药粒粒径和药量对点火管内压力的影响。     

计算绕薄翼型跨音速非定常流的积分方程法

 本文从小扰动方程出发推导出绕薄翼型的跨音速非定常流动的积分方程,为使此方程适用于具有激波的流场引入人工粘性,并对其数值求解。通过算例讨论了人工粘性和计算域大小对计算结果的影响,与实验及其它数值计算结果比较,表明本方法的准确度令人满意。计算量小,收敛性好是此法的特点。     

预处理方法在叶轮机械三维数值模拟中的应用

应用预处理方法求解三维可压缩流动Navier-Stokes方程, 发展了一种适用于叶轮机械的既可以求解跨声速流动也可以求解不可压流动的数值模拟方法, 对NASA Rotor67和北航低速大尺寸压气机算例进行了计算, 计算结果与实验数据相吻合, 表明本文采用的预处理方法具有较大的工程实用性.      

超音进气道的亚音扩压器流场的数值计算和实验研究

本文对二维、亚音扩压管道中的粘性流动进行了数值计算,并对同样的扩压器模型进行了风洞吹风实验。 根据要计算的二维、亚音、壁面扩张角不太大的管道中的粘流,我们对N-S方程进行了抛物化假设,化成可沿流程推进求解的抛物型方程,再用有限差分进行计算。 经与实验结果对比,证明本文采用的结合零方程紊流模型,在几何流线坐标系下求解抛物型方程的方法对解决二维、亚音扩压管中分离区以前的流动问题,可得到满意的结果。     

吹吸气流隔断局部空间的基础研究

本文对用吹吸气流隔断局部空间时的气流运动进行基础研究。本研究的气流流动假定为二维、湍流、等温和稳定流动。动量方程式中的从属变量选用流函数和涡度(ψ—ω),湍流动能方程选用一方程式模型(K—L模型),以此作为吹吸气流流动的数学模型。数值解析结果与实验结果比较表明,用上述方法推导出的数学模型和计算程序能较好的反映实际流动情况。本研究对吹吸气流隔断局部空间的研究、设计、现场调试将具有一定的理论价值和工程应用价值。     

驱动方腔内涡流的数值模拟

本文介绍一种求解涡度流函数形式的定常不可压缩Navier-Stokes方程的半隐式指数型差分格式。涡度方程用上述格式求解,而流函数方程用多层网格法求解。这种组合求解方式具有很好的稳定性及较快的收敛速度。本文对Re=100,400,1000,3000,5000的驱动方腔及Re=400,1000的驱动长方腔进行了数值模拟,结果与现有的计算相吻合。     

激光燃速仪的燃烧室流场分析计算

 Numerical method of flow field of laser servo-controlled strand window bomb was studied and discussed. Stream function and vorticity method of axisymmetric, steady, incompressible flow field was also studied and discussed. Axisymmetric finite difference stream function-vorticity equations on interior points and the points near oblique boundary were derived. Boundary conditions were given and initial interative field was offered before calculation. Cross-bidirectional scanning method was used for increasing parameter transmission rate and relaxation method which could increase convergent stability was applied. At last, the flow field of laser servo-controlled strand window bomb was shown, including motion patterns and vorticity flow scatter diagrams. The solutions were in concordance to the solutions by Mills et al. and the theories by Pan and Acrivos.     

预测轴流压气机叶栅通道端壁附面层及叶片力亏损

本文将二次流理论与端壁附面层理论相关联,提出一种用于预测轴流压气机叶栅通道内的端壁附面层及叶片力亏损的新方法。附面层主流流动采用包含叶片力亏损的动量积分方程,横流流动采用双层速度分布模型,使用二次流分析和有限差分计算方法预测横流外层速度分布,主流流动和横流流动交替计算。用于两种高负荷压气机叶栅的预测结果表明,不仅端壁附面层发展的预测结果与实验结果吻合较好,而且叶片力亏损发展的预测值也与实验值一致性较好;此外,该方法不仅自动,而且还能较准确地预测S形横向流动的发展。     

斜坡增强混合的数值研究

采用隐式方法求解三维可压Ｎ－Ｓ方程及组分方程,数值模拟了后掠、无后掠斜坡的混合增强流场。比较了两种斜坡的混合增强效果及产生的流向涡量的强度。计算结果表明后掠斜坡产生流向涡的强度明显高于无后掠斜坡,其混合增强效果也更好。     

离心压气机S1流面流场的计算

本文对计算轴流压气机 S1流面的流函数矩阵解法的计算机程序作了修改,使之能适合于计算离心压气机的 S1流面。本方法是引入人工密度,积分运动方程求速度场,而后求密度场,进行迭代求解,计算速度非常迅速。对离心压气机的几个算例表明:只要不发生明显分离,计算结果与实验值是符合的。      

波瓣混合流场流向涡的数值分析

借助N-S方程对波瓣喷管混合装置的流场进行了数值分析。采用三维贴体曲线坐标,非交错网格法和K－ε两方程湍流模型,对波瓣喷管内外两侧的气流通道及波瓣喷管下游的混合段流场统一进行了数值计算。得出了速度场和流向涡及正交涡的分布,并对流向涡的发展规律及其对正交涡分布规律的影响作了分析。并对混合管壁面的压力分布和出口速度型与实验值进行了比较。     

弯曲流道内二次流动的运动学分析和 数值可视化方法

通过运动学分析,得到了涡量与二次流之间的关系方程,介绍了弯曲流道内二次流结构的可视化方法.以截面法向涡量为基础,可求解得到与涡量对应的二次流速度场.该方法具有以下几个特点:从运动学求解得到的涡量与二次流流函数之间的方程,直接体现了二次流作为涡的特性,且作为分析流动结构的后处理方法,能够量化地求解出二次流速度以及分布情况.该方法消除了位势流的影响,使二次流显示更加真实清晰,这种方法比速度矢量图法更加准确完善.      

全速位积分方程跨音速绕流的数值解

从全速位方程出发,利用Green公式将其化为激波捕获积分方程和激波装配积分方程,然后离散进行数值解。流场出现激波时,对激波捕获积分方程应用上风技术捕捉到激波,然后应用激波装配技术计算,得到了满意的结果。经算例考核,该方法具有计算区域小,收敛快和CPU时间较少等优点。