Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “eruptions”

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 199 records · Page 11

Bayesian Seismoacoustic Source Location: Acoustic Approach

Seismic waves and infrasound are key technologies in the International Monitoring System (IMS) to monitor explosive events in the solid Earth and atmosphere. Energetic man-made or natural events (e.g., chemical/nuclear explosions, volcanic eruptions, and earthquakes) near the Earth’s surface produce both ground motion and atmospheric pressure disturbances which propagate as seismic waves and infrasound, respectively. Seismic waves have been generally used to detect and identify underground and near-surface events (Myers, et. al., 2007), and infrasound are sensitive to events near the surface or in the atmosphere (Modrak et. al., 2010). Due to their different sensitivities to events, they can complement to each other to improve the event detection and discrimination. The framework of joint seismoacoustic event location has recently reviewed by theoretical research (Koch and Arrowsmith, 2019). Although the early applications showed promising results to improve the accuracy of event location, their application were still limited to a small set of events selected to prove the concepts, and practical capability of the method for operational purpose is not fully evaluated with data. Our final goal is to apply the method of seismoacoustic event location to a larger set of events and evaluate its applicability for operational event location in practice. To that end, we focus on developing and verifying acoustic source location method in this study.

58 GEOSCIENCES↗

Report on Year-4 of Water NSTF Matrix Testing: Facility Maintenance and Accident Testing

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report serves as a summary of maintenance and experimental activities during the program’s fourth year of water-based operation. A planned six-month maintenance period began in August 2021, during which major inspections, repairs, cleaning, and installation of new instrumentation and data acquisition hardware were conducted. Most significantly, two heaters that faulted during Year-3 were repaired, allowing the facility to resume use of the full heated section area and full range of available electric power. The remainder of the year consisted of eight months of test operations, during which the facility logged 211 hours of active heating across one bake-out (following the maintenance period) and seven matrix tests; five classifieds as Accepted per NQA-1, one as Trending, and one as Failed. Testing began by performing two repeat cases to confirm expected facility response and behavior during both single- and two-phase flow conditions, ensuring no changes were introduced during the maintenance period that might have altered the thermal-hydraulic characteristics of the facility. In continuation of the power parametric series initiated in previous years, a high-power test case was then performed examining heat removal performance at a decay heat load equivalent to 2.4 MWt, full-scale, a level exceeding maximum design targets. Additional testing then introduced various blockage scenarios along the network piping, examining the effects of partial and complete blockages of the flow paths on the system behavior and heat removal performance. A study of static boiling tests directed at understanding the geysering two-phase instability was also conducted. The loop was filled only to the bottom of the tank outlet, creating an open loop configuration that prevents any natural circulation flow from occurring, and the heaters were powered on until the facility reached saturation conditions. Following, a series of quasi-steady-state conditions were introduced by adjusting the inventory level in the adiabatic chimney piping at decreasingly lower elevations above the heated region. A strong correlation of geysering characteristics and loop level was observed, with flow and temperature excursions decreasing in intensity, but increasing in frequency, as the fill were reduced to lower elevations along the chimney piping. Once the level fell very low in the chimney, at points near the top of the heated section, the system reached a stable state of continuous boiling without any occurrence of geysering eruptions. A final significant testing accomplishment this year was successful completion of an “accident scenario” test, whose operating conditions were based on a prototypic decay heat curve provided by Framatome and scaled for the NSTF. This test began by establishing steady-state, single-phase “normal operation” conditions, before simulating an accident trip where the availability of active cooling systems was lost. Loop temperatures gradually increased until reaching saturation and subsequent two-phase boiling flow. Over the course of an extended operational period along the defined decay heat curve, steam boil-off caused gradual but continued depletion of liquid inventory until reaching a critically low level causing flow stagnation and cessation of natural circulation heat removal. At this point, after nearly 72 hours of continuous operation, a cold refill was performed to replenish the system inventory and allow the facility to re-establish closed loop natural circulation flow and return to a safe operational state.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Innovative Subsurface Learning and Hawaiian Exploration using Advanced Tomography (ISLAND HEAT) Phase 1 (Final Report)

The Innovative Subsurface Learning and Hawaiian Exploration using Advanced Tomography (ISLAND HEAT) project aims to leverage and enhance existing Hawai'i geothermal play fairway analysis (PFA) results and validate a conceptual model-driven, optimized, least-cost exploration and geophysical suite at a proven geothermal field with known permeability, the Puna Geothermal Venture, located in the Lower East Rift Zone. The abundance of geophysical scrutinization before, during, and after the 2018 Kilauea eruption at and around the broader Puna system offers a rare opportunity to elucidate geophysical expressions within a dynamic magmatic rift setting. Development of the methodology continues with its application at a second prospective site identified by the PFA, Mauna Kea, which serves as confirmation of the integrated approach.

15 GEOTHERMAL ENERGY↗

Distributed Acoustic Sensing as a Monitoring Tool at LANL [Slides]

DAS records a variety of signals, including (1) earthquakes and volcanic eruptions, (2) explosions (DAG, mining, etc.), (3) acoustic waves, (4) ocean (waves, whales), (5) vehicle, ship and pedestrian traffic, (6) glaciology and landslides, and (7) oil and gas. Main questions include: (1) what is the extent of the signal content? (2) how does DAS compare to traditional sensors? (3) how do interrogators, fiber types, and fiber depth burials impact measurements? (4) special applications offshore to overcome the lack of instrumentation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The geology, geochemistry, and geographic distribution of southern Idaho obsidian

The Snake River Plain Volcanic Province of Southern Idaho is host to an abundance of pyroclastic obsidian deposits that constitute the principal lithic raw material utilized by Native Americans of the region throughout the Precontact Period. Archaeologists have long looked to obsidian as a source of information on past patterns of raw material selection, landscape use, and mobility through geochemical analysis of obsidian artifacts and lithic raw materials. Formation of the Snake River Plain through cataclysmic “super-eruptions” of the Yellowstone hotspot, however, presents one of the most complex and poorly understood obsidian source-scapes of North America. To clarify the geographic distribution of Southern Idaho obsidian source groups, the INL CRMO has compiled a comprehensive reference library of geologic obsidian from 138 source locales distributed across the southern half of the state. In this report, we define 28 geochemically distinct obsidian source groups that occur at these deposits through X-ray Fluorescence spectrometry and contextualize the geographic distribution of each source group relative to the history of silicic magmatism on the Snake River Plain. This research facilitates ongoing investigations at the INL CRMO into patterns of obsidian procurement, conveyance, and exchange in Southern Idaho throughout the Precontact Period.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai

The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.

54 ENVIRONMENTAL SCIENCES↗

Multiscale Wildfire Simulation Framework and Remote Sensing

Wildfire as one type of climate extreme events causes huge socioeconomic losses and damages. Large wildfires (i.e., generated pyrocumulonimbus (PyroCb)) can inject tremendous amounts of smoke into the stratosphere, where black carbon and organic carbon aerosols can persist months to years and influence climate by imposing a significant reduction in the radiative forcing like that associated with large volcanic eruptions or proposed via climate interventions such as geoengineering. Both observations and numerical modeling results clearly indicate an increasing trend in wildfire frequency and intensity in many regions during the recent decades with climate change. However, current understanding of wildfire remains largely uncertain owing to limitations of modeling capabilities in representing the multiscale wildfire physics and dynamics and a scarcity of observations constraining important wildfire and environmental variables. This study primarily aims to improve the wildfire simulation capabilities in the state-of-the-art climate model by filling in two major gaps: (1) model resolution is typically too coarse to resolve fine scale processes associated with fires, and (2) chemistry and aerosol processes in fire smoke are poorly represented.

54 ENVIRONMENTAL SCIENCES↗

Elastic Changepoint Detection for Globally-indexed Functional Time Series Data with Climate Applications

Changepoint detection is a vital tool in the application of climate data analysis. Numerous types of climate observation data are most properly represented by functional time series, implying a need for accurate changepoint detection methods applicable to functional time series data. Such data taken at a global scale often contain both spatial heterogeneity and dependence as well as phase (time) misalignment. In this report, we present methods which can detect spatially-dependent changepoints while allowing different estimates of change time and change strength depending on location. Additionally, we provide extensions to this spatially-predicted model which controls for phase variability among observations. Our methods provide the ability to detect a single change, or control for epidemic changes (where a “return-to-normal” change is more likely to be detected than the initial change). We showcase results analyzing the June 1991 eruption of Mt. Pinatubo, where our methods demonstrate the ability to accurately detect both single and epidemic changepoints even in the presence of strong seasonal variability. We find that our spatially-predicted model improves the detection of relevant changepoints versus methods which do not take spatial information into account, and we find that controlling for phase variability helps to control the false discovery rate during the detection process.

54 ENVIRONMENTAL SCIENCES↗

Interactions Between Aerosols, Meteorology, and Early Convective Cloud Lifecycle as Measured During CACTI (Final Technical Report)

The research supported by this award sought to improve the understanding and forecasting of thunderstorms. We did so by using data from the RELAMPAGO-CACTI project, which deployed a suite of instruments around a mountain range in Argentina that sees thunderstorms erupt over the same general area almost daily. The main portion of our research examined the relationship between the concentration of atmospheric particulates (dust, smoke, etc) on the intensity of thunderstorms. Despite prior studies finding that increased particulate concentration corresponded to more intense storms, we found no effect, or if anything a slightly opposite effect.

54 ENVIRONMENTAL SCIENCES↗

Computational Prediction of Infrasound Arrival Times and Directions from Stationary and Moving Impulsive Sources

This report addresses the need to predict infrasound signal arrival times and back azimuths at monitoring stations, enabling more focused and efficient searches within recorded waveform data. The primary challenge is estimating expected signal arrival windows for stationary and moving acoustic sources, such as chemical explosions, volcanic eruptions, meteoroids, and spacecraft re-entry events. To address this challenge, a reproducible methodology is described that uses simplified propagation speeds for boundary layer, tropospheric, stratospheric, and thermospheric atmospheric waveguides. While the Python source code itself is not freely available, this document provides detailed, step-by-step instructions, and equations enabling users to replicate and adapt the method independently. The method reliably predicts signal arrival intervals and back azimuths, thereby supporting rapid detection and accurate interpretation of infrasound events. Results demonstrate that this method effectively identifies plausible signal arrival intervals and directions, facilitating faster event detection and more reliable interpretation. This methodology directly supports atmospheric monitoring, planetary defense, and forensic analysis of explosive atmospheric events.

47 OTHER INSTRUMENTATION↗

The MURAVES Experiment: A Study of the Vesuvius Great Cone with Muon Radiography

The MURAVES experiment aims at the muographic imaging of the internal structure of the summit of Mt. Vesuvius, exploiting muons produced by cosmic rays. Though presently quiescent, the volcano carries a dramatic hazard in its highly populated surroundings. The challenging measurement of the rock density distribution in its summit by muography, in conjunction with data from other geophysical techniques, can help the modeling of possible eruptive dynamics. The MURAVES apparatus consists of an array of three independent and identical muon trackers, with a total sensitive area of 3 square meters. In each tracker, a sequence of 4 XY tracking planes made of plastic scintillators is complemented by a 60 cm thick lead wall inserted between the two downstream planes to improve rejection of background from low-energy muons. The apparatus is currently acquiring data. Preliminary results from the analysis of the first data sample are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Carbon stored in live ponderosa pines in the Sierra Nevada will not return to pre-drought (2012) levels during the 21st century due to bark beetle outbreaks

Outbreaks of several bark beetle species can develop rapidly in response to drought and may result in large transfers of carbon (C) stored in live trees to C stored in dead trees (10s of Tg C yr -1 in the western U.S. alone), which over time will be released back to the atmosphere. The western pine beetle (WPB) outbreak incited by the 2012–2015 mega-drought in the Sierra Nevada, California, U.S., could portend more frequent and/or severe bark beetle outbreaks as the temperature warms and drought frequency and intensity increase in the future. However, changes in the frequency and/or severity (resultant levels of host tree mortality) of beetle outbreaks are difficult to predict as outbreaks are complex with non-linear and eruptive processes primarily driven by interactions among beetle populations, the demography of hosts and other tree species, and climate and weather. Using an insect phenology and tree defense model, we projected the future likelihood of WPB outbreaks in the Sierra Nevada with climate drivers from different Earth System Models. Our goal was to understand how host (ponderosa pine, PIPO) recovery and future warming and drought affect the frequency and severity of WPB outbreaks and their C consequences. Our projections suggested that by 2100 the C stored in live PIPO (mean: 1.98 kg C m -2 , 95% CI: 1.74–2.21 kg C m -2 ) will not return to levels that occurred before the 2012–2015 drought (2012: ~2.30 kg C m –2 ) due to future WPB outbreaks. However, differences in climate models indicate a wide range of possible WPB outbreak frequencies and severities. Our results suggest that total plot basal area is the most significant factor in the mortality rate of PIPO by WPB in any given year, followed by drought severity and temperature. High levels of host basal area, higher temperature, and extreme drought all contribute to the frequency and severity of future WPB outbreaks. While PIPO basal area may decline under increased drought and warming, limiting high-stand basal area (>60 m 2 ha -1 ) may reduce the severity of future WPB outbreaks in the Sierra Nevada.

54 ENVIRONMENTAL SCIENCES↗

Life Expectancy of Evaporating Capillary Bridges Predicted by Tertiary Creep Modeling

The evaporation of capillary bridges is experimentally investigated at the microscale through a three-grain capillary cluster. This setting provides the minimum viable description of Haines jumps during evaporation, that is, capillary instabilities stemming from air entry into a saturated granular material. The displacement profile of a meniscus is obtained via digital image correlation for different grain materials, geometries, and separations. While it is well known that Haines jumps are triggered at the pore throat, we find that these instabilities are of three types depending on the separation. We also provide a temporal characterization of Haines jumps; we find that they are accurately described, as tertiary creep instabilities, by Voight’s relation, similarly to landslides and volcanic eruptions. This finding extends the description of capillary instabilities beyond their onset predicted by Laplace equilibrium. Our contribution also paves the way for a microscopically-informed description of desiccation cracks, of which Haines jumps are the precursors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quantized Information in Spectral Cyberspace

The constant-Q Gabor atom is developed for spectral power, information, and uncertainty quantification from time–frequency representations. Stable multiresolution spectral entropy algorithms are constructed with continuous wavelet and Stockwell transforms. The recommended processing and scaling method will depend on the signature of interest, the desired information, and the acceptable levels of uncertainty of signal and noise features. Selected Lamb wave signatures and information spectra from the 2022 Tonga eruption are presented as representative case studies. Resilient transformations from physical to information metrics are provided for sensor-agnostic signal processing, pattern recognition, and machine learning applications.

74 ATOMIC AND MOLECULAR PHYSICS↗

Lithium in T Coronae Borealis

T Coronae Borealis is a recurrent, symbiotic nova system currently in quiescence between its periodic ≈80 yr cycle of eruptions. Observations during inter-outburst epochs provide an opportunity to study properties of the accretion disk and the M red giant. Here we present new irradiated (blackbody veiling) models, incorporating modern molecular opacities and line lists, of spectra derived from high-resolution (22,000 ≲ R ≲ 120,000) optical echelle observations obtained at two epochs, one prior to and one post the 2015 rebrightening event at similar spectroscopic system phase. We find a lithium abundance in the secondary at both epochs to be comparable. The non-irradiated (classical) model atmospheres yield a lithium abundance, A(Li) = 1.3 ± 0.1. The irradiated model (veiled) atmospheres, which are likely a better representation of the system in which the white dwarf and accretion disk illuminate the red giant, give A(Li) = 2.4 ± 0.1.

79 ASTRONOMY AND ASTROPHYSICS↗

Near-infrared Characterization of Four Massive Stars in Transition Phases

Massive stars typically undergo short-lived post-main-sequence evolutionary phases with strong mass loss and occasional mass eruptions. Many of such massive stars in transition phases have been identified based on their dusty envelopes. The ejected material often veils the stellar photospheres so that the central stars cannot be assigned proper spectral types and evolutionary stages. The infrared spectral range has proved to be ideal for the classification of evolved massive stars and for the characterization of their environments. To improve our knowledge on the central stars of four such dust-enshrouded objects, [GKF 2010] MN 83, [GKF 2010] MN 108, [GKF 2010] MN 109, and [GKF 2010] MN 112, we collect and present their first medium-resolution K-band spectra in the 2.3–2.47 μm region and discuss the location of the stars in the JHK color–color diagram. We find that the emission-line spectra of both MN 83 and MN 112 show characteristics typically seen in luminous blue variable (LBV) stars. In addition, we propose that the presence and strength of the newly reported Mg ii lines might be used as a new complementary criterion to identify LBV candidates. The spectra of the other two objects imply that MN 108 is an O-type supergiant, whereas MN 109 could be an LBV candidate in its active phase. We derive lower limits for the reddening toward the stars and find that three of all dereddened fall into the region of confirmed LBVs.

79 ASTRONOMY AND ASTROPHYSICS↗

Laboratory Study of the Torus Instability Threshold in Solar-relevant, Line-tied Magnetic Flux Ropes

Coronal mass ejections (CMEs) occur when long-lived magnetic flux ropes (MFRs) anchored to the solar surface destabilize and erupt away from the Sun. This destabilization is often described in terms of an ideal magnetohydrodynamic instability called the torus instability. It occurs when the external magnetic field decreases sufficiently fast such that its decay index, ${n}_{}=-z\,\partial (\mathrm{ln}{B}_{})/\partial z$, is larger than a critical value, $n\gt {n}_{\mathrm{cr}}^{}$, where ${n}_{\mathrm{cr}}^{}=1.5$ for a full, large aspect ratio torus. However, when this is applied to solar MFRs, a range of conflicting values for ${n}_{\mathrm{cr}}^{}$ is found in the literature. To investigate this discrepancy, we have conducted laboratory experiments on arched, line-tied flux ropes and applied a theoretical model of the torus instability. Our model describes an MFR as a partial torus with foot points anchored in a conducting surface and numerically calculates various magnetic forces on it. This calculation yields better predictions of ${n}_{\mathrm{cr}}^{}$ that take into account the specific parameters of the MFR. Here, we describe a systematic methodology to properly translate laboratory results to their solar counterparts, provided that the MFRs have a sufficiently small edge safety factor or, equivalently, a large enough twist. After this translation, our model predicts that ${n}_{\mathrm{cr}}^{}$ in solar conditions falls near ${n}_{\mathrm{cr}}^{\mathrm{solar}}\sim 0.9$ and within a larger range of ${n}_{\mathrm{cr}}^{\mathrm{solar}}\sim (0.7,1.2)$, depending on the parameters. The methodology of translating laboratory MFRs to their solar counterparts enables quantitative investigations of CME initiation through laboratory experiments. These experiments allow for new physics insights that are required for better predictions of space weather events but are difficult to obtain otherwise.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiple Sources of Solar High-energy Protons

During the 24th solar cycle, the Fermi Large Area Telescope (LAT) has observed a total of 27 solar flares possessing delayed γ -ray emission, including the exceptionally well-observed flare and coronal mass ejection (CME) on 2017 September 10. Based on the Fermi/LAT data, we plot, for the first time, maps of possible sources of the delayed >100 MeV γ -ray emission of the 2017 September 10 event. The long-lasting γ -ray emission is localized under the CME core. The γ -ray spectrum exhibits intermittent changes in time, implying that more than one source of high-energy protons was formed during the flare–CME eruption. We find a good statistical correlation between the γ -ray fluences of the Fermi/LAT-observed delayed events and the products of corresponding CME speed and the square root of the soft X-ray flare magnitude. Data support the idea that both flares and CMEs jointly contribute to the production of subrelativistic and relativistic protons near the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗