Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “S WAVES”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Identifying Ultra Low Frequency Waves in the Lunar Plasma Environment Using Trajectory Analysis and Resonance Conditions

Abstract Recent studies show that localized crustal magnetic fields on the lunar surface can reflect a significant portion of the incoming solar wind protons. These reflected ions can drive a wide range of plasma waves. It is difficult to determine the intrinsic properties of low-frequency waves with single-spacecraft observations, which can be heavily Doppler shifted. We describe a technique to combine trajectory analysis of reflected protons with the Doppler shift and resonance conditions to identify ultralow-frequency waves at the Moon. On 31 January 2014 plasma waves were detected by one of the Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun (ARTEMIS) probes as it approached the lunar wake; these waves were not detected by the second ARTEMIS probe located upstream in the undisturbed solar wind. The observed waves had a frequency below the local ion cyclotron frequency and had right-hand circular polarization in the reference frame of the Moon. By solving the Doppler shift and the cyclotron resonance equations, we determined the conditions for reflected ions to excite the observed waves. Simulated trajectories of reflected ions correspond to ARTEMIS ion observations and support the hypothesis that reflected ions are the primary driver of the waves. By combining trajectory analysis with the resonance conditions, we identify scenarios where ions that satisfy the resonance conditions are present in the right location to generate the observed waves. Using this method, we can uniquely identify the observed waves as upstream propagating right-hand polarized waves, subject to the assumption that they are generated by cyclotron resonance with ions.

Magnetic↗
NASA NTRS2012

Interpretation of Vector Electric Field Measurements of Bursty Langmuir Waves in the Cusp

An analysis of auroral Langmuir waves in the cusp observed by two high-frequency electric field instruments on the TRICE high-flyer sounding rocket shows many examples of Langmuir wave bursts modulated at approximately 10 kHz. Previous studies have explained these and similar observations as the result of beating between waves with very close frequencies near the Langmuir cutoff, resulting from either wave-wave interactions or independent linear excitations. The unique three-dimensional (3-D) data set provided by the NASA Goddard Space Flight Center TAEFWD instrument shows that up to 25% of waveforms selected from the most intense bursts exhibit anisotropic modulations, i.e., the beat nulls and peaks are not aligned in time across the three perpendicular electric field components. Modulations of this type arise when superposed wave normal modes possess differing polarizations, and simulations using wave modes calculated with the J-WHAMP numerical dispersion code show (1) that waves with differing polarizations can exist in conditions like those observed by TRICE-High and (2) that superpositions of such waves can produce anisotropic modulation. Fourier analysis of the 3-D waveform data suggests that both linear and elliptically polarized waves are present near the Langmuir cutoff at these times. These observations illustrate how 3-D measurements can give valuable insight into the nature of wave interactions in the auroral plasma environment.

Dombrowski, M. P.↗
NASA NTRS2019

Satellite Estimates of Momentum Fluxes from High-Impact Gravity Wave Events in the Stratosphere and Their Effects on Circulation

Recent assessments of chemistry-climate models (CCMs) reveal biases in temperatures and winds in, especially but not limited to, the Southern Hemisphere stratosphere, where winds are generally too strong and temperatures too cold. The reasons for these biases are not completely understood, but it is thought that missing wave drag in models is a major culprit. Observational and modeling studies support this idea by elucidating the role of infrequent but very high-impact gravity wave events in the stratosphere. These highly intermittent gravity wave events with large momentum fluxes are the most important drivers of circulation and transport in the stratosphere, yet they are not treated correctly in most global models. This has implications for the cold pole problem in the Southern Hemisphere and the global Brewer-Dobson circulation in general. In this presentation we show results combining HIRDLS and AIRS to derive detailed gravity wave properties and obtain new quantitative estimates of the local and intermittent gravity wave drag in the stratosphere. The combination of high-vertical resolution (1 km) and near-global (60°S to 80°N), close horizontal sampling (100 km) makes HIRDLS temperatures the best available dataset for retrieving gravity wave properties needed to diagnose gravity wave effects on circulation. We further exploit the close zonal sampling of HIRDLS near the turnaround latitude in the Southern Hemisphere to obtain estimates of the missing drag. We combine the HIRDLS results with AIRS brightness temperature images, which reveal high-spatial resolution detail of long vertical wavelength waves, to obtain 3-D, day-to-day variability in gravity wave properties and attribute the wave events to wave sources. The AIRS and HIRDLS datasets complement each other well since the two instruments have very different resolutions and horizontal sampling.

Holt, Laura↗
NASA NTRS2020

Mechanisms Associated with Daytime and Nighttime Heat Waves over the Contiguous United States

Heat waves are extreme climate events that have the potential to cause immense stress on human health, agriculture, and energy systems, so understanding the processes leading to their onset is crucial. There is no single accepted definition for heat waves, but they are generally described as a sustained amount of time where temperature exceeds a local threshold. Multiple different temperature variables are potentially relevant, as high values of both daily maximum (Tmax) and minimum (Tmin) temperatures can be detrimental to human health. In this study, we focus explicitly on the different mechanisms associated with summertime heat waves manifested during daytime versus nighttime hours over the contiguous United States. Heat waves are examined using the National Aeronautics and Space Administration (NASA) Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2). Over 1980–2018, the increase in the number of heat wave days per summer was generally stronger for nighttime heat wave days than daytime heat wave days, with localized regions of significant positive trends. Processes linked with daytime and nighttime heat waves are identified through composite analysis of precipitation, soil moisture, clouds, humidity and fluxes of heat and moisture. Daytime heat waves are associated with dry conditions, reduced cloud cover, and increased sensible heating. Mechanisms leading to nighttime heat waves differ regionally across the US, but they are typically associated with increased clouds, humidity and/or low-level temperature advection. In the Midwest US, enhanced moisture is transported from the Gulf of Mexico during nighttime heat waves.

Extreme events↗
NASA NTRSDate not supplied

Effect of Injection Dynamics on Wave Propagation in a Linear Detonation Combustor

Injector response to self-excited and sustained detonation wave propagation in an optically-accessible nonpremixed natural gas-oxygen linear detonation combustor is investigated. 100 kHz measurements of planar laser induced fluorescence (PLIF) of trace quantities of acetone injected into the fuel supply and simultaneous OH∗ chemiluminescence measurements were used to correlate the effect on injection dynamics on detonation wave structure and propagation. Measurements over a range of equivalence ratios (0.61-1.48) revealed wave propagation through regions of high fuel stratification in the axial and transverse direction of reactant injection. Consequently, the propagation of the detonation wave along the channel is through a flow-field with significant variation in fuel mass fraction, reactivity, and sound speeds. At lean conditions, wave propagation frequency is generally lower (∼6 kHz) resulting in longer time between wave passages for mixing and pre-heat resulting in steeper detonation wave-fronts anchored near the injector exit plane. Conversely, at rich equivalence ratios the wave frequency is higher (∼10 kHz) resulting in shorter intra-cycle period for reactant fill and mixing resulting waves lifted from the injection plane that are more compact. Higher chemiluminescence intensity in the post-wave region with deflagrative combustion is observed in these cases. Phase averaged representation of the acetone-PLIF measurements indicate that the fraction of time between wave passages available for reactant fill, mixing and pre-heat is consistent across all equivalence ratios.

Propulsion↗
NASA NTRSDate not supplied

Effect of Injection Dynamics on Detonation Wave Propagation in a Linear Detonation Combustor

Injector response to self-excited and sustained detonation wave propagation in an optically-accessible nonpremixed natural gas-oxygen linear detonation combustor is investigated. 100 kHz measurements of planar laser induced fluorescence (PLIF) of trace quantities of acetone injected into the fuel supply and simultaneous OH∗ chemiluminescence measurements were used to correlate the effect on injection dynamics on detonation wave structure and propagation. Measurements over a range of equivalence ratios (0.61-1.48) revealed wave propagation through regions of high fuel stratification in the axial and transverse direction of reactant injection. Consequently, the propagation of the detonation wave along the channel is through a flow-field with significant variation in fuel mass fraction, reactivity, and sound speeds. At lean conditions, wave propagation frequency is generally lower (∼6 kHz) resulting in longer time between wave passages for mixing and pre-heat resulting in steeper detonation wave-fronts anchored near the injector exit plane. Conversely, at rich equivalence ratios the wave frequency is higher (∼10 kHz) resulting in shorter intra-cycle period for reactant fill and mixing resulting waves lifted from the injection plane that are more compact. Higher chemiluminescence intensity in the post-wave region with deflagrative combustion is observed in these cases. Phase averaged representation of the acetone-PLIF measurements indicate that the fraction of time between wave passages available for reactant fill, mixing and pre-heat is consistent across all equivalence ratios.

Propulsion↗
DOE OSTI2024

Wave-powered water pump for upwelling in aquaculture: Numerical model and ocean test

Wave-powered upwelling can increase the productivity and survivability of several aquaculture species. This enhancement is due to transporting cold, nutrient-rich ocean water, typically found lower in the water column, to the surface. Macroalgaes, like kelp, exhibit increased growth from these altered conditions. The University of New Hampshire’s (UNH) wave-powered water pump (wave pump) is a point absorber wave energy converter (WEC) that uses ocean waves to create relative motion between a spar buoy and a concentric float which drives an internal pump. A numerical model of the wave pump was developed using WEC-Sim to predict device performance in the ocean. Wave pump performance was evaluated during a five day ocean test near Appledore Island in Maine in March 2023, where volumetric flow rate, relative distance between spar and float, and wave conditions were measured. These data were then used for numerical model validation. The ocean deployment recorded the device’s performance in a variety of sea states, with average significant wave heights up to 0.7 m. The ocean test data were compared to the WEC-Sim numerical model of the device with favorable results. Average values of device stroke period, stroke height, and flow rate agreed between the ocean test and model data to within approximately 16 to 22%. Furthermore, the validated numerical model provides a valuable tool for improving the design and developing a commercial-scale, wave-powered water pump for use in aquaculture.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Recurrence analysis of friction based dry-couplant ultrasonic Lamb waves in plate-like structures

Here, in this study, the effect of friction on the generation of dry-coupled Lamb waves is experimentally investigated. Recurrence analysis is performed to analyze the complex behavior of friction based dry-coupled Lamb waves. In particular, the effect of the normal force, which is necessary for a stronger dry-coupled Lamb wave generation and the friction, on the transmission of mechanical energy and determinism characteristics of Lamb waves are investigated. The results verify that larger friction coefficient and friction force are crucial for generation and propagation of strong Lamb waves supporting the fact that the main mechanism to transfer mechanical energy using dry-couplant is friction. The sensitivity of Lamb waves to the friction coefficient, highlights the importance of designing specific pads with respect to condition of the surface. Besides, the results show that the normal force and friction coefficient can change the determinism characteristics behavior of multimode Lamb waves. Furthermore, it is shown that the determinism value is sensitive to the friction coefficient and normal force. A similar trend is observed in the determinism values and friction coefficient. In general, a smaller friction coefficient indicates smaller determinism value. Additionally, it is shown that a normal load can change the behavior of a system, as observed from recurrence plots, owing to changes in the Lamb waves trajectories in the phase-space domain. In addition, it is shown that recurrence plots enable the detection of mode transitions in multimode Lamb waves. Recurrence analysis is a complementary tool to frequency domain methods for accurate analysis of multimode Lamb waves behavior.

42 ENGINEERING↗
DOE OSTI2022

Whistler waves generated by nongyrotropic and gyrotropic electron beams during asymmetric guide field reconnection

Using a two-dimensional particle-in-cell simulation of asymmetric reconnection with a guide field whose strength is 0.3 times the reconnecting magnetic field, we study electron distribution functions and wave intensities in the diffusion region, focusing on the electron diffusion region (EDR). Wave activities with frequencies below the electron cyclotron frequency are observed, and these are whistler waves propagating almost anti-parallel to the magnetic field. The waves are concentrated near the magnetospheric separatrix away from the X line, but the wave activity also spreads through the EDR near the X line. The reconnection outflows are asymmetric in the outflow direction in the magnetospheric side, and the wave intensity is stronger in the side of the faster electron outflow. We study the whistler waves using the fast Fourier transform, analyses of electron velocity distribution functions, and the dispersion solver calculation. Along the magnetospheric separatrix in the stronger outflow side, highly anisotropic electron beams exist with super-Alfvénic drift speeds. The dispersion analysis shows that there are two modes: a temperature anisotropy mode and a beam mode. Outside the EDR, the whistler wave intensity is highest near the separatrix, but the wave intensity decreases if we move away from the separatrix toward the magnetic neutral line because of the increase in the electron population near zero parallel velocity. In the EDR, in the velocity plane perpendicular to the magnetic field, ring/crescent electron distribution functions are observed. Near the X-line, the wave power is enhanced where nongyrotropic electrons contribute to increase the perpendicular temperature anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Uncertainty Assessment of CFD Investigation of the Nonlinear Difference-Frequency Wave Loads on a Semisubmersible FOWT Platform

Current mid-fidelity modeling approaches for floating offshore wind turbines (FOWTs) have been found to underpredict the nonlinear, low-frequency wave excitation and the response of semisubmersible FOWTs. To examine the cause of this underprediction, the OC6 project is using computational fluid dynamics (CFD) tools to investigate the wave loads on the OC5-DeepCwind semisubmersible, with a focus on the nonlinear difference-frequency excitation. This paper focuses on assessing the uncertainty of the CFD predictions from simulations of the semisubmersible in a fixed condition under bichromatic wave loading and on establishing confidence in the results for use in improving mid-fidelity models. The uncertainty for the nonlinear wave excitation is found to be acceptable but larger than that for the wave-frequency excitation, with the spatial discretization error being the dominant contributor. Further, unwanted free waves at the difference frequency have been identified in the CFD solution. A wave-splitting and wave load-correction procedure are presented to remove the contamination from the free waves in the results. A preliminary comparison to second-order potential-flow theory shows that the CFD model predicted significantly higher difference-frequency wave excitations, especially in surge, suggesting that the CFD results can be used to better calibrate the mid-fidelity tools.

17 WIND ENERGY↗
DOE OSTI2021

Prototype and Codesign of Nascent Flexible Wave Energy Converter Concepts

Prevailing mainstream ocean wave energy converter (WEC) archetypes are often relatively large monolithic rigid bodies with singular means of power take-off (e.g., relative motion between rigid bodies) and solo prime mover mechanisms (e.g. a rotary generator, a hydraulic piston). These mainstream WEC archetypes have existed for decades with some having histories spanning more than a century of significant investments and interests. As an alternative to the domain of mainstream WEC archetypes is the domain of distributed embedded energy converter technologies (also known as DEEC-Tec). Use of DEEC-Tec to create WECs - often called flexWECs - is a nascent domain of research and development that centers itself upon the interconnection of many relatively small energy converters - very small energy transducers - that interconnect and combine to create a much larger, overall, energy converting structure. The DEEC-Tec domain enables promising possibilities such as: direct in situ ocean wave energy harvesting and conversion throughout the volume of an entire flexWEC; inherent broad-banded ocean wave energy capture and conversion arising from near semi-continuous structural degrees of freedom; an innate absence of focusing and concentrating of ocean wave energy into a single prime mover and/or transmission system; and built-in implicit redundancy. As such, the DEEC-Tec domain could revolutionize how we conceptualize ocean wave energy conversion - both in terms of ocean wave energy converter design and what we view as viable environments and circumstances for ocean wave energy conversion. The research and development opportunities for the DEEC-Tec domain, therefore, are vast compared to the prevailing mainstream ocean wave energy converter domain. Correspondingly, there are ample opportunities to invest time and money into the DEEC-Tec domain with an aim to enable a paradigm shift in how ocean waves can be utilized as a more viable form of renewable energy. To that end, two seedling efforts have been enacted providing a very cursory investigation into the DEEC-Tec domain: (1) a preliminary simple fabrication/prototyping effort and (2) a preliminary co-design effort. The prototype seedling centers upon the design of a fabrication and evaluation rig for individual DEEC-Tec transducers; e.g., dielectric elastomer generators. The co-design seedling centers upon the analysis of those methods suitable for DEEC-Tec, ocean wave energy conversion, and fluid-structure interaction co-design efforts.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2026

Measurement of the Alfvén Wave Parametric Decay Instability Growth Rate

Alfvén waves, a fundamental mode of magnetized plasmas, are ubiquitous in space and laboratory plasmas. The nonlinear behavior of these modes is thought to play a key role in important problems in space plasma, such as the heating of the solar corona and solar wind turbulence. In particular, theoretical predictions show that these Alfvén waves may be unstable to various parametric instabilities, but space observations of these processes are limited. We demonstrate the first measurement of the Alfvén wave parametric decay instability (PDI) growth rate. Experiments are conducted on the Large Plasma Device at UCLA in which a high amplitude 𝛿⁢𝐵/𝐵 0 ∼ 0.7% pump Alfvén wave is launched from one end of the device and a smaller seed Alfvén wave is launched from the other side. When the frequency of the seed wave is chosen to match the backward wave expected from PDI, damping of the seed wave is reduced. We compare this reduction in damping to the theoretically expected PDI growth rate while accounting for acoustic mode damping. Results show agreement between measurements and theoretical predictions. As a result, this not only provides critical validation for PDI theories and simulations that could help interpret future space observations but also suggests a new way of studying similar nonlinear wave phenomena.

Alfvén waves↗
DOE OSTI2021

Exploring wave propagation in heterogeneous metastructures using the relaxed micromorphic model

Metamaterials are artificial structures that can manipulate and control sound waves in ways not possible with conventional materials. While much effort has been undertaken to widen the bandgaps produced by these materials through design of heterogeneities within unit cells, comparatively little work has considered the effect of engineering heterogeneities at the structural scale by combining different types of unit cells. In this paper, we use the relaxed micromorphic model to study wave propagation in heterogeneous metastructures composed of different unit cells. We first establish the efficacy of the relaxed micromorphic model for capturing the salient characteristics of dispersive wave propagation through comparisons with direct numerical simulations for two classes of metamaterial unit cells: namely phononic crystals and locally resonant metamaterials. We then use this model to demonstrate how spatially arranging multiple unit cells into metastructures can lead to tailored and unique properties such as spatially-dependent broadband wave attenuation, rainbow trapping, and pulse shaping. In the case of the broadband wave attenuation application, we show that by building layered metastructures from different metamaterial unit cells, we can slow down or stop wave packets in an enlarged frequency range, while letting other frequencies through. In the case of the rainbow-trapping application, we show that spatial arrangements of different unit cells can be designed to progressively slow down and eventually stop waves with different frequencies at different spatial locations. Finally, in the case of the pulse-shaping application, our results show that heterogeneous metastructures can be designed to tailor the spatial profile of a propagating wave packet. Collectively, these results show the versatility of the relaxed micromorphic model for effectively and accurately simulating wave propagation in heterogeneous metastructures, and how this model can be used to design heterogeneous metastructures with tailored wave propagation functionalities.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Wave energy converter buoy for Arctic observations

Ocean observation buoys are currently powered-constrained by battery storage capacity or available solar power. Power constraints limit the number of measurements that can be made and the lifetime of the buoys. Ocean surface wave energy could be used for power production, but wave energy converters are not yet commercially available for ocean observation buoys. Here we present the design and testing of a drifting wave energy converter buoy using a pendulum transmission system (PTS). The wave energy converter buoy was designed to operate in Arctic temperatures and wave conditions, but it could be used in areas with warmer temperatures and larger waves. Field tests measured a maximum power production of 5 watts in waves with a significant wave height of 0.42 meters and a ten minute average power production of 37 milliwatts. Power production increased with the energy period and significant wave height. The energy harvesting capabilities of the PTS showed the utility of adding a wave energy conversion device to a drifting ocean observation buoy.

Arctic observations↗
DOE OSTI2024

On hydromagnetic wave interactions in collisionless, high- β plasmas

We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high- β collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew–Goldberger–Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfvén wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfvén wave. The difference in time scales of acoustic and Alfvénic fluctuations at high- β means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfvénic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfvén-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfvén-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfvénic turbulence and to the interaction between the compressive and Alfvénic cascades in strong, collisionless turbulence are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Generation of low–latitude seamount–trapped waves: a case study of the Seychelles Plateau

Baroclinic seamount-trapped waves are thought to influence their surrounding ecosystem; however, trapped waves are not well-studied in near-equatorial settings, where stratification is strong and Burger numbers are large. Motivated by observations, we use daily output (2005-2009) from the global 0.1-degree Parallel Ocean Program Model (POP) to examine topographically trapped baroclinic waves around the Seychelles Plateau (S > 400) in the tropical Indian Ocean. These trapped waves are associated with velocity and temperature oscillations at periods of 15–16 days, similar to the dominant period of some equatorial Yanai waves. Energy flux maps using POP output suggest that quasi-biweekly equatorial Yanai waves excite trapped waves on the western and south-western flanks of the Seychelles Plateau, near the surface. Here, the anticyclonic energy flux associated with the trapped wave extends vertically throughout the water column and around most of the plateau circumference, diminishing on the eastern flank of the plateau. This work highlights the role that equatorial planetary waves and trapped waves play in facilitating energy redistribution, dissipation, and mixing in the tropical ocean.

58 GEOSCIENCES↗
DOE OSTI2022

Whistler Waves Associated With Electron Beams in Magnetopause Reconnection Diffusion Regions

Abstract Whistler waves are often observed in magnetopause reconnection associated with electron beams. We analyze seven MMS crossings surrounding the electron diffusion region (EDR) to study the role of electron beams in whistler excitation. Waves have two major types: (a) Narrow‐band waves with high ellipticities and (b) broad‐band waves that are more electrostatic with significant variations in ellipticities and wave normal angles. While both types of waves are associated with electron beams, the key difference is the anisotropy of the background population, with perpendicular and parallel anisotropies, respectively. The linear instability analysis suggests that the first type of wave is mainly due to the background anisotropy, with the beam contributing additional cyclotron resonance to enhance the wave growth. The second type of broadband waves are excited via Landau resonance, and as seen in one event, the beam anisotropy induces an additional cyclotron mode. The results are supported by particle‐in‐cell simulations. We infer that the first type occurs downstream of the central EDR, where background electrons experience Betatron acceleration to form the perpendicular anisotropy; the second type occurs in the central EDR of guide field reconnection. A parametric study is conducted with linear instability analysis. A beam anisotropy alone of above ∼3 likely excites the cyclotron mode waves. Large beam drifts cause Doppler shifts and may lead to left‐hand polarizations in the ion frame. Future studies are needed to determine whether the observation covers a broader parameter regime and to understand the competition between whistler and other instabilities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications

Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind-farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, which contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind-farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches relies heavily on ad hoc processes. For instance, the widely used Rayleigh damping method requires ad hoc tuning to produce a setup that may only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct a large-eddy simulation (LES) study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics. In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This, in turn, will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths (the representative wavelengths of the entire wave spectrum) to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency. Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upward-propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow–outflow boundaries is of equal importance to damping the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.

17 WIND ENERGY↗