对于浅水方程其实有诸多假设,比如(1)忽略垂向及速度、流速等水力参数的垂向变化(2)沿水深方向压力分布符合静水压强分布等等。对于二维浅水方程中,还需要沿水深积分取平均值: 忽略风应力以及柯氏力的二维浅水运动方程为:
一般而言,我们认为求解紊动涡粘度是可以通过k-ε方程求解,但是在Delft3D FM Suite使用二维浅水方程求解水流情况时,
是利用背景值(ambient/background coefficients)求得,也就是在mdu文件中<Physical Parameters>设置,或者在GUI中General-Physical Parameters中设置统一值,或者通过spatial editor设置空间分布值。但是无论怎么设置,都是背景值,并不是通过k-ε求解,最终方程中的
也采用的是背景值。这就导致先前为什么在二维浅水方程模拟水流无法提取k、ε参数,因为二维浅水模拟
并不是通过k-ε求解,而是用户自定义值。(详见D-FLOW用户手册 8.3.5)
水流中扩散包括三方面:(1)分子扩散;(2)紊动扩散(手册中3D-turbulence);(3)离散(手册中2D-turbulence)。在二维模拟中,水平速度沿垂向分布无法求得,因此无法求解离散系数;紊动扩散是sub-grid尺度,因此也无法通过雷诺平均求解。这也就说明了为什么在Delft3D FM Suite使用二维浅水方程模拟水流,即使能够导出水平涡粘度(同样由三部分组成:分子、紊动、离散)以及水平扩散系数
,但是为设定的背景值,因为真实值难以求解。
图源自(D-FLOW User Manual, section 8.3.5)
但是在3D模拟中会通过k-ε或者k-模型求解sub-grid涡结构,以及水平速度沿水深的分布导致的离散,因此在3D模拟中是可以提取k、ε、
、水平涡粘度
以及水平扩散系数
。
而对于背景值的设定会不会影响3D模拟结果,通过我对提取到参数的观察对比,是没有影响的。对于用户手册中的第九章,物质输移扩散系数=分子扩散+背景扩散系数+紊动涡粘度/Prandtl-Schmidt number 也是针对二维模拟的,并不是针对三维模拟。所以对于3D水流模拟中物质输移,背景值会不会产生影响,目前难以给出确切的答复。
本文感谢博主@小归海爱吃小海龟的提供的思路。