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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Elastic Stochastic Full Waveform Inversion (eSFWI)

This collaboration between Lawrence Livermore National Security, LLC (LLNS) as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Chevron USA Inc., acting through its Chevron Technical Center division, aimed at developing next-generation computational methods for the Elastic Stochastic Full Waveform Inversion (eSFWI). Seismic imaging is heavily used in the oil and gas industry for identifying and operating subsurface reservoirs. Improved seismic imaging methods can improve productivity, lower costs, and improve operational and environmental safety. This CRADA demonstrated that new high-performance computing (HPC) architectures being rolled out over the next five years can enable unprecedented seismic imaging resolution when using eSFWI techniques to process active seismic data. An open-source computational mini-application was developed, capable of demonstrating near-peak performance for eSFWI algorithms on CPU and GPU enabled HPC platforms. Performance was demonstrated on LLNL HPC systems such as Lassen, as well as on Chevron systems. This project benefited Chevron USA Inc. by demonstrating the potential computational efficiency of their full waveform inversion capabilities used to characterize oil/gas reservoirs, which in turn benefits the public through potential increases in capabilities to perform analysis of leasing sites.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

Electrochemically induced fracture in LLZO: How the interplay between flaw density and electrostatic potential affects operability

Fracture and short circuit in the Li 7 La 3 Zr 2 O 12 (LLZO) solid electrolyte are two key issues that prevent its adoption in battery cells. In this paper, we utilize phase-field simulations that couple electrochemistry and fracture to evaluate the maximum electric potential that LLZO electrolytes can support as a function of crack density. In the case of a single crack, we find that the applied potential at the onset of crack propagation exhibits inverse square root scaling with respect to crack length, analogous to classical fracture mechanics. Here, we further find that the short-circuit potential scales linearly with crack length. In the realistic case where the solid electrolyte contains multiple cracks, we reveal that failure fits the Weibull model. The failure distributions shift to favor failure at lower overpotentials as areal crack density increases. Furthermore, when flawless interfacial buffers are placed between the applied potential and the bulk of the electrolyte, failure is mitigated. When constant currents are applied, current focuses in near-surface flaws, leading to crack propagation and short circuit. We find that buffered samples sustain larger currents without reaching unstable overpotentials and without failing. Our findings suggest several mitigation strategies for improving the ability of LLZO to support larger currents and improve operability.

25 ENERGY STORAGE↗
DOE OSTI2021

Anisotropic colloidal interactions & assembly in AC electric fields

In this work, we match experimental and simulated configurations of anisotropic epoxy colloidal particles in high frequency AC electric fields by identifying analytical potentials for dipole-field and dipole-dipole interactions. We report an inverse Monte Carlo simulation algorithm to determine optimal fits of analytical potentials by matching simulated and experimental distribution functions for non-uniform liquid, liquid crystal, and crystal microstructures in varying amplitude electric fields. Two potentials that include accurate particle volume and dimensions along with a concentration dependent prefactor quantitatively capture experimental observations. At low concentrations, an effective ellipsoidal point dipole potential works well, whereas a novel stretched point dipole potential is found to be suitable at all concentrations, field amplitudes, and degrees of ordering. The simplicity, accuracy, and adjustability of the stretched point dipole potential suggest it can be applied to model field mediated microstructures and assembly of systematically varying anisotropic particle shapes.

42 ENGINEERING↗
DOE OSTI2025

First-Principles Studies on Sc 2 RuZ (Z = Si, Ge, Sn) Inverse Heusler Alloys: Structural, Electronic, and Transport Properties

The continuous demand for efficient, nontoxic, and thermally stable materials for room-temperature energy conversion motivates the exploration of novel thermoelectric systems beyond the traditional magnetic Heusler alloys. While full and half-Heusler compounds, especially Co-, Ni-, and Mn-based systems, have demonstrated promising thermoelectric properties, their typically high operating temperatures and magnetic complexities limit their applicability in ambient thermal management. In this context, we investigate whether Sc-based inverse Heusler alloys can offer a viable nonmagnetic alternative with competitive thermoelectric performance. In this work, we perform a systematic first-principles study of the inverse Heusler compounds Sc 2 RuZ (Z = Si, Ge, Sn), focusing on their structural, electronic, mechanical, and thermodynamic-thermoelectric properties. Density Functional Theory (DFT) was employed to compute optimized lattice structures and band dispersion, while dynamical stability was assessed via phonon calculations. Thermoelectric transport coefficients, including Seebeck coefficient, electrical conductivity, and thermal conductivity, were estimated using the semiclassical Boltzmann transport theory within the constant relaxation time approximation. Our results show that all Sc 2 RuZ compounds are thermodynamically stable semiconductors with indirect band gaps of 0.12–0.16 eV and exhibit high elastic moduli, especially Sc 2 RuSn, which demonstrates superior stiffness and incompressibility. Importantly, all compounds display promising room-temperature thermoelectric characteristics, including high Seebeck coefficients and power factors. These findings reveal that Sc 2 RuZ alloys represent a rare class of stable, nonmagnetic inverse Heusler semiconductors with intrinsic thermoelectric potential at room temperature, unlike many existing Heusler systems optimized for spintronics or high-temperature operation. This work expands the known design space for Heusler-based thermoelectrics and offers a theoretical basis for experimental realization of efficient, low-temperature, nonmagnetic thermoelectric materials.

alloys↗
DOE OSTI2024

Groundwater Circulation Within the Mountain Block: Combining Flow and Transport Models With Magnetotelluric Observations to Untangle Its Nested Nature

Abstract Mountains are vital water sources for humans and ecosystems, continuously replenishing lowland aquifers through surface runoff and mountain recharge. Quantifying these fluxes and their relative importance is essential for sustainable water resource management. However, our mechanistic understanding of the flow and transport processes determining the connection between the mountain block and the basin aquifer remains limited. Traditional conceptualizations assume groundwater circulation within the mountain block is predominantly shallow. This view neglects the role of deep groundwater flowpaths significantly contributing to the water, solute, and energy budgets. Overcoming these limitations requires a holistic characterization of the multiscale nature of groundwater flow along the mountain‐to‐valley continuum. As a proof‐of‐concept, we use a coupled groundwater flow and transport model to design a series of numerical experiments that explore the role of geology, topography, and weathering rates in groundwater circulation and their resulting resistivity patterns. Our results show that accumulating solutes near stagnation zones create contrasting electrical resistivity patterns that separate local, intermediate, and regional flow cells, presenting a target for magnetotelluric observations. To demonstrate the sensitivity of magnetotelluric data to features in our resistivity models, we use the MARE2DEM electromagnetic modeling code to perform forward and inverse simulations. This study highlights the potential of magnetotelluric surveys to image the resistivity structure resulting from multiscale groundwater circulation through relatively impervious crystalline basement rocks in mountainous terrains. This capability could change our understanding of the critical zone, offering a holistic perspective that includes deep groundwater circulation and its role in conveying solutes and energy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Coherence from multiorbital tunneling ionization of molecules

We present a simple and general coherence model for multiorbital tunnel ionization of molecules, which we incorporate into our previously developed density-matrix approach for sequential double ionization. Here, the influence of this coherence is investigated through simulations of single ionization and sequential double ionization of N 2 and O 2 using few-cycle near-infrared laser pulses. In the case of single ionization, our results reveal the crucial role played by this coherence in generating population inversion in N 2 + , suggesting a potential mechanism for air lasing. Regarding sequential double ionization, we observe only minor changes in the kinetic energy release spectra when the coherence is included, while noticeable differences in the angle-dependent dication yield for both N 2 and O 2 are found. Based on these findings, we recommend the inclusion of multiorbital tunnel ionization coherence in models for single ionization of general molecules, while suggesting that it can be safely neglected in the case of sequential double ionization.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Design, tuning, and blackbox optimization of laser systems

Chirped pulse amplification (CPA) and subsequent nonlinear optical (NLO) systems constitute the backbone of myriad advancements in semiconductor manufacturing, communications, biology, defense, and beyond. Accurately and efficiently modeling CPA+NLO-based laser systems is challenging because of the complex coupled processes and diverse simulation frameworks. Our modular start-to-end model unlocks the potential for exciting new optimization and inverse design approaches reliant on data-driven machine learning methods, providing a means to create tailored CPA+NLO systems unattainable with current models. To demonstrate this new, to our knowledge, technical capability, we present a study on the LCLS-II photo-injector laser, representative of a high-power and spectro-temporally non-trivial CPA+NLO system.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential for specialized applications in segmented inverse beta decay (IBD) neutrino detectors, astronomy, machine learning, and more. Although this method may easily extend to 3D scalar fields, our focus here is on 2D real-valued fields as it directly applies to directionality. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Data Analysis, Statistics and Probability (physics↗
DOE OSTI2020

The significant impact of aerosol vertical structure on lower atmosphere stability and its critical role in aerosol–planetary boundary layer (PBL) interactions

The aerosol–planetary boundary layer (PBL) interaction was proposed as an important mechanism to stabilize the atmosphere and exacerbate surface air pollution. Despite the tremendous progress made in understanding this process, its magnitude and significance still have large uncertainties and vary largely with aerosol distribution and meteorological conditions. In this study, we focus on the role of aerosol vertical distribution in thermodynamic stability and PBL development by jointly using micropulse lidar, sun photometer, and radiosonde measurements taken in Beijing. Despite the complexity of aerosol vertical distributions, cloud-free aerosol structures can be largely classified into three types: well-mixed, decreasing with height, and inverse structures. The aerosol–PBL relationship and diurnal cycles of the PBL height and PM 2.5 associated with these different aerosol vertical structures show distinct characteristics. The vertical distribution of aerosol radiative forcing differs drastically among the three types, with strong heating in the lower, middle, and upper PBL, respectively. Such a discrepancy in the heating rate affects the atmospheric buoyancy and stability differently in the three distinct aerosol structures. Absorbing aerosols have a weaker effect of stabilizing the lower atmosphere under the decreasing structure than under the inverse structure. As a result, the aerosol–PBL interaction can be strengthened by the inverse aerosol structure and can be potentially neutralized by the decreasing structure. Moreover, aerosols can both enhance and suppress PBL stability, leading to both positive and negative feedback loops. This study attempts to improve our understanding of the aerosol–PBL interaction, showing the importance of the observational constraint of aerosol vertical distribution for simulating this interaction and consequent feedbacks.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1997

Mars Surveyor '98 Landers MVACS Robotic Arm Control System Design Concepts

This paper describes the control system design concepts for the Mars Volatiles and Climate Surveyor (MVACS) Robotic Arm which supports the scientific investigations to be conducted as part of the Mars Surveyor '98 Lander project. Solutions are presented to some of the problems encountered in this demanding space application with its tight constraints on mass, power, volume, and computing resources. Problems addressed include 4-DOF forward and inverse kinematics, trajectory planning to minimize potential impact damage, joint drive train protection, Lander tilt prevention, hardware fault monitoring, and collision avoidance.

Bonitz, Robert G.↗
NASA NTRS2000

Mass Spectrometric and Langmuir Probe Measurements in Inductively Coupled Plasmas in Ar, CHF3/Ar and CHF3/Ar/O2 Mixtures

Absolute fluxes and energy distributions of ions in inductively coupled plasmas of Ar, CHF3/Ar, and CHF3/Ar/O2 have been measured. These plasmas were generated in a Gaseous Electronics Conference (GEC) cell modified for inductive coupling at pressures 10-50 mTorr and 100-300 W of 13.56 MHz radio frequency (RF) power in various feedgas mixtures. In pure Ar plasmas, the Ar(+) flux increases linearly with pressure as well as RF-power. Total ion flux in CHF3 mixtures decreases with increase in pressure and also CHF3 concentration. Relative ion fluxes observed in the present studies are analyzed with the help of available cross sections for electron impact ionization and charge-exchange ion-molecule reactions. Measurements of plasma potential, electron and ion number densities, electron energy distribution function, and mean electron energy have also been made in the center of the plasma with a RF compensated Langmuir probe. Plasma potential values are compared with the mean ion energies determined from the measured ion energy distributions and are consistent. Electron temperature, plasma potential, and mean ion energy vary inversely with pressure, but increase with CHF3 content in the mixture.

Kim, J. S.↗
NASA NTRS2006

Using CO2:CO Correlations to Improve Inverse Analyses of Carbon Fluxes

Observed correlations between atmospheric concentrations of CO2 and CO represent potentially powerful information for improving CO2 surface flux estimates through coupled CO2-CO inverse analyses. We explore the value of these correlations in improving estimates of regional CO2 fluxes in east Asia by using aircraft observations of CO2 and CO from the TRACE-P campaign over the NW Pacific in March 2001. Our inverse model uses regional CO2 and CO surface fluxes as the state vector, separating biospheric and combustion contributions to CO2. CO2-CO error correlation coefficients are included in the inversion as off-diagonal entries in the a priori and observation error covariance matrices. We derive error correlations in a priori combustion source estimates of CO2 and CO by propagating error estimates of fuel consumption rates and emission factors. However, we find that these correlations are weak because CO source uncertainties are mostly determined by emission factors. Observed correlations between atmospheric CO2 and CO concentrations imply corresponding error correlations in the chemical transport model used as the forward model for the inversion. These error correlations in excess of 0.7, as derived from the TRACE-P data, enable a coupled CO2-CO inversion to achieve significant improvement over a CO2-only inversion for quantifying regional fluxes of CO2.

Palmer, Paul I.↗
NASA NTRS2009

Monitoring Exchange of CO2 - A KISS Workshop Report 2009

The problem and context: Can top-down estimates of carbon dioxide (CO2) fluxes resolve the anthropogenic emissions of China, India, the United States, and the European Union with an accuracy of +/-10% or better?The workshop "Monitoring Exchange of Carbon Dioxide" was convened at the Keck Institute for Space Studies in Pasadena, California in February 2010 to address this question. The Workshop brought together an international, interdisciplinary group of 24 experts in carbon cycle science, remote sensing, emissions inventory estimation, and inverse modeling. The participants reviewed the potential of space-based and sub-orbital observational and modeling approaches to monitor anthropogenic CO2 emissions in the presence of much larger natural fluxes from the exchange of CO2 between the land, atmosphere, and ocean. This particular challenge was motivated in part by the NRC Report "Verifying Greenhouse Gas Emissions" [Pacala et al., 2010]. This workshop report includes several recommendations for improvements to observing strategies and modeling frameworks for optimal and cost-effective monitoring of carbon exchange

Greenhouse Gas Emissions↗
NASA NTRS2017

Iterative Discrete Ordinates Solution of the Equation for the Surface-Reflected Radiance

This paper presents a new method of numerical solution of the integral equation for the radiance reflected from an anisotropic surface. The equation relates the radiance at the surface level with BRDF and solutions of the standard radiative transfer problems for a slab with no reflection on its surfaces. It is also shown that the kernel of the equation satisfies the condition of the existence of a unique solution and the convergence of the successive approximations to that solution. The developed method features two basic steps: discretization on a 2D quadrature, and solving the resulting system of algebraic equations with successive over-relaxation method based on the Gauss-Seidel iterative process. Presented numerical examples show good coincidence between the surface-reflected radiance obtained with DISORT and the proposed method. Analysis of contributions of the direct and diffuse (but not yet reflected) parts of the downward radiance to the total solution is performed. Together, they represent a very good initial guess for the iterative process. This fact ensures fast convergence. The numerical evidence is given that the fastest convergence occurs with the relaxation parameter of 1 (no relaxation). An integral equation for BRDF is derived as inversion of the original equation. The potential of this new equation for BRDF retrievals is analyzed. The approach is found not viable as the BRDF equation appears to be an ill-posed problem, and it requires knowledge the surface-reflected radiance on the entire domain of both Sun and viewing zenith angles.

Alexander Radkevich↗
DOE OSTI2022

Metal dissolution from first principles: Potential-dependent kinetics and charge transfer

Metal dissolution and its inverse process are integral to both corrosion and electrodeposition; however, many mechanistic details regarding the dissolution process are challenging to decipher. These include how ion dissolution kinetics and charge transfer are influenced by the competition between metal and solvent interactions under an electrode potential. Here, in this work, we introduce a computational framework based on density functional theory with grand-canonical treatment of electrons to directly predict the potential energy landscape for metal dissolution at a constant potential. Using aluminum as an example, we demonstrate that dissolution kinetics is governed by competing kinetics between two physical processes associated with metal–metal bond breaking and ion-migration within the electrical double layer, respectively. We identify a kinetic transition between regimes dominated by each of these processes and show that this transition depends on the operating electrode potential, among other key factors. It is further found that kinetics and thermodynamics of these processes can be described with a simple, one-parameter Marcus-theory-type model. Beyond offering new understanding of charge transfer during dissolution, our simulation protocol provides a recipe for directly predicting other important quantities in electrochemical reactions from first principles that are difficult to measure, such as the symmetry factor.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Transonic cascade flow analysis using viscous/inviscid coupling concepts

Transonic two-dimensional cascade flows have been analyzed using viscous/inviscid coupling concepts. A full potential cascade code is coupled with an inverse integral boundary layer/wake method that permits calculation of separated laminar or turbulent flow. The semi-inverse coupling method of Wigton converges slowly in the case of a strong shock in the region between the shock and the trailing edge. The location of a strong shock is not well predicted by the coupling method, which indicates the need for an entropy correction in the potential code or the inclusion of a shock-boundary layer interaction module.

Olling, C. R.↗
DOE OSTI2023

Double-Atom Catalysts Featuring Inverse Sandwich Structure for CO 2 Reduction Reaction: A Synergetic First-Principles and Machine Learning Investigation

Electrocatalytic CO 2 reduction reactions (CO 2 RR) based on scalable and highly efficient catalysis provide an attractive strategy for reducing CO 2 emissions. Here in this work, we combined first-principles density functional theory (DFT) and machine learning (ML) to comprehensively explore the potential of double-atom catalysts (DACs) featuring an inverse sandwich structure anchored on defective graphene (gra) to catalyze CO 2 RR to generate C 1 products. We started with five homonuclear M 2 ⊥gra (M = Co, Ni, Rh, Ir, and Pt), followed by 127 heteronuclear MM'⊥gra (M = Co, Ni, Rh, Ir, and Pt, M' = Sc–Au). Stable DACs were screened by evaluating their binding energy, formation energy, and dissolution potential of metal atoms, as well as conducting first-principles molecular dynamics simulations with and without solvent water molecules. Based on DFT calculations, Rh 2 ⊥gra DAC was found to outperform the other four homonuclear DACs and the Rh-based single- and double-atom catalysts of noninverse sandwich structures. Out of the 127 heteronuclear DACs, 14 were found to be stable and have good catalytic performance. An ML approach was adopted to correlate key factors with the activity and stability of the DACs, including the sum of radii of metal and ligand atoms (d M–M' , d M–C , and d M'–C ), the sum and difference of electronegativity of two metal atoms (P M + P M' , P M – P M '), the sum and difference of first ionization energy of two metal atoms (I M + I M' , I M – I M '), the sum and difference of electron affinity of two metal atoms (A M + A M' , A M – A M '), and the number of d-electrons of the two metal atoms (Nd). The obtained ML models were further used to predict 154 potential electrocatalysts out of 784 possible DACs featuring the same inverse sandwich configuration. Overall, this work not only identified promising CO 2 RR DACs featuring the reported inverse sandwich structure but also provided insights into key atomic characteristics associated with high CO 2 RR activity.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS1977

Influence of temperature and the role of chromium on the kinetics of sulfidation of 310 stainless steel

The sulfidation of 310 stainless steel was studied over the temperature range from 910 K to 1285 K. By adjusting the ratio of hydrogen sulfide, variations in sulfur potential were obtained. The effect of temperature on sulfidation was determined at three different sulfur potentials: 39/sqNm, 0.014/sqNm, and 0.00015/sqNm. All sulfide scales contained one or two surface layers in addition to a subscale. The second outer layer (OL-II), furthest from the alloy, contained primarily Fe-Ni-S. The first outer layer (OL-I), nearest the subscale, contained FE-Cr-S. The subscale consisted of sulfide inclusions in the metal matrix. At a given temperature and sulfur potential, the weight gain data obeyed the parabolic rate law after an initial transient period. The parabolic rate constants obtained at the sulfur potential of 39/sqNm did not show a break when the logarithm of the rate constant was plotted as a function of the inverse of absolute temperature. Sulfidation carried out at sulfur potentials below 0.02/sqNm, however, did show a break at 1145 K, which is termed as the transition temperature. This break was found to be associated with the changes which had occurred in the Fe:Cr ratio of OL-I. Below the transition temperature the activation energy was found to be approximately 125 kj/mole. Above the transition temperature the rate of sulfidation decreased with temperature but dependent on the Fe:Cr ratio in the iron-chromium-sulfide layers of the OL-I. A reaction mechanism consistent with the experimental results has been proposed.

Rao, D. B.↗