基于风洞模拟的不同出风口高度导风板的聚风效应.pdf

第3 0卷第6期2 0 2 3年1 2月水土保持研究R e s e a r c ho fS o i l a n dW a t e rC o n s e r v a t i o nV o l.3 0,N o.6D e c.,2 0 2 3 收稿日期:2 0 2 2-1 0-2 9 修回日期:2 0 2 2-1 1-0 5 资助项目:内蒙古科技攻关项目“库布齐沙漠沿黄段流动沙丘综合治理技术研究”(2 0 2 1 G G 0 0 7 3);内蒙古农业大学“高层次人才引进科研启动项目”(N D Y B 2 0 1 8-3 6)第一作者:徐铮铮(1 9 9 6),女,辽宁凌源人,硕士研究生,研究方向为荒漠化防治。E-m a i l:1 6 3 7 2 7 8 6 4 4q q.c o m 通信作者:高永(1 9 6 2),男,内蒙古包头人,博士,教授,主要从事荒漠化防治研究。E-m a i l:1 3 9 4 8 8 1 5 7 0 91 6 3.c o mh t t p:s t b c y j.p a p e r o n c e.o r gD O I:1 0.1 3 8 6 9/j.c n k i.r s w c.2 0 2 3.0 6.0 2 2.徐铮铮,高永,韩彦隆,等.基于风洞模拟的不同出风口高度导风板的聚风效应J.水土保持研究,2 0 2 3,3 0(6):2 8 4-2 9 4,3 0 6.X uZ h e n g z h e n g,G a oY o n g,H a nY a n l o n g,e t a l.W i n dG a t h e r i n gE f f e c t o fW i n dG u i d eP l a t ew i t hD i f f e r e n tO u t l e tH e i g h t sB a s e do nW i n dT u n n e lT e s tJ.R e s e a r c ho fS o i l a n dW a t e rC o n s e r v a t i o n,2 0 2 3,3 0(6):2 8 4-2 9 4,3 0 6.基于风洞模拟的不同出风口高度导风板的聚风效应徐铮铮1,高 永1,2,3,韩彦隆1,2,3,张 超1,田晓宁1,王 鹏1(1.内蒙古农业大学 沙漠治理学院,呼和浩特0 1 0 0 1 1;2.内蒙古杭锦荒漠生态系统国家定位观测研究站,内蒙古 鄂尔多斯0 1 7 4 0 0;3.内蒙古自治区风沙物理与防沙治沙工程重点实验室,呼和浩特0 1 0 0 1 1)摘 要:目的 揭示不同出风口高度导风板下的风速流场特征及规律,阐明其聚风效能,进而为导风板配置模式提供理论支撑。方法 应用风洞模拟手段,测定不同风向夹角(=9 0,=6 0,=4 5 和=3 0)条件下,不同出风口高度(h=3c m,h=5c m,h=7c m,h=9c m)导风板的流场特征,分析1/2出风口高度处的相对风速。从而选取风向夹角()、出风口高度(h)、来流风速(V0)以及距离导风板的水平距离(x)等对导风板1/2出风口高度处各点位风速V(x)进行多元回归建模。结果(1)各风向夹角条件下,不同出风口高度导风板均形成一定范围的强风区和增速区,但不同出风口高度及风向夹角下其流场特征存在差异。=9 0 时,h=3c m和h=5c m导风板的增速区和强风区影响范围分别为最大;=6 0 时,h=9c m导风板形成的增速区和强风区影响范围均最大;=4 5 和=3 0 时,h=9c m导风板增速区和集流加速区影响范围均最大。(2)=9 0 和=6 0 时,各导风板1/2出风口高度处相对风速随水平距离的变化近似呈“N”形。由于狭管效应,在0.5H(H为模型高度)处达到相对风速峰值,而后气流发生扩散,到2H处降至最低值。随迎风夹角减小,导风板垂直风向投影面积急剧减小,过境气流在板后汇集,此时对水平测点1H后气流有较强加速效果。(3)1/2出风口高度处导风板风速预测模型为V(x)=-0.1 4 4+0.1 4 5h+(-0.0 2 2)+0.8 0 8V0+0.0 9 2x(R2=0.7 5 8)。结论 风向夹角=6 0 时,h=9c m导风板聚风效果最优;构建的预测模型参数选用较为合理,模型拟合程度较高,可为导风板干扰下1/2出风口高度处各点位风速预测提供理论支撑。关键词:荒漠化防治;风速流场;导风板;风洞模拟;聚风效应中图分类号:S 7 7 5 文献标识码:A 文章编号:1 0 0 5-3 4 0 9(2 0 2 3)0 6-0 2 8 4-1 1W i n dG a t h e r i n gE f f e c t o fW i n dG u i d eP l a t ew i t hD i f f e r e n tO u t l e tH e i g h t sB a s e do nW i n dT u n n e lT e s tX uZ h e n g z h e n g1,G a oY o n g1,2,3,H a nY a n l o n g1,2,3,Z h a n gC h a o1,T i a nX i a o n i n g1,W a n gP e n g1(1.C o l l e g eo fD e s e r tS c i e n c ea n dE n g i n e e r i n g,I n n e rM o n g o l i aA g r i c u l t u r a lU n i v e r s i t y,H u h h o t0 1 0 0 1 1,C h i n a;2.H a n g j i nN a t i o n a lP o s i t i o n i n gO b s e r v a t i o na n dR e s e a r c hS t a t i o no fD e s e r tE c o s y s t e m,O r d o s,I n n e rM o n g o l i a0 1 7 4 0 0,C h i n a;3.I n n e rM o n g o l i aK e yL a b o r a t o r yo fA e o l i a nP h y s i c sa n dD e s e r t i f i c a t i o nE n g i n e e r i n g,H u h h o t0 1 0 0 1 1,C h i n a)A b s t r a c t:O b j e c t i v eT h ea i m so f t h i ss t u d ya r et oe x p l o r et h ec h a r a c t e r i s t i c sa n dr u l e so f t h ew i n ds p e e df l o wf i e l du n d e rd i f f e r e n to u t l e th e i g h tw i n dg u i d ep l a t e s,c l a r i f y i t sw i n dg a t h e r i n ge f f e c t,a n dt h e nt op r o-v i d e t h e o r e t i c a l s u p p o r t f o r t h ec o n f i g u r a t i o nm o d eo fw i n dg u i d ep l a t e s.M e t h o dB ym e a n so fw i n dt u n n e ls i m u l a t i o n,t h e f l o wc h a r a c t e r i s t i c so f t h ew i n dg u i d ep l a t ew i t hd i f f e r e n ta i ro u t l e th e i g h t sw e r em e a s u r e du n d e rd i f f e r e n tw i n dd i r e c t i o na n g l e s,a n dt h er e l a t i v ew i n ds p e e da t 1/2a i ro u t l e th e i g h tw a sa n a l y z e d.T h ew i n ds p e e dV(x)a te a c hp o i n ta t1/2t h eh e i g h to ft h ea i rv e n to ft h ew i n dg u i d ep l a t ew a sm o d e l e db ym u l t i p l e r e g r e s s i o nb ys e l e c t i n g t h e a n g l eb e t w e e n t h ew i n dd i r e c t i o n(),t h eh e i g h t o f t h e a i rv e n t(h),t h ei n c o m i n gw i n ds p e e d(V0),a n dt h eh o r i z o n t a l d i s t a n c e f r o mt h ew i n dg u i d ep l a t e(x).R e s u l t s(1)U n d e rt h ec o n d i t i o no fw i n da n g l e,t h ew i n dg u i d ep l a t e s a l l f o r m e dac e r t a i nr a n g eo f s t r o n gw i n dz o n ea n dg r o w t hz o n e,b u ti t sf l o wf i e l dc h a r a c t e r i s t i c sw e r ed i f f e r e n tu n d e rd i f f e r e n ta i ro u t l e th e i g h ta n dw i n dd i r e c t i o na n g l e.Wh e nw a s 9 0,t h e i n f l u e n c e r a n g eo f t h eg r o w i n gz o n e a n d t h e s t r o n gw i n dz o n eo f t h eh=3c ma n dh=5c m w i n dg u i d ep l a t ew e r e t h e l a r g e s t.Wh e nw a s6 0,t h e i n f l u e n c e r a n g eo f t h eg r o w t hz o n ea n dt h es t r o n gw i n dz o n e f o r m e db yt h eh=9c m w i n dg u i d ep l a t ew a st h e l a r g e s t.Wh e nw a s4 5 a n dw a s3 0,t h e i n f l u e n c e r a n g eo fh=9c mw i n dg u i d ep l a t e i n c r e a s i n gz o n e a n dg a t h e r i n g f l o wa c c e l e r a t i n gz o n ew e r e t h el a r g e s t.(2)Wh e nw a s9 0 a n dw a s6 0,t h er e l a t i v ew i n ds p e e da t1/2o u t l e th e i g h to fe a c hw i n dg u i d ep l a t ec h a n g e dw i t hh o r i z o n t a ld i s t a n c ea p p r o x i m a t e l yi nt h es h a p eo fN.D u et ot h en a r r o wp i p ee f f e c t,m o s to f t h e mc o u l dr e a c ht h ep e a kv a l u eo f r e l a t i v ew i n ds p e e da t0.5H(Hs t a n d sf o rm o d e lh e i g h t),a n dt h e nt h ea i r f l o ws p r e a da n dd r o p p e d t o t h e l o w e s t v a l u eo f r e l a t i v ew i n ds p e e da t 2H.A s t h ea n g l eb e t w e e nw i n d w a r dd e c r e a s e d,t h ev e r t i c a lw i n dd i r e c t i o np r o j e c t i o na r e ao ft h ew i n dg u i d ep l a t ed e c r e a s e ds h a r p l y,a n d t h e t r a n s i t a i r c o l l e c t e db e h i n d t h ew i n dg u i d ep l a t e.A t t h i s t i m e,t h e a i r f l o wo f t h eh o r i z o n t a lm e a s u r i n gp o i n ta f t e r1Hh a da s t r o n ga c c e l e r a t i o ne f f e c t.(3)T h ew i n ds p e e dp r e d i c t i o nm o d e l o f t h ew i n dg u i d ep l a t ea t t h eh e i g h to f1/2a i ro u t l e tw a sV(x)=-0.1 4 4+0.1 4 5h+(-0.0 2 2)+0.8 0 8V0+0.0 9 2x(R2=0.7 5 8).C o n c l u s i o nI nt h ee f f e c to fw i n dg a t h e r i n ge f f e c t,w h e nt h ew i n dd i r e c t i o na n g l e=6 0,h=9c m w i n dg u i d ep l a t eh a s t h eb e s tw i n da c c u m u l a t i o ne f f e c t.T h ep a r a m e t e r so f t h ec o n s t r u c t e dp r e d i c t i o nm o d e lw e r er e a s o n a b l e,a n dt h em o d e lh a sah i g hd e g r e eo f f i t t i n g,w h i c hc a np r o v i d et h e o r e t i c a ls u p p o r t f o rt h ew i n ds p e e dp r e d i c t i o na t e a c hp o i n t a t t h eh e i g h to f 1/2o u t l e tu n d e r t h e i n t e r f e r e n c eo f t h ew i n dg u i d ep l a t e.K e y w o r d s:c o m b a t i n gd e s e r t i f i c a t i o n;w i n ds p e e d f l o wf i e l d;w i n dg u i d ep l a t e;w i n d t u n n e l t e s t;w i n dg a t h e r-i n ge f f e c t 在荒漠化防治进程中,输(导)沙工程措施可以改变风沙流的结构和方向改变沙物质的蚀积状态,包括导风板工程、羽毛排导沙工程和输沙断面工程等。其中导风板工程由栅栏工程发展而来,是一种典型的输沙工程措施,根据原理不同可分为下导风板和侧导风板1。常被应用于铁路2、公路的沙害3及风吹雪防治工作中4。目前导风板应用已深入到生产、生活等5方方面面,而在荒漠化防治领域,对于导风板组成的下导风工程的作用机理及工程体系等问题缺乏系统研究。其主要包括导风板不同倾角6、线路清沙7及其结构的优化配置等方面7-8。在风力强劲、沙源丰富的高大密集型流动沙丘地区,采用固沙、阻沙措施虽在短时间内可以达到控制沙害的目的,但随时间推移会造成更大的沙害9。导风板工程的作用机理是通过聚合加速作用使风沙流体加速通过区域,气流本身携带的沙物质基本不会降落在防护区内,防护区内积沙也会被吹走,达到聚风输沙的目的1 0。并且因其成本低廉、体积小、易于运输以及可实现大规模施工等优势,被广泛应用于输沙工程。因此,利用地形借助风力设置导风板工程措施在高大密集型流动沙丘的开发利用中更为有效。作为导风板工程的主要组成结构,其板面高度和开口高度直接影响工程聚风输沙能力1 1,合理确定导风板板面高度与出风口高度对于获得高聚风输沙效益具有重要意义。作为风沙流常用的3种研究方法之一,相较于野外观测,风洞试验可以有效控制风力条件,而相较于数值模拟过于理想化的模拟环境,风洞试验又几乎可以完美模拟自然状况下的流场分布,并还原沙粒受力运动状态。因此,基于前人研究,运用风洞模拟手段,对与风向呈不同角度设置的不同出风口高度前倾式4 5 导风板的流场特征进行系统研究,实现对其流场规律的认识与把握,而后采用多元回归法建立导风板聚风效应定量模型,为预测不同条件下导风板干扰下风速,进而确定其配置模式提供理论支撑。1 试验设计1.1 试验风洞试验在中国林业科学研究院沙漠林业试验中心可适应野外起伏地形的便携式环境风洞开展。由入口段、动力段、导流段、整流段、过渡段和试验段等组成的开口吹气直流式风洞,建成于2 0 1 3年,洞体总长3 0m(本文安装长度1 8m),横断面面积1.8 0m1.9 0m(内壁尺寸:宽高)。风洞通过变频器调节风速,设计最大风速1 8m/s,气流稳定性小于4%,风洞下边界层厚度2 8 c m。582第6期 徐铮铮等:基于风洞模拟的不同出风口高度导风板的聚风效应该风洞可以进行土壤风蚀、风沙动力地貌、风沙防治工程、种子风力传播等模拟试验研究。1.2 试验材料试验所用导风板原型为封闭式导风板,实际板面规格为2m(长)1m(宽),本试验模型按5 1将导风板缩小为4 0c m 2 0c m,材质为三合板,采用固定支撑架支撑于地面(图1)。模型在风洞试验段中的最大阻塞率为2.3%,满足风洞试验对阻塞率的要求1 2。图1 导风板模型示意图F i g.1 S c h e m a t i cd i a g r a mo fw i n dg u i d ep l a t em o d e l1.3 试验方案本试验选取4种出风口高度的前倾式4 5 导风板进行试验,分别为h=3c m,h=5c m,h=7c m和h=9c m(换算为野外实际规格分别为h=1 5c m,h=2 5c m,h=3 5c m和h=4 5c m)。设置4个风向夹角,=9 0,=6 0,=4 5 和=3 0。风速采用T S IA I R P R O风速计A P 5 0 0进行测定。在净风条件下,选择3种风速梯度(4m/s,6m/s,8m/s)进行空洞和不同模型气流速度场测定,设置数据采集频率为1s,待气流稳定后,记录吹刮时间为1m i n。测风高度分别为:1/2出风口高度处(1/2h)、出风口高度处(h)、导风板板面1/2处(1/2L),3 0c m和4 0c m。如图2所示,测点位置为模型迎风侧3H,模型背风侧:0.5H,1H,2H,3H,4H,5H,7H,9H,其中H表示模型高度(4个导风板模型高度分别为H=1 8.2c m,H=2 0.2c m,H=2 2.2c m和H=2 4.2c m)。测点换算为水平距离见表1。在空洞条件下测定风洞试验段上述各测点的风速作为对照风速。图2 风洞内不同角度导风板测点位置示意图(俯视图)F i g.2 S c h e m a t i cd i a g r a mo fm e a s u r i n gp o i n t so fw i n dg u i d ep l a t ea td i f f e r e n t a n g l e s i nt h ew i n dt u n n e l(t o pv i e w)表1 三种风速梯度下不同出风口高度导风板观测点位T a b l e1 O b s e r v a t i o np o i n t so fw i n dg u i d ep l a t ea td i f f e r e n t a i ro u t l e th e i g h t su n d e r t h r e ew i n ds p e e ds h a v i n g s导风板参数导风板与水平面夹角()/()出风口高度(h)/c m迎风侧测点位置/c m背风侧4 535 4.79.1,1 8.2,3 6.5,5 4.7,7 2.9,9 1.2,1 2 7.6,1 6 4.1 56 0.71 0.1,2 0.2,4 0.5,6 0.7,8 0.9,1 0 1.2,1 4 1.6,1 8 2.176 6.71 1.1,2 2.2,4 4.5,6 6.7,8 8.9,1 1 1.2,1 5 5.6,2 0 0.197 2.71 2.1,2 4.2,4 8.5,7 2.7,9 6.9,1 2 1.2,1 6 9.6,2 1 8.11.4 数据处理方法采用E x c e l 2 0 1 9对1m i n内所测持续稳定(4%)的风速数据进行平均,采用S u r f e r 2 1.1,O r i g i n2 0 2 1软件进行风速流场和相对风速等的制图,利用S P S S2 2.0的线性回归分析功能模块对导风板干扰下风速与各个因素之间的关系进行分析并构建预测模型。采用相对风速对导风板的聚风效能变化进行定量分析。相对风速计算方法如下:u=u(x,y)U(x0,y0)(1)式中:u 为相对风速;u(x,y)为距离导风板x处,高度y处的平均风速(m/s);U(x0,y0)为空洞下试验段水平距离x,高度y处的平均风速(m/s)。2 结果与分析2.1 导风板对风速流场的影响图47分别为=9 0,=6 0,=4 5 与=3 0 4种风向夹角条件下,导风板模型纵剖面风速流场特征。以空洞测得的风速流场作为对照。对比图3可知,在测定范围内,各风速流场大体可分为4个重新组合的风速能区,但能区大小、强弱各有不同。来流受到导风板阻挡后能量强烈衰减,导致其进风口区域风速明显下降,形成降速区并在模型中后部形成较为明显的风影区(后文将降速区及风影区统称为弱风速区)。当气流到达导风682 水 土 保 持 研 究 第3 0卷板处时,在板面某一点产生分离,一部分气流被迫抬升形成绕流,在导风板顶部形成增速区;另一部分气流沿板面向下流动与近地表气流相汇加速通过出风口,在出风口处产生狭管效应,形成一定范围的强风区。而后气流到达相对宽阔地域时发生扩散,形成风速减弱区。随着导风板与主风向夹角的减小,导风板垂直风向投影面积减小,流经导风板处的气流产生平行和垂直导风板方向的分量,使导风板侧导效应加强,流动的区域逐渐合并,划分的区域越来越少。此时,过境气流在背风侧汇集,又产生新的集流加速区。注:水平距离0表示模型摆放位置,“-”表示导风板模型迎风侧。图3 空洞风速流场F i g.3 A i rv e l o c i t y f l o wf i e l do f c a v i t y 当风向夹角=9 0 时,不同风速条件下不同出风口高度导风板流场特征如图4所示。4个出风口高度导风板顶部增速区分布范围依次为-1.4 4.1H,-2.4 1.9H,-1.10.9H,-2.31.6H。增速区内最大风速分别为4.2 8,4.3 5,4.2 1,4.2 8m/s,较C K同位置分别增加8.9 1%,1 0.6 9%,7.1 3%,8.9 1%。随风速增大,顶部增速区逐渐向导风板处移动。当h5c m时,随风速增大导风板顶部增速区影响范围呈减小趋势;h5c m时则呈逐渐增大趋势。h=3c m和h=7c m导风板在4m/s风速条件下无强风区。h=5c m和h=9 c m导风板强风区分布范围分别为-1.6 0.8H和-1.5 0.7H。强风区内最大风速分别为4.0 7,3.9 7m/s,较C K同位置分别增加2 9.2 1%,2 3.2 9%。当h5c m时,随风速的增大导风板强风区影响范围呈减小趋势;h5c m时呈逐渐增大趋势。随导风板出风口高度的增加,强风区的影响范围呈先增大后减小再增大趋势。弱风区的影响范围随风速的增大逐渐增大。且出风口高度越高,弱风区影响范围越小。扩散减速区随风速的增大影响范围逐渐增大。随出风口高度的升高,扩散减速区由原来的1个区逐渐分解为2个区域。当=6 0 时(图5),4个出风口高度导风板顶部增速区分布范围分别为-0.51.9H,-0.50.7H,-1.1 0.8H,-3 3.6H。增速区内最大风速分别为4.1 8,4.0 6,4.1 9,4.3 8m/s,较C K同位置分别增加6.4 1%,3.3 6%,6.6 7%,1 1.5 1%。当h7c m时,随风速增大导风板顶部增速区影响范围呈先增大后减小趋势;h7c m时,则呈减小趋势。h=3c m,h=5c m和h=7c m导风板强风区分布范围分别为-0.9 0.6H,0 0.4H,-1.4 0.8H,h=9c m导风板无强风区。强风区内最大风速分别为3.4 0,3.2 8,3.8 4m/s,较C K同位置分别增加8.2 8%,4.1 3%,2 0.0%。随导风板出风口高度增加,强风区影响范围呈增加趋势。随风速增大,弱风区无明显变化。扩散减速区影响范围随风速增大逐渐增大。随出风口高度增加,扩散减速区逐渐分解为两个区域,且影响范围逐渐减小。当=4 5 时(图6),4个出风口高度导风板顶部增速区分布范围分别为-0.40.7H,-0.30.5H,-1.6 1.8H和-3 9H。增速区内最大风速分别为4.0 7,3.9 6,4.1 9,4.3 8m/s,较C K同位置分别增加3.6 1%,0.8 1%,6.6 7%,1 1.5 1%。随风速增大,h=3c m导风板顶部增速区影响范围呈先增大后减小趋势;h=5c m和h=9c m导风板呈逐渐增大趋势;h=7c m导风板呈先减小后增大趋势。顶部增速区影响范围随出风口高度的增加呈先减小后增大趋势。在该风向夹角条件下,导风板垂直风向投影面积减小,导风板侧导作用加强。各风速条件下,4个出风口高度的导风板均无强风区。过境气流在背风侧汇集,形成新的集流加速区。h=3c m导风板无集流加速区,其余3个出风口高度导风板集流加速区分布范围分别为0.2 1.0H,-0.1 2.9H和1.3 9H。随风速及出风口高度的增加,集流加速区影响范围均逐渐增大。782第6期 徐铮铮等:基于风洞模拟的不同出风口高度导风板的聚风效应图4=9 0 导风板风速流场F i g.4 W i n dv e l o c i t yf l o wf i e l dd i a g r a mo f t h ew i n dg u i d ep l a t ew h e nt h ew i n dd i r e c t i o na n g l e=9 0 如图7所示,当=3 0 时,h=3c m和h=5c m导风板无顶部增速区,其余导风板顶部增速区分布范围分别为-2.54.5H和-33.4H。增速区内最大风速分别为4.2 0,4.2 3m/s,较C K同位置分别增加6.9 2%,7.6 9%。随风速的增大,h=3c m,h=5c m和h=9c m导风板顶部增速区影响范围逐渐增大,h=7c m导风板顶部增速区影响范围呈先减小后增加的趋势。顶部增速区影响范围随出风口高度的增加呈增大趋势。仅h=5c m导风板在6,8m/s风速条件下存在强风区,分布范围分别为-1.31H和-0.50.4H。强风区内最大风速分别为5.4 3,6.9 6m/s,较C K同位置分别增加1 2.5 5%,4 4.2 7%。h=7c m和h=9c m导风板集流加速区分布范围分别为0.42.4H,1.74.5H。随风速的增大,集流加速区影响范围逐渐增大。2.2 导风板相对风速变化特征相对风速反映导风板聚风能力的变化特征。相对风速越大,导风板聚风效能越高。对来流中颗粒物起疏导作用的主要为近地面的强风区,故对导风板1/2出风口高度处相对风速特征进行分析。如图8所示,不同的风向夹角条件下,随着导风板出风口高度的增加,导风板1/2出风口高度处相对风速变化趋势存在差异。当=9 0(图AC)时,导风板1/2出风口高度处相对风速随出风口高度的变化呈现出先升高再降低然后升高的趋势。4个出风口高度导风板相对风速平均值分别为1.1 4,1.3 1,1.1 3,1.2 5;=6 0(图DF)时,呈先升高后降低趋势。4个出风口高度导风板相对风速平均值分别为0.9 3,1.1 3,1.2 3,1.1 4;=4 5(图GI)时,呈逐渐升高趋势。4个出风口高度导风板相对风速平均值分别为0.8 5,0.9 0,0.9 9,1.0 3;=3 0(图JL)时,呈先升高后降低趋势。4个出风口高度导风板相对风速平均值分别为0.8 7,1.0 5,1.0 0,0.9 7。说明不同的风向夹角、不同出风口高度的导风板对相对风速变化882 水 土 保 持 研 究 第3 0卷的影响程度存在差异,=9 0,=6 0 和=4 5 时,h=9c m导风板对相对风速变化影响的波动性更大,=3 0 时,h=5c m导风板对相对风速变化影响的波动性更大。图5=6 0 导风板风速流场F i g.5 W i n dv e l o c i t yf l o wf i e l dd i a g r a mo f t h ew i n dg u i d ep l a t ew h e nt h ew i n dd i r e c t i o na n g l e=6 0 综上所述,不同风向夹角、不同出风口高度的导风板对相对风速变化的影响程度存在差异。当迎风夹角=9 0 和=6 0 时,各导风板在1/2出风口高度处相对风速均随水平距离的增加均呈先增加后减小再增加而后趋于稳定的趋势。由于出风口狭管效应,绝大部分导风板均能在水平测点0.5H处达到相对风速峰值。而后由于气流发生扩散,在水平测点2H(在4m/s风速条件下,h=9c m导风板为3H)处风速减弱为最低值。当迎风夹角=4 5 和=3 0 时,导风板垂直风向投影面积急剧减小,侧导作用加强。此时,过境气流在导风板后汇集,对水平距离1H后的气流有较大的加速效果。2.3 导风板聚风效果的预测模型多元线性回归模型是描述变量的相关性模型,在满足模型决定系数、相关性检验等基础上,可对因变量进行预测1 3。研究表明,导风板倾角、气流扩散角、山坡坡度、风向与下导风走向的夹角以及风速等均能作为下导风防护效果的标志1 4。在前人研究的基础上,本研究选取 风向的夹角()、出 风口高度(h)、来流风速(V0)和距离导风板水平距离(x)4个影响因素作为自变量,运用多元回归分析方法,建立导风板1/2出风口高度处风速多元回归预测模型。利用S P S S软件对3个风速条件下4个风向夹角、4个出风口高度的导风板各9个水平点位的4 3 2个样本数据进行多元线性回归。回归分析结果见表24。此多元回归模型方程调整后的决定系数R2=0.7 5 8(表2),说明自变量(风向的夹角、出风口高度、来流风速以及距离导风板的水平距离)可以解释因变量(导风板干扰下1/2出风口高度处的风速)7 5.8%982第6期 徐铮铮等:基于风洞模拟的不同出风口高度导风板的聚风效应的变化,模型拟合程度较高。对回归方程进行显著性检验(表3),得出p值为0.0 0 0导风板出风口高度距离导风板的水平距离风向夹角。可得最终预测模型为V(x)=-0.1 4 4+0.1 4 5h+(-0.0 2 2)+0.8 0 8V0+0.0 9 2x(2)式中:V(x)为导风板干扰下1/2出风口高度处的风速(m/s);h为出风口高度,本文取3,5,7,9c m;为风向夹角,取9 0,6 0,4 5,3 0;V0为来流风速,取4,6,8m/s;x为距离导风板的水平距离,取-3,0.5,1,2,3,4,5,7,9H(“-”号代表导风板模型迎风侧,H代表导风板模型高,单位为c m)。图6=4 5 导风板风速流场F i g.6 W i n dv e l o c i t yf l o wf i e l dd i a g r a mo f t h ew i n dg u i d ep l a t ew h e nt h ew i n dd i r e c t i o na n g l e=4 5 3 讨 论本研究对不同风速条件下,与风向呈不同夹角的不同出风口高度导风板的流场特性进行风洞模拟,基本掌握了各条件下的流场特征。在自然过程中,风向是不固定的,本研究通过调整导风板与风向的夹角,模拟了自然状况下多变的风向,更具实践意义。从风洞模拟结果来看,与韩彦隆1 5对不同倾角导风板聚风输沙机理的野外试验流场规律相一致,表明导风板风洞试验和野外试验对于流场的作用具有一定的相似性,可为实际应用提供理论参考依据。导风板布设与地面呈4 5 角,将导风板所在流层流体的动能和位能重新分配,加强了导风板上、下端流体运动的能量,使原来流经地表的“一次流”变为人为干扰下的“次生流”1 6-1 7。研究表明,4种不同出风口高度的导风板在不同风向夹角条件下均表现出明显的流速分区特征。大体可以划分为4个典型区域:顶部增速区、底部强风区、弱风区及扩散减速区。092 水 土 保 持 研 究 第3 0卷图7=3 0 导风板风速流场F i g.7 W i n dv e l o c i t yf l o wf i e l dd i a g r a mo f t h ew i n dg u i d ep l a t ew h e nt h ew i n dd i r e c t i o na n g l e=3 0 这些典型区域分别表示气流的不同运移特征,影响了导风板周围沙粒的跃移传输及沉积特征1 8。来流到达导风板处时在板面某一点产生分离,一部分气流被迫抬升,与上方气流汇合,在顶部形成增速区。另一部分气流沿板面向下流动与近地表气流相汇加速通过出风口,在出风口处产生狭管效应,导致气流携沙能力加强。辛林桂7、陈柏羽8等认为导风板与挡沙墙对风速流场的作用具有一定的相似性,都是一种带有尖缘的钝体,也有减速区,涡流区以及上方加速区和下方加速区生成,此与本研究结果一致。且加速区的范围和风速大小可以反映导风板输导效果。导风板与光伏电板结构类似,其对过境气流的作用机理也类似。陈曦1 9、唐国栋1 6、袁方等2 0、孙涛等2 1研究显示光伏电板出风口后方出现一定范围的风速增大,