Search NASA⌕ Search

SEARCH · Search NASA

Results for “non-equilibrium gas 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗
NASA NTRS2004

The Interaction of Hot and Cold Gas in the Disk and Halo of Galaxies

Most of the thermal energy in the Galaxy and perhaps most of the baryons in the Universe are found in hot (log T approximately 5.5 - 7) gas. Hot gas is detected in the local interstellar medium, in supernova remnants (SNR), the Galactic halo, galaxy clusters and the intergalactic medium (IGM). In our own Galaxy, hot gas exists in large superbubbles up to several hundred pc in diameter that locally dominate the interstellar medium (ISM) and determine its thermal and dynamic evolution. While X-ray observations using ROSAT, Chandra and XMM have allowed us to make dramatic progress in mapping out the morphology of the hot gas and in understanding some of its spectral characteristics, there remain fundamental questions that are unanswered. Chief among these questions is the way that hot gas interacts with cooler phase gas and the effects these interactions have on hot gas energetics. The theoretical investigations we proposed in this grant aim to explore these interactions and to develop observational diagnostics that will allow us to gain much improved information on the evolution of hot gas in the disk and halo of galaxies. The first of the series of investigations that we proposed was a thorough exploration of turbulent mixing layers and cloud evaporation. We proposed to employ a multi-dimensional hydrodynamical code that includes non-equilibrium ionization (NEI), radiative cooling and thermal conduction. These models are to be applied to high velocity clouds in our galactic halo that are seen to have O VI by FUSE (Sembach et ai. 2000) and other clouds for which sufficient constraining observations exist.

Slavin, Jonathan↗
NASA NTRSDate not supplied

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗
DOE OSTI2021

Role of surface diffusion in formation of unique reactivity for graphite oxidation: Time-resolved measurements in a pulsed diffusion reactor

Quantification of oxidation kinetics is essential to develop graphitic materials for diverse applications: from refractories found in gas-cooled nuclear reactors to catalysts needed for chemical manufacturing. In this work, using well-defined highly oriented pyrolytic graphite, low-pressure isotopic transient experiments combined with controlled annealing periods, we resolve the role of surface diffusion and quantify oxidation kinetics with nanomole-precision. We observe an unexpected increase in reactivity following annealing which is explained by the role of surface diffusion increasing the probability for trapping mobile oxygen at more reactive edge sites. Here, the locus of adsorption and spillover to the basal plane is distinct from the trapping location creating a more active oxygen species. Isotopic products reflect the population dynamics of oxygen added at the edge and surface diffusion that relocates basal plane oxygen to more reactive edge sites. Since this process proceeds in parallel with direct oxidation reactions, it is not likely to be observed using steady-state or conventional ‘bulk’ characterization techniques. Our unique time-resolved non-equilibrium measurement in a well-defined transport regime, enables observation of three distinct behaviors: short-term deactivation due to the balance of rates in oxygen supply/product formation, reactivity increases due to surface diffusion and longer-term reactivity increase with oxygen accumulation.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Fluid dynamic modeling and numerical simulation of low-density hypersonic flow

The concept of a viscous shock-layer and several related versions of continuum theories/methods are examined for their adequacy as a viable framework to study flow physics and aerothermodynamics of relevance to sustained hypersonic flights. Considering the flat plate at angle of attack, or the wedge, as a generic example for the major aerodynamic component of a hypersonic vehicle, the relative importance of the molecular-transport effects behind the shock (in the form of the 'shock slip') and the wall-slip effects are studied. In the flow regime where the shock-transition-zone thickness remains small compared to the shock radius of curvature, a quasi-one-dimensional shock structure under the Burnett/thirteen-moment approximation, as well as particulate/collisional models, can be consistently developed. The fully viscous version of the shock-layer model is shown to provide the crucial boundary condition downstream the shock in this case. The gas-kinetic basis of the continuum description for the flow behind the bow shock, and certain features affecting the non-equilibrium flow chemistry, are also discussed.

Cheng, H. K.↗
NASA NTRS2004

A Massively Parallel Hybrid Dusty-Gasdynamics and Kinetic Direct Simulation Monte Carlo Model for Planetary Applications

In order to understand the global structure, dynamics, and physical and chemical processes occurring in the upper atmospheres, exospheres, and ionospheres of the Earth, the other planets, comets and planetary satellites and their interactions with their outer particles and fields environs, it is often necessary to address the fundamentally non-equilibrium aspects of the physical environment. These are regions where complex chemistry, energetics, and electromagnetic field influences are important. Traditional approaches are based largely on hydrodynamic or magnetohydrodynamic (MHD) formulations and are very important and highly useful. However, these methods often have limitations in rarefied physical regimes where the molecular collision rates and ion gyrofrequencies are small and where interactions with ionospheres and upper neutral atmospheres are important. At the University of Michigan we have an established base of experience and expertise in numerical simulations based on particle codes which address these physical regimes. The Principal Investigator, Dr. Michael Combi, has over 20 years of experience in the development of particle-kinetic and hybrid kinetichydrodynamics models and their direct use in data analysis. He has also worked in ground-based and space-based remote observational work and on spacecraft instrument teams. His research has involved studies of cometary atmospheres and ionospheres and their interaction with the solar wind, the neutral gas clouds escaping from Jupiter s moon Io, the interaction of the atmospheres/ionospheres of Io and Europa with Jupiter s corotating magnetosphere, as well as Earth s ionosphere. This report describes our progress during the year. The contained in section 2 of this report will serve as the basis of a paper describing the method and its application to the cometary coma that will be continued under a research and analysis grant that supports various applications of theoretical comet models to understanding the inner comae of comets (grant NAGS- 13239 from the Planetary Atmospheres program).

Combi, Michael R.↗
DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

An Open-Source Framework for the Computational Analysis and Design of Autothermal Chemical Processes

The main project object was to develop software tools for simulating non-equilibrium autothermal processes, improving the prospects for identifying and designing such systems. The project demonstrates the use of these tools to simulate autothermal pyrolysis, a process recently developed at the pilot-plant scale at Iowa State University. In such process, instead of externally heating a reactor to pro-vide the enthalpy of pyrolysis, sufficient oxygen in the form of air is introduced into the reactor to support partial oxidation of reactants and products with the exothermic energy released supporting endothermic pyrolysis reactions. A fluidized bed is used to assure good mixing of biomass and oxidant and provide an isothermal reaction environment. The amount of oxygen required depends upon the kind of biomass being pyrolyzed and parasitic heat losses from the reactor. For example, for woody biomass pyrolyzed under conditions that simulate adiabatic operation, equivalence ratios can be as low as 0.06, compared to 0.20 or higher for autothermal gasifiers. By removing the heat transfer bottleneck of conventional pyrolysis, operation in autothermal mode allowed a significant increase in reactor throughput process, approaching five times the throughput of the conventionally operated pyrolyzer. Different simulation strategies were considered and developed: a zero-dimensional chemistry model was used to verify the applicability of kinetic schemes to predict biomass fast pyrolysis in autothermal conditions. Conventional chemical reactor models such as the plug flow reactor and the partially stirred reactor were used to investigate the role of mixing in the fluidized bed pyrolyzer and to establish the impact of mixing time on the gas-phase reactions. A comprehensive multiphase computational fluid dynamics (mCFD) framework, including polydisperse granular phase modeling and detailed chemical kinetics was formulated and used to model the experimental setup for autothermal biomass fast pyrolysis at ISU. Multiphase CFD was also used to investigate the role of biomass feed positioning on the mixing of biomass in the pyrolizer. Finally, a reduced order model (ROM), suitable to be implemented in process simulators was obtained. Both the mCFD and the ROM were validated against experiments.

09 BIOMASS FUELS↗
DOE OSTI2024

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

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

Investigation of post-injection strategies for diesel engine Catalyst Heating Operation using a vapor-liquid-equilibrium-based spray model

Most multidimensional engine simulations spend much time solving for non-equilibrium spray dynamics (atomization, collision, vaporization). However, their accuracy is limited by significant grid dependency, and the need for extensive calibration. This is critical for modeling cold-start diesel fuel post injections, which occur at low temperatures and pressures, far from typical model validation ranges. At the same time, resolving micron-scale spray phenomena would render full Eulerian multiphase calculations prohibitive. In this study, an improved phase equilibrium based approach was implemented and assessed for simulating diesel catalyst heating operation strategies. A phase equilibrium solver based on the model by Yue and Reitz [1] was implemented: a fully multiphase CFD solver is employed with an engineering-size engine grid, and fuel injection is modeled using the standard Lagrangian parcels approach. Mass and energy from the liquid parcels are released to the Eulerian multiphase mixture according to an equilibrium-based liquid jet model. An improved phase equilibrium solver was developed to handle large real-gas mixtures such as those from accurate chemical kinetics mechanisms. The liquid-jet model was improved such that momentum transfer to the Eulerian solver better reproduces the physical spray jet structure. Validation of liquid/vapor penetration predictions showed that the model yields accurate results with very limited tuning and low sensitivity to the few calibration constants. In-cylinder simulations of diesel catalyst heating operation strategies showed that capturing spray structure is paramount when short, transient injection pulses and low temperatures are present. Furthermore, the EP model provides improved predictions of post-injection spray structure and ignitability, while conventional spray modeling does not capture the increase of liquid penetration during the expansion stroke. Finally, the only important EP model calibration constant, Cliq, does not affect momentum transfer, but it changes the local charge cooling distribution through the local energy transfer, which makes it candidate to additional research. The results confirm that non-equilibrium spray processes do not need to be resolved in engineering simulations of high-pressure diesel sprays.

42 ENGINEERING↗
DOE OSTI2025

Laser Measurements and Modeling of Shock Tunnel Freestream Velocity and Multispecies Thermal Nonequilibrium

Coherent anti-Stokes Raman scattering (CARS) and nitric oxide molecular tagging velocimetry (NO-MTV) were used to characterize the freestream in Sandia’s Hypersonic Shock Tunnel (HST) using a burst-mode laser operated at 100 kHz. Experiments were performed at nominal freestream velocities of 3 and 4 km/s using both air and N 2 test gas. The CARS diagnostic provides nonequilibrium characterization of the flow by measuring vibrational and rotational temperatures of N 2 and O 2 , which are compared to NO temperatures from separate laser absorption experiments. Simultaneous, collinear freestream velocities were measured using NO MTV along with pitot pressure measurements. Furthermore, this extensive freestream dataset is compared to nonequilibrium CFD that is capable of modeling species-specific, vibrational temperatures throughout the nozzle expansion. Significant nonequilibrium between vibrational and rotational temperatures was measured at each flow condition. N 2 exhibits the most nonequilibrium, followed by O 2 and NO. The CFD model captures this trend, although it consistently overpredicts N 2 and O 2 vibrational temperatures. At 3 km/s, the modeled NO nonequilibrium is underpredicted, whereas it is overpredicted at 4 km/s. Good agreement is seen between CFD and the velocity and rotational temperature measurements. Experiments with water added to the test gas yielded no discernable difference in vibrational relaxation.

Aerodynamics↗
DOE OSTI2022

Modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge

In this work, we present the results of two-dimensional modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge. A multiscale adaptive reduced chemistry solver for plasma assisted combustion (MARCS-PAC) based on PASSKEy discharge modeling package and compressible multi-component reactive flow solver ASURF+ is developed and validated. This model is applied to simulate the impact of non-equilibrium plasma excitation and electrode geometry and heat loss on the dynamics of the discharge from streamer to spark and ignition kernel development in a H 2 /air mixture with a pair of cylindrical electrodes. The results show that the plasmagenerated species (N 2 (A), N 2 (B), N 2 (a'), N 2 (C), O( 1 D), O and H) in the spark and afterglow significantly accelerate the ignition kernel development. The increase of discharge voltage at the same total discharge energy promotes the non-equilibrium active species production. It is found that the production of electronically excited species at higher reduced electric field strength is more efficient in enhancing ignition in comparison to the vibrational excitation and heating. Moreover, the 2D simulation clearly reveals that the electric field and active species distribution are highly non-uniform. The streamers are initiated at the sharp outer edges of the negative and positive electrodes by a strong electric field while the electric field is much weaker at the centerline of the electrodes. Furthermore, the simulations reveal that the ignition enhancement is sensitive to the variation of electrode shape, diameter, and gap size due to the changes of electric field distribution and location of streamer formation. A cylindrical electrode produces a larger discharge volume and ignition kernel than the parabolic and spherical electrodes, when the discharge is localized near the axis of the gap. It is found that there is a non-monotonic dependence of ignition kernel size on the electrode diameter and inter-electrode distance. The increase of electrode diameter and gap size above the optimal conditions leads to the reduction of ignition kernel volume, due to the decrease of active species concentration and gas temperature. At a larger electrode surface area and electrode diameter as well as smaller electrode gap size, the heat loss to electrode plays a greater role in reducing the ignition kernel size and slowing ignition kernel development. This work provides insights and guidance to understand the kinetic enhancement of non-equilibrium plasma and the effects of electrode geometries on ignition for the optimization ignitors in advanced engines.

42 ENGINEERING↗
DOE OSTI2024

1D Simulation of Avalanche to Streamer to Spark Transition of Plasma Discharge in Ammonia-Air Combustion

The backward problem of plasma assisted combustion emphasizes evaluating the effect of the evolving thermochemical state on the plasma discharge. This paper investigates the dependence of avalanche to streamer to spark formation dynamics and kinetics on the gas composition and temperature at different points in an ammonia-air premixed laminar flame using a self-consistent multigrid-based 1D plasma solver. Different values of ..alpha.., the coefficient for effective ionization events per unit length, have been reported for electron avalanches in air and stoichiometric NH3-air mixtures. The streamer inception has been shown to obey the Meek's criterion. An exponential reduction in streamer and spark formation time has been observed from plasma simulations at different points in the unburnt, pre-heat zone, reaction zone and the fully burnt regions of the premixed flame. While the enhancement of the reduced electric field with increasing temperature affects effective ionization, there exists a minimum breakdown field for streamer formation, which does not vary proportionally with the changing number density of the gas. The change in the mixture from reactants (NH3, O2, N2) to products of complete combustion of ammonia in air (N2, H2O) has also been shown to affect the streamer and spark formation. Finally, the major pathways during the streamer and spark phases which are responsible for producing important radicals used in combustion of NH3 are also discussed.

ADVANCED PROPULSION SYSTEMS,INORGANIC, ORGANIC, PH↗
DOE OSTI2025

Effects of non-equilibrium ionization and excitation on radiation absorption in plasma plumes induced by ablation of metal targets with nanosecond laser pulses

Ionization and radiation absorption in nanosecond laser-induced plasma plumes are routinely modeled using the Saha–Boltzmann equilibrium ionization model (EQM). However, the equilibrium assumption can be inaccurate during the laser pulse when non-equilibrium effects significantly impact radiation absorption. In the present work, the EQM and non-equilibrium collisional-radiative model (CRM) are compared to reveal the effect of plasma non-equilibrium on radiation absorption in non-homogeneous plumes and degree of plasma shielding. Simulations of plume expansion induced by irradiation of a copper target in 1 atm argon background gas with a 10 ns Gaussian pulse at a fluence from 8 Jcm −2 to 14 Jcm −2 are performed with a hybrid computational model that couples a lumped particle direct simulation Monte Carlo method with either CRM or EQM. The simulations show that the EQM strongly underestimates the effects of ionization and radiation absorption compared to CRM and, contrary to the CRM, predict strong ionization of the background gas. The differences between the models are explained by the qualitatively different coupling between plume expansion and dynamics of ionization and exitation processes in the CRM and EQM under conditions when the characteristic times for most radiation- and electron-induced processes are longer than the pulse duration. The CRM-based predictions are also found to agree much better with available experimental data. In conclusion, these results indicate that the model of Saha–Boltzmann equilibrium cannot be used for reliable prediction of the degree of plasma shielding in plumes induced by nanosecond laser pulses or for processing results of spectroscopic measurements at early stages of expansion of such plumes.

Physics - Plasma physics↗
DOE OSTI2022

The effect of elastic anisotropy on the symmetry selection of irradiation-induced void superlattices in cubic metals

We report self-organized microstructures and patterns have been widely observed in non-equilibrium physical systems. In particular, irradiation in metals creates far-from-equilibrium environments, in which the competing dynamics of defect production and annihilation can lead to unique self-organized superlattice structures, e.g., void and gas bubble superlattices. From a physical point of view, the superlattice structures are dictated by the intrinsic symmetry breaking in the metals, i.e., anisotropy caused by the breaking of continuous rotational symmetry. In the literature, two distinctive anisotropies, elastic anisotropy and diffusion anisotropy of interstitials, have been proposed to be the origins of superlattice formation. However, it is still unclear which anisotropy dominates the symmetry selection of superlattice structures. In this paper, we study elastic anisotropy and its effect on the symmetry of void superlattices. By using theoretical analyses and phase field simulations, we show that elastic anisotropy in cubic metals can lead to either face-centered cubic or simple cubic superlattices depending on the Zener anisotropy ratio. The superlattices formed under this elastic anisotropy mechanism must form under the influence of spinodal decomposition, as the mechanism requires perturbations in the vacancy concentration field to develop into spatially-static concentration waves. We compare to existing work on symmetry selection in superlattices via diffusion anisotropy and to experimental observations, and we suggest that concentration wave development under the influence of elastic anisotropy is not the mechanism for symmetry selection during the formation of irradiation-induced void superlattices, but that diffusion anisotropy could be the dominant mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

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

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.↗
NASA NTRS1996

Nonlinear Dynamics and Nucleation Kinetics in Near-Critical Liquids

The objective of our study is to model the nonlinear behavior of a near-critical liquid following a rapid change of the temperature and/or other thermodynamic parameters (pressure, external electric or gravitational field). The thermodynamic critical point is manifested by large, strongly correlated fluctuations of the order parameter (particle density in liquid-gas systems, concentration in binary solutions) in the critical range of scales. The largest critical length scale is the correlation radius r(sub c). According to the scaling theory, r(sub c) increases as r(sub c) = r(sub 0)epsilon(exp -alpha) when the nondimensional distance epsilon = (T - T(sub c))/T(sub c) to the critical point decreases. The normal gravity alters the nature of correlated long-range fluctuations when one reaches epsilon approximately equal to 10(exp -5), and correspondingly the relaxation time, tau(r(sub c)), is approximately equal to 10(exp -3) seconds; this time is short when compared to the typical experimental time. Close to the critical point, a rapid, relatively small temperature change may perturb the thermodynamic equilibrium on many scales. The critical fluctuations have a hierarchical structure, and the relaxation involves many length and time scales. Above the critical point, in the one-phase region, we consider the relaxation of the liquid following a sudden temperature change that simultaneously violates the equilibrium on many scales. Below T(sub c), a non-equilibrium state may include a distribution of small scale phase droplets; we consider the relaxation of such a droplet following a temperature change that has made the phase of the matrix stable.

Patashinski, Alexander Z.↗