海南某滨海核电厂温排水三维数值模拟-广东水利水电科学研究院.PDF

海南某滨海核电厂温排水三维数值模拟 倪培桐，张晓艳 (广东省水利水电科学研究院，广州，510610) 摘 要：目前工程上广泛应用温排水数值模型多为平面二维模型，但平面二维数值模拟不能完全反映 温排水的近区水力、热力效应，同时也不能模拟工程海域的温盐分层环境，制约着平面二维温排水数 学模型的成果应用。本文在海南某滨海核电厂温排水研究基础上建立了温排水三维数值模型，并成功 应用于工程方案比选。 关键词：温排水；三维数值模型 前言 温排水运动是复杂的水力、热力现象，完善的数学模型需要模拟近、远区紊动要素的影响， 反映水体水平、垂向出流掺混，热水回归，反映水面散热影响、水体浮力作用等要素。但由于对 环境水体近区、远区的水力、热力特性认识局限，以及数值模拟计算工作量、边界处理、参数选 取等原因，三维数值模型仍未在工程实践中广泛应用。目前工程上广泛应用温排水数值模型多为 平面二维模型，但平面二维数值模拟不能完全反映温排水的近区水力、热力效应，同时也不能模 拟工程海域的温盐分层环境，制约着平面二维温排水数学模型的成果应用。鉴于此，我国多家单 位在应用平面二维数值模型对工程问题进行研究的同时，开展了温排放三维数学模型研究 [1]2[3][4][5] 。本文在对我国海南某滨海核电厂相关课题研究的基础上建立温排水三维数值模型。由于 三维温排水数学模型研究仍处于发展阶段，且计算工作量大、计算时间长，本次三维温排水数值 模拟重点研究温排水自身作用下的水力、热力特性对不同方案取水温升的影响，为二维数学模型 方案比选工作提供支持，为冷却水工程方案优化布置服务。 1 温排水三维数学模型 1.1 控制方程 （1）潮流方程 ∂u ∂v ∂w + + S (1) ∂x ∂y ∂z ∂u ∂u2 ∂vu ∂wu ∂η 1 ∂P + + + fv =−g − a − ∂t ∂x ∂y ∂z ∂x ρ0 ∂x (2) g η ∂ρ 1 ∂s ∂s  ∂ ∂u xx xy   dz −  + +F + υ +u S ρ0 ∫z ∂x ρ0h  ∂x ∂y  u ∂z  t ∂z  s ∂v ∂uv ∂v2 ∂wv ∂η 1 ∂P + + + =−fu −g − a − ∂t ∂x ∂y ∂z ∂y ρ0 ∂y (3) g η ∂ρ 1 ∂s ∂s  ∂ ∂v yx yy   dz −  + +F + υ +v S ρ0 ∫z ∂y ρ0h  ∂x ∂y  v ∂z  t ∂z  s 55 其中t 是时间，x , y ,z 为笛卡尔坐标，η是水位，d 是静水深，h d +η 是总水深，u(x , y , z ) 、 v(x, y, z) 、w(x, y, z) 是水流速度，f 2Ωsinφ是科势力因子，g 是重力加速度，ρ(x , y )