基于迎风型紧致格式的开闸式浊流在长直水槽中的数值模拟

发布时间：2020-10-15 01:12

　　 本文将一般多相流的思想借鉴到含分散质多相流(本文指浊流)的模拟中,用改进的低色散的迎风紧致算法对长直水槽中的开闸式浊流进行三维数值模拟。首先,对泥沙输运方程中对流项的求解算法进行优化,构造了保色散的迎风结合紧致差分格式(Upwind Dispersion Relation Preserving Combined Compact Difference,UDRPCCD),同时,引入水平集合法的距离函数思想对多相流浓度界面进行平滑处理并用其描述浓度场。改进的迎风DRPCCD算法和其与水平集函数的结合使本文的多相流模型适用于大梯度流场的同时实现对方程中非线性对流项的高精度求解。之后,对改进的低色散算法进行数学分析,包括精度比较、色散性与耗散性比较、各向异性分析及稳定性分析。接下来,将本模型运用在平坡及不同尺寸障碍物存在下的平坡地形中,对浊流演变的外部特征(头部位置、浓度剖面与三维形态)及内部规律(能量转变与速度场)进行模拟和研究。数值结果表明,浊流在平坡上的头部速度呈现三段式趋势,经历加速、匀速和减速阶段,模拟结果与现有的研究理论及实验结果基本一致;泥沙的悬浮与沉降结果表明粗泥沙颗粒具有更快的沉降速度,三维结构图则进一步表明细泥沙颗粒的加入对整体浊流形态的决定性影响;能量方面,可以看到本模型较好的能量守恒性质,各能量分量(动能、势能、粘性损失和颗粒沉降损失)的变化趋势与Nasr-Azadani等的三维直接模拟的结果一致-最初势能快速转化为动能,泥沙沉降损失由初始大于粘性损失逐渐减小;同时,从数模结果中也可观察到出现在物模试验中的其它特性,如异重流头部高于临界点的结构特点,由以上结果验证了本浊流数值模型的合理性和可行性。最后,考虑到浊流通常出现的海底起伏地形,在平坡模型中加入了两种尺寸的高斯型障碍物并进行模拟,结果与Meiburg团队的研究一致,从而进一步验证了本模型在浊流问题上的适用性和可延展性。
【学位单位】：浙江大学
【学位级别】：硕士
【学位年份】：2017
【中图分类】：TV149.2
【部分图文】：

需要考虑密度变化外，其他项密度的变化均可以忽略，从而使我们能够在密度变??化只引起微小修正时把它忽略，同时在密度变化有动力学意义时又能保持它［５９］。??本文模型的初始设置见下图（图２．１），采取全高度的瞬时开闸设置以产生??重力流。打开闸门后，两种密度不同的流体在压力作用下相互侵入，其中，重流??体中含有两种粗细不同的颗粒，颗粒与流体的密度差为定义约化重力＃??图２．１开闸式双组分浊流的初始设置图??对于物理问题中的参数，有两种方法进行处理。一种是直接将物理问題的变??量本身带入到数值模型中进行计算，它的好处在于，在输出结果中可以直接看到??目标量纲的结果。但对于不同单位的输入或输出量，需要人工进行相应量纲换算，??且最终不容易对比不同量纲下的结果。另一种方法是引入各变量的一些特征值，??１１??

需要考虑密度变化外，其他项密度的变化均可以忽略，从而使我们能够在密度变??化只引起微小修正时把它忽略，同时在密度变化有动力学意义时又能保持它［５９］。??本文模型的初始设置见下图（图２．１），采取全高度的瞬时开闸设置以产生??重力流。打开闸门后，两种密度不同的流体在压力作用下相互侵入，其中，重流??体中含有两种粗细不同的颗粒，颗粒与流体的密度差为定义约化重力＃??图２．１开闸式双组分浊流的初始设置图??对于物理问题中的参数，有两种方法进行处理。一种是直接将物理问題的变??量本身带入到数值模型中进行计算，它的好处在于，在输出结果中可以直接看到??目标量纲的结果。但对于不同单位的输入或输出量，需要人工进行相应量纲换算，??且最终不容易对比不同量纲下的结果。另一种方法是引入各变量的一些特征值，??１１??

图２．３泥沙输运方程的推导示意图??泥沙的质量守恒，考虑量纲，泥沙的浓＾＋（ｕ＋ｗ，ｅｇ）－Ｖｃ，?＝ＤｐＶ２ｃ，，?ｉ＝＼，?２浓度，ｉ表示本文模拟的两种不同粒径ｒ从，心，７＼从．分别指玻尔兹曼常数（Ｂ）的粘度。同上节的方法，对泥沙输侧同时除以一。可以得到，原始ＨＩ２ｌｕｂ??由此，无量纲化后系数依次为ｕ久丄ｌｕ＼的泥沙输运方程可具体表示如下。??Ｑｑ?１??－－—２
【相似文献】

相关硕士学位论文 前1条

1 王妲;基于迎风型紧致格式的开闸式浊流在长直水槽中的数值模拟[D];浙江大学;2017年

本文编号：2841465