Search NASA⌕ Search

SEARCH · Search NASA

Results for “Equilibrium programming”

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

Crystalline Silicotitanate (CST) Ion Exchange Media Performance Evaluations to Support TSCR DSA IX Media Equilibrium Contacts

The primary objective of this work is to calculate the maximum loading expected on the Hanford Tank Side Cesium Removal (TSCR) ion exchange columns. A key consideration in the design of the columns is the amount of 137 Cs that loads onto the Crystalline Silicotitanate (CST) and the heat generated by the loaded column during storage. Per request of Washington River Protection Solutions (WRPS), Savannah River National Laboratory (SRNL) has utilized ZAM, a computer program developed by the research group of Professor Rayford G. Anthony of Texas A&M University, to predict the cesium loading on the CST for a variety of waste compositions expected to be processed by TSCR. The study evaluated cesium loadings for the following waste compositions: 1. Seventeen DFLAW campaign batches to cover projected supernate composition ranges within which TSCR may be expected to operate within the first ten years, 2. Hanford tank AP-105 and AP-107 waste solutions that will be processed by the TSCR system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Thermal-hydraulic analyses of MIT reactor LEU transition cycles

The Massachusetts Institute of Technology Reactor (MITR) is a research reactor in Cambridge, Massachusetts designed primarily for experiments using neutron beam and in-core irradiation facilities. In the framework of non-proliferation policy, research and test reactors have started a program to convert highly enriched uranium (HEU) fuel to low-enriched uranium (LEU) fuel. A new type of LEU fuel based on a high-density alloy of uranium and molybdenum (U-10Mo) is expected to allow conversion. A transition core plan, from 22 fresh LEU fuel elements to 24 in an equilibrium state, is proposed. This study evaluates the thermal-hydraulic safety margins of each transition cycle and state (i.e. beginning, middle and end). The STAT7 and RELAP5 codes are used in the study. STAT7 provides an integrated platform for statistical propagation of uncertainties in determining steady-state thermal-hydraulic margins. A RELAP5 model was created to verify MITR's STAT7 model, and provide transient analyses. The steady-state analyses determined that, during the entire transition from a fresh to equilibrium core, the Limiting System Safety Settings (LSSS) has a 10.2% margin to the minimum power at which onset of nucleate boiling (ONB) is precluded at 3σ confidence. The loss of flow (LOF) transient, which is the most limiting anticipated transient for the MITR, was simulated. With a conservative initial state of the reactor (operating at μ + 3σ power level and LSSS mass flow rate), the maximum calculated centerline fuel temperature was 115 degrees C, which is a significantly lower value than the blistering limit of 350 degrees C. The cladding wall temperatures did not exceed the ONB temperature throughout the transient. Therefore, during LOF, no nucleate boiling is expected, precluding any critical heat flux trigger. Overall, this study indicates that the proposed LEU transition core specifications have significant margin to thermal-hydraulic limits, during steady-state operation and the LOF transient.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Atmospheric CO 2 estimates for the Miocene to Pleistocene based on foraminiferal $δ^{11}$B at Ocean Drilling Program Sites 806 and 807 in the Western Equatorial Pacific

Constraints on the evolution of atmospheric CO 2 levels throughout Earth's history are foundational to our understanding of past variations in climate. Despite considerable effort, records vary in their temporal and spatial coverage and estimates of past CO 2 levels do not always converge, and therefore new records and proxies are valuable. Here we reconstruct atmospheric CO 2 values across major climate transitions over the past 16 million years using the boron isotopic composition (δ 11 B) of planktic foraminifera from 89 samples obtained from two sites in the West Pacific Warm Pool, Ocean Drilling Program (ODP) Sites 806 and 807, measured using high-precision multi-collector inductively coupled plasma mass spectrometry. We compare our results to published data from ODP Site 872, also in the Western Equatorial Pacific, that goes back to 22 million years ago. These sites are in a region that today is near equilibrium with the atmosphere and are thought to have been in equilibrium with the atmosphere for the interval studied. We show that δ 11 B data from this region are consistent with other boron-based studies. The data show evidence for elevated pCO 2 during the Middle Miocene and Early to Middle Pliocene, and reductions in pCO 2 of ~200 ppm during the Middle Miocene Climate Transition, ~250 ppm during Pliocene Glacial Intensification and ~50 ppm during the Mid-Pleistocene Climate Transition. During the Mid-Pleistocene Transition there is a minimum pCO 2 at marine isotopic stage (MIS) 30. Our results are consistent with a coupling between pCO 2 , temperature and ice sheet expansion from the Miocene to the late Quaternary.

54 ENVIRONMENTAL SCIENCES↗

Phosphoproteome Response to Dithiothreitol Reveals Unique Versus Shared Features of Saccharomyces cerevisiae Stress Responses

To cope with sudden changes in the external environment, the budding yeast Saccharomyces cerevisiae orchestrates a multifaceted response that spans many levels of physiology. Several studies have interrogated the transcriptome response to endoplasmic reticulum (ER) stress and the role of regulators such as the Ire1 kinase and Hac1 transcription factors. However, less is known about responses to ER stress at other levels of physiology. Here, we used quantitative phosphoproteomics and computational network inference to uncover the yeast phosphoproteome response to the reducing agent dithiothreitol (DTT) and the upstream signaling network that controls it. We profiled wild-type cells and mutants lacking IRE1 or MAPK kinases MKK1 and MKK2, before and at various times after DTT treatment. In addition to revealing downstream targets of these kinases, our inference approach predicted new regulators in the DTT response, including cell-cycle regulator Cdc28 and osmotic-response kinase Rck2, which we validated computationally. Our results also revealed similarities and surprising differences in responses to different stress conditions, especially in the response of protein kinase A targets. Furthermore, these results have implications for the breadth of signaling programs that can give rise to common stress response signatures.

59 BASIC BIOLOGICAL SCIENCES↗

Characteristics of Ultralow Aspect Ratio Toroidal Plasmas and Non-Solenoidal Startup and Edge Stability Studies at Near-Unity Aspect Ratio (Final Report)

This grant supported a long-term research and development program to explore the physics and technology of spherical tokamak (ST) plasma behavior at near-unity aspect ratio ( A ). The PEGASUS Toroidal Experiment is a university-scale ( R ) ~ 0.40 m, α ~ .35 m, I p ≤ 0.3 MA, B T ≤ 0.15 T) experiment that accesses A → 1, to provide large values of the normalized plasma current, I N = I p /α B T , or equivalently the toroidal field utilization factor I p / I TF . The research activities focused on: 1) developing the technique of Local Helicity Injection (LHI) as a viable means of initiating and growing an ST plasma without the use of a central solenoid; 2) exploring the equilibrium and stability properties of spherical tokamak plasmas at extremely low aspect ratio; and 3) obtaining new insights in the properties of the enhanced confinement H-mode regime at A ~ 1. The initial (Phase I) facility employed a unique high-stress central solenoid Ohmic for current drive and heating with simple magnet current waveforms. A subsequent major upgrade (Phase II) offered more flexible control of the plasma evolution along with a powerful LHI system for plasma initiation, heating and sustainment. The LHI technique consists of injection of high-density current streams into vacuum toroidal and vertical magnetic fields. These streams become unstable and support magnetic reconnection to allow the system to relax into a tokamak-like magnetic configuration. This program focused on: 1) development of electron current sources that provide sufficiently high injection currents in the challenging edge plasma region; 2) validation of theoretical limits to achievable I p arising from magnetic helicity and energy conservation plus absolute limits arising from relaxation to a minimal energy states; 3) tests of a power balance model to describe the evolution of the plasma current I p (t); and 4) study of magnetic fluctuations and their role in the I p evolution. This culminated in the production of toroidal currents of I p ≥ 0.2 MA with only ~8 kA of injected current. Operation with Ohmic induction alone in Phase I created plasmas with soft limits on the achievable I p and field utilization ( I p / I TF ~ 1) due to the onset of large low-order (2/1 or 3/2) internal resonant tearing modes. With the expanded capabilities in Phase II operation, methods employed to mitigate those modes and achieve I p / I TF > 1 were: 1) I TF ramp-downs; 2) non-solenoidal startup via LHI; and 3) extremely low I TF Ohmic operations assisted by LHI current injector pre-ionization. Extensive use of LHI-driven startup and drive for non-solenoidal operation produced favorable current profiles to stably operate tokamak plasmas with unprecedentedly high I N = I p / aB T ~ 15 and I p / I TF > 2. With strong reconnection-driven anomalous ion heating, this in turn provided access to tokamak plasmas with B t up to 100%, a world record 2 times higher than that previously achieved in advanced or spherical tokamaks. This high B t regime had a large fraction (~50%) of the plasma volume in a region of an absolute minimum- B configuration. The MHD stability limit in this regime was the ideal external kink mode without a stabilizing wall. These results provided the first demonstration of access to tokamak plasmas with very high B t , high i N , high elongation κ, and low internal inductance ℓ i , and as such served as a proof-ofprinciple for access to this long-sought ST regime. The naturally low B T needed for stability at A ≤ 1.2 reduces the L-H mode power threshold, P LH , so that H-mode plasmas were achieved with Ohmic heating alone. The power threshold for the transition to the H-mode regime was found to be 10× the accepted values from the ITER international scaling, confirming initial trends indicated by experiments on MAST and NSTX, and suggesting that the physics behind this transition is not yet fully understood. The density dependence of the transition and the observed insensitivity to whether the plasma is bounded by a limiter or magnetic divertor are consistent with the recent FM 3 model. Multiple unstable toroidal magnetic modes were measured during ELM (Edge Localized Mode) crashes, with generally lower toroidal mode numbers than seen at higher A , which was attributed to the larger peeling mode instability drive at low A . First-ever edge current profile measurements through an ELM crash showed current-carrying filaments were generated, as seen in ELM simulations and in electromagnetic blob transport theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

From graphical to model-based distillation column design: A McCabe-Thiele-inspired mathematical programming approach

We propose a mixed-integer nonlinear programming (MINLP) model for simple and complex distillation column design and optimization. The model is based upon the concepts and equations underpinning the McCabe-Thiele method. Generalizing this method, we introduce material balances at various locations of the column and employ binary variables to determine the optimal number of trays and optimal feed locations. We model the vapor–liquid equilibrium using continuous piecewise linear approximating functions. The model is extended to account for multicomponent mixtures and non-constant-molar overflow. We also discuss how to estimate the minimum number of trays and the minimum reflux ratio.

42 ENGINEERING↗

Thermally induced transformations of Au@Cu 2 O core–shell nanoparticles into Au–Cu nanoparticles from temperature-programmed in situ powder X-ray diffraction

Temperature-programmed in situ X-ray diffraction with whole-powder-pattern modeling is used to investigate the reaction of Au@Cu 2 O core–shell nanoparticles to form nanocrystalline bimetallic Cu x Au 1− x alloys ( x = 0, 0.25, 0.5, 0.75, 1.0) in a reducing atmosphere. The mechanisms of the reactions are key to informed design of tailored non-equilibrium nanostructures for catalytic and plasmonic materials. The Au@Cu 2 O reaction is initiated by reduction of the Cu 2 O cuprite shell to form nanocrystalline metallic Cu at about 413 K. Alloying begins immediately upon formation of metallic Cu at 413 K, with the nucleation of an Au-rich alloy phase which reaches the nominal Cu content of the overall system stoichiometry by 493 K. All bimetallic alloys form a transient ordered Cu 3 Au intermetallic compound at intermediate temperatures, with the onset of ordering and subsequent disordering varying by composition. No evidence for an ordered Au 3 Cu intermetallic is found for any composition. Significant crystal growth in the bimetallic phase is apparent at higher temperatures, with the onset temperature increasing with Cu concentration and initial Cu-shell thickness. The reduction of the cuprite phase is slowed by the presence of the core–shell interface, and crystal growth in the Cu shell is completely suppressed within the alloy systems.

Chemistry↗

Accelerating discoveries at DIII-D with the Integrated Research Infrastructure

DIII-D research is being accelerated by leveraging high performance computing (HPC) and data resources available through the National Energy Research Scientific Computing Center (NERSC) Superfacility initiative. As part of this initiative, a high-resolution, fully automated, whole discharge kinetic equilibrium reconstruction workflow was developed that runs at the NERSC for most DIII-D shots in under 20 min. This has eliminated a long-standing research barrier and opened the door to more sophisticated analyses, including plasma transport and stability. These capabilities would benefit from being automated and executed within the larger Department of Energy Advanced Scientific Computing Research program’s Integrated Research Infrastructure (IRI) framework. The goal of IRI is to empower researchers to meld DOE’s world-class research tools, infrastructure, and user facilities seamlessly and securely in novel ways to radically accelerate discovery and innovation. For transport, we are looking at producing flux matched profiles and also using particle tracing to predict fast ion heat deposition from neutral beam injection before a shot takes place. Our starting point for evaluating plasma stability focuses on the pedestal limits that must be navigated to achieve better confinement. This information is meant to help operators run more effective experiments, so it needs to be available rapidly inside the DIII-D control room. So far this has been achieved by ensuring the data is available with existing tools, but as more novel results are produced new visualization tools must be developed. In addition, all of the high-quality data we have generated has been collected into databases that can unlock even deeper insights. This has already been leveraged for model and code validation studies as well as for developing AI/ML surrogates. The workflows developed for this project are intended to serve as prototypes that can be replicated on other experiments and can be run to provide timely and essential information for ITER, as well as next stage fusion power plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of Computational Tools for Modeling Thermal and Radiation Effects on Grain Boundary Segregation and Precipitation in Fe-Cr-Ni-based Alloys

This work aims at developing computational tools for modeling thermal and radiation effects on solute segregation at grain boundaries (GBs) and precipitation. This report described two major efforts. One is the development of computational tools on integrated modeling of thermal equilibrium segregation (TES) and radiation-induced segregation (RIS), from which synergistic effects of thermal and radiation, pre-existing GB segregation have been taken into consideration. This integrated modeling was used in describing the Cr and Ni segregation in the Fe-Cr-Ni alloys. The other effort is thermodynamic modeling on the Fe-Cr-Ni-Mo system which includes the major alloying elements in the investigated alloys in the Advanced Radiation Resistant Materials (ARRM) program. Through thermodynamic calculation, we provide baseline thermodynamic stability of the hardening phase Ni2(Cr,Mo) in selected Ni-based super alloys, and contribute knowledge on mechanistic understanding on the formation of Ni2(Cr,Mo) in the irradiated materials. The major outcomes from this work are listed in the following: 1) Under the simultaneous thermal and irradiation conditions, radiation-induced segregation played a dominant role in the GB segregation. The pre-existing GB segregation only affects the subsequent radiation-induced segregation in the short time. For the same element, the segregation tendency of Cr and Ni due to TES is opposite to it from RIS. The opposite tendency can lead to the formation of W-shape profile. These findings are consistent with literature observation of the transitory W-shape profile. 2) While TES only affects the distance of one or two atomic layers from GBs, the RIS can affect a broader distance from GB. Therefore, the W-shape due to pre-existing GB segregation is much narrower than that due to composition gradient formed during the transient state. Considering the measurement resolution of Auger or STEM analysis, the segregation tendency due to RIS should play a dominant role in the measured values. However, The GB segregation due to pre-existing GB segregation may affect the chemical potential of element at GB, and subsequently the corrosion resistance. 3) Based on the newly developed thermodynamic database of Fe-Cr-Ni-Mo, we predicted the Ni2(Cr,Mo) as a thermodynamically stable phase in all investigated low Fe-content Ni-based alloys. The calculated phase amount decreases with the increasing Fe content, being consistent with that observed in the irradiated materials. 4) The formation of the Ni2(Cr,Mo) phase in irradiated materials is due to irradiation enhanced diffusion. The calculated equilibrium Ni2(Cr,Mo) amount is more than that observed in the irradiated materials, suggesting that the amount of Ni2(Cr,Mo) is likely to increase more with further irradiation.

36 MATERIALS SCIENCE↗

The Double Galaxy Cluster A2465. III. X-Ray and Weak-lensing Observations

We report Chandra X-ray observations and optical weak-lensing measurements from Subaru/Suprime-Cam images of the double galaxy cluster A2465 ( z = 0.245). The X-ray brightness data are fit to a β model to obtain the radial gas density profiles of the northeast (NE) and southwest (SW) subcomponents, which are seen to differ in structure. We determine core radii, central temperatures, the gas masses within r {sub 500c}, and the total masses for the broader NE and sharper SW components assuming hydrostatic equilibrium. There is no large X-ray excess between the two components. The central entropy of the NE subcluster is about two times higher than the SW. Along with its structural properties and an apparent radio halo that is a sign of a merger, this suggests that the NE component has undergone merging on its own. The weak-lensing analysis gives virial masses for each substructure, which compare well with earlier dynamical results. The derived outer mass contours of the SW sub-component from weak lensing are more irregular and extended than those of the NE. Although there is a weak enhancement and small offsets between X-ray gas and mass centers from weak lensing, the lack of large amounts of gas between the two subclusters indicates that A2465 is in a pre-merger state. We discuss star formation enhancement in this system resulting from its dynamics and shock-induced star formation scenarios. A dynamical model that is consistent with the observed cluster data, based on the FLASH program and the radial infall model, is constructed, where the subclusters currently separated by ∼1.2 Mpc are approaching each other at ∼2000 km s{sup −1} and will meet in ∼0.4 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗

TOI-1231 b: A Temperate, Neptune-sized Planet Transiting the Nearby M3 Dwarf NLTT 24399

We report the discovery of a transiting, temperate, Neptune-sized exoplanet orbiting the nearby (d = 27.5 pc), M3V star TOI-1231 (NLTT 24399, L 248-27, 2MASS J10265947-5228099). The planet was detected using photometric data from the Transiting Exoplanet Survey Satellite and followed up with observations from the Las Cumbres Observatory and the Antarctica Search for Transiting ExoPlanets program. Combining the photometric data sets, we find that the newly discovered planet has a radius of 3.65{sub −0.15}{sup +0.16} R{sub ⊕} and an orbital period of 24.246 days. Radial velocity measurements obtained with the Planet Finder Spectrograph on the Magellan Clay telescope confirm the existence of the planet and lead to a mass measurement of 15.5 ± 3.3 M {sub ⊕}. With an equilibrium temperature of just 330 K, TOI-1231 b is one of the coolest small planets accessible for atmospheric studies thus far, and its host star’s bright near-infrared brightness (J = 8.88, K {sub s} = 8.07) makes it an exciting target for the Hubble Space Telescope and the James Webb Space Telescope. Future atmospheric observations would enable the first comparative planetology efforts in the 250–350 K temperature regime via comparisons with K2-18 b. Furthermore, TOI-1231's high systemic radial velocity (70.5 km s{sup −1}) may allow for the detection of low-velocity hydrogen atoms escaping the planet by Doppler, shifting the H i Lyα stellar emission away from the geocoronal and interstellar medium absorption features.

47 OTHER INSTRUMENTATION↗

The Hubble Space Telescope PanCET Program: An Optical to Infrared Transmission Spectrum of HAT-P-32Ab

We present a 0.3−5 μm transmission spectrum of the hot Jupiter HAT-P-32Ab observed with the Space Telescope Imaging Spectrograph and Wide Field Camera 3 instruments mounted on the Hubble Space Telescope, combined with Spitzer Infrared Array Camera photometry. The spectrum is composed of 51 spectrophotometric bins with widths ranging between 150 and 400 Å, measured to a median precision of 215 ppm. Comparisons of the observed transmission spectrum to a grid of 1D radiative-convective equilibrium models indicate the presence of clouds/hazes, consistent with previous transit observations and secondary eclipse measurements. To provide more robust constraints on the planet’s atmospheric properties, we perform the first full optical to infrared retrieval analysis for this planet. The retrieved spectrum is consistent with a limb temperature of 1248{sub −92}{sup +92} K, a thick cloud deck, enhanced Rayleigh scattering, and ∼10× solar H{sub 2}O abundance. We find log(Z/Z {sub ⊙}) = 2.41{sub −0.07}{sup +0.06}, and compare this measurement with the mass–metallicity relation derived for the solar system.

79 ASTRONOMY AND ASTROPHYSICS↗

Probing direct carbon-carbon acylation of furans and long-chain acids over H-ZSM-5

The acylation of furan and 2-methylfuran (2-MF) with acetic acid (AA) and 1-hexanoic acid (HA) was studied over an H-ZSM-5 catalyst. Temperature-programmed desorption (TPD) measurements showed evidence that AA adsorbs strongly at Brønsted sites, with reaction to CO 2 and acetone beginning at 523 K, while both furan and 2-MF oligomerize upon adsorption at room temperature. Steady-state, flow-reactor measurements demonstrate that selectivities for the acylation reactions can approach 100% below 523 K when reactions are operated in excess organic acids, with maximum conversion limited by equilibrium. Coking becomes important when the reaction is operated in excess furan or 2-MF and at high reactant concentrations, while ketonization products are observed when the reaction is performed above 523 K. Lastly, potential uses for the reactions are discussed.

1-Hexanoic Acid↗

Multi-Physics Simulations for Molten Salt Reactor Evaluation: Chemistry Modeling and Database Development

To aid in design and licensing of molten salt reactors, a framework integrating the complex interaction of reactor neutronics, thermal hydraulics, and chemistry is being developed within the Department of Energy Advanced Reactor Technology Program’s Molten Salt Reactor (MSR) campaign. The challenges of integrating thermochemical and thermophysical behavior into a multi-physics reactor simulations include the following: (1) population of data needed for refinement of current models and development of nonexistent models through experimental measurements, first principles calculations, and development of a machine learning approach (2) thermochemical and thermophysical model development, (3) further development of Thermochimica, an open-source efficient equilibrium solver used to link thermochemical models to the multi-physics code, (4) a framework for integrating kinetic phenomena: nucleation, precipitation, mass/heat transport, and corrosion models, and (5) a computational environment to efficiently utilize the data and models within a multi-physics modeling tool. These challenges are being addressed through a collaboration among Oak Ridge National Laboratory, Argonne National Laboratory, the University of South Carolina, and the University of Ontario Institute of Technology.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

AdvEP

AdvEP is a code repository which contains PyTorch implementations of various adversarial attacks on a deep neural network trained with Equilibrium Propagation (EP), which is a neuromorphic learning framework. AdvEP allows for the training, testing, and conducting white/black-box attacks of EP models on a wide variety of applications and datasets. AdvEP is based on the open-source code https://github.com/Laborieux-Axel/Equilibrium-Propagation which was developed to train energy models. AdvEP was created by modifying the original code to perform and test against adversarial attacks. AdvEP was developed in Python, a high-level programming language that takes advantage of the Python ecosystem of high-quality open-source packages for machine learning. AdvEP interfaces heavily with the open-source PyTorch Python package as well as the open-source Adversarial Robustness Toolbox (ART) package.

Mansingh, Siddarth↗

Geochemical Changes in Response to CO 2 Injection in a CO 2 -EOR Complex in Northern Michigan (Volume II.1)

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. The report describes the geochemical monitoring program under MRCSP to use stable and radiogenic isotope geochemistry in concert with analysis of general geochemical parameters of fluids and gases and analysis of core samples to determine geochemical processes occurring in the reef structure because of CO 2 injection. Specifically, brine and gas samples were collected and analyzed to determine changes occurring between reefs prior to and following CO 2 injection. The analytical results for general geochemical parameters were modeled with chemical equilibrium models to determine if the injection of CO 2 resulted in the mineral dissolution or precipitation. Finally, core samples were collected and analyzed to determine if there was evidence of dissolution features or mineral precipitation.

01 COAL, LIGNITE, AND PEAT↗

Quantum simulation of charge and exciton transfer in multi-mode models using engineered reservoirs

Quantum simulation enables studies of open-system dynamics in non-perturbative regimes by programming electronic, vibrational, and environmental interactions on comparable energy scales. Trapped ions offer this capability, combining spins, phonons, and tunable dissipation on one platform. We demonstrate an open-system quantum simulation of charge and exciton transfer in a multi-mode linear vibronic coupling model. Using tailored spin-phonon interactions with reservoir engineering, we emulate a system with two dissipative vibrational modes coupled to donor and acceptor sites and track its non-equilibrium dynamics. We continuously tune the system from the charge transfer regime to the vibrationally assisted exciton transfer regime and find that degenerate modes enhance transfer rates at large energy gaps, while non-degenerate modes activate pathways that reduce the energy-gap dependence. Thus, the presence of one additional vibration introduces interfering pathways and reshapes non-perturbative excitation transfer. Our results establish a scalable, hardware-efficient route to simulate vibronic processes with engineered environments.

74 ATOMIC AND MOLECULAR PHYSICS↗

The General Atomics Center for Magnetohydrodynamic Modeling (CEMM) Final Report

The Center for Extended Magnetohydrodynamic Modeling (CEMM) is built around the further development, verification, validation, and other application of two major extended magnetohydrodynamics (MHD) codes NIMROD and M3D-C1. The General Atomics (GA) CEMM Program focuses on three main research areas: extension of the M3D-C1 MHD code plasma response calculations to new regimes, coupling of the DK4D drift kinetic equation solver to a temperature evolution equation, and simulations of disruption mitigation for unstable plasmas with pre-existing magnetic islands using massive gas injection (MGI) with the NIMROD MHD code. A systematic scan of the zero-crossing of the rotation profile in a DIII-D ITER-similar shape equilibrium using linear, time-independent modeling with M3D-C1 was performed. A new Python utility, autoC1, has been developed to automate the running of M3D-C1 for linear simulations. A number of ELM suppression experiments were analyzed including a validation and cross-code verification with the MARS-F MHD code of magnetic signal phase shifts across ELM-suppression bifurcations. A linearized, time-dependent temperature equation was tightly coupled to the electron and ion drift-kinetic equations (DKEs) in the DK4D drift-kinetic equation solver. Various explicit and implicit implementations were attempted, but all methods proved numerically unstable, even with the addition of numerical dissipation. Several simulations were completed investigating the effects of pre-existing islands on disruption mitigation by MGI in DIII-D using NIMROD.

3D MHD Simulation↗