基于数值模拟的防护林规划设计探讨..doc

基于数值模拟的防护林规划设计探讨 摘 要：采用CFD方法，对崇明岛瀛东村鱼塘防风林的规划进行了数值模拟分析。在二维计算流体力学模拟中，采用标准k–ε湍流模型。分别模拟防风林在固定林带宽度下，有叶期和无叶期以及护效果变化。分析了绕林流场,比较了防护距离的差异。结果显示，无叶期的防护距离L50为85m，覆盖了鱼塘的65%。而有叶期可以完整保护鱼塘。若直接栽种树苗有5年左右的培育期，防护距离L20可以覆盖鱼塘，L50也可以达到鱼塘的35%左右。 关键词：防风林；数值模拟；计算流体力学（CFD）；疏透度 中图分类号: 　　　　　　文献标识码: A 引言 人类利用防风林改善局部风环境已经有几百年的历史了。过去人们常常进行实地测试来评价防风林的效果，再对防风林进行修缮[1-2]。后来引入了风洞模型实验的方式。但这些方式成本较高，而且修缮具有滞后性。近些年来，随着计算机技术的发展，许多科研工作者将计算流体力学（简称CFD）技术应用于防风领域来评价防风林的设计效果。它可大大减少人力物力消耗，及时反馈设计效果。结果表明,防风林对风的阻力随林带宽增加而增加,但是整体阻尼系数略微降低。防风林带宽度明显影响最小风速值及其位置[8]。目前国内对林业领域内数值模拟的应用还比较少。本文应用CFD技术崇明岛瀛东村的防护林布局进行分析。 崇明县陈家镇瀛东村，位于崇明岛的最东南端，由滩涂围垦而成。崇明自古有“瀛洲”的美名，瀛东村成为“瀛洲”最东、最早迎来旭日东升的村庄。从1985年开始，在村党支部书记陆文忠的率领下，全村党员和群众历经艰辛，围垦荒滩4200亩，建立起一个以淡水养殖为主的新村庄，1989年定名为瀛东村。长三角地区是我国重要的经济发展区，也是人口、城市最密集地区，其社会和财产易损性因此也较大。生态风险增加，这对长三角地区区域防护林体系的构建标准提出了更高要求[9]。 瀛东村属于亚热带海洋性季风气候。风向以西北风和东南风为主。本次规划设计中在鱼塘的周边布设防风林，保护鱼塘。防护林位置见图1中部分。 Wang H和Takle(1995年)等用数值模型得到了一系列防护林流动状态，并研究了它们随孔隙度的变化。模拟与现场测量基本一致[3]。Wilson等（2004）将户外测量结果和数值模拟结果进行了对照，发现使用雷诺模型模拟的数据与实验数据有较好的一致性，只有10%-15%的误差。应用有限差分格式求解方程时，高阶迎风格式的选用需仔细考虑[4]。Wei Li（2007）进行了在包含防风林和建筑物的城市开放空间中风环境数值模拟和实验研究。使用KE两方程模型也得出了精确的数据[5]。Santiago（2007年）使用了standard k–ε, RNG k–ε和realizable k–ε湍流模型模拟了防护林的防风效果，并将数据与风洞实验的结果对比，意图求出最佳孔隙度[6]。Bourdin, P.等（2008年）研究了数值模拟在防风林模拟方面的适用性，使用realizable k–ε模型模拟，设置二维三维情况下的工况，对比了各个情况下的模拟数据与实验数据[7]。Mulati Yusaiyin（2009年）通过在日本扎晃人学水压力工程实验室开展的埃菲尔型非循环风洞试验,研究了防风林阻尼变化及不同林带宽的整体阻力系数，阐明在迎风向上防风林带宽对风速减弱的作用。以二维雷诺时均N-S方程k-ε闭合模型研究了不同防风林带宽的流场变化。k-ε湍流模型进行计算，该方法应用较广也以得到了实验验证。 2.1求解区域 求解设置两排防风林的二维区域。沿水塘西北到东南前后两排林带，林带高H，宽取0.4H，水平方向上第一层林带带前-20H，第二层林带带后40H，垂直方向上0～10H的二维区域为求解空间。详见图1。数值计算过程中H取为单位长度。H取值为7m。 图2 几何模型示意图 2.3计算模型和控制方程 湍流模型采用标准的模型进行模拟。这种模型的模拟效果已经得到一些学者的验证[4、6、8]，近壁区域采用壁面函数法来处理。 连续方程、动量方程以及描述湍流的两方程模型都可以采用通用形式的控制方程表达： （1） 方程中分别代表u，v，w及湍流动能和湍流耗散率等变量，方程中代表扩散系数，方程最右端的代表源项。 连续方程 1 0 0 X方向动量方程 u Y方向动量方程 v 湍动能 k 耗散率 ε 其中 （2） （3） 在标准模型中，根据Launder等的推荐值及后来的实验验证[12]，模型中各个系数值分别为：方程中各个系数值 1.44 1.92 0