scholarly journals Capturing Day-to-day Diurnal Variations in Stability in the Convective Atmospheric Boundary Layer Using Large Eddy Simulation

2018 ◽  
Vol 12 (1) ◽  
pp. 107-131 ◽  
Author(s):  
Jordan Nielson ◽  
Kiran Bhaganagar

Introduction:Large Eddy Simulation (LES) modelers must begin to answer the question of how to better incorporate large datasets into simulations. This question is important because, at a given location, the diurnal, seasonal, and day-to-day variations of atmospheric stability have significant consequences for the power generated by wind turbines. The following study provides a methodology to obtain discrete values of surface flux, inversion height and geostrophic wind for LES using field data over multiple diurnal cycles (averaged over a month) at 12 Local Time (LT) (during the convective ABL). The methodology will allow the discrete LES to quantify the day-to-day variations over multiple diurnal cycles.Methods:The study tests the hypothesis that LES can capture the mean velocity and TKE profiles from the averaged variations in surface heat flux at 12 LT measured in the field (mean, +1 standard deviation, and -1 standard deviation). The discrete LES from the mean, +1 standard deviation, and -1 standard deviation surface heat flux represent the variations in the ABL due to the day-to-day variations in surface heat flux. The method calculates the surface heat flux for the NREL NWTC M5 dataset. The field data were used to generate Probability Density Functions (PDFs) of surface heat flux for the January and July 12 LT. The PDFs are used to select the surface heat fluxes as inputs into the discrete LES.Results / Conclusion:A correlation function between the surface heat flux and the boundary layer height was determined to select the initial inversion height, and the geostrophic departure function was used to determine the geostrophic wind for each surface heat flux. The LES profiles matched the averaged velocity profiles from the field data to 4% and the averaged TKE profiles to 6% and, therefore, validated the methodology. The method allows for further quantification of day-to-day stability variations using LES.

2019 ◽  
Vol 76 (3) ◽  
pp. 785-799 ◽  
Author(s):  
Cheng Liu ◽  
Evgeni Fedorovich ◽  
Jianping Huang ◽  
Xiao-Ming Hu ◽  
Yongwei Wang ◽  
...  

AbstractEntrainment is critical to the development of the atmospheric convective boundary layer (CBL), but little is known about how entrainment is impacted by the aerosol radiative effect. An aerosol radiation transfer model is used in conjunction with large-eddy simulation (LES) to quantify the impact of aerosol shortwave radiative heating on entrainment and thermodynamics of an idealized dry CBL under aerosol-loading conditions. An entrainment equation is derived within the framework of a zero-order model (ZOM) with the aerosol radiative heating effect included; the equation is then examined against the LES outputs for varying aerosol optical depths (AODs) and free-atmosphere stratification scenarios. The results show that the heat flux profiles become more nonlinear in shape as compared to the case of the clean (no aerosol pollution) CBL, with the degree of nonlinearity being highly dependent on the AOD of the layer for the given type of radiation-absorbing aerosols. As AOD increases, less solar radiation reaches the surface and thus the surface heat flux becomes smaller, and both actual (LES) and ZOM-derived entrainment flux ratios decrease. This trend is opposite to the clean CBL where the LES-predicted flux ratios show an increasing trend with diminishing surface heat flux, while the ZOM-calculated flux ratio remains constant. The modified dimensionless entrainment rate closely follows the −1 power law with a modified Richardson number. The study suggests that including the aerosol radiative effect may improve numerical air quality predictions for heavy-air-pollution events.


1990 ◽  
Vol 27 (8) ◽  
pp. 1048-1060 ◽  
Author(s):  
David W. S. Eaton ◽  
Frederick A. Cook

The Valhalla complex, situated in the Omineca crystalline belt in southeastern British Columbia, is a Cordilleran metamorphic core complex bordering the suture zone between Quesnellia and North American rocks. The region is tectonically interposed between a convergent plate margin along Canada's west coast and the stable North American craton, and is characterized by a crustal thickness of ~ 35 km, high surface heat flux, and elevated lower crustal electrical conductivity. In this study, Lithoprobe deep-crustal seismic-reflection data, potential-field data, and geological constraints have been used to gain a better understanding of crustal structure in the vicinity of the Valhalla complex. Analysis of Bouguer gravity and total-field aeromagnetic data indicates that mafic oceanic rocks and various syn- and post-accretionary granitoid plutonic rocks are not major constituents of the upper crust underlying the complex. The seismic data reveal a moderately reflective upper crust and image several fault zones, including a very high amplitude, west-dipping reflection that is interpreted as a significant Late Cretaceous or Paleocene thrust fault. The fault-zone reflectivity may be related to compositional heterogeneity and (or) seismic anisotropy associated with mylonites. The lower crust appears to be nonreflective, in contrast with other areas of high surface heat flux and elevated lower crustal conductivity. Taken together, the various data show that the Valhalla complex is likely cored by North American metasedimentary rocks and reveal features related to the Jurassic to Paleocene compressional fabric, which has been largely overprinted at the surface by subsequent Eocene extension.


2010 ◽  
Vol 67 (5) ◽  
pp. 1655-1666 ◽  
Author(s):  
David M. Romps ◽  
Zhiming Kuang

Abstract Tracers are used in a large-eddy simulation of shallow convection to show that stochastic entrainment (and not cloud-base properties) determines the fate of convecting parcels. The tracers are used to diagnose the correlations between a parcel’s state above the cloud base and both the parcel’s state at the cloud base and its entrainment history. The correlation with the cloud-base state goes to zero a few hundred meters above the cloud base. On the other hand, correlations between a parcel’s state and its net entrainment are large. Evidence is found that the entrainment events may be described as a stochastic Poisson process. A parcel model is constructed with stochastic entrainment that is able to replicate the mean and standard deviation of cloud properties. Turning off cloud-base variability has little effect on the results, which suggests that stochastic mass-flux models may be initialized with a single set of properties. The success of the stochastic parcel model suggests that it holds promise as the framework for a convective parameterization.


2019 ◽  
Vol 43 (6) ◽  
pp. 625-638 ◽  
Author(s):  
Jordan Nielson ◽  
Kiran Bhaganagar

A novel and a robust high-fidelity numerical methodology has been developed to realistically estimate the net energy production of full-scale horizontal axis wind turbines in a convective atmospheric boundary layer, for both isolated and multiple wind turbine arrays by accounting for the wake effects between them. Large eddy simulation has been used to understand the role of atmospheric stability in net energy production (annual energy production) of full-scale horizontal axis wind turbines placed in the convective atmospheric boundary layer. The simulations are performed during the convective conditions corresponding to the National Renewable Energy Laboratory field campaign of July 2015. A mathematical framework was developed to incorporate the field-based measurements as boundary conditions for the large eddy simulation by averaging the surface flux over multiple diurnal cycles. The objective of the study is to quantify the role of surface flux in the calculation of energy production for an isolated, two and three wind turbine configuration. The study compares the mean value, +1 standard deviation, and −1 standard deviation from the measured surface flux to demonstrate the role of surface heat flux. The uniqueness of the study is that power deficits from large eddy simulation were used to determine wake losses and obtain a net energy production that accounts for the wake losses. The frequency of stability events, from field measurements, is input into the calculation of an ensemble energy production prediction with wake losses for different wind turbine arrays. The increased surface heat flux increases the atmospheric turbulence into the wind turbines. Higher turbulence results in faster wake recovery by a factor of two. The faster wake recovery rates result in lowering the power deficits from 46% to 28% for the two-turbine array. The difference in net energy production between the +1 and −1 standard deviation (with respect to surface heat flux) simulations was 10% for the two-turbine array and 8% for the three-turbine array. An ensemble net energy production by accounting for the wake losses indicated the overestimation of annual energy production from current practices could be corrected by accounting for variation of surface flux from the mean value.


2020 ◽  
Vol 143 (1) ◽  
Author(s):  
Yu Xia ◽  
Patrick Sharkey ◽  
Stefano Orsino ◽  
Mike Kuron ◽  
Florian Menter ◽  
...  

Abstract Accurate numerical prediction of surface heat transfer in the presence of film cooling within aero-engine sub-components, such as blade effusion holes and combustor liners, has long been a goal of the aero-engine industry. It requires accurate simulation of the turbulent mixing and reaction processes between freestream and the cooling flow. In this study, the stress-blended eddy simulation (SBES) turbulence model is used together with the flamelet generated manifold (FGM) combustion model to calculate the surface heat flux upstream and downstream of an effusion cooling hole. The SBES model employs a blending function to automatically switch between Reynolds-averaged Navier–Stokes (RANS) and large eddy simulation (LES) based on the local flow features, and thus significantly reduces the computational cost compared to a full LES simulation. All simulations are run using ansys fluent®, a commercial finite-volume computational fluid dynamics (CFD) solver. The test case corresponds to an experimental rig run at Massachusetts Institute of Technology (MIT), which is essentially a flat plate brushed by a uniform freestream of argon with ethylene seeded inside, and is cooled by either a reacting air or a non-reacting nitrogen jet inclined at 35 deg to the freestream. Calculations are performed for both reacting and non-reacting jet cooling cases across a range of jet-to-stream blowing ratios and compared with the experimental data. The effects of mesh resolution are also investigated. Calculations are also performed across a range of Damköhler number (i.e., flow to chemical time ratio) from zero to 30, with unity blowing ratio, and the differences in the maximum surface heat flux magnitude in the reacting and non-reacting cases at a specific location downstream of the hole are investigated. Results from these analyses show good correlation with the experimental heat flux data upstream and downstream of the cooling hole, including the heat flux augmentation due to local reaction. Results from the Damköhler number sweep also show a good match with the experimental data across the range investigated.


2009 ◽  
Vol 66 (2) ◽  
pp. 412-431 ◽  
Author(s):  
Rob Stoll ◽  
Fernando Porté-Agel

Abstract Large-eddy simulation, with recently developed dynamic subgrid-scale models, is used to study the effect of heterogeneous surface temperature distributions on regional-scale turbulent fluxes in the stable boundary layer (SBL). Simulations are performed of a continuously turbulent SBL with surface heterogeneity added in the form of streamwise transitions in surface temperature. Temperature differences between patches of 6 and 3 K are explored with patch length scales ranging from one-half to twice the equivalent homogeneous boundary layer height. The surface temperature heterogeneity has important effects on the mean wind speed and potential temperature profiles as well as on the surface heat flux distribution. Increasing the difference between the patch temperatures results in decreased magnitude of the average surface heat flux, with a corresponding increase in the mean potential temperature in the boundary layer. The simulation results are also used to test existing models for average surface fluxes over heterogeneous terrain. The tested models fail to fully represent the average turbulent heat flux, with models that break the domain into homogeneous subareas grossly underestimating the heat flux magnitude over patches with relatively colder surface temperatures. Motivated by these results, a new parameterization based on local similarity theory is proposed. The new formulation is found to correct the bias over the cold patches, resulting in improved average surface heat flux calculations.


1977 ◽  
Vol 80 (1) ◽  
pp. 153-177 ◽  
Author(s):  
R. A. Antonia ◽  
H. Q. Danh ◽  
A. Prabhu

Measurements of both the velocity and the temperature field have been made in the thermal layer that grows inside a turbulent boundary layer which is subjected to a small step change in surface heat flux. Upstream of the step, the wall heat flux is zero and the velocity boundary layer is nearly self-preserving. The thermal-layer measurements are discussed in the context of a self-preserving analysis for the temperature disturbance which grows underneath a thick external turbulent boundary layer. A logarithmic mean temperature profile is established downstream of the step but the budget for the mean-square temperature fluctuations shows that, in the inner region of the thermal layer, the production and dissipation of temperature fluctuations are not quite equal at the furthest downstream measurement station. The measurements for both the mean and the fluctuating temperature field indicate that the relaxation distance for the thermal layer is quite large, of the order of 1000θ0, where θ0is the momentum thickness of the boundary layer at the step. Statistics of the thermal-layer interface and conditionally sampled measurements with respect to this interface are presented. Measurements of the temperature intermittency factor indicate that the interface is normally distributed with respect to its mean position. Near the step, the passive heat contaminant acts as an effective marker of the organized turbulence structure that has been observed in the wall region of a boundary layer. Accordingly, conditional averages of Reynolds stresses and heat fluxes measured in the heated part of the flow are considerably larger than the conventional averages when the temperature intermittency factor is small.


2010 ◽  
Vol 40 (5) ◽  
pp. 934-948 ◽  
Author(s):  
Steven J. Lentz

Abstract The mean heat and salt balances over the Middle Atlantic Bight continental shelf are investigated by testing the hypothesis that surface fluxes of heat or freshwater are balanced by along-isobath fluxes resulting from the mean, depth-averaged, along-isobath flow acting on the mean, depth-averaged, along-isobath temperature or salinity gradient. This hypothesized balance is equivalent in a Lagrangian frame to a column of water, for example, warming because of surface heating as it is advected southward along isobath by the mean flow. Mean depth-averaged temperatures increase from north to south along isobath at a rate of 2°C (1000 km)−1 at midshelf, which is consistent with the hypothesized balance and mean surface heat flux estimates from the 50-yr NCEP Reanalysis. However, mean surface heat flux estimates from the higher-resolution 20-yr Objectively Analyzed Air–Sea Fluxes (OAFlux) reanalysis are too small to balance the along-isobath heat flux divergence implying a cross-shelf heat flux convergence. It is unclear which surface heat flux estimate, NCEP or OAFlux, is more accurate. The cross-shelf heat flux convergence resulting from the mean cross-shelf circulation is too small to balance the along-isobath heat flux divergence. Mean depth-averaged salinities increase from north to south along isobath at a rate of 1 (psu) (1000 km)−1 at midshelf. Mean precipitation and evaporation rates nearly balance so that the net freshwater flux is too small by more than an order of magnitude to account for the observed along-isobath increase in salinity. The cross-shelf salt flux divergence resulting from the mean cross-shelf circulation has the wrong sign to balance the divergence in the along-isobath salt flux. These results imply there must be an onshore “eddy” salt flux resulting from the time-dependent current and salinity variability. The along-isobath temperature and salinity gradients compensate for each other so that the mean, depth-averaged, along-isobath density gradient is approximately zero. This suggests that there may be a feedback between the along-isobath density gradient and the onshore salt and heat fluxes that maintains the density gradient near zero.


Sign in / Sign up

Export Citation Format

Share Document