Search NASA⌕ Search

SEARCH · Search NASA

Results for “shear layer turbulence”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Baroclinic Residual Circulation and Mass Transport Due to Internal Tides

Baroclinic (BC) tidal residual circulation due to internal tides is investigated around islands over a shallow ridge using a numerical ocean model. Internal tides enhance vertical mixing over shallow slopes, leading to horizontal density gradients that drive BC residual circulation along the main thermocline. For a strongly stratified summer case, the vertical diffusivity estimated by the Mellor and Yamada turbulence closure model exceeds 1 × 10 -2 m 2 s -1 , and the velocity of BC residual circulations reaches 0.2 m s -1 . The magnitude of BC residual circulation is larger than that of barotropic residual tidal circulation, implying that BC residual circulation due to internal tides plays an important role in coastal ocean circulation. Furthermore, BC residual flow accounts for an equal percentage (5%) of the total tidal kinetic energy as the barotropic residual flow under summer stratification conditions. Results from a coupled sediment resuspension and transport model show the growth of intermediate nepheloid layers formed by strong bottom shear stress and BC residual flow. The residual component contributes largely (23% of the total in summer) to the total suspended sediment flux. Seasonal variability is explored, with weaker winter stratification leading to reduced mixing, and thus weaker BC residual circulation and sediment flux. Lastly, the magnitude of the BC residual circulation is also shown to be proportional to the square of the tidal amplitude.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗
DOE OSTI2022

Development of Explainable Data-Driven Turbulence Models with Application to Liquid Fuel Nuclear Reactors

Liquid fuel nuclear reactors offer innovative possibilities in terms of nuclear reactor designs and passive safety systems. Molten Salts Reactors (MSRs) with a fast spectrum are a particular type of these reactors using liquid fuel. MSFRs often involve large open cavities in their core in which the liquid fuel circulates at a high speed to transport the heat generated by the nuclear reactions into the heat exchangers. This high-speed flow yields a turbulent field with large Reynolds numbers in the reactor core. Since the nuclear power, the neutron precursor’s transport and the thermal exchanges are strongly coupled in the MSFR’s core cavity, having accurate turbulent models for the liquid fuel flow is necessary to avoid introducing significant errors in the numerical simulations of these reactors. Nonetheless, high-accuracy simulations of the turbulent flow field in the reactor cavity of these reactors are usually prohibitively expensive in terms of computational resources, especially when performing multiphysics numerical calculations. Therefore, in this work, we propose a novel method using a modified genetic algorithm to optimize the calculation of the Reynolds Shear Stress Tensor (RST) used for turbulence modeling. The proposed optimization methodology is particularly suitable for advanced liquid fuel reactors such as the MSFRs since it allows the development of high-accuracy but still low-computational-cost turbulence models for the liquid fuel. We demonstrate the applicability of this approach by developing high accuracy Reynolds-Averaged Navier–Stokes (RANS) models (averaged flow error less than 5%) for a low and a large aspect ratio in a Backward-Facing Step (BFS) section particularly challenging for RANS models. The newly developed turbulence models better capture the flow field after the boundary layer tipping, over the extent of the recirculation bubble, and near the boundary layer reattachment region in both BFS configurations. The main reason for these improvements is that the developed models better capture the flow field turbulent anisotropy in the bulk region of the BFS. Then, we illustrate the interest in using this turbulence modeling approach for the case of an MSFR by quantifying the impact of the turbulence modeling on the reactor key parameters.

42 ENGINEERING↗
DOE OSTI2025

The effect of shear flow on the resistive tearing instability

We develop a new scaling theory for the resistive tearing mode instability of a current sheet with a strong shear flow across the layer. The growth rate decreases with increasing flow shear and is completely stabilized as the shear flow becomes Alfvénic: both in the constant-Ψ regime, as in previous results, but we also show that the growth rate is in fact suppressed more strongly in the nonconstant-Ψ regime. As a consequence, for sufficiently large flow shear, the maximum of the growth rate is always affected by the shear suppression, and the wavenumber at which this maximum growth rate is attained is an increasing function of the strength of the flow shear. These results may be important for the onset of reconnection in imbalanced MHD turbulence.

astrophysical plasmas↗
DOE OSTI2022

Rational surfaces, flows and radial structure in the TJ-II stellarator

Abstract In this work, we report on the results obtained by measuring several turbulent quantities well inside the plasma edge by means of a Langmuir probe during dynamical rotational transform scans in the TJ-II stellarator, while applying a radial electric field to the edge plasma using a biasing probe. By calculating the intermittence parameter from floating potential measurements, we are able to identify a major low order rational surface and hence relate the probe measurements to the local value of the rotational transform. Based on the former, we are able to show that the poloidal plasma velocity (and hence radial electric field) has a significant radial structure that is clearly related to the rotational transform profile and in particular the lowest order rational surfaces in the range studied. The poloidal velocity is also affected by the edge biasing. The particle flux Γ was also found to exhibit a radial pattern, as did the flow shear suppression term ω E × B , but the relation of the former to the low-order rational surfaces was less clear. We surmise that this lack of direct correspondence is due to an unknown term in the turbulence evolution equation: the instability growth rate, γ . We make use of a reduced Magnetohydrodynamic turbulence model to interpret the results. Overall, a picture is obtained in which the plasma self-organizes towards a state with a clear radial pattern of the radial electric field, in line with expectations from some numerical studies describing the spontaneous formation of an ‘ E × B staircase’, consisting of alternating layers with fast and slow radial transport. In this state, the radial profiles of various quantities (density, temperature, pressure) will not be smooth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Dynamics of multilayer Rayleigh–Taylor instability at moderately high Atwood numbers

This paper investigates the multilayer Rayleigh–Taylor instability (RTI) using statistically stationary experiments conducted in a gas tunnel. Employing diagnostics such as particle image velocimetry (PIV) and planar laser induced fluorescence (PLIF), we make simultaneous velocity–density measurements to study how dynamics and mixing are linked in this variable density flow. Experiments are conducted in a newly built, blow-down three-layer gas tunnel facility. Mixing between three gas streams is studied, where the top and bottom streams are comprised of air, and the middle stream is an air–helium mixture. Shear is minimized between these streams by matching their inlet velocities. The four experimental conditions investigated here consist of two different density ratios (Atwood numbers 0.3 and 0.6), each investigated at two instability development times (or equivalently, two streamwise locations), and all experiments are with the same middle stream thickness of 3 cm. The growth of the middle layer is measured using laser-based planar Mie scattering visualization. The mixing width is found to grow linearly with time at late times. Various quantitative measures of molecular mixing indicate a very high degree of molecular mixing at late times in the multilayer RTI flow. The vertical turbulent mass flux $a_y$ is calculated. In addition to mostly negative values of $a_y$ , typical of buoyancy-dominated flows due to negative correlation between velocity and density fluctuations, positive regions are also observed in profiles of $a_y$ due to entrainment and erosion at the lower edge of the mixing region. Global energy budgets are calculated for the multilayer RTI flow at late times and it is found that the majority of potential energy released has been dissipated due to viscous effects, and a large value of mixing efficiency ( $\sim$ 60 %) is observed.

Mechanics↗
DOE OSTI2025

Observations of Subduction, Downward Heat Flux and Dense Filament Collapse in the Northern Gulf of Mexico

Submesoscale processes are important contributors to the global heat budget and generally support upward heat transport through restratification. However, in salinity‐stratified regions, such as the northern Gulf of Mexico with its influx of freshwater from the Mississippi‐Atchafalaya river system, temperature can act like a passive tracer and submesoscale processes can contribute to downward heat transport. Oceanic heat content is a factor in many environmental risks the region faces, for example, hurricane intensification, and marine heatwaves. During the 2022 field campaign of the Submesoscales Under Near‐Resonant Inertial Shear Experiment, a sampling plan was developed to study such submesoscale processes in high resolution. Over 31 hr, four assets (two research ships and two remotely controlled boats) drove in parallel across a dense filament, capturing its evolution in time and space. The observations show that surface waters, warmed by daytime solar radiation, were subducted and that the associated overturning circulation transported heat below the surface layer where it was later irreversibly mixed away. The estimated downward heat flux was as strong as the concurrent net air‐sea heat flux into the ocean. The filament was then observed to rapidly collapse which we attribute to boundary layer turbulence and the breakdown of geostrophic balance. The collapsing fronts display behaviors indicative of gravity currents. These observations highlight how in salinity‐stratified regions, frontal dynamics can be associated with downward heat flux and how the submesoscale can play an important role in the oceanic heat budget.

58 GEOSCIENCES↗
DOE OSTI2022

SOL width broadening by spreading of pedestal turbulence

The pedestal turbulence intensity required to convert the thin, laminar H-mode scrape-off layer (SOL) to a broad turbulent SOL is calculated using the theory of turbulence spreading. A lower bound on the pedestal turbulence level to exceed the neoclassical heuristic drift (HD) width is derived. A reduced model of SOL turbulence spreading is used to determine the SOL width as a function of intensity flux from the pedestal to the SOL. The cross-over value for exceeding the HD model width is then calculated. We determine the pedestal turbulence levels—and the critical scalings thereof—required to achieve this level of broadening. Both drift wave and ballooning mode turbulence are considered. Here, a sensitivity analysis reveals that the key competition is that between spreading and linear E × B shear damping. The required pedestal turbulence levels scale with ρ/R.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simultaneous measurements of kinetic and scalar energy spectrum time evolution in the Richtmyer–Meshkov instability upon reshock

The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths , vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn. , vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$ , and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids , vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. ( Phys. Fluids , vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Unstable turbulent channel flow response to spanwise-heterogeneous heat fluxes: Prandtl's secondary flow of the third kind

Turbulent secondary flows are defined as Prandtl's secondary flow of the first or second kind, the former produced by stretching and/or tilting of vorticity, the latter produced via spatial heterogeneity of Reynolds stresses. Both mechanisms are instantaneously active within inertia-dominated wall turbulence; Reynolds stress spatial heterogeneity is required for Reynolds-averaged secondary flows. Spanwise-variable surface roughness can induce turbulent stress spatial heterogeneity in the spanwise–wall-normal plane and provide sustenance for streamwise-aligned mean secondary flows. Herein, we demonstrate that turbulent secondary flows can also be sustained by spanwise variability in the surface heat flux in unstably stratified turbulent channels, defined hereafter as Prandtl's secondary flow of the third kind. Support for this mechanism is established with scaling arguments, while large-eddy simulation is used to model inertia-dominated channel turbulence responding to a lower boundary with uniform aerodynamic/hydrodynamic roughness but spanwise-variable surface heat flux. Transport equations for streamwise vorticity and turbulent kinetic energy, k, outline the conditions needed for third-kind production: shear and buoyancy production over the elevated heat flux regions necessitates lateral entrainment of low- k fluid, inducing mean counter-rotating secondary cells aligned such that upwelling and downwelling occur over the high and low heat flux regions, respectively. Buoyancy-driven production of k alters aggregate flow response and thus is a distinctly different mechanism responsible for sustenance of secondary flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Scalar Power Spectra and Scalar Structure Function Evolution in the Richtmyer–Meshkov Instability Upon Reshock

The Richtmyer–Meshkov instability of a twice-shocked gas interface is studied using high-speed planar laser-induced fluorescence in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. Furthermore, this initial condition is accelerated by a shock of nominal strength M = 1.9, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. Three individual experiments are analyzed, the energy spectrum and the structure functions of the light gas mole fraction field are calculated and compared.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Insights into the relation between vertical cloud structure and dynamics of three polar lows: Observations from COMBLE

The vertical structure of three polar lows is described using profiling radar, lidar, and passive remote sensors deployed at a coastal site in northern Norway. The data were collected as part of the Cold-air Outbreaks in the Marine Boundary Layer Experiment (COMBLE). These data are examined in the context of satellite imagery, operational weather radar reflectivity imagery, and output from the AROME-Arctic model. A comparison to satellite images shows that AROME-Arctic captured the polar lows well. All three polar lows form in marine cold-air outbreaks and are surrounded by open-cellular convection typical of such outbreaks. Two of the lows form in a forward shear environment, and one under reverse shear. Initially, the three polar lows are comma-shaped; two of them transition to be more spirali-form at their mature stage and have mostly cloud-free warm cores. The warm cores are the result of a warm-seclusion process. All three lows have stratiform precipitation bands, marked by little cloud liquid water, and rather high surface precipitation rate. The vertical drafts and turbulence in these stratiform clouds are generally weak. In conclusion, all three lows also have convective clouds, which have stronger vertical drafts, stronger turbulence, and pockets of high liquid water content.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Active Aerodynamic Load Control for Wind Turbines

The goal of this project was to develop and demonstrate an advanced dielectric barrier discharge (DBD) plasma actuator technology. We set out to demonstrate the efficacy and impact of the new actuator technology as a key component of an active load control system for wind turbines. A DBD plasma actuator consists of a thin layer of dielectric material separating a pair of offset electrodes. One of those electrodes is embedded between the dielectric layer and a non-conductive substrate, while the other is exposed to air. When driven by an appropriate high voltage waveform, the device ionizes the air near the surface of the dielectric and adjacent to the exposed electrode. Collisions between the ions in the plasma and neutral air molecules result in a wall-jet – a region of induced air velocity that can be used to modify the flow around a lifting body. Use of the device near the trailing edge of a wind turbine blade, designed in such a way as to amplify the effects of the flow-modifying device, can result in large changes to the global forces experienced by the blade. Because of the fast response time of the device, it can allow the turbine to react in real time to changes in the wind associated with turbulence, wind shear, gusts, and the like, when paired with appropriate sensors and control algorithms. The goal of this project was to increase the capacity of the device to induce velocity on its surface at the levels required by large, utility-scale wind turbines. The first technical goal was to increase the induced velocity from the current industry-best of about 4 m/s to 10 m/s by modifying the electrical waveform used to drive the device and by introducing a semi-conductive surface coating to control electrostatic charge build-up on the surface. The second goal was to use the device to modify the lift on a representative airfoil in a wind tunnel, demonstrating a reduction in lift coefficient of 0.2 or better. The third goal was to produce an actuator-induced velocity of 20 m/s. Alongside our partners at the University of Texas at Dallas, we also applied modern, advanced design methods to optimize the impact of the technology on the design of wind turbines. We also had planned to install a segmented, active load control system on a test turbine to demonstrate the ability to reduce unsteady aerodynamic forces on the turbine associated with changes in the wind. The highest induced velocity achieved during this research was 11 m/s. However, practical design constraints limited the change in lift coefficient to about 0.12 for a representative airfoil in the wind tunnel at a Reynolds number of 400,000. The primary conclusion is that the plasma actuator control authority remains insufficient for practical purposes when extrapolated to Reynolds numbers over 1 million. The potential impact of the active lift control concept was evaluated through detailed simulations for three wind turbine sizes: a 3.4MW onshore turbine, a 10MW offshore turbine, and a 15MW offshore turbine. A feedback control system was designed for each turbine within two scenarios: one where the active lift control is used as a retrofit capability on the baseline design, and the other where the designers were permitted to “upscale” the turbines in order increase annual energy production. The levelized cost of energy was then evaluated for the range of turbine sizes and design configuration. It was found that the LCOE reduction associated with active lift control fell in the range of 0.7% to 7.2%, with the highest reduction associated with upscaling the 3.4MW turbine.

17 WIND ENERGY↗
DOE OSTI2022

Correlation between the relative blob fraction and plasma parameters in NSTX

Here, we present a simple analysis of the fraction of blobs observed in the scrape-off layer of NSTX and calculates the correlation between this blob fraction and various plasma parameters. The measurements were made using the gas puff imaging (GPI) diagnostic on a database of 103 shots from 2010 with neutral beam power P NBI = 0–6 MW. The blob fraction ranged from f b = 0.1% to 4.8%, where f b is defined as the fraction of time the GPI signal was larger than three times its average value within 0–7.5 cm outside the separatrix. The blob fraction generally decreased with increasing neutral beam injection power and was lower on average in H-mode plasmas than Ohmic and L-mode plasmas. The blob fraction had its highest correlation with the measured poloidal turbulence velocity and its radial gradient just inside the separatrix. This is in part consistent with a model in which the blob generation rate increases with the poloidal velocity shear.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗
DOE OSTI2022

Experimental validation of universal plasma blob formation mechanism

Anomalous plasma transport in the boundary region of a tokamak plasma is commonly associated with the formation and evolution of coherent density structures known as blobs. Recently, a theory for a universal mechanism of plasma blob formation has been put forward. It is based on a breaking process of a radially elongated streamer due to poloidal and radial velocity shears. The theory is well supported by two-dimensional and three-dimensional numerical simulation results but lacks experimental validation. Here, we report the first ever experimental validation of this universal criterion by testing it against NSTX data on blobs obtained using the gas-puff imaging diagnostic. It is found that the criterion is widely satisfied in most L-mode discharges and may explain the significantly larger number of blob events. We also validate the theoretical criterion against ADITYA Langmuir probe data taken in the scrape-off layer region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Modeling Offshore Wind Farm Performance in Coastal Low-Level Jets Using Coupled Mesoscale-Microscale Large Eddy Simulations

Accurately predicting wind farm reliability under complex offshore atmospheric conditions remains a key challenge, particularly during noncanonical meteorological events such as coastal low-level jets (LLJs). LLJs, characterized by strong nonmonotonic vertical shear and directional veer, depart significantly from the simplified inflow assumptions embedded in conventional design standards, low-fidelity engineering models, and microscale large eddy simulations of the atmospheric boundary layer. In this work, we use the virtual wind farm framework—an exascale, graphics processing unit–accelerated large eddy simulation platform coupled with high-fidelity aeroservoelastic turbine models and advanced mesoscale-microscale coupling via the ExaWind software stack—to investigate turbine responses under realistic LLJ forcing. Simulations are performed over the U.S. North Atlantic offshore domain with the use of meteorological inputs from New York State Energy Research and Development Authority buoy data, focusing on a representative LLJ case impacting the International Energy Agency 15 MW reference turbine. Our results show that LLJs can cause up to 50% power deficits in downstream turbine rows and significantly amplify low-speed shaft and tower loads through nonlinear coupling between complex inflow characteristics and turbine structural dynamics. Two primary mechanisms drive these load amplifications: (1) unique LLJ inflow features—including veer and vertical/lateral shear—and (2) the downstream evolution of the flow under stable thermal stratification, which suppresses turbulence mixing and alters wake recovery. These mechanisms produce streamwise variations in turbine loading not captured by standard hub height–based metrics or existing design load case (DLC) definitions. This study highlights the critical role of rotor-scale flow gradients in driving fatigue and system-level aeroelastic responses, challenging current DLC and control strategies. We advocate the integration of full-flow field, environment-aware wind inputs into load modeling and control algorithms. By leveraging exascale computing to resolve mesoscale-microscale coupling, this work lays the groundwork for next-generation offshore wind turbine design and operation in meteorologically complex marine environments.

17 WIND ENERGY↗
DOE OSTI2022

Behavior and mechanisms of Doppler wind lidar error in varying stability regimes

Abstract. Wind lidars are widespread and important tools in atmospheric observations. An intrinsic part of lidar measurement error is due to atmospheric variability in the remote-sensing scan volume. This study describes and quantifies the distribution of measurement error due to turbulence in varying atmospheric stability. While the lidar error model is general, we demonstrate the approach using large ensembles of virtual WindCube V2 lidar performing a profiling Doppler-beam-swinging scan in quasi-stationary large-eddy simulations (LESs) of convective and stable boundary layers. Error trends vary with the stability regime, time averaging of results, and observation height. A systematic analysis of the observation error explains dominant mechanisms and supports the findings of the empirical results. Treating the error under a random variable framework allows for informed predictions about the effect of different configurations or conditions on lidar performance. Convective conditions are most prone to large errors (up to 1.5 m s−1 in 1 Hz wind speed in strong convection), driven by the large vertical velocity variances in convective conditions and the high elevation angle of the scanning beams (62∘). Range-gate weighting induces a negative bias into the horizontal wind speeds near the surface shear layer (−0.2 m s−1 in the stable test case). Errors in the horizontal wind speed and direction computed from the wind components are sensitive to the background wind speed but have negligible dependence on the relative orientation of the instrument. Especially during low winds and in the presence of large errors in the horizontal velocity estimates, the reported wind speed is subject to a systematic positive bias (up to 0.4 m s−1 in 1 Hz measurements in strong convection). Vector time-averaged measurements can improve the behavior of the error distributions (reducing the 10 min wind speed error standard deviation to <0.3 m s−1 and the bias to <0.1 m s−1 in strong convection) with a predictable effectiveness related to the number of decorrelated samples in the time window. Hybrid schemes weighting the 10 min scalar- and vector-averaged lidar measurements are shown to be effective at reducing the wind speed biases compared to cup measurements in most of the simulated conditions, with time averages longer than 10 min recommended for best use in some unstable conditions. The approach in decomposing the error mechanisms with the help of the LES flow field could be extended to more complex measurement scenarios and scans.

17 WIND ENERGY↗