Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non-equilibrium processing”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Solidification of eutectic system alloys in space (M-19)

It is well known that in the liquid state eutectic alloys are theoretically homogeneous under 1 g conditions. However, the homogeneous solidified structure of this alloy is not obtained because thermal convection and non-equilibrium solidification occur. The present investigators have clarified the solidification mechanisms of the eutectic system alloys under 1 g conditions by using the in situ observation method; in particular, the primary crystals of the eutectic system alloys never nucleated in the liquid, but instead did so on the mold wall, and the crystals separated from the mold wall by fluid motion caused by thermal convection. They also found that the equiaxed eutectic grains (eutectic cells) are formed on the primary crystals. In this case, the leading phase of the eutectic must agree with the phase of the primary crystals. In space, no thermal convection occurs so that primary crystals should not move from the mold wall and should not appear inside the solidified structure. Therefore no equiaxed eutectic grains will be formed under microgravity conditions. Past space experiments concerning eutectic alloys were classified into two types of experiments: one with respect to the solidification mechanisms of the eutectic alloys and the other to the unidirectional solidification of this alloy. The former type of experiment has the problem that the solidified structures between microgravity and 1 g conditions show little difference. This is why the flight samples were prepared by the ordinary cast techniques on Earth. Therefore it is impossible to ascertain whether or not the nucleation and growth of primary crystals in the melt occur and if primary crystals influence the formation of the equiaxed eutectic grains. In this experiment, hypo- and hyper-eutectic aluminum copper alloys which are near eutectic point are used. The chemical compositions of the samples are Al-32.4mass%Cu (Hypo-eutectic) and Al-33.5mass%Cu (hyper-eutectic). Long rods for the samples are cast by the Ohno Continuous Casting Process and they show the unidirectionally solidified structure. Each flight and ground sample was made of these same rods. The dimensions of all samples are 4.5 mm in diameter and 23.5 mm in length. Each sample is put in a graphite capsule and then vacuum sealed in a double silica ampoule. Then the ampoule is put in the tantalum cartridge and sealed by electron beam welding. For onbard experiments, a Continuous Heating Furnance (CHF) will be used for melting and solidifying samples under microgravity conditions. Six flight samples will be used. Four samples are hypo-eutectic and two are hyper-eutectic alloys. The surface of the two hypo-eutectic alloy samples are covered with aluminum oxide film to prevent Marangoni convection expected under microgravity conditions. Each sample will be heated to 700 C and held at that temperature for 5 min. After that the samples will be allowed to cool to 500 C in the furnace and they will be taken out of the furnace for He gas cooling. The heating and cooling diagrams for the flight experiments are shown. After collecting the flight samples, the solidified structures of the samples will be examined and the mechanisms of eutectic solidification under microgravity conditions will be determined. It is likely that successful flight experiment results will lead to production of high quality eutectic alloys and eutectic composite materials in space.

Ohno, Atsumi↗
DOE OSTI2025

REBOUND: Reverse Engineering Bidirectional Outflow Under Non-Equilibrium Diffusion

Rare-earth elements (REEs) are essential for electronics, renewable energy, and defense technologies. However, the current supply of REEs relies on mining concentrated in a few countries and energy-intensive separations. DOE’s Basic Energy Sciences (BES) program has launched a grand challenge which aims to ensure a sustainable supply of critical REEs by developing innovative and environmentally friendly separation methods. As an alternative to costly and harmful traditional methods, the Non-Equilibrium Transport Driven Separations (NETS) initiative has created a microfluidic Y-channel co-flow method that applies external fields to exploit magneto- and electrohydrodynamic effects for separating dilute REE ions from complex feedstocks. Computational fluid dynamics (CFD) studies have identified a few operating conditions with promising ion selectivity and separation efficiency. However, challenges remain regarding Y-channel versatility across feedstocks and accurate incorporation of physical phenomena into CFD models. In this work, we develop a multi-fidelity modelling approach which integrates experimental results with CFD simulation to build a surrogate model for the dependence of separation efficiency to variation of design parameters. The surrogate model enables a reinforcement learning (RL) method to adaptively launch CFD and experimental runs, improving model fidelity around optimal Y-channel parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Preliminary Kinetic Analysis of Non-Equilibrium Plasma- Assisted Methanol Pyrolysis and Oxidation Experiments

Efforts to enhance power generation efficiency and reduce emissions have driven interest in novel combustion techniques, including non-equilibrium plasma (NEP) ignitors. NEP ignitors show promise in improving energy conversion efficiency, fuel reforming, emission control, and lean-flammability limits. However, their adoption is hindered by a limited understanding of the interplay between plasma-enhanced combustion and thermal chemistry, particularly for complex fuels under engine-relevant conditions. Developing experimentally validated kinetic mechanisms is therefore critical. Additionally, the increasing interest in renewable biofuels like ethanol and methanol, coupled with the desirable qualities of NEP ignitors, presents a compelling opportunity for study. Therefore, this work acts as an extension of a previous work (Bopaiah et al., 2023) pertaining to the experimental results of NEP-assisted methanol pyrolysis and oxidation. Experiments were performed with a custom-built plasma flow reactor at 0.5 atm and temperatures from 523-1203 K. All reactive mixtures are extremely diluted to minimize exothermicity due to reactivity, allowing isothermal assumptions and the isolation of plasma chemistry from thermal chemistry. A dielectric barrier discharge plasma, at 14 kV and 15 ns full-width half maximum, was applied to the reactive mixture at varying frequencies to maintain the number of pulses with increasing temperature. Steady-state product speciation was performed downstream of the reactor with ex-situ GC/MS diagnostics. The attained experimental results were examined through in-depth analysis performed by means of an in-development plasma-coupled kinetic mechanism. As discussed in the previous work, the plasma significantly accelerates methanol pyrolysis, increasing stable intermediate production, including oxygenated and nitrile species. Plasma-assisted oxidation shows even faster fuel consumption compared to pyrolysis and a 200 K ignition shift compared to thermal oxidation. For plasma-assisted pyrolysis, the model demonstrates that accelerated fuel consumption stems from dissociative quenching of excited N2 states with fuel and H2, generating H radicals that react to rapidly form CH3 and CH2OH radicals. At low temperatures, these radicals recombine to produce oxygenates, while CH3 drives nitrile and hydrocarbon formation at higher temperatures. While the model captures pyrolysis trends well, discrepancies in methane, ethylene, and ethanol predictions are present. Similarly, the model faces challenges in accurately representing plasma-assisted oxidation, predicting a much steeper fuel gradient and ignition 100 K earlier than the experiment. While a similar scheme to pyrolysis is nested in the reaction pathway, the enhancement of the O and H radical fluxes and their initiation of the OH and HO2 radical pools dominate fuel and intermediate oxidation. The overestimation of these processes is shown to be responsible for the divergence of model from experiment. While the modelling predictions of this preliminary mechanism are not perfect, they serve as a valuable starting point. Primarily, they elicited new reaction pathways that are not otherwise possible in thermal chemistry induced reaction kinetics. The results also provide a basis for the future work that should be performed. For example, theoretical and experimental studies on excited nitrogen species and fuel/fuel radical interactions, quantification of the NOx production, and the kinetics behind the slow ignition observed in oxidation should be emphasized.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Quantum simulation of charge and exciton transfer in multi-mode models using engineered reservoirs

Quantum simulation enables studies of open-system dynamics in non-perturbative regimes by programming electronic, vibrational, and environmental interactions on comparable energy scales. Trapped ions offer this capability, combining spins, phonons, and tunable dissipation on one platform. We demonstrate an open-system quantum simulation of charge and exciton transfer in a multi-mode linear vibronic coupling model. Using tailored spin-phonon interactions with reservoir engineering, we emulate a system with two dissipative vibrational modes coupled to donor and acceptor sites and track its non-equilibrium dynamics. We continuously tune the system from the charge transfer regime to the vibrationally assisted exciton transfer regime and find that degenerate modes enhance transfer rates at large energy gaps, while non-degenerate modes activate pathways that reduce the energy-gap dependence. Thus, the presence of one additional vibration introduces interfering pathways and reshapes non-perturbative excitation transfer. Our results establish a scalable, hardware-efficient route to simulate vibronic processes with engineered environments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Foundations of magnetohydrodynamics

In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Intense cavitation-assisted electric discharge as a promising tool for water treatment

This study investigates interrelations between one-electrode Cavitation-Assisted Electric Discharge (CAED), two-electrode CAED, and recently discovered Intense CAED (I-CAED). The one-electrode CAED is a self-triggered nanosecond discharge with pulse energy in the micro-Joule range, which can be generated even by a DC voltage. I-CAED consists of a non-equilibrium part within a low-pressure cavitating region and a micro-spark traversing a liquid film. We hypothesize that CAED propagates from the high-voltage electrode as an ionization wave through bubbles of saturated vapor. Subsequently, the streamer-like discharges in the bubbles may form a continuous plasma channel. Inside the cavitating region, the plasma is strongly non-equilibrium, providing an ideal environment for generating chemically unstable species such as hydrogen peroxide (H 2 O 2 ). Plasma of I-CAED spark is characterized by high electron density and near-thermal equilibrium, emitting a continuous ultraviolet spectrum. The combination of these different discharge parts makes I-CAED in water a highly effective tool for the Advanced Oxidation Process, particularly in water disinfection. Experimentally demonstrated Electric Energy per Order value for disinfection of E. coli-contaminated water is as low as 0.135 ± 0.035 kWh/m 3 /order. Estimates show that the implementation of “dry electrodes” configuration reduces the erosion rate of the electrode material by at least one order of magnitude. Spectral analysis reveals that the continuum emission generated by I-CAED in proximity to metal electrodes deviates from the spectra of discharges spatially decoupled from the electrodes. We assume that this spectral divergence is attributable to blackbody-like emission originating from metallic nanoparticles form during the electrode's erosion process.

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

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗
DOE OSTI2025

Applying deep learning methods to develop new models of molecular charge transfer, nonadiabatic dynamics, and nonlinear spectroscopy in the condensed phase

Photon- and field-induced charge transfer has central importance in the generation and storage of electricity, the novel properties of materials, photo-induced catalysis, and electro-optic activity (e.g., photovoltaic cells, fuel cells, and organic chromophores for use in optical fibers and light-emission diodes). These non-equilibrium electronic and chemical transformations are probed by ultrafast, nonlinear spectroscopies. Accurate simulations play a crucial role in our ability to understand, optimize, and control these transformations. This project applies modern deep learning and machine learning (ML) methods to dramatically improve models of electronic dynamics, electronic-nuclear dynamics, and spectroscopic measurements for improved simulations of chemistry in complex environments, far from equilibrium phenomena, and processes in extreme environments, such as materials exposed to strong or resonant fields. This project develops accurate neural net models that go beyond predictive capability to also provide new insight into the fundamental physics underlying electron and nuclear dynamics. To achieve its objectives, this project explores and develops customized versions of high-capacity deep learning algorithms/models. These techniques are developed with an emphasis on fundamental chemical insight, not just predictive accuracy, to assist the development of the next generation of quantum simulation methods.

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

Understanding Mechanisms of Selective Crystal Nucleation and Growth Using In Situ Electron Microscopy

Separation processes are essential for purifying and isolating critical materials, but are challenging due to the natural tendency of systems to mix. This project explores far-from-equilibrium separation methods using external electric and magnetic fields to drive selective nucleation and crystal growth. Using in situ electrochemical transmission electron microscopy (TEM), we observed ion transport and crystallization behavior under applied electric fields. We found that ions migrate toward the working electrode, where nucleation preferentially occurs. These findings advance our understanding of non-equilibrium separation mechanisms and support the development of efficient, sustainable methods for extracting critical materials from unconventional domestic sources.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Technical note: A modified formulation of dynamic energy budget theory for faster computation of biological growth

Abstract. The mass conservation equation in the presence of boundary fluxes and chemical reactions from non-equilibrium thermodynamics is used to derive a modified dynamic energy budget (mDEB) model. Compared to the standard dynamic energy budget (sDEB) model (Kooijman, 2009), this modified formulation does not place the dilution effect in the mobilization kinetics of reserve biomass, and it maintains the partition principle for reserve mobilization dynamics for both linear and non-linear kinetics. Overall, the mDEB model shares most features with the sDEB model. However, for biological growth that requires multiple nutrients, the mDEB model is computationally much more efficient by not requiring numerical iterations for obtaining the specific growth rate. In an example of modeling the growth of Thalassiosira weissflogii in a nitrogen-limiting chemostat, the mDEB model was found to have almost the same accuracy as the sDEB model while requiring almost half of the computing time of the sDEB model. Since the sDEB model has been successfully applied in numerous studies, we believe that the mDEB model can help improve the modeling of biological growth and the associated ecosystem processes in various contexts.

Tang, Jinyun↗
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 NTRS1992

Evolution of attached and detached slabs and their associated mantle dynamics

Over the two years of the NASA grant, this project has produced a significant amount of research results related to the plate subduction process and the surface crustal deformation at convergent boundaries (i.e., above subduction zones). While some research objectives are completely accomplished, other research tasks remain active and continue to be investigated at present. A steady state analytic thermal model for subducting slabs was used to examine the torques acting on a descending slab. It is found that gravitational torque vanishes when a slab is dipping either vertically or horizontally, unlike previous studies indicating that the magnitude of gravitational torque decreases as dip angle increases. Subsequently, a new time-dependent, analytic thermal model for a subducting slab was developed. The new model enables us to study transient phenomena associated with plate subduction analytically. On the basis of this model, the nature of slab dip angles was evaluated. Slab dip angles are found to be transient features. As they penetrate into the mantle and increase their lengths, the associated gravitational torque also increases resulting in a downward pulling of the slab to the steeper dip angle. This is especially true once a slab penetrates the olivine-spinel phase boundary at about 400 km depth. However, if the phase transformation does not follow the equilibrium condition, the gravitational torque may have a different behavior. This problem was investigated. Except for fast descending slabs, non-equilibrium phase transformation can only slow down the transient increase of slab dip angles discussed earlier. Its effect is not sufficiently strong to reverse the downward pulling for most of the slabs. However, when slabs subducting at 10 cm/yr or faster, a sufficient amount of metastable olivine can exist beneath 400 km. Because of its low density compared with the surrounding spinel, an upward buoyancy is produced resulting in an upward bending of the slab and possibly an upward rotation of the slab such that smaller dip angles are formed. Seismic studies of the Japanese Slab seem to support this interpretation. The development of oroclinal geometries at convergent boundaries was also examined to study plate obduction which is an important ingredient to the initiation of plate subduction. Although the study suggests that surface features are better modeled by block models, the large scale deformation can be adequately studied by viscous models. Such a model is now under development to complete our original objective to study the initiation of plate subduction. Finally, a three-dimensional, finite element, spherical convective model is developed to study dynamic plate subductions. The model development is now complete and it is being tested to ensure its proper operation. The model is able to generate convection results with a viscosity contrast of about 100. Our research continues to push the viscosity contrast to a level that is appropriate for a subducting slab.

Hsui, Albert T.↗
NASA NTRS1992

Pulsatile instability in rapid directional solidification: Strongly-nonlinear analysis

In models of rapid directional solidification, non-equilibrium interfacial conditions are employed. As a result, there is an oscillatory mode of instability, as well as the steady cellular mode, found in the equilibrium model of Mullins and Sekerka. When the temperature field is decoupled from the problem, the preferred wave number for the oscillatory mode is zero, and the interface pulsates in time while remaining spatially uniform. Results from multiple scale analyses in the two limiting cases of the parameters are reported. In these limits, it is found that the instability is a bifurcation to relaxation oscillations; these nonlinear oscillations may be related to the observed microstructure that results from rapid solidification processes such as laser surface remelting.

Braun, Richard J.↗
NASA NTRS1993

Study and modeling of finite rate chemistry effects in turbulent non-premixed flames

The development of numerical models that reflect some of the most important features of turbulent reacting flows requires information about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between turbulent and chemical processes is so strong that it is extremely difficult to isolate the role played by one individual physical phenomenon. Direct numerical simulation (hereafter DNS) allows us to study in detail the turbulence-chemistry interaction in some restricted but completely defined situations. Globally, non-premixed flames are controlled by two limiting regimes: the fast chemistry case, where the turbulent flame can be pictured as a random distribution of local chemical equilibrium problems; and the slow chemistry case, where the chemistry integrates in time the turbulent fluctuations. The Damkoehler number, ratio of a mechanical time scale to chemical time scale, is used to distinguish between these regimes. Today most of the industrial computer codes are able to perform predictions in the hypothesis of local equilibrium chemistry using a presumed shape for the probability density function (pdt) of the conserved scalar. However, the finite rate chemistry situation is of great interest because industrial burners usually generate regimes in which, at some points, the flame is undergoing local extinction or at least non-equilibrium situations. Moreover, this variety of situations strongly influences the production of pollutants. To quantify finite rate chemistry effect, the interaction between a non-premixed flame and a free decaying turbulence is studied using DNS. The attention is focused on the dynamic of extinction, and an attempt is made to quantify the effect of the reaction on the small scale mixing process. The unequal diffusivity effect is also addressed. Finally, a simple turbulent combustion model based on the DNS observations and tractable in real flow configurations is proposed.

Vervisch, Luc↗
NASA NTRS1995

GRO: Black hole models for gamma-ray bursts

The Burst and Transient Source Experiment (BATSE) on board the Compton Gamma Ray Observatory (CGRO) has established that the distribution of gamma-ray bursts (GRB's) is isotropic but is bound radially. This finding suggests that the bursts are either cosmological or they originate from an extended Galactic halo. The implied luminosities and the observed variability of the GRB's on time scales as short as one millisecond suggest that they originate from compact objects. We are presently studying black hole models for GRB's. Any such model must produce a non-thermal photon spectrum to agree with the observed properties. For a wide range of burst parameters the assumed bursting source consists of a non-thermal electron-positron-photon plasma of very high density. It seems possible to produce such a plasma in accretion onto black holes. In our on-going work, we are developing the kinetic theory for a non-equilibrium pair plasma. The main new features of our work are as follows: (1) We do not assume the presence of a thermal electron bath. (2) Non-thermal, high-energy pairs are allowed to have an arbitrary concentration and energy distribution. (3) There is no soft photon source in our model; initially all the photons in the plasma are either energetic X-rays or gamma-rays. (4) The initial energy distribution of the pairs as well as photons is arbitrary. (5) We collect the analytical expressions for the kinetic kernels for all relevant processes. And (6) we present a different approach to finding the time-evolution of pair and photon spectra, which is a combination of the kinetic-theory and the non-linear Monte-Carlo schemes. We have developed many Monte-Carlo programs to model various process, to take into account the time evolution, and to incorporate various physical effects which are unique to non-thermal plasmas. The hydrodynamics of fireballs in GRB's was studied before. Applying results from kinetic theory will improve our understanding of these systems.

Ruderman, Malvin↗
NASA NTRS1996

Dynamical Evolution of a Coronal Streamer-Flux Rope System: 2. A Self-Consistent Non-Planar Magnetohydrodynamic Simulation

The dynamical response of a helmet streamer to a flux rope escape from the sub-photosphere is examined in a physically self-consistent manner within the approximation of axisymmetric three-dimensional magnetohydrodynamics (i.e., so-called '2 1/2 D'). In contrast to the previous planar analyses of Paper 1 (Wu, Guo, and Wang), the present study shows, with the inclusion of out-of-plane components of magnetic and velocity fields, that the magnetic configuration represents a helical flux rope instead of a planar bubble as shown in Paper 1. Because of this more physically-realistic configuration, we are able to examine the dynamical evolution of the helical flux rope's interaction with the helmet streamer. This process leads to the formation of two parts of the solar mass ejection: (i) the expulsion of the helmet dome due to eruption of this flux rope, and (ii) the flux rope's eruption itself. When this two-part feature propagates out to the interplanetary space, it exhibits all the physical characteristics of observed interplanetary magnetic clouds. These numerical simulations also show that the dynamical behavior of the streamer-flux rope system has three distinct states: (i) quasi-equilibrium, (ii) non-equilibrium, and (iii) eruptive state depending on the energy level of the flux rope.

Wu, S. T.↗
NASA NTRS1998

Non-Equilibrium Turbulence Modeling for High Lift Aerodynamics

This phase is discussed in ('Non linear kappa - epsilon - upsilon(sup 2) modeling with application to high lift', Application of the kappa - epsilon -upsilon(sup 2) model to multi-component airfoils'). Further results are presented in 'Non-linear upsilon(sup 2) - f modeling with application to high-lift' The ADI solution method in the initial implementation was very slow to converge on multi-zone chimera meshes. I modified the INS implementation to use GMRES. This provided improved convergence and less need for user intervention in the solution process. There were some difficulties with implementation into the NASA compressible codes, due to their use of approximate factorization. The Helmholtz equation for f is not an evolution equation, so it is not of the form assumed by the approximate factorization method. Although The Kalitzin implementation involved a new solution algorithm ('An implementation of the upsilon(sup 2) - f model with application to transonic flows'). The algorithm involves introducing a relaxation term in the f-equation so that it can be factored. The factorization can be into a plane and a line, with GMRES used in the plane. The NASA code already evaluated coefficients in planes, so no additional memory is required except that associated the the GMRES algorithm. So the scope of this project has expanded via these interactions. . The high-lift work has dovetailed into turbine applications.

Durbin, P. A.↗
NASA NTRS2001

Quantitative Species Measurements in Microgravity Combustion Flames

Flame-vortex interactions are canonical configurations that can be used to study the underlying processes occurring in complicated turbulent reacting flows. The elegant simplicity of the flame-vortex interaction permits the study of these complex interactions under relatively controllable experimental configurations, in contrast to direct measurements in turbulent flames. The ability to measure and model the fundamental phenomena that occur in a turbulent flame, but with time and spatial scales which are amenable to our diagnostics, permits significant improvements in the understanding of turbulent combustion under both normal and reduced gravity conditions. In this paper, we report absolute mole fraction measurements of methane in a reacting vortex ring. These microgravity experiments are performed in the 2.2-sec drop tower at NASA Glenn Research Center. In collaboration with Drs. Chen and Dahm at the University of Michigan, measured methane absorbances are incorporated into a new model from which the temperature and concentrations of all major gases in the flame can be determined at all positions and times in the development of the vortex ring. This is the first demonstration of the ITAC (Iterative Temperature with Assumed Chemistry) approach, and the results of these computations and analyses are presented in a companion paper by Dahm and Chen at this Workshop. We believe that the ITAC approach will become a powerful tool in understanding a wide variety of combustion flames under both equilibrium and non-equilibrium conditions.

Silver, Joel A.↗