SEARCH · Engineering Papers
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.
This dataset contains streamwise velocity retrieved from the Range Height Indicator (RHI) scans done by the UC Davis scanning Lidar. The quality control of the radial wind speed is performed following the algorithm of Beck & Kuhn (2017). Instantaneous realizations characterized by wind misalignment greater than 60 degrees (based on the closest VAD scan in time) are discarded.
The Amazon rainforest is a vital component of the global climate system, influencing the hydrological cycle and tropical circulation. However, understanding and modeling the evolution of convection in this region remain a scientific challenge. Here, we assess the environmental conditions associated with shallow, congestus, and isolated deep convection days during the wet season (December to April), employing measurements from the Green Ocean Amazon 2014–2015 (GoAmazon2014/5) experiment and large-scale wind fields from the constrained variational analysis. Composites of deep days show moister than average conditions below 3 km early in the morning. Analyzing the water budget at the surface through observations only, we estimated the water vapor convergence term as a residual of the water balance closure. Convergence remains nearly zero during the deep days until early afternoon (13:00 LST), when it becomes a dominant factor in the water budget. At 14:00 LST, the deep days experience a robust upward large-scale vertical velocity, especially above 4 km, which supports the shallow-to-deep convective transition occurring around 16:00–17:00 LST. In contrast, shallow and congestus days exhibit drier pre-convective conditions, along with diurnal water vapor divergence and large-scale subsidence that extend from the surface to the lower free troposphere. Moreover, afternoon precipitation exhibits the strongest linear correlation (0.6) with large-scale vertical velocity, nearly double the magnitude observed for other environmental factors, even moisture, at different levels and periods of the day. Precipitation also exhibits a moderate increase with low-level wind shear, while upper-level shear has a relatively minor negative impact on convection.
Magnetic holes (MHs) are coherent structures characterized by a strong and localized magnetic field amplitude dip, commonly observed in the heliosphere. These structures come in different sizes, from magnetohydrodynamic to kinetic scales. Subion-scale MHs are usually sustained by an electron current vortex and exhibit a strong electron temperature anisotropy, with higher temperatures perpendicular to the background magnetic field. Magnetospheric multiscale observations (MMSs) have revealed electron-scale MHs to be ubiquitous in the turbulent Earth’s magnetosheath and the solar wind, potentially playing an important role in the energy cascade and dissipation. Despite abundant observations, the origin of electron-scale MHs is still unclear and debated. In this work, we use fully kinetic simulations to investigate the role of plasma turbulence in generating electron-scale MHs. We find that the turbulence spontaneously produces electron-scale MHs via the following mechanism: first, large-scale turbulent velocity shears produce regions with high electron temperature anisotropy; these localized regions become unstable, generating oblique electron-scale whistler waves; as they propagate over the inhomogeneous turbulent background, whistler fluctuations develop an electrostatic component, turning into Bernstein-like modes; the strong electrostatic fluctuations produce current filaments that merge into an electron-scale current vortex; the resulting electron vortex locally reduces the magnetic field amplitude, finally evolving into an electron-scale MH. We show that MHs generated by this mechanism have properties consistent with MMSs and nontrivial kinetic features with a “mushroom”-shaped electron velocity distribution function. Our results have potential implications for understanding the formation and occurrence of electron-scale MHs in astrophysical turbulent and space environments, such as the Earth’s magnetosheath and the solar wind.
We study the evolution of collisionless plasmas that, due to their macroscopic evolution, are susceptible to the firehose instability, using both analytic theory and hybrid-kinetic particle-in-cell simulations. We establish that, depending on the relative magnitude of the plasma β, the characteristic time scale of macroscopic evolution and the ion-Larmor frequency, the saturation of the firehose instability in high-β plasmas can result in three qualitatively distinct thermodynamic (and electromagnetic) states. By contrast with the previously identified ‘ultra-high-beta’ and ‘Alfvén-inhibiting’ states, the newly identified ‘Alfvén-enabling’ state, which is realised when the macroscopic evolution time τ exceeds the ion-Larmor frequency by a β-dependent critical parameter, can support linear Alfvén waves and Alfvénic turbulence because the magnetic tension associated with the plasma’s macroscopic magnetic field is never completely negated by anisotropic pressure forces. We characterise these states in detail, including their saturated magnetic-energy spectra. The effective collision operator associated with the firehose fluctuations is also described; we find it to be well approximated in the Alfvén-enabling state by a simple quasi-linear pitch-angle scattering operator. The box-averaged collision frequency is ν eff ∼ β/τ, in agreement with previous results, but certain subpopulations of particles scatter at a much larger (or smaller) rate depending on their velocity in the direction parallel to the magnetic field. Our findings are essential for understanding low-collisionality astrophysical plasmas including the solar wind, the intracluster medium of galaxy clusters and black hole accretion flows. We show that all three of these plasmas are in the Alfvén-enabling regime of firehose saturation and discuss the implications of this result.
The main goal of the TRacking Aerosol Convection interactions ExpeRiment (TRACER) campaign was to study aerosol–cloud interactions during deep convection over the Houston area. This project deployed a suite of instrumentation with the aim to (1) quantify turbulent vertical particle fluxes during at DOE-ARM sites, including TRACER, (2) assess hygroscopic growth factors and hygroscopicity parameters of the material driving modal aerosol growth during new particle formation and growth events, (3) derive turbulent aerosol mass fluxes using co-located Doppler LIDAR measurements, and (4) create quality-controlled PI data products to support future research utilizing data collected during the TRACER campaign. This report summarized the main findings from the deployments at two DOE-ARM sites. Briefly, we found that new particle formation may occur aloft, in a residual layer, near the top of the boundary layer. Small grown particles appear later due to downward mixing with daytime turbulence. The species that are responsible for aerosol modal growth had hygroscopicity parameters varying between 0.05 and 0.34. These values systematically depended on the wind sector, suggesting that the chemical composition of the precursors differed. This work demonstrated that lidar retrievals of the elastic backscatter and Doppler velocity can be used to obtain surface number emissions of particles with a diameter greater than 0.53 µm. During TRACER, emission particle number fluxes peaked near ∼ 100 cm−2 s−1. Multiple quality-controlled PI data products that will support future TRACER related science were generated and made publically available.
As wind energy deployment grows, interactions between wind plants and the surrounding environment become more prevalent. The current investigation seeks to understand these interactions by characterizing the impact of wind plants on the planetary boundary layer height (PBLH), utilizing observations from the American WAKE experimeNt (AWAKEN) campaign. Given the ambiguity of the definition of PBLH under stable atmospheric conditions, where the impact of wind plants is expected to be strongest, a comparison of different methods for identifying PBLH is first conducted using data collected by multiple types of instruments. The Heffter method is selected as the thermodynamic method because it generates the most consistent results for the radiosonde and infrared spectrometer. A minimum vertical velocity variance method is used for a turbulence-based definition. Using both of these methods, the values of PBLH measured at spatially distributed sites are compared under a range of atmospheric conditions. Both methods show a clear increase in PBLH downstream of a wind plant for stable conditions. These impacts are strongest when the upstream PBLH is shallow (less than 0.25 km), with the thermodynamic method showing a PBLH increase of 35 %–39 % and the turbulence-based method showing a 143 % increase. At a site 20 km downstream of the wind plant, these effects are no longer observed, suggesting PBLH is not influenced by the wind plant at this distance. The results of this investigation show that wind plants can modify PBLH in their vicinity. As PBLH is a key parameter for numerical models, this insight into wind plant–atmosphere interactions can be used to validate and improve the accuracy of wind energy forecasts.
The Kelvin-Helmholtz Instability (KHI), arising from velocity shear across the magnetopause, plays a significant role in the viscous-like transfer of mass, momentum, and energy from the shocked solar wind into the magnetosphere. While the KHI leads to growth of surface waves and vortices, suitable detection methods for these applicable to magnetohydrodynamics (MHD) are currently lacking. A novel method is derived based on the well-established λ-family of hydrodynamic vortex identification techniques, which define a vortex as a local minimum in an adapted pressure field. The J × B Lorentz force is incorporated into this method by using an effective total pressure in MHD, including both magnetic pressure and a pressure-like part of the magnetic tension derived from a Helmholtz decomposition. The λ MHD method is shown to comprise of four physical effects: vortical momentum, density gradients, fluid compressibility, and the rotational part of the magnetic tension. A local three-dimensional MHD simulation representative of near-flank magnetopause conditions (plasma β’s 0.5 – 5 and convective Mach numbers M f ∼ 0.4) under northward interplanetary magnetic field (IMF) is used to validate λ MHD . Analysis shows it correlates well with hydrodynamic vortex definitions, though the level of correlation decreases with vortex evolution. Overall, vortical momentum dominates λ MHD at all times. During the linear growth phase, density gradients act to oppose vortex formation. By the highly nonlinear stage, the formation of small-scale structures leads to a rising importance of the magnetic tension. Compressibility was found to be insignificant throughout. Finally, a demonstration of this method adapted to tetrahedral spacecraft observations is performed.
Abstract We study the magnetospheric evolution of a nonaccreting spinning black hole (BH) with an initially inclined split monopole magnetic field by means of 3D general relativistic magnetohydrodynamic simulations. This serves as a model for a neutron star (NS) collapse or a BH–NS merger remnant after the inherited magnetosphere has settled into a split monopole field creating a striped wind. We show that the initially inclined split monopolar current sheet aligns over time with the BH equatorial plane. The inclination angle evolves exponentially toward alignment, with an alignment timescale that is inversely proportional to the square of the BH angular velocity, where higher spin results in faster alignment. Furthermore, magnetic reconnection in the current sheet leads to exponential decay of event-horizon-penetrating magnetic flux with nearly the same timescale for all considered BH spins. In addition, we present relations for the BH mass and spin in terms of the period and alignment timescale of the striped wind. The explored scenario of a rotating, aligning, and reconnecting current sheet can potentially lead to multimessenger electromagnetic counterparts to a gravitational-wave event due to the acceleration of particles powering high-energy radiation, plasmoid mergers resulting in coherent radio signals, and pulsating emission due to the initial misalignment of the BH magnetosphere.
This dataset contains uncorrected sap velocity measurements using the heat ratio method (HRM) collected using ICT International SFM1x sensors at five urban sites in Chicago, Illinois, as part of the DOE CROCUS project. The data includes continuous monitoring of sap velocity from various tree species, including Maples (Acer spp.): Sugar Maple (Acer saccharum), Silver Maple (Acer saccharinum), and Red Maple (Acer rubrum); Oaks (Quercus spp.): Swamp White Oak (Quercus bicolor); American Elm (Ulmus americana); Honey Locust (Gleditsia triacanthos); Cottonwood (Populus deltoides); and Tree of Heaven (Ailanthus altissima) across Chicago State University (CSU), Northeastern Illinois University (NEIU), Northwestern University (NU), University of Illinois Chicago (UIC), and West Woodlawn "Blacks in Green" (BIG). These include both street trees and those in urban park locations. Measurements were collected at 15-20 minute intervals, depending on the sensor, and transmitted via Long Range Wide Area Network (LoRaWAN) protocols. The wireless data was collected by Sage Network (https://sagecontinuum.org/) nodes. The dataset includes sensor ID, Global Positioning System (GPS) coordinates, tree species (common and scientific names), tree identification number, diameter at breast height (DBH in cm), uncorrected sap velocity measurements (cm/hr) from both inner and outer probes, and Sage Node identifiers so the data can be mapped to related variables such as air quality and wind speed that were collected on the Sage nodes. All timestamps are in local Chicago time (CDT/CST). Quality control flags are provided using a 3-bit binary system indicating physical range violations (< -10 or > 60 cm/hr), step spikes (absolute difference > 36 cm/hr), and stuck sensor conditions (> 10 consecutive identical values). These are raw data, not corrected for wood anatomy or species-specific characteristics. Data is provided in comma separated (CSV) format. This dataset is part of a larger collection of CROCUS environmental monitoring data, including linked datasets from Air Quality Transmitter (AQT) sensors, Weather Transmitter (WXT) sensors, and Multi-Function Research LoRaWAN (MFR) Nodes. DOIs for the supporting data are provided as part of this data package.
We use in situ measurements from the first 19 encounters of Parker Solar Probe and the most recent five encounters of Solar Orbiter to study the evolution of the turbulent sonic Mach number M t (the ratio of the amplitude of velocity fluctuations to the sound speed) with radial distance and its relationship to density fluctuations. We focus on the near-Sun region with radial distances ranging from about 11 to 80 R ⊙ . Our results show that (1) the turbulent sonic Mach number M t gradually moves toward larger values as it approaches the Sun, until at least 11 R ⊙ , where M t is much larger than the previously observed value of 0.1 at and above 0.3 au; (2) transonic turbulence with M t ∼ 1 is observed in situ for the first time and is found mostly near the Alfvén critical surface; (3) Alfvén Mach number of the bulk flow M A shows a strong correlation with the plasma beta, indicating that most of the observed sub-Alfvénic intervals correspond to a low-beta plasma; (4) the scaling relation between density fluctuations and M t gradually changes from a linear scaling at larger radial distances to a quadratic scaling at smaller radial distances; and (5) transonic turbulence is more compressible than subsonic turbulence, with enhanced density fluctuations and slightly flatter spectra than subsonic turbulence. A systematic understanding of compressible turbulence near the Sun is necessary for future solar wind modeling efforts.
This numerical investigation employs Large Eddy Simulation (LES) coupled with Actuator Disk Model (ADM) to evaluate wind farm layout optimization strategies. The study presents a systematic analysis of aligned, horizontal staggering, vertical staggering, and mixed (combination of horizontal and vertical) staggering configurations, aiming to establish optimal design parameters for enhanced power production. The investigation examines key performance metrics including mean velocity distributions, turbulence intensity characteristics, and power generation efficiency. Results demonstrate better performance of both horizontal and vertical staggering patterns compared to conventional aligned configurations, with horizontal staggering exhibiting notably higher power output than vertical arrangements. Our findings also suggest that mixed configurations, incorporating both horizontal and vertical staggering, can offer optimal performance characteristics. As a result, this research advances the understanding of wake interactions in complex wind farm layouts and provides design guidelines for maximizing wind farm power generation efficiency through strategic turbine positioning.
Wind and marine energy are rapidly growing and complementary technologies that share some techniques for simplified modeling and control, particularly in below-rated flow speeds. A turbine operator has several choices of controller for maximizing power in Region 2. The simple and ubiquitous KΩ 2 control law is often effective but limited in its flexibility. Alternative controllers use reference tracking to split the control objectives into a low-bandwidth optimal tip-speed ratio tracking loop to maximize steady-state power and a higher-bandwidth proportional-integral control loop to reject inflow turbulence. Several options exist for identifying the slowly varying optimal set point during operation, based on estimating the inflow velocity or filtering the power or torque signals. This study compares the trade-offs between performance and other design priorities for a few choices of reference-tracking controller in the literature for reference wind and marine turbines. Analysis is performed in the frequency domain using the linearization of each controller, and the impact of turbulent disturbances on the closed-loop system is described. The controllers are simulated in OpenFAST to analyze their performance with higher-order nonlinear turbine dynamics.
The X-ray binary system Cygnus X-3 (4U 2030+40, V1521 Cyg) is luminous but enigmatic owing to the high intervening absorption. High-resolution X-ray spectroscopy uniquely probes the dynamics of the photoionized gas in the system. In this Letter, we report on an observation of Cyg X-3 with the XRISM/Resolve spectrometer, which provides unprecedented spectral resolution and sensitivity in the 2–10 keV band. We detect multiple kinematic and ionization components in absorption and emission whose superposition leads to complex line profiles, including strong P Cygni profiles on resonance lines. The prominent Fe xxv Heα and Fe xxvi Lyα emission complexes are clearly resolved into their characteristic fine-structure transitions. Self-consistent photoionization modeling allows us to disentangle the absorption and emission components and measure the Doppler velocity of these components as a function of binary orbital phase. We find a significantly higher velocity amplitude for the emission lines than for the absorption lines. The absorption lines generally appear blueshifted by ∼−500–600 km s −1 . We show that the wind decomposes naturally into a relatively smooth and large-scale component, perhaps associated with the background wind itself, plus a turbulent, denser structure located close to the compact object in its orbit.
Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.
W49B is a unique Galactic supernova remnant with centrally peaked, “bar”-like ejecta distribution, which was once considered evidence for a hypernova origin that resulted in a bipolar ejection of the stellar core. However, chemical abundance measurements contradict this interpretation. Closely connected to the morphology of the ejecta is its velocity distribution, which provides critical details for understanding the explosion mechanism. We report the first ever observational constraint on the kinematics of the ejecta in W49B using the Resolve microcalorimeter spectrometer on the X-ray Imaging and Spectroscopy Mission (XRISM). Using XRISM/Resolve, we measured the line-of-sight velocity traced by the Fe Heα emission, which is the brightest feature in the Resolve spectrum, to vary by ±300 km s −1 with a smooth east-to-west gradient of a few tens of kilometers per second per parsec along the major axis. Similar trends in the line-of-sight velocity structure were found for other Fe-group elements Cr and Mn, traced by the Heα emission, and also for intermediate-mass elements Si, S, Ar, and Ca, traced by the Lyα emission. The discovery of the east–west gradient in the line-of-sight velocity, together with the absence of a twin-peaked line profile or enhanced broadening in the central region, clearly rejects the equatorially expanding disk model. In contrast, the observed velocity structure suggests bipolar flows reminiscent of a bipolar explosion scenario. An alternative scenario would be a collimation of the ejecta by an elongated cavity sculpted by bipolar stellar winds.
We conduct 3D magnetohydrodynamic simulations of decaying turbulence in the context of the solar wind. To account for the spherical expansion of the solar wind, we implement the expanding box model. The initial turbulence comprises uncorrelated counterpropagating Alfvén waves and exhibits an isotropic power spectrum. Our findings reveal the consistent generation of negative residual energy whenever nonlinear interactions are present, independent of the normalized cross helicity σc and compressibility. The spherical expansion facilitates this process. The resulting residual energy is primarily distributed in the perpendicular direction, with S 2 (b) − S 2 (u) ∝ l ⊥ or equivalently –E r ∝ k$^{–2}_{⟂}$. Here S 2 (b) and S 2 (u) are second-order structure functions of magnetic field and velocity respectively. In most runs, S 2 (b) develops a scaling relation S 2 (b) ∝ $l$$^{1/2}_{⊥}$ (E b ∝ k$^{–3/2}_{⟂}$). In contrast, S 2 (u) is consistently shallower than S 2 (b), which aligns with in situ observations of the solar wind. We observe that the higher-order statistics of the turbulence, which act as a proxy for intermittency, depend on the initial σc and are strongly affected by the expansion effect. Generally, the intermittency is more pronounced when the expansion effect is present. Finally, we find that in our simulations, although the negative residual energy and intermittency grow simultaneously as the turbulence evolves, the causal relation between them seems to be weak, possibly because they are generated on different scales.
Abstract. Scanning lidars enable the collection of spatially distributed measurements of turbine wakes and the estimation of wake properties such as magnitude, extent, and trajectory. Lidar-based characterizations, however, may be subject to distortions due to the observational system. Distortions can arise from the resolution of the measurement points across the wake, the projection of the winds onto the beam, averaging along the beam probe volume, and intervening evolution of the flow over the scan duration. Using a large-eddy simulation and simulated measurements with a virtual lidar model, we assess how scanning lidar systems may influence the properties of the retrieved wake using a case study from the Perdigão campaign. We consider three lidars performing range-height indicator sweeps in complex terrain, based on the deployments of lidars from the Danish Technical University (DTU) and German Aerospace Center (DLR) at the Perdigão site. The unwaked flow, measured by the DTU lidar, is well-captured by the lidar, even without combining data into a multi-lidar retrieval. The two DLR lidars measure a waked transect from different downwind vantage points. In the region of the wake, the observation system reacts to the smaller spatial and temporal variations of the winds, allowing more significant observation distortions to arise. While the measurements largely capture the wake structure and trajectory over its 4–5 D extent, limited spatial resolution of measurement points and volume averaging lead to a quicker loss of the two lobes in the near wake, smearing of the vertical bounds of the wake (< 30 m), wake center displacements up to 10 m, and dampening of the maximum velocity deficit by up to a third. The virtual lidar tool, coupled with simulations, provides a means for assessing measurement capabilities in advance of measurement campaigns.
Acoustic tomography offers path-integrated measurements of atmospheric velocity and temperature fluctuations with high spatial resolution. Classical implementations of time-dependent stochastic inversion rely on homogeneous, isotropic covariance models that are poorly suited to the anisotropic structure of wind turbine wakes. By directly estimating heterogeneous covariances from large-eddy simulations (LESs) into the time-dependent stochastic inversion operator, we relax implicit assumptions in the analytical models used historically. Retrievals using these LES-informed models improve agreement with true fields in variance, turbulent kinetic energy, and spectral content compared to analytical and precursor-based covariance models. The results indicate that LES-informed covariance models can enhance the accuracy of acoustic tomography retrievals in complex, anisotropic flows such as wind turbine wakes in some cases and highlight instances where analytical models still offer competitive performance, despite their simplifying assumptions.