Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium reconstruction”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Measurement of L-shell emission from mid-Z targets under non-LTE conditions using Transmission Grating Spectrometer and DANTE power diagnostics

In this work, we present the measurement of L-band emission from buried Sc/V targets in experiments performed at the OMEGA laser facility. The goal of these experiments was to study non-local thermodynamic equilibrium plasmas and benchmark atomic physics codes. The L-band emission was measured simultaneously by the time resolved DANTE power diagnostic and the recently fielded time integrated Soreq-Transmission Grating Spectrometer (TGS) diagnostic. The TGS measurement was used to support the spectral reconstruction process needed for the unfolding of the DANTE data. The Soreq-TGS diagnostic allows for broadband spectral measurement in the 120 eV–2000 eV spectral band, covering L- and M-shell emission of mid- and high-Z elements, with spectral resolution λ/Δλ = 8–30 and accuracy better than 25%. The Soreq-TGS diagnostic is compatible with ten-inch-manipulator platforms and can be used for a wide variety of high energy density physics, laboratory astrophysics, and inertial confinement fusion experiments.

Fraenkel, M. (ORCID:0000000291939544)↗
DOE OSTI2026

Direct estimation of the density of states for fermionic systems

Simulating time evolution is one of the most natural applications of quantum computers and is thus one of the most promising prospects for achieving practical quantum advantage. Here, we develop quantum algorithms to extract thermodynamic properties by estimating the density of states (DOS), which is a central object in quantum statistical mechanics. We introduce several key innovations that significantly improve the practicality and extend the generality of previous techniques. First, our approach allows one to estimate the DOS only for a specific subspace of the full Hilbert space. This is crucial for fermionic systems, since both canonical and grand canonical ensemble thermal equilibrium properties depend on subspaces of fixed number. Second, in our approach, by time evolving very simple, random initial states, such as randomly chosen computational basis states, we can exactly recover the DOS on average. Third, due to circuit-depth limitations, we only reconstruct the DOS up to a convolution with a Gaussian window—thus all imperfections that shift the energy levels by less than the width of the convolution window will not significantly affect the estimated DOS. For these reasons, we find the approach is a promising candidate for early quantum advantage as even short-time, noisy dynamics can yield a semiquantitative reconstruction of the DOS (convolution with a broad Gaussian window), while early fault-tolerant devices will likely enable higher-resolution DOS reconstruction through longer time evolutions. We demonstrate the practicality of our approach in representative Fermi-Hubbard and spin models and indeed find that our approach is highly robust against algorithmic errors in the time evolution and against gate noise. We further demonstrate that our approach is compatible with noisy intermediate-scale quantum (NISQ) computing NISQ-friendly variational techniques, introducing and leveraging a technique for variational time evolution.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2014

Microstructures of Hibonite From an ALH A77307 (CO3.0) CAI: Evidence for Evaporative Loss of Calcium

Hibonite is a comparatively rare, primary phase found in some CAIs from different chondrite groups and is also common in Wark-Lovering rims [1]. Hibonite is predicted to be one of the earliest refractory phases to form by equilibrium condensation from a cooling gas of solar composition [2] and, therefore, can be a potential recorder of very early solar system processes. In this study, we describe the microstructures of hibonite from one CAI in ALH A77307 (CO3.0) using FIB/TEM techniques in order to reconstruct its formational history.

Han, Jangmi↗
DOE OSTI2025

Single-shot electro-optic sampling with arbitrary terahertz polarization

With the recent development of diversity electro-optic sampling (DEOS), significant progress has been made in the range of applicability of single-shot EOS measurements, allowing broadband THz waveforms to be captured in a single shot over large temporal windows. In addition to the decrease in acquisition time compared to standard multishot data acquisition, this technique allows measurements on systems far from equilibrium with large shot-to-shot noise or with irreversible or poorly repeatable dynamics. Although DEOS has been demonstrated and verified for linearly polarized THz waveforms, we investigate the effects resulting from the presence of a secondary polarization component. This imposes new challenges for accurate waveform reconstruction, and opens the opportunity to measure out complex polarization states such as arbitrary elliptically polarized THz field. We demonstrate a single-shot diversity-electro-optic-sampling-based approach to capture both x- and y-THz fields simultaneously with a single (110)-cut EO crystal for THz polarimetry and ellipsometry over a wide range of frequencies.

Lenz, Maximilian [University of California, Los An↗
DOE OSTI2021

Dynamic in situ imaging of methane hydrate formation in coal media

Fast 3D X-ray imaging has proved to provide crucial insights into multi-phase dynamic processes in various geomaterials. In particular, many in situ imaging experiments have been made to study gas-hydrate formation in porous sandy samples. Such imaging is challenging for the methane gas-hydrate formation in coal samples because of the coal micro-porosity structure and lower X-ray contrast. Here we present results of the first dynamic in situ micro-computed tomography experiment of methane-hydrate formation in coal samples. Synchrotron phase-contrast tomography techniques allowed to achieve necessary contrast levels to separate all the materials of interest (gas, water, coal, and gas hydrate) in reconstructed images with high spatial and temporal resolution. The imaging results are compared to the ones from a similar tomographic experiment with sand samples. Methane-hydrate formation is accompanied by the water movements caused by cryogenic water suction that happens in sequences of short fast movements with longer equilibrium states in between. Only one type of the hydrate formation was observed in coal (growth as shells on grain boundaries) as opposed to three types in sand (shells, growth into gas pockets, and inside water volumes). In particular, this leads to a slower hydrateformation speed in coal. For the coal sample, we also observed water extraction from grains, and interpreted it as competitive sorption of methane. Lastly, we visualized the dynamic behavior of this water extraction via microchannels inside the coal grains, and performed nano-tomography imaging of these channels for a better understanding of this phenomenon.

58 GEOSCIENCES↗
DOE OSTI2025

Comparing quantum annealing and spiking neuromorphic computing for sampling binary sparse coding QUBO problems

We consider the problem of computing a sparse binary representation of an image. Given an image and an overcomplete, non-orthonormal basis, we aim to find a sparse binary vector indicating the minimal set of basis vectors that when added together best reconstruct the given input. We formulate this problem with an L 2 loss on the reconstruction error, and an L 0 loss on the binary vector enforcing sparsity. First, we solve the sparse representation QUBOs by solving them both on a D-Wave quantum annealer with Pegasus chip connectivity, as well as on the Intel Loihi 2 spiking neuromorphic processor using a stochastic Non-equilibrium Boltzmann Machine (NEBM). Second, using Quantum Evolution Monte Carlo with Reverse Annealing and iterated warm starting on Loihi 2 to evolve the solution quality from the respective machines. We demonstrate that both quantum annealing and neuromorphic computing are suitable for solving binary sparse coding QUBOs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Reduced fast-ion transport calculations of m = n = 1 fishbone-like instabilities in MAST-U

Fast-ion transport associated with an m = n = 1 fishbone-like burst in MAST-U discharge 47128 is investigated using a reduced guiding-center-based transport model (ORBIT-Kick) constrained by multi-diagnostic measurements. The two-dimensional beam-emission spectroscopy system provides measurements of the core poloidal mode structure and fluctuation amplitude, while EFIT++ reconstructions constrained by the motional Stark effect diagnostic indicate a flat q-profile with q 0 > 1⁠, indicating the absence of a resonant q = 1 surface and supporting a pressure-driven infernal-mode interpretation. Analytic m = n = 1 displacement profiles consistent with the measured core mode structure and equilibrium constraints are used as the mode structure inputs to ORBIT-Kick. The calculations show that the dominant resonances occur between the mode and co-passing fast ions, producing redistribution localized near the magnetic axis. Synthetic neutron camera signals from TRANSP-Kick recover up to 90% of the experimentally observed neutron deficit at the time of peak mode amplitude, indicating that the measured m = n = 1 mode is a dominant contributor to core fast-ion transport. However, the synthetic neutron signals recover rapidly, whereas the measured neutron emission continues to decrease after the peak amplitude. In conclusion, the remaining discrepancy may arise from contributions not included in the present single-harmonic model, including higher-m and higher-n harmonics, multi-harmonic interactions, and additional transport mechanisms, motivating future diagnostic development and modeling efforts to resolve and incorporate these additional contributions.

Wong, Henry H. [University of California, Los Ange↗
DOE OSTI2024

The Reservoir Temperature Estimator (RTEst): A multicomponent geothermometry tool

The Reservoir Temperature Estimator (RTEst) is a multicomponent geothermometry tool for estimating reservoir geochemical parameters including reservoir temperature, CO 2 fugacity, mass of water lost or gained, and a reaction factor. It estimates these parameters and their associated uncertainties by minimizing an objective function that is the weighted sum of squares of the saturation indices of a user-selected set of minerals believed to be equilibrated with the reservoir fluid. RTEst accomplishes these estimates by combining the geochemical modeling capabilities of The Geochemist’s Workbench® with the optimization/parameter estimation resources of PEST®. An included interface aids the user in selecting plausible mineral phases to comprise the objective function and calculates their weighting factors. The working principles of RTEst are described and its efficacy is illustrated by presenting results of its application to various geothermal fields with known conditions. These examples show RTEst can account for the alteration of ascending reservoir fluid by mineral (calcite) re-equilibration with changes in temperature, reconstruct waters with CO 2 loss, correct for the deficit of water and other volatiles (CO 2,gas , H 2 S gas ) from boiling, and determine the amount of mixing of thermal and non-thermal waters. RTEst can use data with basis species below detection limit, missing, or unreliable either by assuming equilibrium with a controlling mineral (fixed-analyte method) or by treating the analyte concentration as an optimization parameter. The inverse of variance weighting method included in RTEst provides more representative results than either the normalization or unit weighting methods. Finally, the ability of RTEst to calculate reservoir temperatures, gas fugacity, and mixing fractions demonstrates its usefulness as a tool for evaluating geothermal systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Laser-heated diamond anvil cell synthesis and recovery of metastable MnSb 2 and YbZn 2 for post-synthesis transport studies

The creation and exploration of new materials under extreme pressure–temperature conditions has become increasingly reliant on laser-heated diamond anvil cell (LHDAC) techniques, which provide direct access to previously unexplored regions of multinary phase diagrams. Whereas numerous high-pressure phases have been identified in situ, systematic recovery and post-synthesis physical property characterization of these materials remain significant challenges. In this work, we describe the setup and implementation of an LHDAC-based synthesis and recovery workflow and demonstrate its application to metastable MnSb 2 and YbZn 2 phases. Synchrotron x-ray diffraction and spatial mapping confirm dominant formation of the targeted phases, whereas laboratory-based refinement quantifies phase fractions despite intrinsic microstrain and minor secondary phases. High-pressure transport measurements on recovered samples reveal pressure-tunable electronic instabilities in both systems. In MnSb 2 , pressure suppresses two high-temperature magnetic ordering anomalies, observed in transport, by ∼5 GPa and, for higher pressures, induces a new low-temperature feature that increases with further pressure increase. In hexagonal high-pressure YbZn 2 , an electronic reconstruction emerges at ∼11 GPa, characterized by semiconducting-like behavior from ∼30 to 300 K and a broad low-temperature coherence crossover near 30 K. Our results establish LHDAC synthesis not only as a structural discovery tool but also as an experimental platform for investigating correlated quantum states stabilized far from equilibrium thermodynamic conditions.

Huyan, S. [Iowa State University, Ames, IA (United↗
DOE OSTI2024

Plasma pressure profiles in a sheared-flow-stabilized Z-pinch

We report the plasma pressure reached inside the central plasma column of a sheared-flow-stabilized Z-pinch using Thomson scattering measurements. Building on previously reported experimental results and the analysis methods established for the high temperature and moderate density plasmas generated on the FuZE device, we show evidence of a central plasma region with higher electron temperature and density, which is consistent with a pinch behavior. Elevated electron temperatures up to 2.25 ± 0.8 keV and densities up to (4.9±0.2)×1017 cm−3 are observed to temporally coincide with the fusion neutron production from the plasma. Reconstructed plasma pressure profiles highlight the presence of a several millimeter-wide column with elevated pressure whose location varies shot-to-shot. The plasma pressure rises as neutron production from the deuterium plasma increases, reaching a peak value of 2.6 kBar. This peak value is consistent with a radially force-balanced pinch equilibrium model based on the measured ∼320 kA pinch current. Complete datasets were obtained at two axial locations, 10 and 20 cm axial position from the tip of the central electrode, which corroborate the estimated neutron source axial lengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2012

Constraints on the Parental Melts of Enriched Shergottites from Image Analysis and High Pressure Experiments

Martian basalts can be classified in at least two geochemically different families: enriched and depleted shergottites. Enriched shergottites are characterized by higher incompatible element concentrations and initial Sr-87/Sr-86 and lower initial Nd-143/Nd-144 and Hf-176/Hf-177 than depleted shergottites [e.g. 1, 2]. It is now generally admitted that shergottites result from the melting of at least two distinct mantle reservoirs [e.g. 2, 3]. Some of the olivine-phyric shergottites (either depleted or enriched), the most magnesian Martian basalts, could represent primitive melts, which are of considerable interest to constrain mantle sources. Two depleted olivine-phyric shergottites, Yamato (Y) 980459 and Northwest Africa (NWA) 5789, are in equilibrium with their most magnesian olivine (Fig. 1) and their bulk rock compositions are inferred to represent primitive melts [4, 5]. Larkman Nunatak (LAR) 06319 [3, 6, 7] and NWA 1068 [8], the most magnesian enriched basalts, have bulk Mg# that are too high to be in equilibrium with their olivine megacryst cores. Parental melt compositions have been estimated by subtracting the most magnesian olivine from the bulk rock composition, assuming that olivine megacrysts have partially accumulated [3, 9]. However, because this technique does not account for the actual petrography of these meteorites, we used image analysis to study these rocks history, reconstruct their parent magma and understand the nature of olivine megacrysts.

Collinet, M.↗
DOE OSTI2020

In Vivo Thermodynamic Analysis of Glycolysis in Clostridium thermocellum and Thermoanaerobacterium saccharolyticum Using 13 C and 2 H Tracers

Clostridium thermocellum and Thermoanaerobacterium saccharolyticum are thermophilic anaerobic bacteria with complementary metabolic capabilities that utilize distinct glycolytic pathways for the conversion of cellulosic sugars to biofuels. We integrated quantitative metabolomics with 2 H and 13 C metabolic flux analysis to investigate the in vivo reversibility and thermodynamics of the central metabolic networks of these two microbes. We found that the glycolytic pathway in C. thermocellum operates remarkably close to thermodynamic equilibrium, with an overall drop in Gibbs free energy 5-fold lower than that of T. saccharolyticum or anaerobically grown Escherichia coli . The limited thermodynamic driving force of glycolysis in C. thermocellum could be attributed in large part to the small free energy of the phosphofructokinase reaction producing fructose bisphosphate. The ethanol fermentation pathway was also substantially more reversible in C. thermocellummthan in T. saccharolyticum. These observations help explain the comparatively low ethanol titers of C. thermocellum and suggest engineering interventions that can be used to increase its ethanol productivity and glycolytic rate. In addition to thermodynamic analysis, we used our isotope tracer data to reconstruct the T. saccharolyticum central metabolic network, revealing exclusive use of the Embden-Meyerhof-Parnas (EMP) pathway for glycolysis, a bifurcated tricarboxylic acid (TCA) cycle, and a sedoheptulose bisphosphate bypass active within the pentose phosphate pathway. Thermodynamics constitutes a key determinant of flux and enzyme efficiency in metabolic networks. Here, we provide new insights into the divergent thermodynamics of the glycolytic pathways of C. thermocellum and T. saccharolyticum, two industrially relevant thermophilic bacteria whose metabolism still is not well understood. We report that while the glycolytic pathway in T. saccharolyticum is as thermodynamically favorable as that found in model organisms, such as E. coli or Saccharomyces cerevisiae , the glycolytic pathway of C. thermocellum operates near equilibrium. The use of a near-equilibrium glycolytic pathway, with potentially increased ATP yield, by this cellulolytic microbe may represent an evolutionary adaptation to growth on cellulose, but it has the drawback of being highly susceptible to product feedback inhibition. The results of this study will facilitate future engineering of high-performance strains capable of transforming cellulosic biomass to biofuels at high yields and titers.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1999

3D-Stereoscopic Analysis of Solar Active Region Loops: SoHo/EIT Observations at Temperatures of 1.5-2.5 MK - 2

In this paper we study the three-dimensional (3D) structure of hot (T(sub e) approximately equals 1.5 - 2.5 MK) loops in solar active region NOAA 7986, observed on 1996 August 30 with the Extreme-ultraviolet Imaging Telescope (EIT) onboard the Solar and Heliospheric Observatory (SoHO). This complements a first study on cooler (T(sub e) approximately equals 1.0 - 1.5 MK) loops of the same active region, using the same method of Dynamic Stereoscopy to reconstruct the 3D geometry. We reconstruct the 3D-coordinates x(s), y(s), z(s), the density n(sub e)(s), and temperature profile T(sub e)(s) of 35 individual loop segments (as function of the loop coordinate s) using EIT 195 A and 284 A images. The major findings are: (1) All loops are found to be in hydrostatic equilibrium, in the entire temperature regime of T(sub e) = 1.0 - 2.5 MK; (2) The analyzed loops have a height of 2-3 scale heights, and thus only segments extending over about one vertical scale height have sufficient emission measure contrast for detection; (3) The temperature gradient over the lowest scale height is of order dT/ds is approximately 1 - 4 K/km; (4) The radiative loss rate is found to exceed the conductive loss rate by about two orders or magnitude, making thermal conduction negligible to explain the temperature structure of the loops; (5) A steady-state can only be achieved when the heating rate E(sub H) matches the radiative loss rate in hydrostatic equilibrium, requiring a heat deposition length lambda(sub H) of the half density scale height lambda, predicting a scaling law with the loop base pressure, EH varies as p(sub 0 exp 2). This favors coronal heating mechanisms that operate near the loop footpoints; (6) We find a reciprocal correlation between the loop pressure p(sub 0) and loop length L, i.e. p(sub 0) varies as 1/L, implying a scaling law of the steady-state requirement with loop length, i.e. E(sub H ) varies as 1/L(exp 2). The heating rate shows no correlation with the loop-aligned magnetic field component B(sub z) at the footpoints, but is correlated with the azimuthal field B(sub phi) = Bz(RDelta Phi/L) of a twisted loop, and is thus consistent with heating mechanisms based on field-aligned currents.

Aschwanden, Markus J.↗
DOE OSTI2022

Local measures of entanglement in black holes and CFTs

We study the structure and dynamics of entanglement in CFTs and black holes. We use a local entanglement measure, the entanglement contour, which is a spatial density function for von Neumann entropy with some additional properties. The entanglement contour can be calculated in many 1+1d condensed matter systems and simple models of black hole evaporation. We calculate the entanglement contour of a state excited by a splitting quench, and find universal results for the entanglement contours of low energy non-equilibrium states in 2d CFTs. We also calculate the contour of a non-gravitational bath coupled to an extremal AdS _2 2 black hole, and find that the contour only has finite support within the bath, due to an island phase transition. The particular entanglement contour proposal we use quantifies how well the bath’s state can be reconstructed from its marginals, through its connection to conditional mutual information, and the vanishing contour is a reflection of the protection of bulk island regions against erasures of the boundary state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Edge ICRF simulations in 3D geometry: From MHD equilibrium to coupling determination

We present in this work a consistent numerical scheme that allows the computation of 3D magnetic fields a nd 3D density profiles and their usage in ion cyclotron range of frequencies (ICRF) coupling simulations. We first utilize the PARVMEC code to compute the 3D free-boundary plasma equilibrium in the ideal magnetohydrodynamic (MHD) approximation. Since the PARVMEC solution is only defined within the last closed flux surface (LCFS), the magnetic field domain is extended to the scrape-off layer (SOL) via the BMW code, which computes a divergence-free magnetic field solution arising from the external conductors’ vacuum field a nd t he PARVMEC flux surface currents. This magnetic reconstruction is then used in the EMC3-EIRENE transport code in order to compute 3D density profiles. In the last step, the RAPLICASOL code is utilized to compute the ICRF antenna S-matrices resulting from the 3D density profiles. We exemplify this scheme for the ASDEX Upgrade tokamak. A new implementation of a curved model for the ASDEX Upgrade ICRF 2-strap antenna in RAPLICASOL allows simulations in realistic geometry, without any coordinate transformations.

Lopez, Graham↗
DOE OSTI2023

Isotopic Constraints on the Nature of Primary Precipitates in Archean–Early Paleoproterozoic Iron Formations from Determinations of the Iron Phonon Density of States of Greenalite and 2L- and 6L-Ferrihydrite

Iron formations (IFs) are chemical sedimentary rocks that were widely deposited before the Great Oxidation Event (GOE) around 2.4- 2.2 Ga. It is generally thought that IFs precipitated as hydrated Fe 3+ oxides (HFOs) such as ferrihydrite following surface oxidation of Fe 2+ -rich, anoxic deep waters. This model often implicates biological oxidation and underpins reconstructions of marine nutrient concentrations. However, nanoscale petrography indicates that an Fe 2+ silicate, greenalite, is a common primary mineral in well-preserved IFs, motivating an alternative depositional model of anoxic ferrous silicate precipitation. It is unclear, however, if Fe 2+ -rich silicates can produce the Fe isotopic variations in IFs that are well explained by Fe 2+ oxidation. To address this question, we constrain the equilibrium Fe isotopic ( 56 Fe/ 54 Fe) fractionation of greenalite and ferrihydrite by determining the iron phonon densities of states for those minerals. Here, we use ab initio density functional theory (DFT + U) calculations and nuclear resonant inelastic X-ray scattering spectroscopy to show that ferrous greenalite should be isotopically lighter than ferrihydrite by similar to 1-1.2 parts per thousand at equilibrium, and fractionation should scale linearly with increasing Fe 3+ content in greenalite. By anchoring ferrihydrite-greenalite mineral pair fractionations to published experimental Fe isotopic fractionations between HFOs and aqueous Fe 2+ , we show that ferrous greenalite may produce all but the heaviest pre-GOE Fe isotopic compositions and mixed valence greenalites can produce the entire record. Our results suggest that heavy Fe isotope enrichments alone are not diagnostic of primary IF mineralogies, and ferrihydrite and partially oxidized or even purely ferrous greenalite are all viable primary IF mineralogies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Quasiparticle Poisoning of Superconducting Qubits from Resonant Absorption of Pair-Breaking Photons

The ideal superconductor provides a pristine environment for the delicate states of a quantum computer: because there is an energy gap to excitations, there are no spurious modes with which the qubits can interact, causing irreversible decay of the quantum state. As a practical matter, however, there exists a high density of excitations out of the superconducting ground state even at ultralow temperature; these are known as quasiparticles. Observed quasiparticle densities are of order 1 μ m - 3 , tens of orders of magnitude greater than the equilibrium density expected from theory. Nonequilibrium quasiparticles extract energy from the qubit mode and can induce dephasing. Here we show that a dominant mechanism for quasiparticle poisoning is direct absorption of high-energy photons at the qubit junction. We use a Josephson junction-based photon source to controllably dose qubit circuits with millimeter-wave radiation, and we use an interferometric quantum gate sequence to reconstruct the charge parity of the qubit. We find that the structure of the qubit itself acts as a resonant antenna for millimeter-wave radiation, providing an efficient path for photons to generate quasiparticles. A deep understanding of this physics will pave the way to realization of next-generation superconducting qubits that are robust against quasiparticle poisoning.

Liu, Chuan-Hong↗