非均匀来流稳态畸变流场的数值计算

本把TEACH计算机程序推广到三维，用定常，绝热，不可压的时均三维N－S方程和连续方程，并运用紊流动能和紊流耗散率的微分输运方程封闭雷诺方程，模拟计算了非均匀来流稳定畸变流场。数值计算中除了对通用变量进行欠松驰，还对修正压力通过惯量进行了欠松驰，惯量取为压力修正方程中的系数|b|，保证了数值计算的稳定性和收敛性。在进口马赫数小于0．35时，流场畸变指数的计算值和试验值符合很好。     

离心叶轮出口三维流场研究

文发展了一种通过求解不可压缩、定常、时均 N- S方程、连续性方程和 k- ε湍流模型方程分析离心叶轮出口三维流场的方法。该方法的差分计算网格通过代数变换关系生成 ,可方便地应用于离心叶轮设计分析。为检验本文分析方法的有效性 ,对一个三元闭式离心压缩机叶轮出口流场进行了实验测量。计算与实验比较表明 ,本文方法预示出的平均流场与实验数据符合良好     

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

本文采用时均Ｎ－Ｓ方程和Ｂａｌｄｗｉｎ／Ｌｏｍａｘ代数湍流模型计算了典型拉伐尔喷管内正激波与湍流边界层干扰流场。     

轴对称弯道流动的N—S方程数值解

利用涡量—流函数方法求解离心压气机进气弯道中的轴对称粘性流场。本文推导出了q—τ湍流模型方程在正交曲线坐标系下的表达式,并用以预测湍流流场。文中在涡量方程及q—τ方程中引进非定常项,采用时间推进法求解。通过对二维直通道及复杂边界的轴对称弯道内的层流、湍流流动的计算,获得了令人满意的结果。      

考虑传质传热效应的翼型积冰计算

本文采用计算流体力学的方法对NACA0012翼型前缘涉及传质传热效应的glaze ice积冰进行了预测.在流场计算中,用有限体积法对二维定常不可压粘流的时均N-S方程进行离散求解;采用四阶龙戈-库塔(RK)法求解水滴运动方程;通过求解在翼型表面控制体积内积冰过程中所遵循的质量和能量守恒方程,并且假定冰沿着与当地翼型表面法向一致的方向增长,来预测积冰的形状.文中同时对结有glaze ice冰翼型的气动特性进行了计算与分析.     

翼型积冰的数值模拟

本文采用计算流体力学的方法对NACA0012翼型前沿的Rimeice积冰进行了预测。在流场计算中,用有限体积(FV)法对二维定常不可压粘流的时均N-S方程进行离散求解;采用四阶龙戈-库塔(RK)法求解水滴运动方程;假定水滴在与翼型相碰撞的点处完全凝结,并且冰沿着与当地翼型表面法向一致的方向增长,以此来预测积冰的形状。文中同时对迎角为4°时结冰翼型的气动特性进行了分析。     

用粘性／非粘性相互作用法求解机翼跨声速分离流

提供了一种计算机翼的跨声速绕流的粘性／非粘性相互作用的计算方法,包括无粘流场计算,混合边界层计算及两者之间的相互作用,其中三维混合边界的计算包括了层流边界层,转捩区,湍流边界层和分离流的积分方法计算了,特别是在靠近分离的区域采用边界层反方法计算,无粘流场由全速势方程计算得到,通过粘流无粘流耦合迭代求得了Ｍ６机翼跨声速绕流的收敛解,与实验结果比较,吻合得较好,本方法能够计算出激波诱导分离泡和后缘分离     

高超声速钝头双锥体高温流场辐射研究

研究再入高温绕流及尾流流场的辐射特性对高超声速飞行器突防和识别技术的发展具有重要意义。本文基于Navier-Stokes方程和辐射传输方程,建立了沿任意视角方向的再入体全目标区流场(包括绕流及尾迹流场)辐射特性计算方法,计算分析了在0.4~12.0μm光谱范围内钝头双锥体高温化学反应流场辐射在不同视角(0°~180°)方向以及不同再入高度(40~60km)下的变化特征。     

平面压气机叶栅通道分离流动高精度直接数值模拟

采用高精度有限差分格式直接求解二维Navier-Stokes方程组,数值模拟V103平面压气机叶栅分离流动,数值结果表明:在瞬时流场中,吸力面后部发生流动分离,在分离区前端存在大尺度分离涡,分离涡下游是由二次涡和脱落涡交替形成的涡串,直至叶片尾缘,形成以脱落涡为主结构的尾迹;在时均流场中,吸力面后部存在短分离泡,分离区压力分布存在明显压力平台。与逆压力梯度下平板边界层分离流动相比,瞬时和时均流场结构相似;叶栅通道内无量纲涡脱落频率是前者的两倍。与文献计算结果对比表明:叶片表面压力分布除分离区外吻合很好;非定常计算所得分离区轴向长度比定常计算大41%。在分离区内三个二阶统计量均达到最大值,表明流场强非定常性集中在分离区。     

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

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

一种适用于颗粒非规则分布的阻止SPH数值断裂的方法

 〖HT5”H〗摘〓要〖HTSS〗: 光滑粒子流体动力学(SPH)方法已广泛应用于航空航天领域中大变形的碰撞计算。数值断裂是SPH计算中一种常见问题,可严重影响计算结果。把适用于颗粒规则分布的阻止SPH数值断裂的方法推广到颗粒非规则分布的情形,得到了一种适用于颗粒非规则分布的阻止数值断裂的加颗粒方法。然后将该方法应用于存在数值断裂的悬臂梁弯曲仿真计算,验证了其有效性。最后把该方法应用于弹丸撞击飞机蒙皮的非规则颗粒模型计算,并和规则颗粒模型得到的结果及试验结果进行比较。对比结果表明：对于非规则分布的颗粒模型,该加颗粒方法可有效地阻止数值断裂,提高计算精度,并可较为普遍使用。     

大攻角机翼的气动弹性计算方法

 本文利用大攻角带分离流机翼的完全非定常非线性气动力计算方法,通过与机翼运动方程的同时求解,在时间域内实现了大攻角机翼非线性气动弹性的数值模拟,根据弹性机翼各种状态下的运动过程,可以得到大攻角机翼的颤振速度等重要参数以及亚临界、超临界等飞行状态的运动规律。算例结果表明,大攻角下机翼的气动弹性问题需引起足够的重视。     

一种求解全位势方程的跨音松弛法

 本文发展了一种适用于求解二维跨音速叶栅流场的松弛方法。采用流线和平行于y-轴的直线族组成的非正交网格系统,在物理平面中求解全位势方程,沿流线构造相关差分格式。数值结果表明,本方法具有收敛速度快、差分格式稳定的优点,松弛150步以后,便可获得基本满意的结果。计算结果和物理分析一致,和试验数据相近。     

A second-order numerical method for elliptic equations with singular sources using local flter

The presence of Dirac delta function in differential equation can lead to a discontinuity,which may degrade the accuracy of related numerical methods.To improve the accuracy,a secondorder numerical method for elliptic equations with singular sources is introduced by employing a local kernel flter.In this method,the discontinuous equation is convoluted with the kernel function to obtain a more regular one.Then the original equation is replaced by this fltered equation around the singular points,to obtain discrete numerical form.The unchanged equations at the other points are discretized by using a central difference scheme.1D and 2D examples are carried out to validate the correctness and accuracy of the present method.The results show that a second-order of accuracy can be obtained in the fltering framework with an appropriate integration rule.Furthermore,the present method does not need any jump condition,and also has extremely simple form that can be easily extended to high dimensional cases and complex geometry.     

化学非平衡喷管流计算方法的研究

本文提供了考虑化学反应非平衡和振动非平衡的喷管流动计算方法。在忽略扩散、粘性和热传导及辐射损失的假设下，得到完全气体一维定常流动方程，再加上振动能松驰方程和化学反应速率方程即可数值积分求解。采用扰动法得到接近平衡流的一组积分初值。采用“逆积分方法”求解通过喉道。计算结果是令人满意的。     

旋转变厚度圆盘的样条积分方程法

一、前言 对于旋转变厚度圆盘的应力分析,一般可用有限元法或差分法,但相对于边界元来说精度较差,若要提高精度就必须将单元分细,随之而来的是方程阶数的增加和计算时间的增长。为了克服上述缺点,本文提出了一种新的样条积分方程法,直接利用弹性力学二维问题的环基本解建立积分方程并用样条函数进行求解。这样只需划分少量单元,就可获得精度很高的数值解,可大大节约计算机内存和计算时间。     

采用二阶迎风格式的 PISO 算法研究

在ＰＩＳＯ算法基础上,发展了一种新的二阶迎风格式以提高计算精度又保证算法的稳定性。新格式应用于液体火箭发动机燃烧流动过程计算,其数值稳定性很好。     

任意非正交曲线坐标系中湍流的数值计算

本文旨在运用张量分析理论,结合成熟的SIMPLE算法,发展一种适用于复杂几何形状内湍流流场数值求解的方法。假设流动为二维单连通域内的等温、不可压、定常、无化学反应、单组分湍流流动,采用双方程湍流模型。由雷诺平均的湍流控制方程的基本形式[1,2],结合张量分析理论。      

数值求解可压缩N—S方程的增量迭代法

本文在Beam-Warming格式的基础上,引入通量分裂的概念,应用行高斯迭代技术提出了一个二阶隐式的求解BFC坐标变换下的N-S方程的增量迭代法。它消除了因式分解带来的寄生误差,合理地处理了控制方程中的粘性项,具有较好的收敛性和稳定性。二维Poissuill流动和Couette流动的计算结果和精确解吻合良好。S形进气道的数值解反映出的现象和实验结果是一致的。     

静止叶栅三维粘性数值模拟

应用Ｎ－ Ｓ方程压力修正算法和ｋ－ ε两方程湍流模型,对端弯叶片和基准（无端弯）叶片的平面叶栅进行了三维粘性数值模拟。在采用简单交错网格的情况下,通过对计算控制体通量的速度进行修正的方法解决了Ｎ－ Ｓ方程压力修法算中遇到的压力场和速度场呈锯齿状分布的问题。用实验测量值对本数值模拟方法的精度和可靠性进行了检验