后向台阶横向喷流干扰流场的数值研究

 用有限差分法求解层流全N-S方程,对超音速后向台阶绕流与台阶后横向喷流干扰流场进行了数值模拟。不带喷流的超音速后向台阶绕流数值模拟结果与实验及其他计算结果吻合较好。最后给出横向喷流压力比p_j/p_∞=3.033～7.033和喷口宽度w=0.1h～0.2h(h为台阶高度)时的数值模拟结果。得到了合理的流场结构;并可看到喷流压力比和喷口宽度对流场特性有显著影响,在计算参数范围内,底部压力与系数Q(=wp_j/hp_∞)几乎成线性关系。     

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

直接的矩阵分裂计算跨音速喷管流场

给出了分离系数矩阵求解欧拉方程组时系数矩阵的直接分离公式，减少和简化了应用该方法时的矩阵运算，参照Ｂｅａｍ－Ｗａｒｍｉｎｇ显式格式构造与ＳＣＭ方法一致的内点差分格式，并分析了它的数值特性，对有斜激波的跨音速内流场进行了数值仿真。计算中采用变系数的当地时间步长，加快了流场达到的收敛速度，计算无需引入任何人工参数，计算值与实验结果吻合。     

紧致格式数值模拟超音速粘性绕流问题

通过数值求解可压Navier-Stokes方程的方法计算了球锥粘性绕流问题。方程中粘性项按通常办法逼近,为改善收敛速度利用了文[1,2]中之算子附加修正方法。方程中之流向导数利用了四阶精度的紧致差分。这是由于该方向上流动参量光滑和边界条件易于处理的缘故。在法向方向上利用了中心差分。对差分方程之隐式部分利用了近似因式分解法。由于选用了文[3]中之特殊Jacobiall系数矩阵分裂法而使计算工作量大为减少。用这种方法计算了超音速粘性绕流问题。计算结果与实验结果相符很好。     

横向喷流与超音速主流干扰流场的数值研究

 采用欧拉方程与可压缩涡对模型相结合的数值方法,研究了超音速主流与横向喷流干扰流场的气动特性。喷流横截面用一对反向旋转的可压缩涡来模拟。采用有限元Osher格式及非结构网格数值求解二维欧拉方程以得到非等熵外流场的气动特性及弓波位置。将欧拉方程数值解结果与可压缩涡对流场解结果相结合以研究喷流的轨迹及喷流横截面的变化。最后给出了数值模拟结果并与工程估算结果进行了比较,表明本文所述方法是可行的。     

侧向喷流干扰工程估算模型研究

喷流的直接力控制适用于全速域和全空域，具有响应快、效率高等特点，是控制飞行器飞行的有效手段之一。喷流与飞行器流场的相互干扰十分复杂，准确地获得喷流干扰气动力非常困难。笔者调研了国内外大量的喷流实验和局部凸起物干扰的分离流动实验，利用二维平板及轴对称凸起物绕流的分离特性，初步建立了计算喷流干扰的工程模型。该模型可计算喷流穿透高度、分离区长度、二维分离区边界、分离区平台压力系数，并可对分离特性进行三维修正。利用这一模型，笔者进行了某飞行器的超声速喷流的气动干扰特性计算，并与数值计算结果进行了比较。     

定常轴对称超音速喷管和超音速射流的计算

本文叙述轴对称超音速喷管和超音速射流的数值计算。内点用MacCormack有限差分格式,分步加粘性;边界用从特征关系列出的二阶精度格式;喷管口用Prandtl-Meyer流公式;射流边界用一个迭代过程确定。对有激波喷管和无激波喷管做了大量的数值实验,并作了比较和讨论,给出了三个喷管的计算结果。     

火箭级间分离喷流干扰数值模拟与风洞试验研究

针对典型风洞试验状态,采用有限体积法求解三维雷诺平均Navier-Stokes(RANS)方程,数值模拟了给定级间距下火箭级间分离喷流干扰流动情况,获得了清晰的流场结构和有/无喷流干扰时两级气动特性变化规律,并对仅主机喷流和主机与四游机同时喷流的情况进行了计算,数值模拟研究结果为风洞试验方案的设计提供了参考依据。在典型状态数值模拟研究的基础上采用了冷喷流模拟技术对多个状态下火箭级间分离喷流干扰特性进行了风洞试验研究,获得了两级在有/无喷流、同轴变攻角情况下的气动力系数,研究结果表明,级间喷流干扰会对两级气动特性带来较大影响。相同状态下计算得到的两级气动特性变化与试验符合较好,表明了针对风洞试验状态的预先数值模拟对试验方案的设计验证可提供有益的参考。     

壁面喷流可控环量机翼的数值计算

本文提出一种适用于可控环量机翼研究的势流/粘性流的数值解法,包括钝体分离绕流的一种势流解,以及改进的Spalding统一边界层理论的解。通过对具有壁面切向喷流可控环量圆柱绕流和类椭圆翼绕流所作出的数值计算表明,在其计算模型确定后,就不必依赖实验数据来模拟实际问题,计算结果与实验数据符合较好。     

超音速来流中侧向喷流干扰流场的数值模拟

通过求解三维可压缩N S方程组,对尖锥/圆柱组合体模型在超音速绕流中的侧向喷流干扰流场进行了模拟,数值方法采用有限体积法,差分格式采用NND格式。分析了干扰流场的结构,研究了干扰效应对气动性能的影响,得到来流Ma∞=3 3和4 5时力放大因子和干扰力矩系数随攻角变化的规律,计算结果与实验一致。     

含侧喷流弹丸三维绕流干扰流场数值分析

从全Navier-Stokes方程出发，应用三维隐式有限体积TVD格式，在超声速和零攻角条件下，分别对无侧喷流和含侧喷流的低阻远程弹丸三维粘性绕流流场进行了数值模拟，为避免分区计算而用代数方法生成了弹丸绕流场O型网络，计算结果得到了清晰的流场波系结构，并从流场波系结构、弹体和弹底表面压力分布等角度将两者进行了对比与分析，表明侧喷流不仅提供增程推力，而且改变了弹丸阻力特性，其结果可用于研究弹丸的火箭增程技术。     

超声速横喷干扰湍流场数值模拟

本文探讨了超声速平板横向喷流干扰湍流流场的数值模拟方法。分别采用Roe、Harten-Yee及AUSM 格式对三维雷诺平均Navier-Stokes方程进行离散,湍流模拟应用Baldwin-Lomax代数模型,讨论了B-L模型中Fy最大值的不同选取方式。分析了数值格式、网格划分对计算结果的影响。将来流马赫数为5、喷流马赫数为3的实验结果与计算结果进行比较,选择了与实验符合较好的数值模拟方法,为进一步利用数值计算方法研究喷流干扰流场特性奠定基础。     

用于粘性可压缩流动数值计算的SIMPLE方法

对SIMPLE方法加以推广,使之适用于可压缩粘性流动的数值计算。采用非交错网格技术,推出三维任意曲线坐标系下可压缩形式SIMPLE方法的计算方程。对亚音速、跨音速及超音速等5个流场进行了数值计算并与有关文献及实验数据进行比较。      

超音速三维分离系数矩阵差分方法

本文提供了计算超音速无粘、定常、三维非对称体分离系数矩阵(SCM)差分方法,这种方法是基于数学特征理论,在分离系数差分法中,其系数是按照特征线的斜率而分离的,这些系数乘以合理的单向差分,前差与负特征值相关,而后差与正特征值相关。 本文给出了球头钝锥的算例,并与文献[2]作了比较,发现SCM方法在较少的网格点情况下,仍能获得较高的精度。     

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

非定常低雷诺数流数值模拟方法研究

低雷诺数翼型流动的一个重要特征是分离泡的产生,并且分离泡通常是非定常的。基于有限体积法,N-S方程和S-A模型的控制方程非耦合求解,时间均采用牛顿型LU-SGS方法推进,无粘项和粘性项用中心差分方法离散计算,在计算非定常流场时,采用伪时间子迭代方法,湍流模型计算时采用在固定点转捩,研究了二维翼型低雷诺数流粘性流动的数值模拟方法,并利用所开发的非定常可压缩粘性流计算方法,数值模拟了在低雷诺数下两种翼型的非定常流动。计算结果表明：气动力系数及分离泡呈规则的周期性变化,证明方法是成功的、可行的。     

复杂回旋流场的数值模拟和实验研究

本文通过四角喷射炉膛内回旋流场的数值计算和实验测试的对比研究，比较了几种计算方法的预估结果与实验结果的符合程序。结果表明：为克服当流动方向与网格线交角过大时所产和珠过分数值扩散，引入流向平衡粘性方法与κ－ε湍流模型相结合可获得与实验结果相一致的计算结果，而采用混合分格式加κ－ε模型所得到的计算结果比流向平衡粘性方法加等效常湍流粘性系数这样的简化处理的计算结果还要差。     

横向喷流对鸭式导弹滚转特性影响研究

通过数值方法求解三维可压缩雷诺平均N-S方程,对导弹横向喷流的干扰流场进行了数值模拟,计算结果与实验数据吻合较好,基本验证了计算方法在横向喷流复杂流场数值模拟方面的有效性。在此基础上,对固定尾翼鸭式布局导弹亚、跨、超声速流场进行了数值模拟,并且计算分析了横向喷流对鸭式布局导弹滚转控制特性的影响,计算结果表明,横向喷流可以有效提高鸭式布局导弹的滚转控制能力。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

二维喷管分离流线性/非线性湍流模式的研究

采用三个线性和七个非线性涡粘性湍流模式，对二维非对称喷管中的不可压缩分离流进行了研究。通过将数值计算结果和实验结果进行比较，对有关的湍流模式进行了评估和分析。计算结果表明，非线性模式的模化系数与平均流动应变不变量以及旋转不变量有关，反映了湍流的各向异性，比线性模式优越得多。另外，所有线性／非线性涡粘性湍流模式不能够捕捉流动的非定常分离过程。