风沙流环境地表剪切力测量.pdf

文章编号：1672-9897(2023)05-0093-08doi：10.11729/syltlx20230074风沙流环境地表剪切力测量黄宁1,2，杨博1,2，何畏3，张洁1,2,*1.兰州大学土木工程与力学学院，兰州7300132.兰州大学西部灾害与环境力学教育部重点实验室，兰州7300133.温州理工学院建筑与能源工程学院，温州325035摘要：地表剪切力是风沙运动的关键驱动力，在风沙运动研究中非常重要，但目前针对风沙流环境地表剪切力的实验测量却较为匮乏，已经成为风沙运动研究深入发展的瓶颈。本文介绍了一种基于热膜技术的适用于风沙流环境地表剪切力精确测量的方法，使用热膜传感器和二维热线探头对净风条件下的地表剪切力进行了测量，结果显示：热膜传感器对地表剪切力有较好的响应，其与二维热线探头对摩阻风速的测量结果相差小于 6.6%。热膜传感器对方柱绕流条件下地表剪切力的测量结果与 OpenFOAM 数值模拟结果具有较好的一致性，进一步证明了热膜传感器测量地表剪切力的可行性。热膜传感器的沙粒撞击实验结果表明：沙粒撞击对热膜信号的影响可以忽略，在风沙流环境中可以采用热膜传感器对地表剪切力进行有效测量。基于该技术测量了风沙流环境地表剪切力沿流向的空间分布特征，采用实测方式获得了风沙流环境地表剪切力的空间分布规律。关键词：热膜传感器；风沙流；地表剪切力；二维热线探头中图分类号：O359；O368文献标识码：AMeasurement of wall shear stress in wind-blown sand environmentHUANG Ning1,2，YANG Bo1,2，HE Wei3，ZHANG Jie1,2,*1.College of Civil Engineering and Mechanics,Lanzhou University,Lanzhou730013,China2.Key Laboratory of Mechanics on Disaster and Environment in Western China,The Ministry of Education of China,Lanzhou University,Lanzhou730013,China3.College of Architecture and Energy Engineering,Wenzhou University of Technology,Wenzhou325035,ChinaAbstract:Asakeydrivingforceofaeoliansandmovement,wallshearstressisofimportanceinthestudyofaeoliansandmovement.However,atpresent,theexperimentalmeasurementofthe surface shear force in the wind-blown sand environment is extremely scarce,which hasbecomeabottleneckforthefurtherdevelopmentofthewind-blownsandresearch.Thispaperintroducesamethodbasedonthehot-filmtechnologyforaccuratemeasurementofthewallshearstressinthewind-blownsandenvironment.Thewallshearstresswasmeasuredunderconditionswithoutsandusinghot-filmprobesandtwo-dimensionalhot-wireprobes,respectively.Theresultshowsthatthehot-filmprobehasagoodresponsetothewallshearstress,andthedifferencebetweenthemeasurementresultsofthefrictionvelocityofthesetwoprobesislessthan6.6%.Themeasurementresultsofthehot-filmprobeforthewallshearstressofflowaroundthesquarecolumnareingoodagreementwiththeOpenFOAMnumericalsimulationresults,whichfurtherprovestheeffectivenessofthehot-filmprobeforthemeasurementofthewallshearstress.Inaddition,theparticleimpacttestresultsofthehotfilmprobesshowthattheinfluenceofparticleimpact on the hot-film signal is negligible,so the hot-film probe can be used to effectively收稿日期：2023-05-22；修回日期：2023-07-18；录用日期：2023-07-24*通信作者E-mail：zhang-引用格式：黄宁,杨博,何畏,等.风沙流环境地表剪切力测量 J.实验流体力学,2023,37(5):93-100.HUANG N,YANG B,HE W,et al.Measurement of wall shear stress in wind-blown sand environmentJ.Journal ofExperimentsinFluidMechanics,2023,37(5):93-100.第37卷第5期实验流体力学Vol.37，No.52023年10月Journal of Experiments in Fluid MechanicsOct.，2023measurethewallshearstressinthewind-blownsandenvironment.Basedonthistechnology,thespatialdistributioncharacteristicsofthewallshearstressalongtheflowdirectioninthewind-blownsandenvironmentweretested.Andforthefirsttime,thespatialdistributionofthewallshearstressinthewind-blownsandenvironmentwasobtainedbymeansofactualmeasurement.Inviewofthehigh-frequencyresponsecapabilityofthehot-filmprobe,thistechnologywillplayanimportantroleintheexperimentalresearchonthemechanismofaeoliansandmovementintheturbulentboundarylayer.Keywords:hot-filmprobe;wind-blownsand;wallshearstress;two-dimensionalhot-wireprobe 0 引言 风沙运动是一种对人类生活影响巨大的极端环境，不但会破坏建筑、道路等人工设施，还会严重影响人类健康，每年都会造成巨大的经济损失1。在风沙物理学中，风沙运动可分为沙粒释放、沙粒输运、沙粒沉降等 3 个力学过程2，这些过程都与地表剪切力 联系密切（在流体力学中，地表剪切力通常表示为摩阻风速 u*的函数：=u*2）。在纯流体中，研究者已经针对不同尺度湍流结构、雷诺数及地表粗糙度条件下的地表剪切力平均值和脉动值变化情况开展了充分的理论和实验研究3-4；而在气固两相流中，由于测量技术的限制，关于地表剪切力的研究主要针对平均值，对脉动值的研究则相对较少5-8。近年来，研究者开始认识到地表剪切力脉动值对风沙流中沙粒的影响不能忽略9-10。在沙粒起动方面，通过风沙流环境野外实验研究发现：真实风沙运动中的沙粒起动具有间歇性11-12，单纯的地表剪切力平均值不能很好地反映间歇性现象，对输沙率预测会产生影响。在沙粒输运方面，研究者对跃移输沙通量（工程与科学研究共同关注的重要宏观物理量）进行了深入研究，针对不同工况及环境得到了一系列输沙率公式13-15。通过对比研究发现，输沙率与摩阻风速的三次方成正比，两者之间存在密切的联系。对于这种非线性关系，平均值显然不能很好地反映输沙率的变化。另外，在野外环境下，由温度的日变化等带来的风速剧烈变化会导致地表剪切力平均值的研究价值降低。例如，对于地表加热不均匀而产生的尘卷风，地表剪切力平均值就无法真实反映其内在机理。在沙粒沉降方面，沉降末速度是一个重要参数，会对沉积通量产生影响。通过对几种沉积方案的研究发现：地表剪切力对沉降末速度存在重要影响16-17。同样，由于沉降末速度与摩阻风速之间的非线性关系，基于地表剪切力平均值的研究可能不足以实现对沉积通量的准确估计。测量风沙流中的地表剪切力脉动时间序列，对于风沙物理研究具有重要意义。研究者针对风沙流环境地表剪切力已经进行了大量测量。在环境复杂恶劣的野外实验中，一般使用三维超声风速仪测量风速，进而推导出地表剪切力。而在风洞实验中，由于空间分辨率不足，无法使用超声风速仪精确测量地表剪切力脉动时间序列。针对这一问题，Irwin 模仿普雷斯顿管式表面摩擦计18发明了 Irwin 探头19，其测量基于垂直于地表的探头细管顶端与地表之间压差是摩阻风速的指数函数的原理，通过测量压差即可得到摩阻风速。在测量风沙流环境地表剪切力平均值的实验中，Irwin 探头的表现令人满意；但由于其测量的是近地表不同高度处通过气体传播的压差，压力信号存在延迟，使得 Irwin 探头对地表剪切力的瞬时响应速度较低。经过多年发展，基于电信号的测量手段逐渐发展成熟。其中，热膜传感器的测量精度得到明显提升，已经广泛应用于飞行器、潜水器的表面阻力测量20-24。带有保护层的热膜传感器抗干扰能力更强，对沙粒撞击具有一定抵挡作用，在气固两相流地表剪切力测量中应用潜力较大。本文通过实验验证热膜地表剪切力测试仪的测量精度，评估其在风沙流环境下的表现及风沙流对热膜传感器的影响。通过风洞实验，测量方柱绕流分离区内的地表剪切力空间分布情况，并与数值模拟结果进行对比，验证热膜传感器测量复杂流场环境下地表剪切力的准确性。1 测量原理简介以热膜传感器测量地表剪切力，是通过测量壁面与流体之间的热传导实现的。使用热膜传感器需满足 3 个前提条件：1）热膜传感器足够薄且紧贴被测对象表面，其对流场的影响可以忽略；2）工作温度低于 100，可忽略辐射传热，且不会对流体性质及流动特性造成影响；3）热膜传感器的发热单元在流向上足够短，其产生的温度边界层厚度 T远小于速94实验流体力学http:/度边界层厚度。本文使用的热膜传感器以聚酰亚胺箔（PI）作为基板，将 1 m 热敏镍（Ni）磁控溅射于 PI 箔上，再采用光刻及湿法蚀刻工艺将传感器元件进行图案化。传感元件长 3 mm，宽 50 m，热膜总厚度不超过 80 m。该传感器的有效响应频率可达到 4 kHz，结合风沙流研究特点及实际测量环境，将测试频率设置为 2 kHz，该频率已足以反映测量过程中风沙流环境下地表剪切力的相关信息。热膜装置如图 1 所示。（a）柔性热膜电阻丝（b）热膜地表剪切力测试仪主机图 1 热膜装置Fig.1 Illustration of hot-film 1.1 热膜传感器标定关系以二维定常流为例推导标定关系式。式（1）和（2）分别为层流边界层微分方程和能量积分方程：uux+vy=dpdx+2uy2（1）ddxw0u(T T)dy=kCp(Ty)w（2）式中：k 为热传导系数；Cp为定压比热；T 为传感器温度，T为来流温度；p、分别为流体压力和密度；u、v 分别为流向速度、垂向速度；、分别为运动黏性系数、动力黏性系数；下标 w 表示物理量在壁面处的值。由前述第 3 个前提条件（即 T ），可将式（1）简化为：2uy2=dpdx（3）对式（3）积分，代入边界条件 y=0 和 =w（w表示高度为 0 时地表剪切力的边界条件）可得：u=wy+12dpdxy2（4）将式（4）代入式（2）可得：ddxww0y(T T)dy+12dpdxw0y2(T T)dy=2kPrqw(x)（5）式中，Prandtl 数 Pr=Cp/k，热膜单位面积热传导qw(x)=k(T/y)w率。式中积分可用 Curle 法求解。假设：T T=(Tw T)f()（6）式中，Tw为 y=0 时的温度值。由 T 可得：f()=(1 43+34)+(33+24)（7）其 中，高 度 的 无 量 纲 值 =y/T和 =qwT/k(Tw T)为温度边界层的形式参数。由式（6）和（7）可得式（5）的第一个积分结果为：w0y(T T)dy=w102T(Tw T)f()d=ak2(Tw T)3/q3w(x)（8）a=2w102f()d式中，。同理可得第二个积分结果为：w0y2(T T)dy=w1023T(Tw T)f()d=bk3(Tw T)4/q3w(x)（9）b=3w102f（）d式 中，。将 上 述 积 分 结 果 代 入式（5）可得：ddxak2w(Tw T)3q2w(x)bk3(Tw T)42q3w(x)dpdx=2kPrqw(x)（10）假设 Tw和 qw为 x 的阶跃函数，可以得到：ak2w(TQw)312bdpdxk3(TQw)4=2LkPr（11）式中：L 为热膜的有效长度；Qw为热膜与流体之间传递的总热量；T 为来流温度与壁面温度差值。假设静压梯度很小，即 dp/dx 0，可将式（11）简化为：1/3w=(2Lak3Pr)1/3QwT（12）式中，Qw等于热膜电阻发热产生的热能 I2R 减去壁面热损失 Q，则式（12）可写作：1/3w=(2Lak3Pr)1/3I2RT(2Lak3Pr)1/3QT（13）以上推导过程中使用了若干假设，且风洞实验中难以测量地表的热损失，因此通过式（13）无法直接得出地表剪切力，需通过实验进行标定。根据式（13），假设标定关系为如下形式：1/3w=A1I2RT+B（14）式中，A1和 B 为标定时需确定的参数。实际上，热膜地表剪切力测试仪的工作原理与恒温式热线风速仪相同，通过测量电桥电压得到 w，此时 Tw为恒定第 5 期黄宁等：风沙流环境地表剪切力测量95值。若来流温度 T的变化忽略不计，则式（14）可写作：1/3w=A2E2+B（15）式中：A2为标定参数；E 为电桥电压，即式（14）中I2R。在实际测量中，一般采用如下形式：E20 E2w=A1/3w+B（16）式中，E0为无风状态下电桥电压读数，A 和 B 为标定参数。2 热膜地表剪切力测试仪的标定和实验验证在航天飞行器、潜艇等测量环境下，热膜地表剪切力测试仪的准确性已经得到证实20-24，但在风沙流环境的恶劣条件下，其准确性仍需进一步验证。本文对热膜传感器进行了标定，通过风洞实验验证其测量结果的准确性，同时验证其在承受沙粒撞击过程中工作的稳定性。2.1 实验布置如图 2 所示，在风洞实验段铺设木板，形成光滑平坦地表，以进行标定和测试实验。在地表上方相同高度布置皮托管和二维热线探头，并将两者安装于移动坐标架上，以测量不同高度处的流动信息。在地表沿流向布置 3 个热膜传感器，探头间距 5 cm。对皮托管、二维热线探头和热膜传感器进行同步观测，热膜传感器的采样频率为 2 kHz，采样时间为3 min。皮托管热膜二维热线探头图 2 实验布置图Fig.2 Experimental setup 2.2 热膜传感器与热线探头测量结果对比在风沙流环境下，对仪器抗干扰能力有着较高要求，因此需对热膜传感器的测量精度进行验证，并评估风沙中运动沙粒对热膜测量结果的影响程度。实验中，以皮托管测量风速廓线，并通过壁面率计算出 5 种不同风速下的地表剪切力，与热膜传感器测量结果进行对比，最终完成标定。图 3 为实验中 3 个热膜传感器（No.1 No.3）在 5 种不同风速下的标定结果，同时还使用式（16）进行了线性拟合。0.0020.0030.0040.0050.0060.007E2 E020.30.40.50.6w1/3No.1No.2No.3图 3 不同风速下热膜传感器的标定结果Fig.3 Calibration results of hot-film probes at different velociiesuw完成标定之后，还需对热膜传感器的测量结果进行验证。图 4（a）为二维热线探头测量的不同高度处的雷诺应力（u、w分别为流向瞬时脉动风速、垂向瞬时脉动风速；图中 u为来流速度）。随着200150100500z/mm00.050.100.150.20（a）二维热线探头测得的雷诺应力廓线uw/(m2s2)0.50.40.30.20.1u*/(ms1)48126101416（b）热膜传感器与二维热线探头测量结果对比u/(ms1)u=6 m/su=8 m/su=10 m/su=12 m/su=14 m/sHot-film probeHot-wire probe图 4 热膜传感器与二维热线探头测量结果对比Fig.4 Comparison of hot-wire and hot-film probe96实验流体力学http:/uw高度 z 增大，雷诺应力先增大后减小，摩阻风速由的峰值确定25。图 4（b）为热膜传感器与二维热线探头的摩阻风速 u*测量结果的对比，从图中可以看出，两者之间的测量误差在 0.7%6.6%之间。E20E2w由图 4 可知，热膜传感器对地表剪切力的测量结果是准确的。进一步地，对其在风沙流环境下的工作表现进行验证。为评估沙粒撞击对热膜传感器测量结果的影响程度，在无风状态下以沙粒撞击热膜传感器。图 5（a）为沙粒撞击产生的噪声信号与实际风场测量中输出信号的对比（图中 e0=），可以发现，沙粒撞击对信号的影响远小于信号本身的波动。图 5（b）为 3 种风速下沙粒撞击产生的噪声信号与实际风场测量中输出信号的比值（图中 ew=）。经对比发现，沙粒撞击产生的噪声信号比实际测量的信号低一个数量级。因此，在风沙流环境中，沙粒撞击引入的信号噪声相当小。0.5660.5640.5620.5600.5580.556e0024810t/s（a）沙粒撞击热膜传感器产生的噪声101102103104(e0 ew)/e00.5750.5800.5850.5900.595e0（b）噪声信号与输出信号的比值u=0 m/su=0 m/su=9.28 m/sSand hombardmentu=0 m/s,with bombardmentIn airflow,without bombardment(1)(2)(3)(1)u=9.28 m/s(2)u=13.26 m/s(3)u=17.41 m/s图 5 沙粒撞击对热膜传感器测量结果的影响Fig.5 Effect of particle impact on measurement results of hot-filmprobe 3 热膜测量结果与数值模拟对比为进一步验证热膜传感器的测量精度，本节以方柱绕流算例结果与热膜传感器测量结果进行对比。方柱扰流实验设置如图 6 所示，在热膜传感器上游布置一个方柱（5 cm 5 cm 10 cm）。实验中，改变方柱与热膜传感器之间的距离，在方柱的每个位置上测量 3 min，对测量结果进行时间平均，即可得出方柱后方地表剪切力平均值随方柱与热膜传感器之间距离的变化情况。Air flow 图 6 方柱绕流实验布置图Fig.6 Experimental setup of flow around a square column针对上述实验工况，基于 OpenFOAM 开源代码中求解不可压缩流体的瞬态求解器进行了数值模拟。数值模拟的计算域和网格划分如图 7 所示，计算域设置为 1.0 m 1.0 m 2.5 m，在地表和方柱附近进行了网格加密。设置入口风速 u=12 m/s，计算时长为 25 s。湍流模型为 IDDES（k Shear StressTransport (SST)Improved Delayed DetachedEddy Simulation），该模型被广泛应用于复杂流场的湍流计算，准确性已经得到验证。图 7 数值模拟网格Fig.7 Grid division of numerical simulation在求解湍流过程中，近壁面区域使用 RANS 模型进行计算，远离壁面区域则使用 LES 模型，以尽量确保数值模拟的准确性。图 8 为方柱绕流的地表剪切力云图（x、y 分别为流向和展向坐标），从图中可以清晰分辨出方柱后方的湍流结构。数值模拟能够较为精确地反映方柱后方的地表剪切力变化情况，从而对风洞实验结果进行验证。从图 9 可以看到，通过实验得到的方柱后方不同距离（s）处的摩阻风速和数值模拟结果之间能够较好地匹配，这从另一个角度证明了热膜传感器的第 5 期黄宁等：风沙流环境地表剪切力测量97准确性。y/m 0.40.200.20.40.500.51.01.52.0 x/m 1.051 1030.038 0.075 0.112 0.151Wall shear stress/(Nm2)图 8 数值模拟中的地表剪切力分布Fig.8 Spatial distribution of wall shear stress in numerical simulation0.30.20.100.1u*/(ms1)020406080100s/cm50 mm 50 mm 100 mmExperiment resultSimulation result图 9 实验得到的方柱绕流摩阻风速与数值模拟结果对比Fig.9 Comparison between wind tunnel experiment and numericalsimulation results of the friction velocity of flow behind asquare column 4 热膜地表剪切力测试仪在风洞实验中的应用前文介绍了热膜地表剪切力测试仪的工作原理，对其准确性和抗干扰能力进行了分析，证明了其在风沙流环境下测量地表剪切力脉动值的可行性。本节采用热膜地表剪切力测试仪，在风洞实验中对风沙流不同发展阶段的地表剪切力进行测量。图 10 为风洞实验示意图。风洞实验段长度为22 m，工作截面尺寸为 1.30 m 1.45 m。在风洞中，沿流向铺设了宽 0.5 m、厚 0.08 m 的沙床，沿风洞轴线布置了宽为 3 cm 的金属方管，方管上表面与沙面平齐。在方管上表面沿流向粘贴 11 个热膜传感器，对风沙流中的流体剪切力 a进行测量（在风沙流中，地表剪切力 分为沙粒剪切力 p和流体剪切力 a，此处测量的地表剪切力为流体剪切力 a）。实验风速设定为 12.3 m/s。图 11 给出了风沙流发展过程中不同流向位置的地表剪切力，其中 ft表示使沙粒离开沙床（地表）进入空中所需的地表剪切力流体起动临界值，该值可ft=au2ft由式得出，其中 u*ft由下式确定26：uft=ANp aagdp+adp（17）式中：实验参数 AN=0.111；修正参数26=2.9 104 N/m；a为空气密度，dp、p分别为沙粒粒径和密度。u=12.3 m/sNatural sandHot-film probes图 10 风沙流环境地表剪切力测量风洞实验示意图Fig.10 Schematic of measurement of wall shear stress in wind-blown environment0.30.20.1a/Pa0246810 x/mu=12.3 m/s ft图 11 地表剪切力随流向位置的变化Fig.11 Variation of wall shear stress with flow direction position从图 11 可以看出：在沙床的起始位置（x 3 m；沙床前缘起始处为 x=0 m），地表剪切力的平均值高于地表剪切力流体起动临界值 ft，沙粒起动以流体起动为主；随着风沙流的发展，进入空中和落回沙床的沙粒越来越多，沙粒落地冲击地表，沙粒起动变为以冲击起动为主，此时地表剪切力平均值小于地表剪切力流体起动临界值，这与目前风沙物理学中的普遍认知相同。通过上述实验，证明热膜传感器可以对风沙流环境下的地表剪切力进行测量。5 结论本文介绍了热膜传感器的工作原理，推导了标定公式，对其测量风沙流环境地表剪切力的可行性进行了验证，提出了在风洞实验中采用热膜地表剪切力测试仪测量风沙流环境地表剪切力脉动值的技术。通过风洞实验和数值模拟，验证了热膜地表剪切力测试仪的准确性，对沙粒撞击引入测量的噪声信号进行了评估。实验结果显示，热膜传感器与二98实验流体力学http:/维热线探头测量的摩阻风速相差不超过 6.6%，风沙流中沙粒撞击产生的噪声信号比实际测量输出信号低一个数量级，数值模拟结果与热膜传感器测量结果具有较好的一致性。在风洞实验中实现了风沙流环境地表剪切力的实际测量，得出了风沙流环境下地表剪切力沿流向的空间分布规律。以上结果表明，热膜地表剪切力测试仪是一种可靠的、可用于风沙运动研究的地表剪切力测量仪器。参考文献：DUPONTS,BERGAMETTIG,MARTICORENAB,etal.ModelingsaltationintermittencyJ.JournalofGeophysicalResearch:Atmospheres,2013,118(13):71097128.doi：10.1002/jgrd.505281BAGNOLD R A.The physics of blown sand and desertdunesM.NewYork:W.Morrow&company,1941.2AGOSTINI L,LESCHZINER M.The connection betweenthe spectrum of turbulent scales and the skin-frictionstatisticsinchannelflowatRe1000J.JournalofFluidMechanics,2019,871:2251.doi：10.1017/jfm.2019.2973DUANYC,ZHONGQ,WANGGQ,etal.Contributionsof different scales of turbulent motions to the mean wall-shear stress in open channel flows at low-to-moderateReynoldsnumbersJ.JournalofFluidMechanics,2021,918:A40.doi：10.1017/jfm.2021.2364JACOBC,ANDERSONW.Conditionallyaveragedlarge-scale motions in the neutral atmospheric boundary layer:insights for aeolian processesJ.Boundary-Layer Meteo-rology,2017,162(1):2141.doi：10.1007/s10546-016-0183-45RICHTERDH,SULLIVANPP.Modificationofnear-wallcoherentstructuresbyinertialparticlesJ.PhysicsofFluids,2014,26(10):103304.doi：10.1063/1.49005836RANA S,ANDERSON W,DAY M.Turbulence-basedmodel for subthreshold aeolian saltationJ.GeophysicalResearchLetters,2020,47(15):e2020GL088050.doi：10.1029/2020GL0880507LI G,ZHANG J,HERRMANN H J,et al.Study ofaerodynamic grain entrainment in aeolian transportJ.GeophysicalResearchLetters,2020,47(11):e2019GL086574.doi：10.1029/2019GL0865748COMOLA F,KOK J F,CHAMECKI M,et al.Theintermittencyofwind-drivensandtransportJ.GeophysicalResearchLetters,2019,46(22):1343013440.doi：10.1029/2019GL0857399YINX,JIANGC,SHAOYP,etal.Large-eddy-simulationstudyonturbulentparticledepositionanditsdependenceonatmospheric-boundary-layer stabilityJ.Atmospheric Che-mistryandPhysics,2022,22(7):45094522.doi：10.5194/acp-22-4509-202210TANLH,ANZS,ZHANGK,etal.Intermittentaeoliansaltation over a Gobi surface:threshold,saltation layerheight,andhigh-frequencyvariabilityJ.JournalofGeophysicalResearch:EarthSurface,2020,125(1):e2019JF005329.doi：10.1029/2019JF00532911ZHENG X J,JIN T,WANG P.The influence of surfacestress fluctuation on saltation sand transport aroundthresholdJ.Journal of Geophysical Research:EarthSurface,2020,125(5):e2019JF005246.doi：10.1029/2019JF00524612KAWAMURA R.Study on sand movement by windJ.AnnualReportInstituteScienceandTechnology,1951,5(3):95112.13LETTAU K.Experimental and micro-meteorological fieldstudiesofdunemigrationM.LETTAUHH,LETTAUK,ExploringintheWorldsdriestclimate.Madison:UniversityofWisconsin-Madison,1978:110-147.14SRENSEN M.On the rate of aeolian sand transportJ.Geomorphology,2004,59(1-4):5362.doi：10.1016/j.geomorph.2003.09.00515SLINN W G N.Predictions for particle deposition tovegetative canopiesJ.Atmospheric Environment(1967),1982,16(7):17851794.doi：10.1016/0004-6981（82）90271-216ZHANGJ,SHAOY.Anewparameterizationofparticledrydeposition over rough surfacesJ.Atmospheric ChemistryandPhysics,2014,14(22):1242912440.doi：10.5194/acp-14-12429-201417PRESTONJH.Thedeterminationofturbulentskinfrictionby means of pitot tubesJ.The Journal of the RoyalAeronauticalSociety,1954,58(518):109121.doi：10.1017/s036839310009770418IRWINHPAH.Asimpleomnidirectionalsensorforwind-tunnelstudiesofpedestrian-levelwindsJ.JournalofWindEngineering and Industrial Aerodynamics,1981,7(3):219239.doi：10.1016/0167-6105（81）90051-919FAGEA,FALKNERVM.OntherelationbetweenheattransferandsurfacefrictionforlaminarflowR.Aeronautical Research Communications Report Memoran-dum,ReportNo.1408,1931.20LUDWIEGH.Instrumentformeasuringthewallshearingstress of turbulent boundary layersR.NACA TechnicalMemorandum1284,1950.21LIEPMANN H W,SKINNER G T.Shearing-stressmeasurements by use of a heated elementR.NACATechnicalNote3268,1954.22BELLHOUSE B J,SCHULTZ D L.The measurement offluctuating skin friction in air with heated thin-filmgaugesJ.JournalofFluidMechanics,1968,32(4):675680.doi：10.1017/s002211206800095923SUNBY,WANGPB,LUOJ,etal.Aflexiblehot-filmsensor array for underwater shear stress and transitionmeasurementJ.Sensors,2018,18(10):3469.doi：10.3390/s1810346924WALTERB,GROMKEC,LEONARDK,etal.Spatially25第 5 期黄宁等：风沙流环境地表剪切力测量99resolved skin friction velocity measurements using Irwinsensors:A calibration and accuracy analysisJ.Journal ofWindEngineeringandIndustrialAerodynamics,2012,104-106:314321.doi：10.1016/j.jweia.2012.02.018SHAO Y P,LU H.A simple expression for wind erosionthresholdfrictionvelocityJ.Journa