Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonequilibrium 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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electric Field Effects on Water and Ion Structure and Diffusion at the Orthoclase (001)–Water Interface

Understanding the electrochemical properties of mineral–water interfaces tends to rely upon electrical double layer (EDL) models, but these models are based on the assumption that electrostatic equilibrium is constantly maintained. In reality, interfacial reactions, ion diffusion, and their electrochemical signatures are based in nonequilibrium conditions of locally or globally imbalanced electrical fields where current EDL models have limited purview. In this work, we performed molecular dynamics (MD) simulations of the orthoclase (001) surface in contact with a 1 M NaCl aqueous solution under various electric fields, to explore the interplay between EDL structure and dynamics when perturbed by electric fields of different direction and strength, by confinement, and by different distributions of structural surface charge. The simulations showed that confinement between two opposing (001) surfaces led to the development of an induced field when the applied field was perpendicular to the surfaces and, as a result, to ionic diffusion coefficients that were independent of electric field strength. In contrast, when the applied field was parallel to the surfaces, confinement resulted in ionic diffusion coefficients that were more strongly dependent on the magnitude of the electric field than in bulk water. Differences in the density and distribution of aluminol groups on the two surfaces had a significant impact on how the interfacial structure and dynamics varied in the presence of an electric field. Notably, these differences resulted in an electro-osmotic flow with opposite directions at the two surfaces under parallel applied electric field. Overall, the MD simulations highlighted the importance of considering atomic-level structure and heterogeneities when developing models of the electrochemical properties of mineral–water interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

Quantum noise and the threshold of hearing

It is argued that the sensitivity of the ear reaches a limit imposed by the uncertainty principle. This is possible only if the receptor cell holds the detector elements in a special nonequilibrium state which has the same noise characteristics as a ground (T = 0 K) state. To accomplish this 'active cooling' the molecular dynamics of the system must maintain quantum mechanical coherence over the time scale of the measurement.

Bialek, W.↗
NASA NTRS1990

Theory of radiation from low velocity heated air

The application of hypersonic computational fluid dynamics models to low velocity vehicles is examined. Modeling aspects such as chemical kinetics, electronic excitation/deexcitation mechanisms, and existence of equilibrium versus nonequilibrium conditions in the flow were examined. Flowfield properties and in-band radiances in the wavelength region of .25 micron in the vicinity of the stagnation streamline are given for a three-inch hemisphere.

Levin, D. A.↗
NASA NTRS1994

A Class of High-Resolution Explicit and Implicit Shock-Capturing Methods

The development of shock-capturing finite difference methods for hyperbolic conservation laws has been a rapidly growing area for the last decade. Many of the fundamental concepts, state-of-the-art developments and applications to fluid dynamics problems can only be found in meeting proceedings, scientific journals and internal reports. This paper attempts to give a unified and generalized formulation of a class of high-resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock waves, perfect gases, equilibrium real gases and nonequilibrium flow computations. These numerical methods are formulated for the purpose of ease and efficient implementation into a practical computer code. The various constructions of high-resolution shock-capturing methods fall nicely into the present framework and a computer code can be implemented with the various methods as separate modules. Included is a systematic overview of the basic design principle of the various related numerical methods. Special emphasis will be on the construction of the basic nonlinear, spatially second and third-order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. Some perbolic conservation laws to problems containing stiff source terms and terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas-dynamics problems. The use of the Lax-Friedrichs numerical flux to obtain high-resolution shock-capturing schemes is generalized. This method can be extended to nonlinear systems of equations without the use of Riemann solvers or flux-vector splitting approaches and thus provides a large savings for multidimensional, equilibrium real gases and nonequilibrium flow computations.

Yee, H. C.↗
DOE OSTI2024

Thermochemical Nonequilibrium Modeling in a Continuous-Galerkin, Finite-Element Framework

The presented work discusses the implementation, verification, and validation of Park's two-temperature model in a scalable, computational fluid dynamics (CFD) code developed at the US Department of Energy's Oak Ridge National Laboratory (ORNL). The implementation of Park's two-temperature model was verified through 0D test cases involving an adiabatic reactor and a nitrogen thermal bath. The implementation was then validated through comparisons with other validated CFD codes and experimental data on a hypersonic cylinder and double cones. These are standard benchmark test cases for thermochemical non-equilibrium (TCNE) modeling, and all data are shared publicly. The verification and validation results showed that ORNL's in-house CFD code could model complex, high-speed flow problems with and without TCNE modeling. This work is essential for future research involving 3D shock wave/boundary layer interactions (SBLIs).

Nutter, Nicole↗
DOE OSTI2023

Out-of-equilibrium Majorana zero modes in interacting Kitaev chains

Here, we employ a time-dependent real-space local density-of-states method to study the movement and fusion of Majorana zero modes in the one-dimensional interacting Kitaev model, based on the time evolution of many-body states. We analyze the dynamics and both fusion channels of Majoranas using time-dependent potentials, either creating walls or wells. For fast moving Majoranas, we unveil nonequilibrium signatures of the “strong-zero-mode” operator (quasiparity degeneracy in the full spectrum) and its breakdown in the presence of repulsive Coulomb interactions. Focusing on forming a full electron after fusion, we also discuss the upper and lower limits on the Majorana speed needed to reduce nonadiabatic effects and to avoid poisoning due to decoherence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1972

Interrelated structures of the transport shock and collisional relaxation layer in a multitemperature, multilevel ionized gas

The gas dynamic structures of the transport shock and the downstream collisional relaxation layer are evaluated for partially ionized monatomic gases. Elastic and inelastic collisional nonequilibrium effects are taken into consideration. Three electronic levels are accounted for in the microscopic model of the atom. Nonequilibrium processes with respect to population of levels and species plus temperature are considered. By using an asymptotic technique the shock morphology is found on a continuum flow basis. The asymptotic procedure gives two distinct layers in which the nonequilibrium effects to be considered are different. A transport shock appears as the inner solution to an outer collisional relaxation layer in which the gas reaches local equilibrium. A family of numerical examples is displayed for different flow regimes. Argon and helium models are used in these examples.

Vinolo, A. R.↗
NASA NTRS1973

Interrelated structures of the transport shock and collisional relaxation layer in a multitemperature, multilevel ionized gas

The gas dynamic structures of the transport shock and the downstream collisional relaxation layer are evaluated for partially ionized monatomic gases. Elastic and inelastic collisional nonequilibrium effects are taken into consideration. In the microscopic model of the atom, three electronic levels are accounted for. By using an asymptotic technique, the shock morphology is found on a continuum flow basis. This procedure gives two distinct layers in which the nonequilibrium effects to be considered are different. A transport shock appears as the inner solution to an outer collisional relaxation layer. The results show four main interesting points: (1) on structuring the transport shock, ionization and excitation rates must be included in the formulation, since the flow is not frozen with respect to the population of the different electronic levels; (2) an electron temperature precursor appears at the beginning of the transport shock; (3) the collisional layer is rationally reduced to quadrature for special initial conditions, which (4) are obtained from new Rankine-Hugoniot relations for the inner shock.

Vinolo, A. R.↗
DOE OSTI2025

Electronic effects in radiation-induced collision cascades in nickel

The accurate treatment of electronic effects in multi-million-atom simulations of radiation-induced collision cascades is crucial for reliable predictions of primary radiation damage. In this work, we explore the fidelity of a recently developed two-temperature molecular dynamics model implementing an electron density-dependent coupling of electronic and atomic subsystems for cascade simulations in nickel. We show that the parameter-free model realistically captures the instantaneous energy losses during all stages of the highly nonequilibrium cascade process. Our simulations predict two distinct coupling regimes, corresponding to the rapid energy losses through electronic stopping in the early stages of the cascade and to the slow equilibration through the electron-phonon coupling mechanism in the later stages, without the use of separate models or coupling terms. The intermediate stage of the cascade dynamics displays a complex energy transfer between the subsystems, which cannot be interpreted by comparison to either electronic stopping or electron-phonon coupling theories. We therefore compare the predicted atomic mixing, which is sensitive to the energy losses during the intermediate cascade stage, with experimental ion beam mixing measurements. We find good agreement with the experiments, validating the coupling model for the intermediate stage of the cascade. Predictions of final defect numbers and cluster sizes are found in line with the results from conventional electronic stopping-based methods, while significantly reducing the theoretical uncertainty in the outcomes of conventional models stemming from arbitrary choices of thresholds for different coupling terms. Our results represent a notable improvement in cascade damage predictions in nickel, providing validation of the electron density-dependent coupling model for radiation damage simulations in general. The results lead us to propose an interpretation of the electronic energy losses in the intermediate regime of velocities, where we find an effectively nonlinear dissipation.

Crystal defects↗
DOE OSTI2020

Real-Time Time-Dependent Nuclear–Electronic Orbital Approach: Dynamics beyond the Born–Oppenheimer Approximation

The quantum mechanical treatment of both electrons and nuclei is crucial in nona- diabatic dynamical processes such as proton-coupled electron transfer. The nuclear electronic orbital (NEO) method provides an elegant framework for including nuclear quantum effects beyond the Born-Oppenheimer approximation. To enable the study of nonequilibrium properties, we derive and implement a real-time NEO (RT-NEO) approach based on time-dependent Hatree-Fock or density functional theory, in which both the electronic and nuclear degrees of freedom are propagated in a time-dependent variational framework. Nuclear and electronic spectral features can be resolved from the time-dependent dipole moment computed using the RT-NEO method. The test cases demonstrate the dynamical interplay between the quantum nuclei and the electrons through vibronic coupling. Moreover, vibrational excitation in the RT-NEO approach is demonstrated by applying a resonant driving eld, and electronic excitation is demonstrated by simulating excited state intramolecular proton transfer. This work shows that the RT-NEO approach is a promising tool to study nonadiabatic quantum dynamical processes within a time-dependent variational description for the coupled electronic and nuclear degrees of freedom.

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

Dynamic Motions of Topological Defects in Nematic Liquid Crystals under Spatial Confinement

Abstract Topological defects (TDs) have become a sensational topic due to their significant influence on the unusual optical and physiochemical characteristics of the material. To facilitate their application across a wide range of disciplines it is desirable to analyze and gain fundamental understanding of TDs in both equilibrium and nonequilibrium systems. Liquid crystals (LCs) are considered an ideal system for study given the direct visualization of TDs and a straightforward analyzation process. In addition to the equilibrium morphology of LC TDs, it is also of great interest to track and control the formation and annihilation of defects during thermodynamic processes. However, controlling the dynamic behavior of formed defects remains a challenge. Here, nematic LCs confined in a cell exhibiting surfaces with periodic anchoring conditions containing surface topography are relied on. The effects of patterned surface characteristics such as width, periodicity, and degree of curvature on defects dynamic motion, stabilization, and annihilation are explored. The computational experiments recapitulate the TDs transition path and provide free energy‐based predictions of critical distances for defect annihilation. Taken together, this simple approach offers a promising opportunity to control the dynamics of TDs in LCs through chemical patterned surfaces with topography.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Improving noisy free-energy measurements by adding more noise

Estimating free-energy differences using nonequilibrium work relations, such as the Jarzynski equality, is hindered by poor convergence when work fluctuations are large. For systems governed by overdamped Langevin dynamics, we propose the counterintuitive approach of adding noise in order to increase the precision of such calculations. Here, by introducing additional stochastic fluctuations to the system and rescaling its potential energy accordingly, we leave the thermodynamics of the system unchanged while increasing its relaxation rate. For a given time-dependent protocol this modification reduces the dissipated reduced work, leading to more accurate free-energy estimates. The method is designed to be used in experiment, and we illustrate its operation using computer simulations applied to two model systems. However, the regime of applicability of this strategy is likely limited, because it requires control of the system's potential energy in a way that is feasible in only a few experimental settings.

Whitelam, Stephen [Lawrence Berkeley National Labo↗
NASA NTRS1987

Study of non-equilibrium transport phenomena

Nonequilibrium phenomena due to real gas effects are very important features of low density hypersonic flows. The shock shape and emitted nonequilibrium radiation are identified as the bulk flow behavior parameters which are very sensitive to the nonequilibrium phenomena. These parameters can be measured in shock tubes, shock tunnels, and ballistic ranges and used to test the accuracy of computational fluid dynamic (CFD) codes. Since the CDF codes, by necessity, are based on multi-temperature models, it is also desirable to measure various temperatures, most importantly, the vibrational temperature. The CFD codes would require high temperature rate constants, which are not available at present. Experiments conducted at the NASA Electric Arc-driven Shock Tube (EAST) facility reveal that radiation from steel contaminants overwhelm the radiation from the test gas. For the measurement of radiation and the chemical parameters, further investigation and then appropriate modifications of the EAST facility are required.

Sharma, Surendra P.↗
NASA NTRS1987

Upwind and symmetric shock-capturing schemes

The development of numerical methods for hyperbolic conservation laws has been a rapidly growing area for the last ten years. Many of the fundamental concepts and state-of-the-art developments can only be found in meeting proceedings or internal reports. This review paper attempts to give an overview and a unified formulation of a class of shock-capturing methods. Special emphasis is on the construction of the basic nonlinear scalar second-order schemes and the methods of extending these nonlinear scalar schemes to nonlinear systems via the extact Riemann solver, approximate Riemann solvers, and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows is discussed. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas dynamics problems.

Yee, H. C.↗
DOE OSTI2025

Interpreting dynamic-compression experiments to uncover the time dependence of freezing: Application to gallium

Using pulsed-power magnetic field sources to compress gallium to gigapascal pressures on nanosecond timescales, we report here experiments on shockless dynamic compression of a liquid metal. Time-resolved velocimetry data reveal signatures of rapid freezing from a metastable liquid state, and we demonstrate that the kinetics of this nonequilibrium solidification can be accurately simulated with a computational modeling framework we have developed in previous studies, where classical nucleation theory is coupled with hydrodynamics. Notably, velocity traces in some of our experiments show evidence of a phase transition, while others do not, even though other types of evidence suggest that solidification may be occurring in all of them. We explain how predictions made by our models regarding the presence or absence of these phase-transition signatures motivated additional experiments that later confirmed the theoretical predictions. Our analysis shows that due to the rapid, quasi-isentropic nature of the loading path, our experiments were able to compress liquid gallium to metastable states that are undercooled below the equilibrium melt temperature by more than 300 K and exhibit pressures that approach five times the equilibrium melt pressure. The understanding gained in this study should form the basis for future dynamic-compression experiments aimed at interrogating melt curves at high pressures.

gallium↗
NASA NTRS1989

An efficient, explicit finite-rate algorithm to compute flows in chemical nonequilibrium

An explicit finite-rate code was developed to compute hypersonic viscous chemically reacting flows about three-dimensional bodies. Equations describing the finite-rate chemical reactions were fully coupled to the gas dynamic equations using a new coupling technique. The new technique maintains stability in the explicit finite-rate formulation while permitting relatively large global time steps.

Palmer, Grant↗
NASA NTRS1997

Noise and Dynamical Pattern Selection in Solidification

The overall goal of this project was to understand in more detail how a pattern-forming system can adjust its spacing. "Pattern-forming systems," in this context, are nonequilibrium contina whose state is determined by experimentally adjustable control parameter. Below some critical value of the control system then has available to it a range of linearly stable, spatially periodic steady states, each characterized by a spacing which can lie anywhere within some band of values. These systems like directional solidification, where the solidification front is planar when the ratio of growth velocity to thermal gradient is below its critical value, but takes on a cellular shape above critical. They also include systems without interfaces, such as Benard convection, where it is the fluid velocity field which changes from zero to something spatially periodic as the control parameter is increased through its critical value. The basic question to be addressed was that of how the system chooses one of its myriad possible spacings when the control parameter is above critical, and in particular the role of noise in the selection process. Previous work on explosive crystallization had suggested that one spacing in the range should be preferred, in the sense that weak noise should eventually drive the system to that spacing. That work had also suggested a heuristic argument for identifying the preferred spacing. The project had three main objectives: to understand in more detail how a pattern-forming system can adjust its spacing; to investigate how noise drives a system to its preferred spacing; and to extend the heuristic argument for a preferred spacing in explosive crystallization to other pattern-forming systems.

Kurtze, Douglas A.↗
DOE OSTI2020

Transport phenomena in electrolyte solutions: Nonequilibrium thermodynamics and statistical mechanics

The theory of transport phenomena in multicomponent electrolyte solutions is presented here through the integration of continuum mechanics, electromagnetism, and nonequilibrium thermodynamics. The governing equations of irreversible thermodynamics, including balance laws, Maxwell's equations, internal entropy production, and linear laws relating the thermodynamic forces and fluxes, are derived. Green–Kubo relations for the transport coefficients connecting electrochemical potential gradients and diffusive fluxes are obtained in terms of the flux–flux time correlations. Furthermore, the relationship between the derived transport coefficients and those of the Stefan–Maxwell and infinitely dilute frameworks are presented, and the connection between the transport matrix and experimentally measurable quantities is described. To exemplify the application of the derived Green–Kubo relations in molecular simulations, the matrix of transport coefficients for lithium and chloride ions in dimethyl sulfoxide is computed using classical molecular dynamics and compared with experimental measurements.

42 ENGINEERING↗