SEARCH · Engineering Papers
Results for “nonequilibrium thermodynamics”
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.
Flux-split algorithms for flows with non-equilibrium chemistry and vibrational relaxation
The present consideration of numerical computation methods for gas flows with nonequilibrium chemistry thermodynamics gives attention to an equilibrium model, a general nonequilibrium model, and a simplified model based on vibrational relaxation. Flux-splitting procedures are developed for the fully-coupled inviscid equations encompassing fluid dynamics and both chemical and internal energy-relaxation processes. A fully coupled and implicit large-block structure is presented which embodies novel forms of flux-vector split and flux-difference split algorithms valid for nonequilibrium flow; illustrative high-temperature shock tube and nozzle flow examples are given.
Nonequilibrium effects for hypersonic transitional flows
Presented are the results of numerical simulations of hypersonic flow about blunt cones and hemispherical nose configurations for reentry velocities of 7.5 and 10 km/s. Cone half angles 0, 5, and 10 deg are considered at zero angle of incidence; however, the focus is for the 5 deg cone. The body size and altitude ranges considered (70 to 110 km) are such that the flow is in the transitional regime. Translational, thermodynamic, and chemical nonequilibrium effects are considered in the numerical simulation by utilizing the direct simulation Monte Carlo (DSMC) method of Bird. The DSMC results are compared with those obtained with viscous shock-layer and Navier-Stokes methods. Comparisons between the DSMC and continuum calculations show the altitude range where differences in flowfield structure and surface quantities become significant. The current calculations show that the binary scaling similitude provides a means of correlating the blunt body surface quantities in the hypersonic, transitional regime. Furthermore, for the higher velocity entry conditions, the results highlight some of the concerns in the application of multitemperature continuum formulations, particularly the use of some proposed functional relations for the chemical rate constants under thermodynamic nonequilibrium conditions.
NONEQUILIBRIUM GASDYNAMICS RELATED TO PROPULSION SYSTEMS.
Nonequilibrium gas dynamics related to propulsion systems
Foundations of magnetohydrodynamics
In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.
Required toroidal confinement for fusion and omnigeneity
Deuterium–tritium (DT) burning requires a long energy confinement times compared to collision times, so the particle distribution functions must approximate local Maxwellians. Non-equilibrium thermodynamics is applicable, which gives relations among transport, entropy production, the collision frequency, and the deviation from a Maxwellian. The distribution functions are given by the Fokker–Planck equation, which is an advection–diffusion equation. A large hyperbolic operator, the Vlasov operator with the particle trajectories as its characteristics, equals a small diffusive operator, the collision operator. The collisionless particle trajectories would be chaotic in stellarators without careful optimization. This would lead to rapid entropy production and transport—far beyond what is consistent with a self-sustaining DT burn. Omnigeneity is the weakest general condition that is consistent with a sufficiently small entropy production associated with the thermal particle trajectories. Omnigeneity requires that the contours of constant magnetic field strength be unbounded in at least one of the two angular coordinates in magnetic surfaces and that there be a symmetry in the field-strength wells along the field lines. Even in omnigenous plasmas, fluctuations due to microturbulence can produce chaotic particle trajectories and the gyro-Bohm transport is seen in many stellarator and tokamak experiments. The higher the plasma temperature above 10 keV, the smaller the transport must be compared to gyro-Bohm for a self-sustaining DT burn. The hot alphas of DT fusion heat the electrons. When the ion–electron equilibration time is long compared to the ion energy confinement time, a self-sustaining DT burn is not possible, which sets a limit on the electron temperature.
On principles of emergent organization
After more than a century of concerted effort, physics still lacks basic principles of spontaneous organization. To appreciate why, we first state the problem, outline historical approaches, and survey the present state of the physics of self-organization. This frames the particular challenges arising from mathematical intractability and the resulting need for computational approaches, as well as those arising from a chronic failure to define structure. Then an overview of two modern mathematical formulations of organization—intrinsic computation and evolution operators—lays out a way to overcome these challenges. Together, the vantage point they afford shows how to account for the emergence of structured states via a statistical mechanics of systems arbitrarily far from equilibrium. The result is a constructive path forward to principles of organization that builds on mathematical identification of structure.
Investigation of flow regimes in arc plasma–gas interactions using a two-temperature arc in crossflow model
The perpendicular impingement of a gas stream on an electric arc, a configuration known as the arc in crossflow, is of primary relevance in the study of plasma-gas interactions, as well as in industrial applications such as circuit breakers and wire-arc spraying. The flow dynamics in the arc in crossflow are the result of coupled fluid-thermal-electromagnetic phenomena accompanied by large property gradients, which can produce significant deviations from Local Thermodynamic Equilibrium (LTE) among electrons and gas species. These characteristics can lead to the establishment of distinct flow regimes depending on the relative values of the controlling parameters of the system, such as inflow velocity, arc current, and inter-electrode spacing. A two-temperature non-LTE (NLTE) model is used to investigate the arc dynamics and the establishment of flow regimes in the arc in crossflow. The plasma flow model is implemented within a nonlinear Variational Multiscale (VMS) numerical discretization approach that is less dissipative, and hence better suited to capture unstable behavior, than traditional VMS methods commonly used in computational fluid dynamics simulations. The Reynolds and the Enthalpy dimensionless numbers, characterizing the relative flow strength and arc strength, respectively, are chosen as the controlling parameters of the system. Simulation results reveal the onset of dynamic behavior and the establishment of steady, periodic, quasi-periodic, and chaotic or potentially turbulent, regimes, as identified by distinct spatiotemporal fluctuations. The computational findings reveal the role of increasing the relative arc strength on enhancing flow stability by delaying the growth of fluctuating and unstable flow behavior.
Stratification of polymer mixtures in drying droplets: Hydrodynamics and diffusion
In this work, we study the evaporation-induced stratification of a mixture of short and long polymer chains in a drying droplet using molecular simulations. We systematically investigate the effects of hydrodynamic interactions (HI) on this process by comparing hybrid simulations accounting for HI between polymers through the multiparticle collision dynamics technique with free-draining Langevin dynamics simulations neglecting the same. We find that the dried supraparticle morphologies are homogeneous when HI are included but are stratified in core–shell structures (with the short polymers forming the shell) when HI are neglected. The simulation methodology unambiguously attributes this difference to the treatment of the solvent in the two models. We rationalize the presence (or absence) of stratification by measuring phenomenological multicomponent diffusion coefficients for the polymer mixtures. The diffusion coefficients show the importance of not only solvent backflow but also HI between polymers in controlling the dried supraparticle morphology.
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.
Coupling solidification kinetics with phase-behavior computations in hydrodynamic simulations of high-pressure, dynamic-compression processes
In this study, we report a numerical scheme to integrate models for the kinetics of solidification processes together with phase-behavior computations in the context of continuum-scale hydrodynamic simulations. The objective of the phase-behavior computations is to determine the pressure and temperature, given the following three sets of inputs: (1) an appropriate equation of state to describe our system, (2) the phase fraction(s) produced by the kinetic models, (3) and the volume and internal energy obtained by solving the conservation equations that govern the hydrodynamic behavior. The kinetics are assumed to be governed by the Kolmogorov–Johnson–Mehl–Avrami equation, and the nucleation and growth rates that enter into that equation are functions of the pressure and temperature produced by the phase-behavior computations. Our formulation allows for the fluid and solid phases to be at different temperatures (thermal nonequilibrium) and pressures (arising from surface-tension-induced Laplace contributions). The formulation is presented in a fairly general setting that is independent of any particular material, although we demonstrate it in some examples related to high-energy-density science applications where materials are rapidly compressed to pressures exceeding several gigapascals in less than a microsecond. We conclude with a critical evaluation of our approach and provide suggestions for future work to improve the predictive capabilities and generality of the models.
Analytic and numerical vibronic spectra from quasi-classical trajectory ensembles
The truncated Wigner approximation to quantum dynamics in phase space is explored in the context of computing vibronic line shapes for monomer linear optical spectra. We consider multiple model potential forms including a shifted harmonic oscillator with both equal and unequal frequencies on the ground and excited state potentials as well as a shifted Morse potential model. For the equal-frequency shifted harmonic oscillator model, we derive an analytic expression for the exact vibronic line shape that emphasizes the importance of using a quantum mechanical distribution of phase space initial conditions. For the unequal-frequency shifted harmonic oscillator model, we are no longer able to obtain an exact expression for the vibronic line shape in terms of independent deterministic classical trajectories. We show how one can rigorously account for corrections to the truncated Wigner approximation through nonlinear responses of the line shape function to momentum fluctuations along a classical trajectory and demonstrate the qualitative improvement in the resulting spectrum when the leading-order quantum correction is included. Finally, we numerically simulate absorption spectra of a highly anharmonic shifted Morse potential model. We find that, while finite quantization and the dissociation limit are captured with reasonable accuracy, there is a qualitative breakdown of the quasi-classical trajectory ensemble’s ability to describe the vibronic line shape when the relative shift in Morse potentials becomes large. The work presented here provides clarity on the origin of unphysical negative features known to contaminate absorption spectra computed with quasi-classical trajectory ensembles.
Tuning nonequilibrium phase transitions with inertia
In striking contrast to equilibrium systems, inertia can profoundly alter the structure of active systems. Here, we demonstrate that driven systems can exhibit effective equilibrium-like states with increasing particle inertia, despite rigorously violating the fluctuation–dissipation theorem. Increasing inertia progressively eliminates motility-induced phase separation and restores equilibrium crystallization for active Brownian spheres. This effect appears to be general for a wide class of active systems, including those driven by deterministic time-dependent external fields, whose nonequilibrium patterns ultimately disappear with increasing inertia. The path to this effective equilibrium limit can be complex, with finite inertia sometimes acting to accentuate nonequilibrium transitions. The restoration of near equilibrium statistics can be understood through the conversion of active momentum sources to passive-like stresses. Unlike truly equilibrium systems, the effective temperature is now density dependent, the only remnant of the nonequilibrium dynamics. Further, this density-dependent temperature can in principle introduce departures from equilibrium expectations, particularly in response to strong gradients. Our results provide additional insight into the effective temperature ansatz while revealing a mechanism to tune nonequilibrium phase transitions.
Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas
Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.
Energy-conserving contact dynamics of nonspherical rigid-body particles
Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.
Inducing thermodynamically blocked atomic ordering via strongly driven nonequilibrium kinetics
Ultrafast light-matter interactions enable inducing exotic material phases by promoting access to kinetic processes blocked in equilibrium. Despite potential opportunities, actively using nonequilibrium kinetics for material discovery is limited by the poor understanding on intermediate states of driven systems. Here, using single-pulse time-resolved imaging with x-ray free-electron lasers, we found intermediate states of photoexcited bismuth nanoparticles that showed kinetically reversed surface ordering during ultrafast melting. This entropy-lowering reaction was further investigated by molecular dynamics simulations to reveal that observed kinetics were thermodynamically buried in equilibrium, which emphasized the critical role of electron-mediated ultrafast free-energy modification in inducing exotic material phases. This study demonstrated that ultrafast photoexcitations of electrons provide an efficient strategy to induce hidden material phases by overcoming thermodynamic barriers via nonequilibrium reaction pathways.
Analysis of nonequilibrium air streams in the Ames 1-foot shock tunnel
Aerodynamic and thermodynamic calibration of nonequilibrium air flow in Ames 1-foot shock tunnel
Unified, Geometric Framework for Nonequilibrium Protocol Optimization
Controlling thermodynamic cycles to minimize the dissipated heat is a long-standing goal in thermodynamics, and more recently, a central challenge in stochastic thermodynamics for nanoscale systems. Here, we introduce a theoretical and computational framework for optimizing nonequilibrium control protocols that can transform a system between two distributions in a minimally dissipative fashion. These protocols optimally transport a system along paths through the space of probability distributions that minimize the dissipative cost of a transformation. Furthermore, we show that the thermodynamic metric—determined via a linear response approach—can be directly derived from the same objective function that is optimized in the optimal transport problem, thus providing a unified perspective on thermodynamic geometries. As a result, we investigate this unified geometric framework in two model systems and observe that our procedure for optimizing control protocols is robust beyond linear response.