In-situ Data Analysis System for High Resolution Meteorological Large Eddy Simulation Model

Author(s):  
Takuto Sato ◽  
Osamu Tatebe ◽  
Hiroyuki Kusaka
2020 ◽  
Author(s):  
Antti Hellsten ◽  
Klaus Ketelsen ◽  
Matthias Sühring ◽  
Mikko Auvinen ◽  
Björn Maronga ◽  
...  

Abstract. Large-eddy simulation provides a physically sound approach to study complex turbulent processes within the atmospheric boundary layer including urban boundary layer flows. However, such flow problems often involve a large separation of turbulent scales, requiring a large computational domain and very high grid resolution near the surface features, leading to prohibitive computational costs. To overcome this problem, an online LES-LES nesting scheme is implemented into the PALM model system 6.0. The hereby documented and evaluated nesting method is capable of supporting multiple child domains which can be nested within their parent domain either in a parallel or recursively cascading configuration. The nesting system is evaluated by simulating first a purely convective boundary layer flow system and then three different neutrally-stratified flow scenarios with increasing order of topographic complexity. The results of the nested runs are compared with corresponding non-nested high- and low-resolution results. The results reveal that the solution accuracy within the high-resolution nest domain is clearly improved as the solutions approach the non-nested high-resolution reference results. In obstacle-resolving LES, the two-way coupling becomes problematic as anterpolation introduces a regional discrepancy within the obstacle canopy of the parent domain. This is remedied by introducing canopy-restricted anterpolation where the operation is only performed above the obstacle canopy. The test simulations make evident that this approach is the most suitable coupling strategy for obstacle-resolving LES. The performed simulations testify that nesting can reduce the CPU time up to 80 % compared to the fine-resolution reference runs while the computational overhead from the nesting operations remained below 16 % for the two-way coupling approach and significantly less for the one-way alternative.


2021 ◽  
Vol 14 (6) ◽  
pp. 3185-3214
Author(s):  
Antti Hellsten ◽  
Klaus Ketelsen ◽  
Matthias Sühring ◽  
Mikko Auvinen ◽  
Björn Maronga ◽  
...  

Abstract. Large-eddy simulation (LES) provides a physically sound approach to study complex turbulent processes within the atmospheric boundary layer including urban boundary layer flows. However, such flow problems often involve a large separation of turbulent scales, requiring a large computational domain and very high grid resolution near the surface features, leading to prohibitive computational costs. To overcome this problem, an online LES–LES nesting scheme is implemented into the PALM model system 6.0. The hereby documented and evaluated nesting method is capable of supporting multiple child domains, which can be nested within their parent domain either in a parallel or recursively cascading configuration. The nesting system is evaluated by first simulating a purely convective boundary layer flow system and then three different neutrally stratified flow scenarios with increasing order of topographic complexity. The results of the nested runs are compared with corresponding non-nested high- and low-resolution results. The results reveal that the solution accuracy within the high-resolution nest domain is clearly improved as the solutions approach the non-nested high-resolution reference results. In obstacle-resolving LES, the two-way coupling becomes problematic as anterpolation introduces a regional discrepancy within the obstacle canopy of the parent domain. This is remedied by introducing canopy-restricted anterpolation where the operation is only performed above the obstacle canopy. The test simulations make evident that this approach is the most suitable coupling strategy for obstacle-resolving LES. The performed simulations testify that nesting can reduce the CPU time up to 80 % compared to the fine-resolution reference runs, while the computational overhead from the nesting operations remained below 16 % for the two-way coupling approach and significantly less for the one-way alternative.


Energies ◽  
2019 ◽  
Vol 12 (10) ◽  
pp. 1925 ◽  
Author(s):  
Takanori Uchida ◽  
Susumu Takakuwa

The first part of the present study investigated the relationship among the number of yaw gear and motor failures and turbulence intensity (TI) at all the wind turbines under investigation with the use of in situ data. The investigation revealed that wind turbine #7 (T7), which experienced a large number of failures, was affected by terrain-induced turbulence with TI that exceeded the TI presumed for the wind turbine design class to which T7 belongs. Subsequently, a computational fluid dynamics (CFD) simulation was performed to examine if the abovementioned observed wind flow characteristics could be successfully simulated. The CFD software package that was used in the present study was RIAM-COMPACT, which was developed by the first author of the present paper. RIAM-COMPACT is a nonlinear, unsteady wind prediction model that uses large-eddy simulation (LES) for the turbulence model. RIAM-COMPACT is capable of simulating flow collision, separation, and reattachment and also various unsteady turbulence–eddy phenomena that are caused by flow collision, separation, and reattachment. A close examination of computer animations of the streamwise (x) wind velocity revealed the following findings: As we predicted, wind flow that was separated from a micro-topographical feature (micro-scale terrain undulations) upstream of T7 generated large vortices. These vortices were shed downstream in a nearly periodic manner, which in turn generated terrain-induced turbulence, affecting T7 directly. Finally, the temporal change of the streamwise (x) wind velocity (a non-dimensional quantity) at the hub-height of T7 in the period from 600 to 800 in non-dimensional time was re-scaled in such a way that the average value of the streamwise (x) wind velocity for this period was 8.0 m/s, and the results of the analysis of the re-scaled data were discussed. With the re-scaled full-scale streamwise wind velocity (m/s) data (total number of data points: approximately 50,000; time interval: 0.3 s), the time-averaged streamwise (x) wind velocity and TI were evaluated using a common statistical processing procedure adopted for in situ data. Specifically, 10-min moving averaging (number of sample data points: 1932) was performed on the re-scaled data. Comparisons of the evaluated TI values to the TI values from the normal turbulence model in IEC61400-1 Ed.3 (2005) revealed the following: Although the evaluated TI values were not as large as those observed in situ, some of the evaluated TI values exceeded the values for turbulence class A, suggesting that the influence of terrain-induced turbulence on the wind turbine was well simulated.


2013 ◽  
Vol 20 (6) ◽  
pp. 1095-1112 ◽  
Author(s):  
A. Petronio ◽  
F. Roman ◽  
C. Nasello ◽  
V. Armenio

Abstract. In the present paper a state-of-the-art large eddy simulation model (LES-COAST), suited for the analysis of water circulation and mixing in closed or semi-closed areas, is presented and applied to the study of the hydrodynamic characteristics of the Muggia bay, the industrial harbor of the city of Trieste, Italy. The model solves the non-hydrostatic, unsteady Navier–Stokes equations, under the Boussinesq approximation for temperature and salinity buoyancy effects, using a novel, two-eddy viscosity Smagorinsky model for the closure of the subgrid-scale momentum fluxes. The model employs: a simple and effective technique to take into account wind-stress inhomogeneity related to the blocking effect of emerged structures, which, in turn, can drive local-scale, short-term pollutant dispersion; a new nesting procedure to reconstruct instantaneous, turbulent velocity components, temperature and salinity at the open boundaries of the domain using data coming from large-scale circulation models (LCM). Validation tests have shown that the model reproduces field measurement satisfactorily. The analysis of water circulation and mixing in the Muggia bay has been carried out under three typical breeze conditions. Water circulation has been shown to behave as in typical semi-closed basins, with an upper layer moving along the wind direction (apart from the anti-cyclonic veering associated with the Coriolis force) and a bottom layer, thicker and slower than the upper one, moving along the opposite direction. The study has shown that water vertical mixing in the bay is inhibited by a large level of stable stratification, mainly associated with vertical variation in salinity and, to a minor extent, with temperature variation along the water column. More intense mixing, quantified by sub-critical values of the gradient Richardson number, is present in near-coastal regions where upwelling/downwelling phenomena occur. The analysis of instantaneous fields has detected the presence of large cross-sectional eddies spanning the whole water column and contributing to vertical mixing, associated with the presence of sub-surface horizontal turbulent structures. Analysis of water renewal within the bay shows that, under the typical breeze regimes considered in the study, the residence time of water in the bay is of the order of a few days. Finally, vertical eddy viscosity has been calculated and shown to vary by a couple of orders of magnitude along the water column, with larger values near the bottom surface where density stratification is smaller.


2015 ◽  
Vol 8 (8) ◽  
pp. 2515-2551 ◽  
Author(s):  
B. Maronga ◽  
M. Gryschka ◽  
R. Heinze ◽  
F. Hoffmann ◽  
F. Kanani-Sühring ◽  
...  

Abstract. In this paper we present the current version of the Parallelized Large-Eddy Simulation Model (PALM) whose core has been developed at the Institute of Meteorology and Climatology at Leibniz Universität Hannover (Germany). PALM is a Fortran 95-based code with some Fortran 2003 extensions and has been applied for the simulation of a variety of atmospheric and oceanic boundary layers for more than 15 years. PALM is optimized for use on massively parallel computer architectures and was recently ported to general-purpose graphics processing units. In the present paper we give a detailed description of the current version of the model and its features, such as an embedded Lagrangian cloud model and the possibility to use Cartesian topography. Moreover, we discuss recent model developments and future perspectives for LES applications.


2019 ◽  
Vol 12 (6) ◽  
pp. 2523-2538 ◽  
Author(s):  
Sadiq Huq ◽  
Frederik De Roo ◽  
Siegfried Raasch ◽  
Matthias Mauder

Abstract. Large-eddy simulation (LES) has become a well-established tool in the atmospheric boundary layer research community to study turbulence. It allows three-dimensional realizations of the turbulent fields, which large-scale models and most experimental studies cannot yield. To resolve the largest eddies in the mixed layer, a moderate grid resolution in the range of 10 to 100 m is often sufficient, and these simulations can be run on a computing cluster with a few hundred processors or even on a workstation for simple configurations. The desired resolution is usually limited by the computational resources. However, to compare with tower measurements of turbulence and exchange fluxes in the surface layer, a much higher resolution is required. In spite of the growth in computational power, a high-resolution LES of the surface layer is often not feasible: to fully resolve the energy-containing eddies near the surface, a grid spacing of O(1 m) is required. One way to tackle this problem is to employ a vertical grid nesting technique, in which the surface is simulated at the necessary fine grid resolution, and it is coupled with a standard, coarse, LES that resolves the turbulence in the whole boundary layer. We modified the LES model PALM (Parallelized Large-eddy simulation Model) and implemented a two-way nesting technique, with coupling in both directions between the coarse and the fine grid. The coupling algorithm has to ensure correct boundary conditions for the fine grid. Our nesting algorithm is realized by modifying the standard third-order Runge–Kutta time stepping to allow communication of data between the two grids. The two grids are concurrently advanced in time while ensuring that the sum of resolved and sub-grid-scale kinetic energy is conserved. We design a validation test and show that the temporally averaged profiles from the fine grid agree well compared to the reference simulation with high resolution in the entire domain. The overall performance and scalability of the nesting algorithm is found to be satisfactory. Our nesting results in more than 80 % savings in computational power for 5 times higher resolution in each direction in the surface layer.


Sign in / Sign up

Export Citation Format

Share Document