为解决山区峡谷风场数值模拟过程中入口边界难以合理给定的问题,在中尺度气象模式(WRF)基础上利用多尺度耦合方法对山区峡谷桥址风场进行了精细化数值模拟。模拟过程中,首先基于WRF利用多尺度耦合的方法得到山区峡谷入口处的中尺度速度场信息,然后以入口边界位置风场波动情况为原则,对模拟的中尺度风场信息在耦合面进行分块,分块后分别运用多项式方法将其风速进行拟合,并通过UDF程序将拟合的速度赋给大涡模拟的入口边界。最后以张家界澧水大桥所在峡谷为研究背景,在桥址位置安装风速实时监测系统,将现场实测结果与所提方法的数值模拟结果进行了对比。研究结果表明:WRF的运行结果通过降尺后能较好地运用在山区峡谷风场CFD数值模拟的入口边界上,这种处理较好地解决了山区峡谷风场数值模拟过程中入口平均风的给定问题;分块多项式拟合插值方法解决了以往数值模拟过程中出现的“人为峭壁”问题;利用该方法可以较为准确地得到桥址处平均风速、风向角和风攻角等参数的分布情况。
In order to resolve the difficulty in providing reasonable inlet boundary conditions for the numerical simulation of wind field in mountainous area, the multiscale coupling method was adopted for accurately simulating the wind field on the bridge site at mountainous area based on a meteorological software of weather research and forecast (WRF). Firstly, at the process of simulation, the multiscale coupling method was adopted to obtain the mesoscale velocity field information by the WRF. Secondly, the coupling surface was divided into several pieces based on the principle of wind fluctuating degree in the inlet boundary. Block polynomial interpolation (BPI) fitting method was used to fit the wind velocity, and then the fitting velocities were given to the inlet boundary of the large eddy simulation (LES) by the UDF program. Finally, taking the I.ishui Bridge in Zhangjiajie as the research background, the wind velocities real-time observation stations were established at the bridge site for a comparison purpose. The measured results were compared with the numerical simulation results by the proposed method. The results show that the computational results of WRF can be well applied to the inlet boundary of wind field CFD numerical simulation in mountainous area after downsealing. The method can provide a reasonable value of average wind in inlet boundary for the numerical simulation in mountainous gorge area. Meanwhile, the BPI method provides a good solution to the irrational inlet problems caused by the "artificial cliff". The distributions of parameters including the average wind velocity, wind direction and wind attack angle can be accurately obtained by the proposed method.