三维不可压缩层流、转捩、湍流边界层的数值计算

本文提出了一种计算三维不可压缩边界层的数值解法,其中包括对层流边界层的微分法,对湍流边界层的积分法以及利用Orr-Sommerfeld方程计算转捩的数值解,当扰动振幅达到e~(13.5)时发生转捩。 按照上述方法,本文计算了ONERA M6后掠角为30°的机翼,计算结果分别和Cebeci和Cousteix的结果作了对比,二者吻合较好。     

跨声速翼型绕流的Euler／边界层方程干扰数值解

本文利用Euler方程和可压缩湍流边界层积分方程研究绕跨声速翼型的有粘与无粘强干扰流动。应用有限差分法在贴体的网格上求解时间相关的Euler方程,以剪功积分方法求解翼面贴附和分离湍流边界层流动,并引入一个松弛方程描述剪应力对上游湍流历程的延迟响应。有粘/无粘干扰采用表面源模型。计算结果表明,对翼面存在强干扰流动情况,获得了与实验值基本吻合的结果。     

跨声速机翼流场的N—S方程计算

本文给出了三维机翼粘性流场的数值模拟过程，控制方程为时均化的三维可压缩薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，湍流模型采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数湍流模型，空间离散采用中心有限体积格式，时间采用Ｒｕｎｇｅ－Ｋｕｔｔａ多齿格式进行式积分。     

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

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

绕扁长椭球体三维分离的积分法

本文利用三维积分边界层法并结合流线法决定涡层型流动分离。由三维不可压湍流边界层方程求得边界层参数后,按流线法计算粘性流线。作为算例,研究了一带迎角的扁长椭球体在不可压缩湍流里的绕流,经与水洞里所观察到的α=15°流型的对比,二者基本一致。     

弯曲进气道三维边界层预测

本文研究了弯曲进气道内三维可压缩湍流边界层的预测方法。采用非正交曲线坐标系中三维边界层积分方程,并引入掺混系数微分方程,用滞后掺混法求解,数值方法为预估-校正法,用三次样条函数解决了离散数据的插值、求导的精度问题。所编程序可在微机上运行。实例计算清楚地揭示了弯曲进气道内强横向流动对边界层发展的影响,计算所得的进气道总压恢复系数与实验值的误差为1～2％。     

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

垫块法半模型风洞试验的数值模拟研究

以M6机翼半模型为算例,采用数值方法模拟垫块法风洞试验洞壁边界层及垫块厚度对半模型试验结果的影响.控制方程采用全场三维可压、雷诺平均的Navier-Stokes方程;计算格式采用Jameson中心有限体积、多步Runge-Kutta时间步长格式进行求解;湍流模型采用Baldwin-Lomax模型.通过本文的研究,可以为跨声速风洞半模型试验垫块厚度的选取提供计算方法和设计依据.     

可压缩紊流分离边界层的计算

本文介绍了一种求解二维可压缩紊流分离边界层的正反积分方法。该方法以边界层动量积分方程和平均流动能积分方程作为基本方程,采用适用于贴附和分离紊流的Swafford速度型解析式,对比计算了四种代数涡流粘性模型对沿流向各点的耗散积分及诸边界层特征参数的影响。计算结果和存在激波边界层较强干扰的跨音速翼型边界层的测量数据吻合。该方法可扩展用于计算翼型的尾迹特性。     

微可压缩模型(SCM)与可压缩N-S方程数值计算对比研究

本文开展了微可压缩模型与可压缩N-S方程之间的对比研究.通过对层流平板边界层和椭球高雷诺数湍流绕流的数值模拟,比较了微可压缩模型和可压缩N-S方程之间的计算精度和效率.计算结果表明,微可压缩模型和可压缩N-S方程具有相同的计算精度,但是微可压缩模型具有更高的计算效率,内存和CPU时间都节省20%左右.     

后掠翼可压缩三维边界层稳定性计算

采用可压缩线化稳定性理论对后掠机翼的层流边界层稳定性进行了计算研究。对三维边界层中三维扰动波的可压缩稳定性方程组特征问题采用相关两步格式求解。与常用的初值法相比,具有不需预估特征值和快速收敛的特点。为得到渐近边界条件的准确表达式,导出了外边界上方程的解析解,这对精确的特征值计算和分析至关重要。为避免出现伪不稳定模态,选择了合适的稳定性方程组形式。最后,通过层流控制机翼和自然层流机翼的算例,研究了后掠翼上典型Tollmien-Schlichting(T-S)型扰动和横流(C-F)扰动稳定性问题,其结果与已有计算结果符合甚好。     

三维跨声速粘性-无粘相互作用的半逆方法

本文为计算三维跨声速粘性-无粘相互作用提供一半逆方法版本。湍流边界层采用积分反方法。证明了对于跨声速流,选取外部流线角α和等价的不可压缩形参作为输入量,边界层积分方程总是双曲型的。无粘外流采用FL027程序。通过粘性-无粘迭代求得小分离区情况下的与实验符合的收敛解。     

三维跨声速流粘性—无粘相互作用问题的分析和计算

本文对跨声速粘性-无粘相互作用问题进行了分析,论证了如何选择湍流边界层反模式的输入量。在FL027多重网格基础上作了局部加密。计算例子表明,本文建议的半逆方法适用于弱激波情形,可用于实际设计之中。     

层流控制和自然层流机翼可压缩边界层计算和应用

本文对后掠层流机翼可压缩边界层进行了计算研究。引入新的相似变换形式,避免了Kaups-Cebeci方法在变换中可能出现的奇性,导出了在新的相似变换下的偏微分方程组。对层流控制(LFC)和自然层流(NLF)机翼进行了计算,得到了满意的结果。     

二维可压粘性叶栅流场计算的预测修正迭代多重网格法

 本文中提出了一种解可压N-S方程的快速数值法。它由预测步,修正步,迭代步以及多重网格法构成。叶栅出口采用了一种使数值解收敛性进一步加快的技术,我们称之为定常无反射边界条件技术。本文应用了k-ε双方程紊流模型。 对若干叶栅在-13°～8°的攻角范围内进行了数值试验,并与对应的物理实验结果进行了比较。结果是满意的。     

运用GAO-YONG湍流模型对扩压器内跨声速流动的数值模拟

运用GAO-YONG可压缩湍流方程组,采用同位网格SIMPLE算法,对扩压器跨声速流动中的二维激波/湍流边界层干扰现象进行了数值模拟。将计算得到的流场的时均参数与实验值进行比较,数值模拟结果在激波强度、壁面压力分布以及分离点和再附点位置等方面,与实验值吻合较好,表明GAO-YONG可压缩湍流方程组能够比较准确的模拟较强激波/湍流边界层干扰流动,从而进一步为GAO-YONG湍流模型的正确性及其在可压缩流场模拟方面的适用性提供了佐证。     

一维可压缩流的有限元状态空间模型

通过引入位置交错的两种有限控制体积单元——状态单元和速度单元,对其分别应用质量、能量守恒方程和动量守恒方程, 建立了一维可压缩流的有限元状态空间模型, 模型单元边界条件的处理方式简单有效, 适用于处理任意的物理边界条件。与分析解的比较表明该模型是精确的     

Euler方程与边界层积分方程耦合求解跨音速翼型绕流

 应用Euler方程求解跨音速翼型特性时考虑了粘性影响,粘性影响是通过边界层动量和能量积分方程求解的,即粘流/无粘流迭代方法。其中Euler方程采用LU-ADI方法求解;边界层方程均由正解法过渡到反解法,以解决强激波干扰区出现小分离泡的计算问题。计算中使用了贴体C网格,通过一定变换使其保持基本正交。计算结果表明,压力分布、摩阻系数分布与实验结果符合较好。