Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time Evolution”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Fourier-MIONet: Fourier-enhanced multiple-input neural operators for multiphase modeling of geological carbon sequestration

Geologic carbon sequestration (GCS) is a safety-critical technology that aims to reduce the amount of carbon dioxide in the atmosphere, which also places high demands on reliability. Multiphase flow in porous media is essential to understand CO 2 migration and pressure fields in the subsurface associated with GCS. However, numerical simulation for such problems in 4D is computationally challenging and expensive, due to the multiphysics and multiscale nature of the highly nonlinear governing partial differential equations (PDEs). It prevents us from considering multiple subsurface scenarios and conducting real-time optimization. Here, we develop a Fourier-enhanced multiple-input neural operator (Fourier-MIONet) to learn the solution operator of the problem of multiphase flow in porous media. Fourier-MIONet utilizes the recently developed framework of the multiple-input deep neural operators (MIONet) and incorporates the Fourier neural operator (FNO) in the network architecture. Once Fourier-MIONet is trained, it can predict the evolution of saturation and pressure of the multiphase flow under various reservoir conditions, such as permeability and porosity heterogeneity, anisotropy, injection configurations, and multiphase flow properties. Compared to the enhanced FNO (U-FNO), the proposed Fourier-MIONet has 90% fewer unknown parameters, and it can be trained in significantly less time (about 3.5 times faster) with much lower CPU memory (<15%) and GPU memory (<35%) requirements, to achieve similar prediction accuracy. In addition to the lower computational cost, Fourier-MIONet can be trained with only 6 snapshots of time to predict the PDE solutions for 30 years. Furthermore, we observed that Fourier-MIONet can maintain good accuracy when predicting out-of-distribution (OOD) data. The excellent generalizability of Fourier-MIONet is enabled by its adherence to the physical principle that the solution to a PDE is continuous over time. Furthermore, the developed Fourier-MIONet makes it possible to solve the long-time evolution of geological carbon sequestration in a large-scale three-dimensional space accurately and efficiently.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Chemical decomposition in shock-compressed 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals: Time-resolved Raman spectroscopy

Understanding the evolution of chemical decomposition—at the molecular level—in shock-compressed insensitive high explosive (IHE) single crystals is an important need. Toward this goal, time-resolved Raman spectra were measured in 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals shock compressed to stresses above 20 GPa. At 22 GPa and higher stresses, the Raman peak intensities were significantly reduced, showing the onset and extent of chemical changes. At 33 GPa, no Raman peaks were observed, suggesting complete decomposition of FOX-7. The 22 GPa onset from Raman data is in marked contrast to the much higher (32 GPa) onset determined previously from wave profile measurements—showing the significantly greater sensitivity of Raman measurements for gaining insight into shock-induced molecular-level changes in IHE single crystals. The close match between the complete loss of Raman peaks at 33 GPa and the 32 GPa reaction threshold determined from continuum data makes a good case that the two results are related and they reflect significant energy release due to decomposition.

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

Optical emission spectroscopy and imaging of low-pressure N2 plasmas generated by intense fast-pulsed electron beams

An optical emission spectroscopic (OES) and imaging characterization is conducted on N2 plasmas generated by a 100 keV fast-pulsed electron beam. The electron beams are injected into an N2 gas filled volume with a current of 4.5 kA (300 A/cm2) and a 100 ns pulse width. The characterization is conducted at the pressures, 1 Torr and 0.1 Torr, corresponding to two distinct regimes that exhibit significantly different plasma dynamics. Beam impact ionization is shown to be a primary mechanism for producing low temperature plasmas at 1 Torr during beam output. After beam termination, ionization by an inductive electric field becomes the primary mechanism for plasma formation later in time for both pressures. OES and plasma imaging are used in this work as a diagnostic tool to track the distribution of electronic, vibrational, and rotational state transitions and ionized species. This is achieved with the use of a multi-resolution suite of spectrometers capable of acquiring time-resolved spectra. Vibrational and rotational bands of the N2 second positive system (C3Πu→B3Πg) and the N2+ first negative system (B2Σu+→X2Σg+) are identified in both regimes as well as N+ states exclusively in the lower pressure regime. Vibrational and rotational spectra are shown to track the evolution of the time-varying plasma current. Plasma imaging also reveals spatially nonuniform plasma emission at 0.1 Torr. A few hundred shots are recorded to fully characterize the emissions, and results are shown to be highly reproducible (≤±1% with 2σ confidence).

Kaiser, E. R. (ORCID:0000000336493938)↗
DOE OSTI2018

Tracing Phase Transformation and Lattice Evolution in a TRIP Sheet Steel under High-Temperature Annealing by Real-Time In Situ Neutron Diffraction

Real-time in situ neutron diffraction was used to characterize the crystal structure evolution in a transformation-induced plasticity (TRIP) sheet steel during annealing up to 1000 °C and then cooling to 60 °C. Based on the results of full-pattern Rietveld refinement, critical temperature regions were determined in which the transformations of retained austenite to ferrite and ferrite to high-temperature austenite during heating and the transformation of austenite to ferrite during cooling occurred, respectively. The phase-specific lattice variation with temperature was further analyzed to comprehensively understand the role of carbon diffusion in accordance with phase transformation, which also shed light on the determination of internal stress in retained austenite. These results prove the technique of real-time in situ neutron diffraction as a powerful tool for heat treatment design of novel metallic materials.

Yu, Dunji [ORNL] (ORCID:0000000189467851)↗
DOE OSTI2024

Skeletal reaction models for methane combustion

A local-sensitivity-analysis technique is employed to generate new skeletal reaction models for methane combustion from the foundational fuel chemistry model (FFCM-1). Here, the sensitivities of the thermo-chemical variables with respect to the reaction rates are computed via the forced-optimally time dependent (f-OTD) methodology. In this methodology, the large sensitivity matrix containing all local sensitivities is modeled as a product of two low-rank time-dependent matrices. The evolution equations of these matrices are derived from the governing equations of the system. The modeled sensitivities are computed for the auto-ignition of methane at atmospheric and high pressures with different sets of initial temperatures, and equivalence ratios. These sensitivities are then analyzed to rank the most important (sensitive) species. A series of skeletal models with different number of species and levels of accuracy in reproducing the FFCM-1 results are suggested. The performances of the generated models are compared against FFCM-1 in predicting the ignition delay, the laminar flame speed, and the flame extinction. The results of this comparative assessment suggest the skeletal models with 24 and more species generate the FFCM-1 results with an excellent accuracy.

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

Reaction-diffusion modeling provides insights into biophysical carbon-concentrating mechanisms in land plants

Carbon-concentrating mechanisms (CCMs) have evolved numerous times in photosynthetic organisms. They elevate the concentration of CO 2 around the carbon-fixing enzyme rubisco, thereby increasing CO 2 assimilatory flux and reducing photorespiration. Biophysical CCMs, like the pyrenoid-based CCM (PCCM) of Chlamydomonas reinhardtii or carboxysome systems of cyanobacteria, are common in aquatic photosynthetic microbes, but in land plants appear only among the hornworts. To predict the likely efficiency of biophysical CCMs in C 3 plants, we used spatially resolved reaction-diffusion models to predict rubisco saturation and light use efficiency. We found that the energy efficiency of adding individual CCM components to a C 3 land plant is highly dependent on the permeability of lipid membranes to CO 2 , with values in the range reported in the literature that are higher than those used in previous modeling studies resulting in low light use efficiency. Adding a complete PCCM into the leaf cells of a C 3 land plant was predicted to boost net CO 2 fixation, but at higher energetic costs than those incurred by photorespiratory losses without a CCM. Two notable exceptions were when substomatal CO 2 levels are as low as those found in land plants that already use biochemical CCMs and when gas exchange is limited, such as with hornworts, making the use of a biophysical CCM necessary to achieve net positive CO 2 fixation under atmospheric CO 2 levels. This provides an explanation for the uniqueness of hornworts' CCM among land plants and the evolution of pyrenoids multiple times.

Kaste, Joshua A. M.↗
DOE OSTI2024

A Spatial-Temporal Analysis of Travel Time Gap and Inequality between Public Transportation and Personal Vehicles

The increased use of personal vehicles presents environmental challenges, prompting the exploration of public transportation as an affordable, eco-friendly alternative. However, obstacles like fixed schedules, limited routes, and extended travel times impede widespread adoption. This study investigates the temporal evolution of spatial inequality in the travel time gap between public transportation and personal vehicles, reflecting disparities across states and time periods. Analyzing Census Transportation Planning Program data for six northeastern states in 2010 and 2016 reveals no significant increase in the travel time gap, but notable growth in inequality in a few urban and disadvantaged communities. Comprehending these trends is vital for fostering equitable advancements in transportation infrastructure and enhancing public transportation competitiveness.

Pan, Melrose↗
DOE OSTI2024

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

2D analysis of tokamak divertor-plasma detachment-bifurcation with operational parameters and geometries

UEDGE simulations with density scans for various input power, transport coefficients and outer poloidal leg length are performed to study the conditions for the existence of a bifurcation-like drop of T e at the outer strike point, commonly referred to as a detachment cliff, when transitioning to a detached plasma from an attached plasma in the outer divertor as the upstream density increases (McLean et al., 2015). The simulation results show that a detachment cliff tends to occur with a higher power input regardless of diffusivities and leg length. Further analysis of change of plasma profiles at a cliff indicate that, in addition to the sharp reduction of the E x B drift fluxes in the outer divertor studied in Jaervinen et al., (2018), the substantial change of the Mach number in the outer divertor and the decrease of the outer mid-plane T e due to the radiation front moving across the separatrix into the confinement region above the X-point consistently occur for all UEDGE density scans that have a detachment cliff. UEDGE time-dependent simulation of the evolution of a detachment cliff shows that the rapid increase of radiation above the X-point occurs in a time scale of ~0.3–0.5, which could possibly be the trigger for the formation of a detachment cliff, quicker than the Mach number change in a time scale of ~1 ms and the drop of T e in a time scale of ~2–3 ms in the outer divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Shear Response of Ionizable Polymer Melts at the Crossover from Ionomers to Polyelectrolytes

Nonlinear shear response of polymers is affected by inherent barriers for diffusive motion, including entanglements and topology. In ionizable polymers, ionic clusters further constrain the intrinsic dynamics of the polymers, significantly enhancing their viscosity. Here, using fully atomistic molecular dynamics simulations, the nonlinear shear response of ionizable polymers is presented, across the transition from the ionomer regime where distinctive clusters dominate the structure to the polyelectrolyte regime where clusters percolate, in polystyrene randomly sulfonated with fractions of SO 3 − groups of f = 0.20 and 0.35, in pristine and tetrahydrofuran (THF) swollen polyelectrolyte melts. For f = 0.20, the ionic clusters first fracture into smaller clusters followed by splitting into individual ionic groups and eventually reform. At higher f, the clusters morph in shape but do not break under high shear. At very high shear rates, all of the chains stretch and recoil rapidly. As the shear rate is reduced, some chains stretch and recoil, while others remain largely unaffected by the shear. Macroscopically, for all systems, the shear viscosity displays initially an elastic response, followed by nonlinear shear stress overshoot and, eventually, a steady state. The evolution of viscosity with time and shear reflects that of the ionizable domains.

cluster chemistry↗
DOE OSTI2025

Footprints of Worldwide Adaptation in Structured Populations of Drosophila melanogaster Through the Expanded DEST 2.0 Genomic Resource

Abstract Large-scale genomic resources can place genetic variation into an ecologically informed context. To advance our understanding of the population genetics of the fruit fly Drosophila melanogaster, we present an expanded release of the community-generated population genomics resource Drosophila Evolution over Space and Time (DEST 2.0; https://dest.bio/). This release includes 530 high-quality pooled libraries from flies collected across six continents over more than a decade (2009 to 2021), most at multiple time points per year; 211 of these libraries are sequenced and shared here for the first time. We used this enhanced resource to elucidate several aspects of the species' demographic history and identify novel signs of adaptation across spatial and temporal dimensions. For example, we showed that the spatial genetic structure of populations is stable over time, but that drift due to seasonal contractions of population size causes populations to diverge over time. We identified signals of adaptation that vary between continents in genomic regions associated with xenobiotic resistance, consistent with independent adaptation to common pesticides. Moreover, by analyzing samples collected during spring and fall across Europe, we provide new evidence for seasonal adaptation related to loci associated with pathogen response. Furthermore, we have also released an updated version of the DEST genome browser. This is a useful tool for studying spatiotemporal patterns of genetic variation in this classic model system.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Microbial Evolution Drives Adaptation of Substrate Degradation on Decadal to Centennial Time Scales Relevant to Global Change

ABSTRACT Understanding microbial adaptation is crucial for predicting how soil carbon dynamics and global biogeochemical cycles will respond to climate change. This study employs the DEMENT model of microbial decomposition, along with empirical mutation and dispersal rates, to explore the roles of mutation and dispersal in the adaptation of soil microbial populations to shifts in litter chemistry, changes that are anticipated with climate‐driven vegetation dynamics. Following a change in litter chemistry, mutation generally allows for a higher rate of litter decomposition than dispersal, especially when dispersal predominantly introduces genotypes already present in the population. These findings challenge the common idea that mutation rates are too low to affect ecosystem processes on ecological timescales. These results demonstrate that evolutionary processes, such as mutation, can help maintain ecosystem functioning as the climate changes.

Abs, Elsa↗
DOE OSTI2024

Electron Influence on the Parallel Proton Firehose Instability in 10-moment, Multifluid Simulations

Instabilities driven by pressure anisotropy play a critical role in modulating the energy transfer in space and astrophysical plasmas. For the first time, we simulate the evolution and saturation of the parallel proton firehose instability using a multifluid model without adding artificial viscosity. These simulations are performed using a 10-moment, multifluid model with local and gradient relaxation heat-flux closures in high-β proton–electron plasmas. When these higher-order moments are included and pressure anisotropy is permitted to develop in all species, we find that the electrons have a significant impact on the saturation of the parallel proton firehose instability, modulating the proton pressure anisotropy as the instability saturates. Even for lower β's more relevant to heliospheric plasmas, we observe a pronounced electron energization in simulations using the gradient relaxation closure. Our results indicate that resolving the electron pressure anisotropy is important to correctly describe the behavior of multispecies plasma systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Constraints on Dynamical Dark Energy from Multiple Probes in the Full Dark Energy Survey

We present results on dark energy evolution, assuming a time-dependent equation of state $w(a)=w_0+w_a(1-a)$, from growth and geometric probes using the full six-year Dark Energy Survey dataset: type Ia supernovae, baryon acoustic oscillations, and weak gravitational lensing and galaxy clustering (3$\times$2pt). The combination yields $w_0=-0.84^{+0.10}_{-0.10}$ and $w_a=-0.44^{+0.60}_{-0.55}$, the tightest constraints ever obtained from a single survey, with $2.2σ$ deviation from a cosmological constant. Adding the DESI DR2 BAO data yields $w_0=-0.84^{+0.06}_{-0.07}$ and $w_a=-0.53^{+0.33}_{-0.28}$, representing the most stringent low-redshift-only test of dynamical dark energy to date, with a $2.3σ$ deviation. In this combination, adding 3$\times$2pt doubles the constraining power. Finally, when combined with primary CMB information, we obtain $w_0=-0.82^{+0.05}_{-0.05}$, $w_a=-0.63^{+0.21}_{-0.18}$, with a $3.0σ$ deviation. We find that including 3$\times$2pt in the previously studied SN + DESI BAO + CMB combination leaves the significance essentially unchanged ($3.2 σ$ to $3.0σ$) while improving the figure of merit by $\sim$10%. We systematically investigate the impact of leaving out each one of the probes and find that the significance of the deviation from a cosmological constant ranges from 2.3 to 3.2$σ$, with best-fit parameters consistently in the region $w_0 >-1$ and $w_a <0$. Excluding SN from the all data combination yields a $2.6σ$ departure from $Λ$CDM, providing a cross-check independent of supernova photometric calibration. These results support the weak preference for evolving dark energy reported by several recent cosmological analyses. By combining growth and geometric probes from a single survey, this work realizes the multi-probe dark energy program envisioned at the inception of DES.

Abbott, T. M.C. [Cerro-Tololo InterAmerican Obs.] ↗
DOE OSTI2025

Characteristics of the IBEX Ribbon and Their Implications for a Source Region Outside the Heliopause

This paper presents a comprehensive exploration of the Interstellar Boundary Explorer energetic neutral atom (ENA) ribbon, focusing on its spatial and temporal variations over 14 yr. Methodological advancements, including a refined map modeling procedure and a new ribbon separation technique with appropriate error propagation, enable a detailed investigation of the ribbon’s features. Utilizing statistically robust metrics, this study reveals details of the ribbon across energy and time. Key findings include energy- and time-dependent variations in flux, angular radius, ribbon profile width, and higher moments. By applying these metrics, we reveal new complexity to the evolution of the ribbon over time, highlighting the nuanced relationship between it and the solar wind. Furthermore, the study examines for the first time the ribbon as it passes through the starboard/heliotail region (Lon EC 120°–180°), revealing properties distinct from other portions of the ribbon. The analysis uncovers an anticorrelation between ribbon width and flux, which provides quantitative support for a multisource ribbon created by a combination of solar wind neutrals that generate a spatiall narrow ribbon component and heliosheath neutrals giving rise to a broad component. Finally, differences in the temporal evolution of the ENA flux at different energies provide additional support that the location of the ribbon source region is beyond the heliopause.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Immersion Freezing in Particle-Based Aerosol-Cloud Microphysics: A Probabilistic Perspective on Singular and Time-Dependent Models

Cloud droplets containing immersed ice-nucleating particles (INPs) may freeze at temperatures above the homogeneous freezing threshold temperature in a process referred to as immersion freezing. In modeling studies, immersion freezing is often described using either so-called “singular” or “time-dependent” parameterizations. Here, we compare both approaches and discuss them in the context of probabilistic particle-based (super-droplet) cloud microphysics modeling. First, using a box model, we contrast how both parameterizations respond to idealized ambient cooling rate profiles and quantify the impact of the polydispersity of the immersed surface spectrum on the frozen fraction evolution. Presented simulations highlight that the singular approach, constituting a time-integrated form of a more general time-dependent approach, is only accurate under a limited range of ambient cooling rates. The time-dependent approach is free from this limitation. Second, using a prescribed-flow two-dimensional cloud model, we illustrate the macroscopic differences in the evolution in time of ice particle concentrations in simulations with flow regimes relevant to ambient cloud conditions. The flow-coupled aerosol-budget-resolving simulations highlight the benefits and challenges of modeling cloud condensation nuclei activation and immersion freezing on insoluble ice nuclei with super-particle methods. The challenges stem, on the one hand, from heterogeneous ice nucleation being contingent on the presence of relatively sparse immersed INPs, and on the other hand, from the need to represent a vast population of particles with relatively few so-called super particles (each representing a multiplicity of real particles). We discuss the critical role of the sampling strategy for particle attributes, including the INP size, the freezing temperature (for singular scheme) and the multiplicity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Ultrafast Production of NiCO and Ni Following 197 nm Photodissociation of Nickel Tetracarbonyl

Herein, we report on the ultrafast photodissociation of nickel tetracarbonyl—a prototypical metal–ligand model system—at 197 nm. Using mid-infrared transient absorption spectroscopy to probe the bound C≡O stretching modes, we find evidence for the picosecond time scale production of highly vibronically excited nickel dicarbonyl and nickel monocarbonyl, in marked contrast with a prior investigation at 193 nm. Further spectral evolution with a 50 ps time constant suggests an additional dissociation step; the absence of any corresponding growth in signal strongly indicates the production of bare Ni, a heretofore unreported product from single-photon excitation of nickel tetracarbonyl. Thus, by probing the deep UV-induced photodynamics of a prototypical metal carbonyl, this Letter adds time-resolved spectroscopic signatures of these dynamics to the sparse literature at high excitation energies.

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

Feedforward equilibrium trajectory optimization with GSPulse

One of the common tasks required for designing new plasma scenarios or evaluating capabilities of a tokamak is to design the desired equilibria using a Grad-Shafranov (GS) equilibrium solver. However, most standard equilibrium solvers are time-independent and do not include dynamic effects such as plasma current flux consumption, induced vessel currents, or voltage constraints. Another class of tools, plasma equilibrium evolution simulators, do include time-dependent effects. These are generally structured to solve the forward problem of evolving the plasma equilibrium given feedback-controlled voltages. In this work, we introduce GSPulse, a novel algorithm for equilibrium trajectory optimization, that is more akin to a pulse planner than a pulse simulator. GSPulse includes time-dependent effects and solves the inverse problem: given a user-specified set of target equilibrium shapes, as well as limits on the coil currents and voltages, the optimizer returns trajectories of the voltages, currents, and achievable equilibria. This task is useful for scoping performance of a tokamak and exploring the space of achievable pulses. The computed equilibria satisfy both Grad-Shafranov force balance and axisymmetric circuit dynamics. The optimization is performed by restructuring the free-boundary equilibrium evolution equations into a form where it is computationally efficient to optimize the entire dynamic sequence. GSPulse can solve for hundreds of equilibria simultaneously within a few minutes. GSPulse has been validated against NSTX-U and MAST-U experiments and against SPARC feedback control simulations, and is being used to perform scenario design for SPARC. The computed trajectories can be used as feedforward inputs that are connected to the feedback controller to inform and improve feedback performance. The code for GSPulse is available open-source at github.com/jwai-cfs/GSPulse_public.

equilibrium↗