针对 Delft 3D FM Suite 中二维浅水方程求解的k与ε提取难题的回应

自然选择号，前进四

已于 2024-05-03 21:32:59 修改

阅读量443

点赞数 6

文章标签： 3d

于 2024-05-03 21:21:40 首次发布

本文链接：https://blog.csdn.net/weixin_59408766/article/details/138423705

版权

对于浅水方程其实有诸多假设，比如（1）忽略垂向及速度、流速等水力参数的垂向变化（2）沿水深方向压力分布符合静水压强分布等等。对于二维浅水方程中，还需要沿水深积分取平均值：忽略风应力以及柯氏力的二维浅水运动方程为：

$\frac{\partial H}{\partial t}+\frac{\partial HU}{\partial x} +\frac{\partial HV}{\partial y}=0$

$\frac{\partial HU}{\partial t}+\frac{\partial HUU}{\partial x} +\frac{\partial HUV}{\partial y}=-gH\frac{\partial z_{s} }{\partial x} -gn^2\frac{U\sqrt{U^2+V^2} }{H^{1/3}}+\nu _{t}( \frac{\partial^2}{\partial^2 x}HU+\frac{\partial^2}{\partial^2 y }HU )$

$\frac{\partial HV}{\partial t}+\frac{\partial HUV}{\partial x} +\frac{\partial HVV}{\partial y}=-gH\frac{\partial z_{s} }{\partial y} -gn^2\frac{V\sqrt{U^2+V^2} }{H^{1/3}}+\nu _{t}( \frac{\partial^2}{\partial^2 x}HV+\frac{\partial^2}{\partial^2 y }HV )$

一般而言，我们认为求解紊动涡粘度 $\nu _{t}$ 是可以通过k-ε方程求解，但是在Delft3D FM Suite使用二维浅水方程求解水流情况时， $\nu _{t}$ 是利用背景值（ambient/background coefficients）求得，也就是在mdu文件中<Physical Parameters>设置，或者在GUI中General-Physical Parameters中设置统一值，或者通过spatial editor设置空间分布值。但是无论怎么设置，都是背景值，并不是通过k-ε求解，最终方程中的 $\nu _{t}$ 也采用的是背景值。这就导致先前为什么在二维浅水方程模拟水流无法提取k、ε参数，因为二维浅水模拟 $\nu _{t}$ 并不是通过k-ε求解，而是用户自定义值。（详见D-FLOW用户手册 8.3.5）

水流中扩散包括三方面：（1）分子扩散；（2）紊动扩散（手册中3D-turbulence）；（3）离散（手册中2D-turbulence）。在二维模拟中，水平速度沿垂向分布无法求得，因此无法求解离散系数；紊动扩散是sub-grid尺度，因此也无法通过雷诺平均求解。这也就说明了为什么在Delft3D FM Suite使用二维浅水方程模拟水流，即使能够导出水平涡粘度 $\nu _{H}$ (同样由三部分组成：分子、紊动、离散)以及水平扩散系数 $D _{H}$ ，但是为设定的背景值，因为真实值难以求解。

图源自（D-FLOW User Manual, section 8.3.5）

但是在3D模拟中会通过k-ε或者k- $\tau$ 模型求解sub-grid涡结构，以及水平速度沿水深的分布导致的离散，因此在3D模拟中是可以提取k、ε、 $\nu _{t}$ 、水平涡粘度 $\nu _{H}$ 以及水平扩散系数 $D _{H}$ 。

而对于背景值的设定会不会影响3D模拟结果，通过我对提取到参数的观察对比，是没有影响的。对于用户手册中的第九章，物质输移扩散系数=分子扩散+背景扩散系数+紊动涡粘度/Prandtl-Schmidt number 也是针对二维模拟的，并不是针对三维模拟。所以对于3D水流模拟中物质输移，背景值会不会产生影响，目前难以给出确切的答复。

本文感谢博主@小归海爱吃小海龟的提供的思路。

自然选择号，前进四

关注

6
点赞
踩
7

收藏

觉得还不错? 一键收藏
0
评论
针对 Delft 3D FM Suite 中二维浅水方程求解的k与ε提取难题的回应

（3）离散（手册中2D-turbulence）。是利用背景值（ambient/background coefficients）求得，也就是在mdu文件中设置，或者在GUI中General-Physical Parameters中设置统一值，或者通过spatial editor设置空间分布值。但是无论怎么设置，都是背景值，并不是通过k-ε求解，最终方程中的。模型求解sub-grid涡结构，以及水平速度沿水深的分布导致的离散，因此在3D模拟中是可以提取k、ε、
复制链接

扫一扫