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.
High frequency (HF) radar systems measure the speed and direction of ocean surface currents in near real time provided by Integrated Ocean Observing System (IOOS). The HFR data along the Pacific Coast are at a spatial resolution of 6 km. These data are collected as part of an observational database developed to support the floating offshore wind energy research within the ORACLE project funded by DOE WETO.
Aerosols influence fog properties such as visibility and lifetime by affecting fog droplet number concentrations (N d ). Numerical weather prediction (NWP) models often represent aerosol–fog interactions using highly simplified approaches. Incorporating prognostic size-resolved aerosol microphysics from climate models could allow them to simulate N d and aerosol–fog interactions without incurring excessive computational expense. However, microphysics code designed for coarse spatial resolution may struggle with sub-kilometer-scale grid spacings. Here, we test the ability of the UK Met Office Unified Model to simulate aerosol and fog properties during case studies from the ParisFog field campaign in 2011. We examine the sensitivity of fog properties to variations in N d caused by modifications to simulated aerosol activation. Our model, with a 500 m horizontal resolution and interactive aerosol and cloud microphysics, significantly underpredicts N d , although it only slightly underestimates the cloud condensation nuclei concentration. With an updated version of the Abdul-Razzak and Ghan (2000) activation scheme, we produce N d that are more consistent with those predicted by a cloud parcel model under fog-like conditions. We activate droplets only by adiabatic cooling. We incorporate more realistic hygroscopicities for sulfate and organic aerosols and explore the sensitivity of simulated N d to unresolved updrafts. We find that both N d and simulated fog liquid water content are very sensitive to the updated activation scheme but remain less affected by the update to hygroscopicities. Our improvements offer insights into the physical processes regulating N d in stable conditions, potentially laying foundations for improved operational fog forecasts that incorporate interactive aerosol simulations or aerosol climatologies.
Crystal-strain variation imposes significant limitations on many quantum sensing and information applications for solid-state defect qubits in diamond. Thus, the precision measurement and control of diamond crystal strain is a key challenge. Here, we report diamond strain measurements with a unique set of capabilities, including micron-scale spatial resolution, a millimeter-scale field of view, and a 2-order-of-magnitude improvement in volume-normalized sensitivity over previous work, reaching 5(2)×10 -8 /√Hzμm -3 (with spin-strain coupling coefficients representing the dominant systematic uncertainty). We use strain-sensitive spin-state interferometry on ensembles of nitrogen-vacancy (N-V) color centers in single-crystal bulk diamond with low strain gradients. This quantum interferometry technique provides insensitivity to magnetic-field inhomogeneity from the electronic and nuclear spin bath, thereby enabling long N-V–ensemble electronic spin dephasing times and enhanced strain sensitivity, as well as broadening the potential applications of the technique beyond isotopically enriched or high-purity diamond. We demonstrate the strain-sensitive measurement protocol first on a confocal scanning laser microscope, providing quantitative measurement of sensitivity as well as three-dimensional strain mapping; and second on a wide-field-imaging quantum diamond microscope. Our strain-microscopy technique enables fast, sensitive characterization for diamond material engineering and nanofabrication; as well as diamond-based sensing of strains applied externally, as in diamond anvil cells or embedded diamond stress sensors, or internally, as by crystal damage due to particle-induced nuclear recoils.
Many critical phenomena in soft matter occur at large length scales, necessitating the resolution of their structure and dynamics at low wavenumbers. However, resolving wavenumber-dependent dynamics computationally via molecular dynamics simulations presents significant challenges, as these phenomena span several orders of magnitude in both time and length scales, resulting in high computational costs and memory demands. Here, this work highlights the computational and memory challenges associated with analyzing molecular trajectories in reciprocal space and demonstrates a method to address them. We introduce FLAMESFast, Low-storage, Accurate, and Memory-Efficient adaptive Sampling, which is a direct method for calculation of structure factors, allowing us to select only the required number of wavevectors for binning. We also use wavenumber-dependent time steps to extract dynamics. Our FLAMES approach effectively mitigates computational and memory/storage bottlenecks. We demonstrate the method using simulations of a model system, liquid octane, at various temperatures. Comparisons with experimental data and real space computation show that the FLAMES technique achieves high accuracy in resolving temperature- and spatially dependent dynamics while being significantly more computationally efficient and requiring less memory and storage than methods based on a uniform wavevector grid and fixed temporal spacing.
The Wigner crystal has fascinated condensed matter physicists for nearly 90 years. Signatures of two-dimensional (2D) Wigner crystals were first observed in 2D electron gases under high magnetic field, and recently reported in transition metal dichalcogenide moiré superlattices. Direct observation of the 2D Wigner crystal lattice in real space, however, has remained an outstanding challenge. Conventional scanning tunnelling microscopy (STM) has sufficient spatial resolution but induces perturbations that can potentially alter this fragile state. Here we demonstrate real-space imaging of 2D Wigner crystals in WSe 2 /WS 2 moiré heterostructures using a specially designed non-invasive STM spectroscopy technique. This employs a graphene sensing layer held close to the WSe 2 /WS 2 moiré superlattice. Local STM tunnel current into the graphene layer is modulated by the underlying Wigner crystal electron lattice in the WSe 2 /WS 2 heterostructure. Different Wigner crystal lattice configurations at fractional electron fillings of n = 1/3, 1/2 and 2/3, where n is the electron number per site, are directly visualized. The n = 1/3 and n = 2/3 Wigner crystals exhibit triangular and honeycomb lattices, respectively, to minimize nearest-neighbour occupations. The n = 1/2 state spontaneously breaks the original C3 symmetry and forms a stripe phase. Here our study lays a solid foundation for understanding Wigner crystal states in WSe 2 /WS 2 moiré heterostructures and provides an approach that is generally applicable for imaging novel correlated electron lattices in other systems.
Here, we present comprehensive phononic and charge density wave properties (CDW) of rare-earth van der Waals tritellurides through temperature dependent angle-resolved Raman spectroscopy measurements. All the possible rare-earth tritellurides (RTe3) ranging from R = La–Nd, Sm, Gd–Tm were synthesized through a chemical vapor transport technique to achieve high quality crystals with excellent CDW characteristics. Raman spectroscopy studies successfully identify the emergence of the CDW state and transition temperature (TCDW), which offers a non-destructive method to identify their CDW response with micron spatial resolution. Temperature dependent Raman measurements further correlate how the atomic mass of metal cations and the resulting chemical pressure influence its CDW properties and offer detailed insight into the strength of CDW amplitude mode-phonon coupling during the CDW transition. Angle-resolved Raman measurements offer the first insights into the CDW-phonon symmetry interplay by monitoring the change in the symmetry of phonon mode across the CDW transition. Overall results introduce the library of RTe3 CDW materials and establish their characteristics through the non-destructive angle-resolved Raman spectroscopy technique.
Liquid argon time projection chambers (LArTPCs) are a class of detectors that produce high resolution images of charged particles within their sensitive volume. In these images, the clustering of distinct particles into superstructures is of central importance to the current and future neutrino physics program. Electromagnetic (EM) activity typically exhibits spatially detached fragments of varying morphology and orientation that are challenging to efficiently assemble using traditional algorithms. Similarly, particles that are spatially removed from each other in the detector may originate from a common interaction. Graph neural networks (GNNs) were developed in recent years to find correlations between objects embedded in an arbitrary space. The graph particle aggregator (GrapPA) first leverages GNNs to predict the adjacency matrix of EM shower fragments and to identify the origin of showers, i.e., primary fragments. On the PILArNet public LArTPC simulation dataset, the algorithm achieves a shower clustering accuracy characterized by a mean purity of 99.4%, a mean efficiency of 99.6% and a primary identification accuracy of 99.8%. It yields a relative shower energy uncertainty of (4.1 + 1.4 / $\sqrt{\text{E(GeV)})}$% and a shower direction uncertainty of (2.1/ $\sqrt{\text{E(GeV)})}$°. Finally, the optimized algorithm is then applied to the related task of clustering particle instances into interactions and yields a mean purity of 99.8% and a mean efficiency of 99.5% for an interaction density of $\mathcal{O}(1)$ m –3 .
Fast neutrons enable a nondestructive examination of dense, large, and highly attenuating samples due to their lower interaction probability compared to thermal neutrons. However, this also creates a challenge in fast neutron imaging, as the thicker sensors necessary to detect fast neutrons degrade an image’s spatial resolution due to scattering within the sensor and the indeterminate depth of interaction in the sensor. This work explores the advantages of a fast neutron imaging screen with a layered polymer-phosphor screen approach as opposed to a mixed polymer-phosphor screen typically used in fast neutron imaging. Proton recoil is the primary conversion mechanism for fast neutron imaging. Simulations showed that the recoil proton range of typical fast neutrons is approximately 200 µm, however, tests at Idaho National Laboratory revealed that the light output of these screens increased at much greater polymer thicknesses. The NECTAR fast neutron beamline at FRM II was used to test the imaging performance of layered fast neutron imaging screens. Distinguishing between the fast-neutron and γ-ray signals is a major challenge in fast neutron imaging because all fast neutron sources also produce γ-rays. A relative comparison between a control plate and the fast neutron screen was made to distinguish between a γ-ray and fast neutron signals. MCNP modeling quantified the γ-ray and fast neutron contributions to the images measured at NECTAR, which were approximately a 75% γ-ray image.
Effective flood simulation capabilities can tremendously support early warning and disaster prevention. To examine the applicability of a fully physics-based and high-performance flood simulation and forecasting modeling framework for a flood-prone region in Taiwan, we conduct a numerical experiment that couples the Weather Research and Forecasting (WRF) Model, WRF-Hydrological modeling system (WRF-Hydro), and the Two-Dimensional Runoff Inundation Toolkit for Operational Needs (TRITON) to perform integrated rainfall, streamflow, and flood simulations. Furthermore, we first use the coupled WRF and WRF-Hydro (WWH) to predict rainfall and streamflow and then drive TRITON with the predicted streamflow hydrographs to simulate flood depth and inundation area. With the refined spatial resolution and parameterization, this framework can better predict rainfall with reasonable spatial patterns. Although WWH could overestimate the amount of rainfall in some areas, the uncertain rainfall–streamflow predictions produce reasonable flood maps able to pinpoint regions at risk of flooding. In terms of model efficiency, the graphics processing unit–based computation can yield a speed-up factor as high as ∼13 compared to the central processing unit–based computation, promoting the efficacy of the coupled modeling framework in practical real-time flood forecasting.
Uncovering drivers of risk is crucial to understanding interactions between wildlife and wind turbines, and identifying options for impact minimization. These drivers tie to co-variates linked to behavior and movement patterns that allow us to estimate locations and periods of risk. For volant species, atmospheric flow can have significant influence on flight patterns. For obligate soaring birds, like golden eagles, updraft velocities can inform where eagles are likely to travel, at what altitude, and where conditions are not likely sufficient to sustain soaring flight. This has been an active area of study in recent years, using relatively coarse atmospheric data generally at the 20km x 20km scale or larger. Leveraging a 20-year dataset from the Weather Research and Forecasting Model (WRF) (https://www.mmm.ucar.edu/weather-research-and-forecasting-model), we are quantifying vertical velocities across the continental United States at a 2km x 2km resolution. Wind resource data sets originally were static maps showing the mean annual wind speed over an area. However, for these data sets to be optimally used for various applications they must be high-resolution time series, seamlessly span large geographic contexts, and account for uncertainty in wind speed. The National Renewable Energy Laboratory is producing public available datasets that meet these criteria and will be bias corrected to yield the most accurate wind resource data. This effort is an augment to the current WIND Toolkit which houses a high resolution data set. In the new iteration of the WIND Toolkit, a 20-year dataset will be used to improve the accuracy and estimate uncertainty using ensemble and machine-learning techniques. The resulting product will be the most accurate dataset of its size and at a 2km x 2km spatial and 5-minute temporal resolution. Through this work, a mesoscale vertical velocity layer will be produced by calculating the likelihood of orographic updraft and thermal updraft conditions across the continental United States. Specifically, we will use WRF model output combined with digital elevation maps to predict updrafts and then determine if vertical velocities are sufficient to support Golden Eagle soaring and gliding. Ultimately this data layer will be made available as a GIS layer in the Wind Prospector (maps.nrel.gov/wind-prospector/) tool or a similar framework. Data that will be incorporated include wind speed, direction temperature, relative humidity, barometric pressure, air density, precipitation rate, solar radiation, atmospheric stability, skin temperature, and upward heat flux. These products will advance research on interactions between volant species and wind energy by providing open access to highly resolved data with uncertainty quantification not previously available at this scale.
Association of organic matter (OM) with mineral phases via co-precipitation is expected to be a widespread process in environments with high OM input and frequent mineral dissolution and re-precipitation. In contrast to surface area-limited adsorption processes, co-precipitation may allow for greater carbon (C) accumulation. However, the potential sub-micrometer scale structural and compositional differences that affect the bioavailability of co-precipitated C are largely unknown. In this work, we used a combination of high-resolution analytical electron microscopy and bulk spectroscopy to probe interactions between a mineral phase (ferrihydrite, nominally Fe 2 O 3 •0.5H 2 O) and organic soil-derived water-extractable OM (WEOM). In co-precipitated WEOM-Fe, nanometer-scale scanning transmission electron microscopy with electron energy loss spectroscopy (STEM-EELS) revealed increased Fe(II) and less Fe aggregation relative to adsorbed WEOM-Fe. Spatially distinct lower- and higher-energy C regions were detected in both adsorbed and co-precipitated WEOM-Fe. In co-precipitates, lower-energy aromatic and/or substituted aromatic C was spatially associated with reduced Fe(II), but higher-energy oxidized C was enriched at the oxidized Fe(III) interface. Therefore, we show that co-precipitation does not constitute a non-specific physical encapsulation of C that only affects Fe chemistry and spatial distribution, but may cause a bi-directional set of reactions that lead to spatial separation and transformation of both Fe and C forms. In particular, we propose that abiotic redox reactions between Fe and C via substituted aromatic groups (e.g., hydroquinones) play a role in creating distinct co-precipitate composition, with potential implications for its mineralization.
The injection of multi-species fluids into high-pressure and high-temperature environments beyond the species' critical points is commonly found in engineering applications. At these conditions, for immiscible species, both subcritical interfacial dynamics and supercritical mixing can coexist due to variations in temperature around the mixture critical point. The modeling of these complex transcritical phenomena for large-scale configurations is so far not possible. To address this issue, we propose the Regularized-Interface Method (RIM) as a unified formulation that can describe both sub- and supercritical processes as well as the transition between them. The proposed method is derived via filtering of the nanoscale interface-resolving formulation based on van der Waals' linear gradient theory. Thus, this approach allows for the consistent modeling of interfacial dynamics that vanishes at supercritical conditions, while significantly reducing the temporal and spatial resolution constraints of the original nanoscale formulation. The resulting RIM formulation is examined in interface-capturing simulations of sub-, trans-, and supercritical fuel injection processes, involving droplets and jets. Furthermore, these results highlight the importance of resolving spatio-temporal transitions from subcritical interfacial dynamics to supercritical mixing in high-pressure multiphase simulations, in contrast to commonly employed diffused-interface methods, where interfacial dynamics are often neglected.
A simple repeating step-filtered structure has been designed and deployed for the detection of MeV-energy proton beams produced from high-intensity short-pulse laser-driven experiments. Due to the diagnostic's ability to recover energy-resolved images of proton beams, it is called the proton beam imager and energy spectrometer (or PROBIES). The flexible design has two useful applications. First, it can be fielded in a configuration with a scintillator and imaging setup in order to collect proton beam information in high-repetition rate experiments. Here we show a configuration that collects the equivalent of a nine-film radiochromic film (RCF) stack, and a technique for using a neural network (NN) to rapidly analyze the images on millisecond time-scales. Second, when using RCF stacks as the detector, the energy resolution capability is greatly increased compared to traditional RCF stacks, while simultaneously reducing the total number of films required to make each measurement. We demonstrate a configuration that provides up to 90 spatially- and energy-resolved datapoints using just ten films. Both designs presented here can be easily customized for anticipated experimental outputs and manufactured at low cost. Such designs allow for rapid prototyping using modern high-accuracy 3D printers or on-demand computer numerical control machining to reduce the reliance on time-consuming film media, analysis, or bulky electro-magnetically dispersing spectrometers. The analysis of data is discussed and a methodology for constructing a NN that is capable of rapidly analyzing data on-the-fly with greater than 95% accuracy is presented.
Previous studies have documented how forests influence snow at fine spatial scales, but none have documented the influence that existing forest-snow variability has on streamflow. Here to test how much forest-controlled snow variability influences streamflow, a tiling parameterization based on classifications from high-resolution (1–3 m) vegetation maps was incorporated into the Distributed Hydrology and Soil Vegetation Model (DHSVM). Within each grid cell (90–150 m), the tiling parameterization simulated forest-snow variability with four independently evolving snowpacks. Each tile had unique radiation conditions to represent conditions underneath the canopy, in exposed areas, and along north- and south-facing forest edges. This tiled parameterization was used to test where and when detailed forest-snow modelling should be considered further and where and when the impacts are too small to be worth the effort. To test this, tiled model simulations of streamflow were compared to non-tiled model simulations in the Sierra Nevada, CA, the Jemez Mountains, NM, and the Eastern Cascades, WA. In Tuolumne, CA, the tiled model simulated little difference in grid cell average SWE, and late-season streamflow decreased by only 3%–4% compared to the non-tiled model. In Jemez, NM, the tiled model decreased late-season streamflow by 18% due to increased sublimation. In Chiwawa, WA, the tiled model increased late-season streamflow by 15% due to high shortwave radiation attenuation and less longwave radiation enhancement from the forest. Furthermore, within the Chiwawa, a substantial silvicultural practice was synthetically implemented to increase the north-facing edge's fractional area. This silvicultural experiment, which used the same fractional forest area in all simulations increased late-season streamflow by 35% compared to tiled model simulations that did not represent forest edges. In conclusion, representing forest-SWE variability had an effect on late-season streamflow in some watersheds but not in others based on the fractional area of the forest edges, forest characteristics, and climate conditions.
Tomographic imaging of time-evolving samples is a challenging yet important task for various research fields. At the nanoscale, current approaches face limitations of measurement speed or resolution due to lengthy acquisitions. We developed a dynamic nanotomography technique based on sparse dynamic imaging and 4D tomography modeling. We demonstrated the technique, using ptychographic x-ray computed tomography as its imaging modality, on resolving the in situ hydration process of polymer electrolyte fuel cell (PEFC) catalyst. The technique provides a 40-time increase in temporal resolution compared to conventional approaches, yielding 28 nm half-period spatial and 12 min temporal resolution. The results allow a quantitative characterization of the water intake process inside PEFC catalysts with nanoscale resolution, which is crucial for understanding their electrochemical mechanisms and optimizing their performance. Our technique enables high-speed operando nanotomography studies and paves the way for wider application of dynamic tomography at the nanoscale.
This manuscript describes a software suite that can be used for objective evaluation of tropical cyclones (TCs) in gridded climate data. Using cyclone trajectories derived from 6-hourly data, a comprehensive set of metrics is defined to systematically compare and contrast products to one another. In addition to annual TC climatologies, attention is paid to spatial and temporal patterns of storm occurrence and intensity. Assessment can be performed either on the global scale or regional domains. Simple to visualize ‘scorecards’ allow for rapid credibility assessment. We showcase three key findings enabled by this suite. First, we compare the representation of TCs in seven current-generation global reanalyses and conclude that higher resolution models and those with TC-specific assimilation contain more accurate storm climatologies. Second, using a free-running Earth system model (ESM) we find that full basin refinement is required in variable-resolution configurations to adequately simulate North Atlantic TC frequency. Furthermore, upstream refinement over northern Africa offers little benefit in simulating storm occurrence but spatial genesis patterns are improved. Finally, we show that TCs simulated by ESMs can be highly sensitive to individual parameterizations in climate models, with North Atlantic TC metrics varying greatly depending on version of the Morrison-Gettelman microphysics package.
Neutron reflectometers are a class of instruments that employ neutron diffraction for measuring the structure of thin films. The technique provides valuable information over a wide variety of scientific and technological applications including chemical aggregation, polymer and surfactant adsorption, structure of thin film magnetic systems, biological membranes, etc. Such instruments are considered to be one of the most important in the neutron science field. Newer facilities place growing demands on the count rate capability of the specialized neutron detectors required. Proportional Technologies completed the Phase I effort successfully with a fully operational prototype detector built at PTI and tested at the CG1A neutron beam of the High Flux Isotope Reactor (HFIR, Oak Ridge, TN). The prototype incorporated 3 panels of 14 copper-walled cells each (total of 42 cells), each lined with 1.3 μm of 10B4C. The detector was tested for efficiency, spatial resolution, count rate linearity, uniformity, overlap and scatter. Testing revealed stable operation, uniform response across cells, and successful resolution of individual cells. The maximum count rate tested in the 3 panels combined was 590 kHz, limited by a deadtime of about 1 μs in the digital electronics. The local rate measured in a single cell showed no deadtime losses up to 810 cps/mm2, indicating that a rate of 20 Mhz can be supported in the 308-cell detector planned for Phase II.
CdZnTe-based detectors are highly valued because of their high spectral resolution, which is an essential feature for nuclear medical imaging. However, this resolution is compromised when there are substantial defects in the CdZnTe crystals. In this study, we present a learning-based approach to determine the spatially dependent bulk properties and defects in semiconductor detectors. This characterization allows us to mitigate and compensate for the undesired effects caused by crystal impurities. We tested our model with computer-generated noise-free input data, where it showed excellent accuracy, achieving an average RMSE of 0.43% between the predicted and the ground truth crystal properties. In addition, a sensitivity analysis was performed to determine the effect of noisy data on the accuracy of the model.