Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion couples”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Updates on the Predictive Materials Modeling Software Tools

Updates on NASA‘s efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. The PMM effort is part of the Entry Systems Modeling (ESM) project under NASA’s Game Changing Development (GCD) program. To reduce the need for extensive testing and accelerate the design cycle process, ESM is developing simulation and modeling tools that enable the characterization of the properties of thermal protection materials and their response to extremely hot plasma. The Porous Microstructure Analysis (PuMA) software has been developed to compute effective material properties and perform material response simulations on digitized microstructures of porous media. PuMA is able to import three-dimensional digital images obtained from X-ray microtomography or to generate artificial microstructures that mimic real materials. PuMA also provides a module for interactive 3D visualizations. Version 3, which was recently released as open-source, includes modules to compute simple morphological properties such as porosity, volume fractions, pore diameter, and specific surface area. Additional capabilities include the determination of effective thermal and electrical conductivity (both radiative and solid conduction - including the ability to simulate local anisotropy for the latter); effective diffusivity and tortuosity from the continuum to the rarefied regime; techniques to determine the local material orientation, as well as mechanical properties (elasticity coefficients), and permeability. Computed properties are then used to inform a macro-scale material response model, such as those implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software developed within ESM. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale, assuming Local Thermal Equilibrium. Recent efforts include the development of a mechanical erosion model as well as a unified model allowing an intrinsic coupling between fluid and material. Comparison to flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI] and Mars 2020 MEDLI2) is critical in order to validate these computational tools. Examples of ablative material response using the code will be presented, including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrated the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains through the use of massively parallel computations.

material modeling↗
NASA NTRS1985

Critical Transport Phenomena in Fluid Helium Under Low Gravity

The feasibility of carrying out measurements of certain critical transport properties of pure fluid under conditions of low gravity was studied. These properties are the thermal conductivity, kappa, the shear viscosity zeta and the diffusive relaxation time tau, which are predicted to diverge (tend to infinity) as the liquid-vapor critical point is approached. However, in this critical region, the Earth's gravity effect becomes very important. As the critical point is approached, the gravity effects increasingly distort the results. The reason for this is that the compressibility of the fluid also diverges and under the influence of gravity causes a vertical density gradient in the fluid, which is significant even when very thin fluid layers (typically 1 mm high) are being used. The result is that the temperature dependence of kappa, zeta, and tau tends to flatten off as T sub c is approached instead of continuing to increase, and therefore the predictions from the renormalization group and mode coupling theories cannot be subjected to a satisfactory test.

Meyer, H.↗
NASA NTRS1994

Numerical Simulations of High-Speed Chemically Reacting Flow

The Essentially NonOscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacting flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.

Ton, V. T.↗
NASA NTRS1994

Numerical Simulations of High-Speed Chemically Reacting Flow

The essentially nonoscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacting flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.

Ton, V. T.↗
NASA NTRS2005

The current development status of the Orbiting Carbon Observatory (OCO) instrument optical design

The status of the OCO instrument optical design is presented in this paper. The optical bench assembly comprises three cooled grating spectrometers coupled to an all-reflective telescope/relay system. Dichroic beam splitters are used to separate the light from a common telescope into the three spectral bands. The three bore sighted spectrometers allow the total column CO2 absorption path to be corrected for optical path and surface pressure uncertainties, aerosols, and water vapor. The design of the instrument is based on classic flight proven technologies.

global carb↗
NASA NTRS2001

Quantum Mechanical Study of Nanoscale MOSFET

The steady state characteristics of MOSFETS that are of practical Interest are the drive current, off-current, dope of drain current versus drain voltage, and threshold voltage. In this section, we show that quantum mechanical simulations yield significantly different results from drift-diffusion based methods. These differences arise because of the following quantum mechanical features: (I) polysilicon gate depletion in a manner opposite to the classical case (II) dependence of the resonant levels in the channel on the gate voltage, (III) tunneling of charge across the gate oxide and from source to drain, (IV) quasi-ballistic flow of electrons. Conclusions dI/dV versus V does not increase in a manner commensurate with the increase in number of subbands. - The increase in dI/dV with bias is much smaller then the increase in the number of subbands - a consequence of bragg reflection. Our calculations show an increase in transmission with length of contact, as seen in experiments. It is desirable for molecular electronics applications to have a small contact area, yet large coupling. In this case, the circumferential dependence of the nanotube wave function dictates: - Transmission in armchair tubes saturates around unity - Transmission in zigzag tubes saturates at two.

Svizhenko, Alexei↗
NASA NTRS2005

Control of Combustion-Instabilities Through Various Passive Devices

It is well known that under some operating conditions, rocket engines (using solid or liquid fuels) exhibit unstable modes of operation that can lead to engine malfunction and shutdown. The sources of these instabilities are diverse and are dependent on fuel, chamber geometry and various upstream sources such as pumps, valves and injection mechanism. It is believed that combustion-acoustic instabilities occur when the acoustic energy increase due to the unsteady heat release of the flame is greater than the losses of acoustic energy from the system [1, 2]. Giammar and Putnam [3] performed a comprehensive study of noise generated by gasfired industrial burners and made several key observations; flow noise was sometimes more intense than combustion roar, which tended to have a characteristic frequency spectrum. Turbulence was amplified by the flame. The noise power varied directly with combustion intensity and also with the product of pressure drop and heat release rate. Karchmer [4] correlated the noise emitted from a turbofan jet engine with that in the combustion chamber. This is important, since it quantified how much of the noise from an engine originates in the combustor. A physical interpretation of the interchange of energy between sound waves and unsteady heat release rates was given by Rayleigh [5] for inviscid, linear perturbations. Bloxidge et al [6] extended Rayleigh s criterion to describe the interaction of unsteady combustion with one-dimensional acoustic waves in a duct. Solutions to the mass, momentum and energy conservation equations in the pre- and post-flame zones were matched by making several assumptions about the combustion process. They concluded that changes in boundary conditions affect the energy balance of acoustic waves in the combustor. Abouseif et al [7] also solved the one-dimensional flow equations, but they used a onestep reaction to evaluate the unsteady heat release rate by relating it to temperature and velocity perturbations. Their analysis showed that oscillations arise from coupling between entropy waves produced at the flame and pressure waves originating from the nozzle. Yang and Culick [8] assumed a thin flame sheet, which is distorted by velocity and pressure oscillations. Conservation equations were expressed in integral form and solutions for the acoustic wave equations and complex frequencies were obtained. The imaginary part of the frequency indicated stability regions of the flame. Activation energy asymptotics together with a one-step reaction were used by McIntosh [9] to study the effects of acoustic forcing and feedback on unsteady, one-dimensional flames. He found that the flame stability was altered by the upstream acoustic feedback. Shyy et al [10] used a high-accuracy TVD scheme to simulate unsteady, one-dimensional longitudinal, combustion instabilities. However, numerical diffusion was not completely eliminated. Recently, Prasad [11] investigated numerically the interactions of pressure perturbations with premixed flames. He used complex chemistry to study responses of pressure perturbations in one-dimensional combustors. His results indicated that reflected and transmitted waves differed significantly from incident waves.

Frendi, Kader↗
NASA NTRSDate not supplied

Microbial Communities in Microgravity: Simulation in Lab and on the Computer

Microorganisms grow differently in spaceflight than they do on Earth. While much remains unexplained about how microgravity affects microbial growth, one dominant hypothesis is that the lack of density-driven convection in the liquid growth environment makes mixing diffusion-limited, and therefore slower. This is supported by evidence that individual microbial strains experience starvation and acid stress in microgravity. However, if it is true, then microgravity would also have measurable effects on microbes in mixed communities, because many interspecies interactions involve the exchange of soluble metabolites through the medium (cross-feeding). Specifically, cooperative cross-feeding communities would grow more slowly in microgravity, and cooperation would be less stable on evolutionary timescales. Here we describe our efforts to test this hypothesis by simulating microgravity in silico and in the lab, using a model system of Eschericia coli and Salmonella enterica that grow only when they can exchange methionine and acetate. We created CAMDLES (CFD-DEM Artificial Microgravity Developments for Living Ecosystem Simulation) as an extension of CFDEM®coupling software, to carry out computational modeling of biological flows, growth, and mass transfer in microgravity and also in laboratory artificial microgravity devices (rotating wall vessels, RWV). Using CAMDLES, we found distinct differences in growth rates between RWV and true microgravity, and we were able to identify several features, such as spatial distribution, biofilm formation, and product yield parameters, that influence the degree to which RWV growth recapitulates microgravity growth. In addition, we report on the development of a laboratory system for monitoring growth rates and species ratios of the community in RWVs, using fluorescent strains. Pairing CAMDLES with the laboratory model system allows us to generate quantitative predictions about the effects of spaceflight on organisms that will be essential to sustaining human space exploration in the long term.

microbiology↗
NASA NTRS1999

Quantitative Interpretation of Optical Emission Sensors for Microgravity Experiments

Microgravity flight experiments uniquely test our knowledge and understanding of the coupling between chemistry and fluid mechanics. However, compared to ground based laboratory experiments the number of useful diagnostic tools suitable for microgravity environments is severely limited by the space, weight, power consumption, and operator complexity requirements. One of the available tools is the observation of optical emission, and total emission has already proven useful for observations of stable "flame balls" on the space shuttle. Wavelength resolved tomographic measurements of flame emission offer the promise of diagnostics to test our understanding of flame chemistry and structure. Individual emissions from electronically excited radicals, e.g., CH*, OH*, and C2*, can be identified in a methane/air flame. Spatially resolved measurements of the intensity of this resolved optical emission from a specific excited molecule enable chemically resolved flame structure studies. Wavelength resolved emission measurements to determine such structure in diffusion flames are being readied for flight experiments by a group at Yale headed by Profs. Smooke and Long. A quantitative relationship between emission intensity and flame properties, as expressed by a flame model, is needed for species specific optical emission measurements to fulfill its promise. The Yale group compared models and measurements of optical emission in laboratory tests at 1-g. Unfortunately, these experiments show disagreement between measurement and state-of-the-art flame models by over a factor of 50. Therefore, an improved chemical mechanism for optical emission from flames is needed to enable quantitative tests of microgravity flame models. The connection between excited state emission and flame chemistry is not yet adequate.

Jeffries, Jay B.↗
NASA NTRS2001

Experiments and Model Development for the Investigation of Sooting and Radiation Effects in Microgravity Droplet Combustion

Today, despite efforts to develop and utilize natural gas and renewable energy sources, nearly 97% of the energy used for transportation is derived from combustion of liquid fuels, principally derived from petroleum. While society continues to rely on liquid petroleum-based fuels as a major energy source in spite of their finite supply, it is of paramount importance to maximize the efficiency and minimize the environmental impact of the devices that burn these fuels. The development of improved energy conversion systems, having higher efficiencies and lower emissions, is central to meeting both local and regional air quality standards. This development requires improvements in computational design tools for applied energy conversion systems, which in turn requires more robust sub-model components for combustion chemistry, transport, energy transport (including radiation), and pollutant emissions (soot formation and burnout). The study of isolated droplet burning as a unidimensional, time dependent model diffusion flame system facilitates extensions of these mechanisms to include fuel molecular sizes and pollutants typical of conventional and alternative liquid fuels used in the transportation sector. Because of the simplified geometry, sub-model components from the most detailed to those reduced to sizes compatible for use in multi-dimensional, time dependent applied models can be developed, compared and validated against experimental diffusion flame processes, and tested against one another. Based on observations in microgravity experiments on droplet combustion, it appears that the formation and lingering presence of soot within the fuel-rich region of isolated droplets can modify the burning rate, flame structure and extinction, soot aerosol properties, and the effective thermophysical properties. These observations led to the belief that perhaps one of the most important outstanding contributions of microgravity droplet combustion is the observation that in the absence of asymmetrical forced and natural convection, a soot shell is formed between the droplet surface and the flame, exerting an influence on the droplet combustion response far greater than previously recognized. The effects of soot on droplet burning parameters, including burning rate, soot shell dynamics, flame structure, and extinction phenomena provide significant testing parameters for studying the structure and coupling of soot models with other sub-model components.

Choi, Mun Young↗
NASA NTRS1975

Effect of electric fields on the daytime high-latitude E and F regions

We have obtained solutions of the coupled continuity, momentum, and energy equations for NO(+), O(+), and O2(+) ions for conditions appropriate to the daytime high-latitude E and F regions. Owing to the rapid increase of the reaction O(+) + N2 yielding NO(+) + N with ion energy, high-latitude electric fields and consequent perpendicular-E x B drifts deplete O(+) in favor of NO(+). For electric field strengths less than about 10 mV/m the depletion of O(+) is small, and the altitude profiles of ion density are similar to those found at mid-latitudes. However, for moderate electric field strengths (50 mV/m), NO(+) is substantially increased in relation to O(+) and becomes an important ion throughout the F region. For large electric fields (200 mV/m), NO(+) completely dominates the ion composition to at least 600 km, decreasing at high altitudes with a diffusive equilibrium scale height. Since the overall F region electron density decreases markedly with increasing electric field strength, it appears that high-latitude, daytime electron density troughs are directly related to the presence of ionospheric electric fields.

Schunk, R. W.↗
NASA NTRSDate not supplied

Predicting Thermal Response in a Li-ion Cell on a UAV Fight Profile

As the energy storage devices continue to "pack" more energy in a small space, any damage, battery component failure, manufacturing defect, or electrically abusing the battery can lead to catastrophic thermal runaway events. A catastrophic thermal event in a cell leads to high temperature, in some instances spewing of battery materials due to gas development from side reactions initiated due to high internal temperatures. Also, a thermal runaway event can propagate from a single "failed" cell to the pack in a battery pack, leading to a more significant event. Mitigating a thermal runaway event is important in the commercial and automotive sectors. However, preventing such events in an electric aircraft (or air taxis) is paramount due to the lack of alternatives in the event of a failure. Battery prognostics algorithms allow the prediction of state-of-charge (SOC) and end-of-life (EOL) of a Li-ion battery in a UAV (unmanned air vehicle) [1]. For this presentation, we will extend this two-level battery predictive algorithm to predict SOC, EOL, and estimated maximum temperature during a simulated flight. The model is extended by integrating a lumped physics-driven thermal model for high current densities [2]. The parameters used to control SOC and EOL are maximum storable charge, time constant for Li-ion diffusivity in the carbon particles, and internal cell resistance. Cycling leads to an increase in the heat generated by an aged Li-ion cell with a LiyCoO2 (LCO) cathode and a LixC6 (MCMB) anode. The aging of a cell leads to increase in SEI layer thickness, the diffusion time for the lithium ions inside the electrodes, and the local reaction rates, in addition to the thermodynamic abuse caused by fixed cycling voltages controlled by a Battery Management System. As the battery ages, the cell resistance increases, while the onset temperature of the thermal runaway decreases (depends on the cell chemistry and cell abuse history). Any large deviation of the cell temperature from the estimated (expected) value can identify a faulty cell. Since SEI decomposition has the lowest onset temperature in the series of reactions leading to thermal runaway, the model considers the self-heating rate of the SEI decomposition as onset temperature (similar to Ref. [3]). The parameters in the Arrhenius equation for the SEI heating rate depend on the number of cycles, the cell's operating temperature, and the cell's abuse history [4,5]. Coupling the electrochemical, thermal, and aging model allow the prognostic algorithm to estimate a typical cell voltage and temperature as a function of age (cycling and calendar), whose departure from measured values from the BMS is used to identify a safety event. In addition, we will present the results from two simulated flight scenarios for a UAV: typical and extreme, since the power requirements vary significantly during take-off, landing, and changing altitudes, while the power requirements remain low during the cruise. For this presentation, the power requirement for a battery pack in a UAV is scaled to a single cell. This cell is cycled through a simulated profile, and the data is collected and used to predict a safety event.

Li-ion↗
NASA NTRS1995

The double oval UV auroral distribution. 1: Implications for the mapping of auroral arcs

During the later stages of the auroral substorm the luminosity distribution frequently resembles a double oval, one oval lying poleward of the normal or main UV auroral oval. We interpret the double oval morphology as being due to the plasma sheet boundary layer becoming active in the later stages of the substorm process. If the disturbance engulfs the nightside low-latitude boundary layers, then the double oval configuration extends into the dayside ionospheric region. The main UV oval is associated with the inner portion of the central plasma sheet and can rapidly change its auroral character from being diffuse to discrete. This transition is associated with the substorm process and is fundamental to understanding the near-Earth character of substorm onset. On the other hand, the poleward arc system in the nightside ionosphere occurs adjacent to or near the open-closed field line boundary. This system activates at the end of the optical expansion phase and is a part of the recovery phase configuration in substorms where it occurs. These two source regions for nightside discrete auroral arcs are important in resolving the controversy concerning the mapping of arcs to the magnetosphere. The dayside extension of this double oval configuration is also investigated and shows particle signatures which differ considerably from those on the nightside giving clues to the magnetospheric source regions of the aurora in the two local time sectors. Near-Earth substorm onsets are shown to be coupled to processes occurring much further tailward and indicate the importance of understanding the temporal development of features within the double oval. Using 'variance images,' a new technqiue for the investigation of these dynamics is outlined.

Elphinstone, R. D.↗
NASA NTRS2000

Thermographic Imaging of Defects in Anisotropic Composites

Composite materials are of increasing interest to the aerospace industry as a result of their weight versus performance characteristics. One of the disadvantages of composites is the high cost of fabrication and post inspection with conventional ultrasonic scanning systems. The high cost of inspection is driven by the need for scanning systems which can follow large curve surfaces. Additionally, either large water tanks or water squirters are required to couple the ultrasonics into the part. Thermographic techniques offer significant advantages over conventional ultrasonics by not requiring physical coupling between the part and sensor. The thermographic system can easily inspect large curved surface without requiring a surface following scanner. However, implementation of Thermal Nondestructive Evaluations (TNDE) for flaw detection in composite materials and structures requires determining its limit. Advanced algorithms have been developed to enable locating and sizing defects in carbon fiber reinforced plastic (CFRP). Thermal Tomography is a very promising method for visualizing the size and location of defects in materials such as CFRP. However, further investigations are required to determine its capabilities for inspection of thick composites. In present work we have studied influence of the anisotropy on the reconstructed image of a defect generated by an inversion technique. The composite material is considered as homogeneous with macro properties: thermal conductivity K, specific heat c, and density rho. The simulation process involves two sequential steps: solving the three dimensional transient heat diffusion equation for a sample with a defect, then estimating the defect location and size from the surface spatial and temporal thermal distributions (inverse problem), calculated from the simulations.

Plotnikov, Y. A.↗
NASA NTRS1992

Columnar and subsurface silicide growth with novel molecular beam epitaxy techniques

We have found novel growth modes for epitaxial CoSi2 at high temperatures coupled with Si-rich flux ratios or low deposition rates. In the first of these modes, codeposition of metal and Si at 600-800 C with excess Si leads to the formation of epitaxial silicide columns surrounded by single-crystal Si. During the initial stages of the deposition, the excess Si grows homoepitaxially in between the silicide, which forms islands, so that the lateral growth of the islands is confined. Once a template layer is established by this process, columns of silicide form as a result of selective epitaxy of silicide on silicide and Si on Si. This growth process allows nanometer control over silicide particles in three dimensions. In the second of these modes, a columnar silicide seed layer is used as a template to nucleate subsurface growth of CoSi2. With a 100 nm Si layer covering CoSi2 seeds, Co deposited at 800C and 0.01 nm/s diffuses down to grow on the buried seeds rather than nucleating surface silicide islands. For thicker Si caps or higher deposition rates, the surface concentration of Co exceeds the critical concentration for nucleation of islands, preventing this subsurface growth mode from occurring. Using this technique, single-crystal layers of CoSi2 buried under single-crystal Si caps have been grown.

Fathauer, R. W.↗
NASA NTRS1998

Coupling Processes between Atmospheric Chemistry and Climate

This is the fourth semi-annual report for NAS5-97039, covering the time period July through December 1998. The overall objective of this project is to improve the understanding of coupling processes between atmospheric chemistry and climate. Model predictions of the future distributions of trace gases in the atmosphere constitute an important component of the input necessary for quantitative assessments of global change. We will concentrate on the changes in ozone and stratospheric sulfate aerosol, with emphasis on how ozone in the lower stratosphere would respond to natural or anthropogenic changes. The key modeling tools for this work are the Atmospheric and Environmental Research (AER) two-dimensional chemistry-transport model, the AER two-dimensional stratospheric sulfate model, and the AER three-wave interactive model with full chemistry. For this six month period, we report on a modeling study of new rate constant which modify the NOx/NOy ratio in the lower stratosphere; sensitivity to changes in stratospheric water vapor in the future atmosphere; a study of N2O and CH4 observations which has allowed us to adjust diffusion in the 2-D CTM in order to obtain appropriate polar vortex isolation; a study of SF6 and age of air with comparisons of models and measurements; and a report on the Models and Measurements II effort.

Ko, Malcolm K. W.↗
NASA NTRS1993

Ultra low-noise charge coupled device

Special purpose CCD designed for ultra low-noise imaging and spectroscopy applications that require subelectron read noise floors, wherein a non-destructive output circuit operating near its 1/f noise regime is clocked in a special manner to read a single pixel multiple times. Off-chip electronics average the multiple values, reducing the random noise by the square-root of the number of samples taken. Noise floors below 0.5 electrons rms are possible in this manner. In a preferred embodiment of the invention, a three-phase CCD horizontal register is used to bring a pixel charge packet to an input gate adjacent a floating gate amplifier. The charge is then repeatedly clocked back and forth between the input gate and the floating gate. Each time the charge is injected into the potential well of the floating gate, it is sensed non-destructively. The floating gate amplifier is provided with a reference voltage of a fixed value and a pre-charge gate for resetting the amplifier between charge samples to a constant gain. After the charge is repeatedly sampled a selected number of times, it is transferred by means of output gates, back into the horizontal register, where it is clocked in a conventional manner to a diffusion MOSFET amplifier. It can then be either sampled (destructively) one more time or otherwise discarded.

Janesick, James R.↗
NASA NTRS1993

A model of concurrent flow flame spread over a thin solid fuel

A numerical model is developed to examine laminar flame spread and extinction over a thin solid fuel in lowspeed concurrent flows. The model provides a more precise fluid-mechanical description of the flame by incorporating an elliptic treatment of the upstream flame stabilization zone near the fuel burnout point. Parabolic equations are used to treat the downstream flame, which has a higher flow Reynolds number. The parabolic and elliptic regions are coupled smoothly by an appropriate matching of boundary conditions. The solid phase consists of an energy equation with surface radiative loss and a surface pyrolysis relation. Steady spread with constant flame and pyrolysis lengths is found possible for thin fuels and this facilitates the adoption of a moving coordinate system attached to the flame with the flame spread rate being an eigen value. Calculations are performed in purely forced flow in a range of velocities which are lower than those induced in a normal gravity buoyant environment. Both quenching and blowoff extinction are observed. The results show that as flow velocity or oxygen percentage is reduced, the flame spread rate, the pyrolysis length, and the flame length all decrease, as expected. The flame standoff distance from the solid and the reaction zone thickness, however, first increase with decreasing flow velocity, but eventually decrease very near the quenching extinction limit. The short, diffuse flames observed at low flow velocities and oxygen levels are consistent with available experimental data. The maximum flame temperature decreases slowly at first as flow velocity is reduced, then falls more steeply close to the quenching extinction limit. Low velocity quenching occurs as a result of heat loss. At low velocities, surface radiative loss becomes a significant fraction of the total combustion heat release. In addition, the shorter flame length causes an increase in the fraction of conduction downstream compared to conduction to the fuel. These heat losses lead to lower flame temperatures, and ultimately, extinction. This extinction mechanism differs from that of blowoff, where the flame is unable to be stabilized due to the high flow velocity.

Ferkul, Paul V.↗