Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonequilibrium chemistry”

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 163 records · Page 9

Tailoring Lifetimes and Properties of Carbodiimide-Fueled Covalently Cross-linked Polymer Networks

The mechanical properties of nonequilibrium polymer hydrogels obtained from the transient cross-linking of polymer chains by a chemical fuel were investigated. Aqueous polymers featuring pendant carboxylic acids were treated with a carbodiimide to give anhydride-cross-linked gels. The anhydrides spontaneously hydrolyze back to the polymer solution, and the cycle can be repeated multiple times. Oscillatory rheology was employed to study the effects of temperature, fuel concentration, chain length, and polymer composition on the storage and loss moduli of the polymeric materials as well as the time taken for the polymers to undergo decross-linking. Regardless of the temperature used, at constant carbodiimide concentration, degelation times are more sensitive to experimental temperature than are the peak storage moduli. Decross-linking times decrease with increasing temperature. As carbodiimide concentration decreases, there is a decrease in moduli and decross-linking times. Within the scope of materials studied, the polymer structure was found to have a relatively small impact on the transient properties of gel networks compared to the fuel concentration and temperature. Furthermore, these findings facilitate the design of tunable on-demand networks and gels.

36 MATERIALS SCIENCE↗

On the molecular correlations that result in field-dependent conductivities in electrolyte solutions

Employing recent advances in response theory and nonequilibrium ensemble reweighting, we study the dynamic and static correlations that give rise to an electric field-dependent ionic conductivity in electrolyte solutions. We consider solutions modeled with both implicit and explicit solvents, with different dielectric properties, and at multiple concentrations. Implicit solvent models at low concentrations and small dielectric constants exhibit strongly field-dependent conductivities. We compare these results to Onsager-Wilson theory of the Wien effect, which provides a qualitatively consistent prediction at low concentrations and high static dielectric constants but is inconsistent away from these regimes. The origin of the discrepancy is found to be increased ion correlations under these conditions. Explicit solvent effects act to suppress nonlinear responses, yielding a weakly field-dependent conductivity over the range of physically realizable field strengths. By decomposing the relevant time correlation functions, we find that the insensitivity of the conductivity to the field results from the persistent frictional forces on the ions from the solvent. Our findings illustrate the utility of nonequilibrium response theory in rationalizing nonlinear transport behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Equilibrium and disequilibrium chemistry of adiabatic, solar-composition planetary atmospheres

The impact of atmospheric and cloud-structure models on the nonequilibrium chemical behavior of the atmospheres of the Jovian planets is discussed. Quantitative constraints on photochemical, lightning, and charged-particle production of organic matter and chromophores are emphasized whenever available. These considerations imply that inorganic chromophore production is far more important than that of organic chromophores, and that lightning is probably a negligibly significant process relative to photochemistry on Jupiter. Production of complex molecules by gas-phase disequilibrium processes on Saturn, Uranus, and Neptune is severely limited by condensation of even simple intermediates.

Lewis, J. S.↗

Spin-orbit-proximitized ferromagnetic metal by monolayer transition metal dichalcogenide: Atlas of spectral functions, spin textures, and spin-orbit torques in Co/MoSe 2 , Co/WSe 2 and Co/TaSe 2 heterostructures

In this work, the heterostructures composed of an ultrathin ferromagnetic metal (FM) and a material hosting strong spin-orbit (SO) coupling are the principal resource for SO torque and spin-to-charge conversion nonequilibrium effects in spintronics. The key quantity in theoretical description of these effects is nonequilibrium spin density, which can appear on any monolayer of the heterostructure through which the current is flowing as long as the monolayer bands are affected by the native or proximity induced SO coupling. Here we demonstrate how hybridization of wave functions of Co layer and a monolayer of transition metal dichalcogenides (TMDs)—such as semiconducting MoSe 2 and WSe 2 or metallic TaSe 2 —can lead to dramatic transmutation of electronic and spin structure of Co within some distance away from its interface with TMD, when compared to the bulk of Co or its surface in contact with a vacuum. This is due to proximity induced SO splitting of Co bands encoded in the spectral functions and spin textures on its monolayers, which we obtain using noncollinear density functional theory (ncDFT) combined with equilibrium Green function (GF) calculations. In fact, SO splitting is present due to structural inversion asymmetry of the bilayer—i.e., just the presence of the interface—even if SO coupling within TMD monolayer is artificially switched off in ncDFT calculations. However, switching it on makes the effects associated with proximity SO coupling within Co layer about five times larger. Injecting spin-unpolarized charge current through SO-proximitized monolayers of Co generates nonequilibrium spin density over them, so that its cross product with the magnetization of Co determines SO torque. The SO torque computed via first-principles quantum transport methodology, which combines ncDFT with nonequilibrium GF calculations, can be used as the screening parameter to identify optimal combination of materials and their interfaces for applications in spintronics. In particular, we identify heterostructure two-monolayer-Co/monolayer-WSe 2 as the most optimal one, at least in the clean limit.

2-dimensional systems↗

Stratification of polymer mixtures in drying droplets: Hydrodynamics and diffusion

In this work, we study the evaporation-induced stratification of a mixture of short and long polymer chains in a drying droplet using molecular simulations. We systematically investigate the effects of hydrodynamic interactions (HI) on this process by comparing hybrid simulations accounting for HI between polymers through the multiparticle collision dynamics technique with free-draining Langevin dynamics simulations neglecting the same. We find that the dried supraparticle morphologies are homogeneous when HI are included but are stratified in core–shell structures (with the short polymers forming the shell) when HI are neglected. The simulation methodology unambiguously attributes this difference to the treatment of the solvent in the two models. We rationalize the presence (or absence) of stratification by measuring phenomenological multicomponent diffusion coefficients for the polymer mixtures. The diffusion coefficients show the importance of not only solvent backflow but also HI between polymers in controlling the dried supraparticle morphology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A code validation strategy and facility for nonequilibrium, reacting flows

The rotational and vibrational population distributions in a hypersonic flow of hot iodine vapor were investigated using a pilot wind tunnel. The tunnel provides run times of about 20 min, with a scale-up to a larger facility appearing to be possible. It will provide a 20-cm-diameter, Mach 9 flow of iodine vapor with only small amounts of free stream nonequilibrium. An analysis of the flow response to a normal shock wave indicates that significant chemistry will occur in model flow fields in the larger facility.

Muntz, E. P.↗

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗

Chemically Fueled Reinforcement of Polymer Hydrogels

Carbodiimide-fueled anhydride bond formation has been used to enhance the mechanical properties of permanently crosslinked polymer networks, giving materials that exhibit transitions from soft gels to covalently reinforced gels, eventually returning to the original soft gels. Temporary changes in mechanical properties result from a transient network of anhydride crosslinks, which eventually dissipate by hydrolysis. Over an order of magnitude increase in the storage modulus is possible through carbodiimide fueling. Here, the time-dependent mechanical properties can be modulated by the concentration of carbodiimide, temperature, and primary chain architecture. Because the materials remain rheological solids, new material functions such as temporally controlled adhesion and rewritable spatial patterns of mechanical properties have been realized.

36 MATERIALS SCIENCE↗

Treating random sequential addition via the replica method

While many physical processes are non-equilibrium in nature, the theory and modeling of such phenomena lag behind theoretical treatments of equilibrium systems. The diversity of powerful theoretical tools available to describe equilibrium systems has inspired strategies that map non-equilibrium systems onto equivalent equilibrium analogs so that interrogation with standard statistical mechanical approaches is possible. In this work, we revisit the mapping from the non-equilibrium random sequential addition process onto an equilibrium multi-component mixture via the replica method, allowing for theoretical predictions of non-equilibrium structural quantities. We validate the above approach by comparing the theoretical predictions to numerical simulations of random sequential addition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Programming Function via Soft Materials (PFvSM) (Final Report)

The PFvSM cluster at the University of Illinois greatly deepened the science of soft materials and expanded the means for embedding within them unprecedented functionalities. We pioneered a mesoscopic approach to materials design that combined leading competencies in fundamental experimental and theoretical soft matter science and fabrication that powerfully establishes equilibrium and nonequilibrium assembly of novel and complex building blocks as a preeminent contributor to progress in materials chemistry. Our cluster’s vision was an integrated program of work centered on soft matter science with the ultimate goal of impacting important energy technologies. As elementary units we employed function-encoded colloids, nanoparticles, large mesh fibrillar gels, deterministically fabricated networks and hybrid objects. These materials were employed as building blocks for diverse model systems that allowed fundamental soft matter principles of wide relevance to be established. Their assembly into materials with useful functional behavior provided new approaches to challenging problems in basic energy science including energy harvesting, transport and storage. To achieve these goals the cluster brought together a team of recognized world-leaders in materials synthesis and processing, in-situ characterization, theory and modeling, and materials design, integration and fabrication

36 MATERIALS SCIENCE↗

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

Preconditioning for the Navier-Stokes equations with finite-rate chemistry

The extension of Van Leer's preconditioning procedure to generalized finite-rate chemistry is discussed. Application to viscous flow is begun with the proper preconditioning matrix for the one-dimensional Navier-Stokes equations. Eigenvalue stiffness is resolved and convergence-rate acceleration is demonstrated over the entire Mach-number range from nearly stagnant flow to hypersonic. Specific benefits are realized at the low and transonic flow speeds typical of complete propulsion-system simulations. The extended preconditioning matrix necessarily accounts for both thermal and chemical nonequilibrium. Numerical analysis reveals the possible theoretical improvements from using a preconditioner for all Mach number regimes. Numerical results confirm the expectations from the numerical analysis. Representative test cases include flows with previously troublesome embedded high-condition-number areas. Van Leer, Lee, and Roe recently developed an optimal, analytic preconditioning technique to reduce eigenvalue stiffness over the full Mach-number range. By multiplying the flux-balance residual with the preconditioning matrix, the acoustic wave speeds are scaled so that all waves propagate at the same rate, an essential property to eliminate inherent eigenvalue stiffness. This session discusses a synthesis of the thermochemical nonequilibrium flux-splitting developed by Grossman and Cinnella and the characteristic wave preconditioning of Van Leer into a powerful tool for implicitly solving two and three-dimensional flows with generalized finite-rate chemistry. For finite-rate chemistry, the state vector of unknowns is variable in length. Therefore, the preconditioning matrix extended to generalized finite-rate chemistry must accommodate a flexible system of moving waves. Fortunately, no new kind of wave appears in the system. The only existing waves are entropy and vorticity waves, which move with the fluid, and acoustic waves, which propagate in Mach number dependent directions. The nonequilibrium vibrational energies and species densities in the unknown state vector act strictly as convective waves. The essential concept for extending the preconditioning to generalized chemistry models is determining the differential variables which symmetrize the flux Jacobians. The extension is then straight-forward. This algorithm research effort will be released in a future version of the production level computational code coined the General Aerodynamic Simulation Program (GASP), developed by Walters, Slack, and McGrory.

Godfrey, Andrew G.↗

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↗

Calculation of self-diffusion coefficients in supercritical carbon dioxide using mean force kinetic theory

This paper presents an application of mean force kinetic theory (MFT) to the calculation of the self-diffusivity of CO 2 in the supercritical fluid regime. Two modifications to the typical application of MFT are employed to allow its application to a system of molecular species. Therefore, the first is the assumption that the inter-particle potential of mean force can be obtained from the molecule center-of-mass pair correlation function, which in the case of CO 2 is the C–C pair correlation function. The second is a new definition of the Enskog factor that describes the effect of correlations at the surface of the collision volume. The new definition retains the physical picture that this quantity represents a local density increase, resulting from particle correlations, relative to that in the zero density homogeneous fluid limit. These calculations are facilitated by the calculation of pair correlation functions from molecular dynamics (MD) simulations using the FEPM2 molecular CO2 model. The self-diffusivity calculated from theory is in good agreement with that from MD simulations up to and slightly beyond the density at the location of the Frenkel line. The calculation is compared with and is found to perform similarly well to other commonly used models but has a greater potential for application to systems of mixed species and to systems of particles with long range interatomic potentials due to electrostatic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparison of vibration dissociation coupling and radiative heat transfer models for AOTV/AFE flowfields

A series of detailed studies comparing various vibration dissociation coupling models, reaction systems and rates, and radiative heating models has been conducted for the nonequilibrium stagnation region of an AFE/AOTV vehicle. Atomic and molecular nonequilibrium radiation correction factors have been developed and applied to various absorption coefficient step models, and a modified vibration dissociation coupling model has been shown to yield good vibration/electronic temperature and concentration profiles. While results indicate sensitivity to the choice of vibration dissociation coupling model and to the nitrogen electron impact ionization rate, by proper combinations accurate flowfield and radiative heating results can be obtained. These results indicate that nonequilibrium effects significantly affect the flowfield and the radiative heat transfer. However, additional work is needed in ionization chemistry and absorption coefficient modeling.

Carlson, Leland A.↗

Multiple Boundary Layer Instability Modes with Nonequilibrium and Wall Temperature Effects Using LASTRAC

Prediction and control of boundary layer transition from laminar to turbulent is important to many flow regimes and vehicle designs, including vehicles operating at hypersonic conditions where nonequilibrium effects may be encountered. Wall cooling is known to affect the instability characteristics of the boundary layer and subsequently the transition location. Design considerations, including material failure and fuel chemistry, require the use of actively cooled walls in hypersonic vehicles, further motivating the study of wall temperature effects on top of the considerations of reducing heat flux, drag, and uncertainty. In this work, we analyze the stability of a boundary layer with chemical and thermal nonequilibrium on a Mach 20, 6 wedge. We investigate the effects of wall temperature on multiple unstable modes individually and on the integrated growth of disturbances along the surface. We use the LAngley Stability and TRansition Analysis Code (LASTRAC) to evaluate boundary layer stability, using capabilities implemented by the authors. Included are results that address chemical nonequilibrium with both thermal equilibrium and nonequilibrium.

Kline, H. L.↗

Entropy production in model colloidal suspensions under shear via the fluctuation theorem

Dissipative systems often exhibit novel and unexpected properties. This is, for instance, the case of simple liquids, which, when subjected to shear and after reaching a steady state, can exhibit a negative entropy production over finite length scales and timescales. This result, among others, is captured by nonequilibrium relations known as fluctuation theorems. Using nonequilibrium molecular dynamics simulations, we examine how, by fine-tuning the properties of the components of a complex fluid, we can steer the nonequilibrium response of the fluid. More specifically, we show how we control the nonequilibrium probability distribution for the shear stress and, in turn, how often states with a negative entropy production can occur. To achieve this, we start by characterizing how the size for the liquid matrix impacts the probability of observing negative entropy states, as well as the timescale over which these can be observed. We then measure how the addition of larger particles to this liquid matrix, i.e., simulating a model colloidal suspension, results in an increase in the occurrence of such states. As a result, this suggests how modifications in the composition of the mixture and in the properties of its components lead to an increase in the probability of observing states of negative entropy production and, thus, for the system to run in reverse.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An analysis of combustion studies in shock expansion tunnels and reflected shock tunnels

The effect of initial nonequilibrium dissociated air constituents on the combustion of hydrogen in high-speed flows for a simulated Mach 17 flight condition was investigated by analyzing the results of comparative combustion experiments performed in a reflected shock tunnel test gas and in a shock expansion tunnel test gas. The results were analyzed and interpreted with a one-dimensional quasi-three-stream combustor code that includes finite rate combustion chemistry. The results of this study indicate that the combustion process is kinetically controlled in the experiments in both tunnels and the presence of the nonequilibrium partially dissociated oxygen in the reflected shock tunnel enhances the combustion. Methods of compensating for the effect of dissociated oxygen are discussed.

Jachimowski, Casimir J.↗