Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Injection data analysis”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Intense Equatorial Electrojet and Counter Electrojet caused by the 15 January 2022 Tonga Volcanic Eruption: Space and Ground-based Observations

Abstract Text: We present space and ground-based multi-instrument observations demonstrating the impact of the 2022 Tonga volcanic eruption on dayside equatorial electrodynamics. A strong counter electrojet (CEJ) was observed by Swarm and ground-based magnetometers on 15 January after the Tonga eruption and during the recovery phase of a moderate geomagnetic storm. Swarm also observed an enhanced equatorial electrojet (EEJ) preceding the CEJ in the previous orbit. The observed EEJ and CEJ exhibited complex spatiotemporal variations. We combine them with the Ionospheric Connection Explorer (ICON) neutral wind measurements to disentangle the potential mechanisms. Our analysis indicates that the geomagnetic storm had minimal impact; instead, a large-scale atmospheric disturbance propagating eastward from the Tonga eruption site was the most likely driver for the observed intensification and directional reversal of the equatorial electrojet. The CEJ was associated with strong eastward zonal winds in the E-region ionosphere, as a direct response to the lower atmosphere forcing. Plain-Language Summary: The Earth's E-region ionosphere (~100-150 km altitude) consists of both ionized and neutral gasses, and the two components are coupled through ion-neutral collisions. The state of this region is closely influenced by neutral atmospheric activities from the lower atmosphere and the variability of the solar drivers. On 15 January 2022, the Tonga volcano had a massive eruption and injected an enormous amount of mass and energy into the atmosphere causing disturbances in the E-region ionosphere or even higher. There was also a moderate geomagnetic storm that started one day before the eruption and ended days after. These conditions offer a unique opportunity to understand the different roles they play in controlling the ionosphere. Coordinated observations including the atmosphere, ionosphere and magnetosphere were made from both space and on the ground during this event. We analyzed the magnetic field and neutral wind data and found that a large-scale atmospheric disturbance generated by the volcano eruption was responsible for the observed directional reversal of the dayside equatorial electric field and electric current.

Guan Le↗
DOE OSTI2020

Ultrafast Relaxations in Ruthenium Polypyridyl Chromophores Determined by Stochastic Kinetics Simulations

Maximizing the efficiency of solar energy conversion using dye assemblies rests on understanding where the energy goes following absorption. Transient spectroscopies in solution are useful for this purpose, and the time-resolved data are usually analyzed with a sum of exponentials. This treatment assumes that dynamic events are well separated in time, and that the resulting exponential prefactors and phenomenological lifetimes are related directly to primary physical values. Such assumptions break down for coincident absorption, emission, and excited state relaxation that occur in transient absorption and photoluminescence of tris(2,2'-bipyridine)ruthenium(2+) derivatives, confounding the physical meaning of the reported lifetimes. Here, we use inductive modeling and stochastic chemical kinetics to develop a detailed description of the primary ultrafast photophysics in transient spectroscopies of a series of Ru dyes, as an alternative to sums of exponential analysis. Commonly invoked three-level schemes involving absorption, intersystem crossing (ISC), and slow nonradiative relaxation and incoherent emission to the ground state cannot reproduce the experimentally measured spectra. The kinetics simulations reveal that ultrafast decay from the singlet excited state manifold to the ground state competes with ISC to the triplet excited state, whose efficiency was determined to be less than unity. The populations predicted by the simulations are used to estimate the magnitudes of transition dipoles for excited state excitations and evaluate the influence of specific ligands. The mechanistic framework and methodology presented here are entirely general, applicable to other dye classes, and can be extended to include charge injection by molecules bound to semiconductor surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

Radionuclide measurements by accelerator mass spectrometry at Arizona

Over the past years, Tandem Accelerator Mass Spectrometry (TAMS) has become established as an important method for radionuclide analysis. In the Arizona system the accelerator is operated at a thermal voltage of 1.8MV for C-14 analysis, and 1.6 to 2MV for Be-10. Samples are inserted into a cesium sputter ion source in solid form. Negative ions sputtered from the target are accelerated to about 25kV, and the injection magnet selects ions of a particular mass. Ions of the 3+ charge state, having an energy of about 9MeV are selected by an electrostatic deflector, surviving ions pass through two magnets, where only ions of the desired mass-energy product are selected. The final detector is a combination ionization chamber to measure energy loss (and hence, Z), and a silicon surface-barrier detector which measures residual energy. After counting the trace iosotope for a fixed time, the injected ions are switched to the major isotope used for normalization. These ions are deflected into a Faraday cup after the first high-energy magnet. Repeated measurements of the isotope ratio of both sample and standards results in a measurement of the concentration of the radionuclide. Recent improvements in sample preparation for C-14 make preparation of high-beam current graphite targets directly from CO2 feasible. Except for some measurements of standards and backgrounds for Be-10 measurements to date have been on C-14. Although most results have been in archaeology and quaternary geology, studies have been expanded to include cosmogenic C-14 in meteorites. The data obtained so far tend to confirm the antiquity of Antarctic meteorites from the Allan Hills site. Data on three samples of Yamato meteorites gave terrestrial ages of between about 3 and 22 thousand years.

Jull, A. J. T.↗
DOE OSTI2024

Advancing spatiotemporal forecasts of CO 2 plume migration using deep learning networks with transfer learning and interpretation analysis

Accurate and timely forecasts of CO 2 plume distribution throughout the injection and post-injection phases are crucial for detecting plume migration, assessing leakage risks, and supporting operational decisions in geologic carbon storage (GCS). Current convolutional neural network-based approaches primarily focus on spatial information and overlook temporal dependencies in plume distributions, thus limiting their ability to capture dynamic movement effects and provide accurate predictions of plume migration. In this work, we propose two deep learning models, Auto-Encoder (AE)-LSTM and Encoder-Decoder (ED)-ConvLSTM, each uniquely designed to capture both spatial and temporal features. We apply the proposed methods to forecast the dynamic distribution of CO 2 plumes based on 108 reservoir simulations over a 30-year injection and a 30-year post-injection period. The results indicate that the ED-ConvLSTM model outperforms the AE-LSTM model in accurately predicting the spatiotemporal dynamics of CO 2 plume migration, achieving R 2 values above 0.99. To provide a deeper understanding of these model predictions, we employ a gradient-based explanation method on the trained models. This approach provides insights into the influence of input variables on plume migration forecasts and uncovers the underlying prediction mechanisms of the proposed models. Furthermore, we introduce a transfer learning technique, enabling fast and accurate plume migration forecasting in the post-injection phase by leveraging the trained model during the injection phase. This reduces the necessity for extensive data collection or re-training. In conclusion, the methods proposed in our work enhances the performance and interpretability of CO 2 plume migration forecasts, thereby facilitating informed decision-making throughout the entire lifecycle of GCS applications.

58 GEOSCIENCES↗
DOE OSTI2025

Control Mechanisms for Self‐Sealing in Activated Clay‐Rich Faults Through Controlled Hydraulic Injection Experiment

Abstract In a high‐pressure injection fault activation experiment conducted at the Mont Terri underground research laboratory in Switzerland, the transmissivity of the Opalinus Clay fault significantly increased due to opening and shearing. The fluid injection, spanning a few hours, generated a 10 m radius fault activation patch. Subsequent pressure pulse tests conducted bi‐weekly for a year revealed the gradual return of fault transmissivity to its initial state. The study utilized fluid pressure decay analysis, optical fiber monitoring, continuous active source seismic measurements and borehole displacement sensors for measuring fault displacements. The fault zone exhibited a dilation of approximately 1.4 mm, associated with both normal and tangential movements during activation, resulting in a sudden transmissivity increase from 1 × 10 −12 to 3.2 × 10 −7 m 2 /s. Early post‐activation, transient compaction and the subsequent slow compaction were observed, transitioning to an extension regime. The pressure pulse tests demonstrated a rapid transmissivity drop by more than two orders of magnitude within the first 10 days, followed by a gradual and less pronounced decrease. Plastic shear and compaction dominated the transmissivity evolution until 70 days after injection ended, followed by a period where additional factors, such as clay mineral swelling, influenced the behavior. Extrapolation suggested a sealing process taking at least 50 years after the initial activation. Plain Language Summary A field‐scale fault activation experiment offers valuable insights into the elasto‐plastic processes governing the sealing of shale faults. The experiment reveals a rapid increase in the fault's transmissivity by approximately five orders of magnitude during activation. Subsequent observations show a gradual transmissivity decrease by about three orders of magnitude post‐activation, with slow long‐term plastic shear and compaction of the fault competing against secondary processes, notably clay mineral swelling. All conceptual models employed to interpret these field data converge on the estimation that the fault's return to its initial low transmissivity state would require a minimum of 50 years. Key Points High‐pressure injection fault activation experiment at the Mont Terri underground research laboratory Continuous transmissivity measurements record self‐sealing inside a clay‐rich fault zone Transmissivity undergoes a phase of domination by slow plastic compaction and shearing during the initial post‐activation period, with mineral swelling exerting its influence over the long term

Guglielmi, Yves↗
NASA NTRS1996

Grant Proposal for the Continuation of the Voyager Interstellar Mission: LECP Investigation

This proposal documents the plans of the Low Energy Charged Particle (LECP) investigation team for participation in NASA's Voyager Interstellar Mission (VIM) as the Voyager 1 and 2 spacecraft explore the outer reaches of the heliosphere and search for the termination shock and the heliopause. The proposal covers the four year period from 1 January 1997 to 31 December 2000. The LECP instruments on Voyager 1 and 2 measure in situ intensities of charged particles with energies from about 30 keV to 100 MeV for ions, and about 20 keV to greater than 10 MeV for electrons. The instruments provide detailed spectral, angular, and compositional information about the particles. Composition is available for greater than 200 keV/nuc using multi-parameter measurements. Angular information is obtained by a mechanically scanned platform that rotates at various commanded rates. Measurements of low energy ion and electron intensities versus time and spatial location within the heliosphere contain an abundance of information regarding various transport and acceleration processes on both local (approx. 1 hr, approx. 0.01 AU) and global (approx. 11 yrs, approx. 100 AU) scales. The LECP instruments provide unique observations of such dynamical processes, and we anticipate that it will return critical information regarding the boundaries of the heliosphere. Several recent and exciting discoveries based on LECP measurements emphasize the important role that low energy charged particle distributions play in physical processes in the interplanetary medium. Yet, at the same time, these discoveries also underscore the fact that our understanding of processes in the outer heliosphere is, in most cases, incomplete, and in others, only rudimentary at best. Among the discoveries referred to above are the following: (1) Shocks: Examination of greater than 30 keV ion intensities have revealed: (a) a total absence of acceleration beyond only -100-200 keV at a strong transient shock in May 1991 at 35 AU, despite an enhanced level of seed particles; (b) a large transient shock in September 1991 of global scale, with intensities of shock-accelerated ions greater than or equal to 30 keV to approx. 30 MeV showing complex, highly energy-dependent spatial evolution, and small-scale (approx. few gyroradii), often anisotropic, micro-structures; (c) recurrent intensity increases in greater than or equal to 30 keV to -few MeV ions, with structures that, in some cases, show no correlation with the associated corotating shock. (2) Superthermal ion pressure: A global merged interaction region with a leading shock, downstream of which the superthermal ion (greater than or equal to 30 keV to approx. 4 MeV) pressure is comparable to that of the thermal plasma, and the total particle pressure yields a plasma beta of order unity. (3) Pickup ions: Measurements of the C/O ratio within transient structures at 35-45 AU showing the first clear evidence that transient shocks can pre-accelerate interstellar pickup ions from approx. 1 keV/nuc to at least 1 MeV/nuc. (4) Seed particles: Injection of ions for acceleration to high energies at the termination shock is unlikely to be a problem, since interplanetary transient and recurrent shocks are continually accelerating ions, of solar wind or interstellar origin, to highly superthermal energies. (5) Precursor electrons: Ambient solar electrons (greater than or equal to few tens of keV) that exist in the outer heliosphere ca form a broad precursor, several days wide, that is upstream of the termination shock and potentially observable a few months prior to the shock crossing. (6) Solar wind velocity at Voyager 1: We can use LECP ion data to obtain the solar wind velocity at Voyager 1, enabling us to provide critical measurement of the plasma flow as we approach and encounter the termination shock and other regions (necessary due to the partial failure of the Voyager 1 PLS experiment). The work of the LECP investigator team during the VIM will include: (1) Continuing operations with regard to the receipt, processing, verification, cataloging, display, and distribution of the data from the LECP instruments on Voyager 1 and 2, (2) Monitoring the health and performance of the LECP instruments, and evaluating and characterizing the response of the LECP instruments to various energetic particle and plasma environments, (3) Participating in, and supporting Voyager Project planning exercises and other coordinated activities relevant to exploration of the outer heliosphere, (4) Developing analysis techniques and operational procedures suitable for searching for and characterizing the boundaries and unique regions of the outher heliosphere, (5) Continuing the preparation of data sets appropriate for submission to the National Space Sciences Data Center (NSSDC) and, where appropriate, the Planetary Data System (PDS), (6) Maintaining direct Web access to online LECP data through the JHU/APL Voyager LECP home page, (7) Performing scientific evaluations of the Voyager 1 and 2 LECP data sets in conjunction with other data sets and other investigators, with particular focus on the outer regions of the heliosphere, and (8) Publishing the results of these evaluations in the scientific literature and presenting the results in scientific conferences.

Krimigis, Stamatios M.↗
DOE OSTI2024

Improving Robustness of Spectrogram Classifiers with Neural Stochastic Differential Equations

Signal analysis and classification is fraught with high levels of noise and perturbation. Computer-vision-based deep learning models applied to spectrograms have proven useful in the field of signal classification and detection; however, these methods aren't designed to handle the low signal-to-noise ratios inherent within non-vision signal processing tasks. While they are powerful, they are currently not the method of choice in the inherently noisy and dynamic critical infrastructure domain, such as smart-grid sensing, anomaly detection, and non-intrusive load monitoring. Currently, these models can be brittle, which makes them susceptible to noisy input. This also means they have sub-optimal stability of explanation outputs. Experts and technicians using these models to make decisions in real world scenarios need assurance that a model is performing as it is supposed to. The classification or prediction outputs it generates should be sound and grounded, not likely to change in the presence of shifting noise landscapes. In this work, we explore the idea of Neural Stochastic Differential Equations (NSDE's) to improve the robustness of models trained to classify time series data and the effect of NSDE's on the explainability of outputs. We then test the effectiveness of these approaches by applying them to a non-intrusive load monitoring (NILM) dataset that consists of simulated harmonic signals injected into a real building.

Brogan, Joel↗
NASA NTRS1996

Application of the kappa-omega Turbulence Model to Quasi-Three-Dimensional Turbomachinery Flows

Many computational fluid dynamics codes for turbomachinery use the Baldwin-Lomax (B-L) turbulence model. It is easy to implement in two dimensions and works well for predicting overall turbomachinery performance. However, it is awkward to implement in three dimensions, often has difficulty finding the length scale, has a crude transition model, and neglects freestream turbulence, surface roughness, and mass injection. The kappa-omega model developed by Wilcox is an appealing alternative for several reasons. First, it is the only two-equation model that can be integrated to the wall without requiring damping functions or the distance to the wall, and hence, should behave well numerically. Second, the effects of freestream turbulence, surface roughness, and mass injection are easily included in the model. Finally, transition can be simulated using the low Reynolds number version of the model. Menter applied the kappa-model to external flows and showed very good results for flows with adverse pressure gradients. Liu and Zheng described their implementation of the kappa-model in a cascade code that included an area change term to account for endwall convergence. They validated the model for a flat plate, and compared the B-L and kappa-models to measured surface pressures for a low-pressure turbine cascade. Since they did not use the low Reynolds number version of the model, their results showed problems resulting from early transition. In this Note the low Reynolds number kappa-model was incorporated in the author's quasi-three-dimensional turbomachinery analysis code. The code includes the effects of rotation, radius change, and stream-surface thickness variation, and also includes the B-L turbulence model. The kappa-omega model was implemented using many of Menter's recommendations and an implicit approximate-factorization scheme described by Baldwin and Barth. The model was tested for a transonic compressor with rotation and variable stream-surface radius and height, and for a transonic turbine vane with transition and heat transfer. Results were compared to the B-L model and to experimental data.

Chima, Rodrick V.↗
DOE OSTI2025

Benchmark Exercise Report for Experimental Study of Bubble Scrubbing in Sodium Pool

Mechanistic source term (MST) analyses are likely to be an important part of advanced reactor licensing applications. For the purpose, an MST analysis code applicable to newly introduced advanced reactors, called SRT (Simplified Radionuclide Transport) code, has been developed by Argonne National Laboratory. SRT can track overall behaviors of radionuclides especially in metal fuel-based sodium fast reactors (SFRs) and microreactors. Throughout the simulation, migration inside fuel pins before failure, interaction with coolant (for SFR), removal/leakage in cover gas and containment (or confinement), and environmental dose impacts are considered alongside radioactive decay for short-lived nuclides. Among the postulated process, pool scrubbing phenomenon, especially under sodium pool condition, has been identified as high importance with limited supportive data. The phenomenon plays a crucial role in assessing the degree of radiological impacts as radioactive aerosols or vapors are efficiently and effectively removed during the process. To provide validation basis for SRT in assessing pool scrubbing performance inside sodium pools, the University of Wisconsin-Madison performed tests including extensive parametric effects. Separate effect tests were conducted to directly evaluate the SRT models and to estimate degree of contribution by each contributing factor. Specifically, bubble size, aerosol size, aerosol density, aerosol concentration, pool depth, system temperature, and bubble swarm effects were considered. According to the parametric effects, decontamination performance enhances with decreasing bubble size, large density, and deeper pool height. Aerosol concentration provides no effect for the whole range of interest, and pool temperature variation shows minor effects under the considered temperature condition. When multiple bubbles are injected generating a bubble swarm condition, DF performance further enhances by bubble interactions and turbulence characteristics. The measurement shows the exceptional importance of aerosol size range considered, with the lowest decontamination, where most radionuclides are expected to escape.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Review of Design Elements within Power Infrastructure Cyber–Physical Test Beds as Threat Analysis Environments

Cyber–physical systems (CPSs) are an integral part of modern society; thus, enhancing these systems’ reliability and resilience is paramount. Cyber–physical testbeds (CPTs) are a safe way to test and explore the interplay between the cyber and physical domains and to cost-effectively enhance the reliability and resilience of CPSs. Here a review of CPT elements, broken down into physical components (simulators, emulators, and physical hardware), soft components (communication protocols, network timing protocols), and user interfaces (visualization-dashboard design considerations) is presented. Various methods used to validate CPS performance are reviewed and evaluated for potential applications in CPT performance validation. Last, initial simulated results for a CPT design, based on the IEEE 33 bus system, are presented, along with a brief discussion on how model-based testing and fault–injection-based testing (using scaling and ramp-type attacks) may be used to help validate CPT performance.

14 SOLAR ENERGY↗
DOE OSTI2023

Assessment of Sodium Thermal Stratification Models Utilizing the TSTF Benchmark

As a result of certain transient scenarios, a thermally stratified layer of liquid sodium can develop in the bulk coolant volumes of a sodium-cooled fast reactor (SFR). In addition to the effects a stratification layer has on the temperature of the heat transport system, a stratification layer can also influence the transition to and establishment of natural circulation flow, which plays an important role in passive cooling and the inherent safety of a pool-type SFR. Therefore, the ability to accurately capture thermal stratification phenomena is important when demonstrating the safety basis of a pool-type SFR during transient sequences. The present work assesses various computational models with different fidelities in their ability to predict thermal stratification in the upper plenum of an SFR. Each computational model will be assessed using the data generated at the Thermal Stratification Test Facility (TSTF) located at the University of Wisconsin-Madison. Using measured flow rate and inlet temperature data, the measured temperature distributions of the tests are compared to the predictions of the lumped volume-based models in SAS4A/SASSYS-1, a 1D-based model in SAM, and a 3-D computational fluid dynamics (CFD) model using STAR-CCM+. The relative performance of the various computational methods is assessed with respect to key metrics such as bulk coolant temperature distribution and plenum exit temperature. A total of eight tests are analyzed, covering different combinations of flow rates (3 and 10 GPM) and upper internal structure (UIS) configurations (none, solid, porous, and open) The perfect mixing model of SAS4A/SASSYS-1 provides the highest accuracy when the flow rate is high and there is no UIS in the test vessel, as high flow rate injection promotes thermal mixing of the sodium in the test vessel. For most of the analyzed tests, the stratified volume model of SAS4A/SASSYS-1 is able to predict the delay in the outlet temperature drop and temperature distribution in the test vessel by a small number of layers to represent thermal stratification. However, the stratified volume model can only simulate a maximum of three temperature layers within a volume and when a layer approaches the elevation of the outlet, the predicted outlet temperature can demonstrate rapid, non-physical changes. The 1-D axial mixing model of SAM provides results that agree reasonably well with the measured data in the prediction of the temporal evolution of the outlet temperature with the exception of the case with a high flow rate and no UIS. The SAM 1-D model has a similar level of accuracy to CFD results when it comes to predicting the outlet temperature. CFD shows overall good agreement in predicting the temperature distribution in the test vessel and outlet temperature. As CFD can model the test vessel geometry in detail, it performs well in the cases of complex geometries such as tests that included a UIS and internal flow through the UIS resulting in active mixing of the coolant in the test vessel. Each of the models discussed in the present work has the potential to be useful during the various stages of reactor design, analysis, and licensing. The lumped-volume approach can be applied for fast turnaround safety calculations to obtain overall reactor behavior during transients. The 1-D models provide improved accuracy when stratification is expected for a relatively low increase in the computational cost. The CFD model can be utilized for confirmatory analysis of the 1-D model, when experimental measurements are not available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

L-H transition trigger physics in ITER-similar plasmas with applied n = 3 magnetic perturbations

The L–H transition power threshold PLH is observed to increase with applied n = 3 resonant magnetic perturbations (RMP) in ITER-similar-shape plasmas with balanced neutral beam torque injection in DIII-D. The increase is most pronounced with added electron–cyclotron heating: PLH increases with decreasing edge plasma collisionality as PLH/PLH-08 ~ (ν*)−0.5, where PLH-08 is the 2008 ITPA multi-machine power threshold scaling. This result raises concerns for H-mode access at low edge collisionality in ITER, where RMP may have to be applied before the L–H transition to safely suppress the first edge-localized mode. Non-axisymmetric modifications with RMP include a simultaneous reduction of the radial electric field (Er) well depth and E × B shear. This can be attributed to increasing edge toroidal co-current rotation, and is consistent with substantially increased local long-wavelength turbulence (measured via beam emission spectroscopy). At high RMP perturbation strength the edge electric field Er reverses sign locally (becomes positive), with changes in dominant turbulence modes. Edge magnetic stochasticity provides an attractive explanation of the observed modifications, and the observed changes in toroidal rotation and Er are consistent with a simple fluid model describing radial electron current flow along stochastic fieldlines. The observed collisionality dependence of the L-mode edge electric field with applied RMP is also qualitatively consistent with this model. Reflectometry data indicate a significant reduction of the normalized L-mode radial density gradient a/Ln at high RMP field with simultaneous increase in radial particle flux and electron thermal flux from power balance analysis. We conjecture that the increase of PLH with RMP results from the combined effects of reduced E × B flow shear (increasing turbulent transport levels) and toroidal/poloidal flow modulation due to edge stochasticity. Initial experiments indicate that non-resonant n = 3 magnetic perturbations lead only to relatively small changes in Er, E × B shear and fluctuation characteristics, and have less impact on the L–H transition power threshold. This motivates further exploration of the RMP spectrum dependence of PLH for possible mitigation of the observed threshold increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Finite Element Analysis (FEA) for Water-Foam Fracturing of Granite Rock

In addition to the foam data that were obtained from literature and that were collected from the current study, simulation data was also generated from finite element analysis (FEA) conducted in this study using COMSOL Multiphysics software. The FEA models were built to simulate the experiments conducted at Oak Ridge National Laboratory (ORNL) on cement and granite samples. In these FEA models, temperature was kept at ambient while the pressure profile resembled the loading conditions during the ORNL experiments, where pressure was either monotonically increased or applied cyclically. The cement material was used as a model material and was used to study Von Mises stress and tensile stress distribution for different bore hole length geometry using a parametric sweep with water as fracturing fluid using solid-fluid interaction module. For the granite material, FEA models were developed for stress analysis of cylindrical samples with water or foam fluids. The solid mechanics module in COMSOL was implemented to solve for Von Mises stress and tensile stress. The fluid-structure interaction module was implemented to solve for water-foam interaction on granite cylinder with addition of fluid-loading on structure, i.e., large deformation in solid mechanics with no impact on fluid deformation. Foam was considered as a pseudo single-phase compressible fluid for which material properties were calculated from water and gas (nitrogen) phases. The density of foam is calculated as a function of the densities of water and nitrogen, while viscosity is a function of temperature. Four types of FEA analyses were modelled: 1. Monotonic injection with water 2. Monotonic injection with foam 3. Cyclic injection with water 4. Cyclic injection with foam All the COMSOL files are converted to a zip file which is save in .mph.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗
DOE OSTI2020

Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 1—Theory

Summary Core-scale measurements are considered the ground truth that oil and gas or electric utility companies use to predict the migration of fluids such as oil, natural gas, carbon dioxide, or brine deep underground during their extraction or injection operations. To provide a greater understanding of petrophysical properties of low-permeability geologic formations such as shales and tight gas sandstones, this study introduces a novel core-analysis procedure. The technique follows conventional pressure-pulse-decay permeametry, where the pressure in an inlet chamber adjacent to a cylindrical core plug undergoes a rapid pressurization, the system is shut in, and the pressure reaches a new equilibrium. However, unlike a standard unidirectional pressure-pulse decay, the full-immersion pressure-pulse decay (Hannon 2019) applies a pressure disturbance to the entire outer surface area of the sample. This article covers the numerical simulator designed to model flow through an anisotropic porous sample in this scenario. The model assumes distinct but uniform permeabilities along the radial and axial directions of the cylindrical plug sample. When extracting a plug vertically (or perpendicular to bedding), the permeability along the radial direction associates with the horizontal permeability (i.e., parallel to bedding), whereas in the axial direction, flow occurs perpendicular to bedding (a vertical permeability). Investigations of these model outputs demonstrate an approximately 20-fold decrease in time to complete a full-immersion experiment compared with conventional pressure-pulse decay. Furthermore, the pressure-decay curves resulting from the full-immersion method have slightly different shapes than those resulting from other unidimensional transient methods. These differences begin to demonstrate that, under achievable experimental conditions, the analysis of pressure data from one full-immersion test could enable the simultaneous estimation of the apparent permeabilities parallel and perpendicular to bedding of a cylindrical sample in addition to its porosity. A follow-up article finalizes the proof of this capability with a parameter-estimation procedure and presents experimental verification through a proof-of-concept study. Described in greater detail in that article, the parameter estimator requires multiple forward simulations to analyze the data, which behooves minimizing the compute time for each simulation. By using an alternating direction implicit time-marching scheme and a structured but variably spaced grid, the numerical simulator built for this purpose provides a forward model output with suitable accuracy in approximately 0.5 seconds.

Energy & Fuels↗
NASA NTRS1994

Simultaneous observations of VLF ground transmitter signals on the DE 1 and COSMOC 1809 satellites: Detection of a magnetospheric caustic and a duct

Khabarovsk transmitter signals (15.0 kHz, 48 deg N, 135 deg E) were observed on the high-altitude (approximately 15000 km) Dynamic Explorer 1 (DE 1) and the low-altitude (approximately 960) km COSMOS 1809 satellites during a 9-day period in August 1989. On 7 out of 9 days the linear wave receiver (LWR) on the DE 1 satellite also detected signals from the Alpha transmitter (11.9-15.6 kHz) in Russia and an Omega transmitter (10.2-13.6 kHz) in Australia, as well as natural VLF emissions such as hiss, chorus, whistlers, and wideband impulsive signals. On two days, August 23 and 27, 1989, observations of the Khabarovsk transmitter signals were simultaneously carried out at high altitude on thre DE 1 satellite and at low altitude on the COSMOS 1809 satellite. Analysis of data from these 2 days has led to several new results on the propagation of whistler mode signals in the Earth's magnetosphere. New evidence was found of previously reported propagation phenomena, such as (1) confinement of transmitter signals in the conjugate hemisphere at ionospheric heights (approximately 1000 km), (2) observation of direct multipath propagation on both DE 1 and COSMOS 1809, (3) detection of ionospheric irregularities of is less than or equal to 100 km scale size with a few percent enhancement in electron density. We also report the first direct detection of a magnetospheric dust at L = 2.94 which was believed to be responsible for the ducted propagation of Khabarovsk signals observed in the COSMOS 1809 satellite. Our study also indicates that duct end points can extend down to at least approximately 1000 km. The peak electric and magnetic fields of ducted Khabarovsk transmitter signals at approximately 1000 km were 520 micro V/m and 36 pT respectively. Estimated field strengths of these signals inside the dust at the geomagnetic equator were 57 micro V/m and 12 pT for electric and magnetic field respectively. The results of two-dimensional ray tracing simulations were consistent with the observations of the nonducted whistler-mode propagation of Khabarovsk (15 kHz) and Alpha (11.9 kHz) signals from the transmitter location to the DE 1 and COSMOS 1809 satellites. Our results have direct implications for the question of accessibility of waves injected from the ground to various regions of the ionosphere and the magnetosphere.

Chmyrev, V. M.↗
NASA NTRS2010

Altair Descent and Ascent Reference Trajectory Design and Initial Dispersion Analyses

The Altair Lunar Lander is the linchpin in the Constellation Program (CxP) for human return to the Moon. Altair is delivered to low Earth orbit (LEO) by the Ares V heavy lift launch vehicle, and after subsequent docking with Orion in LEO, the Altair/Orion stack is delivered through translunar injection (TLI). The Altair/Orion stack separating from the Earth departure stage (EDS) shortly after TLI and continues the flight to the Moon as a single stack. Altair performs the lunar orbit insertion (LOI) maneuver, targeting a 100-km circular orbit. This orbit will be a polar orbit for missions landing near the lunar South Pole. After spending nearly 24 hours in low lunar orbit (LLO), the lander undocks from Orion and performs a series of small maneuvers to set up for descending to the lunar surface. This descent begins with a small deorbit insertion (DOI) maneuver, putting the lander on an orbit that has a perilune of 15.24 km (50,000 ft), the altitude where the actual powered descent initiation (PDI) commences. At liftoff from Earth, Altair has a mass of 45 metric tons (mt). However after LOI (without Orion attached), the lander mass is slightly less than 33 mt at PDI. The lander currently has a single descent module main engine, with TBD lb(sub f) thrust (TBD N), providing a thrust-to-weight ratio of approximately TBD Earth g's at PDI. LDAC-3 (Lander design and analysis cycle #3) is the most recently closed design sizing and mass properties iteration. Upgrades for loss of crew (LDAC-2) and loss of mission (LDAC-3) have been incorporated into the lander baseline design (and its Master Equipment List). Also, recently, Altair has been working requirements analyses (LRAC-1). All nominal data here are from the LDAC-3 analysis cycle. All dispersions results here are from LRAC-1 analyses.

Kos, Larry D.↗
DOE OSTI2021

Design and Experimental Demonstration of an Additive Manufactured Smart Oxy-Methane Burner for High Pressure and Supercritical Combustion

Pressurized oxy combustion-based systems can improve efficiency by recovering latent heat of the steam in the Flue Gas and achieving 90% CO2 capture. In addition, the novel Directly Heated Supercritical Carbon Dioxide (DH-SCO2) power cycles can achieve high thermal efficiencies and provide nearly full carbon capture. Additionally, due to the reduction of flue gas at higher pressure, smaller system size and capital cost reductions are also possible. Recent thermodynamic analysis of the DH-SCO2 cycle performed by the UTEP research team shows that combustion conditions in the vicinity of 300 bar pressure and 1000-1400 K temperature allow for relatively high system efficiencies while operating within the limit of available combustor materials. However, the realization of a directly heated supercritical power cycle requires combustion systems to operate in supercritical conditions and at temperature far below the blowout limit of conventional flames (above 1500 K). The thermodynamic properties along with the combustion properties and kinetics are unexplored at such conditions. Additionally, the interaction of the supercritical environment with energy components is also unknown. High-pressure combustion tests are performed using a smart burner at some intermediate pressure ranges (<20 bar) to help minimize these knowledge gaps. The knowledge obtained from the high-pressure test will assist in understanding the combustion chamber pressurization mechanism, ignition and flame behavior at the elevated pressure. The obtained data will act as a systematic first step in testing at higher pressures of 100 and 300 bar pressures.The primary purpose of this Dissertation is to demonstrate the operability of a low NOx smart burner for oxy-methane combustion at high pressure (< 20 bar) and scalable up to supercritical conditions. A shear co-axial smart burner is designed with a real-time temperature monitoring capability to understand the burner face interaction at high-pressure conditions. In addition, the burner has four independent injection ports to allow the independent injection of fuel and dilution gases in the combustor. The maximum operating capability of the burner is 575 kWth. The burner was fabricated using a Laser Powder Bed Fusion process with Nickel Alloy 718. A powder removal technology was invented comprising ultrasonic vibration, liquid nitrogen exposure and media blasting to remove powders from internal channels of the burner. The burner operability tests are performed in the high-pressure combustor and the swirl combustor. The high-pressure combustor was used to investigate the burner operability and thermal soak back at different pressurized conditions. The experimental tests in high-pressure combustor up to 275 kWth input resulted in 16.5 bar chamber pressure and 198°C thermal soaks back to the burner. The burner was capable of providing the required thermal input within a 3% deviation range. In addition, soot formation occurred at high-pressure tests. The swirl combustor was used to observe the flame stability of the burner. Flame lift-off was observed for jet velocities above 450 m/s. Additionally, lift-off decreased for low co-flow velocities. CO2 dilution experiments showed increased flame instabilities for all conditions above 50% dilution ratios. At lower thermal inputs, partial flame blow-offs occurred for dilution ratios above 50%. All conditions significantly reduced the flame temperature and increased the flame lift-off height. Finally, a 2nd generation AM smart burner was designed using the knowledge from the 1st generation burner experiments. The 2nd generation burner incorporated two sets of swirlers with 0.9 swirl no. A cooling system was also designed for long-duration tests at higher pressures. The thermal input and division of the burner's power are kept the same as the 1st generation burner. The burner is to be fabricated using nickel-alloy 718 for high-pressure handling capability. The design can sustain at high-pressure conditions up to 100 bar.

Islam, Md Nawshad Arslan↗