资源描述：
《论文黄河小浪底水库三维水沙数学模型初步研究》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、黄河小浪底水库三维水沙数学模型初步研究一建模和率定摘要本文首先通过转换得到贴体坐标系下三维动量方程和k-£方程；然后，采用基于交错网格的有限体积法对其离散，再采用SIMPLE-C方法对速度场和压力场进行耦合求解，所有代数方程都采用TDMA方法求解；通过求解基于二维水深方程的Poisson方程来确定自由面，固体边界则采用壁面函数法处理；最后，以三维方腔流为算例，对程序的可靠性进行验证。关键词贴体坐标；有限体积法；壁面函数；三维数学模型；方腔流1前言冃前小浪底水库调水调沙调度运用亟待解决的工程泥沙是不同泄水建筑物组合运用条件下坝区输沙流态和出库水流、含沙量、级配过程预测分析等。该问题具有强烈的
2、三维性，一维、平面二维和立面二维模型能反映断面平均、垂线平均或横向平均的水流泥沙运动情况，但不能充分反映具有强烈三维性的小浪底水库坝区水沙运动变化规律。因此，耍全面准确预测分析不同泄水建筑物运用条件下出库水沙过程，满足小浪底水库调水调沙实时调度运用，必须尽快建立小浪底库区三维水沙模型。2控制方程离散与求解2.1贴体坐标下的控制方程通过贴体坐标系统(Body-FittedCoordinates)⑴可以将复杂的物理域的流动问题转换到规则、简单的计算域中进行计算。在计算域中，实现复杂边界条件下的水沙流数值模拟，相应贴体坐标下控制方程可以表示为下列通用的张量形式:O(p0)8(p如)_d(d(/>
3、牴臨jI£严戸二通用控制方程中主要符号的具体形式见表lo表1通用控制方程中各符号的具体形式0s.g』・G10“、Uf"-xrv.k“齐）・（*乡）7*务）1•（护仔♦护许）•（护治•护许）♦（护仔♦册许）］—•冷—宙#・卜［（芳■务）・（*夸）•（&务）］•（扑％•护冷）•（斥忖•册忖）•（&冷•/务）卜—4••哙务畑分方程S0式中：f为时间；号为物理坐标；灯为计算坐标；P为水流密度；竹为流速；0分别代表w,V,w,k.£o2.2控制方程离散采用有限体积法离散控制方程，其中变量布置采用交错网格。取图1所示的控制体积，通用控制方程在节点P的离散方程为：A伸卩=4血+九如+A、如+4©+44+
4、A神+S©采用混合格式⑵离散对流项，该格式虽然只有一阶精度，但是一种绝对稳定的离散格式，精度可以满足工程需要。扩散项和源项采用中心差分离散，最终得到各系数表达式如下：图1控制体积示意图Ap=Ag+人”+Ay++Aj+A”+A；+(Fe-F^Fn-FAF-Fh)其中，F为对流强度，表达式如下:SU,，Fw=(pU)w，Fn=(pV-,九=（pu）$，巧=（pW）,，Fh=（pW）h；Q是扩导系数，其表达式如下：2=（厂汹）+025（「為2）”—0.25（「為2）$+0.25（厂拒/一0.25（「為彳人Ds=（r諾22）$-0.25（厂局必+0.25（『砧2Jw-0・25（「砧23）（+0.
5、25（「祸3）/>。几烧d+0.25（r屈几—0.25（r>J“，+2=(厂諾九一°・25(厂為2)”+0.25(厂為2)$-0.25(r；g/+0.25(r；g/D,=(厂略22)$-°・25(「諾2丄+0.25(r祸仇-0.25(r^23),+0.25(「局3)〃0.25（厂也2）“一°・25（厂砧32）$Db=（「諾33）6-0.25（厂砧31）e+0.25（「屈1）“-°・25（厂后2）”+0・25（厂恋2）$其中，S,nk=(m=l,2,3；k=l,2,3；j=l,2,3)oox.dx.J•»2.4初边界条件2.4.1自由面边界自由表面的运动对非恒定流的泥沙输移起着不可忽略的作用
6、，泥沙的自由表面条件应根据自由表面的上浮通量与沉降通量相平衡给出的，仅当自由表面不变化，或变化十分缓慢的情况下，自由表而法向才与垂向重合。事实上，若床而形态复杂、地形变化较大时，自由表面均会产生运动，其至出现剧烈的波动，导致自由表面发生一定的弯曲和变形；再者，自由表面的运动，会直接引起水流的静水压力的变化，进而对水体的流场，特别是二次流的精细结构产生影响。进行水流的数值计算时，由于自由表面未知，不能直接将自由表面的大气压力作为一个主要的自由水面边界条件。所以,得到的压力场类似于一个封闭管中的压力场，在计算过程中无法直接得到自由水面的位置。近十年来，处理自由表面问题主要有标记结点法、空隙比法
7、和标高函数法等方法，本模型采用自由面位置的Poisson方程求解自由面位置⑵。2.4.2壁面边界采用壁面函数法⑸，即在粘性底层内不布置任何节点，把靠近壁面的第一个节点布置在粘性底层之外的完全湍流区，要求第一个计算节点与壁面间的无因次距离z+=zp^/v（Zp为第一个节点距离壁面的距离）在30〜100Z间。设计算壁面相邻的第一个节点到壁面的无量纲距离Z—pUjv,定义摩阻速度仏二弓咲"2,则计算边界上平行于壁面的流速况满足
