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.

SourceResearch connectionFrom yearThrough year

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

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS1995

A comparison of two algorithms for simulating collisionless systems

Two completely different simulation algorithms are compared by applying them to the same stellar dynamical problems: one is a self-consistent field (SCF) method for solving Poisson's equation and the other is a phase-space method for integrating the collisionless Boltzmann equation. We consider simulations of spherical stellar systems which are initially far from equilibrium and relax to their final states by gravitational collapse. The initial conditions consist of either uniform-density spheres or nonequilibrium models having Plummer density profiles, in which velocity dispersions are assigned according to given virial ratios. If a few tens of radial expansion terms with hundreds of thousands of particles are used in the SCF code, excellent agreement is found between the results it generates and those obtained with the phase-space solver, provided that a sufficiently large number of grid cells are employed with the latter. These findings imply that for simulating collisionless systmes over many dynamical times, the SCF approach based on sampling phase space is competitive with the approach treating phase space as a continuous fluid. The results of our tests make it possible to estimate the number of particles and basis functions required in situations like those modeled. Limitations of the SCF method and the choice of an optimal set of basis functions are also discussed.

Hozumi, Shunsuke↗

DOE OSTI2023

Fluctuation-driven excess noise near superconducting phase transition

We discuss intrinsic mechanisms of nonequilibrium excess noise in superconducting devices and transition edge sensors. In particular, we present an overview of fluctuation-driven contributions to the current noise in the vicinity of the superconducting transition. We argue that sufficiently close to the critical temperature fluctuations of conductivity may become correlated provided that the rate of quasiparticle relaxation is slow as compared to dynamics of superconducting fluctuations. In this regime, fluctuations of conductivity adiabatically follow the fluctuations of the electron distribution. This leads to a substantial enhancement of current noise. Here, the corresponding spectral power density of noise has a Lorentzian shape in the frequency domain while its magnitude scales proportionally to the inelastic relaxation time. It also sensitively depends on the dephasing and Ginzburg–Landau timescales. Further estimates suggest that this mechanism dominates over the conventional temperature fluctuations in the same range of parameters. To describe these effects microscopically, we use the nonequilibrium Keldysh technique in the semiclassical approximation of superconductivity with Boltzmann–Langevin random forces to account for correlations of fluctuations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Plasmonic Energetic Electrons Drive CO 2 Reduction on Defective Cu 2 O

Plasmonic photoreduction of CO 2 is valuable for decarbonization and producing value-added chemicals. However, insights into the mechanisms of this reaction remain elusive, particularly regarding the roles of structural defects and their interplay with the nonequilibrium charge carriers. Here, we report density functional theory calculations on Cu 2 O, a prototype photocatalyst, through which we investigate CO 2 reduction over three defected facets to reveal the interfacial charge transfer and bond dynamics under plasmonic excitation. We find that the activation barrier of C–O bond cleavage decreases from 3.2 to about 1 eV, assisted by oxygen vacancies, and that the remaining barrier can be further reduced or eliminated at the plasmon-excited states when Cu 2 O is integrated with plasmonic metals. The regeneration of oxygen vacancies (by H 2 to form water) on Cu 2 O to complete the catalysis cycle is feasible and not affected by the energetic electrons. Furthermore, our calculations thus show the important synergistic effect of energetic electrons and point defects to promote CO 2 reduction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Concerted Electron-Ion Transport by Polyacrylonitrile Elucidated with Reactive Deep Learning Potentials

Charge transport in polymers, such as polyacrylonitrile (PAN), is crucial for electronics and energy storage. For instance, PAN can transport cations e.g., Li + , by facilitating dynamic cation-nitrile coordination in batteries. However, little is known regarding the underlying role of complex reactive polymer configurations. Herein, we develop a deep-learning potential, trained on ab initio energies and forces of nonequilibrium reactive PAN configurations, to unravel the kinetics of PAN cyclization initiated by a nucleophile (OH – dissociated from LiOH) attacking the terminal nitrile carbon. We find, based on the reaction free-energetics, rates, and charge analysis, that the nucleophile attack producing the first ring is the rate-limiting step, which subsequently triggers Li + -coupled electron transfer along the PAN backbone, causing ∼10 4 times faster sequential ring-formation of the remaining nitriles. PAN’s extended configurations, where dipolar and H-bonding interactions are minimal, enable such rapid kinetics. By validating our computational findings with IR and NMR experiments, we establish a pathway for designing reactive polymers with enhanced charge transport for energy applications.

Chahal-Crockett, Rajni [Oak Ridge National Laborat↗

DOE OSTI2026

Inductance meets memory in a quantum magnet

Orbital degrees of freedom offer a largely untapped route to emergent dynamical phenomena in correlated quantum materials. However, it remains unclear whether collective orbital states can intrinsically generate both reactive and memory functionalities in a bulk system. Here we show that in the ferrimagnet Mn₃Si₂Te₆, nonequilibrium reconfiguration of chiral orbital currents produces both emergent inductance and nonvolatile memristance as intrinsic properties of a single crystal. At low frequency and under a magnetic field along the c axis, coherent orbital-current domains generate robust clockwise inductive I-V loops. At higher frequency and low field, current-driven first-order reconfiguration leads to incomplete reversal and metastable trapping, producing an intrinsic electromotive force and a finite remanent voltage at zero current. These results establish orbital currents as a class of quantum state variables that encode both reactive and memory functionalities, opening routes toward intrinsically reconfigurable and energy-efficient electronic systems.

Cao, Tristan [University of Colorado, Boulder]↗

NASA NTRS1989

Review of a code development and calibration program in support of the aeroassist flight experiment

The code development and calibration program for the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) is reviewed, with emphasis directed toward support of the Aeroassist Flight Experiment (AFE). The flight project is designed to obtain data which will be used in the validation of computational fluid dynamic approximation methods. Comparisons between experimental data and numerical simulation focus on perfect-gas tests over a scale model of the AFE and on flight and ground tests which challenge some aspect of the thermochemical nonequilibrium model. In the first case, the gas model is simple, but the grid-related problems of defining the real vehicle are present. In the second case, the vehicle geometries are simple, but thermochemical processes must be modeled correctly in order to compare with the experimental data. These comparisons are described as calibration runs because they test elements of the numerical simulation, but no single data set adequately simulates the full-scale AFE flight conditions. Comparisons between computation and experiment over a broad range of data sets show generally good agreement, though some aspects of the numerical model require further development.

Gnoffo, Peter A.↗

DOE OSTI2025

Unravelling Microstructure Selection in an Additively Manufactured Eutectic High‐Entropy Alloy

High-entropy alloys (HEAs) are promising candidates for advanced structural applications due to their excellent mechanical properties. Additive manufacturing (AM), with its rapid solidification conditions, enables the creation of unique nonequilibrium microstructures. To fully leverage the synergy between AM and HEAs, understanding how processing affects structure and properties is essential. Here, how solidification rate influences microstructure evolution and phase transformation pathway in laser additively manufactured AlCrFe2Ni2 eutectic HEAs is investigated. By increasing the laser scan speed and hence the solidification rate, distinct solidification modes evolving from coupled eutectic to anomalous eutectic and eventually to single-phase solidification are revealed. These transitions result in distinct microstructures and a wide range of mechanical properties. Thermodynamic modeling and molecular dynamics simulations reveal that low cooling rates allow for sufficient atomic diffusion and phase separation, facilitating coupled eutectic growth. In contrast, rapid cooling suppresses diffusion and destabilizes the solid–liquid interface, promoting anomalous or single-phase solidification. This integrated experimental and computational approach provides a multiscale understanding of solidification mechanisms in HEAs and underscores how kinetic effects can over-ride thermodynamic predictions under nonequilibrium conditions. Furthermore, these results demonstrate that AM can serve as a powerful tool to design HEAs with tailored microstructures and properties.

36 MATERIALS SCIENCE↗

NASA NTRS2001

Turbulent Flame Processes Via Diffusion Flame-Vortex Ring Interactions

Flame-vortex interactions are canonical configurations that can be used to study the underlying processes occurring in turbulent reacting flows. This configuration contains many of the fundamental aspects of the coupling between fluid dynamics and combustion that could be investigated with more controllable conditions than are possible under direct investigations of turbulent flames. Diffusion flame-vortex ring interaction contains many of the fundamental elements of flow, transport, combustion, and soot processes found in turbulent diffusion flames. Some of these elements include concentrated vorticity, entrainment and mixing, strain and nonequilibrium phenomena, diffusion and differential diffusion, partial premixing and diluent effects, soot formation and oxidation, and heat release effects. Such simplified flowfield allows the complex processes to be examined more closely and yet preserving the physical processes present in turbulent reacting flows. Furthermore, experimental results from the study of flame-vortex interactions are useful for the validation of numerical simulations and more importantly to deepen our understanding of the fundamental processes present in reacting flows. Experimental and numerical results obtained under microgravity conditions of the diffusion flame-vortex ring interaction are summarized in this paper. Results are obtained using techniques that include Flame Luminosity Imaging (FLI), Laser Soot-Mie Scattering (LSMS), Computational Fluid Dynamics and Combustion (CFDC), and Diode Laser Spectroscopy/Iterative Temperature with Assumed Chemistry (DLS/ITAC).

Dahm, Werner J. A.↗

NASA NTRS1988

Computational fluid dynamics and aerothermodynamics

Approximations applicable to the radiating, reacting, and conducting stagnation region of a hypervelocity vehicle were incorporated into a method for rapidly obtaining approximate solutions. This solution utilizes a coordinate system based upon the origin of the radiative losses and includes in a phenomenologically correct manner the effects of chemical and thermal nonequilibrium, and nonequilibrium, nongray radiative transfer. Results were presented which demonstrate the usefulness of the method and indicate which radiation parameters require further study and definition. Excellent comparisons were obtained with published results for the Fire2 data. An axisymmetric nonequilibrium inverse method was modified and extended and used to investigate and compare various vibration dissociation chemistry coupling models and radiative heat transfer approximations. The similarities, differences, and consequences of using these models in the Aero-assist Orbital Transfer Vehicles flight regime will be discussed.

Carlson, Leland A.↗

NASA NTRS1990

A three-dimensional upwind parabolized Navier-Stokes code for chemically reacting flows

A new upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the three-dimensional flow of chemically reacting air around hypersonic vehicles. The code is a modification of the perfect gas, three-dimensional UPS code of Lawrence et al. (1986) which has been extended in the present study to permit the calculation of hypersonic, viscous flows in chemical nonequilibrium. The algorithm solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that has been modified to account for real gas effects. The present code solves the fluid dynamic and species continuity equations in a loosely-coupled manner. The fluid medium is assumed to be a chemically reacting mixture of thermally perfect (but calorically imperfect) gases in thermal equilibrium. Results are presented for the hypersonic laminar flow over a cone at 0- and 10-deg angles of attack and for a generic hypersonic vehicle. Calculations are performed assuming either perfect gas, equilibrium air, or finite-rate chemistry.

Buelow, Philip E.↗

DOE OSTI2024

Influence of the Dielectric Constant on the Ionic Current Rectification of Bipolar Nanopores

In this paper, we investigate how the dielectric constant, ϵ, of an electrolyte solvent influences the current rectification characteristics of bipolar nanopores. It is well recognized that bipolar nanopores with two oppositely charged regions rectify current when exposed to an alternating electric potential difference. Here, we consider dilute electrolytes with NaCl only and with a mixture of NaCl and charged nanoparticles. These systems are studied using two levels of description, all-atom explicit water molecular dynamics (MD) simulations and coarse-grained implicit solvent MD simulations. The charge density and electric potential profiles and current-voltage relationship predicted by the implicit solvent simulations with ϵ = 11.3 show good agreement with the predictions from the explicit water simulations. Under nonequilibrium conditions, the predictions of the implicit solvent simulations with a dielectric constant closer to the one of bulk water are significantly different from the predictions obtained with the explicit water model. Further, these findings are closely aligned with experimental data on the dielectric constant of water when confined to nanometric spaces, which suggests that ϵ decreases significantly compared to its value in the bulk. Moreover, the largest electric current rectification is observed in systems containing nanoparticles when ϵ = 78.8. Using enhanced sampling, we have shown that this larger rectification arises from the presence of a significantly deeper minimum in the free energy of the system with a larger ϵ, and when a negative voltage bias is applied. Since implicit solvent models and mean-field continuum theories are often used to design Janus membranes based on bipolar nanopores, this work highlights the importance of properly accounting for the effects of confinement on the dielectric constant of the electrolyte solvent. The results presented here indicate that the dielectric constant in implicit solvent simulations may be used as an adjustable parameter to approximately account for the effects of nanometric confinement on aqueous electrolyte solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1995

Formation of Hydrocarbons in the Outflows from Red Giants

The formation of hydrocarbons in the oxygen-rich outflows from red giants was studied. The existence of organic molecules in such outflows has been known for several years; however, their surprisingly high abundances has been a mystery since all of the carbon had been thought to be irretrievably locked up in CO, the most strongly bound molecule. CO is the first molecule to form from the atoms present in the star's extended atmosphere, and as strong stellar winds drive a cooling outflow, dust grains condense out. In oxygen-rich outflows, the dust is thought to be composed mainly of silicates and other metal oxides. Perhaps the noble metals can condense out in metallic form, in particular the relatively abundant transition metals iron and nickel. We proposed that perhaps the carbon reservoir held as CO can be accessed through a catalytic process involving the chemisorption of CO and H2 onto grains rich in metallic iron. CO and H2 are the two most abundant molecules in circumstellar outflows, and they both are known to dissociate on transition metal surfaces at elevated temperatures, freeing carbon to form organic molecules such as methane. We believe methane is a precursor molecule to the organics observed in oxygen-rich red giants. We have developed a nonequilibrium numerical model of a surface chemical (catalytic) process. Based on this model, we believe that methane can be formed under the conditions present in circumstellar outflows. Although the methane formation rates are exceptionally low under these conditions, over dynamical timescales, a significant amount of CO can be converted to methane and driven further out in the envelope, explaining the presence of organics there.

Roberge, Wayne↗

NASA NTRSDate not supplied

Aerothermal Analysis of the Dragonfly Titan Entry

The Dragonfly mission will send a rotorcraft lander to the surface of Saturn’s moon Titan as part of the New-Frontiers program. This will be the first spacecraft to land on Titan since the Huygens probe’s descent in January of 2005 and only the second spacecraft to enter Titan’s atmosphere. The Dragonfly entry capsule is significantly larger than the Huygens probe and will experience higher aerothermal environments. This poster will provide an overview of the aerothermodynamic models used in the Dragonfly aeroshell design process. Titan provides a unique entry environment that has several fundamental differences from the environments at the more common entry destinations of Earth and Mars. The lack of atmospheric oxygen significantly reduces heatshield recession (expected to be negligible), which simplifies material response analysis but also reduces the efficiency of the thermal protection system. The atmos-phere is composed of mostly nitrogen and trace amounts of methane. During hypersonic entry this methane dissociates and leads to the formation of the molecule CN, which is known to radiate strongly in the shock layer. Current Dragonfly heating predictions estimate that radiation contributes about 50% of the total heat flux along the forebody and up to 90% on the aftbody, making radiation a key element of aerothermal analysis. To properly address the importance of radiation, a unique radiative-heating correlation was developed that incorporates normal-shock equilibrium chemistry predictions based on work from the Mars 2020 program. These new heating correlations agree well with Computational Fluid Dynamics (CFD) predictions and enable large-scale trajectory analyses by providing rapid heating predictions along candidate trajectories. Dragonfly aerothermal environments are generated with the Data Parallel Line Relaxation (DPLR) finite-volume Navier-Stokes solver alongside the Nonequilibrium Radiative Transport and Spectra (NEQAIR) radiation transport code. DPLR produces shock-aligned flowfields that incorporate the effects of both chemical and thermal nonequilibrium while NEQAIR solves the radiation transport equations along lines-of-sight through the flowfield solutions to predict the radiative heating on the aeroshell. Titan’s atmosphere produces complex thermal and chemical nonequilibria and so classical radiation approximations that extrapolate from a single line-of-sight, such as the Tangent-Slab method, were found to overpredict peak radiative heat fluxes by 10-20%. Therefore, all radiative heating is directly calculated using 3-D transport for both the forebody and aftbody to reduce unnecessary conservatism. A 21-species finite-rate chemistry model, including electrons and ionized species, is used in the CFD to ensure accurate vehicle heating predictions (the level of ionization is particularly for accurate radiative flux predictions). This poster will cover: 1) The models and assumptions used to simulate entry into Titan’s atmosphere, 2) Dragonfly aerothermal engineering correlations alongside statistical heating values, and 3) Driving aerothermal considerations for a Titan entry such as CN formation and radiation-flowfield coupling.

EDL↗

NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

NASA NTRSDate not supplied

Improved Heat Transfer Prediction for High-Speed Flows over Blunt Bodies using Adaptive Mixed-Element Unstructured Grids

Aerothermodynamic prediction is important for the design and analysis of many aerospace vehicles. Computational Fluid Dynamics (CFD) tools used for these predictions commonly rely on structured shock-aligned grids due to the strong shocks exhibited at high speeds. This work details the incorporation of an improved HLLE++ scheme for both perfect gas and thermochemical nonequilibrium flows which robustly captures shock waves on both non-shock-aligned and shock-aligned grids. The HLLE++ scheme reduces the impact of eigenvalue limiting on flow solutions which is required for Roe's scheme for high Mach number flows. An adapted grid approach is also demonstrated using mixed-element unstructured grids to improve heat prediction. Thin prismatic boundary layers are used with tetrahedra in the farfield and for capturing shock waves. Various results are presented for hypersonic flows over blunt bodies including cylinders, spheres, and capsules.

CFD↗

NASA NTRSDate not supplied

Improved Heat Transfer Prediction for High-Speed Flows over Blunt Bodies using Adaptive Mixed-Element Unstructured Grids

CFD↗

DOE OSTI2025

Ultrafast simultaneous manipulation of multiple ferroic orders through nonlinear phonon excitation

Recent experimental studies have demonstrated the possibility of utilizing strong terahertz pulses to manipulate individual ferroic orders on pico- and femtosecond timescales. Here, we extend these findings and showcase the simultaneous manipulation of multiple ferroic orders in BiFeO 3 , a material that is both ferroelectric and antiferromagnetic at room temperature. We find a concurrent enhancement of ferroelectric and antiferromagnetic second-harmonic generation (SHG) following the resonant excitation of a high-frequency fully-symmetric phonon mode. Based on first-principles calculations and phenomenological modeling, we ascribe this observation to the inherent coupling of the two ferroic orders to the nonequilibrium distortions induced in the crystal lattice by nonlinearly driven phonons. Our finding highlights the potential of nonlinear phononics as a technique for manipulating multiple ferroic order parameters at once. In addition, this approach provides a promising avenue to studying the dynamical magnetic and polarization behavior, as well as their intrinsic coupling, on ultrashort timescales.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Detecting low-energy interactions and the effects of energy accumulation in materials

Elusive as dark matter particles are, they are not the only entities that can produce small energy releases deep inside the most sensitive detectors. Cosmogenic and residual radioactivity as well as other factors can produce and slowly accumulate long-living excitations and defect configurations in materials. Unsteady and avalanche-like releases of accumulated energy can limit the dynamic range of detectors and the sensitivity of experiments. This type of mechanism, while not widely discussed in the dark matter and neutrino detectors community, got a lot of attention in discussions about systems with energy flow in the context of nonequilibrium thermodynamic. Here this paper explores the avalanche hypothesis using published data on low-energy background in detectors and available condensed-matter physics information on excitations, defects, and energetic molecules which can be produced in detector materials and examine ways to use these models to understand the origin of low-energy background in dark matter and other detectors and for suppression of these backgrounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗