Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dissipative particle 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.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Kinetic Simulations of Radiative Magnetic Reconnection in the Coronae of Accreting Black Holes

We perform 2D and 3D particle-in-cell simulations of reconnection in magnetically dominated e ± plasmas subject to strong Compton cooling. Magnetic reconnection under such conditions can operate in accretion disk coronae around black holes, which produce hard X-rays through Comptonization. Our simulations show that most of the plasma in the reconnection layer is kept cold by Compton losses and locked in magnetically dominated plasmoids with a small thermal pressure. Compton drag clears cavities inside plasmoids and also affects their bulk motions. These effects, however, weakly change the reconnection rate and the plasmoid size distribution from those in nonradiative reconnection. This demonstrates that the reconnection dynamics is governed by similar magnetic stresses in both cases and weakly affected by thermal pressure. We examine the energy distribution of particles energized by radiative reconnection and observe two distinct components: (1) A mildly relativistic peak, which results from bulk motions of cooled plasmoids. This component receives most of the dissipated reconnection power and dominates the output X-ray emission. The peak has a quasi-Maxwellian shape with an effective temperature of ~100 keV. Thus, it mimics thermal Comptonization used previously to fit hard-state spectra of accreting black holes. (2) A high-energy tail, which receives ~20% of the dissipated reconnection power. It is populated by particles accelerated impulsively at X-points or "picked up" by fast outflows from X-points. The high-energy particles immediately cool, and their inverse Compton emission explains the MeV spectral tail detected in the hard state of Cyg X-1. Our first-principle simulations support magnetic reconnection as a mechanism powering hard X-ray emission from magnetically dominated regions of accreting black holes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Subphotospheric Emission from Short Gamma-Ray Bursts: Protons Mold the Multimessenger Signals

Abstract The origin of the observed Band-like photon spectrum in short gamma-ray bursts (sGRBs) is a long-standing mystery. We carry out the first general relativistic magnetohydrodynamic simulation of an sGRB jet with initial magnetizationσ 0 = 150 in dynamical ejecta from a binary merger. From this simulation, we identify regions along the jet of efficient energy dissipation due to magnetic reconnection and collisionless subshocks. Taking into account electron and proton acceleration processes, we solve for the first time the coupled transport equations for photons, electrons, protons, neutrinos, and intermediate particle species up to close to the photosphere (i.e., up to 1 × 10 12 cm), accounting for all relevant radiative and cooling processes. We find that the subphotospheric multimessenger signals carry strong signatures of the hadronic interactions and their resulting particle cascades. Importantly, the spectral energy distribution of photons is significantly distorted with respect to the Wien one, commonly assumed below the photosphere. Our findings suggest that the bulk of the nonthermal photon spectrum observed in sGRBs can stem from hadronic processes occurring below the photosphere and previously neglected, with an accompanying energy flux of neutrinos peaking in the GeV energy range.

Astronomy & Astrophysics↗
DOE OSTI2020

Self-sustaining sound in collisionless, high- β plasma

Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $δn/n_0≳2/β$ (where $β\gg1$ is the ratio of thermal to magnetic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly to reduce the pressure anisotropy by pitch-angle scattering and trapping particles, respectively, thereby impeding the maintenance of Landau resonances that enable such waves’ otherwise potent collisionless damping. The result is wave dynamics that evince a weakly collisional plasma: the ion distribution function is near-Maxwellian, the field-parallel flow of heat resembles its Braginskii form (except in regions where large-amplitude magnetic mirrors strongly suppress particle transport), and the relations between various thermodynamic quantities are more ‘fluid-like’ than kinetic. A nonlinear fluctuation–dissipation relation for self-sustaining IAWs is obtained by solving a plasma-kinetic Langevin problem, which demonstrates suppressed damping, enhanced fluctuation levels and weakly collisional thermodynamics when IAWs with $δn/n_0≳2/β$ are stochastically driven. We investigate how our results depend upon the scale separation between the wavelength of the IAW and the Larmor radius of the ions, and discuss briefly their implications for our understanding of turbulence and transport in the solar wind and the intracluster medium of galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Active Colloids as Models, Materials, and Machines

Active colloids use energy input at the particle level to propel persistent motion and direct dynamic assemblies. We consider three types of colloids animated by chemical reactions, time-varying magnetic fields, and electric currents. For each type, we review the basic propulsion mechanisms at the particle level and discuss their consequences for collective behaviors in particle ensembles. These microscopic systems provide useful experimental models of nonequilibrium many-body physics in which dissipative currents break time-reversal symmetry. Freed from the constraints of thermodynamic equilibrium, active colloids assemble to form materials that move, reconfigure, heal, and adapt. Colloidal machines based on engineered particles and their assemblies provide a basis for mobile robots with increasing levels of autonomy. This review provides a conceptual framework for understanding and applying active colloids to create material systems that mimic the functions of living matter. We highlight opportunities for chemical engineers to contribute to this growing field.

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

Kramers' degeneracy for open systems in thermal equilibrium

Kramers' degeneracy theorem underpins many interesting effects in quantum systems with time-reversal symmetry. We show that the generator of dynamics for Markovian open fermionic systems can exhibit an analogous degeneracy, protected by a combination of time-reversal symmetry and the microreversibility (detailed balance) property of systems at thermal equilibrium—the degeneracy is lifted if either condition is not met. We provide simple examples of this phenomenon and show that the degeneracy is reflected in the single-particle Green's functions. Furthermore, we show that certain experimental signatures of topological edge modes in open many-body systems can be protected by microreversibility in the same way. Our results highlight the importance of detailed balance in characterizing open topological matter.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Kinetic Simulations of the Kruskal–Schwarzschild Instability in Accelerating Striped Outflows: Dynamics and Energy Dissipation

Astrophysical relativistic outflows are launched as Poynting-flux dominated, yet the mechanism governing efficient magnetic dissipation, which powers the observed emission, is still poorly understood. We study magnetic energy dissipation in relativistic "striped" jets, which host current sheets separating magnetically dominated regions with opposite field polarity. The effective gravity force g in the rest frame of accelerating jets drives the Kruskal–Schwarzschild instability (KSI), a magnetic analog of the Rayleigh–Taylor instability. By means of 2D and 3D particle-in-cell simulations, we study the linear and nonlinear evolution of the KSI. The linear stage is well described by linear stability analysis. The nonlinear stages of the KSI generate thin (skin-depth-thick) current layers, with length comparable to the dominant KSI wavelength. There, the relativistic drift-kink mode and the tearing mode drive efficient magnetic dissipation. The dissipation rate can be cast as an increase in the effective width Δ eff of the dissipative region, which follows dΔ eff /dt ≃ 0.05 $\sqrt{Δ_{eff} g}$. Our results have important implications for the location of the dissipation region in gamma-ray burst and active galactic nuclei jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dynamics of nonequilibrium magnons in gapped Heisenberg antiferromagnets

Nonequilibrium dynamics in spin systems is a topic currently under intense investigation as it provides fundamental insights into thermalization, universality, and exotic transport phenomena. While most of the studies have been focused on ideal closed quantum many-body systems such as ultracold atomic quantum gases and one-dimensional spin chains, driven-dissipative Bose gases in steady states away from equilibrium in classical systems also lead to intriguing nonequilibrium physics. Here, in this work, we theoretically investigate out-of-equilibrium dynamics of magnons in a gapped Heisenberg quantum antiferromagnet based on Boltzmann transport theory. We show that, by treating scattering terms beyond the relaxation-time approximation in the Boltzmann transport equation, energy and particle number conservation mandate that nonequilibrium magnons cannot relax to equilibrium, but decay to other nonequilibrium stationary states. The only decay channel for these stationary states back to equilibrium is through the nonconserving interactions (i.e., changing particle number and/or energy within the magnon system) such as boundary or magnon-phonon scattering. At low temperatures, these nonconserving interactions are much slower processes than intrinsic magnon-magnon interaction in a gapped spin system. Using magnon-phonon interaction as a quintessential type of nonconserving interaction, we then propose that nonequilibrium steady states of magnons can be maintained and tailored using periodic driving at frequencies faster than relaxation due to phonon interactions. These findings reveal a class of classical material systems that are suitable platforms to study nonequilibrium statistical physics and macroscopic phenomena such as classical Bose-Einstein condensation of quasiparticles and magnon supercurrents that are relevant for spintronic applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A particle finite element-based model for droplet spreading analysis

Here, a particle finite element method-based model is proposed to analyze droplet dynamics problems, particularly droplet spreading on solid substrates (wetting). The model uses an updated Lagrangian framework to formulate the governing equations of the liquid. The curvature of the liquid surface is tracked accurately using a deforming boundary mesh. In order to predict the spreading rate of the droplet on the solid substrate and track the corresponding contact angle evolution, dissipative forces at the contact line are included in the formulation in addition to the Navier-slip boundary conditions at the solid–liquid interface. The inclusion of these boundary conditions makes it possible to account for the induced Young’s stress at the contact line and for the viscous dissipation along the solid–liquid interfacial region. These are found to be essential to obtain a mesh-independent physical solution. The temporal evolution of the contact angle and the contact line velocity of the proposed model are compared with spreading droplets and micro-sessile droplet injection experiments and are shown to be in good agreement.

42 ENGINEERING↗
DOE OSTI2024

Collaborative Research: Vlasov–Maxwell Simulations to resolve electron heating and dissipation, in quasiperpendicular shocks (Final Technical Report)

The project is to study electron heating at quasi-perpendicular shocks. We have analyzed measurements from the Magnetosphere Multiscale (MMS) spacecraft and particle-in-cell simulations to address the question. One key discovery is that waves and are excited in the shock transition region, which further develop into complicated structures and thin current layers, which contribute to accelerate and heat electrons. We further investigated the excitation and evolution of these waves and structures, and their association with ion and electron dynamics. In addition, the analysis technique of Vlasov equation terms is further applied to study electron dynamics at the dayside magnetopause, to understand the effect of density, velocity and temperature variations on electron distributions and kinetics.

58 GEOSCIENCES↗
DOE OSTI2020

Connecting energetics to dynamics in particle growth by oriented attachment using real-time observations

The interplay between crystal and solvent structure, interparticle forces and ensemble particle response dynamics governs the process of crystallization by oriented attachment (OA), yet a quantitative understanding is lacking. Using as a model system ZnO, we combine in situ TEM observations of single particle and ensemble assembly dynamics with simulations of interparticle forces and responses to relate experimentally derived interparticle potentials to the underlying interactions. We show that OA is driven by forces and torques due to a combination of electrostatic ion-solvent correlations and dipolar interactions that act at separations well beyond 5 nm. Importantly, coalignment is achieved before particles reach separations at which strong attractions drive the final jump to contact. The observed barrier to attachment is negligible, while dissipative factors in the quasi-2D confinement of the TEM fluid cell lead to abnormal diffusivities with timescales for rotation much less than for translation, thus enabling OA to dominate.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Stochastic mean-field theory and applications to multinucleon transfer and kinetic energy dissipation processes in heavy-ion collisions

In this Review article, a brief description of the stochastic mean-field (SMF) theory for describing reaction dynamics in low-energy heavy-ion collisions at bombarding energies in the vicinity of the Coulomb barrier is presented. In these collisions, as a result of strong Pauli blocking, binary nucleon collisions do not have a significant effect on the dissipation and fluctuations. At low energies, the mean-field fluctuations, due to initial correlations, have a dominant effect on fluctuations of macroscopic variables. The SMF theory proposes the determination of an ensemble of single-particle density matrices by specifying random initial fluctuations according to a distribution law. Employing an ensemble of single-particle density matrices, not only the mean values but also the distribution functions of the one-body observables can be determined. If the di-nuclear structure is maintained in heavy-ion collisions, such as deep inelastic collisions and fast quasi-fission reactions, a much simpler description of the reaction mechanism can be derived in terms of several macroscopic variables such as mass and charge asymmetry, and relative linear and relative angular momentum. In this case, by geometric projection of the SMF equations, it is possible to derive the quantal Langevin equations for macroscopic variables. As an application of quantal transport description, an analysis of multinucleon transfers and kinetic energy dissipation and fluctuations is presented for selected quasi-fission reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Structure preservation using discrete gradients in the Vlasov-Poisson-Landau system

We present a novel structure-preserving framework for solving the Vlasov-Poisson-Landau system of equations using a particle in cell (PIC) discretization combined with discrete gradient time integrators. The Vlasov-Poisson-Landau system is an accurate model for studying hot plasma dynamics at a kinetic scale where small-angle Coulomb collisions dominate. Our scheme guarantees conservation of mass, momentum and energy as well as preservation of the monotonicity of entropy production in both the time-continuous and discrete systems. We employ the conservative integrator for both the Hamiltonian Vlasov-Poisson equations and the dissipative Landau equation using the PETSc library (www.mcs.anl.gov/petsc) to showcase structure-preserving properties.

Discrete gradients↗
DOE OSTI2021

First evidence of dominant influence of E × B drifts on plasma cooling in an advanced slot divertor for tokamak power exhaust

Addressing power exhaust in tokamaks is presently recognized as one of the major remaining open issues for the development of fusion reactors. At the forefront of this endeavor is the effort to develop an advanced divertor by maximizing dissipation of plasma power and momentum inside the divertor. Here, we demonstrate, for the first time, that the electromagnetic (E × B) drifts exert a key influence on plasma and gas dynamics in a new advanced slot divertor in the DIII-D tokamak, named the small angle slot (SAS). SAS leverages the effect of drifts to achieve a highly dissipative divertor with electron temperature T e ≲ 10 eV over a wide range of plasma densities, for ion B × ∇B away from the divertor, as used for advanced tokamak operation on DIII-D. Modeling with the SOLPS-ITER code shows that for this drift direction, the E × B flow carries particles toward the outer common flux region (CFR) via the private flux region (PFR), reinforcing neutral recycling and enhancing divertor dissipation. In contrast, for the opposite field direction, the E × B flow carries particles away from the outer CFR into the PFR, offsetting the anticipated SAS geometric effects. Furthermore, this finding is an important step in the understanding of the behavior of advanced divertors for achieving a power exhaust solution for fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Properties of Hot Nuclear Matter

A fully quantitative description of the equilibrium and dynamical properties of hot nuclear matter will be needed for the interpretation of the available and forthcoming astrophysical data, providing information on the post-merger phase of a neutron star coalescence. We discuss the results of a recently developed theoretical model, based on a phenomenological nuclear Hamiltonian including two- and three-nucleon potentials, to study the temperature dependence of average and single-particle properties of nuclear matter relevant to astrophysical applications. The potential of the proposed approach for describing dissipative processes leading to the appearance of bulk viscosity in neutron star matter is also outlined.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Microstructure effects on high velocity microparticle impacts of copper

Constitutive models can fail to predict high-rate deformation behavior due to their inability to account for microstructural effects. In part, this is because of a dearth of experimental benchmarking data in the high strain-rate, low pressure regime, since many high-rate experiments also probe a region of strong shockwaves, at which point microstructure effects no longer play a primary role. This work uses laser-induced particle impact testing to quantitatively study high velocity impacts of small, rigid alumina microspheres on flat copper substrates with varying amounts of initial cold work in the weak shock regime, but at very high strain rates up to ~10 7 s –1 . Through paired experiments and numerical simulations, this work shows that the initial microstructure condition can have significant influence on dynamical mechanical properties in this range. Specifically, prior work hardening of the copper substrate leads to increased rebounding of the microparticles (i.e., less plastic dissipation in the impact) as well as smaller craters. Each of these experimental measurables can be converted into a strength measure, i.e., the dynamic yield strength or dynamic hardness, respectively, neither of which is well predicted consistently by existing constitutive laws. The general trend of hardening can be captured by such models by incorporating an existing “pre-strain,” suggesting that future calibration of the materials parameters may yield a good fit over a broader range of conditions. Our results emphasize the importance of reporting the microstructural condition in dynamic studies, as well as the necessity of accounting for these factors when formulating and optimizing constitutive models.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dissipative ground state preparation in ab initio electronic structure theory

Dissipative engineering is a powerful tool for quantum state preparation, and has drawn significant attention in quantum algorithms and quantum many-body physics in recent years. In this work, we introduce a novel approach using the Lindblad dynamics to efficiently prepare the ground state for general ab initio electronic structure problems on quantum computers, without variational parameters. These problems often involve Hamiltonians that lack geometric locality or sparsity structures, which we address by proposing two generic types of jump operators for the Lindblad dynamics. Type-I jump operators break the particle number symmetry and should be simulated in the Fock space. Type-II jump operators preserves the particle number symmetry and can be simulated more efficiently in the full configuration interaction space. For both types of jump operators, we prove that in a simplified Hartree-Fock framework, the spectral gap of our Lindbladian is lower bounded by a universal constant. For physical observables such as energy and reduced density matrices, the convergence rate of our Lindblad dynamics with Type-I jump operators remains universal, while the convergence rate with Type-II jump operators only depends on coarse grained information such as the number of orbitals and the number of electrons. To validate our approach, we employ a Monte Carlo trajectory-based algorithm for simulating the Lindblad dynamics for full ab initio Hamiltonians, demonstrating its effectiveness on molecular systems amenable to exact wavefunction treatment.

Quantum chemistry↗
DOE OSTI2020

Ion irradiation and modification: The role of coupled electronic and nuclear energy dissipation and subsequent nonequilibrium processes in materials

Understanding material responses to energy deposition from energetic charged particles is important for defect engineering, ion-beam processing, ion-beam analysis and modification, geologic aging, space exploration, and nuclear applications. As an incident ion penetrates a solid, its energy is transferred to electrons and to atomic nuclei of the solid. Much of this electronic energy deposition is subsequently transferred to the atomic structure via electron–phonon (e–ph) coupling, leading to local inelastic thermal spikes in which energy dissipation is influenced by the local environment. In addition, intense ionization can lead to high densities of localized electronic excitations in wide-bandgap materials and ceramics that can affect defect dynamics and atomic mobility. Specifically, energy exchange between electrons and atomic nuclei, along with localized electronic excitations, can lead to substantial competitive (ionization-induced annealing), additive (both electronic and nuclear energy depositions contributing to damage production), and synergistic (more damage than the sums of separate damage processes) effects. Although nonmonotonic effects of the e–ph coupling strength and athermal processes are demonstrated for pre-existing defects and residual damage during ion–solid interactions, there is limited understanding of when such electronic effects must be considered in atomic-scale models of damage production and evolution in a broad variety of materials. Complex ceramics and chemically disordered solid solution alloys with different constituent elements allow a systematic evaluation of defect dynamics and irradiation performance with increasing complexity. Current knowledge regarding tuning of bonding characteristics and chemical disorder to control atomic-level dynamics is reviewed. Although a lack of fundamental understanding obstructs the advancement of reliable predictions for ion beam material modification, it highlights challenges and opens research opportunities. Insights into the complex electronic and atomic correlations with extreme energy deposition will strengthen our ability to design materials and predict ion-irradiation-induced damage in a radiation environment, and they may pave the way to better control fundamental processes and design new material functionalities for advanced technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Collisional-energy-cascade model for nonthermal velocity distributions of neutral atoms in plasmas

Nonthermal velocity distributions with much greater tails than a Maxwellian have been observed for radical atoms in plasmas for a long time. Historically, such velocity distributions have been modeled by a two-temperature Maxwell distribution. In this paper, I propose a model based on collisional energy cascade, which has been studied in the field of granular materials. In the collisional energy cascade, a particle ensemble undergoes energy input at the high-energy region, entropy production by elastic collisions among particles, and energy dissipation. For radical atoms, energy input may be caused by the Franck-Condon energy of molecular dissociation or charge-exchange collision with hot ions, and the input energy is eventually dissipated by collisions with the walls. Here, I show that the steady-state velocity distribution in the collisional energy cascade is approximated by the generalized Mittag-Leffler distribution, which is a one-parameter extension of the Maxwell distribution. This parameter indicates the degree of the nonthermality and is related to the relative importance of energy dissipation over entropy production. This model is compared with a direct molecular dynamics simulation for simplified gaseous systems with energy input, as well as some experimentally observed velocity distributions of light radicals in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗