EFFECTS OF WIND SPEED AND DURATION ON PROCESS OF SNOW DISTRIBUTION AROUND A CUBE
-
摘要: 建筑物周边风致积雪是冬季多雪大风地区最常见的自然灾害之一。由于风吹雪的机理复杂，虽然数值模拟方法近几年取得了丰硕的成果，但是仍缺乏在稳定风雪条件下的广泛验证，限制了其预测方法的进一步提升和推广。为了揭示风速和持续时间对立方体周边积雪演化过程的影响，该文在西北严寒地区的稳定吹雪条件下，对不同风速下立方体模型周边不同时刻的积雪深度进行实测。讨论了不同工况下堆积区域最大积雪深度、积雪剖面的变化以及立方体迎风区两侧的积雪堆积形态。结果表明：在不同工况下，立方体周边风致积雪均在迎风和背风壁面处堆积，在模型前缘和后方出现侵蚀，随着风速的变化，积雪深度、积雪范围、积雪堆积位置都会发生相应的改变；随着时间的变化，积雪最大深度位置也会转移。Abstract: Snow accumulation around buildings is one of the most common natural disasters in snowy and windy areas in winter. Although numerical simulation has achieved fruitful results in recent years, it still lacks an extensive verification due to its complex mechanism, which limits the effectiveness and accuracy of the prediction methods. To reveal the influences of wind speed and duration on the process of snow accumulation around the cube, the field measurements of a stable blowing snow condition in severely cold northwestern region were carried out by measuring the depth of snow at different moments around the cube model at different wind speeds. The maximum snow depth and snow profile change in the accumulation area under different conditions and the snow accumulation patterns on both sides of the windward region were discussed. The results show that under different conditions, the snow is accumulated against the wall of the cube, and erosion occurs at the front and rear of the model. With the change of the wind speed, snow depth, snow deposition range, snow accumulation location will change accordingly. As time changes, the position of maximum snow depth will also change.
-
Keywords:
- snowdrift /
- cube /
- field measurement /
- wind speed /
- duration
-
在多雪大风的寒冷地区，风致建筑物周边积雪堆积会引发室外设施和通行问题，另外，长期积雪的覆盖和融化也会加速附属设施的腐蚀。建筑物周边的积雪与其复杂的三维风致绕流状态密切相关，如壁面处气流的分离、建筑物的迎风侧和背风侧形成的涡流等都会对建筑物附近的积雪堆积情况产生巨大影响。为了探究建筑物附近的积雪现象，使用雪颗粒模拟物降雪的风洞试验(Kind[1]、Smedley等[2]、Nemoto等[3]、Liu等[4]、Nikolas等[5]、Ding等[6]、Yan等[7])及数值模拟方法(Liston等[8]、Tominaga[9]、周晅毅等[10-13]、Jiang等[14]、孙晓颖等[15]、徐枫等[16]、刘多特等[17-18])进行了很多研究，这些研究成果往往需要以现场实测(Oikawa等[19]、Tsuchiya等[20]、Thiis等[21])的结果作为其验证依据。
影响建筑周边积雪堆积的因素很多，如雪颗粒的粒径、雪颗粒之间的粘聚力、雪密度、种类、地面粗糙度、风速、风吹雪时长，还有温度、湿度等其他环境因素[22-25]。由于风吹雪引起的积雪堆积是一个变化的过程，因此风速和持续时间对积雪的堆积形态会产生明显的影响。
Oikawa等[26]在日本札幌对立方体周边的积雪分布进行了现场实测，风速分别为1.7 m/s、4.6 m/s和5.2 m/s时，立方体周围的积雪分布有明显差异。Liu等[4]的风雪联合试验中也发现了类似的现象，其试验结果表明迎风区积雪堆积有两种模式，当风速为1.5 m/s和2.5 m/s时，积雪在迎风侧为一个整体；当风速为3.5 m/s和4.5 m/s时，积雪分为两部分，一部分在迎风侧底部，另一部分在迎风侧前方。周晅毅等[13]采用拉格朗日方法模拟立方体周围风致积雪漂移现象时发现，风速越高，计算域中雪颗粒总数的减少幅度就越大，计算域内的侵蚀或沉积质量通量也越大。刘多特[27]对比5 m/s、15 m/s、25 m/s风速下桥面雪浓度和积雪重分布情况，结果表明随着来流风速的增大，桥面过剪切范围随之增大，积雪面积将有所减小，而当风速超过某一限值时，特征吹雪浓度可能发生轻微减小，此时，局部积雪深度将不再持续增大。赵雷等[28]通过研究风雪流对低矮建筑屋盖的作用时发现，随着风速变化，积雪漂移使积雪形态发生了改变。以上的研究均表明，风速是影响积雪分布的一个重要参数，但是其作用机理和影响特征还不明确。
Beyers等[29]和Tominaga等[30]研究表明积雪的累积过程是随时间变化的，这意味着积雪深度在不同位置以不同的模式随时间不断变化。在李雪峰[31]的风洞试验中，480 s、960 s、1440 s时立方体周边的积雪分布有明显的差别。赵有栓等[32]采用RSM湍流模型和k-kl-ω湍流模型分别模拟了立方体模型周围的风致积雪运动，对比12 h、24 h、36 h三个时间节点的数值模拟结果可以看出，在风致积雪运动中，立方体周围的积雪深度随时间不断变化。
Michael等[33]选定粉碎的胡桃壳模拟雪颗粒，在纽约伦斯勒理工学院水洞实验室进行的高低屋盖屋面积雪重分布试验，试验得到了风速和持续时间对高屋盖屋面积雪传输率的影响。张玉杰等[34]基于Euler-Euler方法，采用Mixture多相流模型，对轻型四坡房屋表面的积雪分布情况进行模拟，研究表明，随着风速和模拟时间的增大，屋面积雪整体减少。
从已有的研究来看，积雪堆积的演化过程受到风速、持续时间等各种环境因素的影响，且其实测建筑周边积雪堆积情况与风洞试验和数值模拟结果吻合较差。为了确保多雪地区冬季的生存环境安全、舒适，精确预测不同地区建筑周边积雪分布的演化过程是解决该问题的主要策略，而这又受到风吹雪现象复杂机理的限制。目前尚没有成熟的预测方法，更缺少以自然雪为对象，且条件较为稳定的实测数据对各种预测方法进行验证，所以建筑物附近积雪现象的再现性还没有被很好地阐明。
本文选择在冬季新疆北疆的玛依塔斯风区，对不同风速下立方体模型的不同积雪堆积阶段进行现场实测，测量期间无降雪且环境因素稳定。分别对不同风速同一时刻和不同时刻同一风速的立方体周边积雪分布进行对比研究，给出了风速和持续时间对立方体周边积雪演化过程的影响，为建筑周边积雪堆积预测方法提供了基础的对比和验证数据。
1 现场模型试验概况
1.1 实测场地
实测时间为2019年12月−2020年1月，地点为新疆北疆九大风区之一的玛依塔斯风区。该区域地处吾尔喀什尔山和加依尔山之间，东西两侧分别为准格尔盆地与塔额盆地。高山、盆地高差大，走廊型山间盆地贯穿于山地之间，形成了狭长的大气环流通道，导致玛依塔斯风力强劲且风速相对稳定，测试区域如图1所示。
此处冬季盛行东风，在降雪天气时可能会变为西风，每次刮风时的风向较为稳定[35]，本次实测采用风杯式风速风向仪测量了实测场地处测试时段内的风向及风速，风向为WSW，风速时程如图2所示。
玛依塔斯气温较低，冬季平均气温约为−12 ℃，最高气温不超过0 ℃，平均湿度约为71%，冬季降水量为109.4 mm，冬春两季有大范围降雪，提供了较为稳定且十分充足的新雪颗粒。
1.2 风吹雪特性
雪颗粒落在放有标尺的黑色硬纸板上，拍照记录多组样本，并采取图像处理获得雪颗粒的投影外形，根据等效面积法，得到所采样本雪颗粒平均粒径约为295 μm，如图3所示。
通过筛雪收集一桶雪样本测量质量和体积的方式测量雪的密度，如表1所示。
表 1 雪的堆积密度Table 1. Bulk density of snow雪种类 密度/(g/cm3) 新雪 0.03~0.12 >24 h 0.18 >48 h 0.29 >72 h 0.30 沉积雪 0.29~0.42 通过筛雪的方式使雪颗粒在圆桶上形成圆锥且再次降落在圆锥上的雪颗粒可以自然下落时，测量圆锥侧边与底面形成的夹角为其休止角，如表2所示。
表 2 雪的休止角Table 2. The angle of repose雪种类 休止角/(°) 新雪 60.78 >24 h 58.14 >30 h 56.97 表层积雪 40.54 底层积雪 32.03 表3为利用自制的雪通量仪测量不同高度雪质量和记录测量时间的方式得到的不同风速下不同高度的雪通量，仪器如图4所示，左侧为进风口，右侧为出风口，中间采用设置隔板的方式存雪。
表 3 雪通量Table 3. Mass flux高度/cm 6.7 m/s雪通量/(g·m−2·s−1) 11.4 m/s雪通量/(g·m−2·s−1) 2.5 68.75 127.78 7.5 10.42 23.61 12.5 4.17 11.81 17.5 2.08 6.94 22.5 2.08 6.25 27.5 0.00 2.78 结果表明：绝大部分的雪颗粒在距离地面约2.5 cm高区域内运动，由此可见新雪颗粒处于跃移层近地面位置，这与现场观察到的雪颗粒运动方式一致，当风速增大时，雪通量随之增大，这也是导致积雪量增大的原因。
1.3 试验概况
为了便于在恶劣的环境中根据风向变化调整模型，立方体由边长为B=0.3 m、厚5 mm的木板制成。将其放在长6.6 m、宽3.51 m、厚5 mm的木底板上，最外圈测点位置距离模型中心前、后、左、右的距离分别为5.5B、15.5B、5.5B、5.5B，使得风雪流在模型周围造成的积雪都落在平坦的底板上，从而观测更为准确。风速采用风杯式风速仪测量，风速测试点位于模型前侧来流方向1 m高度处。
为了描述方便，将整体测试区域划分为迎风区、尾侧区和背风区，在模型周边布置490个测点，每个测点插入高18 cm，直径2.3 mm的钢签，钢签上每间隔5 mm做不同颜色的标识作为雪深量度，测点布置及区域划分如图5(a)所示。在模型底板四周以不同时间间隔拍摄积雪分布的照片，每组照片将全部测点包括在内，读取钢签标尺刻度作为积雪深度，实测照片如图5(b)所示。
在风雪流中，采用风速风向仪测量风向及风速，并将模型调整为正面迎风，在模型前侧放置挡风板，清扫干净底板上安置模型时留下的积雪，避免残留雪对试验结果造成影响，随后撤去挡风板，通过不同的角度拍照记录立方体周边的积雪分布，由于立方体周围积雪会随着时间变化，因此要尽可能在积雪发生显著变化且最短时间内完成这一过程的记录。风速选取1 m高度处平均风速v =6.7 m/s、v =11.4 m/s。10 min内风速风向稳定，风速时程如表4所示。
表 4 风速测量数据Table 4. Wind speed measurement data工况 1 min 2 min 4 min 8 min 10 min 工况1：平均风速v/(m/s) 6.2 6.8 7.1 6.8 6.8 工况2：平均风速 v/(m/s) 11.0 10.8 12.1 12.1 10.8 时间选取t=1 min、2 min、4 min、8 min、10 min作为记录时刻，且每组拍摄时间在30 s内完成，确保每组照片中第一张和最后一张的拍摄时间相差最小。选取积雪变化较大的t=2 min、4 min、8 min作为三个阶段进行分析，积雪深度采用立方体边长B无量纲化进行整理。
2 研究结果
2.1 最大积雪深度
如表5所示，在迎风区，6.7 m/s风速下积雪最大深度位置保持在X/B=−1.3，Y/B=0.17的位置，11.4 m/s风速下积雪最大深度位置相对6.7 m/s风速下稍微靠近模型了一些，在X/B=−1.1，Y/B=0的位置；在尾侧区，6.7 m/s风速下积雪最大深度位置保持在X/B=1.83，Y/B=1.83的位置，11.4 m/s风速下积雪最大深度位置相对6.7 m/s风速下稍微远离模型了一些，在X/B=2.17，Y/B=1.5的位置；在背风区，积雪最大深度位置一直在壁面X/B=0.5，Y/B=0处。在迎风区和尾侧区，积雪最大深度随时间增长在逐步增大，风速越大，积雪最大深度越大。在背风区积雪最大深度一直保持稳定，随时间变化不大。
表 5 积雪最大深度及其位置Table 5. Maximum snow depth and location堆积区域 平均风速v/(m/s) 时间t/min X/B Y/B D/B 迎风区 6.7 2 −1.30 0.17 0.05 4 −1.30 0.17 0.11 8 −1.30 0.00 0.16 11.4 2 −1.10 0.00 0.08 4 −1.10 0.00 0.13 8 −1.10 0.17 0.18 尾侧区 6.7 2 1.83 1.83 0.05 4 1.83 1.83 0.10 8 1.83 −1.17 0.14 11.4 2 2.17 1.50 0.07 4 1.83 1.50 0.10 8 2.17 1.50 0.16 背风区 6.7 2 0.50 0.00 0.08 4 0.50 0.00 0.09 8 0.50 0.00 0.09 11.4 2 0.50 0.17 0.08 4 0.50 0.17 0.08 8 0.50 0.17 0.07 注：X/B为顺风向位置； Y/B为横风向位置； D/B为积雪深度。 如图6所示，在2 min时，6.7 m/s风速下积雪最大深度位置出现在背风区，而11.4 m/s风速下积雪最大深度位置在背风区和迎风区同时存在；在4 min时，6.7 m/s风速下积雪最大深度位置转移到迎风区，11.4 m/s风速下积雪最大深度位置只出现在迎风区；在8 min时，两种风速下积雪最大深度位置均维持在迎风区不再变化。可以看出，立方体周边风致积雪最大深度位置不受风速影响，但随时间变化。
2.2 积雪堆积形状
图7中图7(a)~图7(c)分别给出了t=2 min、t=4 min、t=8 min在6.7 m/s风速下和11.4 m/s风速下的顺风向中轴线Y/B=0处的积雪堆积情况，在图7(c)中，还给出了Oikawa等[26]的部分实测数据作为参考。积雪深度根据其实测模型边长进行无量纲化处理，需要说明的是，本文的测试过程中并未出现降雪，Oikawa等[26]的测试是在降雪条件下进行的，因此，存在70 mm的基础雪深，不对积雪深度做比较。另外，Oikawa等[26]实测风速为1 m高度处5.2 m/s，与本文两种风速下的工况形成不同大小风速的积雪形态对比。
如图7所示，在靠近立方体迎风侧壁面处，由于驻涡的产生，使得风速在立方体近壁面处减小，无法带动雪颗粒运动从而导致雪颗粒堆积，驻涡造成的回流区域风速和来流风速方向相反且风速较大，使其发生侵蚀，在靠近立方体背风侧壁面处，气流绕过立方体，在两侧产生马蹄涡，在模型后侧有回流形成，在马蹄涡和回流涡的交界处，风速减小，积雪从壁面中心处到两个角部后方产生堆积，由于回流涡带动的反向气流使风速加强，在背风区也有侵蚀产生，但效果不是太明显。
在t=2 min时，6.7 m/s风速下和11.4 m/s风速下迎风区积雪最大深度位置相同，均在X/B=1.2处，11.4 m/s风速下积雪最大深度为6.7 m/s风速下积雪最大深度的2倍，背风区无太大差异。在t=4 min时，6.7 m/s风速下，迎风区积雪在X/B=−1.3处积雪达到最大深度，在X/B=0.8左右处侵蚀，背风区积雪在壁面处积雪达到最大深度，在X/B=1.3左右处侵蚀。11.4 m/s风速下，迎风区积雪深度增大，积雪最大深度位置由于雪颗粒惯性力较大，向靠近模型侧偏移，达到X/B=−1.1处，背风区积雪深度增大，积雪范围增长，当壁面处风速增大，绕过模型产生的马蹄涡和流过模型产生的回流涡交界形成的积雪堆积区域减小，积雪深度也减小。在t=8 min时，对比Oikawa等[26]的实测数据，当风速减小为5.2 m/s时，迎风区积雪最大深度位置由于雪颗粒惯性力较小，向远离模型侧偏移，达到X/B=−1.5处，背风侧绕过模型产生的马蹄涡和流过模型产生的回流涡交界形成的积雪堆积区域增大，壁面处积雪深度加大。风速越大，迎风区积雪最大深度位置会更接近模型，背风侧靠近壁面处积雪堆积区域减小，积雪深度减小，背风侧后侧积雪深度增大，范围增长。背风侧远离立方体的位置，受到风致积雪作用较小，积雪堆积逐渐恢复到基础雪深，Oikawa等[26]实测数据基础雪深为70 mm，本实测数据基础雪深为0 mm，所以存在较大差异。
为了对比不同风速在不同时间的尾迹变化情况，选取Y/B=−2.17处的积雪分布情况进行了单尾侧区分析，由于测试条件限制，此处虽偏离尾迹中轴线，但也能反映出尾迹变化趋势。
对比图8中6.7 m/s风速下和11.4 m/s风速下尾侧区积雪堆积情况发现，尾侧区积雪呈现出流线型，风速越大，积雪深度越大，积雪范围也越广。相对于6.7 m/s风速下的积雪最大深度位置，11.4 m/s风速下的积雪最大深度位置离模型位置更远。对比图8(a)~图8(c)三个时间节点发现，在实测的8 min内，随着时间的增长，两种风速下各自的积雪深度在不断增加，但是积雪范围并没有改变。11.4 m/s风速下的积雪范围要比6.7 m/s风速下积雪范围更大，由此推测，积雪范围越广，积雪深度最大位置相对于模型距离越远。为了观测积雪范围是否发生变化以及验证该结论，选取11.4 m/s风速的工况进行了长达32 min的堆积测量，如图9所示，尾侧区积雪范围明显增大，积雪深度最大位置也随积雪范围的增大而后移。同时发现尾侧区积雪至少在8 min内积雪范围为一固定值，不随时间改变而改变，超过8 min后，尾侧区积雪范围随时间增大而增大，积雪深度也相应的增加。
![]()
图 9 11.4 m/s风速下尾侧区积雪对比Figure 9. Comparison of snow distribution in the lateral region at 11.4 m/s2.3 迎风区两侧的积雪情况
图10中图10(a)~图10(c)为X/B=−0.5处t=2 min、4 min、8 min时不同风速下的积雪分布情况，在模型横风向的两侧，6.7 m/s风速下，可以看到雪堆积，且随着时间的增大，积雪深度也在逐渐加深，而11.4 m/s风速下，8 min前该区域无积雪堆积，在8 min时，有少量堆积。可以见得，在风速不同的情况下，立方体迎风区两侧的积雪分布有明显的差异，当风速较小时，在迎风区两侧很快就有积雪产生，且随时间增长积雪深度逐渐增大，当风速较大时，该区域短时间内很难形成积雪堆积，经过长时间的累积才会出现积雪堆积现象。迎风区两侧是否存在积雪堆积现象，会引起周围流场变动，从而对立方体周边积雪分布造成影响，尤其是迎风区壁面以及背风区的积雪情况。
在迎风区壁面处，由于驻涡和马蹄涡的影响，积雪从壁面中心到两个角部堆积，形成三角形堆积形态，11.4 m/s风速下的积雪深度要比6.7 m/s风速下的积雪深度大。6.7 m/s风速下，该区域积雪较为稳定，随时间变化不大，11.4 m/s风速下，积雪深度随着时间的增长一直在增大。当风速较小时，迎风区两侧已经存在积雪，对迎风区壁面处风场有稳固的效果，所以此处积雪堆积情况较为稳定，当风速较大时，迎风区两侧无积雪堆积，迎风区壁面处积雪持续增加。由图9(b)中11.4 m/s风速下的积雪分布情况来看，迎风区两侧已经产生积雪堆积，可见随着时间的发展，11.4 m/s风速下迎风区壁面的积雪也会趋于稳定。
图11中图11(a)~图11(c)为背风区X/B=1.83处t=2 min、4 min、8 min时不同风速下的积雪分布情况，在立方体模型两侧Y/B<−1，Y/B>1的位置，积雪随着时间的增大深度不断加深，且具有一定的不对称性，根据在野外的观测，这种差异是因为来流风的不稳定性以及风吹雪的不均匀性引起的，在背风区模型后方的位置，积雪深度随时间的增长不断加深，风速越大，积雪深度越大。11.4 m/s风速下在−1<Y/B<1的范围内都有积雪堆积，而6.7 m/s风速下仅在−0.5<Y/B<0.5范围内存在积雪堆积现象，11.4 m/s风速下的积雪范围要比6.7 m/s风速下的积雪范围更大，且与两侧尾迹连在一起，这是因为风速大的时候，迎风区两侧区域没有积雪堆积，从此区域吹来的风经过模型和尾侧区积雪之间的狭窄通道时，由于狭管效应风速加大，当模型的阻挡作用消失后，风速急剧下降，此时加上绕过模型的马蹄涡产生的风速回流形成更强的减速区，不能带动雪颗粒运动，所以在狭管出风口处形成较大范围的积雪，等到迎风区两侧区域积雪堆积后，狭管入风口被积雪阻挡，狭管效应消失，狭管出风口减速区效果减弱，有可能在发生少量侵蚀后保持稳定，如图12所示。
3 结论
通过对立方体周边风致积雪演化过程的现场测量，揭示了风速及持续时间对立方体周边积雪堆积的影响，为数值模拟或风洞试验验证提供了参考。根据对实测结果的对比分析，可以得出以下结论：
(1) 在无降雪且雪源充足的风吹雪环境下，立方体周边主要在迎风区、背风区和两尾侧区域出现积雪，其中最大积雪深度出现在迎风区壁面附近或者两侧区域。风速和持续时间对立方体周边积雪的深度和分布范围都有明显的影响；
(2) 风速较大时，迎风前壁积雪深度增长较快，尾侧区积雪范围更长，背风区靠近壁面处积雪堆积区域较小，积雪深度也较小，背风区后侧积雪范围增长，积雪深度增大。较大风速对应的积雪最大深度在迎风区和尾侧区也比低风速下的对应值更大，且迎风区积雪最大深度位置会更接近模型；
(3) 随着时间的推移，积雪最大深度位置会从背风区转移到迎风区，然后稳定在迎风区，尾侧区积雪范围在初期不变(本测试中的8 min以内)，随着时间的增长，积雪深度逐渐增大，积雪范围也进一步增大，积雪深度最大位置随积雪范围增大而向远离模型的方向移动；
(4) 立方体两侧的积雪会引起流场变化从而对立方体周边积雪分布造成影响，较大风速时，迎风区两侧不易出现积雪。迎风区两侧无积雪时，随着时间推移，迎风侧前壁积雪深度逐渐加大，背风区模型后方积雪深度不断加深，积雪范围也比较大；迎风区两侧有积雪时，随着时间的推移，迎风侧前壁积雪稳定，背风区模型后方在发生少量侵蚀后保持稳定。
-
![]()
图 9 11.4 m/s风速下尾侧区积雪对比
Figure 9. Comparison of snow distribution in the lateral region at 11.4 m/s
表 1 雪的堆积密度
Table 1 Bulk density of snow
雪种类 密度/(g/cm3) 新雪 0.03~0.12 >24 h 0.18 >48 h 0.29 >72 h 0.30 沉积雪 0.29~0.42 
下载: 导出CSV
表 2 雪的休止角
Table 2 The angle of repose
雪种类 休止角/(°) 新雪 60.78 >24 h 58.14 >30 h 56.97 表层积雪 40.54 底层积雪 32.03
下载: 导出CSV
表 3 雪通量
Table 3 Mass flux
高度/cm 6.7 m/s雪通量/(g·m−2·s−1) 11.4 m/s雪通量/(g·m−2·s−1) 2.5 68.75 127.78 7.5 10.42 23.61 12.5 4.17 11.81 17.5 2.08 6.94 22.5 2.08 6.25 27.5 0.00 2.78
下载: 导出CSV
表 4 风速测量数据
Table 4 Wind speed measurement data
工况 1 min 2 min 4 min 8 min 10 min 工况1：平均风速v/(m/s) 6.2 6.8 7.1 6.8 6.8 工况2：平均风速 v/(m/s) 11.0 10.8 12.1 12.1 10.8
下载: 导出CSV
表 5 积雪最大深度及其位置
Table 5 Maximum snow depth and location
堆积区域 平均风速v/(m/s) 时间t/min X/B Y/B D/B 迎风区 6.7 2 −1.30 0.17 0.05 4 −1.30 0.17 0.11 8 −1.30 0.00 0.16 11.4 2 −1.10 0.00 0.08 4 −1.10 0.00 0.13 8 −1.10 0.17 0.18 尾侧区 6.7 2 1.83 1.83 0.05 4 1.83 1.83 0.10 8 1.83 −1.17 0.14 11.4 2 2.17 1.50 0.07 4 1.83 1.50 0.10 8 2.17 1.50 0.16 背风区 6.7 2 0.50 0.00 0.08 4 0.50 0.00 0.09 8 0.50 0.00 0.09 11.4 2 0.50 0.17 0.08 4 0.50 0.17 0.08 8 0.50 0.17 0.07 注：X/B为顺风向位置； Y/B为横风向位置； D/B为积雪深度。
下载: 导出CSV
-
[1] Kind R J. A critical examination of the requirements for model simulation of wind-induced erosion/deposition phenomena such as snow drifting [J]. Atmospheric Environment, 1976, 10(3): 219 − 227. doi: 10.1016/0004-6981(76)90094-9
[2] Smedley D J, Kwok K C S, Kim D H. Snowdrifting simulation around davis station workshop, Antarctica [J]. Journal of Wind Engineering Industrial Aerodynamics, 1993, 50: 153 − 162.
[3] Nemoto M, Sato T, Kosugi K, et al. Effects of snowfall on drifting snow and wind structure near a surface [J]. Boundary-Layer Meteorology, 2014, 152(3): 395 − 410. doi: 10.1007/s10546-014-9924-4
[4] Liu M, Zhang Q, Fan F, et al. Experiments on natural snow distribution around simplified building models based on open air snow-wind combined experimental facility [J]. Journal of Wind Engineering & Industrial Aerodynamics, 2018, 173: 1 − 13.
[5] Nikolas O A, John W P. The effect of coherent structures in the atmospheric surface layer on blowing-snow transport [J]. Boundary-Layer Meteorology, 2018, 167(2): 211 − 233.
[6] Ding Minghu, Zhang Tong, Xiao Cunde, et al. Snowdrift effect on snow deposition: Insights from a comparison of a snow pit profile and meteorological observations in east Antarctica [J]. Science China Earth Sciences, 2017(4): 672 − 685.
[7] Yan Min, Zuo Hejun, Wang Haibing, et al. Snow resisting capacity of caragana microphylla and achnatherum splendens in a typical steppe region of inner mongolia, China [J]. Journal of Arid Land, 2020, 12(2): 294 − 302. doi: 10.1007/s40333-019-0021-x
[8] Liston G E, Haehnel R B, Sturm M, et al. Simulating complex snow distributions in windy environments using snowtran-3D [J]. Journal of Glaciology, 2007, 181: 241 − 256.
[9] Tominaga Y. Computational fluid dynamics simulation of snowdrift around buildings: Past achievements and future perspectives [J]. Cold Regions Ence and Technology, 2018, 150: 2 − 14.
[10] 周晅毅, 顾明. 风致积雪漂移堆积效应的研究进展[J]. 工程力学, 2008(7): 5 − 10, 17. Zhou Xuanyi, Gu Ming. Simulation of the wind-induced snowdrift state of the art [J]. Engineering Mechanics, 2008(7): 5 − 10, 17. (in Chinese)
[11] 周晅毅, 张运清, 顾明. 建筑屋面滑移雪荷载的模拟方法研究[J]. 工程力学, 2014, 31(6): 190 − 196. doi: 10.6052/j.issn.1000-4750.2012.12.1023 Zhou Xuanyi, Zhang Yunqing, Gu Ming. Method for simulating sliding snow loads on roofs [J]. Engineering Mechanics, 2014, 31(6): 190 − 196. (in Chinese) doi: 10.6052/j.issn.1000-4750.2012.12.1023
[12] 周晅毅, 刘长卿, 顾明, 谭敏海. 拉格朗日方法在风雪运动模拟中的应用[J]. 工程力学, 2015, 32(1): 36 − 42. doi: 10.6052/j.issn.1000-4750.2013.07.0607 Zhou Xuanyi, Liu Changqing, Gu Ming, Tan Minhai. Application of lagrangian method to snowdrift model [J]. Engineering Mechanics, 2015, 32(1): 36 − 42. (in Chinese) doi: 10.6052/j.issn.1000-4750.2013.07.0607
[13] 周晅毅, 胡学富, 顾明. 考虑摩擦速度修正的风致积雪重分布模拟方法研究[J]. 工程力学, 2017, 34(10): 19 − 25. doi: 10.6052/j.issn.1000-4750.2016.06.0432 Zhou Xuanyi, Hu Xuefu, Gu Ming. Study on wind-induced snow redistribution considering modification of friction velocity [J]. Engineering Mechanics, 2017, 34(10): 19 − 25. (in Chinese) doi: 10.6052/j.issn.1000-4750.2016.06.0432
[14] Jiang X, Yin Z, Cui H. Wind Tunnel tests and numerical simulations of wind-induced snow drift in an open stadium and gymnasium [J]. Advances in Civil Engineering, 2020(3): 1 − 15.
[15] 孙晓颖, 洪财滨, 武岳, 等. 建筑物周边风致雪漂移的数值模拟研究[J]. 工程力学, 2014(4): 141 − 146. doi: 10.6052/j.issn.1000-4750.2012.11.0839 Sun Xiaoying, Hong Caibin, Wu Yue, et al. Numerical simulation of snow drifting around building model [J]. Engineering Mechanics, 2014(4): 141 − 146. (in Chinese) doi: 10.6052/j.issn.1000-4750.2012.11.0839
[16] 徐枫, 周高照, 肖仪清, 等. 建筑屋盖表面及其周围积雪分布研究[J]. 振动与冲击, 2018, 37(13): 92 − 99. Xu Feng, Zhou Gaozhao, Xiao Yiqing, et al. Snow distribution on surface of building roof and its surrounding [J]. Journal of Vibration and Shock, 2018, 37(13): 92 − 99. (in Chinese)
[17] 刘多特, 李永乐, 汪斌. 风雪绕流数值模拟的积雪预测模型研究[J]. 工程力学, 2016, 33(8): 122 − 131. doi: 10.6052/j.issn.1000-4750.2015.01.0019 Liu Duote, Li Yongle, Wang Bin. A numerical prediction model for snow accumulation caused by ambient snowdrift [J]. Engineering Mechanics, 2016, 33(8): 122 − 131. (in Chinese) doi: 10.6052/j.issn.1000-4750.2015.01.0019
[18] 刘多特, 李永乐. 风雪流滞空颗粒受力及运动行为的数值模拟[J]. 工程力学, 2018, 35(12): 15 − 24. doi: 10.6052/j.issn.1000-4750.2017.09.0672 Liu Duote, Li Yongle. A numerical simulation on aerodynamic forces and motion behaviors of snow granules in air [J]. Engineering Mechanics, 2018, 35(12): 15 − 24. (in Chinese) doi: 10.6052/j.issn.1000-4750.2017.09.0672
[19] Oikawa S, Tomabechi T. Formation processes of the deposition and erosion of snow around a model building [J]. Seppyo, 2003, 65(3): 207 − 218.
[20] Tsuchiya M, Tomabechi T, Hongo T, et al. Wind effects on snowdrift on stepped flat roofs [J]. Journal of Wind Engineering & Industrial Aerodynamics, 2002, 90(12): 1881 − 1892.
[21] Thiis T K. Large scale studies of development of snowdrifts around buildings [J]. Journal of Wind Engineering & Industrial Aerodynamics, 2003, 91(6): 829 − 839.
[22] Schmidt R A. Properties of blowing snow [J]. Reviews of Geophysics and Space Physics, 1982, 20: 39 − 44. doi: 10.1029/RG020i001p00039
[23] Kind R J. Snowdrifting: A review of modelling methods [J]. Cold Regions ence & Technology, 1986, 12(3): 217 − 228.
[24] Pomeroy J W, Gray D M. Saltation of snow [J]. Water Resources Research, 1990, 26(7): 1583 − 1594. doi: 10.1029/WR026i007p01583
[25] Schmidt R A. Transport of drifting snow and the mean wind speed profile [J]. Boundary-Layer Meteorology, 1986, 34(3): 213 − 241. doi: 10.1007/BF00122380
[26] Oikawa S, Tomabechi T, Ishihara T. One-day observation of snowdrifts around a model cube (in Japanese) [J]. Journal of snow Engineering of Japan, 1999, 15: 283 − 291. doi: 10.4106/jsse.15.4_283
[27] 刘多特. 风致积雪漂移沉降特性及线路与桥梁积雪预测研究[D]. 成都: 西南交通大学, 2017. Liu Duote. Study on characteristics of snowdrift and redistribution prediction for snowcover on traffic line and bridge [D]. Chengdu: Southwest Jiaotong University, 2017. (in Chinese)
[28] 赵雷, 余志祥, 齐欣, 等. 低矮建筑屋盖风雪流作用场地实测与数值模拟[J]. 振动与冲击, 2017(22): 225 − 231. Zhao Lei, Yu Zhixiang, Qi Xin, et al. Field measurement and numerical simulation of snowdrift on low-rise buildings [J]. Journal of Vibration and Shock, 2017(22): 225 − 231. (in Chinese)
[29] Beyers J H M, Sundsb P A, Harms T M. Numerical simulation of three-dimensional, transient snow drifting around a cube [J]. Journal of Wind Engineering & Industrial Aerodynamics, 2004, 92(9): 725 − 747.
[30] Tominaga Y, Okaze T, Mochida A. CFD modeling of snowdrift around a building: An overview of models and evaluation of a new approach [J]. Building & Environment, 2011, 46(4): 899 − 910.
[31] 李雪峰. 风致建筑屋盖表面及其周边积雪分布研究[D]. 上海: 同济大学, 2011. Li Xuefeng. Research on snowdrifting on building roof and around building [D]. Shanghai: Tongji University, 2011. (in Chinese)
[32] 赵有栓, 颜卫亨, 姬明辉. 风致积雪运动在立方体周围表现的数值模拟 [C]. 天津: 天津大学, 天津市钢结构学会: 全国现代结构工程学术研讨会学术委员会, 2017. Zhao Youshuan, Yan Weiheng, Ji Minghui. Numerical simulation of wind-induced snow movement around a cube [C]. Tianjin: Tianjin University, Tianjin Steel Structure Society: Academic Committee of the National Symposium on Modern Structural Engineering, 2017. (in Chinese)
[33] Michael O'Rourke, Arthur De Gaetano, Janine Derjue Tokarczyk. Snow drifting transport rates from water flume simulation [J]. Journal of Wind Engineering & Industrial Aerodynamics, 2004, 92(14/15): 1245 − 1264.
[34] 张玉杰, 颜卫亨, 张雯. 轻型四坡房屋表面风致雪漂的数值模拟[C]. 天津: 天津大学, 天津市钢结构学会: 全国现代结构工程学术研讨会学术委员会, 2018. Zhang Yujie, Yan Weiheng, Zhang Wen. Numerical simulation of snow drifting on light hipped roof [C]. Tianjin: Tianjin University, Tianjin Steel Structure Society: Academic Committee of the National Symposium on Modern Structural Engineering, 2018. (in Chinese)
[35] 杨旭. 玛依塔斯交通走廊风吹雪特点研究[J]. 铁道工程学报, 2018, 35(12): 1 − 6. Yang Xu. Research on the characteristics of drifting snow disaster in Xinjiang Mayitas traffic corridor [J]. Journal of Railway Engineering Society, 2018, 35(12): 1 − 6. (in Chinese)
-
期刊类型引用(4)
1. 辛林桂，周晅毅，顾明. 平屋面摩擦速度的数值模拟研究. 工程力学. 2024(02): 151-159 . 
本站查看
2. 杨树瑚，王斌，张云，洪中华，韩彦岭，王静. 地形信息对利用Sentinel-1数据反演南北极积雪密度的影响. 测绘科学. 2024(12): 114-123 . 百度学术
3. 孟少伟，冯涛，张广泽，徐正宣，刘建国，林之恒，宋章，伊小娟，欧阳吉，谢荣强. 川藏铁路雅安至昌都段风吹雪空间分异特征及易发性评价. 地球科学. 2022(03): 987-994 . 百度学术
4. 骆颜，马文勇，孙元春. 高寒地区交通线路风速及风向预测研究. 工程力学. 2022(S1): 195-201 . 本站查看
其他类型引用(8)
计量
- 文章访问数: 474
- HTML全文浏览量: 157
- PDF下载量: 44
- 被引次数: 12
