The traditional RANS turbulence parameterization scheme is modified on the basis of the Computational Fluid Dynamics model (Fluent), and the large eddy model (LES) turbulence scheme is selected, then the buoyancy model of temperature stratification is introduced, finally the modified Fluent model is obtained.The Weather Research and Forecasting model (WRF) and WRF coupled Fluent before and after modified are used to simulate the wind speed and direction of the 4 months from December 2010 to March 2011 in Jishan wind station in Poyang Lake region, Jiangxi and the evolution of the flow field and the temperature field near Jishan station in March 6, 2011.According to the comparison of the observation data, it is found that: After coupling Fluent, the better WRF simulation result of the wind speed and wind direction in Jishan station is obtained.Compared with coupling original Fluent (Example 1), WRF coupling modified Fluent (Example 2) can provide more accurate simulation result of the wind speed and wind direction at each height in Jishan station.The average error of WRF, Example 1 and Example 2 for the 4 months 70 m height average wind speed are 3.963, 2.727 and 2.224 m·s-1 respectively.The Matching ratio to the observed wind direction rose of Example 2 is higher than Example 1 at 70 m height mostly.Example 2 not only can simulate the thermal buoyancy effect under the daytime unstable state and the spatial inhomogeneous and higher turbulent kinetic energy, but also can simulate the temperature inversion of the atmospheric boundary layer at night, while Example 1 can not.All the above results indicate that after the modification of Fluent model, the simulation result of the method that WRF coupled Fluent to the wind field near Jishan Station is improved.To sum up, this research provides a new method and basis for the prediction of the atmospheric flow field in the bottom of atmospheric boundary layer and the multiscale coupling simulation, and has also laid the foundation for the development and improvement of Fluent model in the future .
Jiarong ZHANG
,
Xueling CHENG
. Numerical Simulation of Wind Field in Complex Terrain based on CFD Downscaling[J]. Plateau Meteorology, 2020
, 39(1)
: 172
-184
.
DOI: 10.7522/j.issn.1000-0534.2019.00005
[1]Bechmann A, Berg J, Courtney M S, al et, 2009.The Bolund experiment: Overview and background.[R].Ris?-R-1658(EN).
[2]Beljaars A C M, Walmsley J L, Taylor P A, 1987.A mixed spectral finite difference model for neutrally stratified boundary-layer flow over roughness change and topography[J].Boundary-Layer Meteorology, 38(3): 273-303
[3]Benchmann A, 2006.Large eddy simulation of atmospheric flow over complex terrain[D].Technical University of Denmark, 1-107.
[4]Chow F, Street R, 2009.Evaluation of turbulence closure models for large-eddy simulation over complex terrain: Flow over Askervein Hill[J].Journal of Applied Meteorology & Climatology, 48 (5): 1050-1065.
[5]Coppin P, Bradley E F, Finnigan J J, 1994.Measurements of flow over an elongated ridge and its thermal stability dependence: The mean field[J].Boundary-Layer Meteorology, 69(1/2): 173-199.
[6]Hunt J, Tampier F, Weng W S, al et, 1991.Air flow and turbulence over complex terrain: A colloquium and a computational workshop[J].Journal of Fluid Mechanics, 227(6): 667-688.
[7]Kim H G, Patel P C, Lee C M, 2000.Numerical simulation of wind flow over hilly terrain[J].Journal of Wind Engineering and Industrial Aerodynamics, 87: 45-60.
[8]Koblitz T, Bechmann A, Sogachev A, al et, 2011.Modification of CFD code to model the non-neutral atmospheric boundary layer[R].13th International Conference on Wind Engineering,Amsterdam,Netherlands.
[9]Kristof G, Racz N, Balogh M, 2007.Application of anasys-fluent on mesoscale atmospheric flow simulations[R].ANSYS Conference & 25th CADFEM User’s Meeting,November,Congress Center Dresden,Germany
[10]Launder B E, Spalding D B, 1972.Lectures in mathematical models of turbulence[M].London: Academic Press, 1-7169.
[11]Leblebici E, Ahmety G, Tuncer I H, 2011.Terrain modeling and atmospheric turbulent flow solutions based on meteorological weather forecast data[R].6th Ankara International Aerospace Conference,METU,Turkey Ankara
[12]Li L, Hu F, Cheng X L, al et, 2006.Numerical simulation of the flow within and over an intersection model with Reynolds-averaged Navier-Stokes method[J].Chinese Physics B, 15(1): 149-155.