Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transport phenomena”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Combustion instability analysis for liquid propellant rocket engines

The multi-dimensional numerical model has been developed to analyze the nonlinear combustion instabilities in liquid-fueled engines. The present pressure-based approach can handle the implicit pressure-velocity coupling in a non-iterative way. The additional scalar conservation equations for the chemical species, the energy, and the turbulent transport quantities can be handled by the same predictor-corrector sequences. This method is time-accurate and it can be applicable to the all-speed, transient, multi-phase, and reacting flows. Special emphasis is given to the acoustic/vaporization interaction which may act as the crucial rate-controlling mechanism in the liquid-fueled rocket engines. The subcritical vaporization is modeled to account for the effects of variable thermophysical properties, non-unitary Lewis number in the gas-film, the Stefan flow effect, and the effect of transient liquid heating. The test cases include the one-dimenisonal fast transient non-reacting and reacting flows, and the multi-dimensional combustion instabilities encountered in the liquid-fueled rocket thrust chamber. The present numerical model successfully demonstrated the capability to simulate the fast transient spray-combusting flows in terms of the limiting-cycle amplitude phenomena, correspondence between combustion and acoustics, and the steep-fronted wave and flame propagation. The investigated parameters include the spray initial conditions, air-fuel mixture ratios, and the engine geometry. Stable and unstable operating conditions are found for the liquid-fueled combustors. Under certain conditions, the limiting cycle behavior of the combusting flowfields is obtained. The numerical results indicate that the spray vaporization processes play an important role in releasing thermal energy and driving the combustion instability.

Kim, Y. M.↗
DOE OSTI2025

Laser-driven flash x-ray radiography of a shocked metallic foil

Characterizing hydrodynamic instability evolution in millimeter-scale, high-Z foils is crucial for understanding complex phenomena in high-energy-density physics. Here, we demonstrate a proof-of-concept, laser-driven flash x-ray radiography platform tailored for two-dimensional linear density mapping in shocked high-Z foils. Using chromium (Cr) foils with internal shockwaves (∼100 μm width), our platform achieves a spatial resolution of 59.8 ± 1.4 μm by employing a broadband x-ray source extending into the hundreds of keV range. The setup combines a compound parabolic concentrator cone with a tantalum wire target, a magnetic field to deflect residual transmitted electrons, and a copper casing to shield the sides and rear of the image plate pack. By varying the delay of the short-pulse beam driving the flash x-ray source, we resolve shockwave dynamics, specifically the velocity, position, width, and density profile, within the Cr foil. Reported experimental results are consistent with the corresponding hydrodynamics and radiation transport simulations, which accurately reproduce the measured electron and x-ray source terms. These developments enable the conversion of shockwave radiographs into two-dimensional density maps, enhancing interpretability for hydrodynamic instability evolution applications and validating the simulation approach.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Correlating Tomographic Chemical Inhomogeneity and Low Energy Electronic Structure in Layered Quantum Materials

Photoemission spectroscopy (PES) is a suite of experimental tools to learn about the electronic and chemical structure of materials and surfaces. Normally implemented in a surface-sensitive manner, the research performed under this grant focused on pushing PES into less explored regimes, to reveal bulk electronic structure, to reveal tomographic (layer-resolved) chemistry and electronic structure of layered materials and heterostructures, and to reveal interface phenomena at the junction of two different materials. Standing wave (SW) spectroscopies have also been applied to PdCoO2, a material of interest due to its high conductivity and electron-hydrodynamic tendencies. This material can be modeled as an alternating layered structure consisting of metallic Pd layers and insulating CoO2 layers. Using SW XPS, the total electronic structure has been decomposed into contributions from the two layers, and computations highlighted the different many-body interactions in the two layers (Comm. Phys. 4, 143 (2021)). We have also used hard-x-ray angle-resolved photoemission spectroscopy (ARPES), to investigate LaB6, a technologically important material with widespread application as a cathode material for electron microscopes. We measured the bulk electronic structure of this material and found that the one-step model of photoemission better captured the electronic structure and correlations. This model treats the three steps of the photoemission process—excitation, transport of the photoelectron to the crystal surface, and escape into the vacuum—as a single quantum mechanically coherent process (Phys. Rev. Mater. 5, 055002 (2021)). We also applied x-ray photoelectron spectroscopy, implemented in a near total reflection grazing incidence geometry to elucidate technologically relevant interfaces, such as those between a substrate and photoresist (J. Phys. D: Appl. Phys. 54 464002 (2021)).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies (M2SF-21SN010301072)

The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heatgenerating nuclear waste disposition in deep-seated clay/shale/argillaceous rock. International collaboration activities such as heater tests, continuous field data monitoring, and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum and discrete) to tackle issues related to flow and transport at various scales of the host-rock, its interactions with barrier materials, and EBS design concept.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Chemical Sensing and Chemoresponsive Pumping with Conical-Pore Polymeric Membranes

Synthetic membranes containing asymmetrically shaped pores have been shown to rectify the ionic current flowing through the membrane. Ion-current rectification means that such membranes produce nonlinear current–voltage curves analogous to those observed with solid-state diode rectifiers. In order to observe this ion-current rectification phenomenon, the asymmetrically shaped pores must have pore-wall surface charge. Pore-wall surface charge also allows for electroosmotic flow (EOF) to occur through the membrane. We have shown that, because ion-current is rectified, EOF is likewise rectified in such membranes. This means that flow through the membrane depends on the polarity of the voltage applied across the membrane, one polarity producing a higher, and the opposite producing a lower, flow rate. As is reviewed here, these ion-current and EOF rectification phenomena are being used to develop new sensing technologies. Results obtained from an ion-current-based sensor for hydrophobic cations are reviewed. In addition, ion-current and EOF rectification can be combined to make a new type of device—a chemoresponsive nanofluidic pump. This is a pump that either turns flow on or turns flow off, when a specific chemical species is detected. Results from a prototype Pb 2+ chemoresponsive pump are also reviewed here.

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

Synthesis and Observation of Emergent Phenomena in Epitaxial Heusler Compound Heterostructures

The proposal was on the synthesis and observation of emergent phenomena in epitaxial Heusler compound heterostructures. The large range of properties and number of Heusler compounds opens up a wide number of potential compounds that will exhibit emergent phenomena. The similarity, large range of relatively inexpensive, large area, high crystal quality, III-V bulk substrates, lattice parameters and the ability to tune the lattice parameters through ternary or quaternary III-V compound semiconductor epitaxial growth, makes III-V semiconductors an ideal choice for substrates for epitaxial growth of Heusler compounds. A number of half Heusler compounds have been predicted to exhibit band inversion, making them topological and are therefore expected to exhibit spin-momentum locked topological surface states with linear dispersion. Others are predicted to be semimetals with Weyl points and others semiconducting and magnetic. During the course of this grant, emphasis has been on investigating Heusler compounds with emergent phenomena and demonstrating the ability to tune their properties through alloying and strain. We have grown toplogical semimetal (PtLuSb, PtMnBi), Weyl (Co 2 MnAl, Co 2 TiGe), half metal (PtMnSb, Co 2 MnSi, Co 2 MnAl x Si 1-x , Co2FeAl), and semiconducting (CoTiSb, NiTiSn) and tuned their properties through alloying and epitaxial strain. We also investigated the closely related materials of rare-earth monopnictide, some of which have also been predicted to be topological. During the attempts to grow the PtMnBi, it was discovered that Bi, another predicted topological material when ultrathin, could be grown epitaxially on InSb, results for which are also reported here. The main focus for this effort has been on using variable photon energy and spin-dependent angle resolved photoemission (ARPES) to determine bulk band structure and surface states of pristine epitaxial films grown on III-V semiconductor and MgO substrates and correlate results with theory and transport measurements. Theory has been critical to interpretation of experimental results and has been essential in guiding experiments. The research benefited from several strong collaborations between the PIs and the beamline scientists at the Advanced Light Source at Lawrence Berkeley Laboratory, the Stanford Linear Accelerator Center (SLAC) at Stanford and at the Max Lab at Lund University in Sweden. The strong experiment - theory collaboration between the PI’s groups, the Palmstrøm group at UCSB and the Janotti group at the University of Delaware, has been critical for interpreting the experimental ARPES and magnetotransport measurements results and making predictions to guide experiments. Weekly interactive Zoom meetings made this work well. A collaboration between the Palmstrøm group and Dr. Alexei Fedorov at the Advanced Light Source (ALS) resulted in significant modifications to his end chamber to accommodate the vacuum suitcase that was designed and constructed in the Palmstrøm group at UCSB. In collaboration with beamline scientists, Drs. Makoto Hashimoto and Donghui Lu at SLAC, Palmstrøm made modifications to the vacuum suitcase and developed special sample holders that allowed samples to be grown in the Palmstrøm MBE systems at UCSB and transported in the UHV vacuum suitcase to SLAC for ARPES measurements. The development of the vacuum suitcase was essential for this grant as it has allowed variable photon energies to be used to identify surface versus bulk states on samples that could not be capped and decapped using As- or Sb-capping layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1989

Saltation thresholds and entrainment of fine particles at Earth and Martian pressures

An open circuit wind tunnel designed to operate in a large vacuum chamber was built at NASA-Ames to investigate saltation threshold, flux, deflation rates, and other aeolian phenomena on the planet Mars. The vacuum chamber will operate at pressures as low as 4 mbar, and the tunnel operates at windspeeds as high as 150 m/sec. Either air or CO2 can be used as a working fluid. It was found that, to a first order approximation, the same dynamic pressure was required at Martian pressure to entrain or saltate particles as was required on Earth, although wind and particle speed are considerably higher at Martian pressure. A 2nd wind tunnel, designed to operate aboard the NASA KC-135 0-g aircraft to obtain information on the effect of gravity on saltation threshold and the interparticle force at 0-g, is also described and test data presented. Some of the experiments are summarized and various aspects of low pressure aeolian entrainment for particles 12 to 100 micron in diameter are discussed, some of them unique to low pressure testing and some common in Earth pressure particle transport testing. The facility, the modes of operation, and the materials used are described.

Leach, Rodman↗
DOE OSTI2023

Modeling and Simulation of Xe-100-type Pebble Bed Gas-Cooled Reactor with SCALE

The US Department of Energy (DOE) announced the Advanced Reactor Demonstration Program (ARDP) to accelerate the deployment of advanced reactor concepts. Awardees of ARDP funds are expected to demonstrate the operation of an advanced reactor within 7 years of receiving the award. X-Energy’s advanced reactor concept, the Xe-100, was selected as one of two advanced reactor concepts to receive funding to demonstrate the operation of its high-temperature gas-cooled pebble-bed reactor before the end of this decade. As a result of this push to bring advanced reactors to maturation and commercialization, transition and deployment scenario studies are being performed under the Systems Analysis and Integration (SA&I) campaign within the DOE Office of Nuclear Energy (DOE-NE) to evaluate the transition of the current US commercial fleet of light-water reactors (LWRs) to a future fleet of advanced reactors consisting of a mix of ARDP type reactor concepts and advanced LWRs. To accurately evaluate the front- and back-end resource requirements, it is important to perform reactor physics calculations to determine the discharge burnup and isotopic content, fuel residence time, as well as other parameters. For this purpose, a summer project funded by the SA&I campaign allowed for the setup of SCALE models for full-core Xe-100 type high-temperature gas-cooled pebble-bed reactor and a Xe-100 type slice using publicly available information. The core-averaged equilibrium compositions and zone-wise equilibrium compositions for the slice and 3D models, respectively, were obtained following an iterative depletion method developed by Bostelmann et al. using SCALE’s reactor physics sequence TRITON. The slice model was used with TRITON to generate burnup-dependent cross section libraries at different temperatures which can be used with SCALE’s ORIGAMI code to rapidly determine fuel inventory and therefore to perform quick sensitivity studies on parameters such as the pebble location in the core. The SCALE/TRITON transport and depletion calculation for the Xe-100 type slice model indicates that the isotopic concentrations are in good agreement at 1,300 effective full power days (EFPD) for 235 U. An analysis of 236 U results match 239 Pu results would seem to indicate a typographical error in Mulder and Boyes wherein the reported results of 236 U and 239 Pu are reversed. In addition to SCALE/TRITON calculations, a new capability within SCALE/ORIGAMI for the simulation of pebble-bed reactors was used to study the burnup sensitivity with respect to the pebble pathway through the core. The SCALE/ORIGAMI results show that pebbles that travel closer to the reflector for the entire depletion history have a higher burnup than pebbles that travel through the middle of the core because of the higher thermal to fast flux ratio near the reflector. Consequently, a pebble’s burnup is strongly affected by the pebble’s pathway for each pass. Additional phenomena such as temperature distributions in the core and different travel times of the pebbles in the individual radial zones further affect the burnup distribution. The sensitivity of the discharge vector to the pebble pathways taken during each pass can be evaluated in the future using SCALE/ORIGAMI now that the SCALE inputs have been established.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Extended quantum anomalous Hall states in graphene/hBN moiré superlattices

Electrons in topological flat bands can form new topological states driven by correlation effects. The pentalayer rhombohedral graphene/hexagonal boron nitride (hBN) moiré superlattice was shown to host fractional quantum anomalous Hall effect (FQAHE) at approximately 400 mK, triggering discussions around the underlying mechanism and role of moiré effects. In particular, new electron crystal states with non-trivial topology have been proposed. Here, in this paper, we report electrical transport measurements in rhombohedral pentalayer and tetralayer graphene/hBN moiré superlattices at electronic temperatures down to below 40 mK. We observed two more fractional quantum anomalous Hall (FQAH) states and smaller R xx values in pentalayer devices than those previously reported. In the new tetralayer device, we observed FQAHE at moiré filling factors v = 3/5 and 2/3. With a small current at the base temperature, we observed a new extended quantum anomalous Hall (EQAH) state and magnetic hysteresis, where R xy = h/e 2 and vanishing R xx spans a wide range of v from 0.5 to 1.3. At increased temperature or current, EQAH states disappear and partially transition into the FQAH liquid. Furthermore, we observed displacement field-induced quantum phase transitions from the EQAH states to the Fermi liquid, FQAH liquid and the likely composite Fermi liquid. Our observations established a new topological phase of electrons with quantized Hall resistance at zero magnetic field and enriched the emergent quantum phenomena in materials with topological flat bands.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Characterization of the inner edge of the neutron star crust

The poorly known crustal equation of state plays a critical role in many observational phenomena associated with a neutron star. Here using semiclassical Monte Carlo simulations, we explore the possible configurations of the inner edge of the neutron star crust for a variety of baryon densities and proton fractions. Applying the Kirkwood-Buff theory to these two-component systems, we observe how the isothermal compressibility reaches a maximum when isolated nonsymmetric (or quasispherical) clusters are formed in an extremely dilute neutron gas. To determine the neutron fraction, we suggest a geometrical model based on the behavior of the proton-neutron pair correlation function. Accordingly, the equation of state of the inner crust is calculated, illustrating that the nuclear energy in β equilibrium follows a power-law behavior with baryon density. As a possible astrophysical outcome of this study, our results could help refine the mass-radius relation. Finally, our results pave the way towards further investigations of the impact of the proton-neutron pair correlation function on transport properties within the neutron star crust.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Probing Electron Correlations in 1D Electronic Materials Using Single Quantum Channels (Final Technical Report)

This research focuses on the investigation of strong electron correlations in quasi-one-dimensional oxide nanostructures. The archetypical system consisting of heterostructures formed from an ultrathin layer of LaAlO 3 and TiO 2 -terminated SrTiO 3 offers a unique laboratory for probing strong correlations due to intrinsic richness of the physical system and its ability to be modulated at extreme nanoscale dimensions. Under suitable conditions, SrTiO 3 is known to exhibit a variety of important properties associated with strong electronic correlations, especially superconductivity, magnetism. Understanding the nature of these electronic correlations at fundamental scales, and in reduced spatial dimensions, will provide the scientific basis for new energy technologies, information technologies, and materials. The research program is motivated by recent discoveries and technical advances in the PIs' lab, specifically, (1) pioneering work to create conductive nanostructures to be created in LaAlO 3 /SrTiO 3 heterointerfaces with 2 nm spatial resolution, (2) the discovery of a novel correlated electronic phases in which electrons pair without forming a superconductor, and (3) the demonstrated ability to form high-quality quantum point contacts (QPCs) and achieve ballistic electron transport over micrometer-scales in LaAlO 3 /SrTiO 3 nanowires. The program has three major components. First, we have press forward to develop a full microscopic understanding of the mechanism for electron pairing in LaAlO 3 /SrTiO 3 nanostructures. Understanding the “glue” that leads to pairing and superconductivity in a doped semiconductor will help to form a scientific basis for the correlated oxide nanoelectronics platform being developed. Second, we have used their ability to create well-defined quantum channels in nanowires allows to allow novel “scattering” experiments to be performed that probe fundamentally new quasiparticle excitations using an unprecedented combination of experimental methodologies. Experiments already underway with artificially defined superlattices (5 nm period) demonstrate the feasibility of this approach. Third, we have investigated 1D correlations in proximally coupled nanowires, which are predicted exhibit Majorana-like zero modes in the absence of superconductivity or spin-orbit coupling. These studies provide insight not only into new and emergent phenomena in the rich class of complex oxides, but also fundamental mechanisms of superconductivity and magnetism. Thematically, the research aims to combine some of the most important and interesting challenges in the physics of semiconductor nanostructures and correlated electronic materials. The increased understanding of electron correlations can pave the way to the development of new materials with outstanding or novel properties, which could be relevant for energy-related technologies such as low-power high performance computing, quantum computation, or the development of novel materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

A Comparison of GPU-Accelerated Multiphase CFD Solvers on the Polaris Supercomputer: Part 1

This report is in support of the Innovative and Novel Computational Impact on Theory and Experiment (INCITE) program sponsored by the U.S. Department of Energy (USDOE). With INCITE-level resources, one project, titled BubblyFlow, was granted computational resources for the 2025 calendar year on the Polaris supercomputer at the Argonne Leadership Computing Facility (ALCF). The project aims to conduct simulations to understand the fundamental characteristics of turbulent bubbly flow phenomena in nature. Staff at the ALCF and Argonne’s Computational Science division, along with collaborators at the City College of New York and University of Illinois at Chicago, helped a summer student to assess the accuracy and performance of two high performance computing (HPC) codes. Both codes, ImExLBM and FluTAS, are fundamentally different in their mathematical and numerical modeling. However, both may be used to solve the same physical problem. The collaboration sought to better understand the differences between both codes in terms of accuracy and efficiency. This would ultimately help the BubblyFlow project better utilize resources and establish a knowledge-base of code capabilities in future simulation campaigns. We compare ImExLBM and FluTAS, two high-performance multiphase computational fluid dynamics (CFD) solvers, in terms of physical fidelity, time-to-solution, and parallel efficiency. We validate ImExLBM (Implicit-Explicit Lattice Boltzmann Method) against a canonical benchmark and assess it’s performance relative to FluTAS (Fluid Transport Accelerated Solver), a well-established open-source CFD code.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Berry curvature memory through electrically driven stacking transitions

In two-dimensional layered quantum materials, the interlayer stacking order determines both crystalline symmetry and quantum electronic properties such as Berry curvature, topology and electron correlation. Electrical stimuli can strongly influence quasi-particle interactions and the free energy landscape, thus making it possible to access hidden stacking orders with novel quantum properties and enabling dynamic engineering of these attributes. In this paper, we demonstrate electrically driven stacking transitions and a new type of nonvolatile memory based on Berry curvature in few-layer WTe 2 . The interplay of out-of-plane electric fields and electrostatic doping controls in-plane interlayer sliding and creates multiple polar and centrosymmetric stacking orders. In-situ nonlinear Hall transport reveals such stacking rearrangements result in a layer-parity-selective Berry curvature memory in momentum space, where the sign reversal of the Berry curvature only occurs in odd layer crystals. Our findings open an avenue towards exploring coupling between topology, electron correlations, and ferroelectricity in hidden stacking orders and demonstrate a new low energy cost, electrically-controlled topological memory in the atomically thin limit.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Universality of Collective Density Fluctuations in Liquids at and away from Equilibrium – An Integrated Neutron Scattering, Theoretical, and Computational Study (Final Technical Report)

Liquids, ubiquitous on earth, are prototypical disordered condensed matter. As the mediate phase between non-structured gases and ordered solids, liquids p ossess intriguing complexities that are absent on either extreme counter phase, especially when they are driven away from equilibrium. Despite decades of studies, the physics of liquids still keeps surprising us and challenges our understanding of condensed matter, emergent phenomena, and complex systems. Among the many interesting properties, the nature of collective dynamics in liquids remains elusive because of both the presence of strong interactions and the absence of translational invariance. This project studied the universal features of collective density fluctuations (or collective modes), which include both phonon-like excitations and collective relaxations, in the Q-dependent generalized hydrodynamic regime in liquids and liquid-like matter both at and away from equilibrium. We performed systematic studies of the fast collective dynamics (including both phonon-like excitations and collective relaxations) of four classes of liquids characterized by different fragilities and interactions using synergistically integrated coherent neutron scattering experiments (Inelastic Neutron Scattering (INS), Quasi-Elastic Neutron Scattering (QENS), and Neutron Spin Echo (NSE)), a ViscoElastic Hydrodynamics theory we have been developing, andMolecular Dynamics (MD) simulations and our analysis package LiquidLib. The transformative knowledge of the collective dynamics in liquids is important to understand and ultimately control the transport of mass, energy, and charge in liquids, both at and away from thermal equilibrium, which may will likely impact a wide range of disciplines.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Capturing the ground state of uranium dioxide from first principles: Crystal distortion, magnetic structure, and phonons

Uranium dioxide (UO 2 ) remains a formidable challenge for first-principles approaches due to the complex interplay among spin-orbit coupling, Mott physics, magnetic ordering, and crystal distortions. Here we use DFT+U to explore UO 2 at zero temperature, incorporating all the aforementioned phenomena. The technical challenge is to navigate the many metastable electronic states produced by DFT+U, which is accomplished using f-orbital occupation matrix control to search for the ground state. We restrict our search to the high-symmetry ferromagnetic phase, including spin-orbit coupling, which produces a previously unreported occupation matrix. This newfound occupation matrix is then used as an initialization to explore the broken symmetry phases. We find the oxygen cage distortion of the 3k antiferromagnetic state to be in excellent agreement with experiments, and both the spin-orbit coupling and the Hubbard U are critical ingredients. Further, we demonstrate that only select phonon modes have a strong dependence on the Hubbard U, whereas magnetic ordering has only a small influence overall. We perform measurements of the phonon dispersion curves using inelastic neutron scattering, and our calculations show good agreement when using reasonable values of U. The quantitative success of DFT+U warrants exploration of thermal transport and other observables within this level of theory.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Anatomy of Continuous Mars SEIS and Pressure Data from Unsupervised Learning

The seismic noise recorded by the Interior Exploration using Seismic Investigations, Geodesy, and Heat Transport (InSight) seismometer (Seismic Experiment for Interior Structure [SEIS]) has a strong daily quasi-periodicity and numerous transient microevents, associated mostly with an active Martian environment with wind bursts, pressure drops, in addition to thermally induced lander and instrument cracks. That noise is far from the Earth’s microseismic noise. Quantifying the importance of nonstochasticity and identifying these microevents is mandatory for improving continuous data quality and noise analysis techniques, including autocorrelation. Cataloging these events has so far been made with specific algorithms and operator’s visual inspection. We investigate here the continuous data with an unsupervised deep-learning approach built on a deep scattering network. This leads to the successful detection and clustering of these microevents as well as better determination of daily cycles associated with changes in the intensity and color of the background noise. We first provide a description of our approach, and then present the learned clusters followed by a study of their origin and associated physical phenomena. We show that the clustering is robust over several Martian days, showing distinct types of glitches that repeat at a rate of several tens per sol with stable time differences. We show that the clustering and detection efficiency for pressure drops and glitches is comparable to or better than manual or targeted detection techniques proposed to date, noticeably with an unsupervised approach. Finally, here we discuss the origin of other clusters found, especially glitch sequences with stable time offsets that might generate artifacts in autocorrelation analyses. We conclude with presenting the potential of unsupervised learning for long-term space mission operations, in particular, for geophysical and environmental observatories.

58 GEOSCIENCES↗
DOE OSTI2023

Theory and Simulation of Metal–Insulator–Semiconductor (MIS) Photoelectrodes

A metal–insulator–semiconductor (MIS) structure is an attractive photoelectrode-catalyst architecture for promoting photoelectrochemical reactions, such as the formation of H 2 by proton reduction. The metal catalyzes the generation of H 2 using electrons generated by photon absorption and charge separation in the semiconductor. The insulator layer between the metal and the semiconductor protects the latter element from photo-corrosion and, also, significantly impacts the photovoltage at the metal surface. Understanding how the insulator layer determines the photovoltage and what properties lead to high photovoltages is critical to the development of MIS structures for solar-to-chemical energy conversion. Herein, we present a continuum model for charge-carrier transport from the semiconductor to the metal with an emphasis on mechanisms of charge transport across the insulator. The polarization curves and photovoltages predicted by this model for a Pt/HfO 2 /p-Si MIS structure at different HfO 2 thicknesses agree well with experimentally measured data. The simulations reveal how insulator properties (i.e., thickness and band structure) affect band bending near the semiconductor/insulator interface and how tuning them can lead to operation closer to the maximally attainable photovoltage, the flat-band potential. This phenomenon is understood by considering the change in tunneling resistance with insulator properties. The model shows that the best MIS performance is attained with highly symmetric semiconductor/insulator band offsets (e.g., BeO, MgO, SiO 2 , HfO 2 , or ZrO 2 deposited on Si) and a low to moderate insulator thickness (e.g., between 0.8 and 1.5 nm). Beyond 1.5 nm, the density of filled interfacial trap sites is high and significantly limits the photovoltage and the solar-to-chemical conversion rate. These conclusions are true for photocathodes and photoanodes. This understanding provides critical insight into the phenomena enhancing and limiting photoelectrode performance and how this phenomenon is influenced by insulator properties. The study gives guidance toward the development of next-generation insulators for MIS structures that achieve high performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗