绿化对室外热环境影响模拟研究.docVIP

喷泉对周围热环境影响的模拟研究 林波荣1，李晓锋1，朱颖心1 改善热环境、提高热舒适有重要作用。本文通过 关键字：喷泉，CFD，室外热环境Nobuya Nishimura 等人[1]对公园中喷泉的降温效果进行了实测。实验结果证明了水景的降温效果，同时通过风洞实验证实了可用通过人工方法来控制水景的降温效果。但是这一研究没有对喷泉建立理论分析模型，而只是通过测量定性分析了喷泉周围热环境的特点，无法进行推广和应用。 因此，本文通过喷泉模型水珠 (1) 其中为水珠位置矢量，为水珠速度矢量，水珠的速度通过求解水珠的动量方程得到。 计算水珠的速度的动量方程如下： (2) 拉格朗日时间步长为： (3) 在拉格朗日模型中，湍流的扩散效应是通过随机模型进行模拟的。假定湍流是各向同性的，流体速度由时均速度和波动速度组成： (4) 由连续态的欧拉方程计算得到，为波动部分。用来分析流体速度的波动速度的湍流信息是来自流体方程中的湍动能。假设湍流是各向同性的，其速度的波动部分满足Gauss分布，标准偏差为： (5) 在任何位置波动部分的速度是通过随机采样得到的。 数值求解水珠运动轨迹的求解步骤如下：（1）计算拉格朗日时间步长；（2）水珠按照原来的速度在该时间步长内运动到新的位置；（3）计算得到水珠在新地点的物性（速度、温度等）；（4）计算相间源项。 第四步完成后重新进行第一步计算，直到水珠蒸发完或者运动出模拟区域。 1.2水珠与空气的热质交换过程水珠与空气的热质交换过程可以视为一个准稳态的蒸发过程，也就是在一定的时间间隔内，周围空气的物性不发生变化，同时水珠的物性也没有改变。考虑球对称的蒸发水珠模型，水珠和空气的热质交换过程参见文献[3]。 1.3 PSI-CELL模型 由于空气与水珠之间存在强烈的耦合关系，因此需要采用合适的模型来处理水珠在空气中温度、质量以及运动速度和位置的耦合关系。本文采用的是Particle-Source-In Cell（PSI-CELL）模型[5]，该模型的主要思想是将水珠对空气的影响考虑成空气的质量、动量和能量源项。整个连续态流场通过欧拉方法求解离散的差分方程得到。水珠的轨迹、尺寸和温度变化过程通过对水珠在流场中运动方程进行积分得到，采用粒－气热质交换率，计算采用拉格朗日方法。 计算过程如图1所示。首先在不考虑水珠影响条件下计算初始的空气流场，然后根据流场计算水珠的轨迹、直径和温度变化，由此可以得到所有水珠在每个网格产生的质量、动量和能量源项，然后将这些源项带入空气方程中，再次计算出新的空气流场，新的流场用于计算新的水珠轨迹、直径和温度变化，因而既包括空气对水珠的影响，也包括水珠对空气的影响，由此反复计算得到最终结果。 连续态的运动采用欧拉方程对进行计算，公式如下： (6) 蒸发量、动量源、能量源的计算公式如下： (7) (8) (9) 采用文献[4～6]的方法给出的水珠平均直径参数、水珠直径分布函数、粒径分布函数、平均水珠速度的计算方法来描述喷头特性，即喷出水珠的直径大小与分布。 1.4直管射流边界条件的确定 在建立喷泉数值模拟方法之前还需要得到水珠射流边界条件，即喷泉射流出口水珠的粒径分布、质量分布和速度分布等[7]。选择在喷泉中大量使用的圆口压力直管射流进行分析。通过分析直管压力圆口射流的破碎机理和模式，将射流分为一次破碎和二次破碎两个阶段，分别给出了各阶段的粒径分布、SMD计算公式以及速度分布。从而可以得到水珠射流边界条件，然后可采用CFD软件模拟得到不同边界条件下喷泉周围各种参数的分布情况。 二次破碎后SMD计算公式为： (10) (11) 二次破碎以后直径为d的粒子速度为： (12) 一次破碎初始SMD公式为： (13) (14) 其中为径向积分尺度，对于直径为D的完全发展的湍流管流，轴向积分尺度为。综合考虑一次破碎和二次破碎的水珠直径公式，最后粒径的水珠直径SMD公式如下： (15) 利用上述方法，模拟了喷泉，并与喷泉热2 喷泉对周围热环境影响的联合模拟方法 在小区热环境联合模拟中，首先假定空气流场温度均匀，等于来流空气