1. 引言
阻沙措施作为机械防沙体系的重要组成部分,在防沙工程中得到广泛应用,其工作原理是在风场内产生一种拖曳力,致使来流动量发生净损失,使气流携载沙能力下降,造成沙粒在阻沙措施附近沉积,从而起到阻沙作用 [1] 。防护效益决定于阻沙措施的设计要素,即形状、高度、疏透度、疏透度分布与方位 [2] ,以及防护区风沙流粒度组成等 [3] 。众多研究者通过测量、分析措施附近的风场,如气流速度、湍流性能、剪切力、压力等参数来阐述阻沙措施的防护效益 [4] [5] [6] [7] ,并用栅栏对风速降低的幅度及其影响范围来判定阻沙效果。然而,疏透度为零的紧密栅栏无论从对风速的衰减幅度还是从形成的低风速区范围角度均优于带有一定疏透度的孔隙栅栏,但工程实践及有关风洞实验研究结果均指向疏透度在0.2~0.4的栅栏阻沙效果最佳 [8] ,单纯从气流场角度不能很好的解释这一现象。地表的风蚀过程本身就是一个非常复杂的现象,对沙粒相的动态测量十分困难,尽管现有的PIV技术已实现对沙粒相的动态测量,但观测的尺度十分有限,单次数据采集仅能反映沙粒运动轨迹的很小一部分。当前,在栅栏阻沙机理方面,尚缺少沙粒运动与风场特性相结合方面的研究。风沙活动本质上属于两相流问题,本文依托现有气固两相流理论,采用欧拉–拉格朗日模型对沙粒流经栅栏的物理过程进行数值模拟,综合流场信息和沙粒的运动轨迹对影响沙粒沉积的主要因素进行了分析,尝试对栅栏的阻沙机制进行更全面的剖析。
2. 研究方法
2.1. 数学模型
本文旨在获取沙粒在流经阻沙措施时的动态信息,采二维欧拉–拉格朗日模型,即气相作为连续相,沙粒相为离散相,模型中仅考虑重力和拖曳力 [9] ,沙粒相视为球型稀相,不考虑沙粒之间的相互作用及沙粒与气流之间的耦合作用。
1) 气相连续方程和动量方程分别为:
(1)
(2)
式中:
、
、
分别为气相密度、速度和压力;
为气相粘性应力张量;
为气相体积分数。
2) 湍流模型采用标准
双方程模型进行描述:
(3)
(4)
(5)
式中:
、
分别为湍动耗散率和湍动粘度;
是由于平均速度梯度引起的湍动能产生;
是由于浮力影响引起的湍动能产生;
为可压缩湍流脉动膨胀对总的耗散率的影响;
、
、
为经验常数;
、
分别为湍动能和湍动耗散率对应的普朗特数;
为速度在
方向的分量;
为热膨胀系数;
为马赫数;
为声数;
为湍动普朗数。
3) 颗粒运动方程:
(6)
(7)
式中:
为颗粒质量;
为颗粒直径;
为颗粒速度;
为重力加速度;
为单颗粒所受的流体粘性阻力;
为流体阻力系数;
为颗粒雷诺数。
2.2. 计算问题描述
选取计算区域为高20 m,长80 m的二维矩形区域,阻沙栅栏位于距进口15 m处,高度H为2 m,分疏透度为0和40%两种,阻沙栅栏结构示意图见图1。计算区域入口、出口、顶部分别采用velocity-inlet、outflow和symmetry边界条件进行描述,计算区域底部和栅栏表面采用wall边界条件行进描述。
入口风速廓线取对数率分布,2 m高度处风速取10、15和20 m∙s−1等3种。沙粒密度为2650 kg∙m−3,分0.1、0.25和0.5 mm等3种粒径,在距计算区域入口10 m处0.1 m~1.0 m高度范围内释放,每种粒径沙粒各释放10个,沙粒初速度方向为水平,速度大小详见表1。
![](Images/Table_Tmp.jpg)
Table 1. Sand grain ejection velocity
表1. 沙粒入射初速度
3. 结果与分析
3.1. 阻沙栅栏结构与风场特征
阻沙措施形成的二次流场是促使沙粒沉积的根源,疏透度与风场关系是当前有关阻沙栅栏方面的研究热点。对于特定的阻沙措施,风场结构与风速无关 [10] ,本文以20 m∙s−1风速为例,给出了疏透度为0%和40%两种阻沙栅栏附近的风速等值线、速度矢量图(见图2)。
从图2(a)可以看出,受栅栏阻挡气流在栅栏前端被迫加速抬升(无孔隙的全部抬升,有孔隙的部分抬升),在栅栏顶部形成一个速度等值线密集的剪切层,且沿风向逐渐减弱并向地面附着,在迎风侧和背风侧分别形成一个小的和大的弱风区,这是有无孔隙的栅栏所形成二次流场的相同之处。两者区别在于紧密型栅栏后弱风区为一回漩涡区,而疏透度为40%的孔隙栅栏后则未出现涡区,且前者对风速衰减更多、低风速区域更大。上述风场计算结果与罗万银等 [11] 的风洞实验结果基本一致,他们的实验结果表明当疏透度大于30%后栅栏后端将不再出现回流现象。
3.2. 沙粒流经阻沙措施时的运动轨迹
图3为沙粒流经疏透度为0栅栏时的运动轨迹及速度变化图。由图3可看出,沙粒在遭遇紧密型栅栏阻挡时会发生两种行为,一种是在栅栏前端1 H范围内沉积,一种是从栅栏顶端越过栅栏。沙粒的行为与其粒径和风速密切相关,总体上有如下特点:越过栅栏的能力随粒径减小而增强,随风速增大而加强,如0.1 mm沙粒三种风速下均能越过栅栏,0.25 mm沙粒仅在20 m∙s−1风速下部分越过栅栏,而0.5 mm沙粒三种风速下均未能越过栅栏;越过栅栏后沙粒落回地面的附着点距栅栏的距离随风速增大而增大,低风速下小粒径沙粒会随回漩气流发生反向运动,随风速增大这种现象则消失。
图4为沙粒流经疏透度为40%栅栏时的运动轨迹及速度变化图。由图4可见,沙粒在遭遇栅栏时一部分会在前端沉积,一部分会穿过栅栏在后端沉积,沙粒穿过栅栏的能力与紧密型类似,不同之处在于沙粒是在栅栏孔隙而不是其顶端穿过,同时附着点与栅栏的距离也远远小于紧密型栅栏。
![](//html.hanspub.org/file/9-1770462x47_hanspub.png)
Figure 2. The wind field characteristics around the fence ((a) The distribution of wind speed, (b) The vectors of wind speed)
图2. 阻沙栅栏结构示意图两种疏透度栅栏附近流场特征((a) 速度分布,(b) 速度矢量分布)
![](//html.hanspub.org/file/9-1770462x48_hanspub.png)
Figure 3. Tracks and velocity of different diameter around the fence with η = 0 under different wind velocity (m∙s−1)
图3. 沙粒流经紧密型栅栏时的轨迹及其速度变化(m∙s−1)
![](//html.hanspub.org/file/9-1770462x49_hanspub.png)
Figure 4. Tracks and velocity of different diameter around the fence with η = 40% under different wind velocity (m∙s−1)
图4. 沙粒流经疏透度为40%栅栏时的轨迹及其速度变化(m∙s−1)
4. 讨论
综合图2和图3不难看出,紧密型栅栏对沙粒的拦截作用主要发生在前端,通过顶端越过栅栏的沙粒运动轨迹主要在气流的剪切加速区,沙粒在这里会被气流加速,风速越大时加速效果越明显,除低风速下小颗粒会随回漩气流反向运动沉积在栅栏的后端,其余情况下会随剪切气流回落到栅栏后部远端的风速恢复区,栅栏对这部分的沙粒拦截作用十分有限。同时也可发现,虽然紧密型栅栏后端对风速衰减效果很好,低风速区域也很大,但沙粒并为流经该区域,而是处于风速远高于来流风速的剪切层内。结合图2和图4可以发现,遭遇具有一定疏透度的孔隙栅栏时沙粒是从栅栏孔隙通过,虽然其间也会经历一个加速过程,但这个作用非常微弱,更重要的是穿越栅栏以后,沙粒进入了栅栏后端的低风速区,能量得到进一步衰减,因此很快发生大量沉积。这很好解释了对风速衰减更低低风速区域更大的紧密型栅栏阻沙效果却不如孔隙栅栏这一现象。由此可见,若使栅栏具备良好的阻沙能力必须同时满足两个条件:一是能够对来流风速进行大幅衰减,形成尽可能大的弱风区;二是促使来流沙粒进入弱风区。
5. 结论
本文采用欧拉–拉格朗日方法在栅栏流场中引入沙粒相,同时计算了流场和沙粒运动轨迹及速度,较单纯的以气流场论述栅栏的防护更为直观有效,深化了对栅栏阻沙机理的认识。同时也获得如下重要认识,阻沙措施若具有良好的阻沙性能必须同时满足两个条件,一是能够对来流风速进行大幅衰减,形成尽可能大的弱风区;二是促使来流沙粒进入弱风区,对阻沙措施的设计具有重要参考意义。
资助信息
本文得到国防环境安全专项资金资助(项目号:WZ20160402)。