可压缩流动体积粘性的连续介质理论及其数值研究

来源 :中国科学院大学 | 被引量 : 0次 | 上传用户:qq231582
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
Navier-Stokes方程组是可压缩流动的数理基础,已被广泛作为研究连续介质流体运动规律的物理模型,其解通常在Stokes假设条件下求得,即体积粘性系数ζ=0。然而,在高速高温流动领域中,常伴随有剧烈的膨胀、压缩和非平衡热力学过程,此时的体积粘性效应表现较为显著。近60年来,分子运动论和声波实验已表明:除了单原子气体外,Stokes假设在高温高频条件下是不成立的。体积粘性是可压缩介质流动的重要物性参数,与流体运动因不可逆过程引起的耗散机制有关,它和剪切粘性分别表征了流体微团所受法向的压缩效应和切向的剪切效应,并完整描述了应力与应变速率之间的本构关系,其理论研究、数值模拟和实验测量至今仍有诸多不足。本文围绕可压缩流动的体积粘性效应开展了基础理论和数值研究,主要内容和结果如下:  (1)基于本构方程和局部热力学平衡条件,建立了体积粘性的连续介质理论,即ζ=Kτ=f(A)pτ。新理论首次揭示了体积粘性系数是体积模量K与弛豫时间τ的乘积,它和分子运动论分别是宏观和微观角度的推导结果,二者仅相差一个系数。在300K附近,考察N2、O2和空气的ζ发现:与分子运动论相比,连续介质理论的结果和实验数据符合得更好。体积粘性系数在转动和振动模态下分别随温度缓慢增加和指数衰减,振动模态下的弛豫时间τ、分子碰撞数Z、ζ和ζ/μ要比转动模态下的高出几个数量级。通过引入过渡函数g(T),从数学上解决了N2、O2和空气体积粘性系数的间断问题,并建立了ζ的连续表达式。  (2)发展了一种新的混合迎风格式——K-CUSP格式。通过K-CUSP分裂思想、矩阵特征分析、数值构造形式、界面声速新定义、数值耗散分析和数值测试等几个方面的研究,结果表明:K-CUSP格式既具有FDS格式的低耗散性和高分辨率,适用于捕捉物理间断、剪切层和粘性边界层等流动结构;又具有FVS格式的计算效率和稳健性,不会产生非物理的马赫干扭曲和“Carbuncle”现象。K-CUSP格式能够缓解过多数值耗散掩盖物理耗散的问题,它可作为体积粘性效应数值研究的一种有效计算方法。  (3)正激波结构的体积粘性效应研究。根据体积粘性系数ζ的分子运动论和连续介质理论模型,采用四阶Runge-Kutta法和打靶法数值研究了氮气的激波结构,结果表明:体积粘性使激波层内的密度、速度等物理量的空间分布变得更为平缓;体积粘性增大激波层内的粘性输运耗散,减小热流和熵流的峰值;与分子运动论相比,连续介质理论给出的激波厚度倒数λ/δ与实验结果符合得更好;当1.2≤Ma≤10时,激波层内的Kn≤0.4。与无体积粘性的情形相比,考虑体积粘性的连续流N-S方程组能给出与实验值更接近的计算结果。  (4)环形激波聚焦的体积粘性效应研究。采用体积粘性系数ζ的连续表达式和3阶K-CUSP格式,数值研究了激波马赫数为2.0的环形激波聚焦问题,结果表明:环形激波会聚中心点的压力、温度和密度随着时间先增加而后缓慢减小;与转动模态相比,振动模态下的体积粘性效应更为显著,体积粘性可使会聚中心点的压力峰值减小23.5%,温度峰值增加13.2%,密度峰值减小32.9%(连续介质理论的结果);当粘性系数比ξ/μ增大时,会聚点的压力峰值减小、温度峰值增加、密度峰值减小。  (5)激波/层流边界层干扰的体积粘性效应研究。根据体积粘性系数ζ的分子运动论和连续介质理论,采用3阶K-CUSP格式数值研究了激波/层流边界层干扰问题,其中Ma∞=2.0,激波入射角θ=32°。数值结果表明:在激波与边界层的粘性干扰区,体积粘性使该区域容易出现逆流速度和回流区,改变了流场特征;体积粘性对绝热壁面的流动变量影响较大,当ζ增加时,边界层内的温度减小、密度增加、压力减小,边界层厚度δ随之减小;体积粘性会增大了粘性干扰区的摩阻系数cf。当体积粘性不可忽略时,可压缩层流边界层的方程组可以推广至包含体积粘性的边界层方程组。
其他文献
盐城发电有限公司110KV系统当有隔离开关操作时,母线CVT二次电压绕组出现击穿保险动作,B相熔丝熔断现象。为此,本文在国内外相关研究的基础上详细分析了隔离开关分合时电弧的产生机理、电弧的辐射特性、电磁干扰在空间的传播过程,最后如何耦合到二次回路中。在理论分析的基础上,通过实际测量得到了具体的数据和电压波形,从而证实了理论分析的正确性。综合理论与实验的研究,我们根据盐城电网的实际情况找到可行性措施
现今,电力行业随着科学技术的革新也在迅速的发展。伴随电能的生产、输送、配电及用电的过程中产生的谐波问题引起了多方面的关注。由谐波所引起的电力系统事故也随之增多。这
该文在Fredholm积分方程组解法统一表达式的基础上着重研究了核函数的构成方法,为方法的具体实施提供了各有关表达式.根据Fredholm积分方程组解法的理论,编写了二维和三维的
随着我国国民经济的迅速发展,社会用电需求的不断增加,高压输电线路发生各种故障的几率也随之增大,这给整个电力系统的正常运行和可靠供电带来了严重的危害。因此研究如何能在广
结构可以通过建筑物单元的能量耗散来抵御地震作用.因此,对于结构的某些单元,允许存在一定的变形与破坏.现在有一种迅速被全世界所接受的减震、抗震方法,那就是基础隔震方案.
该文用有限元方法全面深入地研究了复合材料层合结构的一些非线性热相关问题,其中包括稳态和瞬态热传导、非线性热静动态响应、热环境下大振幅自由振动、热屈曲和热后屈曲以
当今世界,资源枯竭、能源消耗过快与环境的严重破坏是人类所要面对的日益严重的国际问题,经济的快速发展对能源的消耗使得这些问题更加严峻。火力发电在当今的世界范围内仍然是
该文采用二维浅水波方程理论把行驶船舶对水面的作用力直接加在二维浅水波方程组,以此来计算水面三维波高和深度平均的水平流速.对于二维浅水波方程,运用Galerkin加权余量法
该文针对这一断裂力学的薄弱环节,以实验与有限元计算相结合,宏观与微观相结合的方法进行了探索研究.文中首先发展了Ⅱ型裂纹动态应力强度因子的实验测定方法.一方面改进Hopk
该文首先对大柔度相间合成绝缘间隔棒受压稳定性态进行了线性和非线性分析,论证了相间合成绝缘间隔棒的后屈曲性态是稳定的,导出了前屈曲和后屈曲状态下载荷-位移关系曲线.其