Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential inversion”

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 415 records · Page 23

STILGAR End-of-Project Report

The Subsurface Tunnel Imaging LeveraGed by Analysis of Rayleigh wave ellipticity (STILGAR) project demonstrated an integrated geophysical approach for detecting, locating, and characterizing underground structural changes using dense seismic arrays and advanced inversion techniques. Field campaigns were conducted at two operational mines—the Redmond salt mine (Utah) and Graymont Pleasant Gap limestone mine (Pennsylvania)—providing real-world testbeds for monitoring anthropogenic subsurface activity. At the Redmond salt mine, seismic interferometry combined with back-projection inversion successfully identified continuous, low-amplitude signals from mining operations. The approach differentiated stationary from migrating anthropogenic sources, captured daily operational cycles, and validated the potential of passive seismic monitoring for remote detection of underground activity. At the Graymont Pleasant Gap mine, two dense seismic deployments in the spring and fall of 2023 generated over 4 TB of high-resolution data. Key outcomes included the relocation of 199 underground and 8 surface explosions with accuracies within tens of meters and the development of a 3D P-wave velocity model using the triple-difference tomography algorithm (tomoTD) that resolved major structural features such as the mine entrance, low-velocity tunnels, and roof-collapse areas. Ambient noise cross-correlation and back-projection analyses revealed persistent sources linked to ongoing mining activity, whereas horizontal-to-vertical spectral ratio (HVSR) and ellipticity studies confirmed stable site responses across seasons and identified soil thickness trends consistent with regional erosional and depositional processes. Checkerboard and sensitivity tests further validated the robustness of the tomographic results. Overall, the findings emphasize that although significant progress has been made in subsurface imaging, further work is needed to enhance the detection and localization of underground structures. Accurate imaging requires higher frequencies, yet anthropogenic sources tend to dominate the seismic record at those frequencies, and high-frequency surface waves are affected by higher modes that complicate interpretation. The improved detection and localization of human-induced signals enabled detailed temporal and spatial mapping of daily mine operations, demonstrating the feasibility of continuous anthropogenic source monitoring. Sensitivity to signals from nontraditional sources, such as fan operations, highlights the broader applicability of this approach to other industrial environments in which continuous and impulsive signals are present. The field campaigns produced a substantial volume of high-quality seismic data, supporting the development and testing of new methods for seismic source characterization and subsurface imaging. Future deployments should include sensors capable of recording lower frequencies to probe deeper structures, increase bandwidth to enhance resolution and sensitivity to both shallow and deep targets, and collect additional large-scale datasets to refine imaging and source characterization techniques. Moreover, conducting 3D modeling studies of seismic wavefields at higher frequencies will provide a better understanding of wave scattering and cavity–wavefield interactions in complex underground environments. In conclusion, the STILGAR project demonstrated that integrated seismic monitoring can effectively characterize underground operations, capturing both natural and anthropogenic signals. The approaches developed provide a foundation for improved detection, localization, and imaging of subsurface structures and are directly transferable to broader industrial monitoring applications.

58 GEOSCIENCES↗

Charge-exchange plasma generated by an ion thruster

The use of high voltage solar arrays greatly reduces or eliminates power processing requirements in space electric propulsion systems. This use also requires substantial areas of solar array to be at high positive potential relative to space and most of the spacecraft. The charge exchange plasma conducts electrons from the ion beam to such positive surfaces, and thereby electrically load the high voltage solar array. To evaluate this problem, the charge-exchange plasma generated by an ion beam was investigated experimentally. Based upon the experimental data, a simple model was derived for the charge-exchange plasma. This model is conservative in the sense that both the electron/ion density and the electron current density should be equal to, or less than, the preducted value for all directions in the hemisphere upstream of the ion beam direction. Increasing the distance between a positive potential surface (such as a high voltage solar array) and the thruster is the simplest way to control interactions. Both densities and currents vary as the inverse square of this distance.

Kaufman, H. R.↗

Application of a viscous-inviscid interaction method to predict transonic separated flows

A viscous-inviscid interaction numerical procedure is used to predict turbulent transonic separated flows over axisymmetric bodies. The viscous flow solutions are obtained by solving the boundary-layer equations inversely in a coupled manner by a finite-difference scheme; the inviscid flow is determined by a finite-difference solution to the full-potential equation in conservative form. Parameters affecting the convergence rate of the interaction procedure are discussed. The predictions of a new turbulence model proposed by Johnson and King are compared with experimental data and numerical solutions to the Navier-Stokes equations for flows with large regions of separation.

Lee, D. S.↗

The Next Generation Marshall Grazing Incidence X-Ray Spectrometer (MaGIXS) Sounding Rocket Experiments

The Marshall Grazing Incidence X-ray Spectrometer (MaGIXS) is a sounding rocket experiment that carries an X-ray slitless imaging spectrograph to observe spectrally dispersed soft X-ray (SXR) images of the solar corona over a wide field-of-view. The first flight of MaGIXS occurred on 30 July 2021, during which several emission lines from coronal structures including X-ray bright points were observed. Analysis of MaGIXS-1 data demonstrated the application of new inversion techniques and opened a new arena of inverting complex spectroheliogram data, which has spatial-spectral information overlapped. The second flight of MaGIXS, with a simplified optical design, is scheduled for 2024. The goal of MaGIXS-2 is to observe high temperature diagnostic emission lines within an active region core. MaGIXS-1 discovered dominant missing emission lines near 15Å wavelength region, which hosts several closely spaced satellite lines of Fe XVII, Fe XVI, and Fe XV ions, arising from relatively cool plasma that peaks near 2MK. These lines are expected to be enhanced at lower temperatures and offers one of the most unique diagnostics to measure electron temperature, as well as signatures of equilibrium state of the plasma in active regions, which has not been studied so far. The science goal of the third sounding rocket flight of MaGIXS is to determine coronal heating parameters, such as the spatial and temporal properties of coronal heating events, by measuring discriminating observations, such as electron temperature, effective temperature, density, abundance, and departures from thermal equilibrium of the plasma, by observing the Sun in the SXR wavelength range. For this, the MaGIXS instrument will be upgraded with a new X-ray telescope mirror that will provide higher spatial and spectral resolution and throughput, allowing for spectral lines to be observed at the relevant spatial and temporal scales. In addition, MaGIXS-3 will also carry The Resolving Inversion Context X-ray Spectrometer (TRICXS), a high dispersion Bragg crystal spectrometer to spectrally resolve the lines near 15Å, critical to unlock the full diagnostic potential of this wavelength range. Here we will present the preliminary concept design of MaGIXS-3 and discuss the potential observations.

X-ray Imaging↗

Surface Charge Distribution Mapping Using an Electrostatic Probe

Undesired accumulation of charge on dielectric materials causing electrostatic discharges can be an issue in pulsed power systems and electrical systems in general. As such, an understanding of surface charge distribution, charge accumulation, and decay is required. An in-house no-contact electrostatic probe designed with a 2 mm resolution was fabricated to measure and map surface charge distribution. Designed as a contactless instrument, the probe consisting of a metal rod and low leakage amplifier circuit probes the electric potential near the surface. A scan of the surface then provides the raw data, a potential distribution. The actual probe response - i.e., the surface charge to voltage transfer function, is accounted for in post-processing. This is accomplished via an Inverse Wien Filter - a technique often applied in image processing - to deconvolve the probe response from the measured data. A commercially available electrostatic probe, the Trek 341B meter with a 3455ET probe, capable of measuring +/- 20 kV was compared to the in- house probe that is designed to cover a wider range from +/- 35 kV. A resolution better than the simple probe resolution is achieved through the distribution’s scanning voltage method and deconvolution. Applying repeated scans, surface charge decay was tracked on various polymer materials to determine the material and environment dependence; materials included were PA6, PTFE, and others. As an example of material dependence, samples charged to 20 kV at 65% humidity experience full charge decay in approximately 45 seconds for PA5 (152 mm dia.), while it took about 100 times longer for PTFE (51 mm dia.) to observe the same decay/redistribution of charge.

36 MATERIALS SCIENCE↗

Characterization of mechanical deformation in aluminum by optical second harmonic generation

Second harmonic generation (SHG), a laser-based surface technique, was used to characterize near-surface defect evolution associated with mechanical deformation in 2024-T3 aluminum due to tensile loads. Electron backscatter diffraction (EBSD) measurements were also performed to determine the geometrically necessary dislocation density before and after tensile deformation. Comparison of the EBSD and SHG results found that dislocation density and SHG signal are inversely proportional; tensile deformation in the plastic regime led to an increase in dislocation density and a decrease in SHG response. A similar trend was seen in the SHG results on as-received (unpolished) Al samples; SHG response decreases with extent of plastic deformation. Furthermore, this study suggests that SHG has potential as a nondestructive evaluation technique.

36 MATERIALS SCIENCE↗

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗

Toward a Comprehensive Carbon Budget for North America: Potential Applications of Adjoint Methods with Diverse Datasets

A detailed mechanistic understanding of the sources and sinks of CO2 will be required to reliably predict future COS levels and climate. A commonly used technique for deriving information about CO2 exchange with surface reservoirs is to solve an "inverse problem," where CO2 observations are used with an atmospheric transport model to find the optimal distribution of sources and sinks. Synthesis inversion methods are powerful tools for addressing this question, but the results are disturbingly sensitive to the details of the calculation. Studies done using different atmospheric transport models and combinations of surface station data have produced substantially different distributions of surface fluxes. Adjoint methods are now being developed that will more effectively incorporate diverse datasets in estimates of surface fluxes of CO2. In an adjoint framework, it will be possible to combine CO2 concentration data from long-term surface monitoring stations with data from intensive field campaigns and with proposed future satellite observations. A major advantage of the adjoint approach is that meteorological and surface data, as well as data for other atmospheric constituents and pollutants can be efficiently included in addition to observations of CO2 mixing ratios. This presentation will provide an overview of potentially useful datasets for carbon cycle research in general with an emphasis on planning for the North American Carbon Project. Areas of overlap with ongoing and proposed work on air quality/air pollution issues will be highlighted.

Andrews, A.↗

Inverse prediction of PuO2 processing conditions using Bayesian seemingly unrelated regression with functional data

Over the past decade, a variety of innovative methodologies have been developed to better characterize the relationships between processing conditions and the physical, morphological, and chemical features of special nuclear material (SNM). Different processing conditions generate SNM products with different features, which are known as “signatures” because they are indicative of the processing conditions used to produce the material. These signatures can potentially allow a forensic analyst to determine which processes were used to produce the SNM and make inferences about where the material originated. This article investigates a statistical technique for relating processing conditions to the morphological features of PuO 2 particles. We develop a Bayesian implementation of seemingly unrelated regression (SUR) to inverse-predict unknown PuO 2 processing conditions from known PuO 2 features. Model results from simulated data demonstrate the usefulness of the technique. Applied to empirical data from a bench-scale experiment specifically designed with inverse prediction in mind, our model successfully predicts nitric acid concentration, while results for Pu concentration and precipitation temperature were equivalent to a simple mean model. Our technique compliments other recent methodologies developed for forensic analysis of nuclear material and can be generalized across the field of chemometrics for application to other materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A combined direct/inverse three-dimensional transonic wing design method for vector computers

A three-dimensional transonic-wing design algorithm for vector computers is developed, and the results of sample computations are presented graphically. The method incorporates the direct/inverse scheme of Carlson (1975), a Cartesian grid system with boundary conditions applied at a mean plane, and a potential-flow solver based on the conservative form of the full potential equation and using the ZEBRA II vectorizable solution algorithm of South et al. (1980). The accuracy and consistency of the method with regard to direct and inverse analysis and trailing-edge closure are verified in the test computations.

Weed, R. A.↗

Chemokinetic motility responses of the cyanobacterium oscillatoria terebriformis

Oscillatoria terebriformis, a gliding, filamentous, thermophilic cyanobacterium, exhibited an inhibition of gliding motility upon exposure to fructose. The observed response was transient, and the duration of nonmotility was directly proportional to the concentration of fructose. Upon resumption of motility, the rate of motility was also inversely proportional to the concentration of fructose. Sulfide caused a similar response. The effect of sulfide was specific and not due to either anoxia or negative redox potential. Exposure to glucose, acetate, lactate, or mat interstitial water did not elicit any motility response.

Richardson, Laurie L.↗

Error in predicting hydrosol backscattering from remotely sensed reflectance

Monte Carlo simulations are carried out to determine the error in the inversion of backscattering from remotely sensed reflectance when geometrical shape factors of the light field are assumed to be unity. The results show that error in backscattering inversion can vary from a 40% overestimation to a 20% underestimation and is dependent on the solar angle and the hydrosol consituents contributing to backscattering. The simulations also demostrate that for chlorophyll concentrations ranging from 0.05 to 20 mg/cu m the most dramatic change in the geometrical shape factors occurs near 1.0 to 1.5 mg/cu m chlorophyll. The potential importance of bacteria in influencing the shape factor and the subsequent effect of bacteria on the retrieval of the backscattering from remote sensing reflectance are shown. Quartzlike material's strong impact on geometrical shape factors and errors of retrieval of backscattering at low chlorophyll concentrations are also demonstrated. Remote sensing reflectance inversion schemes must include information about the backscattering function to be successful.

Weidemann, A. D.↗

Relocating microseismicity from downhole monitoring of the Decatur CCS site using a modified double-difference algorithm

SUMMARY The injection of CO2 at the Decatur carbon capture and storage site has generated significant microseismic activity, which occurs in distinct spatial clusters up to approximately 2.2 km from the primary injection well. Accurate and precise event locations are vital for the characterization of the microseismicity to help understand the reservoir response to the CO2 injection, whilst enabling the identification of minor faults and fractures below the resolution of conventional active seismic imaging. However, microseismic monitoring of fluid injection sites, such as Decatur, is often performed using a network of borehole sensors often from a single well. While these downhole sensors have excellent detection capabilities, their poor azimuthal coverage limits the ability to precisely determine event locations. We have developed a modified double-difference relocation algorithm suitable for both 1-D and 3-D velocity models, and which incorporates differential back azimuth observations to allow the benefits of the original double-difference algorithm to be applicable to a downhole microseismic monitoring setting. Applying the modified double-difference algorithm to the microseismicity at Decatur, we have successfully relocated 4293 events. The relocation included over 59 million observations for 757 285 event pairs, split across seven geographic regions. Despite the majority of observations being recorded in only two boreholes, with an almost identical azimuthal coverage, the results have shown to be reliable with significantly reduced residuals and low uncertainties associated with the final locations. We have analysed the residuals in terms of their association with each geographic region, data type, station and individual events, to fully appreciate their influence in the inversion and the fit of the data to the final set of event locations. For each region, the relocated seismicity has become less diffuse with improved clustering, and with newly visible linear features often orientated in a NE–SW direction. These results show the potential improvements that can be made to microseismic event locations recorded by a borehole network with a limited and variable azimuthal distribution.

58 GEOSCIENCES↗

Creep behavior of tungsten/niobium and tungsten/niobium-1 percent zirconium composites

The creep behavior and microstructural stability of tungsten fiber reinforced niobium and niobium 1 percent zirconium was determined at 1400 and 1500 K in order to assess the potential of this material for use in advanced space power systems. The creep behavior of the composite materials could be described by a power law creep equation. A linear relationship was found to exist between the minimum creep rate of the composite and the inverse of the composite creep rupture life. The composite materials had an order of magnitude increase in stress to achieve 1 percent creep strain and in rupture strength at test temperatures of 1400 and 1500 K compared to unreinforced material. The composite materials were also stronger than the unreinforced materials by an order of magnitude when density was taken into consideration. Results obtained on the creep behavior and microstructural stability of the composites show significant potential improvement in high temperature properties and mass reduction for space power system components.

Petrasek, Donald W.↗

Multimodal spectroscopy of order parameter distributions

We present a multimodal spectroscopic paradigm that enables independent measurement of charge and spin degrees of freedom (DOF) in strongly correlated materials. This spin-based technique probes symmetry-specific Hamiltonian parameters by analyzing how the time delay between applied pulses (τ) affects the response. We demonstrate ways in which charge DOF that couple through the quadrupolar interaction (inversion symmetric) can be independently measured even in the presence of large magnetic noise (inversion asymmetric). The method quantifies both the strength of the interactions and their distribution (noise). We provide protocols to directly and independently measure the distribution of interaction strengths, even when the average value of the interaction is zero. By independently measuring distributions of different forms of disorder, this methodology can elucidate which microscopic symmetry drives a phase transition. Here, we discuss potential applications to study complex phase transitions in strongly interacting quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quenching of single-particle strength inferred from nucleon-removal transfer reactions on 15 C

The difference in the proton and neutron separation energies (ΔS) of the weakly bound 15 C ground state is -19.86 MeV, an extreme value. Data from intermediate-energy heavy-ion induced (HI-induced) knockout reactions on nuclei spanning -20 ≲ ΔS ≲ +20 MeV, suggest that the degree to which single-particle strength is quenched, R s , has a negative correlation with ΔS, decreasing from unity around -20 MeV to around 0.2 at +20 MeV. For the 15 C ground state (R s = 0.96 (4) in HI-induced knockout), contrasting results have recently been obtained via the neutron-adding transfer reaction, which reveal a value of R s = 0.64 (15), similar to the value observed at modest and more extreme values of ΔS with reaction probes other than HI knockout. In order to explore the any potential differences between adding and removing processes in transfer reactions at extreme ΔS, single-neutron removal transfer reactions on 15 C were performed at 7.1 MeV/u in inverse kinematics. The removal of a valence neutron in 2s 1/2 orbit using both (p, d) and (d, t) reactions shows consistent quenching factors and agrees with those from the neutron-adding reaction. The present results, which can be compared with neutron knockout reaction, suggest that correlations, represented by the quenching factor, show limited dependence on neutron-proton asymmetry under the most extreme asymmetry conditions so far achieved in transfer reactions.

quenching factor↗

Correlation-driven electron-hole asymmetry in graphene field effect devices

Electron-hole asymmetry is a fundamental property in solids that can determine the nature of quantum phase transitions and the regime of operation for devices. The observation of electron-hole asymmetry in graphene and recently in twisted graphene and moiré heterostructures has spurred interest into whether it stems from single-particle effects or from correlations, which are core to the emergence of intriguing phases in moiré systems. Here, we report an effective way to access electron-hole asymmetry in 2D materials by directly measuring the quasiparticle self-energy in graphene/Boron Nitride field-effect devices. As the chemical potential moves from the hole to the electron-doped side, we see an increased strength of electronic correlations manifested by an increase in the band velocity and inverse quasiparticle lifetime. These results suggest that electronic correlations intrinsically drive the electron-hole asymmetry in graphene and by leveraging this asymmetry can provide alternative avenues to generate exotic phases in twisted moiré heterostructures.

36 MATERIALS SCIENCE↗

Advanced Inversion Algorithms for Scientific Data Analysis [Slides]

Accurate subsurface characterization is crucial for all subsurface energy exploration. Accurate characterization of uncertain subsurface properties is also critical for monitoring storage of CO 2 , estimating pathways of subsurface contaminant transport, and monitoring potential nuclear explosions for treaty verification. This research will advance our world-leading subsurface sensing capabilities that are crucial for LANL missions in energy security (geothermal energy, oil/gas resource exploration, geologic carbon storage) and nuclear security (facility monitoring, detonation detection).

47 OTHER INSTRUMENTATION↗