2018年工程热力学.docVIP

工程热力学 　　　　5.1 引言 　　根据经典的传热概念，自由水面上的传热通常可分为辐射，对流，传导和蒸发。为了对水面热通量进行数值建模，可根据传热的透水能力对其进行分类。蒸发、传导和长波辐射都是只发生在水面的表面传热作用，而短波辐射则是一种穿透性作用，能将其热量在散布在水柱中相当大的范围内。 ELCOM模拟的两个表面质量通量为降雨与蒸发。 　　5.2 长波辐射 　　对于大气和水面的交互作用，我们通常采用以下系数：(1)辐射系数，（2）云量，（3）水面反射率（参照TVA）。文献提供了几种不同的方法来应用这些经验系数。在数值模型DYRESM（Imberger和Patterson）中，长波辐射问题被分为：(1)水体的辐射散发，和（2）到达并穿透水面的大气净辐射。我们令这两种长波辐射分别为Q(emitted) 和 Q(absorbed)。描述这两种长波辐射的方程如下： 　　　　和 　　　　其中C(cloud)为相对云量，T(water)为水面温度，T(air 2)为水面上两米处测得的大气温度，而Rt(lw)为水面对长波辐射的总反射率。大气的辐射系数（emissivity）为： 　　　　长波辐射的净传热（QR）是水体散发量和对大气辐射的吸收量之和： 　　　　联立以上方程组，可得出以下计算长波辐射的公式： 　　　　依照旧版本的DYRESM模型，我们采用以下常量： 　　C? : 0.937 x 10-5 C-2 （Imberger和Patterson，方程9） 　　? ?(water) : 0.96 （Imberger和Patterson，方程8下方） 　　? Rt(lw) : 0.03（Imberger和Patterson，方程9下方） 　　-8-2-4? ? : 5 .669 x 10 W m C （Holman，第387页） ? 　　5.3 蒸发潜热 　　ELCOM模型中的蒸发热通量项（或蒸发潜热通量项）可采用三种方法来模拟，用户可根据已知的输入数据来选择模拟方法。当用户已知大气比湿和水面比湿时，适当的蒸发热通量模型为： 　　　　式中，设定蒸发潜热为DYRESM中采用的值：2.453x106 J/kg （水体温度等于20.4℃）。 　　若仅已知参考基准面上的比湿，我们即可假定水面湿度接近饱和，则蒸发通量模型可写为： 　　　　饱和比湿度由以下公式估算： 　　　　饱和水汽压可由对实测数据和实验数据进行拟和得到的曲线（Gill，1982，第606页）近似估算出： 　　　　上式的计算结果是单位为Pa的饱和水汽压。温度的单位为℃（有效范围为-40℃～+40℃）。联立以上三个方程，即得第二个蒸发热通量模型为： 　　　　如果我们仅有相对湿度值，比湿的估计值为： 　　　　故第三个蒸发热通量模型为： 　　　　因此，方程5.6、5.10和5.12分别为蒸发热通量模型的三种不同实现方法。 　　5.4 显热通量 　　显热通量（如传导、对流作用）所引起的能量传输可模拟为（文献Imberger和Patterson 1981中的方程7）： 　　式中CH为（无量纲的）显热传递的容积系数，Cp为恒压下的比热容，为水面上方某一参考基准面的干球温度。在该方程中仍保持传输到水体以外的能量为负值的惯例。Imberger和Patterson并不清楚以上方程中的比热项和密度项应取空气的值还是水体的值。Fischer等人采用的是空气的比热和密度（方程6.20），Gill在方程2.4.5中也采用的是空气的值。 　　根据Imberger和Patterson的论述，在ELCOM中，取CH的值为1.4x10-3，则可实现方程5.13。在水面附近常温下，空气比热容约为1003 J/kg/℃。 　　5.5总水面传热 　　水面传热可模拟为净辐射传热、由蒸发引起的潜热传输、由传导和对流引起的显热传输之和（即除了短波辐射外所有通过自由水面的传热）： 　　　　可假定水面能量传输发生在厚度为0.6～1.0米的层中，且能量（沿深度）服从指数衰减，则： 　　式中QS(z)为在高度z处（以能确保在向上方向上z为正值的坐标基线为测量起点）尚未被吸收的表面热能，QS(S)为穿透水面的表面传热，?s为表面传热的总体消光系数，S为水面的高度（测量起点同z）。如果我们假设90% 　　的表面 　　传热被水面至深度0.6米范围内的水域所吸收，则在深度1.0米处，尚未被吸收的表面热的输入量（增量或损量）为表面热能总输入量的2%。 　　表面热能的模拟是通过对整个水柱进行积分来实现的，表面热能决定着每一层的吸收量。若定义第k层顶部处的有效表面传热能为QS(k)，定义最上层为第n层，则方程5.15可写为： 　　　　其中?z(i)为第i层的厚度，??s为总体消光系数（设为3.84，使得98%的表面传热能量被水面至深度1m范围内的