8、TK-2DC软件原理（可下载word文件形式）.doc

8、TK-2DC软件 8.1 模型网格 模型采用贴体正交曲线网格系统来克服边界复杂及计算域尺度悬殊所引起的困难。 (a)物理域 (b)计算域 贴体正交曲线网格坐标变换 所谓贴体正交曲线网格，就是将沿着边界的变化趋势，借助函数变换，将复杂的不规则物理求解域转化为规则的计算域，将复杂区域上的定解问题转化为规则变换区域上的定解问题，使边界条件得到准确满足。在此规则变换域上，无论采用有限差分法、分析法还是有限体积法等进行离散计算，其网格的选取以及边界条件的满足都将很方便，计算精度得以提高。物理域边界ABCD是不规则的，通过坐标变换后，转换成了规则的边界的计算域。边界拟合坐标中坐标转换实际上就是已知计算平面内的坐标，通过计算平面和物理平面之间的函数表达式，找出与已知计算平面相对应的物理平面上的坐标。对生成的边界拟合坐标有以下几个要求：（1）物理平面上的节点应与计算平面上的节点一一对应，同一族中的曲线不能相交，不同族中的两曲线仅能相交一次；（2）在边界拟合坐标中的每一个节点应当是一系列的曲线坐标的交点，而不是一群三角元素的顶点或一个无序的点群，以便设计有效、经济的算法及程序，要做到这一点，只要在计算平面中采用矩形网格即可。物理平面求解区域内部的网格疏密程度要易于控制。在边界拟合坐标的边界上，网格线最好与边界线正交或接近正交，以便于边界条件的离散化。 .2 控制方程 坐标转换关系的基本方程 平面直角坐标系下的任意形状的区域，通过边界贴体坐标，可转化为新坐标系下的规则区域。新旧坐标采用：函数关系联系起来。假定变换关系满足Poisson方程与Dirichllet边界条件。 Poisson方程： （4-1） （4-2） Dirichllet边界条件： 新旧坐标的变换关系为： （4-3） （4-4） 式中，P，Q是与有关的某一函数，反映了（）平面上等值线在（）平面上的疏密程度，适当选择P，Q函数，可使坐标变换为正交变换。根据水流势函数与流函数的性质及水流等势线与等流线的正交性，可导出生成正交曲线网格的转换方程： （4-5） 式中，；；；，，。 拟合坐标系下水流运动控制方程 水流连续方程 （4-6） 方向动量方程 （4-7） 方向动量方程 （4-8） 式中，、分别表示正交曲线坐标系中二个正交曲线坐标；u、v分别表示沿、方向的流速；ｈ表示水深；H表示水位；、表示正交曲线坐标系中的拉梅系数： ，； 式中，、、、表示紊动应力： 式中，表示紊动粘性系数，一般情况下，，=0.5~1.0，表示摩阻流速；对于不规则岸边、整治建筑物、桥墩作用引起的回流，可采用k-ε紊流模型，k表示紊动动能，表示紊动动能耗散率。正交曲线坐标系下， 紊动动能输运方程 （4-9） 紊动动能耗散率输运方程 （4-10） 紊动动能产生项 式中，、表示因床底切应力所引起的紊动效应，它们与摩阻流速间的关系为：；；；；；、、、、——经验常数，采用Rodi建议的值：，，，，，。 推移质不平衡输移方程 推移质输沙守恒示意图 设推移层的厚度为，为系数，为床沙粒径；河宽为；推移质平均运动速度为；推移层在河段内平均含沙量为，则推移层的平均输沙率为 （4-11） 进口断面单位时间的来沙量为 （4-12） 出口断面单位时间的排沙量为 （4-13） 下界面落至河床的单位时间落淤量为 （4-14） 式中，为下界面的含沙量与平均含沙量的比值，为泥沙沉速。 下界面单位时间从河床上冲起的沙量为 （4-15） 式中，为推移质饱和含沙量，为推移质有效输沙率，为饱和条件下，推移层下界面含沙量与平均含沙量的比值。考虑到推移层厚度很小，在下文里假定。 下界面净流入沙量为 (4-16