高立式阻沙网风场扰动特征的风洞模拟
1.
2.
3.
Wind field disturbance induced by tall vertical sand barriers: a wind tunnel simulation
1.
2.
3.
通讯作者:
收稿日期: 2025-06-09 修回日期: 2025-07-07
| 基金资助: |
|
Received: 2025-06-09 Revised: 2025-07-07
关键词:
Keywords:
本文引用格式
马茜茜, 肖建华, 姚正毅, 薛娴.
Ma Xixi, Xiao Jianhua, Yao Zhengyi, Xue Xian.
0 引言
在开阔沙地上,风沙流运动较为稳定。高立式HDPE阻沙网的设置对近地层气流结构产生扰动,使其周围风速、风向及动能空间分布发生显著变化,进而在局部尺度上形成区别于原始地表风沙环境的“次生气流”结构。根据气流运动特征,对阻沙网两侧流场进行分区是解释其空气动力学的重要步骤。首次对阻沙网两侧次生流流场做出分区的是Plate[5]。而后基于风洞实验和数值模拟,国内外学者对Plate的分区进行验证和扩展。相关研究逐步构建了较为成熟的功能分区模式,适用于低疏透度高立式阻沙网的近地层气流结构特征研究[6-8](图1)。高立式阻沙网的布设导致气流运动发生深刻的变化,这种变化无疑会导致其附近沙粒沉积和运移规律的不同。Dong等[9]通过风洞实验确定了高立式阻沙网附近沙粒运动的3种模式:顺风运动、逆风运动和振荡运动。在这种模式下,产生5个沙粒运动区域——顺风区、振荡区、逆风区、振荡区、顺风区(图2)。
图1
图2
高立式HDPE阻沙网通过改变近地层风场结构,削弱地表风速,减少风蚀起沙过程,起到了控风、固沙的双重作用。然而,阻沙网在风沙流中的扰动机制较为复杂,其风场扰动特征受网高、网孔密度、排列方式以及来流风速等多因素影响[10]。阻沙网高度是影响阻沙网重要的设计指标。来流风速是影响风沙活动与阻沙效能的首要因素,而阻沙网高度直接决定其扰动范围与阻风能力。二者的组合效应决定了风场结构的空间调控效果,是高立式阻沙网性能发挥的关键参数。尽管高立式阻沙网在风沙防护工程中已有广泛应用,然而,在风沙灾害日益严重的沙漠公路沿线,尤其是在工程实际应用背景下,现有参数配置多基于经验或通用推荐,缺乏结合特定区域风沙特征的实验验证与科学依据。因此,本文以S315策达公路典型风沙区段为背景,选取实地采集的沙样,开展系列风洞模拟实验,系统分析不同来流风速(8、10、12、16、20 m·s-1)与阻沙网高度(5、8、12、15 cm)组合下的风场扰动特征,以期为该区域防沙工程的科学布局和参数优化提供理论支撑和数据基础。
1 材料与方法
1.1 实验装置
实验在中国科学院沙漠与沙漠化重点实验室风沙环境风洞中进行。直流开口吹气式风洞由动力段、整流段、实验段和扩散段4部分组成。洞体长37.78 m,实验段长16.23 m,实验段截面积为0.6 m×1.0 m(高×宽)。风速在1~40 m·s-1连续可调。边界层的厚度为120 mm。
1.2 实验设计
图3为室内风洞实验布设示意图。为了模拟戈壁地表,在实验段风洞地板上连续铺设4 m长的粗砂纸(阻沙网迎风测1 m,背风侧3 m),以增加地表粗糙度。在洞内垂直架设毕托管,使其平行于横截面的中轴线,分别测量空洞和阻沙网布设条件下的风速。垂直测点h共10个:0.3、0.6、1.2、2、4、8、12、16、20、24 cm;水平测点共17个:-16H、-8H、-4H、-2H、-1H、-0.5H、0、0.5H、1H、2H、3H、4H、6H、8H、12H、16H、20H(H为阻沙网高度,迎风测为-H,背风侧为+H)。阻沙网模型宽1.0 m,疏透度为42%,孔隙均匀布局且为方形。
图3
实验选取了4个高度阻沙网(5、8、12、15 cm)和5组风速(8、10、12、16、20 m·s-1)对阻沙网附近的风速廓线、流场进行测定,获得不同水平位置、不同高度的风速减弱率,以确定阻沙网阻滞作用,从而便于评估阻沙网的有效防护范围。
1.3 实验条件
风洞的雷诺数Re=105~106,气流处于完全湍流状态,可达到自模拟要求。风洞底板铺设了粗砂纸以模拟戈壁下垫面,边界层的厚度为20 cm左右,因此,可以假设风洞的流动深度几乎没有限制阻沙网模型上方边界层的发展。通过绘制空洞条件下不同自由流风速廓线图,可以发现遵从对数分布规律(图4):
式中:Z是距离地面的高度(cm);
图4
自由流速度为6、8、10、16、20 m·s-1的风速廓线分布符合对数规律,决定系数R2均超过0.9726(表1),对应的地表粗糙度分别为0.00017、0.00020、0.00025、0.00020、0.00014 cm。
表1 无阻沙网布设条件下高度 Z 与风速 UZ 的线性拟合
Table 1
| U/(m·s-1) | Z与UZ 的线性拟合 | R2 | Z0/cm |
|---|---|---|---|
| 8 | lnZ = 1.2696UZ - 8.6856 | 0.9726 | 0.00017 |
| 10 | lnZ = 1.0221UZ - 8.5251 | 0.9852 | 0.00020 |
| 12 | lnZ = 0.8365UZ - 8.2878 | 0.9861 | 0.00025 |
| 16 | lnZ = 0.6468UZ - 8.5006 | 0.9854 | 0.00020 |
| 20 | lnZ = 0.5359UZ - 8.8871 | 0.9864 | 0.00014 |
1.4 防风效能
阻沙网干扰下的风速变化可以用归一化无量纲相对风速
式中:
2 结果与分析
2.1 流场特征
近地表流场是反映阻沙网功能的重要指标。由图5可以看出,当阻沙网布施于平坦沙地时,其前后一定范围内的流场受到强烈干扰,气流结构不再保持平衡状态,气流运移规律将会大大改变。从阻沙网前后的流场分布来看,气流在经过时存在明显的功能分区,主要分为迎风侧减速区、网顶加速区、背风侧减速沉降区3个功能区。气流在接近阻沙网的过程中,遇阻减弱,在迎风侧紧邻阻沙网并低于阻沙网高度的一定区域出现能量的损失,形成一个小的气流减速区,范围在-4H至阻沙网布设位置。在此区域内,因底层风速降幅较大,沙粒以沉降方式堆积,久而久之,迎风侧形成沙丘[14-15]。由于阻沙网阻挡,流线被抬升,疏透度越低,抬升作用越强,实验中所采用的阻沙网疏透度为42%,介于最优疏透度30%~60%,故而流线抬升作用较明显。但研究已经指出,当疏透度大于60%时,流线提升幅度将大大减小[14,16]。阻沙网对气流的抬升作用使近地层风动量向垂直方向传递,导致其背风侧顶部区域气流速度增强,形成上方加速区。与此同时,迎风侧气流在通过阻沙网时受到阻挡,部分动能损耗并转化为湍流能量,从而在背风侧近地面区域形成明显的减速与沉降区。阻沙网迎风侧和背风侧的减速区是沙粒最主要的沉积区[15,17-18],但分析各组流场发现,随着阻沙网高度的增加,两个减速区面积都在增加,背风侧减速区面积始终大于迎风侧,其风速的减弱速率和幅度也明显高于迎风侧,例如当自由流风速为8 m·s-1时(图5A),背风侧减速区内最低风速低于2 m·s-1,而迎风侧减速区内风速为4~6 m·s-1。气流经过阻沙网背风侧的减速沉降区后,随着与阻沙网水平距离的增大,气流逐渐摆脱阻沙网的影响,风速逐渐恢复到自由流。阻沙网越高,迎风侧气流抬升强度越大,导致网顶的加速区有向上风向移动的趋势。
图5
图5
不同来流风速下不同高度阻沙网流场分布
Fig.5
Flow fields around fences of different heights under under different inflow wind speeds
不同高度的阻沙网背风侧最小流速出现的起始位置不同,以风速10 m·s-1为例(图5B),4个高度阻沙网背风侧最小风速均为4 m·s-1,但最小位置有很大不同。当网高为5 cm时,最小风速在1.85H,即9.25 cm;当网高为8 cm时,最小风速在1.75H,即14 cm;当网高分别为12、15 cm时,最小风速分别在4H、10H,分别相当于背风侧距阻沙网布设位置48 cm和150 cm。因此,阻沙网越高,速度降低越慢,也可以认为由最小风速包围的区域范围越大,减速沉降区的范围越大。相反,当阻沙网较低时,风速在背风侧很短的距离内达到最小值,继而增大,很快恢复到自由流风速,减速沉降区面积较小,不利于沙粒的沉降。
2.2 风速廓线
风速廓线是气流结构的重要特性。受地面摩擦阻力的影响,在近地表层稳定床面上的风速廓线随高度呈对数分布[19]。但在阻沙网布设条件下,风速廓线并不完全服从对数分布规律。当来流风速分别为8 m·s-1和16 m·s-1时,在-4H处,4种高度阻沙网风速廓线的形式仍与稳定床面基本一致,基本服从对数分布规律;在-0.5H处,阻沙网高度以下风速逐渐减弱,阻沙网越高,风速减弱越强。在下风向0.5H~6H或0.5H~8H在阻沙网高度以下的风速减弱趋势达到最大,这归结于阻沙网前的气流拥塞效应[20]。自阻沙网背风侧12H处继续向下风向延伸,风速衰减程度逐渐减弱。在距网体背风侧20H处,风速剖面形态逐渐接近入流剖面,说明流场扰动效应趋于减弱,气流逐步恢复至原始流动状态。然而,此处风速仍低于入流风速,表明阻沙网的残余影响依然存在(图6)。
图6
图6
8 m·s-1和16 m·s-1来流风速下不同高度阻沙网两侧风速廓线
Fig.6
Wind speed profiles around fences of different heights under wind velocity of 8 m·s-1 and 16 m·s-1
2.3 风速水平分布
图7是在3种来流风速(8、12、16 m·s-1)条件下,在距离风洞地板一定高度h(0.3、0.6、1.2、2、4、8、12、16、20、24 cm)水平方向的风速变化,可以看出,在3个自由流风速下,阻沙网附近的风速分布大致相似,表明自由流风速的影响很小。共同特征是,风速减弱程度在低于阻沙网高度时明显。当阻沙网高度为5 cm时,h=4 cm的风速已基本不受其影响。当阻沙网高度为8 cm时,h=8 cm的风速也没有明显的减弱趋势。当阻沙网高度为12 cm时,h=12 cm的风速也基本恢复到进口风速大小。在距风洞地表h=12 cm的高度,15 cm高度阻沙网背风侧气流扰动还很强;但在h=16 cm高度,风速受阻沙网扰动程度下降,在h=20 cm高度,风速已基本不受阻沙网干扰。根据阻沙网对气流的干扰在超过网高时较弱[20-22],这里主要关注h=12 cm高度以下的风速,据此分析风速减弱率,即防风效能。
图7
图7
不同来流风速下阻沙网水平方向风速分布特征
Fig.7
Distribution of flow velocity in horizontal direction of fences different inflow wind speeds
2.4 防风效能
在低于阻沙网高度范围内,有最强的风速减弱效应。因此,根据图6,对于4个高度阻沙网,分别选择等于及低于2.4、4、8、12 cm的风速,采用
图8
图8
风速为12 m·s-1时,不同高度阻沙网的平均防风效能
Fig.8
The average reduction coefficient of horizontal wind velocity of different height fences, when the wind speed is 12 m·s-1
3 结论
从阻沙网前后的流场分布来看,气流在经过时存在明显的功能分区,主要分为迎风侧减速区、网顶加速区、背风侧减速沉降区3个功能区。阻沙网迎风侧和背风侧的减速区是沙粒最主要的沉积区。背风侧减速区面积始终大于迎风侧,风速的减弱速率和幅度也明显高于迎风侧。
迎风侧减速区从-4H开始,自阻沙网布设位置直至6H或8H,防风效能达到最大,此范围内沙粒大量沉积。由于16H位置风速达到来流风速56%~67%,已基本接近70%,故认为阻沙网在下风向的防护距离至少可延伸至16H。
参考文献
Evaluation of the effectiveness of an expressway sand protection system in a gobi region:case study of the Ceke-Ejina expressway,Ejina Banner,China
[J].
华南沿海风沙危害防护体系及其效益分析
[J].
A wind tunnel study of turbulent flow around single and multiple windbreaks,part I:velocity fields
[J].
两种形状镁水泥板沙障防风积沙效益的风洞试验
[J].
Threshold velocity for wind erosion:the effects of porous fences
[J].
Aerodynamics and morphodynamics of sand fences:a review
[J].
Relationship between shelterbelt structure and mean wind reduction
[J].
Evaluating the optimal porosity of fences for reducing wind erosion
[J].
Optimal windbreak design for wind-erosion control
[J].
A wind tunnel simulation of the mean velocity fields behind upright porous fences
[J].
The blocking effect of the sand fences quantified using wind tunnel simulations
[J].
A wind tunnel simulation of the turbulence fields behind upright porous wind fences
[J].
Wind tunnel test on the effect of metal net fences on sand flux in a Gobi Desert,China
[J].
Wind tunnel observations about the shelter effect of porous fences on the sand particle movements
[J].
Disturbance of the inclined inserting-type sand fence to wind-sand flow fields and its sand control characteristics
[J].
Influence of Golmud-Lhasa section of Qinghai-Tibet Railway on blown sand transport
[J].
Characteristics of the disastrous wind-sand environment along railways in the gobi area of Xinjiang,China
[J].
Effects of windbreak structure on wind flow
[J].
/
| 〈 |
|
〉 |

甘公网安备 62010202000688号