防风网设置的数值模拟-中国力学学会.docVIP

防风网设置的数值模拟 陈素琴1李萍2詹水芬3 （1同济大学数学系 2上海超级计算中心 3 交通部天津水运工程科学研究所） 1、引言 堆料场在有风的情况下经常会产生大量扬尘，既造成了原料的损失，也造成了严重的空气污染。防风网是大型堆料场防尘的重要措施，在国外已普遍应用。随着环境保护的日益重视，防风网在国内也逐渐开始应用。 防风网的基本原理是在大型堆料场周围设置多孔透风屏障，降低来流风速并减弱大气湍流中的漩涡结构，达到抑制扬尘的目的。 目前，风洞模拟实验是研究防风网设置的有效方法，但是实验耗费大量的人力、物力和财力。近年来，基于流体力学的数值模拟成了研究该类问题的重要手段，相比较风洞实验而言，该方法具有费用低、可重复性强等优点；但数值计算结果的有效性还需经过实验的验证。 2、研究对象及目的 研究对象：某岛港煤及矿石堆场防风网设置作为研究对象。该港是中国最大的煤炭输出港之一，煤堆场现有的防尘设施不足以确保实现大气环保指标，需要建设防风网来满足需求。 研究目的：通过数值模拟计算取得相关技术参数，并与风洞实验进行比较，科学地评估防风网设置的数值模拟和风洞实验结果的以一致性。为将来用数值模拟方法取代风洞实验来确定不同防风网布局方案的降尘、抑尘效果提供科学依据，保证工程设计的快速顺利进行。 3、模型建立和仿真结果 3.1数值方法 该模拟使用Fluent软件，Fluent是一个用于模拟和分析在复杂几何区域内的流体流动与热交换问题的专用计算流体动力学（computational fluid dynamics, CFD）软件。 Fluent并行算法介绍：Fluent在高性能计算技术方面，Fluent提供从模型设置，模拟计算结果到结果后处理等一整套完整的并行计算环境。在网格处理技术方面，Fluent在并行环境下可以实现网格的自动分区并在计算过程中进行动态计算载荷平衡。在求解方法上所有网格模式，求解器算法和各种物理模型均可以实现并行。 3.2 模型建立 模型介绍：堆料场原型南北向长1945米，东西向宽1340米；在堆料场内按一定的规律放置26排煤堆，每排煤堆包括6个或7个棱台型煤堆；这些棱台型煤堆的下底的南北向长51米，东西向宽154米，上底的南北向长10米，东西向宽113米，煤堆的高度为17米。为模拟实际的防风效果，将堆场周围的建筑物也作了考虑，作为模型的一部分。计算范围为底半径为4500米，高为150米的16边型棱柱。 几何建模和网格生成使用Gambit软件，由于计算区域较大，网格生成采用分区方法，在煤堆和建筑物附近采用边界层网格，最细网格尺寸为0.1，这样保证了计算所需要的最小网格，同时控制了网格总数，保证计算的可靠性及计算的效率。 3.3仿真模拟 在数值模拟中，由于防风网的防风范围较大，实际的防风网的弯曲尺寸在数值的计算区域中相对较小，在全范围防风网的数值模拟中，由于网格布置的要求，将防风网处理成带有一定孔隙率的平板网。为研究这种处理的有效性，首先对实际的防风网和平板网进行了防风效果的比较，比较结果如下： 当来流速度为5m/s时，平板网后的速度降为1.73m/s，真实形状网后的速度降为1.69m/s，其它的流场定性形态基本一致。基于此，在以下的数值模拟中，将防风网处理成平板是合理的近似。 在数值模拟中，使用膜边界条件来代替防风网，膜的开孔率的参数的设置应保证与风洞实验所用的防风网的防风效果一致，为此对不同开孔率的膜的防风效果进行了研究，并将防风结果与风洞测试的结果进行了对比，结果如下： Fig. 1 Shelter effect of porous wind fences with different porosity 从图形上观察， 开孔率为47%的膜边界与实验结果吻合较好，故在本数值模拟过程中所使用的开孔率取为47%。 由于流动为不可压流，在此采用Fluent中分离式求解器（segregated solver）的定常计算，控制方程中的扩散项及对流项都采用二阶迎风格式进行离散，根据收敛的状况修改方程的松弛因子，边界条件根据实测结果给定。 计算一开始采用1500万网格单元，经过加密使得计算结果不依赖网格尺寸，最终计算采用的网格单元为1600万单元。采用16CPU并行计算12个小时完成一个工况，本项目要对16个不同方向下有无防风网分别计算，共计32个工况，每次可以同时提交4个作业，整个计算过程不到一个星期。 4、结果分析与比较 本文对上述计算模型的16个风向下有无防风网的32个工况采用RNG k-ε模型进行了计算。下面仅对北风向的RNG k-ε模型的结果进行分析和比较。 图2是煤堆的摆放方式，编号从左到右，最左边的为第1行，共26行。 图3是一个煤堆表面的测点位置，共有11个测点，其中的尺寸与实际尺寸之比为1：200。将每个煤堆