SEARCH · Engineering Papers

Results for “transition path”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

503 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Using Membranes with Internal Microchannels to Prevent Drying-out during CO 2 Electrolysis

Scaling up CO 2 electrolysis is a vital aspect in the transition to manufacturing sustainable fuels and chemical compounds, satisfying the demand for chemicals and demand for storing renewable electricity. Depending on the employed catalyst, different products can be produced, such as carbon monoxide and ethylene, by applying a voltage on a CO 2 and H 2 O fed electrolyzer. CO is a desirable product according to techno-economic analysis, because it can be produced selectively using a silver catalyst and is a precursor for hydrocarbons in the Fischer–Tropsch process. In a state-of-the-art CO 2 electrolyzer, two electrodes are directly pressed against an ion-exchange membrane – this is called a zero-gap configuration. Therefore, the membrane is a crucial component for the system, since it has the role of providing a conductive medium between the electrodes. One of the challenges in CO 2 electrolysis is that water is consumed in the reaction. At high current density, this may cause the membrane's surface near the cathode to dry out, lowering efficiency or perhaps stopping the process entirely, since there is no longer a conductive medium. As a result, water management is critical for this process. In this work, we approach the drying-out challenge by studying the novel concept of a membrane with internal microchannels. These channels allow the circulation of water or an electrolyte inside the membrane, which reduces the water diffusion path and affects the membrane’s conductivity. The effects of channel geometry, location, and concentrations of electrolyte inside on water content, conductivity and overall performance are studied in a 2D COMSOL model. In addition, the effect of internal concentration of electrolyte on the membrane’s resistance, on the process performance and the K + cross-over to cathode side were investigated experimentally. Our modeling results prove that the presence of the channels can keep the membrane hydrated. The highest current densities are observed when the channel is closest to the cathode, and with smaller pores. Smaller pores are advantageous due to the trade-off between enhanced membrane conductivity and the lower conductivity of the liquid itself. If water is circulated in a large channel, it increases the membrane’s ionomer conductivity due to hydration but the overall conductivity is decreased since water is not highly conductive. Nonetheless, the results also show that a higher concentration of electrolyte inside the microchannels can significantly increase the total conductivity of the membrane, and therefore the energy efficiency of the process. These effects are most significant at higher current densities. In the experimental results, we’ve observed similar effects in terms of membrane conductivity and current density of the process – the higher the electrolyte concentration the higher the current density. Furthermore, a low concentration decreases the amount of potassium which crosses over to the cathode side, inhibiting salt deposition. We’ve concluded that a small channel, up to 90 µm wide, close to the pore with an electrolyte with a concentration of up to 10 mM could be very beneficial for the water management and energy efficiency of the process. This helps to keep the membrane hydrated at higher current densities, improves the conductivity of the membrane, and it doesn’t have significant impacts on the salt deposition.

Petrov, Kostadin Veselinov↗

DOE OSTI2024

Observation of high-pressure polymorphs in bulk silicon formed at relativistic laser intensities

Silicon polymorphs with exotic electronic and optical properties have recently attracted significant attention due to their wide range of useful band gap characteristics. They are typically formed by static high-pressure techniques, which limits the crystal structures that can be made. This constitutes a major obstacle to study these polymorphs and their incorporation into existing technology. Approaches have attempted to address this shortcoming through using dynamic conditions and chemical precursor materials. Here, we report on an approach to create unusual crystal structures deep in the bulk of a silicon crystal by irradiating it with a laser pulse at ultrarelativistic intensity of up to 7.5 × 10 19 W/ c m 2 . Laser-generated electrons with MeV energy swiftly penetrate the target with speed close to the speed of light and deposit their energy into a large volume across the whole thickness of the sample. The relativistic electron current creates, via branching propagation and ionization, high-energy-density conditions for thermodynamically nonequilibrium phase transformation paths into new crystal polymorphs. X-ray microdiffraction and synchrotron x-ray diffraction analyses indicate, along with conventional dc-Si, the presence of exotic silicon structures in the bulk of the laser intact target volume. These structures are identified as body-centered bc8-Si, rhombohedral r8-Si, hexagonal-diamond hd-Si, and the tetragonal Si-VIII, all phases of Si that have previously been made through static techniques. Additionally, simple-tetragonal st12-Si and body-centered tetragonal bt8-Si were observed along with signatures of not yet identified diffraction spots. Both st12-Si and bt8-Si have only been observed in ultrafast laser microexplosion conditions at much lower laser intensity ∼ 10 14 W/ c m 2 and within a micron-thin surface layer. The findings here are supported by direct observation of nanoparticles with high-resolution transmission electron microscopy and corresponding fast Fourier transform analysis of their interatomic distances. The presented analyses of absorbed laser energy, generation of the MeV electron current, and deposition of energy across the whole target thickness provide a solid basis for drawing the conclusion that the observed silicon polymorphs were produced because of laser-generated high-energy electrons fast-penetrating deeply into the bulk of silicon. In contrast to solid-solid transformations, the plasma-solid transitions offer a paradigm for the creation of exotic, high-energy density materials inside the bulk of the sample by using laser pulses at relativistic intensities. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

The concept of spin ice graphs and a field theory for their charges

Originally detected in rare earth pyrochlores, spin ice physics is now being artificially extended to a variety of geometries that control collective behavior and exotic properties, making graph theory their proper framework. We relate spin ice notions, such as ice rule, ice manifold, Coulomb phases, charges, and monopoles, to graph-theoretical notions, such as balance, in/out-degrees, and Euler paths. We then propose a field-theoretical treatment in which topological charges and monopoles are the degrees of freedom, while the binary spins are subsumed in an entropic interaction among charges. We show that for a spin ice on a graph in a Gaussian approximation, the kernel of the entropic interaction is the inverse of the graph Laplacian, and we compute screening functions from the graph spectra as Green operators for the screened Poisson problem on a graph. We then apply the treatment to star graphs, tournaments, cycles, and regular spin ice in different dimensions. Our aim is twofold: to set spin ice physics in a proper graph setting, where only topological rather than geometrical notions hold, and to invite graph theorists to contribute their powerful tools to the field of spin ice.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Role of chemical disorder on radiation-induced defect production and damage evolution in NiFeCoCr

Understanding chemical disorder in many concentrated solid solution alloys (CSAs) at the levels of electrons and atoms has attracted increasing attention as a path forward to reveal and identify underlying mechanisms for extraordinary mechanical properties and improved radiation tolerance. Single-phase NiFeCoCr CSA is a common base for many high-entropy alloys (HEAs) that have shown improved mechanical strength and radiation tolerance. In this study, defect production and damage evolution in NiFeCoCr under ion irradiation at room temperature to dose over 20 dpa are determined using ion channeling technique along both <100> and <110> directions utilizing multiple probing beam energies. The results obtained from the multi-axial and multi-energy channeling analysis are compared with those previously obtained for Ni crystals irradiated under similar conditions. The influence of chemical complexity on defect production and clustering at early-stage under room temperature irradiation up to dose of 1 dpa is discussed based on positron annihilation spectroscopy results. Defect structure evaluation in Ni and NiFeCoCr is also discussed based on transmission electron microscopy results over a prolonged irradiation at both room and elevated temperatures. Compared with chemically complex NiFeCoCr, larger dislocation loops thus less lattice strain are expected to form in pure Ni. Moreover, the role of chemical disorder in this CSA is also investigated based on ab initio calculations using large supercells. Finally, to understand the impact of chemical complexity effect on defect structure evolution, this integrated research effort attempts to link the relatively large charge redistribution due to difference in valence electron counts resulting from alloying different 3d transition metal elements, moderate lattice distortion arising from similar adaptable atomic size, and notable suppressed or delayed damage evolution in NiFeCoCr.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Automated Generation of Integrated Digital and Spiking Neuromorphic Machine Learning Accelerators

The growing numbers of application areas for artificial intelligence (AI) methods have led to an explosion of domain-specific accelerators that could support every new machine learning (ML) algorithm advancement, clearly highlighting the need for a capability to quickly and automatically transition from algorithm definition to hardware implementation and explore design space along a variety of SWaP (size, weight and Power). The software defined architectures (SODA) synthesizer implements a compiler-based modular infrastructure for the end-to-end generation of machine learning accelerators from high-level frameworks to hardware description language. At the same time, neuromorphic computing, by mimicking how the brain operates, promises to perform artificial intelligence tasks at efficiencies orders of magnitude higher than the current conventional tensor-processing based accelerators, as demonstrated by a variety of specialized designs leveraging Spiking Neural Networks (SNNs). Nevertheless, the mapping of an artificial neural network (ANN) to solutions supporting SNNs is still a non-trivial and very device-specific task, and completely lack the possibility to design hybrid systems that integrate conventional and spiking neural models. In this paper we discuss the support for such an integrated generation leveraging the SODA Synthesizer framework and its modular structure. In particular, we present a new MLIR dialect (part of the SODA frontend) that allows expressing spiking neural network features (e.g., available resources, spiking sequences, analog signal reading, etc.) and illustrate how it enables mapping to Spiking Neurons and deployment to the related specialized hardware (which, in the digital domain, could be generated through the other existing layers of the SODA Synthesizer). We then discuss the opportunities for even deeper integration afforded by the hardware compilation infrastructure, providing a path towards the generation of complex heterogeneous artificial intelligence systems.

Curzel, Serena↗

DOE OSTI2023

Spectroscopic features of dissolved iodine in pristine and gamma-irradiated nitric acid solutions

While iodine speciation is important for a wide range of nuclear activities, understanding the mechanisms of the transformations of iodine between chemical forms and the sensitivity of these transitions to solution conditions and exposure to radiation remains an active area of research. This work curates spectroscopic data from several experimental techniques and establishes their sensitivity and limitations in detecting changes in iodine speciation in both neutral and acidic regimes. The techniques include Raman spectroscopy, Fourier Transform Infrared (FTIR) spectroscopy, 127 I NMR spectroscopy, and ultraviolet-visible (UV-Vis) spectroscopy. Analysis of these data indicates that these commonly accessible spectroscopies often have dynamic ranges of measurable concentrations that do not always overlap between all techniques. The experimental techniques are disparately sensitive to iodide (I - ), molecular iodine (I 2 ), iodate (IO 3 - ), and periodate (IO 4 - ) species. Raman, FTIR, and NMR spectra were subsequentially analyzed using two-dimensional correlation analyses to generate high-resolution autocorrelation spectra. Here, the use of these spectroscopies is then extended to tracking acidification-induced and gamma irradiation-induced transformations of dissolved sodium iodate in deionized water and concentrated nitric acid. Both dissolution into nitric acid and irradiation with a gamma source are demonstrated to perturb the iodine speciation promoting their assembly into molecular iodine (I 2 ) and/or triiodide (I 3 - ). While I 2 and I 3 - species are undetectable with FTIR spectroscopy and 127 I NMR spectroscopy, the species can be detected with UV-Vis spectroscopy, and in some instances, I 3 - can be detected with Raman spectroscopy in the low wavenumber region. Ultimately, the results of this work provide a path to designing optimal combinations of techniques to detect forms of iodine across a wide range of concentrations and conditions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2026

Blueprint for DOE Quantum Supercomputing: Ensuring U.S. Leadership in the Quantum Decade

Quantum computing stands at the threshold of a transformative decade, where the field will evolve from small-scale demonstrations toward practical scientific computing at scale. This Blueprint identifies fault-tolerant quantum computers (FTQCs) as a viable, scalable, and broadly applicable path to achieving “quantum scientific utility,” defined as solving scientifically valuable problems beyond the reach of conventional, classical computers. This capability is expected to show scientific demonstrations in the late 2020s and to mature in the early-to-mid 2030s. This Blueprint outlines a strategy to prepare the U.S. Department of Energy (DOE) for FTQCs and their integration into the U.S. national scientific computing infrastructure. Its purpose is to identify the steps, milestones, and research directions necessary for DOE to enable initial deployment of FTQCs in 2028 as a scientific tool for the nation and mature this capability into the 2030s. DOE has a long history of supporting quantum information science and technology, contributing significantly to research advancements, training a quantum-ready workforce, and providing access to early small-scale quantum hardware. Given recent demonstrations of logical operations on error-corrected logical qubits and the advancement of commercial hardware roadmaps, DOE should begin preparations for large-scale, fault-tolerant quantum computing deployment for DOE science missions. This Blueprint proposes that DOE focus on (1) deploying first-generation scientifically relevant quantum computers with at least 100 logical qubits and performing at least 10,000 to 100,000 hard logical operations in scientifically relevant calculations; (2) developing essential FTQC programming competencies, system software, and facility readiness; and (3) investing in cutting edge focused R&D that fosters breakthroughs in scientific applications, algorithms, and logical architectures needed to accelerate the advent of scientific utility. This effort will position DOE to transition to larger systems: production-scale quantum computers that comprise 1,000 to 10,000 logical qubits, perform 1 to 10 billion hard logical operations, and execute scientifically useful computations at scale. Achieving these goals will require DOE facilities to evolve with urgency to support scientific campaigns that integrate quantum and classical computing resources into efficient workflows, novel software and firmware environments for compiling and routing quantum programs on FTQC machines, and suitable infrastructure for quantum hardware. It will also require further development and optimization of scientific applications from the fields of materials science, quantum chemistry, and high-energy and nuclear physics. The Blueprint calls for transformative R&D and collective action to accelerate the advent of scientific quantum utility and bring it within reach by 2028.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Reducing module soiling with scalable and robust photocatalytic coatings

The air-glass interface at the front of a photovoltaic (PV) module reflects approximately 4% of incident light, decreasing the potential power output of the module by the same amount. Today’s modules reduce this loss by adding a low-refractive-index (1.25-1.30) SiO2 coating to the sunward side of the module glass; this antireflection coating recovers approximately 3% of the 4% light that would otherwise be lost. While such antireflection coatings work very well on clean, new modules, they do not inhibit soiling—the accumulation of soilants such as dust, pollen, soot, or other foreign material—on the module glass, and soilants reflect and scatter incident light. An improved coating would serve provide not only an antireflection effect, but also an anti-soiling effect. The goal of this project was to develop such a coating and provide a path for it to be manufactured in the U.S. The project successfully designed and fabricated coatings that provided >3% transmittance gain compared to bare glass (matching the performance of commercial antireflection coatings) and displayed anti-soiling behavior in standard laboratory soiling effects. This was achieved by using a Swift Coat proprietary coating deposition technique, aerosol impact-driven assembly (AIDA), to control the porosity and thus refractive index of coatings of photocatalytic materials—such as TiO2—that would otherwise increase (instead of decrease) reflection. These combined antireflection/anti-soiling coatings passed PV industry standard module reliability tests as well as coating-specific abrasion tests, showing that they have the durability needed for decades in the field. Swift Coat scaled the AIDA hardware and deposition process to make mini-modules that were monitored for nearly two years during field tests administered by a third party, as well as demonstrated scaling to the widths of full-sized modules. The fielded mini-modules outperformed reference modules (with commercial antireflection coatings) in two locations, providing a 1% absolute average performance boost and larger increases during periods of heavier soiling. Swift Coat’s cost analysis indicated a coating manufacturing cost below the sales price of today’s antireflection coatings. More than five module manufacturers sampled and assessed the coatings, and three provided letters of support. The coating developed in this project increases the energy output of PV modules, thereby decreasing the cost per kilowatt-hour of solar energy generated. Cheaper solar electricity benefits the public by accelerating the transition to a stable, affordable, carbon-free energy economy. In addition, for select applications in which PV modules are highly visible—such as on residential rooftops—the coating provides an aesthetic benefit because it stays cleaner than today’s modules. Finally, Swift Coat and its prospective customers are U.S. companies, and successful commercialization of this technology will provide U.S. jobs and a secure solar supply chain.

14 SOLAR ENERGY↗

DOE OSTI2022

Li 8 MnO 6 : A Novel Cathode Material with Only Anionic Redox

In Li-excess transition-metal-oxide cathode materials, anionic oxygen redox can offer high capacity and high voltages, although peroxo and superoxo species may cause oxygen loss, poor cycling performance, and capacity fading. Previous work showed that undesirable formation of peroxide and superoxide bonds is controlled to some extent by Mn substitution, and the present work uses density functional calculations to examine the reasons for this by studying the anionic redox mechanism Li 8 MnO 6 . This material is obtained by substituting Mn for Sn in Li 8 SnO 6 or for Zr in Li 8 ZrO 6 , and we also compare to previous work on those materials. The calculations predict that Li 8 MnO 6 is stable at room temperature (with a band gap of 3.19 eV as calculated HSE06 and 1.82 eV as calculated with the less reliable with PBE+U), and they elucidate the chemical and structural effects involved in the inhibition of oxygen release in this cathode. Throughout the whole delithiation process, only O 2- ions are oxidized. The directional Mn-O bonds formed from unfilled 3d orbitals effectively inhibit the formation of O-O bonds, and the layered structure is maintained even after removing 3 Li per Li 8 MnO 6 formula unit. The calculated average voltage for removal of 3 Li is 3.69 V by HSE06, and the corresponding capacity is 389 mAh/g. The high voltage of oxygen anionic redox and the high capacity result in a high energy density of 1436 Wh/kg. The Li-ion diffusion barrier for the dominant interlayer diffusion path along the c-axis is 0.57 eV by PBE+U. Finally, these results help us to understand the oxygen redox mechanism in a new lithium-rich Li 8 MnO 6 cathode material and contribute to the design of high-energy-density lithium-ion-battery cathode materials with favorable electrochemical properties based on anionic oxygen redox.

Li8MnO6↗

DOE OSTI2021

Burning velocities of R-32/O 2 /N 2 mixtures: Experimental measurements and development of a validated detailed chemical kinetic model

This work entails characterizing the flammability of the refrigerant R-32 (CH 2 F 2 ) by both experimental measurements and modeling. Burning velocities S u were measured using a constant-volume spherical-flame method for R-32/O2/N 2 mixtures with O 2 /N 2 ratios ranging from 21% (synthetic air) to 40%, pressures of (1 to 3) bar, and equivalence ratios $\phi$ of (0.8 to 1.3). Based on a critical assessment of available data, and extended by our own calculations, a detailed chemical kinetic model was developed and key reactions determined using reaction path and sensitivity analyses. Initiation and combustion were identified as distinct kinetic regimes and burning velocities were found to be controlled by two primary reactions: unimolecular decomposition of CH 2 F 2 → CHF + HF and the subsequent reaction, CHF + O 2 → CHFO + O, the latter reaction initiating the radical chain propagating and branching by producing O atoms. Sensitive rate constants in the kinetic model were critically adjusted within their uncertainties and current knowledge bounds to best fit the experimental burning velocities. We found that rate constants in the model could be adjusted to match a given experimental S u for specific conditions (O 2 loading, P, T, $\phi$). This, however, then fixes predicted burning velocities for other all conditions within (3 to 4)% if physically realistic rate parameters are maintained. Thus, the entire set of experimental data is fit, not just to particular conditions. Relative random uncertainties in the experimental Su measurements were (4 to 6)%, but assumptions made for thermal radiation lost by the burned gas in the spherical-flame experiments add an additional systematic uncertainty. Systematic differences between the limiting cases of adiabatic (no thermal radiation lost) and optically-thin (all thermal radiation lost) varied significantly with conditions and ranged from (4 to 30)% at high to low velocities, respectively, translating into uncertainties of (2 to 15)% considering the average of two limiting cases. Comparison of experimental and kinetically modeled Su values suggests that the burned gas tends towards the optically-thin limit at the lowest pressures and fuel loadings and toward the adiabatic limit at the highest pressures and loadings. We tested and found support for this conclusion with a detailed analysis as a function of all the conditions (T, P, % O 2 , $\phi$). This behavior appears to transition from optically-thin to adiabatic as the density of the initial fuel increases, which results in increased CO 2 in the burned gas and thus increased absorption of the thermal radiation (consistent with the Beer-Lambert Law). The validated detailed model based on evaluated kinetics is shown to accurately predict burning velocities for R-32 O 2 /N 2 mixtures over a wide range of conditions and provides a reliable basis for extrapolation to other conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Revisiting the pseudo-supercritical path method: An improved formulation for the alchemical calculation of solid–liquid coexistence

Alchemical free energy calculations via molecular dynamics have been applied to obtain thermodynamic properties related to solid–liquid equilibrium conditions, such as melting points. In recent years, the pseudo-supercritical path (PSCP) method has proved to be an important approach to melting point prediction due to its flexibility and applicability. In the present work, we propose improvements to the PSCP alchemical cycle to make it more compact and efficient through a concerted evaluation of different potential energies. The multistate Bennett acceptance ratio (MBAR) estimator was applied at all stages of the new cycle to provide greater accuracy and uniformity, which is essential concerning uncertainty calculations. In particular, for the multistate expansion stage from solid to liquid, we employed the MBAR estimator with a reduced energy function that allows affine transformations of coordinates. Free energy and mean derivative profiles were calculated at different cycle stages for argon, triazole, propenal, and the ionic liquid 1-ethyl-3-methyl-imidazolium hexafluorophosphate. Comparisons showed a better performance of the proposed method than the original PSCP cycle for systems with higher complexity, especially the ionic liquid. Furthermore, a detailed study of the expansion stage revealed that remapping the centers of mass of the molecules or ions is preferable to remapping the coordinates of each atom, yielding better overlap between adjacent states and improving the accuracy of the methodology.

25 ENERGY STORAGE↗

DOE OSTI2025

Molecular axis distribution moments in ultrafast transient absorption spectroscopy: A path toward ultrafast quantum state tomography

In ultrafast time-resolved experiments with gas phase molecules, the alignment of the molecular axis relative to the polarization of the interacting laser pulses plays a crucial role in determining the dynamics following this light–matter interaction. The molecular axis distribution is influenced by the interacting pulses and is intrinsically linked to the electronic coherences of the excited molecules. However, in typical theoretical calculations of such interactions, the signal is either calculated for a single molecule in the molecular frame or averaged over all possible molecular orientations to compare with the experiment. Such averaging removes information about anisotropy in the molecular-axis distribution, even though anisotropic contributions can play a significant role in the measured experimental signal. Here, we calculate the laboratory frame transient electronic first-order polarization [P (1) ] spectra in terms of separated molecular frame and laboratory frame quantities. The laboratory frame polarizations are compared with orientation-averaged quantum master equation calculations, demonstrating that orientation-averaging captures only the isotropic contributions. We show that our formalism also allows us to evaluate the anisotropic contributions to the spectrum. Lastly, we discuss the application of this approach to achieve ultrafast quantum state tomography using transient absorption spectroscopy and field observables in nonlinear spectroscopy.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Phenomenological R -Matrix parameterization of direct, doorway, and compound nuclear reactions [Abstract]

Although formal expressions for scattering matrix accounting for direct, doorway, and compound nuclear (CN) resonant reactions have been derived several decades ago in both the transition ( T -)matrix formalism and the reactance ( K -)matrix formalism, the absence of corresponding expressions in phenomenological R -matrix formalism has limited the application of the latter to CN resonant reactions only. We remove this limitation by parameterizing direct, doorway, and CN resonant reactions in a phenomenological R -matrix scattering matrix, and provide a parameterization for a corresponding Reich-Moore approximation of eliminated capture channels. Direct reactions induce (previously neglected) mixing among the incoming or outgoing R -matrix channel wave functions, parameterized by real and orthonormal channel-rotation matrix, M , whereby the original scattering matrix U is transformed into M T UM . Any real and orthonormal matrix, M , can be equivalently expressed as e η , where η is a real and skew-symmetric 2 rotation-generating matrix that subsequently yields a more intuitive parameterization of eliminated direct capture reactions in Reich-Moore approximation. A phenomenological R -matrix parameterization of doorway reactions is inferred by equating the expression for reactance ( K -)matrix, given in terms of Brune’s alternative R -matrix parameterization, to a corresponding expression derived using Feshbach’s projection operator formalism. Assuming that all doorway states, just like CN states, are confined within spheres defined by R -matrix channel radii, a new R -matrix-like term induced by doorway states is gleaned, wherein each doorway state is parameterized by its energy, width, and the strength of its coupling to each CN state. Since a Reich-Moore approximation for retained-channel scattering matrix ought to approximate the effect of eliminated capture channels taking place via direct, doorway, or CN reactions, each of the three kinds of reactions contributing to the capture entails a corresponding Reich-Moore parameterization in a first-order approximation: direct contribution is parameterized by introducing finite diagonal elements of a retained-channel rotation-generating matrix, doorway contribution is parameterized by doorway capture widths, while CN contribution is parameterized by conventional Reich-Moore capture widths. We will present evidence of direct and doorway reactions observed in recent measurements of resolved resonance cross sections at the Gaerttner LINAC Center at Rensselaer Polytechnic Institute, and will outline a path for implementing this new R -matrix parameterization into the SAMMY nuclear data evaluation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗

DOE OSTI2025

High-Resolution Full-Field Structural Microscopy of the Voltage-Induced Filament Formation in VO 2 -Based Neuromorphic Devices

In order to make neuromorphic functions in memristive devices more efficient, information about the structural properties of filament formation at the micro- and mesoscopic scales is necessary. Despite extensive research on VO 2 , a key material due to its filament formation, local operando structural measurements remain challenging and often involve destructive specimen preparation and long rastering times, greatly limiting the scope of experimental studies. Utilizing dark-field X-ray microscopy (DFXM), a fullfield imaging modality, structural signatures of the filament formation process operando are revealed in VO 2 devices. DFXM experiments illustrate that rutile filaments contain isolated monoclinic clusters, indicating structural nonuniformity interior to the filament. The formation of the rutile phase beneath device electrodes was shown to precede filament development, followed by the formation of filament paths guided by nucleation sites within the device. Finally, a medium-term (<30 min) memory mechanism is observed in VO 2 , mediated by sites within the device gap that tend to switch at significantly lower voltages after electrical cycling, a tendency that persists through a brief thermal reset. High spatial resolution, large field-of-view, structure selectivity, and fast signal acquisition of DFXM provided insight into structural features of the filamentary channel and surrounding regions during voltage cycling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Drive-pressure optimization in ramp-wave compression experiments through differential evolution

Ramp-wave dynamic-compression experiments are used to examine quasi-isentropic loading paths in materials. The gradual and continuous increase in pressure created by ramp waves make these types of experiments ideal for studying nonequilibrium material behavior, such as solidification kinetics. In ramp-wave compression experiments, the input drive pressure to the experimental setup may be exerted through one of a number of different mechanisms (e.g., magnetic fields, gas-gun-driven impactors, or high-energy lasers) and is generally required for simulating such experiments. Yet, regardless of the specific mechanism, this drive pressure cannot be measured directly (measurements are generally taken at a location near the back of the experimental setup through a transparent window), leading to an inverse problem where one must determine the drive pressure at the front of the experimental setup (i.e., the input) that corresponds to the particle velocity (the output) measured near the back of the experimental setup. Furthermore, we solve this inverse problem using a heuristic optimization algorithm, known as differential evolution, coupled with a multiphysics, hydrodynamics code that simulates the compression of the experimental setup. By running many rounds of forward simulations of the experimental setup, our optimization process iteratively searches for a drive pressure that is optimized to closely reproduce the experimentally measured particle velocity near the back of the experimental setup. While our optimization methodology requires a significant number of hydrodynamics simulations to be conducted, many of these can be performed in parallel, which greatly reduces the time cost of our methodology. One novel aspect of our method for determining the drive pressure is that it does not require physical modeling of the drive mechanism and can thus be broadly applied to many types of ramp-compression experiments, regardless of the drive mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗