Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion length”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

In-Situ Partial Pressure Measurements and Visual Observation during Crystal Growth of ZnSe by Seeded Physical Vapor Transport

An in-situ monitoring furnace was constructed with side windows to perform partial pressure measurements by optical absorption and visual observation of the growing crystal. A fused silica -rowth ampoule with a 4.5 cm long square tube between the source and the seed was prepared for the optical absorption measurements. A ZnSe crystal was grown by the seeded physical vapor transport (PVT) technique in the horizontal configuration. The growth temperature was 1120 C and the furnace translation rate was 3nmVday. Partial pressures of Se2, P(sub Se2), at three locations along the length of the growth ampoule were measured at 90 min intervals during the growth process. The measured P (sub Se2) were in the range of 2.0 to 6.5 x 10(exp -3) atm. The P(sub Se2) results indicated that the partial pressure profile was inconsistent with the results of the one-dimensional diffusion mass transport model and that the source composition shifted toward Se-rich during the run, i.e. the grown crystal was more Zn-rich than the source. The visual observation showed that the seed crystal first etched back, with greater thermal etching occurring along the edges of the seed crystal. Once the growth started, the crystal crew in a predominately contactless mode and facets were evident during growth. The crystal did not grow symmetrically which is believed to be due to the unintentional asymmetry of the radial thermal profile in the furnace.

Su, Ching-Hua↗
NASA NTRS2013

Self Induced Buoyant Blow Off in Upward Flame Spread on Thin Solid Fuels

Upward flame spread experiments were conducted on a thin fabric cloth consisting of 75% cotton and 25% fiberglass. The sample is sandwiched symmetrically with stainless steel plates with the exposed width varying between 2 to 8.8 cm from test to test and >1.5m tall. The bottom edge was ignited resulting in a symmetric two sided flame. For the narrower samples (. 5cm), two sided flame growth would proceed until reaching some limiting value (15-30 cm depending on sample width). Fluctuation or instability of the flame base on one side would initially become visible and then the flame base would retreat downstream and cause extinguishment on one side. Detailed examination of the still images shows that the fuel continues to vaporize from the extinguished side due to the thermally thin nature of the fuel. But, due to the remaining inert fiberglass mesh, which acts as a flashback arrestor, the extinguished side was not able to be reignited by the remaining flame. The remaining flame would then shrink in length due to the reduced heat transfer to the solid to a shorter length. The one-sided flame will spread stably with a constant speed and a constant flame length to the end of the sample. A constant length flame implies that the pyrolysis front and the burnt out fronts move at the same speed. For the wider samples (. 7cm), no one-sided extinction is observed. Two-sided flames spread all the way to the top of the sample. For these wider widths, the flames are still growing and have not reached their limiting length if it exists. Care was taken to minimize the amount of non-symmetries in the experimental configuration. Repeated tests show that blow-off can occur on either side of the sample. The flame growth is observed to be very symmetric during the growth phase and grew to significant length (>10cm) before extinction of the flame on one side. Our proposed explanation of this unusual phenomenon (i.e. stronger two ]sided flame cannot exist but weaker one-sided flame can) is as follows: The observed one-sided extinction is a blow- off induced by buoyant entrainment. It is known that the flammable diffusion flame regime is bounded by quenching and blow ]off limits when varying incoming air velocity. The narrowest samples tested (between 2 and 5 cm) begin within the flammable range, but as the flame grows, the buoyancy driven air velocity increases at the neighborhood of the flame base. The initially stable flame crosses the extinguishment boundary resulting in a flame blow-off. When one-side of the flame extinguishes, the remaining side shrinks due to the reduced heat transfer to the solid. This reduces the induced velocity and the flame becomes stable. It is proposed that this may have implications to upward flame growth beyond this experiment.

Johnston, Michael C.↗
DOE OSTI2021

The Corrected Distortion model for Lagrangian spray simulation of transcritical fuel injection

In this work, we present a detailed implementation and validation of the droplet modeling framework proposed by Dahms and Oefelein (2016) into the engine commercial CFD software CONVERGE using the User Defined Function (UDF) interface. The model accounts for the nonlinear deformation and oscillation experienced by liquid spray droplet injected into high pressure and temperature. Lagrangian spray simulations of Engine Combustion Network (ECN) Spray A are performed. Model validation against standard experimental measurements of liquid velocity, vapor mixture fraction is conducted. To perform more rigorous model validation, new experimental measurements based on Diffused Back Illumination (DBI) are introduced. The new measurements are processed for Projected Liquid Volume (PLV), which offers as close as possible one-to-one model validation for liquid penetration while offering new insights into the spray physics. Comparison with a One-D model based on adiabatic mixing theory by Siebers (1999) and Desantes et al. (2007) are also conducted. Through these model validation exercises, it is shown that the new framework improves liquid-phase penetration predictions, following a tendency for enhanced evaporation, compared to the standard approach for both Reynolds Average Navier Stokes (RANS) and Large Eddy Simulation (LES). Additionally, at the liquid length, maximum mixture fraction values predicted by the new approach are in good agreement those of an adiabatic mixing model. Qualitative analysis of the spray behaviors during the early stage of the injection process reveals that the proposed framework predicts significant increase in droplet evaporation rate with lower drop drag compared to the current standard approach.

42 ENGINEERING↗
DOE OSTI2024

Tunable Growth of Layered Double Hydroxide Nanosheets through Hydrothermal Conversion of Atomic Layer Deposition Seed Layers

To enable the design and manufacturing of hierarchical nanomaterial architectures, there is a need for synthesis and processing methods that can enable tunable geometric control at the nanoscale while maintaining conformality on complex 3-D templates. Here, in this study, we explore the programmable control of vertically oriented Zn-Al layered double hydroxide (LDH) nanosheet arrays using atomic layer deposition (ALD) to deposit a seed layer of Al 2 O 3 , which is subsequently consumed and converted into the LDH phase under hydrothermal growth conditions. We demonstrate tunable control over the spacing and length of the nanosheets by varying the thickness of the initial ALD seed layer with subnanometer precision. This can be viewed as a nanoscale titration reaction, where Al acts as the limiting reagent during the hydrothermal synthesis of the nanosheets. Elemental mapping demonstrates the dynamic evolution of the resulting morphology, which is driven by surface diffusion and nucleation processes. The conformal nature of ALD allows for hierarchical growth of nanosheets on the surface of a variety of nonplanar substrate geometries, including microposts, paper fibers, and porous ceramic supports. This illustrates the power of ALD to enable bottom-up growth of 3-D nanoarchitectures with tunable geometries by controlling nucleation and growth in subsequent solution reactions.

36 MATERIALS SCIENCE↗
NASA NTRS1973

Thermal expansion method for lining tantalum alloy tubing with tungsten

A differential-thermal expansion method was developed to line T-111 (tantalum - 8 percent tungsten - 2 percent hafnium) tubing with a tungsten diffusion barrier as part of a fuel element fabrication study for a space power nuclear reactor concept. This method uses a steel mandrel, which has a larger thermal expansion than T-111, to force the tungsten against the inside of the T-111 tube. Variables investigated include lining temperature, initial assembly gas size, and tube length. Linear integrity increased with increasing lining temperature and decreasing gap size. The method should have more general applicability where cylinders must be lined with a thin layer of a second material.

Watson, G. K.↗
DOE OSTI2022

Smooth Modified Surfaces of Silicon for the Study of Ionic Liquid Interfaces by Neutron Reflectometry

Silicon surfaces with termination of hydrogen (H–Si), quartz, native silicon (native-Si), and chemically deposited (Au–Si) cd and sputtered gold (Au–Si) s were fabricated with roughness less than 6 Å to study the solid–liquid interface of the ionic liquid (IL) methyltrioctylammonium bis(trifluoromethylsulfonyl)imide, [N 1888 ][TFSI], by neutron reflectometry (NR). The modifications made to the silicon surface were to obtain a native positive or negative charge similar to a charged surface in electrochemical systems. The differences in scattering length density (SLD) based on mass density and composition of IL were examined. Perpendicular to the substrate surface, a distinct transition between the first ion layer and the bulk liquid is resolved out of the reflectivity data for the gold surfaces and in particular the chemically deposited Au–Si. Other interfaces are best described by a diffuse transition between regions of different SLDs. The diffuse transitions between SLDs observed here suggest that the interfacial structure of IL, while ordered, is best described by a continuous variation in composition as opposed to distinct layers of ions.

25 ENERGY STORAGE↗
DOE OSTI2024

Complexes of tubulin oligomers and tau form a viscoelastic intervening network cross-bridging microtubules into bundles

Abstract The axon-initial-segment (AIS) of mature neurons contains microtubule (MT) fascicles (linear bundles) implicated as retrograde diffusion barriers in the retention of MT-associated protein (MAP) tau inside axons. Tau dysfunction and leakage outside of the axon is associated with neurodegeneration. We report on the structure of steady-state MT bundles in varying concentrations of Mg 2+ or Ca 2+ divalent cations in mixtures containing αβ-tubulin, full-length tau, and GTP at 37 °C in a physiological buffer. A concentration-time kinetic phase diagram generated by synchrotron SAXS reveals a wide-spacing MT bundle phase (B ws ), a transient intermediate MT bundle phase (B int ), and a tubulin ring phase. SAXS with TEM of plastic-embedded samples provides evidence of a viscoelastic intervening network (IN) of complexes of tubulin oligomers and tau stabilizing MT bundles. In this model, αβ-tubulin oligomers in the IN are crosslinked by tau’s MT binding repeats, which also link αβ-tubulin oligomers to αβ-tubulin within the MT lattice. The model challenges whether the cross-bridging of MTs is attributed entirely to MAPs. Tubulin-tau complexes in the IN or bound to isolated MTs are potential sites for enzymatic modification of tau, promoting nucleation and growth of tau fibrils in tauopathies.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1994

Electron injection and acceleration at nonlinear shocks: Results of numerical simulations

We present results of numerical simulations of electron injection and acceleration at nonlinear high Mach number shocks. The electrons are assumed to be heated at the thermal subshock to an energy E(sub inj), which is treated as a free parameter, above which they are injected by self-generated whistlers to momentum m(sub p) x V(sub A). This injection mechanism requires Mach numbers greater than (43/(beta(sub -)))(((k T(sub e))/(E(sub inj)))(exp 1/2)), where T(sub e) and beta (sub -) are the upstream electron temperature and plasma beta parameter. Above m(sub p) x V(sub A) electrons are trapped in the shock by Alfven waves. In the proton precursor region the Alfven waves are assumed to be generated by protons accelerated at the shock, and have nonlinear intensities. Below GeV, however, electrons of a given rigidity propagate faster than protons with a similar rigidity and therefore diffuse to regions ahead of the proton precursor. In those regions the Alfven waves are generated by the electrons themselves. The diffusion coefficient appears to increase with decreasing acceleration efficiency. As a result, the number of electrons accelerated to energies GeV and above and, hence, the electron to proton ratio, depend only weakly on the extent of electron heating at the subshock. The negative feedback also renders the electron spectra insensitive to shock compression ratio and smoothing length scale. The estimated e/p ratio at GeV is between approximately 1%-10%.

Levinson, Amir↗
NASA NTRSDate not supplied

Influence of Thermophysical Property Variability on Thermal Predictions of Ti-6Al-4V

The additive manufacturing (AM) industry has experienced rapid growth in recent decades as industrial interest in the process has grown. Because of this, there are many groups interested in simulations of the AM process. However, the influence of thermophysical property variability on the melt pool geometry during processing of various AM alloys is unclear. The goal of my work at NASA Langley Research Center (LaRC) was to characterize this influence on Ti-6Al-4V during AM processing alongside developing a tool to enable equivalent studies on other AM materials. To facilitate this process, a database tool was developed to contain and organize thermophysical data previously reported by primary sources. The specific thermophysical properties of interest for this study were density, specific heat capacity, and conductivity. Thermal diffusivity was calculated using the other three properties. This data was then fit with a Gaussian distribution and then randomly sampled from using Monte Carlo random value sampling. Using the Rosenthal equation, this sample data was used to simulate the temperature field during additive manufacturing and extract the melt pool geometry. Specifically, the melt pool’s length, width, and depth. This process was repeated an arbitrary number of times, with the default being 1000. Once this information was obtained, histograms were made showing the distributions of the sizes of the melt pool’s length, width, and depth. Representative statistical metrics of the distributions were calculated (mean, standard deviation, and coefficient of variance). This work found that at the simulated processing values, the melt pool’s width and depth had a 5.7% variation while the length had a 3.2% variation. Additionally, a tool was constructed to graph a thermal color map representation of the melt pool for specific, arbitrary values of density, specific heat, and conductivity. This tool allowed for more efficient plotting of individual queries of the thermophysical properties. The results that were observed in the research were that values reported in literature vary and this variance can have a significant impact on the simulation. Predictions of the melt pool’s dimensions show that length has a standard deviation of approximately 6 μm, width has a standard deviation of approximately 5 μm, and depth can vary by approximately 3 μm. Considering the mean sizes of length, width, and depth are 185 μm, 88 μm, and 44 μm respectively, such a deviation is significant. This shows that depth, for example, could be more than 10% larger or smaller than expected. This demonstrates that variability in the reported thermophysical properties are not negligible and should be expected to have an influence on the results of the laser powder bed fusion additive manufacturing process. When simulating this process in the future, measures should be taken to account for this discrepancy and the uncertainty involved

Justin Martin↗
DOE OSTI2020

Phonon-mediated lipid raft formation in biological membranes

Short-wavelength collective molecular motions, also known as phonons, have recently attracted much interest in revealing dynamic properties of biological membranes through the use of neutron and x-ray scattering, infrared and Raman spectroscopies, and molecular dynamics simulations. Experimentally detecting unique vibrational patterns such as, shear phonon excitations, viscoelastic crossovers, transverse acoustic phonon gaps, and continuous and truncated optical phonon modes in cellular membranes, to name a few, has proven non-trivial. Here, we review recent advances in liquid thermodynamics that have resulted in the development of the phonon theory of 1 liquids. The theory has important predictions regarding the shear vibrational spectra of fluids, namely the emergence of viscoelastic crossovers and transverse acoustic phonon gaps. Furthermore, we show that these vibrational patterns are common in soft (non-crystalline) materials, including, but not limited to liquids, colloids, liquid crystals (mesogens), block copolymers, and biological membranes. The existence of viscoelastic crossovers and acoustic phonon gaps define the self-diffusion properties of cellular membranes and provide a molecular picture of the transient nature of lipid rafts.1 Importantly, the timescales (picoseconds) for the formation and dissolution of transient lipid rafts match the lifetime of the formation and breakdown of interfacial water hydrogen bonds. Apart from acoustic propagating phonon modes, biological membranes can also support more energetic non-propagating optical phonon excitations, also known as standing waves or breathing modes. Importantly, optical phonons can be truncated due to the existence of finite size nanodomains made up of strongly correlated lipid-cholesterol molecular pairs. These strongly coupled molecular pairs can serve as nucleation centers for the formation of stable rafts at larger length scales, due to correlations of spontaneous fluctuations (Onsager’s regression hypothesis). Finally and importantly, molecular level viscoelastic crossovers, acoustic phonon gaps, and continuous and truncated optical phonon modes may offer insights as to how lipid-lipid and lipid-protein interactions enable biological function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Understanding Viscoelasticity of an Entangled Silicone Copolymer via Coarse-Grained Molecular Dynamics Simulations

Entangled dynamics is important for understanding rheological properties of long-chain polymers. For entangled homopolymers, the classic tube-reptation model and its refinements have been successfully applied to quantify properties like diffusion coefficient and zero-rate viscosity. However, the application of such models to copolymers has been limited despite scientific and industrial importance. Here, we study the entangled melt dynamics of poly(dimethyl-co-diphenyl)siloxane random copolymer for a range of mean-composition-ratio ϕ of the diphenyl component via long-term molecular dynamics simulation with a recently developed coarse-grained model. We found that the segmental relaxation is heterogeneous at the monomeric level because of compositional fluctuations. However, at the chain-entanglement level and higher length scales, the viscoelastic response is homogeneous with compositional dependence only through the overall diphenyl fraction ϕ. The relaxation modulus of the entangled copolymer melt conforms to the Likhtman–McLeish model, and the viscosity predicted using our current coarse-grained parameters is in good quantitative agreement with experimental data.

Copolymers↗
DOE OSTI2022

Microstructural pattern formation during liquid metal dealloying: Phase-field simulations and theoretical analyses

In recent years, liquid metal dealloying has emerged as a promising material processing method to generate micro- and nanoscale bicontinuous or porous structures. Most previous studies focused on the experimental characterization of the dealloying process and on the properties of the dealloyed materials, leaving the theoretical study incomplete to fully understand the fundamental mechanisms of the liquid metal dealloying process. In this paper, we use theoretical models and phase-field simulations to clarify the kinetics and pattern formation during liquid metal dealloying. Our investigation starts from a theoretical analysis of the 1D dissolution of a binary precursor alloy, which reveals that the 1D dissolution process involves two regimes. In the first regime, due to the low solubility of one of the elements in the melt, it accumulates at the solid-liquid interface, which reduces the dissolution kinetics. In the second regime, the interface kinetics reaches a stationary regime where both elements of the precursor alloy dissolve into the melt. Previous works revealed that in the early dealloying stage, the dealloying front is destabilized by an interfacial spinodal decomposition, which triggers the formation of interconnected ligaments. Here we extend this line of work by proposing a linear stability analysis able to predict the initial length scale of the ligaments formed in the initial stage of the dealloying. Combining this analysis with the 1D dissolution model proposed here enables us to better understand the initial conditions (composition of the precursor alloy and the melt) leading to a planar dissolution without interface destabilization. Finally, we report a strong influence of solid-state diffusion on dealloying. Although the solid-state diffusivity is four to five orders of magnitude smaller than in the liquid phase, it is found to affect both dissolution kinetics and ligament morphologies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

MMMnet: A Neural Network Surrogate for Real-Time Transport Prediction Based on the Updated Multi-Mode Model

The Multi-Mode Model (MMM) is a physics-based anomalous transport model integrated into TRANSP for predicting electron and ion thermal transport, electron and impurity particle transport, and toroidal and poloidal momentum transport. While MMM provides valuable predictive capabilities, its computational cost, although manageable for standard simulations, is too high for real-time control applications. MMMnet, a neural network-based surrogate model, is developed to address this challenge by significantly reducing computation time while maintaining high accuracy. Trained on TRANSP simulations of DIII-D discharges, MMMnet incorporates an updated version of MMM (9.0.10) with enhanced physics, including isotopic effects, plasma shaping via effective magnetic shear, unified correlation lengths for ion-scale modes, and a new physics-based model for the electromagnetic electron temperature gradient mode. A key advancement is MMMnet’s ability to predict all six transport coefficients, providing a comprehensive representation of plasma transport dynamics. MMMnet achieves a two-order-of-magnitude speed improvement while maintaining strong correlation with MMM diffusivities, making it well-suited for real-time tokamak control and scenario optimization.

DIII-D↗
NASA NTRS1978

On the propagation of cosmic rays in the Galaxy

Cosmic ray propagation in the Galaxy is examined within the framework of the leaky-box model, whereby phenomenological aspects of diffusion and propagation are related to the physical properties of the interstellar medium. The model is parameterized by data on composition, spectra, spectral change as a function of energy, and cosmic ray particle lifetime. It is assumed that electrons and nucleons have a common source within the Galaxy (probably associated with supernova remnants). Lifetime data and leakage length are compared with radio synchrotron and gamma-ray data to establish a place of origin for cosmic rays observed at earth, i.e., at a radius of 0.5-1.0 kpc from earth. It is noted that if the escape length becomes energy independent, all nuclear spectra will flatten, while the electron spectrum remains unchanged.

Ormes, J.↗
NASA NTRS2003

Effect of Crucible Diameter Reduction on the Convection, Macrosegregation, and Dendritic Morphology during Directional Solidification of Pb-2.2 Wt Pct Sb Alloy

The Pb-2.2 wt pct Sb alloy has been directionally solidified in 1-, 2-, 3-, and 7-mm-diameter crucibles with planar and dendritic liquid-solid interface orphology. For plane front solidification, the experimentally observed macrosegregation along the solidified length follows the relationship proposed by Favier. Application of a 0.4 T transverse magnetic field has no effect on the extent of convection. Reducing the ampoule diameter appears to decrease the extent of convection. However, extensive convection is still present even in the 1-mm-diameter crucible. An extrapolation of the observed behavior indicated that nearly diffusive transport conditions require ampoules that are about 40 microns in diameter. Reduction of the crucible diameter does not appear to have any significant effect on the primary dendrite spacing. However, it results in considerable distortion of the dendrite morphology and ordering. This is especially true for the 1-mm diameter samples.

Chen, Jun↗
DOE OSTI2021

Investigation of Scale-Dependent Groundwater/Surface-water Exchange in Rivers by Gradient Self-Potential Logging: Numerical Modeling and Field Experiments

Exchanges of groundwater and surface-water are fundamental to a wide range of water-supply and water-quality management issues but challenging to map beyond the reach scale. Waterborne gradient self-potential (SP) measurements are directly sensitive to water flow through riverbed sediments and can be used to infer exchange locations, direction (gain versus loss), scale, and relative changes, but to date applications to river corridor hydrology are limited. Numerical modeling and field experiments were therefore performed herein, each emphasizing waterborne gradient SP logging for identifying and locating focused vertical groundwater discharge (surface-water gain) and recharge (surface-water loss) in a river. Two and three-dimensional numerical models were constructed to simulate the polarities, appearances, and peak amplitudes of streaming-potential and electric-field anomalies on a riverbed and in the surface-water that were attributable to steady-state vertical fluxes of groundwater through high-permeability conduits in the riverbed. Effects of varied hydraulic length-scale of exchange and surface-water depth were tested through numerical modeling. Modeling results aided in data acquisition and interpretation for three repeated field experiments performed along a 1.5–2.0 km reach of the Quashnet River in Cape Cod, Massachusetts, where focused, meter-scale groundwater discharges occur at discrete locations within otherwise ubiquitous and more diffuse groundwater upwelling conditions. Strong gradient SP anomalies were repeatedly measured in the Quashnet River at previously confirmed locations of focused groundwater discharge, showing the efficacy of waterborne gradient SP logging in identifying and characterizing groundwater/surface water exchange dynamics at multiple river network scales.

Engineering↗
DOE OSTI2022

Numerical and Experimental Investigations on the Ignition Behavior of OME

On the path towards climate-neutral future mobility, the usage of synthetic fuels derived from renewable power sources, so-called e-fuels, will be necessary. Oxygenated e-fuels, which contain oxygen in their chemical structure, not only have the potential to realize a climate-neutral powertrain, but also to burn more cleanly in terms of soot formation. Polyoxymethylene dimethyl ethers (PODE or OMEs) are a frequently discussed representative of such combustibles. However, to operate compression ignition engines with these fuels achieving maximum efficiency and minimum emissions, the physical-chemical behavior of OMEs needs to be understood and quantified. Especially the detailed characterization of physical and chemical properties of the spray is of utmost importance for the optimization of the injection and the mixture formation process. The presented work aimed to develop a comprehensive CFD model to specify the differences between OMEs and dodecane, which served as a reference diesel-like fuel, with regards to spray atomization, mixing and auto-ignition for single- and multi-injection patterns. The simulation results were validated against experimental data from a high-temperature and high-pressure combustion vessel. The sprays’ liquid and vapor phase penetration were measured with Mie-scattering and schlieren-imaging as well as diffuse back illumination and Rayleigh-scattering for both fuels. To characterize the ignition process and the flame propagation, measurements of the OH* chemiluminescence of the flame were carried out. Significant differences in the ignition behavior between OMEs and dodecane could be identified in both experiments and CFD simulations. Liquid penetration as well as flame lift-off length are shown to be consistently longer for OMEs. Zones of high reaction activity differ substantially for the two fuels: Along the spray center axis for OMEs and at the shear boundary layers of fuel and ambient air for dodecane. Additionally, the transient behavior of high temperature reactions for OME is predicted to be much faster.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1987

Analysis of non-premixed turbulent reacting flows

Studies of chemical reactions occurring in turbulent flows are important in the understanding of combustion and other applications. Current numerical methods are limited in their applications due to the numerical resolution required to completely capture all length scales, but, despite the fact that realistic combustion cannot be solved completely, numerical simulations can be used to give insight into the interaction between the processes of turbulence and chemical reaction. The objective was to investigate the effects of turbulent motion on the effects of chemical reaction to gain some insight on the interaction of turbulence, molecular diffusion, and chemical reaction to support modeling efforts. A direct turbulence simulation spectral code was modified to include the effects of chemical reaction and applied to an initial value problem of chemical reaction between non-premixed species. The influence of hydrodynamics on the instantaneous structure of the reaction was investigated. The local scalar dissipation rates and the local reaction rates were examined to determine the influence of vorticity or rate of strain on the reaction and the structure of the scalar field.

Leonard, Andy D.↗