Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum mechanics”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Field theories with a vector global symmetry

Motivated by recent discussions of fractons, we explore nonrelativistic field theories with a continuous global symmetry, whose charge is a spatial vector. We present several such symmetries and demonstrate them in concrete examples. They differ by the equations their Noether currents satisfy. Simple cases, other than the translation symmetry, are an ordinary (relativistic) one-form global symmetry and its nonrelativistic generalization. In the latter case the conserved charge is associated with a codimension-one spatial manifold, but it is not topological. More general examples involve charges that are integrated over the entire space. We also discuss the coupling of these systems to gauge fields for these symmetries. We relate our examples to known continuum and lattice constructions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Coarse-grained methods for heterogeneous vesicles with phase-separated domains: Elastic mechanics of shape fluctuations, plate compression, and channel insertion

We develop coarse-grained particle approaches for studying the elastic mechanics of vesicles with heterogeneous membranes having phase-separated domains. We perform simulations both of passive shape fluctuations and of active systems where vesicles are subjected to compression between two plates or subjected to insertion into narrow channels. Analysis methods are developed for mapping particle configurations to continuum fields with spherical harmonics representations. Heterogeneous vesicles are found to exhibit rich behaviors where the heterogeneity can amplify surface two-point correlations, reduce resistance during compression, and augment vesicle transport times in channels. The developed methods provide general approaches for characterizing the mechanics of coarse-grained heterogeneous systems taking into account the roles of thermal fluctuations, geometry, and phase separation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Complications in the interpretation of the charge-asymmetry-dependent π flow for the chiral magnetic wave

The charge asymmetry ($A_ch$) dependence of the $π^₋$ and $π^+$ elliptic flow difference, $Δv_2(A_{\text{ch}})$, has been regarded as a sensitive observable for the possible chiral magnetic wave (CMW) in relativistic heavy ion collisions. In this work, we first demonstrate that, due to nonflow backgrounds, the flow measurements by the Q-cumulant method using all charged particles as reference introduce a trivial linear term to $Δv_2(A_{\text{ch}})$. The trivial contribution to the slope of triangle flow difference $Δv_3( A_{\text{ch}})$ can be negative if the nonflow is dominated by back-to-back pairs. After eliminating the trivial term, we find that the nonflow between like-sign pairs gives rise to an additional positive slope to $Δv_2(A_{\text{ch}})$ because of the larger multiplicity dilution to $π^+$ ($π^₋$) at positive (negative) $A_{\text{ch}}$. We further find that the competition between different π sources can introduce another nontrivial linear-$A_{\text{ch}}$ term due to their different multiplicity fluctuations and anisotropic flows. We then study the effect of neutral cluster (resonance) decays as a mechanism for local charge conservation on the slope parameter of $Δv_2(A_{\text{ch}})$. We find that the slope parameter is sensitive to the kinematics of those neutral clusters. Light resonances give positive slopes while heavy resonances give negative slopes. Local charge conservation from continuum cluster mass distribution can give a positive slope parameter comparable to experimental data. Our studies indicate that many non-CMW physics mechanisms can give rise to a $A_{\text{ch}}$-dependent $Δv_2(A_{\text{ch}})$ and the interpretation of $Δv_2(A_{\text{ch}})$ in terms of the CMW is delicate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Modeling the formation of Sedan Crater using the FLAG and HOSS codes

Numerical modeling of explosion crater formation requires accounting for complex physical processes. Numerical validation of explosion cratering is an important step in modeling and requires experimental data for comparison. Models using discrete elements and continuum models have both benefits and drawbacks to their approaches. In this work, we consider both an arbitrary Lagrangian–Eulerian (ALE) hydrocode and a finite discrete element method (FDEM) approach to modeling the formation of the Sedan crater, the largest human-made crater in the United States. The Sedan crater formed from an underground nuclear detonation in the Nevada desert as part of Project Plowshare. Our models show that the continuum approach of the hydrocode matched well compared to early test time prior to the mound rupture and subsequent fireball venting, when most of the alluvium exhibited fluid behavior. Our FDEM approach matched the final crater dimensions well, after material had settled back into the crater, when material strength and solid mechanics play key roles. Our work shows how leveraging the benefits of multiple numerical approaches can lead to better understanding of complex physical problems, especially problems with limited experimental data. By using a continuum approach to early-time hydrodynamics and an FDEM approach to later-time solid mechanics, we can better understand the different physical regimes of explosion crater formation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A data-driven multiscale model for reactive wetting simulations

Here, we describe a data-driven, multiscale technique to model reactive wetting of a silver–aluminum alloy on a Kovar™ (Fe-Ni-Co alloy) surface. We employ molecular dynamics simulations to elucidate the dependence of surface tension and wetting angle on the drop’s composition and temperature. A design of computational experiments is used to efficiently generate training data of surface tension and wetting angle from a limited number of molecular dynamics simulations. The simulation results are used to parameterize models of the material’s wetting properties and compute the uncertainty in the models due to limited data. The data-driven models are incorporated into an engineering-scale (continuum) model of a silver–aluminum sessile drop on a Kovar™ substrate. Model predictions of the wetting angle are compared with experiments of pure silver spreading on Kovar™ to quantify the model-form errors introduced by the limited training data versus the simplifications inherent in the molecular dynamics simulations. The paper presents innovations in the determination of “convergence” of noisy MD simulations before they are used to extract the wetting angle and surface tension, and the construction of their models which approximate physio-chemical processes that are left unresolved by the engineering-scale model. Together, these constitute a multiscale approach that integrates molecular-scale information into continuum scale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Γ valley transition metal dichalcogenide moiré bands

Significance At small twist angles, moiré patterns form in transition metal dichalcogenide bilayers. Using sophisticated ab initio calculations, we develop a model able to describe the band structure of these systems. We found that the moiré pattern acts as a long-wavelength modulating potential, with a sixfold rotational symmetry not present in the underlying lattice. The modulation potential reduces hole band widths and localizes electrons or holes around honeycomb or kagome lattice sites.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Investigation of mode activity in NBI-heated experiments of Wendelstein 7-X

The 2018 operation phase (OP 1.2b) of the stellarator Wendelstein 7-X (W7-X) included, for the first time, neutral beam injection (NBI) to heat the plasma. Since the injection geometry at W7-X is not parallel, this generates both passing and trapped fast particles. During longer phases of NBI injection, with the primary purpose to study the heating efficiency of this system, Alfvén eigenmodes (AEs) were observed by a number of diagnostics, including the phase contrast imaging (PCI) system, the magnetic pick-up coils (Mirnov coils), and the soft x-ray multi-camera tomography system (XMCTS). Alfvén eigenmodes are of great interest for future fusion reactors as it has been shown that the resonant interaction of fast ions with self-excited AEs can lead to enhanced transport of fast ions and potentially to energy losses. This is especially true for so-called gap-modes, Alfvén eigenmodes with frequencies in gaps of the continuous spectrum, since they lack continuum damping. These modes are commonly known to be excited by fast ions, but other destabilizing mechanisms, e.g. the electron-pressure gradient are also possible. In this article we present a first analysis of the experimentally observed frequencies from the theoretical side. The calculation of shear Alfvén wave continua for selected cases and the assignment of observed frequencies to the gaps of the continuous spectra are presented. Using the ideal-MHD code CKA (Könies A. 2007 10th IAEA TM on Energetic Particles in Magnetic Confinement System), we find gap modes that match the experimental measurements in terms of the observed frequencies. We emphasize the crucial roles played by the coupling of sound and Alfvén waves as well as of the Doppler shift arising as a consequence of the radial electric field in W7-X. We employ the perturbative gyrokinetic code CKA-EUTERPE (Feh´er 2013 Simulation of the interaction between Alfv´en waves and fast particles), using a slowing-down distribution function for the fast ions as calculated by the Monte-Carlo particle following code ASCOT (Hirvijoki et al 2014 Comput. Phys. Commun. 185 1310–21) to assess the fast-ion drive. Here we find that the fast-ion drive is insufficient to overcome the background-plasma damping. The fact that unstable modes were observed experimentally may point to problems with the modelling or indicate the existence of other destabilizing mechanisms, e.g. associated with the electron-pressure gradient (Windisch et al 2017 Plasma Phys. Control. Fusion 59 105002) that sensitively depend on the profiles of the background plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quantum Spin Dynamics Due to Strong Kitaev Interactions in the Triangular-Lattice Antiferromagnet CsCeSe 2

The extraordinary properties of the Kitaev model have motivated an intense search for new physics in materials that combine geometrical and bond frustration. In this Letter, we employ inelastic neutron scattering, spin wave theory, and exact diagonalization to study the spin dynamics in the perfect triangular-lattice antiferromagnet (TLAF) CsCeSe 2 . This material orders into a stripe phase, which is demonstrated to arise as a consequence of the off-diagonal bond-dependent terms in the spin Hamiltonian. By studying the spin dynamics at intermediate fields, we identify an interaction between the single-magnon state and the two-magnon continuum that causes decay of coherent magnon excitations, level repulsion, and transfer of spectral weight to the continuum that are controlled by the strength of the magnetic field. Our results provide a microscopic mechanism for the stabilization of the stripe phase in TLAF and show how complex many-body physics can be present in the spin dynamics in a magnet with strong Kitaev coupling even in an ordered ground state. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Spontaneously broken boosts and the Goldstone continuum

The spontaneous breaking of boost invariance is ubiquitous in nature, yet the associated Goldstone bosons are nowhere to be seen. We discuss why some subtleties are to be expected in the Goldstone phenomenon for spontaneously broken boosts, and derive the corresponding quantum mechanical, non-perturbative Goldstone theorem. Despite similarities with more standard Goldstone theorems, we show by examples that ours can be obeyed by quite unusual spectra of low-energy excitations. In particular, for non- relativistic Fermi liquids, we prove that it is obeyed by the particle-hole continuum . To the best of our knowledge, this is the first example of a Goldstone theorem obeyed by a continuum rather than by (approximately stable) single-particle Goldstone boson states in an interacting theory at zero temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Insights into the Li Diffusion Mechanism in Si/C Composite Anodes for Lithium-Ion Batteries

Recently, Si/C composite materials have attracted enormous research interest as the most promising candidates for the anodes of next-generation lithium-ion batteries, owing to their high energy density and mechanical buffering property. However, the fundamental mechanism of Li diffusion behavior in various Si/C composite materials remains unclear, with our understanding limited by experimental techniques and continuum modeling methodologies. Herein, the atomic behavior of Li diffusion in the Si/C composite material is studied within the framework of density functional theory. Two representative structural mixing formats, that is, simple mixture mode and core-shell mode, are modeled and compared. We discover that the carbon material increases Li diffusion in silicon from 7.75 x 10 -5 to 2.097 x 10 -4 cm 2 /s. The boost is about 50% more obvious in the mixture mode, while the core-shell structure shows more dependence on the atomic structures of the carbon layer. These results offer new insights into Li diffusion behavior in Si/C composites and unlock the enhancing mechanism for Li diffusion in Si/C. Furthermore, this understanding facilitates the modeling of batteries with composite anodes and will guide the corresponding structure designs for robust and high-energy-density batteries.

25 ENERGY STORAGE↗
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↗
DOE OSTI2023

Results of the Micromorphic Upscaling Workflow for the PSAAP III Year 3 Report

Predicting the mechanical response and failure of heterogeneous materials has proven difficult. Multiscale numerical methods based in higher order continuum theories attempt to bridge the gap between microscale and macroscale structural behavior. Micromorphic continuum theories have shown promise. The Tardigrade software package is an implementation of Eringen’s micromorphic continuum theory with capabilities to support multiscale material modeling workflows. These include homogenization through the Micromorphic Filter, calibration of micromorphic material models, and macroscale simulation in Tardigrade-MOOSE. This work discusses micromorphic upscaling efforts of the University of Colorado Boulder PSAAP III multidisciplinary simulation center (MSC). Verification studies are presented that compare the accuracy of the upscaling workflow with analytical solutions for a trivial stress state and homogeneous material using direct numerical simulations (DNS) conducted in the Ratel finite element method (FEM) and GEOS material point method (MPM) codes. These verification studies consider upscaling using the Micromorphic Filter for a “single filter domain”. Finally, DNS of a heterogeneous composite material is upscaled using a “multiple filter domain” method.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Physics of the Coronal-line Region for Galaxies in Mapping Galaxies at Apache Point Observatory

The fundamental nature and extent of the coronal-line region (CLR), which may serve as a vital tracer for active galactic nucleus (AGN) activity, remain unresolved. Previous studies suggest that the CLR is produced by AGNdriven outflows and occupies a distinct region between the broad-line region and the narrow-line region, which places it tens to hundreds of parsecs from the galactic center. Here, we investigate 10 coronal line (CL; ionization potential ≳100 eV) emitting galaxies from the Sloan Digital Sky Survey IV Mapping Galaxies at Apache Point Observatory (MaNGA) catalog with emission from one or more CLs detected at ≥5$σ$ above the continuum in at least 10 spaxels—the largest such MaNGA catalog. We find that the CLR is far more extended, reaching out to 1.3–23 kpc from the galactic center. We crossmatch our sample of 10 CL galaxies with the largest existing MaNGA AGN catalog and identify seven in it; two of the remaining three are galaxy mergers and the final one is an AGN candidate. Further, we measure the average CLR electron temperatures as ranging between 12,331 and 22,530 K, slightly above the typical threshold for pure AGN photoionization (~20,000 K) and indicative of shocks (e.g., merger induced or from supernova remnants) in the CLR. In conclusion, we reason that ionizing photons emitted by the central continuum source (i.e., AGN photoionization) primarily generate the CLs, and that energetic shocks are an additional ionization mechanism that likely produce the most extended CLRs we measure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Protonic solid-state electrochemical synapse for physical neural networks

Physical neural networks made of analog resistive switching processors are promising platforms for analog computing. State-of-the-art resistive switches rely on either conductive filament formation or phase change. These processes suffer from poor reproducibility or high energy consumption, respectively. Herein, we demonstrate the behavior of an alternative synapse design that relies on a deterministic charge-controlled mechanism, modulated electrochemically in solid-state. The device operates by shuffling the smallest cation, the proton, in a three-terminal configuration. It has a channel of active material, WO 3 . A solid proton reservoir layer, PdH x , also serves as the gate terminal. A proton conducting solid electrolyte separates the channel and the reservoir. By protonation/deprotonation, we modulate the electronic conductivity of the channel over seven orders of magnitude, obtaining a continuum of resistance states. Proton intercalation increases the electronic conductivity of WO 3 by increasing both the carrier density and mobility. This switching mechanism offers low energy dissipation, good reversibility, and high symmetry in programming.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Situating the Vector Density Approach Among Contemporary Continuum Theories of Dislocation Dynamics

Abstract For the past century, dislocations have been understood to be the carriers of plastic deformation in crystalline solids. However, their collective behavior is still poorly understood. Progress in understanding the collective behavior of dislocations has primarily come in one of two modes: the simulation of systems of interacting discrete dislocations and the treatment of density measures of varying complexity that are considered as continuum fields. A summary of contemporary models of continuum dislocation dynamics is presented. These include, in order of complexity, the two-dimensional statistical theory of dislocations, the field dislocation mechanics treating the total Kröner–Nye tensor, vector density approaches that treat geometrically necessary dislocations on each slip system of a crystal, and high-order theories that examine the effect of dislocation curvature and distribution over orientation. Each of theories contain common themes, including statistical closure of the kinetic dislocation transport equations and treatment of dislocation reactions such as junction formation. An emphasis is placed on how these common themes rely on closure relations obtained by analysis of discrete dislocation dynamics experiments. The outlook of these various continuum theories of dislocation motion is then discussed.

Engineering↗
DOE OSTI2023

Modeling shock-induced void collapse in single-crystal Ta systems at the mesoscales

Understanding the role of microstructural heterogeneities on the shock wave propagation and defect evolution behavior is essential to predicting the dynamic response of metals. Heterogeneities, such as voids, provide challenges to understanding the wave propagation behavior as the shock-void interaction can collapse the void and result in large plastic strains and significant heating (hot spot formation) in the metal. Accurate modeling of this phenomenon requires predicting the void collapse mechanisms and the related heat generation and dissipation mechanisms (hotspot formation) that determine the microstructure evolution. While molecular dynamics (MD) simulations can model the void collapse behavior, the time/length scale capabilities pose a challenge to connect with continuum models or the experimental scales. Here, this study presents the capability of the newly developed quasi-coarse-grained dynamics (QCGD) method that extends the MD simulations to larger system sizes and longer times to model this phenomenon. This study uses QCGD simulations to investigate the mechanisms of shock wave interactions with pre-existing voids in single-crystal Ta microstructures for variations in shock pressures, void size, and loading orientations. For a given orientation, the plasticity contributions and rates of void collapse are observed to vary with shock pressures and void size. The larger void sizes and higher pressures result in increased temperatures (hot spots) and faster void collapse rates and unravel the variations in the plasticity contributions. In addition, QCGD simulations investigate the post-collapse microstructure evolution as a release wave travels through the hot spot region. The simulations reveal that the reduction in temperatures due to heat dissipation initiates the dynamic recrystallization behavior in the hotspot regions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Magic Angles and Fractional Chern Insulators in Twisted Homobilayer Transition Metal Dichalcogenides

We explain the appearance of magic angles and fractional Chern insulators in twisted K-valley homobilayer transition metal dichalcogenides by mapping their continuum model to a Landau level problem. Here, our approach relies on an adiabatic approximation for the quantum mechanics of valence band holes in a layer-pseudospin field that is valid for sufficiently small twist angles and on a lowest Landau level approximation that is valid for sufficiently large twist angles. It provides a simple qualitative explanation for the nearly ideal quantum geometry of the lowest moir´e miniband at particular twist angles, predicts that topological flat bands occur only when the valley-dependent moir´e potential is sufficiently strong compared to the interlayer tunneling amplitude, and provides a convenient starting point for the study of interactions.

36 MATERIALS SCIENCE↗