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.
HOSS (Hybrid Optimization Software Suite) is a fracture mechanics code that properly conserves mass, energy and momentum. HOSS-FSIS was utilized to conduct exploratory calculations (2D plane strain/3D) for the G-tunnel event 2020 where a pseudo NPE 1.1kt HE equivalent source is used. The purpose of this analysis is to determine if the HOSS fluid-solid coupling solver can capture relevant down-tunnel phenomenology.
Summary The objective of this study is to understand how microseismic events are generated during hydraulic fracturing, as well as the role of geomechanical conditions (i.e., stress and mechanical stratigraphy) in this process. In the industry, microseismic event clouds have been generally used as an “outer boundary” of the “stimulated reservoir volume” (SRV). However, by comparing with other surveillance data (low-frequency distributed acoustic sensing, or LF-DAS, strain) in the Hydraulic Fracturing Test Site (HFTS) 2 experiment, we show that this assumption is fundamentally flawed. The HFTS 2 data have three unique observations that have not been commonly observed in other datasets: (1) Due to the influence of offset pad depletion, microseismic data show that hydraulic fractures from the child well can propagate over 3,000 ft into the depleted low-stress zone. (2) By comparing microseismic and horizontal fiber LF-DAS strain data, we observe that the microseismic event cloud does not necessarily reflect the created hydraulic fracture volume. Particularly, the extent of microseismic event clouds near heel stages is much shorter than what is shown with LF-DAS strain data. (3) Microseismic event magnitudes are larger in the depleted regions. Through geomechanical analysis, we demonstrate that the “bedding-plane-slip” model is likely the mechanism for microseismic generation during hydraulic fracturing. This model successfully explains the above field observations from the HFTS 2 experiment. We also provide a quantitative relationship connecting the microseismic event magnitude with fracture width increment and layer mechanical property contrast.
Characterization of induced microseismicity at a carbon dioxide (CO 2 ) storage site is critical for preserving reservoir integrity and mitigating seismic hazards. We apply a multilevel machine learning (ML) approach that combines the nonnegative matrix factorization and hidden Markov model to extract spectral representations of microseismic events and cluster them to identify seismic patterns at the Illinois Basin-Decatur Project. Unlike traditional waveform correlation methods, this approach leverages spectral characteristics of first arrivals to improve event classification and detect previously undetected planes of weakness. By integrating ML-based clustering with focal mechanism analysis, we resolve small-scale fault structures that are below the detection limits of conventional seismic imaging. Our findings reveal temporal bursts of microseismicity associated with brittle failure, providing insights into the spatio-temporal evolution of fault reactivation during CO 2 injection. This approach enhances seismic monitoring capabilities at CO 2 injection sites by improving fault characterization beyond the resolution of standard geophysical surveys.
Microlensing surveys have discovered thousands of events, with almost all events discovered within the Galactic bulge or toward the Magellanic Clouds. The Zwicky Transient Facility (ZTF), while not designed to be a microlensing campaign, is an optical time-domain survey that observes the entire northern sky, including the Galactic plane, every few nights. The ZTF observes ~10 9 stars in the g and r bands and can significantly contribute to the observed microlensing population. In this paper, we predict that the ZTF will observe ~1100 microlensing events in 3 yr of observing within 10° latitude of the Galactic plane, with ~500 events in the outer Galaxy (ℓ ≥ 10°). This yield increases to ~1400 (~800) events by combining every three ZTF exposures, ~1800 (~900) events if the ZTF observes for a total of 5 yr, and ~2400 (~1300) events for a 5 yr survey with postprocessing image stacking. Using the microlensing modeling software PopSyCLE, we compare the microlensing populations in the Galactic bulge and the outer Galaxy. We also present an analysis of the microlensing event ZTF18abhxjmj to demonstrate how to leverage these population statistics in event modeling. The ZTF will constrain Galactic structure, stellar populations, and primordial black holes through photometric microlensing.
Faceted grain boundaries can migrate in interesting and unexpected ways. For example, faceted Σ11 <110> tilt grain boundaries were observed to exhibit mobility values that could be strongly dependent on the direction of migration. In order to understand whether this directionally-anisotropic mobility is a general phenomenon and to isolate mechanistic explanations for this behavior, molecular dynamics simulations of bicrystals evolved under an artificial driving force are used to study interface migration for a range of boundary plane inclination angles and temperatures in multiple face centered cubic metals (Al, Ni, and Cu). We find that directionally anisotropic mobility is active in a large fraction of these boundaries in Ni and Cu and should therefore impact the coarsening of polycrystalline materials. On the other hand, no such anisotropy is observed in any of the Al boundaries, showing that this behavior is material-dependent. Migration of the faceted boundaries is accomplished through transformation events at facet nodes and incommensurate boundary plane facets, which are termed shuffling modes. Three major shuffling modes have been identified, namely Shockley shuffling, slip plane shuffling, and disordered shuffling. A shift from the first two ordered modes to the third disordered mode is found to be responsible for reducing or removing directionally-anisotropic mobility, especially at the highest temperatures studied.
In this paper, we conducted an analysis of the heliospheric current sheets (HCSs) in the heliosheath (HS), utilizing observations by Voyager 2 between 2008 and 2018. Employing rigorous criteria, we identified a total of 34 HCSs that indicate significant changes in magnetic polarity. These occurrences were more prevalent during solar maximum periods when the HCS expanded to higher latitudes, coinciding with the spacecraft positioned at an average latitude of –31° from the solar equator. We determined certain features of the large-scale structures of the identified HCSs. Most importantly, employing two distinct methods indicates that the thickness of the HCSs within the HS ranges from ~0.003 to ~0.4 au with an average thickness of ~0.03 to ~0.1 au, depending on methods of event selection and fitting techniques. This thickness surpasses that known near 1 au or other heliospheric distances. It is also notably thicker than the typical proton inertial length, implying unfavorable conditions for magnetic reconnection. Additionally, our analysis reveals a frequent tilt of HCS planes relative to the solar equatorial plane by a varying angle up to several tens of degrees, likely implying a common occurrence of a warped structure of the HCS within the HS. Longitudinally, the HCS planes closely align with the Parker spiral field direction expected in the HS. Finally, for a large fraction of the identified HCS events, the HCS planes are likely characterized by a rotational discontinuity. These findings are valid within the limits of the 1 hr resolution data used in this study.
NOvA is a long-baseline neutrino oscillation experiment based at Fermilab, IL, USA, that observes $\nu_\mu \ (\bar{\nu}_\mu) \to \nu_\mu\ (\bar{\nu}_\mu)$ disappearance and $\nu_\mu \ (\bar{\nu}_\mu) \to \nu_e \ (\bar{\nu}_e)$ appearance oscillations from a beam of muon neutrinos (anti-neutrinos) provided by the Fermilab's NuMI beamline. The experiment consists of two functionally identical active liquid scintillattor tracking calorimeters, both situated 14.6 mrad off-axis to the beam direction. The detectors are made out of extruded PVC cells arranged in alternating horizontal and vertical planes for 3-dimensional reconstruction of neutrino events. The NOvA experiment has a wide-ranging scientific program that includes studying standard 3-flavor neutrino oscillations, resolution of neutrino mass orderings, measuring the CP-violating phase, $\delta_{CP}$, beyond standard model (BSM) phenomenon such as non-standard interactions (NSI) and sterile neutrino searches, neutrino-nucleus cross-section measurements, exotics, astrophysics and more. The NOvA experiment leverages its identical detector technology to mitigate systematic uncertainties for its neutrino oscillation analysis. This talk will provide an overview of the experiment's detector design and the data-driven techniques used by the experiment in its neutrino oscillation analysis.
Solid-state batteries can suffer from catastrophic failure at high current densities due to solid electrolyte fracture, interface decomposition, or lithium filament growth. Failure is linked to chemomechanical material transformations that can manifest during electrochemical cycling. We systematically investigate how solid electrolyte microstructure and interfacial decomposition (e.g., interphase) affect failure mechanisms in lithium thiophosphates (Li3PS4, LPS) electrolytes. Kinetically metastable interphases are engineered with iodine doping, and microstructural control is achieved using milling and annealing processing techniques. In situ transmission electron microscopy reveals iodine diffusion to the interphase, and upon electrochemical cycling, pores are formed in the interphase region. In situ synchrotron tomography reveals that interphase pore formation drives edge fracture events, which are the origin of through-plane fracture failure. Fractures in thiophosphate electrolytes actively grow toward regions of higher porosity and are affected by heterogeneity in microstructure (e.g., porosity factor). This report provides fundamental design guidelines for high-performance solid-state batteries.
Making high-precision measurements of neutrino oscillation parameters requires an unprecedented understanding of neutrino-nucleus scattering. In this presentation, we present the first muon neutrino charged current double-differential cross sections in kinematic imbalance variables. These variables characterize the imbalance in the plane transverse to an incoming neutrino. We use events with a single muon above 100 MeV/c, a single final state proton above 300 MeV/c, and no recorded final state pions. Thus, these variables act as a direct probe of nuclear effects such as final state interactions, Fermi motion, and multi-nucleon processes. We also present a complementary ongoing analysis using electron neutrinos. This channel is of the utmost importance for the extraction of neutrino oscillation parameters by making high-precision measurements. Our measurements allow us to constrain systematic uncertainties associated with neutrino oscillation results performed by near-future experiments of the Short Baseline Neutrino (SBN) program, as well as by future large-scale experiments like DUNE.
Making high-precision measurements of neutrino oscillation parameters requires an unprecedented understanding of neutrino-nucleus scattering. In this presentation, we present the first muon neutrino charged current double-differential cross sections in kinematic imbalance variables. These variables characterize the imbalance in the plane transverse to an incoming neutrino. We use events with a single muon above 100 MeV/c, a single final state proton above 300 MeV/c, and no recorded final state pions. Thus, these variables act as a direct probe of nuclear effects such as final state interactions, Fermi motion, and multi-nucleon processes. We also present a complementary ongoing analysis using electron neutrinos. This channel is of the utmost importance for the extraction of neutrino oscillation parameters by making high-precision measurements. Our measurements allow us to constrain systematic uncertainties associated with neutrino oscillation results performed by near-future experiments of the Short Baseline Neutrino (SBN) program, as well as by future large-scale experiments like DUNE.
Making high-precision measurements of neutrino oscillation parameters requires an unprecedented understanding of neutrino-nucleus scattering. In this presentation, we present the first muon neutrino charged current double-differential cross sections in kinematic imbalance variables. These variables characterize the imbalance in the plane transverse to an incoming neutrino. We use events with a single muon above 100 MeV/c, a single final state proton above 300 MeV/c, and no recorded final state pions. Thus, these variables act as a direct probe of nuclear effects such as final state interactions, Fermi motion, and multi-nucleon processes. We also present a complementary ongoing analysis using electron neutrinos. This channel is of the utmost importance for the extraction of neutrino oscillation parameters by making high-precision measurements. Our measurements allow us to constrain systematic uncertainties associated with neutrino oscillation results performed by near-future experiments of the Short Baseline Neutrino (SBN) program, as well as by future large-scale experiments like DUNE.
Here, we present methods for analyzing Beam Emission Spectroscopy (BES) data to obtain the plasma density evolution associated with rapid sawtooth crash events at the DIII-D tokamak. BES allows coverage over a 2D spatial plane, inherently local measurements, with fast time responses, and, therefore, provides a valuable new channel for data during sawtooth events. A method is developed to remove sawtooth-induced edge-light pulses contained in the BES data. The edge light pulses appear to be from the D α emission produced by edge recycling during sawtooth events, and are large enough that traditional spectroscopic filtering and data analysis techniques are insufficient to deduce physically meaningful quantities. A cross-calibration of 64 BES channels is performed by using a novel method to ensure accurate measurements. For the large-amplitude density oscillations observed, we discuss and use the non-linear relationship between the BES signal δI/ I 0 and the plasma density variation δn e / n e0 . The 2D BES images cover an 8 × 20 cm 2 region around the sawtooth inversion layer and show large-amplitude density oscillations, with additional significant spatial variations across the inversion layer that grows and peaks near the time of the temperature crash. The edge light removal technique and method of converting large-amplitude δI/ I 0 to δn e / n e0 presented here may help analyze other impulsive MHD phenomena in tokamaks.
We use publicly available data to perform a search for correlations of high energy neutrino candidate events detected by IceCube and high-energy photons seen by the HAWC Collaboration. Our search is focused on unveiling such correlations outside of the Galactic plane. This search is sensitive to correlations in the neutrino candidate and photon skymaps which would arise from a population of unidentified point sources. We find no evidence for such a correlation, but suggest strategies for improvements with new datasets. Published by the American Physical Society 2024
ABSTRACT We present our 500 pc distance-limited study of stellar flares using the Dark Energy Camera as part of the Deeper, Wider, Faster programme. The data were collected via continuous 20-s cadence g-band imaging and we identify 19 914 sources with precise distances from Gaia DR2 within 12, ∼3 deg2, fields over a range of Galactic latitudes. An average of ∼74 min is spent on each field per visit. All light curves were accessed through a novel unsupervised machine learning techniques designed for anomaly detection. We identify 96 flare events occurring across 80 stars, the majority of which are M dwarfs. Integrated flare energies range from ∼1031–1037 erg, with a proportional relationship existing between increased flare energy with increased distance from the Galactic plane, representative of stellar age leading to declining yet more energetic flare events. In agreement with previous studies we observe an increase in flaring fraction from M0 to M6 spectral types. Furthermore, we find a decrease in the flaring fraction of stars as vertical distance from the galactic plane is increased, with a steep decline present around ∼100 pc. We find that $\sim 70{{\ \rm per\ cent}}$ of identified flares occur on short time-scales of <8 min. Finally, we present our associated flare rates, finding a volumetric rate of 2.9 ± 0.3 × 10−6 flares pc−3 h−1.
The Swift Burst Alert Telescope (BAT) is a coded-aperture gamma-ray instrument with a large field of view that primarily operates in survey mode when it is not triggering on transient events. The survey data consist of 80- channel detector plane histograms that accumulate photon counts over periods of at least 5 minutes. These histograms are processed on the ground and are used to produce the survey data set between 14 and 195 keV. Survey data comprise >90% of all BAT data by volume and allow for the tracking of long-term light curves and spectral properties of cataloged and uncataloged hard X-ray sources. Until now, the survey data set has not been used to its full potential due to the complexity associated with its analysis and the lack of easily usable pipelines. Here, we introduce the BatAnalysis Python package, a wrapper for HEASoftpy, which provides a modern, opensource pipeline to process and analyze BAT survey data. BatAnalysis allows members of the community to use BAT survey data in more advanced analyses of astrophysical sources, including pulsars, pulsar wind nebula, active galactic nuclei, and other known/unknown transient events that may be detected in the hard X-ray band. We outline the steps taken by the Python code and exemplify its usefulness and accuracy by analyzing survey data of the Crab Nebula, NGC 2992, and a previously uncataloged MAXI transient. The BatAnalysis package allows for ~18 yr of BAT survey data to be used in a systematic way to study a large variety of astrophysical sources.
The Swift Burst Alert Telescope (BAT) is a coded-aperture gamma-ray instrument with a large field of view that primarily operates in survey mode when it is not triggering on transient events. The survey data consist of 80-channel detector plane histograms that accumulate photon counts over periods of at least 5 minutes. These histograms are processed on the ground and are used to produce the survey data set between 14 and 195 keV. Survey data comprise >90% of all BAT data by volume and allow for the tracking of long-term light curves and spectral properties of cataloged and uncataloged hard X-ray sources. Until now, the survey data set has not been used to its full potential due to the complexity associated with its analysis and the lack of easily usable pipelines. Here, we introduce the BatAnalysis Python package, a wrapper for HEASoftpy, which provides a modern, open-source pipeline to process and analyze BAT survey data. BatAnalysis allows members of the community to use BAT survey data in more advanced analyses of astrophysical sources, including pulsars, pulsar wind nebula, active galactic nuclei, and other known/unknown transient events that may be detected in the hard X-ray band. We outline the steps taken by the Python code and exemplify its usefulness and accuracy by analyzing survey data of the Crab Nebula, NGC 2992, and a previously uncataloged MAXI transient. The BatAnalysis package allows for ~18 yr of BAT survey data to be used in a systematic way to study a large variety of astrophysical sources.
The High Altitude Water Cherenkov (HAWC) observatory and the High Energy Stereoscopic System (H.E.S.S.) are two leading instruments in the ground-based very-high-energy γ-ray domain. HAWC employs the water Cherenkov detection (WCD) technique, while H.E.S.S. is an array of Imaging Atmospheric Cherenkov Telescopes (IACTs). The two facilities therefore differ in multiple aspects, including their observation strategy, the size of their field of view, and their angular resolution, leading to different analysis approaches. Until now, it has been unclear if the results of observations by both types of instruments are consistent: several of the recently discovered HAWC sources have been followed up by IACTs, resulting in a confirmed detection only in a minority of cases. With this paper, we go further and try to resolve the tensions between previous results by performing a new analysis of the H.E.S.S. Galactic plane survey data, applying an analysis technique comparable between H.E.S.S. and HAWC. Events above 1 TeV are selected for both data sets, the point-spread function of H.E.S.S. is broadened to approach that of HAWC, and a similar background estimation method is used. This is the first detailed comparison of the Galactic plane observed by both instruments. H.E.S.S. can confirm the γ-ray emission of four HAWC sources among seven previously undetected by IACTs, while the three others have measured fluxes below the sensitivity of the H.E.S.S. data set. Remaining differences in the overall γ-ray flux can be explained by the systematic uncertainties. Therefore, we confirm a consistent view of the γ-ray sky between WCD and IACT techniques.
Liquid Argon Time Projection Chamber (LArTPC) technology is commonly utilized in neutrino detector designs. It enables detailed reconstruction of neutrino events with high spatial precision and low energy threshold. Its field response (FR) model describes the time-dependent electric currents induced in the anode-plane electrodes when ionization electrons drift nearby. An accurate and precise FR is a crucial input to LArTPC detector simulations and charge reconstruction. Established LArTPC designs have been based on parallel wire planes. It allows accurate and computationally economic two-dimensional (2D) FR models utilizing the translational symmetry along the direction of the wires. Recently, novel LArTPC designs utilize electrodes formed on printed circuit board (PCB) in the shape of strips with through holes. The translational symmetry is no longer a good approximation near the electrodes and a new FR calculation that employs regions with three dimensions (3D) has been developed. Extending the 2D models to 3D would be computationally expensive. Fortuitously, the nature of strips with through holes allows for a computationally economic approach based on the finite-difference method (FDM). In this paper, we present a new software package pochoir that calculates LArTPC field response for these new strip-based anode designs. This package combines 3D calculations in the volume near the electrodes with 2D far-field solutions to achieve fast and precise field response computation. We apply the resulting FR to simulate and reconstruct samples of cosmic-ray muons and 39 Ar decays from a Vertical Drift (VD) detector prototype operated at CERN. We find the difference between real and simulated data within 5%. Current state-of-the-art LArTPC software requires a 2D FR which we provide by averaging over one dimension and estimate that variations lost in this average are smaller than 7%.