水利工程论文-异重流潜入运动的剖面二维数值模拟.doc

《水利工程论文-异重流潜入运动的剖面二维数值模拟.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、水利工程论文-异重流潜入运动的剖面二维数值模拟摘要：本文采用变密度流基本方程和混合有限分析法，求解了突扩边界下的异重流潜入运动。通过对数值计算结果分析，讨论了进口密度弗汝德数Fre对潜入运动的影响，加深了对异重流潜入规律的了解。关键词：异重流潜入运动数值计算混合有限分析法 1前言　　异重流是自然界中常见的一种流动现象，如水库中的浑水异重流，河口地区出现的盐水楔，工业中遇到的温差异重流，底部突然下沉的小型沉沙池中的泥沙异重流等，在水库、船闸引航道、河港、河口、沉沙池等实际工程中，浑水异重流的淤积问题往往是必须考虑的问题[1]。由于异重流潜入发生处的边界条件复杂，具有与一般明渠流不同的水流泥沙运

2、动规律，且潜入处的水力泥沙因素与异重流的发生、发展及运行状况密切相关[2]，历来广为研究者们所关注。用一维的方法，再加上一些经验关系解异重流运动，无疑具有其实际意义，但异重流运动具有较强的三维性，不是一维流动理论基础所能较好解决的，并且在用一维方法求解异重流运动时，常常碰到诸如上下水层交界面位置及交界面上水体之间的动量和质量交换系数等难以确定的困难。近年来，逐步采用数值模拟技术对异重流潜入运动进行研究，如YoonT.H。采用流函数涡量法对异重流的潜入运动进行了模拟和研究[3]；方春明、韩其为用立面二维对异重流潜入流动进行数值模拟，探讨了影响异重流潜入条件的因素[4]。但由于问题本身的复杂性，

3、许多问题仍有待进一步研究。　　考虑到异重流潜入运动的过程中，不仅具有初始动量，而且还受到因密度差引起的浮力作用，与静止环境中的浮力射流运动有相似之处，不同之处在于异重流潜入发生处的边界条件复杂，因此，本文将异重流的潜入运动视为一种具有复杂边界的浮力射流运动，从NavierStokes方程出发，引入满足布辛涅斯克假定的一组变密度方程，并采用混合有限分析法对该运动进行一系列的数值模拟，通过对数值计算结果进行讨论和分析，进一步加深对该流动的认识。2数学模型图1二维异重流潜入运动的示意图Schematicdiagramof2Ddensitycurrentplungingflow　　图1所示：密度为(

4、ρ0+Δρ)的浑水从水平河道流入密度为ρ0的清水矩形水库中，由于密度分布不均匀，流场中就出现了浮力，在几何形状突然扩大的流道中，形成了有浮力回流的流动。图中x、z分别表示水平和竖直方向，进口浑水流速u0满足抛物线分布。2.1控制方程　　该模型遵循以下假定：忽略泥沙与水流之间的速度差异；根据布辛涅斯克近似，仅重力项计及密度变化；认为泥沙在水平和竖直方向上扩散系数相同；仅对异重流潜入段进行模拟，此时浮力作用较大，且计算区间不太长，因此在这里只考虑异重流的对流和扩散，不考虑其沿程沉降[4]。　　为使计算和讨论具有一般性，以水库初始水深H、进口平均流速u0为特征尺度，引入无量纲量X=x/H,Z=z/

5、H,U=u/u0,W=w/u0,T=tu0/H,P*=p/ρu20　　则控制方程的无量纲形式为　　连续方程 (1)　　运动方程 (2)(3)　　浓度方程 (4) 式中Γ=ρ-ρ0/ρ0为浑水的有效密度，Re=U0H/v为雷诺数，Fr=U0/为弗汝德数，Sc=v/D为斯密特数，v、D分别为流体的运动粘性系数和泥沙扩散系数。2.2初始和边界条件2.2.1初始条件　　初始时刻，水库水体静止，进口断面的浑水流速满足抛物线分布，且认为此断面的密度分布均匀；　　初始压力按静水压力给出：P(X,Z)=1-Z。2.2.2边界条件　　研究的流动区域共有进口、出口、自由表面、上游固壁及底部五个不同的边界，如图1

6、所示。　　进口：浑水流速满足抛物线分布，断面密度均匀分布，即U=7.5(6-5Z)(Z-H-H0/H),W=0,Γ=Δρ/ρ0;　　出口：认为流动已充分发展，则;　　壁面：给定无滑移边界条件，即U=0,W=0,(上游)或(底部)；　　自由表面：对自由面作刚盖假定，并采用对称平面条件，即。3计算方法3.1求解区域及网格划分　　该流动问题是开边界问题中的一种，当出口采用充分发展假定时，求解区域应保持足够的长度，以便在我们感兴趣的区域准确地模拟出流体的运动规律。依文献[5]，本文计算区域长度取为L=20。因水库进口附近的流动是主要的，远离进口的流动是渐变的，故在数值计算中，X方向上采用非均匀网格，

7、进口附近的网格密，越远离进口，网格越疏；在Z方向上，除了0、1附近的网格密些（与边界条件有关），其余部分采用均匀网格。依控制方程及边界区域等特点，采用原始变量法求解为便，故采用交错网格来离散并求解整个流场和密度场[6]。计算时，动量方程和浓度方程的控制精度为10-7，连续方程的控制精度为10-4。3.2计算格式　　采用混合有限分析法的五点格式进行离散。　　在局部单元格上，将非线性方程（2）～（4）线性化，写成