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 379 records · Page 21

A fresh perspective on the role of band bending, and related contributors, in light-driven production of electricity and chemicals

It is widely known that semiconductor-based solar energy conversion could power our planet. This is in part because high-quality semiconductor structures are unrivalled in their ability to separate photogenerated electrons and holes. One effective approach to achieving this photoinduced charge separation relies on a phenomenon known as “band bending”. But details to justify why band bending results in photoinduced charge separation are more complex than often appreciated. This underappreciation is an impediment to the rational, hypothesis-driven design of next-generation approaches to solar energy conversion. Herein we show, by means of derivations rooted in physical chemistry, that several phenomena – not just band bending – can facilitate photoinduced charge separation, and that each is influenced by nonequilibrium species concentration and a parameter, such as diffusion coefficient or rate coefficient, that introduces dynamics. To help visualize the impact of each phenomenon, we introduce plots that depict their contributions as free energy, force, flux, force constant, and rate. We reveal that spatial dopant distributions that define band bending are predictors of initial photogenerated species transport rates. But charge separation alone does not guarantee high-efficiency operation. A photogenerated change in energy that is freely available to do useful work is also essential, and is strongly dependent on semiconductor optical properties and reaction kinetics. Notably, this information reveals that specificity of interfacial chemical reactions – even when they are not preceded by charge separation elsewhere – can result in efficient solar energy conversion. We expect that this tutorial will guide researchers in their pursuit to uncover new mechanisms for light to perform useful work.

14 SOLAR ENERGY↗

Driven-dissipative quantum mechanics on a lattice: Simulating a fermionic reservoir on a quantum computer

The driven-dissipative many-body problem remains one of the most challenging unsolved problems in quantum mechanics. The advent of quantum computers may provide a unique platform for efficiently simulating such driven-dissipative systems. But, there are many choices for how one can engineer the reservoir. One can simply employ ancilla qubits to act as a reservoir and then digitally simulate them via algorithmic cooling. A more attractive approach, which allows one to simulate an infinite reservoir, is to integrate out the bath degrees of freedom and describe the driven-dissipative system via a master equation, that can also be simulated on a quantum computer. In this work, we consider the particular case of noninteracting electrons on a lattice driven by an electric field and coupled to a fermionic thermostat. Then, we provide two different quantum circuits: the first one reconstructs the full dynamics of the system using Trotter steps, while the second one dissipatively prepares the final nonequilibrium steady state in a single step. We run both circuits on the IBM quantum experience. For circuit (i), we achieved up to five Trotter steps. When partial resets become available on quantum computers, we expect that the maximum simulation time can be significantly increased. Lastly, the methods developed here suggest generalizations that can be applied to simulating interacting driven-dissipative systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Adaptive time stepping for the two-time integro-differential Kadanoff-Baym equations

The nonequilibrium Green's function gives access to one-body observables for quantum systems. Of particular interest are quantities such as density, currents, and absorption spectra which are important for interpreting experimental results in quantum transport and spectroscopy. We present an integration scheme for the Green's function's equations of motion, the Kadanoff-Baym equations (KBE), which is both adaptive in the time integrator step size and method order as well as the history integration order. We analyze the importance of solving the KBE self-consistently and show that adapting the order of history integral evaluation is important for obtaining accurate results. To examine the efficiency of our method, we compare runtimes to a state-of-the-art fixed time step integrator for several test systems and show an order of magnitude speedup at similar levels of accuracy. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗

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.

Grossman, B.↗

Nonequilibrium radiative heating prediction method for aeroassist flowfields with coupling to flowfield solvers

A method for predicting radiation adsorption and emission coefficients in thermochemical nonequilibrium flows is developed. The method is called the Langley optimized radiative nonequilibrium code (LORAN). It applies the smeared band approximation for molecular radiation to produce moderately detailed results and is intended to fill the gap between detailed but costly prediction methods and very fast but highly approximate methods. The optimization of the method to provide efficient solutions allowing coupling to flowfield solvers is discussed. Representative results are obtained and compared to previous nonequilibrium radiation methods, as well as to ground- and flight-measured data. Reasonable agreement is found in all cases. A multidimensional radiative transport method is also developed for axisymmetric flows. Its predictions for wall radiative flux are 20 to 25 percent lower than those of the tangent slab transport method, as expected, though additional investigation of the symmetry and outflow boundary conditions is indicated. The method was applied to the peak heating condition of the aeroassist flight experiment (AFE) trajectory, with results comparable to predictions from other methods. The LORAN method was also applied in conjunction with the computational fluid dynamics (CFD) code LAURA to study the sensitivity of the radiative heating prediction to various models used in nonequilibrium CFD. This study suggests that radiation measurements can provide diagnostic information about the detailed processes occurring in a nonequilibrium flowfield because radiation phenomena are very sensitive to these processes.

Hartung, Lin C.↗

Model-Free Measurement of Local Entropy Production and Extractable Work in Active Matter

Time-reversal symmetry breaking and entropy production are universal features of nonequilibrium phenomena. Despite its importance in the physics of active and living systems, the entropy production of systems with many degrees of freedom has remained of little practical significance because the high dimensionality of their state space makes it difficult to measure. Here, in this work, we introduce a local measure of entropy production and a numerical protocol to estimate it. We establish a connection between the entropy production and extractability of work in a given region of the system and show how this quantity depends crucially on the degrees of freedom being tracked. We validate our approach in theory, simulation, and experiments by considering systems of active Brownian particles undergoing motility-induced phase separation, as well as active Brownian particles and E.coli in a rectifying device in which the time-reversal asymmetry of the particle dynamics couples to spatial asymmetry to reveal its effects on a macroscopic scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Examination of nonequilibrium effects in an ionized nitrogen flow

A thermo-chemical nonequilibrium ionized nitrogen flow in a shock tube is examined. A one-dimensional computational fluid dynamics code has been developed to study the flowfield incorporating a novel treatment of the electron energy equation. In the present approach, the electron pressure is included in the electron energy flux. In contrast, previous work has removed this term from the flux to facilitate the flux splitting. The code uses a five-species gas model for ionized nitrogen characterized by translational-rotational, vibrational, and electron-electronic temperatures. The results give good agreement with experimental data except for an inadequacy in the electronic source terms. A small difference is observed between the new and old splittings which increases with higher levels of ionization.

Hatfield, John A.↗

CFD Simulations of the IHF Arc-Jet Flow: Compression-Pad/Separation Bolt Wedge Tests

This paper reports computational analyses in support of two wedge tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using two different wedge models, each placed in a free jet downstream of a corresponding different conical nozzle in the Ames 60-MW Interaction Heating Facility. Each panel test article included a metallic separation bolt imbedded in Orion compression-pad and heatshield materials, resulting in a circular protuberance over a flat plate. The protuberances produce complex model flowfields, containing shock-shock and shock-boundary layer interactions, and multiple augmented heating regions on the test plate. As part of the test calibration runs, surface pressure and heat flux measurements on water-cooled calibration plates integrated with the wedge models were also obtained. Surface heating distributions on the test articles as well as arc-jet test environment parameters for each test configuration are obtained through computational fluid dynamics simulations, consistent with the facility and calibration measurements. The present analysis comprises simulations of the nonequilibrium flowfield in the facility nozzle, test box, and flowfield over test articles, and comparisons with the measured calibration data.

CFD↗

Characterizing skyrmion flow phases with principal component analysis

Principal component analysis (PCA) is a powerful method that can identify patterns in large, complex data sets by constructing low-dimensional order parameters from higher-dimensional feature vectors. There are increasing efforts to use space-and-time-dependent PCA to detect transitions in nonequilibrium systems that are difficult to characterize with equilibrium methods. Here, we demonstrate that feature vectors incorporating the position and velocity information of driven skyrmions moving through random disorder permit PCA to resolve different types of disordered skyrmion motion as a function of driving force and the ratio of the Magnus force to the dissipation. Since the Magnus force creates gyroscopic motion and a finite Hall angle, skyrmions can exhibit a greater range of flow phases than what is observed in overdamped driven systems with quenched disorder. We show that in addition to identifying previously known skyrmion flow phases, PCA detects several additional phases, including different types of channel flow, moving fluids, and partially ordered states. Guided by the PCA analysis, we further characterize the disordered flow phases to elucidate the different microscopic dynamics and show that the changes in the PCA-derived order parameters can be connected to features in bulk transport measures, including the transverse and longitudinal velocity-force curves, differential conductivity, topological defect density, and changes in the skyrmion Hall angle as a function of drive. We discuss how asymmetric feature vectors can be used to improve the resolution of the PCA analysis, and how this technique can be extended to find disordered phases in other nonequilibrium systems with time-dependent dynamics.

36 MATERIALS SCIENCE↗

Lightwave terahertz quantum manipulation of nonequilibrium superconductor phases and their collective modes

In this work, we present a gauge-invariant density matrix description of nonequilibrium superconductor (SC) states with spatial and temporal correlations driven by intense terahertz (THz) lightwaves. We derive superconductor Bloch-Maxwell equations of motion that extend Anderson pseudospin models to include the Cooper pair center-of-mass motion and electromagnetic propagation effects. We thus describe quantum control of dynamical phases, collective modes, quasiparticle coherence, and high nonlinearities during cycles of carrier wave oscillations, which relates to our recent experiments. Coherent photogeneration of a nonlinear supercurrent with a dc component, achieved via condensate acceleration by an effective lightwave field, dynamically breaks the equilibrium inversion symmetry. Experimental signatures include high harmonic light emission at equilibrium-symmetry-forbidden frequencies, Rabi-Higgs collective modes and quasiparticle coherence, and nonequilibrium moving condensate states tuned by few-cycle THz fields. We use such lightwaves as an oscillating accelerating force that drives strong nonlinearities and anisotropic quasiparticle populations to control and amplify different classes of collective modes, e.g., damped oscillations, persistent oscillations, and overdamped dynamics via Rabi flopping. Recent phase-coherent nonlinear spectroscopy experiments can be modeled by solving the full nonlinear quantum dynamics including self-consistent light-matter coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Plasma Assisted Combustion and Chemical Processing: Chapter 9 - Plasma Diagnostics

Plasma dynamics and chemistry have a broad range of timescales from picoseconds to milliseconds. In addition, it involves nonequilibrium energy transfer between electrons, ions, electronically and vibrationally excited states, radicals, intermediate species, and reactants and products as well as surface charges and chemistry. To understand plasma physics and chemistry, it is essential to conduct time and space resolved, quantitative detection of nonequilibrium temperature distributions, electron energy and number density, electric field, and species concentrations. There are enormous publications and review articles on this subject. The focus of this chapter is to be placed on the most recent progress in gas phase plasma properties and chemistry, especially on optical emission spectroscopy, laser absorption spectroscopy, Faraday rotational spectroscopy, Raman and Thompson scattering, femtosecond and picosecond (fs/ps) coherent anti-Stokes Raman scattering (CARS) spectroscopy, and electric field-induced second harmonic generation methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reduced Dynamics of Full Counting Statistics

We present a theory of modified reduced dynamics in the presence of counting fields. Reduced dynamics techniques are useful for describing open quantum systems at long emergent timescales when the memory timescales are short. However, they can be difficult to formulate for observables spanning the system and its environment, such as those characterizing transport properties. A large variety of mixed system--environment observables, as well as their statistical properties, can be evaluated by considering counting fields. Given a numerical method able to simulate the field-modified dynamics over the memory timescale, we show that the long-lived full counting statistics can be efficiently obtained from the reduced dynamics. We demonstrate the utility of the technique by computing the long-time current in the nonequilibrium Anderson impurity model from short-time Monte Carlo simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Full Facility Shock Frame Simulations of the Electric Arc Shock Tube

Radiative heating computations are performed for high speed lunar return experiments conducted in the Electric Arc Shock Tube (EAST) facility at NASA Ames Research Center. The nonequilibrium radiative transport equations are solved via NASA's in-house radiation code NEQAIR using flow field input from US3D flow solver. The post-shock flow properties for the 10 km/s Earth entry conditions are computed using the stagnation line of a blunt-body and a full facility CFD (Computational Fluid Dynamics) simulation of the EAST shock tube. The shocked gas in the blunt-body flow achieves a thermochemical equilibrium away from the shock front whereas EAST flow exhibits a nonequilibrium behavior due to strong viscous dissipation of the shock by boundary layer. The full-tube flow calculations capture the influence of the boundary layer on the shocked gas state and provide a realistic fluid dynamic input for the radiative predictions. The integrated radiance behind the shock is calculated in NEQAIR for wavelength regimes from Vacuum-UltraViolet (VUV) to InfraRed (IR), which are pertinent to the emission characteristics of high enthalpy shock waves in air. These radiance profiles are validated against corresponding EAST shots. The full-tube simulations successfully predict a sharp radiance peak at the shock front which gets smeared in the test data due to the spatial resolution in the measurements. The full facility based radiance behind the shock shows a slightly better match with the test data in the VUV and Red spectral regions, as compared to that from a blunt-body based predictions. The UV radiance is very similar for both geometries and under-predicts the test behavior. The IR test data matches better with the blunt-body based predictions where the full-tube simulations show a significant over-prediction.

Radiation↗

Implicit Thermochemical Nonequilibrium Flow Simulations on Unstructured Grids using GPUs

Thermochemical nonequilibrium flow simulation capabilities have been previously implemented, verified, and validated for central processing unit (CPU) systems in NASA’s unstructured-grid computational fluid dynamics solver FUN3D. Many exascale-class high-performance computing systems will rely on graphics processing unit (GPU) architectures for high throughput and energy efficiency; thus, CPU-based scientific computing software unable to effectively utilize these systems must be updated. In this work, we present a CUDA C++ implementation of FUN3D’s thermochemical nonequilibrium flow simulation capabilities targeting NVIDIA Tesla GPUs. An overview of the porting and optimization strategy is described and performance comparisons with other recent architectures are presented. Scaling to thousands of GPUs is demonstrated, yielding computational performance equivalent to that of several million CPU cores. The implementation enables efficient, high-fidelity, scale-resolving simulations of thermochemical nonequilibrium flows for many applications including atmospheric entry, hypersonics, and combustion.

GPU↗

Hydrodynamic generators in relativistic kinetic theory

Here, we resum the nonequilibrium gradient corrections to a single-particle distribution function evolved by the Boltzmann equation in the relaxation time approximation (RTA). We first study a system undergoing Bjorken expansion and show that, for a constant relaxation time, the exact solution of the RTA Boltzmann equation at late times (i.e., after the decay of nonhydrodynamic modes) generates the Borel resummed Chapman-Enskog series. Extending this correspondence to systems without Bjorken symmetry, we construct a ( 3 + 1 )-dimensional hydrodynamic generator for RTA kinetic theory, which is an integral representation of the Chapman-Enskog series in the limit of vanishing nonhydrodynamic modes. Relaxing this limit, we discover at earlier times a set of nonhydrodynamic modes coupled to the RTA Chapman-Enskog expansion. Including the dynamics of these nonhydrodynamic modes is shown to control the emergence of hydrodynamics as an effective field theory description of nonequilibrium fluids, which works well even for far-off-equilibrium situations where the Knudsen number is initially large.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analysis of thermo-chemical nonequilibrium models for carbon dioxide flows

The aerothermodynamics of thermochemical nonequilibrium carbon dioxide flows is studied. The chemical kinetics models of McKenzie and Park are implemented in separate three-dimensional computational fluid dynamics codes. The codes incorporate a five-species gas model characterized by a translational-rotational and a vibrational temperature. Solutions are obtained for flow over finite length elliptical and circular cylinders. The computed flowfields are then employed to calculate Mach-Zehnder interferograms for comparison with experimental data. The accuracy of the chemical kinetics models is determined through this comparison. Also, the methodology of the three-dimensional thermochemical nonequilibrium code is verified by the reproduction of the experiments.

Rock, Stacey G.↗