抛物化稳定性方程在曲面边界层中的应用

为了研究抛物化稳定性方程在曲面边界层中的应用,选取了等曲率圆环管道、等曲率板和NACA0012翼型3种典型的曲面边界层.通过计算小幅值扰动波的演化,抛物化稳定性方程的计算结果和线性稳定性理论、数值模拟扰动方程结果相符,说明可以使用抛物化稳定性方程研究曲面边界层的小幅值扰动波的演化和稳定性分析;通过计算有限幅值扰动波的演化,线性阶段抛物化稳定性方程计算结果和扰动方程相符,在非线性很强的阶段,抛物化稳定性方程计算发散,发散的位置可作为转捩位置.      

湍流通用物理方程的推导及讨论

湍流的通用方程已经获得。利用变形率张量与摄动变形位移的点积分解湍流脉动速度,并使用密度加权的侧偏系综平均方法取代Reynolds平均方法,对Navier-Stokes方程进行平均。在二阶精度上没有出现未知相关量,从而获得了具有二阶精度的封闭的湍流通用物理方程。方程中无相关量,不引入任何经验系数,无需构造任何湍流模型,即可描述任意复杂湍流场并解释湍流物理本质。方程在长时标上描述湍流平均行为,在短时标上描述相应尺度层次的湍流结构现象,从而统一了湍流的统计及结构二大学派。方程简化后得到湍流的规范方程,是混合的Burgers-Korteweg-de Vries方程,对该方程的分析给出了湍流转捩准则。文中对湍流的间歇性、拟序结构、各向异性、湍流能量逆转等现象的物理本质进行了初步讨论。      

扩压器内二维湍流流动的时间推进数值解

ＡＦ有限差分理论。指出薄层Ｎ－Ｓ方程抛物化后右端粘性项是主元素变量矢的一个简单线性函数,这对实际计算很有用。应用隐式ＡＦ方法计算了为发动机进气道设计的典型扩压器内的复杂湍流流场。研究表明,计算结果与实验结果基本一致,这为三维非定常湍流流动进一步研究提供了基础。     

改进的小波均匀化方法求解小周期系数抛物型方程

对传统的小波数值均匀化方法通过加入逐层修正的步骤进行了改进,并将改进后的方法用于求解具有小周期系数的抛物型方程。数值算例结果表明：在保持较高计算效率的同时,改进后的方法比改进前方法在精度方面有明显的提高。因此,改进后方法是小周期系数抛物型方程的一种高效率、高精度的数值求解方法。     

用二维N—S方程计算汽车流场的地面效应

采用有限体积方法求解了Navier-Stokes方程,计算了绕二维类汽车物体的绕流并研究了地面效应。在非正交网格系统上以强守恒形式构造出有限体积方程,采用压力校正算法来求解所得到的方程。对于湍流情况,采用了双方程模型的κ_ε模型。 为了验证采用本程序计算移动壁面效应的正确性,本文给出了Couette流动的计算结果。 本文也给出了汽车绕流的计算结果,它包括层流情况,从汽车头部开始由层流转为湍流的混合流动情况及湍流情况。对于每一种情况都比较了移动壁面的计算结果和固定壁面的计算结果。     

使用GAO-YONG方程组对湍流边界层拟序结构特性的研究

与以往的湍流模型不同,GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。其级数形式的能量方程与非线性现象多尺度层次相对应,具备了描述湍流统计平均和拟序结构的双重功能。本文采用交错网格和SIMPLE方法求解GAO-YONG不可压湍流方程组,考察了GAO-YONG方程组中能量方程阶数以及计算网格对拟序结构的影响,证明了该方程在描述拟序结构方面具有一定的适定性。      

湍流模型源项的数值处理方法对计算结果的影响研究

涡粘性湍流模型广泛应用于各种工程湍流问题的计算。对于大多数湍流模型,在湍流控制方程的右端会出现源项,其刚性给数值计算带来很大的影响。从源项的物理意义出发分析了源项导致求解困难的原因,并在通用的求解湍流模型方程的点隐法基础上,以源项弱刚性的S-A湍流模型和源项强刚性的k-ωSST湍流模型为例,论述了生成项和耗散项之间的平衡关系对数值模拟的重要意义,给出了针对不同源项的具体处理方法。对RAE2882翼型跨声速流动算例的模拟结果表明,这些处理方法有效的提高了控制方程组的计算稳定性。     

桥梁断面静风荷载的格子Boltzmann方法数值计算

通过引入反映湍流涡粘性的湍流松弛,得到了模拟高雷诺数湍流的BGK方程.在速度相空间、物理空间和时间上对BGK方程进行离散得到了三维十九速离散速度模型;结合分区计算技术,设计了格子Boltzmann并行算法;根据亚格子Smagorinsky模型,提出了直接从粒子分布函数计算湍流松弛时间的方法.用开发的并行计算程序对分体双箱截面和闭口箱梁截面的静风荷载进行了数值识别,得到的静力三分力系数和流场压力分布与风洞试验结果及CFD宏观方法计算结果吻合,并从表面压力分布入手分析了两种桥梁截面的绕流特点.     

关于湍流问题的注记

指出N-S方程是研究湍流所遵循的基本方程, Reynolds方程反映了在湍流中随机性与确定性共存;Burgers-KdV方程与湍流机理无关, 不能作为湍流研究的“新方程”      

二维带动力吸气式高超声速飞行器绕流的PNS-NS混合求解

为了提高吸气式高超声速飞行器绕流的求解效率,采用空间推进求解抛物化PNS(parabolizedNavier-Stokes)方程和时间迭代求解Navier-Stokes(N-S)方程的混合计算流体动力学(CFD)方法来求解高超声速飞行器整机绕流.在超声速占主导的流动区域采用空间推进求解抛物化N-S方程的方法,在亚声速和分离区采用时间迭代求解N-S方程的方法.对于求解二维带化学反应的吸气式高超声速飞行器绕流,混合CFD方法和完全时间迭代方法相比,可得到同等准确的数值模拟结果,并且求解效率提高了数倍.      

半波GLM方法封闭的湍流方程对转捩过程的数值模拟

ＴＲＡＮＳＩＴＩＯＮＳＩＭＵＬＡＴＩＯＮＵＳＩＮＧＧＬＭ－ＥＭＢＡＳＥＤＣＬＯＳＥＤＥＱＵＡＴＩＯＮＳＯＦＴＵＲＢＵＬＥＮＣＥＧａｏＧｅ（４ｔｈＤｅｐｔ．ＢｅｉｊｉｎｇＵｎｉｖｅｒｓｉｔｙｏｆＡｅｒｏｎａｕｔｉｃｓａｎｄＡｓｔｒｏｎａｕｔｉｃｓ，Ｂｅ...     

本构关系对连续介质损伤力学基本方程适定性的影响——一维弹性损伤情况

许多工程材料当变形超过某一极限时由于损伤都会出现应变软化行为。首先从非局部理论出发,推导了应变梯度损伤本构方程;然后利用一阶拟线性偏微分方程组的特征理论,在一维弹性损伤情况下分析了两种不同的本构模型,即 Kachanov损伤本构方程与应变梯度损伤本构方程,对连续介质损伤力学基本方程适定性的影响。结果表明,当损伤发展时,与 Kachanov损伤本构模型相关的连续介质损伤力学的基本方程是不适定的;而与应变梯度损伤本构模型相关的,则无论损伤是否发展,其基本方程始终是适定的。这说明在连续介质损伤力学的本构方程中必须考虑材料内部微结构的尺度效应。     

复杂外形飞行器多块结构贴体网格生成

通过求解椭圆型方程生成三维贴体坐标网格,采用分块技术生成多块结构网格体系。每个子区域转换成计算空间内的立体域,具有自身的曲线坐标系,从而所形成的整个网格系统不会过分扭曲和不规则,这将有利于对于复杂边界的真实飞行器绕流流场的求解。文中讨论了控制函数对网格生成的影响。     

钝锥大攻角超声速分离流场的数值模拟

最近,张涵信等人在传统的Beam-Warming隐式、无迭代、空间推进技术的基础上,根据边界层方程的性质,设计了一种可用小步长推进求解抛物化NS(PNS)方程、而不会引起解的“漂移现象”发生的方法。这种方法对轴对称流动的计算是成功的。本文就是将这一思想推广应用于大攻角有周向分离的流场计算。求解的区域为具有薄亚声速层的有粘与无粘干扰的整个激波层内的流场。在对攻角α=0°和α=20°的球钝锥的计算中,关于壁面上的压力、热流率及流场的涡旋结构均得到了满意的结果。文中特别研究了钝锥大攻角绕流的流动分离图象。 为了增强块三对角矩阵的主对角优势,通常在差分方程的左端附加二阶增量项。本文以选取适当小的推进步长的方法来达到增强主对角优势的目的,不需再附加二阶增量项,从而提高了解的精度。     

物面突变引起的二维分离流动的数值模拟

本文用一个单步二阶精度格式,对由物体截面突变引起的二维超声速分离流动进行了数值模拟。计算了近尾迹流场,并详细研究了来流M数和Re数对绕后台阶分离流场结构的影响,得到了与理论和前人计算一致的结果。另外,对超声速绕竖舨流动进行了数值计算,给出了细致的流场结构图画。     

跨音速操纵面嗡鸣Hopf分叉分析及结构参数对嗡鸣特性影响的研究

采用Ｈｏｐｆ分叉分析方法,对跨音速操纵面嗡鸣问题进行了研究。用二维守恒型非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ 方程计算了舵面振动时的非定常气动力,计算了出现Ｈｏｐｆ分叉的Ｍａ 数,研究了一些参数对嗡鸣的影响。将Ｈｏｐｆ分叉的计算结果与嗡鸣的时间历程计算做了对比,两者相当吻合,与飞行试验结果相比,也比较接近。     

旋翼桨叶绕流Euler方程数值模拟

本研究应用格心格式有限体积法求解欧拉方程,模拟绕旋翼浆叶的流动。旋翼尾流的作用通过自由尾流分析来求解局部的诱导下洗速度以修正翼型攻角。选择O—H型式生成弦向和展向网格。通过算例计算为进一步寻求有效的Euler或NS方程数值解积累了经验。     

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

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

三维不可压边界层抛物化稳定性方程的椭圆特性研究

针对三维平行不可压边界层流动的抛物化稳定性方程,研究其剩余的椭圆特性,提出了消除该特性的具体方法,并对抛物化稳定性方程解与全椭圆问题解的一致性进行了证明。