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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Dispersion and the speed-limited particle-in-cell algorithm

This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. Here, these results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Hamiltonian structure of 2D fluid dynamics with broken parity

Isotropic fluids in two spatial dimensions can break parity symmetry and sustain transverse stresses which do not lead to dissipation. Corresponding transport coefficients include odd viscosity, odd torque, and odd pressure. We consider an isotropic Galilean invariant fluid dynamics in the adiabatic regime with momentum and particle density conservation. We find conditions on transport coefficients that correspond to dissipationless and separately to Hamiltonian fluid dynamics. The restriction on the transport coefficients will help identify what kind of hydrodynamics can be obtained by coarse-graining a microscopic Hamiltonian system. Interestingly, not all parity-breaking transport coefficients lead to energy conservation and, generally, the fluid dynamics is energy conserving but not Hamiltonian. We show how this dynamics can be realized by imposing a nonholonomic constraint on the Hamiltonian system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

Magnetic blowout during collisionless reconnection

Analytic calculations and simulations of the steady-state structure of the dissipation region during collisionless reconnection are described. Ions are treated as particles and electrons as a dissipative fluid. A limiting reconnection rate results as the ions 'blowout' the magnetic field in the dissipation area. The implications for comprehending substorm dynamics are discussed.

Drake, J. F.↗
DOE OSTI2021

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

In this work, we employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant–matrix pair correlation function and the penetrant–matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

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

An Assessment of the Drag Models in the Case of a Shock Interacting With a Fixed Bed of Point Particles

Abstract In this work, three-dimensional Euler–Lagrange (EL) point-particle simulations of a shock wave interacting with a fixed bed of particles are carried out. The results from the particle-resolved (PR) simulations are used to assess the performance of the point-particle drag models during short time scales. We demonstrate that in a one-way coupled regime, the point-particle simulations recover the dominant gas dynamic features of the flow and are in a good agreement with the exact Riemann solution of a shock traveling through a sudden area contraction. Although the PR simulations are inviscid, we show that a dissipative drag is necessary to predict the mean behavior of the gas. As a model for the inviscid shock-induced (SI) drag two different models are presented in lieu of the quasi-steady drag. Finally, two-way coupled simulations are performed at four different particle volume fractions {0.10, 0.15, 0.20, 0.25} and three different incident shock Mach numbers {1.22, 1.66, 3.0} and compared against the data from PR inviscid simulations. At a lower Mach number (1.22), averaged flow quantities from the two-way coupled simulations agree well with the PR simulations. As the Mach number increases, we observe that the discrepancies between the point-particle and the PR simulations grow. A sensitivity analysis of the drag models involved reveals a strong influence of the inviscid-unsteady force on the gas quantities especially in the case of a strong shock interacting with a dense bed of particles. The use of Mach correlation beyond the subcritical regime coupled with the model for volume fraction correction is identified as a probable cause for the additional drag.

Engineering↗
DOE OSTI2025

Dynamical Complexity of Non-Gaussian Many-Body Systems with Dissipation

We characterize the dynamical state of many-body bosonic and fermionic many-body models with intersite Gaussian couplings, on-site non-Gaussian interactions, and local dissipation comprising incoherent particle loss, particle gain, and dephasing. We first establish that, for fermionic systems, if the dephasing noise is larger than the non-Gaussian interactions, irrespective of the Gaussian coupling strength, the system state is a convex combination of Gaussian states at all times. Furthermore, for bosonic systems, we show that if the particle loss and particle gain rates are larger than the Gaussian intersite couplings, the system remains in a separable state at all times. Building on this characterization, we establish that at noise rates above a threshold, there exists a classical algorithm that can efficiently sample from the system state of both the fermionic and bosonic models. Finally, we show that, unlike fermionic systems, bosonic systems can evolve into states that are not convex Gaussian even when the dissipation is much higher than the on-site non-Gaussianity. Similarly, unlike bosonic systems, fermionic systems can generate entanglement even with noise rates much larger than the intersite couplings.

Computational complexity↗
NASA NTRS2004

Complexity Induced Anisotropic Bimodal Intermittent Turbulence in Space Plasmas

The "physics of complexity" in space plasmas is the central theme of this exposition. It is demonstrated that the sporadic and localized interactions of magnetic coherent structures arising from the plasma resonances can be the source for the coexistence of nonpropagating spatiotemporal fluctuations and propagating modes. Non-Gaussian probability distribution functions of the intermittent fluctuations from direct numerical simulations are obtained and discussed. Power spectra and local intermittency measures using the wavelet analyses are presented to display the spottiness of the small-scale turbulent fluctuations and the non-uniformity of coarse-grained dissipation that can lead to magnetic topological reconfigurations. The technique of the dynamic renormalization group is applied to the study of the scaling properties of such type of multiscale fluctuations. Charged particle interactions with both the propagating and nonpropagating portions of the intermittent turbulence are also described.

Chang, Tom↗
DOE OSTI2022

Extracting forces from noisy dynamics in dusty plasmas

Extracting environmental forces from noisy data is a common yet challenging task in complex physical systems. Machine learning (ML) represents a robust approach to this problem, yet is mostly tested on simulated data with known parameters. Here we use supervised ML to extract the electrostatic, dissipative, and stochastic forces acting on micron-sized charged particles levitated in an argon plasma (dusty plasma). By tracking the sub-pixel motion of particles in subsequent images, we successfully estimated these forces from their random motion. The experiments contained important sources of non-Gaussian noise, such as drift and pixel-locking, representing a data mismatch from methods used to analyze simulated data with purely Gaussian noise. Our model was trained on simulated particle trajectories that included all of these artifacts, and used more than 100 dynamical and statistical features, resulting in a prediction with 50\% better accuracy than conventional methods. Lastly, in systems with two interacting particles, the model provided non-contact measurements of the particle charge and Debye length in the plasma environment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Formation of giant molecular clouds in global spiral structures: The role of orbital dynamics and cloud-cloud collisions

The different roles played by orbital dynamics and dissipative cloud-cloud collisions in the formation of giant molecular clouds (GMCs) in a global spiral structure are investigated. The interstellar medium (ISM) is simulated by a system of particles, representing clouds, which orbit in a spiral-perturbed, galactic gravitational field. The overall magnitude and width of the global cloud density distribution in spiral arms is very similar in the collisional and collisionless simulations. The results suggest that the assumed number density and size distribution of clouds and the details of individual cloud-cloud collisions have relatively little effect on these features. Dissipative cloud-cloud collisions play an important steadying role for the cloud system's global spiral structure. Dissipative cloud-cloud collisions also damp the relative velocity dispersion of clouds in massive associations and thereby aid in the effective assembling of GMC-like complexes.

Roberts, W. W., Jr.↗
DOE OSTI2023

Fast particles in drift wave turbulence

This study aims to incorporate the effects of fast particles into our present fluid model for tokamak transport. The parameter ε f = ω / ω f, where ω is the mode frequency and ω f is the typical frequency of the fast particles, which enters as a factor in front of the fast particle response. Thus, for trapped fast particles, where ω f = ω pres the precession frequency of the fast particles, this parameter is of order 10 – 2 for drift waves, and thus, the fast particle response can be neglected. However, ε f will be of order 1 for fast particle modes such as in the fishbone instability. An important turbulence property, affecting both these limits, is resonance broadening. Effects of resonance broadening have recently been considered for fast particle instabilities, often coupled directly to the linear growth rate, while we here consider the original Dupree formulation where the turbulence directly drives a nonlinear frequency shift. Resonance broadening has a general tendency to counteract dissipative wave particle resonances. This has been observed for fast particle instabilities. Here, there is a resonant external source for the fast particles, so the instability survives if this source is dominant over the resonance broadening. For drift waves, however, external sources are not resonant since ε f << 1. Furthermore, the resonance broadening is able to remove the dissipative wave particle resonance completely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Acoustic waves in the solar atmosphere. VII - Non-grey, non-LTE H(-) models

The propagation and shock formation of radiatively damped acoustic waves in the solar chromosphere are studied under the assumption that H(-) is the only absorber; the opacity is non-grey. Deviations from local thermodynamic equilibrium (LTE) are permitted. The results of numerical simulations show the depth dependence of the heating by the acoustic waves to be insensitive to the mean state of the atmosphere. After the waves have developed into shocks, their energy flux decays exponentially with a constant damping length of about 1.4 times the pressure scale height, independent of initial flux and wave period. Departures from LTE have a strong influence on the mean temperature structure in dynamical chromosphere models; this is even more pronounced in models with reduced particle density - simulating conditions in magnetic flux tubes - which show significantly increased temperatures in response to mechanical heating. When the energy dissipation of the waves is sufficiently large to dissociate most of the H(-) ions, a strong temperature rise is found that is reminiscent of the temperature structure in the transition zone between chromosphere and corona; the energy flux remaining in the waves then drives mass motions.

Schmitz, F.↗
DOE OSTI2022

Designing Stress-Adaptive Dense Suspensions Using Dynamic Covalent Chemistry

The non-Newtonian behaviors of dense suspensions are central to their use in technological and industrial applications and arise from a network of particle–particle contacts that dynamically adapt to imposed shear. Reported herein are studies aimed at exploring how dynamic covalent chemistry between particles and the polymeric solvent can be used to tailor such stress-adaptive contact networks, leading to their unusual rheological behaviors. Specifically, a room temperature dynamic thia-Michael bond is employed to rationally tune the equilibrium constant (K eq ) of the polymeric solvent to the particle interface. It is demonstrated that low Keq leads to shear thinning, while high K eq produces antithixotropy, a rare phenomenon where the viscosity increases with shearing time. It is proposed that an increase in K eq increases the polymer graft density at the particle surface and that antithixotropy primarily arises from partial debonding of the polymeric graft/solvent from the particle surface and the formation of polymer bridges between particles. Thus, the implementation of dynamic covalent chemistry provides a new molecular handle with which to tailor the macroscopic rheology of suspensions by introducing programmable time dependence. These studies open the door to energy-absorbing materials that not only sense mechanical inputs and adjust their dissipation as a function of time or shear rate but also can switch between these two modalities on demand.

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

Tracking Optical Phonon Dynamics in InP Nanocrystals via Transient Absorption and Femtosecond Stimulated Raman Spectroscopy

Semiconductor nanocrystals (NCs) offer actualized and prospective utility in optoelectronic technologies, yet key aspects of their thermal and vibrational behaviors remain unresolved. Compared with bulk crystals, heat dissipation involving NCs can differ substantially owing to pervasiveness of interfacial scattering, phonon confinement, and relaxed phonon momentum selection rules, as well as influences of ligands. Here, in addition to transient absorption, we use femtosecond stimulated Raman spectroscopy to track nonequilibrium optical-phonon dynamics in InP NCs capped with myristic acid ligands and compare the same particles coated with a ZnS shell. We identify distinct phonon decay pathways in ligand-capped versus ZnS-shelled NCs. Increasing the excitation density and introducing the ZnS shell measurably modify optical-phonon lifetimes, whereas phonon formation times remain largely unchanged. In addition, we observe phonon mode softening in ligand-capped nanocrystals, consistent with lattice expansion. Together, these results demonstrate that core–shell structures can strongly govern nanocrystal thermal dissipation pathways and should be considered a key design parameter for optoelectronic operation.

dissipation↗
DOE OSTI2024

Theories of Relativistic Dissipative Fluid Dynamics

Relativistic dissipative fluid dynamics finds widespread applications in high-energy nuclear physics and astrophysics. However, formulating a causal and stable theory of relativistic dissipative fluid dynamics is far from trivial; efforts to accomplish this reach back more than 50 years. In this review, we give an overview of the field and attempt a comparative assessment of (at least most of) the theories for relativistic dissipative fluid dynamics proposed until today and used in applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Collaborative Research: Vlasov-Maxwell Simulations to Resolve Electron Heating and Dissipation, in Quasi-Perpendicular Shocks (Final Report)

Collisionless shocks are a grand challenge problem in plasma physics and have been the subject of study for more than six decades. Shocks are ubiquitous phenomena in the universe and are responsible for transforming high energy and momentum flows into thermal energy and energetic particles. Understanding how shocks operate is of primary importance to understand the sun-earth coupling, protecting manned missions and spacecrafts from high energy particles, inertial confinement fusion, and radiation observed from astrophysical plasmas, such as supernova remnants and astrophysical jets. A major unanswered question on this frontier is, how does a collisionless plasma transform flow energy and momentum into electron thermal energy? Many potential mechanisms have been proposed to perform the conversion between flow and thermal energy in a collisionless plasma, but the answer has been elusive. To solve the mystery, non-relativistic perpendicular and quasi-perpendicular shocks are studied via a coordinated approach employing cutting-edge numerical simulations with the continuum kinetic code Gkeyll and spacecraft observations supplied by NASA's Magnetospheric Multiscale (MMS) mission. In both sets of data, we employ novel techniques to identify how and where energy is exchanged between electromagnetic fields and plasma particles. By combining cutting-edge spacecraft data and numerical simulations with these diagnostics, we have been able to identify some of the primary mechanisms responsible electron heating and dissipation in heliospheric shocks. In addition to electron dynamics, we have complete access to proton simulation and in situ data, enabling us to characterize proton heating and acceleration in quasi-perpendicular shocks. Many astrophysical quasi-perpendicular shocks have similar parameters to those in this study, including supernova remnants and shocks in galaxy clusters. Additionally, the simulation code, Gkeyll, and associated analysis scripts are openly available for community use. This award has supported a collaborative effort between Princeton University, the University of Arizona - Tucson, and the University of Maryland - College Park. The award funded three early career scientists, including the funding and training of one female postdoctoral researcher, who was involved in all aspects of the MMS data analysis and some portions of the simulation analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

The Dissipation Mechanism of Magnetic Reconnection

Magnetic reconnection is arguably the most efficient transport and energy conversion mechanism in almost ideal plasmas. Reconnection controls the overall dynamics in space and astrophysics plasmas, as well as in many laboratory plasma systems. Reconnection operates by means of a localized diffusion region, where deviations from the plasma idealness condition generate electric fields and permit plasma transport even far away from the diffusion region itself. Recent advances in analytic theory and computer modeling have begun to shed light on the internal dynamics of the diffusion region. In particular, we begin to understand the delicate nature of the force balance in the inner diffusion region, where particles can become unmagnetized and where electric field forces are important. This presentation will provide a brief introduction of the reconnection process and its applications. This introduction will be followed by a detailed analysis of the current understanding of dissipation region physics, and by an outlook toward future research.

Hesse, Michael↗
NASA NTRS2008

The Dissipation Mechanism of Magnetic Reconnection

Magnetic reconnection is arguably the most efficient transport and energy conversion mechanism in almost ideal plasmas. Reconnection controls the overall dynamics in space and astrophysics plasmas, as well as in many laboratory plasma systems. Reconnection operates by means of a localized diffusion region, where deviations from the plasma idealness condition generate electric fields and permit plasma transport even far away from the diffusion region itself. Recent advances in analytic theory and computer modeling have begun to shed light on the internal dynamics of the diffusion region. In particular, we begin to understand the delicate nature of the force balance in the inner diffusion region, where particles can become unmagnetized and where electric field forces are important. This presentation will provide a brief introduction of the reconnection process and its applications. This introduction will be followed by a detailed analysis of the current understanding of dissipation region physics, and by an outlook toward future research.

Hesse, Michael↗