Validation of Parallel WRF Downscaling Methodology using OpenFOAM

The main objective of this study is to obtain real-time atmospheric flow solutions using open source CFD solver OpenFOAM coupled with Numerical Weather Prediction (NWP) model; Weather Research Forecast (WRF). NWP can take moist convection, land surface parameterization, atmospheric boundary layer physics into account, but wind flow features finer than 1 km aren't captured by the turbulence physics of such models. CFD simulations, however, have proved to be useful at capturing the details of smaller scales due to a finer scale topography. Moreover, using the WRF weather prediction data as unsteady and spatially varying BCs for the CFD solution may prove to be one of the most realistic representations for the atmospheric flow field, and also allows daily power production estimations. Coupling the mesoscale weather prediction model WRF (Weather Research and Forecast) with the open source CFD solver OpenFOAM is done via using low resolution WRF data as unsteady and spatially varying boundary conditions for the OpenFOAM domain.For this purpose, a new unsteady and spatially varying boundary condition class (timeVaryingMixed) that switches between Neumann and Dirichlet depending on the flow is entering or exiting the domain to use the WRF data as boundary conditions without convergence issues for continuity, is developed.Due to real-time prediction requirement, parallelization of the process is of utmost importance. But the developed boundary condition class 'timeVaryingMixed' cannot be run in parallel using OpenFOAM's domain decomposition tool decomposePar as the indexes of cells change when the domainis decomposed. Parallelization of the process is done and made automatic using METIS to optimize the number of partitionboundaries, even when all the cells that arein neighbourhoodof the developed boundary condition timeVaryingMixed, are owned by 1 processor. Details about the methodology and parallelization of process will be given in the final paper.Unsteady OpenFOAM solutions coupled with WRF are performed using the methodology on high resolution stretching structured grids seen in Figure 2. High resolution (1.5 arcsec) ASTER GDEM topographical data is used to create the topography in order to capture the viscous effects which dominates the flow characteristics at the surface layer of the atmosphere where majority of the wind turbines reside. Simulations in Alaiz Mountain (Spain) are carried out and validation studies using the met-mast data from the region are done at the met-mast location at 5 different heights(118, 102, 90, 78, 40 meters) above the ground. As a preliminary result, time-series wind speed data at118and 40meters above ground is given in Figure 1.Results show a drastic improvement over the WRF results especially in thevicinity of the ground.
Citation Formats
E. LEBLEBİCİ and İ. H. TUNCER, “Validation of Parallel WRF Downscaling Methodology using OpenFOAM,” presented at the Wind Energy Science Conference (26 - 29 Haziran 2017), Technical University of Denmark, Lyngby, Danimarka, 2017, Accessed: 00, 2021. [Online]. Available: