Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reacting”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Velocity Measurement in a Dual-Mode Supersonic Combustor using Particle Image Velocimetry

Temporally and spatially-resolved, two-component measurements of velocity in a supersonic hydrogen-air combustor are reported. The combustor had a single unswept ramp fuel injector and operated with an inlet Mach number of 2 and a flow total temperature approaching 1200 K. The experiment simulated the mixing and combustion processes of a dual-mode scramjet operating at a flight Mach number near 5. The velocity measurements were obtained by seeding the fuel with alumina particles and performing Particle Image Velocimetry on the mixing and combustion wake of the ramp injector. To assess the effects of combustion on the fuel air-mixing process, the distribution of time-averaged velocity and relative turbulence intensity was determined for the cases of fuel-air mixing and fuel-air reacting. Relative to the mixing case, the near field core velocity of the reacting fuel jet had a slower streamwise decay. In the far field, downstream of 4 to 6 ramp heights from the ramp base, the heat release of combustion resulted in decreased flow velocity and increased turbulence levels. The reacting measurements were also compared with a computational fluid dynamics solution of the flow field. Numerically predicted velocity magnitudes were higher than that measured and the jet penetration was lower.

Goyne, C. P.↗
NASA NTRS2002

Science Goal Driven Observing: A Step Towards Maximizing Science Returns and Spacecraft Autonomy

In the coming decade, the drive to increase the scientific returns on capital investment and to reduce costs will force automation to be implemented in many of the scientific tasks that have traditionally been manually overseen. Thus, spacecraft autonomy will become an even greater part of mission operations. While recent missions have made great strides in the ability to autonomously monitor and react to changing health and physical status of spacecraft, little progress has been made in responding quickly to science driven events. The new generation of space-based telescopes/observatories will see deeper, with greater clarity, and they will generate data at an unprecedented rate. Yet, while onboard data processing and storage capability will increase rapidly, bandwidth for downloading data will not increase as fast and can become a significant bottleneck and cost of a science program. For observations of inherently variable targets and targets of opportunity, the ability to recognize early if an observation will not meet the science goals of variability or minimum brightness, and react accordingly, can have a major positive impact on the overall scientific returns of an observatory and on its operational costs. If the observatory can reprioritize the schedule to focus on alternate targets, discard uninteresting observations prior to downloading, or download them at a reduced resolution its overall efficiency will be dramatically increased. We are investigating and developing tools for a science goal monitoring (SGM) system. The SGM will have an interface to help capture higher-level science goals from scientists and translate them into a flexible observing strategy that SGM can execute and monitor. SGM will then monitor the incoming data stream and interface with data processing systems to recognize significant events. When an event occurs, the system will use the science goals given it to reprioritize observations, and react appropriately and/or communicate with ground systems - both human and machine - for confirmation and/or further high priority analyses.

Koratkar, Anuradha↗
NASA NTRS2002

Laboratory Anion Chemistry: Implications for the DIBs, and a Potential Formation Mechanism for a Known Interstellar Molecule

Due to recent interest in molecular anions as possible interstellar species, we have carried out several laboratory studies of anion chemistry. The reactions of the series C(sub n)(sup -); and C(sub n)H(sup -) with H and H2 were studied to address the viability of such species in the diffuse interstellar medium and to address their ability to be carriers of the diffuse interstellar bands (DIBs). These same molecules were also reacted with N and O to show possible heteroatomic products. C(sub m)N(sup - was a particularly stable product from the reaction of C(sub n)(sup -) + N. C3N(sup -) was further reacted with H to study chemistry that could produce HC3N, a known interstellar species. The reactions were done in a flowing afterglow selected ion flow tube apparatus (FA-SIFT). The anions were generated in an electron impact or cold cathode discharge source and the anion of interest was then selected by a quadrupole mass filter. The selected ion was then reacted with the atomic or molecular species in the flow tube and products were detected by another quadrupole. While the C(sub n)(sup -) species do not appear to be viable DIB carriers, their possible presence could provide a mechanism for the formation of known heteroatomic neutral molecules detected in the interstellar medium (ISM).

Eichelberger, B.↗
NASA NTRS2003

A 10 Kelvin 3 Tesla Magnet for Space Flight ADR Systems

Many future space flight missions are expected to use adiabatic demagnetization refrigerators (ADRs) to reach detector operating temperatures well below one Kelvin. The goal is to operate each ADR with a mechanical cooler as its heat sink, thus avoiding the use of liquid cryogens. Although mechanical coolers are being developed to operate at temperatures of 6 Kelvin and below, there is a large efficiency cost associated with operating them at the bottom of their temperature range. For the multi-stage ADR system being developed at Goddard Space Flight Center, the goal is to operate with a 10 Kelvin mechanical cooler heat sink. With currently available paramagnetic materials, the highest temperature ADR stage in such a system will require a magnetic field of approximately three Tesla. Thus the goal is to develop a small, lightweight three Tesla superconducting magnet for operation at 10 Kelvin. It is important that this magnet have a low current/field ratio. Because traditional NbTi magnets do not operate safely above about six Kelvin, a magnet with a higher Tc is required. The primary focus has been on Nb3Sn magnets. Since standard Nb3Sn wire must be coated with thick insulation, wound on a magnet mandrel and then reacted, standard Nb,Sn magnets are quite heavy and require high currents Superconducting Systems developed a Nb3Sn wire which can be drawn down to small diameter, reacted, coated with thin insulation and then wound on a small diameter coil form. By using this smaller wire and operating closer to the wire s critical current, it should be possible to reduce the mass and operating current of 10 Kelvin magnets. Using this "react-then-wind" technology, Superconducting Systems has produced prototype 10 Kelvin magnets. This paper describes the development and testing of these magnets and discusses the outlook for including 10 Kelvin magnets on space-flight missions.

Tuttle, Jim↗
NASA NTRS1995

Large steric effect in the substitution reaction of amines with phosphoimidazolide-activated nucleosides

Aliphatic amines react with phosphoimidazolide-activated derivatives of guanosine and cytidine (ImpN) by replacing the imidazole group. The kinetics of reaction of guanosine 5'-phospho-2-methylimidazolide (2-MeImpG) with glycine ethyl ester, glycinamide, 2-methoxyethylamine, n-butylamine, morpholine, dimethylamine (Me2NH), ethylmethylamine (EtNHMe), diethylamine (Et2NH), pyrrolidine, and piperidine were determined in water at 37 degrees C. With primary amines, a plot of the logarithm of the rate constant for attack by the amine on the protonated substrate, log kSH(A), versus the pKa of the amine exhibits a good linear correlation with a Bronsted slope, beta nuc = 0.48. Most of the secondary amines tested react with slightly higher reactivity than primary amines of similar pKa. Interestingly, some secondary amines show substantially lower reactivity than might be expected: EtNHMe reacts about eight times, and Et2NH at least 100 times, more slowly than Me2NH although all three amines are of similar basicity. For comparison, the kinetics of reaction of guanosine 5'-phosphoimidazolide (ImpG) and cytidine 5'-phosphoimidazolide (ImpC) were determined with Me2NH, EtNHMe, and Et2NH, and similar results were obtained. These results establish that the increased steric hindrance observed with the successive addition of ethyl groups are not due to any special steric requirements imposed by the guanosine or the methyl on the 2-methylimidazole leaving group of 2-MeImpG. It is concluded that addition of ethyl and, perhaps, groups larger than ethyl dramatically increases the kinetic barrier for addition of aliphatic secondary amines to the P-N bond of ImpN. This study supports the observation that the primary amino groups on the natural polyamines are at least 2 orders of magnitude more reactive than the secondary amino groups in the reaction with ImpN.

Non-NASA Center↗
NASA NTRS1992

Evidence of thyroxine formation following iodine administration in Sprague-Dawley rats

Iodine (I2) has been proposed to be used as a water disinfectant on the manned space station. Previous work has shown that subchronic administration of I2 to Sprague-Dawley rats in drinking water significantly increases plasma thyroxine/triiodothyronine (T4/T3) levels. This is not observed with iodide (I-) treatment. The present study addresses the possibility that I2 reacts with deiodinated T4 metabolites in the gastrointestinal tract to resynthesize T4. Incubation of diiodothyronine (T2), T3, or reverse T3 with I2 in phosphate-buffered saline resulted in the formation of T4 as measured by radioimmunoassay. Washes from the initial segments of the small intestine of the rat show that substrates are present that react with I2 to produce T4. Single oral doses of I2 to rats produced significant dose-related increases in serum T4 and decreases in T3 concentrations after 2 h. Administration of an equivalent dose of I- did not alter significantly plasma T4 concentrations. Higher concentrations of a radioactive substance that bound a T4-specific antibody are present in plasma of animals treated with 125I2 compared to 125I-. These data support the hypothesis that I2 reacts with metabolites of thyroid hormone in the gastrointestinal tract to resynthesize T4 and elevate its levels in blood.

NASA Discipline Environmental Health↗
NASA NTRS1957

Steady Nuclear Combustion in Rockets

The astrophysical theory of stationary nuclear reactions in stars is applied to the conditions that would be met in the practical engineering cases that would differ from the former, particularly with respect to the much lower combustion pressures, dimensions of the reacting volume, and burnup times. This application yields maximum rates of hear production per unit volume of reacting gas occurring at about 10(exp 8) K in the cases of reactions between the hydrogen isotopes, but yields higher rates for heavier atoms. For the former, with chamber pressures of the order of 100 atmospheres, the energy production for nuclear combustion reaches values of about 10(exp 4) kilocalories per cubic meter per second, which approaches the magnitude for the familiar chemical fuels. The values are substantially lower for heavier atoms, and increase with the square of the combustion pressure. The half-life of the burnup in the fastest reactions may drop to values as low as those for chemical fuels so that, despite the high temperature, the radiated energy can remain smaller than the energy produced, particularly if an inefficiently radiating (i.e., easily completely ionized reacting material like hydrogen), is used. On the other hand, the fraction of completely ionized particles in the gases undergoing nuclear combustion must not exceed a certain upper limit because the densities (approximately 10(exp -10) grams per cubic centimeter)) lie in the range of high vacua and only for the previously mentioned fraction of nonionized particles can mean free paths be retained small enough so that the chamber diameters of several dozen meters will suffice. Under these conditions it appears that continuously maintained stable nuclear reactions at practical pressures and dimensions are fundamentally possible and their application can be visualized as energy sources for power plants and propulsion units.

Saenger, E.↗
NASA NTRS1997

Volatile Reaction Products From Silicon-Based Ceramics in Combustion Environments Identified

Silicon-based ceramics and composites are prime candidates for use as components in the hot sections of advanced aircraft engines. These materials must have long-term durability in the combustion environment. Because water vapor is always present as a major product of combustion in the engine environment, its effect on the durability of silicon-based ceramics must be understood. In combustion environments, silicon-based ceramics react with water vapor to form a surface silica (SiO2) scale. This SiO2 scale, in turn, has been found to react with water vapor to form volatile hydroxides. Studies to date have focused on how water vapor reacts with high-purity silicon carbide (SiC) and SiO2 in model combustion environments. Because the combustion environment in advanced aircraft engines is expected to contain about 10-percent water vapor at 10-atm total pressure, the durability of SiC and SiO2 in gas mixtures containing 0.1- to 1-atm water vapor is of interest. The reactions of SiC and SiO2 with water vapor were monitored by measuring weight changes of sample coupons in a 0.5-atm water vapor/0.5-atm oxygen gas mixture with thermogravimetric analysis.

Opila, Elizabeth J.↗
NASA NTRS2004

Thermal characteristics of Lithium-ion batteries

Lithium-ion batteries have a very promising future for space applications. Currently they are being used on a few GEO satellites, and were used on the two recent Mars rovers Spirit and Opportunity. There are still problem that exist that need to be addressed before these batteries can fully take flight. One of the problems is that the cycle life of these batteries needs to be increased. battery. Research is being focused on the chemistry of the materials inside the battery. This includes the anode, cathode, and the cell electrolyte solution. These components can undergo unwanted chemical reactions inside the cell that deteriorate the materials of the battery. During discharge/ charge cycles there is heat dissipated in the cell, and the battery heats up and its temperature increases. An increase in temperature can speed up any unwanted reactions in the cell. Exothermic reactions cause the temperature to increase; therefore increasing the reaction rate will cause the increase of the temperature inside the cell to occur at a faster rate. If the temperature gets too high thermal runaway will occur, and the cell can explode. The material that separates the electrode from the electrolyte is a non-conducting polymer. At high temperatures the separator will melt and the battery will be destroyed. The separator also contains small pores that allow lithium ions to diffuse through during charge and discharge. High temperatures can cause these pores to close up, permanently damaging the cell. My job at NASA Glenn research center this summer will be to perform thermal characterization tests on an 18650 type lithium-ion battery. High temperatures cause the chemicals inside lithium ion batteries to spontaneously react with each other. My task is to conduct experiments to determine the temperature that the reaction takes place at, what components in the cell are reacting and the mechanism of the reaction. The experiments will be conducted using an accelerating rate calorimeter (ARC), which uses a heat-wait-search mode until an exothermic reaction is detected. After an exotherm is found the calorimeter maintains an adiabatic environment around a bomb which holds the test sample. The ARC will help identify important reactions and what temperature these exothermic reactions take place at. In order fully understand the battery, we are first going to take apart the battery and test the individual components of the battery using the ARC. I will first conduct a test on the electrolyte solution by itself. We will then test the electrolyte solution with the anode. We would like to see how the electrolyte solution reacts with the anode and its binder material. The next would be the same test using the cathode instead of the anode. By comparing the results of the electrolyte, electrolyte with anode, and the electrolyte with the cathode we can determine the reactions that are taking place due to each component. Using the heat capacity of the each individual sample and the temperature by which the sample increases, kinetic and thermo-dynamical information can then be found. A Gas chromatograph could be used to help with the task of identifying the by-products at the end of each test. One way of increasing the cycle life is to increase the stability of the materials inside the

Hauser, Dan↗
NASA NTRS2004

Non-Additive Voltametric Currents From a Mixture of Two, Three and Four Redox-Active Compounds and Electroanalytical Implications

We have published recently the effect of dissimilar diffusion coefficients on the size of the voltammetric waves from a mixture of two redox-active compounds. Similarly, at the potential range where three redox-active species, decamethylferrocene (dMeFc), ferrocene (Fc) and N-methylphenothiazine (MePTZ), are oxidized simultaneously with rates controlled by linear diffusion, electrogenerated radicals diffusing outwards from the electrode react with the original species diffusing towards the electrode from the bulk; thus, Fc(+) reacts with dMeFc producing Fc and dMeFc(+), while MePTZ(+) reacts both with dMeFc producing MePTZ and dMeFc(+), and with Fc producing MePTZ and Fc(+). These reactions replace dMeFc with Fc at the second plateau, and both dMeFc and Fc with MePTZ at the third plateau. Since the diffusion coefficients of the three species are not equal, the mass-transfer limited currents of the second and the third oxidation wave plateaus change by approx. 10%. Numerical simulations of the experimental voltamograms support this mechanism. Similar results were also obtained for a mixture of four redoxactive compounds. The implications of this non-additive nature of currents on: (a) the use of internal voltammetric standards for quantitative analysis of a mixture of redox-active compounds; and, (b) the half wave potentials (E1/2) of the 2nd, 3rd and 4th waves for qualitative analysis, will be discussed.

Dass, Amala↗
NASA NTRS2005

Incorporation of a Chemical Equilibrium Equation of State into LOCI-Chem

Renewed interest in development of advanced high-speed transport, reentry vehicles and propulsion systems has led to a resurgence of research into high speed aerodynamics. As this flow regime is typically dominated by hot reacting gaseous flow, efficient models for the characteristic chemical activity are necessary for accurate and cost effective analysis and design of aerodynamic vehicles that transit this regime. The LOCI-Chem code recently developed by Ed Luke at Mississippi State University for NASA/MSFC and used by NASA/MSFC and SSC represents an important step in providing an accurate, efficient computational tool for the simulation of reacting flows through the use of finite-rate kinetics [3]. Finite rate chemistry however, requires the solution of an additional N-1 species mass conservation equations with source terms involving reaction kinetics that are not fully understood. In the equilibrium limit, where the reaction rates approach infinity, these equations become very stiff. Through the use of the assumption of local chemical equilibrium the set of governing equations is reduced back to the usual gas dynamic equations, and thus requires less computation, while still allowing for the inclusion of reacting flow phenomenology. The incorporation of a chemical equilibrium equation of state module into the LOCI-Chem code was the primary objective of the current research. The major goals of the project were: (1) the development of a chemical equilibrium composition solver, and (2) the incorporation of chemical equilibrium solver into LOCI-Chem. Due to time and resource constraints, code optimization was not considered unless it was important to the proper functioning of the code.

Cox, Carey F.↗
NASA NTRS2005

Process Developed for Generating Ceramic Interconnects With Low Sintering Temperatures for Solid Oxide Fuel Cells

Solid oxide fuel cells (SOFCs) have been considered as premium future power generation devices because they have demonstrated high energy-conversion efficiency, high power density, and extremely low pollution, and have the flexibility of using hydrocarbon fuel. The Solid-State Energy Conversion Alliance (SECA) initiative, supported by the U.S. Department of Energy and private industries, is leading the development and commercialization of SOFCs for low-cost stationary and automotive markets. The targeted power density for the initiative is rather low, so that the SECA SOFC can be operated at a relatively low temperature (approx. 700 C) and inexpensive metallic interconnects can be utilized in the SOFC stack. As only NASA can, the agency is investigating SOFCs for aerospace applications. Considerable high power density is required for the applications. As a result, the NASA SOFC will be operated at a high temperature (approx. 900 C) and ceramic interconnects will be employed. Lanthanum chromite-based materials have emerged as a leading candidate for the ceramic interconnects. The interconnects are expected to co-sinter with zirconia electrolyte to mitigate the interface electric resistance and to simplify the processing procedure. Lanthanum chromites made by the traditional method are sintered at 1500 C or above. They react with zirconia electrolytes (which typically sinter between 1300 and 1400 C) at the sintering temperature of lanthanum chromites. It has been envisioned that lanthanum chromites with lower sintering temperatures can be co-fired with zirconia electrolyte. Nonstoichiometric lanthanum chromites can be sintered at lower temperatures, but they are unstable and react with zirconia electrolyte during co-sintering. NASA Glenn Research Center s Ceramics Branch investigated a glycine nitrate process to generate fine powder of the lanthanum-chromite-based materials. By simultaneously doping calcium on the lanthanum site, and cobalt and aluminum on the chromium site, we could sinter the materials below 1400 C. The doping concentrations were adjusted so that the thermal expansion coefficient matched that of the zirconia electrolyte. Also, the investigation was focused on stoichiometric compositions so that the materials would have better stability. Co-sintering and chemical compatibility with zirconia electrolyte were examined by X-ray diffraction, scanning electron microscopy, and energy dispersive spectroscopy (line scanning and dot map). The results showed that the materials bond well, but do not react, with zirconia electrolyte. The electric conductivity of the materials measured at 900 C in air was about 20 S/cm.

Zhong, Zhi-Min↗
NASA NTRS2000

Tough soluble aromatic thermoplastic copolyimides

Tough, soluble, aromatic, thermoplastic copolyimides were prepared by reacting 4,4'-oxydiphthalic anhydride, 3,4,3',4'-biphenyltetracarboxylic dianhydride and 3,4'-oxydianiline. Alternatively, these copolyimides may be prepared by reacting 4,4'-oxydiphthalic anhydride with 3,4,3',4'-biphenyltetracarboxylic dianhydride and 3,4'-oxydiisocyanate. Also, the copolyimide may be prepared by reacting the corresponding tetra acid and ester precursors of 4,4'-oxydiphthalic anhydride and 3,4,3',4'-biphenyltetracarboxylic dianhydride with 3,4'-oxydianiline. These copolyimides were found to be soluble in common amide solvents such as N,N'-dimethyl acetamide, N-methylpyrrolidinone, and dimethylformamide allowing them to be applied as the fully imidized copolymer and to be used to prepare a wide range of articles.

Bryant, Robert G.↗
NASA NTRS2008

Numerical Simulations of a Co-Axial Supersonic-Combusting Free-Jet Experiment

CFD calculations using the Reynolds-averaged Navier-Stokes equations coupled with species continuity equations have been made for a supersonic coaxial-jet CFD-validation experiment to determine the sensitivity of the external flowfield to the main-nozzle exit profile. Four different nozzle exit profiles were used in the study: a uniform profile, one computed using only the nozzle geometry, one computed using the nozzle geometry and part of the upstream facility combustor, and one using the nozzle and the full facility combustor. Two cases were examined using the four profiles: a non-reacting case without coflow and a reacting case with hydrogen coflow. Results show that the nozzle exit profile has a significant effect on the external flowfield. The uniform profile produced the longest jet while the profile created with the full combustor produced the shortest jet. The nozzle-only and part-combustor profiles fell between the other two profiles. The reacting flow was found to be more sensitive to the nozzle exit profile since it affects the downstream mixing and combustion. These calculations indicate the importance of properly setting the nozzle-exit profile for this type of calculation.

Gaffney, Richard L.↗
NASA NTRS2011

Numerical Dissipation and Wrong Propagation Speed of Discontinuities for Stiff Source Terms

In compressible turbulent combustion/nonequilibrium flows, the constructions of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and (b) obtaining correct propagation speed of discontinuities for stiff reacting terms on coarse grids share one important ingredient - minimization of numerical dissipation while maintaining numerical stability. Here coarse grids means standard mesh density requirement for accurate simulation of typical non-reacting flows. This dual requirement to achieve both numerical stability and accuracy with zero or minimal use of numerical dissipation is most often conflicting for existing schemes that were designed for non-reacting flows. The goal of this paper is to relate numerical dissipations that are inherited in a selected set of high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities for two representative stiff detonation wave problems.

Yee, H. C.↗
NASA NTRS2012

Damage Tolerance Assessment of Friction Pull Plug Welds in an Aluminum Alloy

Friction stir welding is a solid state welding process used in the fabrication of cryogenic propellant tanks. Self-reacting friction stir welding is one variation of the friction stir weld process being developed for manufacturing tanks. Friction pull plug welding is used to seal the exit hole that remains in a circumferential self-reacting friction stir weld. A friction plug weld placed in a self-reacting friction stir weld results in a non-homogenous weld joint where the initial weld, plug weld, their respective heat affected zones and the base metal all interact. The welded joint is a composite plastically deformed material system with a complex residual stress field. In order to address damage tolerance concerns associated with friction plug welds in safety critical structures, such as propellant tanks, nondestructive inspection and proof testing may be required to screen hardware for mission critical defects. The efficacy of the nondestructive evaluation or the proof test is based on an assessment of the critical flaw size. Test data relating residual strength capability to flaw size in an aluminum alloy friction plug weld will be presented.

McGill, Preston↗
NASA NTRS2012

Numerical Dissipation and Wrong Propagation Speed of Discontinuities for Stiff Source Terms

In compressible turbulent combustion/nonequilibrium flows, the constructions of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and (b) obtaining correct propagation speed of discontinuities for stiff reacting terms on coarse grids share one important ingredient - minimization of numerical dissipation while maintaining numerical stability. Here coarse grids means standard mesh density requirement for accurate simulation of typical non-reacting flows. This dual requirement to achieve both numerical stability and accuracy with zero or minimal use of numerical dissipation is most often conflicting for existing schemes that were designed for non-reacting flows. The goal of this paper is to relate numerical dissipations that are inherited in a selected set of high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities as a function of stiffness of the source term and the grid spacing.

Yee, H. C.↗
NASA NTRS2013

Spurious Behavior of Shock-Capturing Methods: Problems Containing Stiff Source Terms and Discontinuities

The goal of this paper is to relate numerical dissipations that are inherited in high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities. For pointwise evaluation of the source term, previous studies indicated that the phenomenon of wrong propagation speed of discontinuities is connected with the smearing of the discontinuity caused by the discretization of the advection term. The smearing introduces a nonequilibrium state into the calculation. Thus as soon as a nonequilibrium value is introduced in this manner, the source term turns on and immediately restores equilibrium, while at the same time shifting the discontinuity to a cell boundary. The present study is to show that the degree of wrong propagation speed of discontinuities is highly dependent on the accuracy of the numerical method. The manner in which the smearing of discontinuities is contained by the numerical method and the overall amount of numerical dissipation being employed play major roles. Moreover, employing finite time steps and grid spacings that are below the standard Courant-Friedrich-Levy (CFL) limit on shockcapturing methods for compressible Euler and Navier-Stokes equations containing stiff reacting source terms and discontinuities reveals surprising counter-intuitive results. Unlike non-reacting flows, for stiff reactions with discontinuities, employing a time step and grid spacing that are below the CFL limit (based on the homogeneous part or non-reacting part of the governing equations) does not guarantee a correct solution of the chosen governing equations. Instead, depending on the numerical method, time step and grid spacing, the numerical simulation may lead to (a) the correct solution (within the truncation error of the scheme), (b) a divergent solution, (c) a wrong propagation speed of discontinuities solution or (d) other spurious solutions that are solutions of the discretized counterparts but are not solutions of the governing equations. The present investigation for three very different stiff system cases confirms some of the findings of Lafon & Yee (1996) and LeVeque & Yee (1990) for a model scalar PDE. The findings might shed some light on the reported difficulties in numerical combustion and problems with stiff nonlinear (homogeneous) source terms and discontinuities in general.

Yee, Helen M. C.↗