SEARCH · Engineering Papers

Results for “phase space dynamics”

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.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Intrinsic thermal conductivity of ZrC from low to ultrahigh temperatures: A critical revisit

Current phonon transport theory based on ground-state calculations has been successful in predicting thermal conductivity at room and medium temperatures but may misrepresent behavior at high temperatures. Here, in this work, we predict the thermal conductivity (κ) of ZrC including electronic and phonon contributions from 300 to 3500 K, by including high-order phonon scattering; lattice expansion; temperature-dependent (TD) second-, third-, and fourth-order force constants (2FC, 3FC, and 4FC); and interband phonon conduction by using first principles. For the phonon transport, we find that four-phonon scattering (4ph) significantly reduces the phonon thermal conductivity (κ ph ), by as much as ~75% at 3500 K. After including 4ph scattering and all other factors, κ ph shows a ~T –1.5 rather than ~T –1 dependence. TD 2FC decreases three-phonon scattering rates but increases 4ph rates by decreasing and increasing the scattering phase spaces, respectively. For 4ph phase space, the TD 2FC flattens phonon bands, and allows more redistribution-4ph processes (1 + 2 → 3 + 4) to happen. The combination effect of TD 2FC and TD 4FC reduces 4ph rates of acoustic modes but increases those of optical modes. The TD 3FC and 4FC decrease the phonon scattering cross section and increase the κ ph significantly (by 52% at 3500 K). The contribution from interband (Wigner) phonon conduction is small, even at ultrahigh temperatures. For electronic thermal transport, we find that it is sensitive to and can be changed by 20% by the TD lattice constants. The Lorenz number varies from 1.6 to 3.3 × 10 –8 WΩK –2 at different temperatures. The theoretical prediction in the literature overpredicts κph (e.g., ~28%) and underpredicts the κ el (e.g., ~38%), resulting in an overall underprediction of κ (~26% at 1500 K). The impacts of grain size and defects are found to be strong, leading to the lower observed thermal conductivity in experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

A Conceptual Design of the Reactor Cavity Cooling System for the Horizontal Compact High Temperature Gas Reactor (HC-HTGR)

Horizontal Compact High Temperature Gas Reactor (HC-HTGR) is being designed by a multi-disciplinary team of nuclear, mechanical, and structural engineers under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed HC-HTGR in 3 years and support its commercialization as a safe, low-cost HTGR. Argonne National Laboratory (Argonne) is responsible for the design and analysis of the reactor cavity cooling system (RCCS) as a safety system for passive decay heat removal of the reactor concept. This report documents the design study to derive a conceptual design study of the RCCS for the HC-HTGR. It includes the identification of the functions and requirements of the HC-HTGR RCCS, design analyses including high-level design consideration and the calculations for optimizing design space of the system with supporting component-level analysis to inform the material selection and performance of the water panel, the description of the conceptual design of the HC-HTGR RCCS derived based on the analyses results, and performance evaluation of the conceptual RCCS for the HC-HTGR. A detailed concept of the RCCS has been identified and high-level system requirements has been developed for the HC-HTGR. Design space focusing on the natural circulation loop portion of the RCCS has been investigated to optimize the system performance. The initial baseline dimensions were firstly derived based on the scoping calculations. A component level design analysis was conducted for the water panel to inform the material selection and to assess its conduction performance. A preliminary system-level performance analysis was performed for the 1/8th of the compartment of the initial baseline design of the RCCS using RELAP5-3D. To improve the system thermal performance, the RCCS design has been updated by exploring various design options by design parametric analyses. Based on the results, the conceptual design of the RCCS for the HC-HTGR has been derived, which satisfies the target performance of ~1 MWt at the elevated vessel wall temperature conditions. Transient simulations were conducted for the conceptual RCCS design for the HC-HTGR under various operation modes and heat load conditions using RELAP5-3D. The system dynamics in different operating states was investigated and the system performance under transients of interest was evaluated. The results demonstrated the overall system feasibility that the RCCS design maintains structures temperatures lower than maximum allowable temperature with sufficient system inventory without any active heat removal in the design process with certain transients addressed. The HC-HTGR RCCS will have additional design updates of subsystems or optimization of the system components during the preliminary and final design phases. Since the entire plant has not been integrated yet, this delivered conceptual design is subject to changes for integration, that require additional conceptual design activities and Quality and Assurance implementation (Q&A). The performance assessment of the RCCS for the HC-HTGR will be then revisited and optimized to finalize the system design, and the RCCS integrated primary system analysis will be utilized to simulate selective accident scenarios of interest where efforts are currently undergoing in the project.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Dynamics of dark matter misalignment through the Higgs portal

A light singlet scalar field feebly coupled through the super-renormalizable Higgs portal provides a minimal and well-motivated realization of ultra-light bosonic dark matter. We study the cosmological production of dark matter in this model by elucidating the dynamics of two sources of scalar field misalignment generated during the radiation era. For large scalar masses (above $\mathcal{O}$ (10 −3 eV)), dark matter is produced through thermal misalignment, by which the scalar field is driven towards large field values as a result of the finite-temperature effective potential. The dominance of thermal misalignment in this mass range leads to a sharp relic abundance prediction which is, to a significant extent, insensitive to the initial conditions of the scalar field. On the other hand, for low mass scalars (below $\mathcal{O}$ (10 −5 eV)), dark matter is produced via VEV misalignment, which is caused by the induced scalar field vacuum expectation value triggered by the electroweak phase transition. We show that the relic abundance in this low mass range is sensitive to the scalar field initial conditions. In the intermediate mass range, the relic abundance is a consequence of a competition between thermal misalignment and VEV misalignment, which can potentially lead to novel forced resonance effects which cause a recurring enhancement and suppression in the late time oscillation amplitude as a function of the scalar mass. We compare our relic abundance predictions with constraints and projections from equivalence principle and inverse square law tests, stellar cooling, resonant molecular absorption, and observations of extra-galactic background light and diffuse X-ray backgrounds. New experimental ideas are needed to probe most of the cosmologically motivated regions of parameter space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Learning with Delayed Rewards—A Case Study on Inverse Defect Design in 2D Materials

Defect dynamics in materials are of central importance to a broad range of technologies from catalysis to energy storage systems to microelectronics. Material functionality depends strongly on the nature and organization of defects-their arrangements often involve intermediate or transient states that present a high barrier for transformation. The lack of knowledge of these intermediate states and the presence of this energy barrier presents a serious challenge for inverse defect design, especially for gradient-based approaches. Here, we present a reinforcement learning (RL) [Monte Carlo Tree Search (MCTS)] based on delayed rewards that allow for efficient search of the defect configurational space and allows us to identify optimal defect arrangements in low-dimensional materials. Using a representative case of two-dimensional MoS 2 , we demonstrate that the use of delayed rewards allows us to efficiently sample the defect configurational space and overcome the energy barrier for a wide range of defect concentrations (from 1.5 to 8% S vacancies)-the system evolves from an initial randomly distributed S vacancies to one with extended S line defects consistent with previous experimental studies. Detailed analysis in the feature space allows us to identify the optimal pathways for this defect transformation and arrangement. Comparison with other global optimization schemes like genetic algorithms suggests that the MCTS with delayed rewards takes fewer evaluations and arrives at a better quality of the solution. The implications of the various sampled defect configurations on the 2H to 1T phase transitions in MoS 2 are discussed. In this study, we introduce a RL strategy employing delayed rewards that can accelerate the inverse design of defects in materials for achieving targeted functionality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

The kinematic properties of Milky Way stellar halo populations

ABSTRACT In the Gaia era, stellar kinematics are extensively used to study Galactic halo stellar populations, to search for halo structures, and to characterize the interface between the halo and hot disc populations. We use distribution function-based models of modern data sets with 6D phase space data to qualitatively describe a variety of kinematic spaces commonly used in the study of the Galactic halo. Furthermore, we quantitatively assess how well each kinematic space can separate radially anisotropic from isotropic halo populations. We find that scaled action space (the ‘action diamond’) is superior to other commonly used kinematic spaces at this task. We present a new, easy to implement selection criterion for members of the radially anisotropic Gaia–Enceladus merger remnant. Assuming a 1:1 ratio of Gaia–Enceladus stars to more isotropic halo, we find our selection achieves a sample purity of 86 per cent in our models with respect to contamination from the more isotropic halo. We compare this criterion to literature criteria, finding that it produces the highest purity in the resulting samples, at the expense of a modest reduction in completeness. We also show that selection biases that underlie nearly all contemporary spectroscopic data sets can noticeably impact the E−Lz distribution of samples in a manner that may be confused for real substructure. We conclude by providing recommendations for how authors should use stellar kinematics in the future to study the Galactic stellar halo.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Analysis of Two Models for the Angular Structure of the Outflows Producing the Swift/XRT “Larger-angle Emission” of Gamma-Ray Bursts

The quasi-instantaneous emission from a relativistic surface endowed with a Lorentz factor that decreases away from the outflow symmetry axis can naturally explain the three phases observed by Swift X-Ray Telescope (XRT) in gamma-ray bursts (GRBs) and their afterglows (GRB tail, afterglow plateau, and postplateau) based only on the angular change of the relativistic Doppler boost across the outflow surface. We develop further the analytical formalism of the “larger-angle emission” model for the case of “n-exponential” outflows (where the Lorentz factor Γ dependence of the angular location θ is Γ ~ exp{-(θ/θ c ) n }), and compare its ability to account for the X-ray emission of XRT afterglows relative to that of “power-law” outflows (Γ ∼ θ −g ). Power-law outflows yield longer afterglow plateaus, followed by slower postplateau flux decays than n-exponential outflows, features which may be used in identifying which type of angular structure is at work in a given afterglow. Identifying the Γ(θ) angular structure that accommodates XRT light curves is slightly complicated by the fact that the afterglow X-ray light curve is also determined by how two characteristics of the comoving-frame emission spectrum (peak energy $E'_p$ and peak intensity $i'_p$) change with the angular location or, equivalently, with the Lorentz factor. Here, we assume power-law Γ dependences of those spectral characteristics and find that, unlike power-law outflows, n-exponential outflows cannot account for plateaus with a temporal dynamical range larger than 100 (2 dex in logarithmic space). To capture all the information contained in XRT afterglow measurements (0.3–10 keV unabsorbed flux and effective spectral slope), we calculate 0.3 and 10 keV light curves using a broken-power-law emission spectrum of peak energy and low- and high-energy slopes that are derived from the effective slope measured by XRT. This economical peak energy determination is found to be consistent with the results of more expensive spectral fits. The angular distributions of the Lorentz factor, comoving frame peak energy, and peak intensity (Γ(θ), $E'_p$(θ), $i'_p$(θ)) constrain the (yet-to-be determined) convolution of various features of the production of relativistic jets by solar-mass black holes and of their propagation through the progenitor/circumburst medium, while the $E'_p$(Γ) and $i'_p$(Γ) dependences may constrain the GRB dissipation mechanism and the GRB emission process.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Krylov complexity in mixed phase space

We investigate the Krylov complexity of thermofield double states in systems with mixed phase space, uncovering a direct correlation with the Brody distribution, which interpolates between Poisson and Wigner statistics. Our analysis spans two-dimensional random matrix models featuring (I) GOE-Poisson and (II) GUE-Poisson transitions and extends to higher-dimensional cases, including a stringy matrix model (GOE-Poisson) and the mass-deformed SYK model (GUE-Poisson). Krylov complexity consistently emerges as a reliable marker of quantum chaos, displaying a characteristic peak in the chaotic regime that gradually diminishes as the Brody parameter approaches zero, signaling a shift toward integrability. These results establish Krylov complexity as a powerful diagnostic of quantum chaos and highlight its interplay with eigenvalue statistics in mixed phase systems.

chaos & nonlinear dynamics↗

DOE OSTI2023

On the interaction of a wind turbine wake with a conventionally neutral atmospheric boundary layer

In this work, we investigate the dynamics of wind turbine tip-vortex breakdown in a conventionally neutral atmospheric boundary layer (ABL). To this end, high-resolution data are collected from large-eddy simulations of a wind turbine operating within a neutral ABL and studied by means of proper orthogonal decomposition (POD) and Fourier analysis. The high resolution of the generated data in both space and time allows us to gain insight into the tip-vortex breakdown mechanisms by (i) capturing the energy modes of the coherent structures, (ii) studying their contribution to the tip-vortex breakdown through their power spectra functions and mean kinetic energy (MKE) flux, and (iii) analysing the growth rate of each contributing perturbation frequency along tip vortices. Our analysis shows that under a fully turbulent scenario, the growth rate of perturbations along the tip vortices is largest for low wave numbers, i.e. long-wave perturbations. Additionally, the MKE flux reaches its highest value at two diameters downstream of the rotor plane, a behaviour that can be attributed to the coexistence of multiple interacting POD modes, with the streamwise vortex roller mode being the primary contributor to the total MKE flux budget, contributing approximately . Finally, comparisons with a laminar, uniform flow scenario subject to a single-frequency perturbation highlight the differences between the two ambient flow conditions. In the non-turbulent, uniform flow scenario, the growth rate attains its maximum value at a wave number corresponding to the out-of-phase mutual-inductance mechanism, whereas the MKE flux exhibits local minima and maxima along the wake and at different downstream locations depending on the perturbation frequency. Our analyses suggest that the breakdown of the wind turbine tip vortices under a fully turbulent neutral ABL inflow is due to complex interactions across a range of excitation frequencies, in which the mutual-inductance instability may not be the dominant one.

17 WIND ENERGY↗

DOE OSTI2021

Robust Molecular Predictive Methods for Novel Polymer Discovery and Applications

Polymeric materials are ubiquitous in modern society and they play an instrumental role in almost all industries, undoubtedly including the energy and environment sectors. Increased demand of energy and awareness to sustainability both necessitates the development of novel polymers with enhanced properties. Unfortunately, their structural and behavioral complexity render such discovery challenging and impeded. To address this problem, scientists are developing various computational modeling techniques and leveraging their power to depict the relationship between structural characteristics of polymers and their properties (such as rheological behaviors), and use such prediction to guide the design and syntheses of novel polymeric materials with enhanced performances. Unfortunately, predicting the relationships between polymer structure and composition with rheological properties via atomistic modeling is still a major challenge because of the extended time and length scales involved. Studying dynamic shear viscosity and linear viscoelasticity using molecular models requires capabilities that have been elusive, including representation of large molecular weight chains with an effective internal scale capable of describing entanglement, shear-rates that are in the s-1 scale with accurate quantitative stresses, and chemically-realistic combinations of both homogeneous and heterogeneous systems. Motivated by these unmet challenges, the overall technical objective of this DOE-STTR Phase II project is to develop robust molecular predictive methods for advanced polymer discovery and applications and especially for designing and demonstrating the “smart” polymer-based waterflooding enhanced oil recovery (EOR) process. In particular, we apply state-of-the-art molecular modeling methods developed by our academic partner, Materials Stimulation Center (MSC) at California Institute of Technology (Caltech), to facilitate and accelerate the experimental discovery processes. During the Phase I of this project, we had focused on development and demonstration of the molecular modeling methods to describe rheological properties of non-Newtonian polymer fluids, and to improve our fundamental understandings of shear-thickening mechanism and kinetics. In Phase II, we further apply the theoretical models to guide our experimental programs to improve our design of smart rheology modifier (SRM) polymers and their optimization for EOR. Specifically, we have three objectives in the Phase II study: (1) to further improve out computational modeling methods, coupling with the advanced machine learning algorithms; (2) to develop cost-effective and efficient SRM-flooding process suitable for EOR applications under typical reservoir conditions; and (3) to further explore the application of our molecular predictive models for innovative material discovery in other industrial applications. The recent development of our multiscale predictive framework allows the successful prediction of rheological properties from the chemical structure for polymers of experimentally relevant molecular weights, and provides an in-silico machine learning engine for screening novel compositions and structures with optimized non-Newtonian response, required for both shear-thinning and shear-thickening applications. Our framework provides: (1) procedures and tools for systematic coarsening from atomistic models and reverse mapping of coarse-grain models to atomistic, (2) unique ab initio methods to characterize the atomistic origin of colloidal and interfacial interactions and phenomena, (3) systematic structure and composition builders based on practical descriptors that drive rheological changes in polymer melts and diluted polymer mixtures, (4) a rheological properties engine capable of predicting viscosity in the zero-shear limit and under realistic dynamic conditions (for shear-rates commensurate with experiments) for large heterogeneous systems, (5) coarse-grain force fields with improved non-bond descriptions based on accurate quantum mechanics, (6) an in-silico screening machine learning engine that feeds from the systematic model builders to cover the descriptors search space, computes the rheological properties from converged trajectories spanning sub-milliseconds and ranks them for each structure/composition using an automated viscosity-vs-shear rate fitness function that can be tuned for shear-thickening, shear-thinning and other rheological responses.

02 PETROLEUM↗

DOE OSTI2025

Rethinking 𝛼−RuCl 3 : Parameters, models, and phase diagram

RuCl 3 was likely the first ever deliberately synthesized ruthenium compound, following the discovery of the 44 Ru element in 1844. For a long time it was known as an oxidation catalyst, with its physical properties being discrepant and confusing, until a decade ago when its allotropic form 𝛼−RuCl 3 rose to exceptional prominence. This “rediscovery” of 𝛼−RuCl 3 has not only reshaped the hunt for a material manifestation of the Kitaev spin liquid, but it has opened the floodgates of theoretical and experimental research in the many unusual phases and excitations that the anisotropic-exchange magnets as a class of compounds have to offer. Given its importance for the field of Kitaev materials, it is astonishing that the low-energy spin model that describes this compound and its possible proximity to the much-desired spin-liquid state is still a subject of significant debate ten years later. In the present study, we argue that the existing key phenomenological observations put strong natural constraints on the effective microscopic spin model of 𝛼−RuCl 3 , and specifically on its spin-orbit-induced anisotropic-exchange parameters that are responsible for the nontrivial physical properties of this material. These constraints allow one to focus on the relevant region of the multidimensional phase diagram of the 𝛼−RuCl 3 model, suggest an intuitive description of it via a different parametrization of the exchange matrix, offer a unifying view on the earlier assessments of its parameters, and bring closer together several approaches to the derivation of anisotropic-exchange models. We explore extended phase diagrams relevant to the 𝛼−RuCl 3 parameter space using quasiclassical, Luttinger-Tisza, exact diagonalization, and density-matrix renormalization-group methods, demonstrating a remarkably close quantitative accord between them on the general structure and hierarchy of the phases, with the zigzag, ferromagnetic, and incommensurate phases that are proximate to each other. As a result, one of the highlights is the detailed agreement on the nature of the incommensurate phases that realize two distinct counterrotating helical states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Dynamics of almost strong edge modes in spin chains away from integrability

Results are presented for the dynamics of an almost strong edge mode which is the quasistable Majorana edge mode occurring in nonintegrable spin chains. The dynamics of the edge mode is studied using exact diagonalization, and compared with time evolution with respect to an effective semi-infinite model in Krylov space obtained from the recursion method. The effective Krylov Hamiltonian is found to resemble a spatially inhomogeneous Su-Schrieffer-Heeger model where the hopping amplitude increases linearly with distance into the bulk, typical of thermalizing systems, but also has a staggered or dimerized structure superimposed on it. Here, the nonperturbatively long lifetime of the edge mode is shown to be due to this staggered structure which diminishes the effectiveness of the linearly growing hopping amplitude. On taking the continuum limit of the Krylov Hamiltonian, the edge mode is found to be equivalent to the quasistable mode of a Dirac Hamiltonian on a half line, with a mass which is nonzero over a finite distance, before terminating into a gapless metallic bulk. The analytic estimates are found to be in good agreement with the numerically obtained lifetimes of the edge mode.

1-dimensional systems↗

DOE OSTI2020

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗

DOE OSTI2023

Understanding Impacts of Frequency Calculations on Underfrequency Load Shedding

Power system fundamental frequency-dependent protection decisions and control decisions are increasingly common in the distribution system space. Understanding the algorithms used in frequency measurements is critical to understanding the decisions made by the many intelligent electronic devices in the power system. These decisions are critically important for reliability analysis for interconnected power systems as they bear directly on resource planning to mitigate under-frequency conditions that would result in under-frequency load shedding (UFLS). This work aims to improve the understanding of frequency measurement algorithms, evaluate the intelligent electronic devices that use frequency measurements for decision making, and understand the algorithms' direct impact on frequency related decision making. First, this work presents background information on the use of frequency measurements in protection logic, specifically UFLS. Second, the paper presents a relay hardware evaluation test bed used to detect and protect systems from under-frequency events. Finally, the paper presents the dynamic events used to evaluate commercially available, off-the-shelf relay equipment and the results from the relay evaluation.

COMTRADE↗

DOE OSTI2023

Density variance dynamics in disparate shock tubes

This report discusses two experiments which investigate thin layer, heavy curtain fragmentation from the perspective of a Reynolds-Averaged mix model, in drastically disparate experimental regimes. The first, the centimeter/millisecond-scale “Horizontal Shock Tube” (HST) is a compressed-gas piston-driven shock tube experiment. The second, the “Multishock thin layer” (Mshock) experiment performed at the National Ignition Facility, is a micrometer/nanosecond-scale laser-driven shock tube experiment. Both are situations in which a heavy plane layer (a ‘curtain’) is initially suspended in a lighter medium. After being shocked from at least one side, the layer translates while its interfaces evolve due to the excitation of the Richtmyer-Meshkov instability at its surfaces. The evolution of density variance, which initially exists only on the surface of the layer, as it comes to encompass the whole layer interior is used as a description of layer fragmentation and dissolution. These experiments have each been simulated in the Los Alamos National Laboratory multi-physics code xRAGE, which includes fundamental hydrodynamics, extended plasma physics and radiation effects which are important to drive the high-energy density experiment, and the Besnard-Harlow-Rauenzahn (BHR) turbulence model. In each, the principal diagnostic for comparison is an experimental metric for the density (co)variance, b, which tracks the moments of the density field at the curtain interfaces and body. Due to experimental constraints in different regimes (i.e. optical diagnostics can be deployed on conventional shock tubes, while the plasma shock tubes must be imaged by x-rays; interfaces can be imposed to specification on laser-driven experiments, which are stored in the solid phase, while conventional experiments have imperfect control of the flow fields which separate the layer, etc.) the experiments are not designed to be perfect scaled cognates of one another. However, despite the separation of six orders of magnitude of scaling in time, and four in space, we are able to demonstrate that the same turbulence model, operating in the same fashion in the same computer code, is able to reproduce results in each experiment, by tracking evolution due to common relevant physics. Additionally, we will present preliminary work toward density variance comparisons in a single-interface Richtmyer-Meshkov configuration, the conventional fluid “Vertical Shock Tube” (VST) experiment, and the Modal Initial Conditions (ModCons) campaign fielded at the OMEGA-EP laser facility.

42 ENGINEERING↗

DOE OSTI2023

Robust clustering of the local Milky Way stellar kinematic substructures with Gaia eDR3

Understanding local stellar kinematic substructures in the solar neighbourhood helps build a complete picture of the formation of the Milky Way, as well as an empirical phase space distribution of dark matter that would inform detection experiments. We apply the clustering algorithm HDBSCAN on the Gaia early third data release to identify a list of stable clusters in velocity space and action-angle space by taking into account the measurement uncertainties and studying the stability of the clustering results. We find 1405 (497) stars in 23 (6) robust clusters in velocity space (action-angle space) that are consistently not associated with noise. We discuss the kinematic properties of these structures and study whether many of the small clusters belong to a similar larger cluster based on their chemical abundances. They are attributed to the known structures: the Gaia Sausage-Enceladus, the Helmi Stream, and globular cluster NGC 3201 are found in both spaces, while NGC 104 and the thick disc (Sequoia) are identified in velocity space (action-angle space). Although we do not identify any new structures, we find that the HDBSCAN member selection of already known structures is unstable to input kinematics of the stars when resampled within their uncertainties. We therefore present the stable subset of local kinematic structures, which are consistently identified by the clustering algorithm, and emphasize the need to take into account error propagation during both the manual and automated identification of stellar structures, both for existing ones as well as future discoveries.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Crystallization of binary nanocrystal superlattices and the relevance of short-range attraction

The synthesis of binary nanocrystal superlattices (BNSLs) enables the targeted integration of orthogonal physical properties, such as photoluminescence and magnetism, into a single superstructure, unlocking a vast design space for multifunctional materials. However, the formation mechanism of BNSLs remains poorly understood, restricting the prediction of the structure and properties of superlattices. Here we use a combination of in situ scattering and molecular simulation to elucidate the self-assembly of two common BNSLs (AlB 2 and NaZn 13 ) through emulsion templating. Our self-assembly experiments reveal that no intermediate structures precede the formation of the final binary phases, indicating that their formation proceeds through classical nucleation. Using simulations, we find that, despite the formation of AlB 2 and NaZn 13 typically being attributed to entropy, their self-assembly is most consistent with the nanocrystals possessing short-range interparticle attraction, which we find can accelerate nucleation kinetics in BNSLs. Here, we also find homogeneous, classical nucleation in simulations, corroborating our experiments. These results establish a robust correspondence between experiment and theory, paving the way towards prediction of BNSLs.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Same-sided successive-shock HED instability experiments

Inertial confinement fusion (ICF) and high-energy density (HED) physics experiments experience complicated forcing for instability growth and mix due to the ubiquitous presence of multiple shocks interacting with perturbations on multiple material interfaces. One common driver of instability growth is successive shocks from the same direction. However, there is a severe lack of analytic work and modeling validation for same-sided successive shocks since they are extremely difficult to achieve with conventional (non-HED) drivers. Successive shocks access a large instability parameter space; idealized fluid theory [K. O. Mikaelian, Phys. Rev. A 31, 410 (1985)] predicts 15 different interface evolution scenarios for a sinusoidal perturbation. Growth becomes more complex for multi-mode, compressible HED systems. The Mshock campaign is the first experiment in any fluid regime to probe a wide portion of successive shock parameter space. This is enabled by our development of a hybrid direct/indirect drive platform capable of creating independently controllable successive shocks on the National Ignition Facility. These experiments have delivered the first data capable of rigorously challenging our models and their ability to accurately capture Richtmyer–Meshkov growth under successive shocks. Single-mode and two-mode experiments have successfully demonstrated the ability to access and control the various growth scenarios of the shocked interface, including re-inversion, freeze out, and continued growth. Simulations and theoretical modeling are shown to accurately capture the experimental observations in the linear growth phase, giving us confidence in our ICF/HED design codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

(2+1)-dimensional compact Lifshitz theory, tensor gauge theory, and fractons

The (2+1)-dimensional continuum Lifshitz theory of a free compact scalar field plays a prominent role in a variety of quantum systems in condensed matter physics and high energy physics. It is known that in compact space, it has an infinite ground-state degeneracy. In order to understand this theory better, we consider two candidate lattice regularizations of it using the modified Villain formalism. We show that these two lattice theories have significantly different global symmetries (including a dipole global symmetry), anomalies, ground-state degeneracies, and dualities. In particular, one of them is self-dual. Given these theories and their global symmetries, we can couple them to corresponding gauge theories. Here, these are two different U(1) tensor gauge theories. The resulting models have excitations with restricted mobility, i.e., fractons. Finally, we give an exact lattice realization of the fracton and lineon-elasticity dualities for the Lifshitz theory, and scalar and vector charge gauge theories.

2-dimensional systems↗