Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non-linear heat conduction”

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.

17 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

LOOPREF: A Fluid Code for the Simulation of Coronal Loops

This report documents the code LOOPREF. LOOPREF is a semi-one dimensional finite element code that is especially well suited to simulate coronal-loop phenomena. It has a full implementation of adaptive mesh refinement (AMR), which is crucial for this type of simulation. The AMR routines are an improved version of AMR1D. LOOPREF's versatility makes is suitable to simulate a wide variety of problems. In addition to efficiently providing very high resolution in rapidly changing regions of the domain, it is equipped to treat loops of variable cross section, any non-linear form of heat conduction, shocks, gravitational effects, and radiative loss.

deFainchtein, Rosalinda↗
NASA NTRS2003

New Propagating Mode near the Superfluid Transition in He-4

We have observed a new temperature-entropy wave that propagates opposite to the direction of a steady heat flux Q when the helium column is heated from above. This propagating mode is due to non-linear thermo-conductance of the helium sample in the self-organized state. Such a mode had been predicted to exist on the self-organized heat transport state for Q less than about 100 nW/sq cm. We confirm that this mode exists in this regime. However, we also observe that it propagates even when the helium is pushed away from the self-organized heat transport state into the normal state.

Sergatskov, Dmitri A.↗
NASA NTRS1998

Thermal Property Parameter Estimation of TPS Materials

Accurate knowledge of the thermophysical properties of TPS (thermal protection system) materials is necessary for pre-flight design and post-flight data analysis. Thermal properties, such as thermal conductivity and the volumetric specific heat, can be estimated from transient temperature measurements using non-linear parameter estimation methods. Property values are derived by minimizing a functional of the differences between measured and calculated temperatures. High temperature thermal response testing of TPS materials is usually done in arc-jet or radiant heating facilities which provide a quasi one-dimensional heating environment. Last year, under the NASA-ASEE-Stanford Fellowship Program, my work focused on developing a radiant heating apparatus. This year, I have worked on increasing the fidelity of the experimental measurements, optimizing the experimental procedures and interpreting the data.

Maddren, Jesse↗
NASA NTRS1998

Frequency-Domain Analysis of Diffusion-Cooled Hot-Electron Bolometer Mixers

A new theoretical model is introduced to describe heterodyne mixer conversion efficiency and noise (from thermal fluctuation effects) in diffusion-cooled superconducting hot-electron bolometers. The model takes into account the non-uniform internal electron temperature distribution generated by Wiedemann-Franz heat conduction, and accepts for input an arbitrary (analytical or experimental) superconducting resistance-versus- temperature curve. A non-linear large-signal solution is solved iteratively to calculate the temperature distribution, and a linear frequency-domain small-signal formulation is used to calculate conversion efficiency and noise. In the small-signal solution the device is discretized into segments, and matrix algebra is used to relate the heating modulation in the segments to temperature and resistance modulations. Matrix expressions are derived that allow single-sideband mixer conversion efficiency and coupled noise power to be directly calculated. The model accounts for self-heating and electrothermal feedback from the surrounding bias circuit.

Skalare, A.↗
NASA NTRS1975

Finite-element solution for a combined radiative-conductive analysis with mixed diffuse-specular surface characteristics

The finite-element method has been applied to solve a combined radiative-conductive heat transfer problem for a large space telescope similar to those used in orbiting satellites. The derivation of the underlying matrices and associated solution algorithm for a 2-D triangular element is presented in detail. The resulting expressions for this triangular element typify such an analysis, which yields constitutive matrices when the heat equation is cast in the matrix form. The relevant matrices include those pertaining to thermal conductance, internal heat generation, radiative exchanges, and all possible external thermal loadings. Emphasis is placed on the treatment of non-linear radiative interchange between surfaces in an enclosure having mixed diffuse-specular surface characteristics. Essential differences in governing equations describing these distinctive surface characteristics are identified. Concluding remarks are drawn from an example simulating a Cassegrainian space telescope.

Lee, H.-P.↗
NASA NTRS2017

Numerical Simulation of Illumination and Thermal Conditions at the Lunar Poles Using LOLA DTMs

We are interested in illumination conditions and the temperature distribution within the upper two meters of regolith near the lunar poles. Here, areas exist receiving almost constant illumination near areas in permanent shadow, which were identified as potential exploration sites for future missions. For our study a numerical simulation of the illumination and thermal environment for lunar near-polar regions is needed. Our study is based on high-resolution, twenty meters per pixel and 400 x 400 km large polar Digital Terrain Models (DTMs), which were derived from Lunar Orbiter Laser Altimeter (LOLA) data. Illumination conditions were simulated by synthetically illuminating the LOLA DTMs using the horizon method considering the Sun as an extended source. We model polar illumination for the central 50 x 50 km subset and use it as an input at each time-step (2 h) to evaluate the heating of the lunar surface and subsequent conduction in the sub-surface. At surface level we balance the incoming insolation with the subsurface conduction and radiation into space, whereas in the sub-surface we consider conduction with an additional constant radiogenic heat source at the bottom of our two-meter layer. Density is modeled as depth-dependent, the specific heat parameter as temperature-dependent and the thermal conductivity as depth- and temperature-dependent. We implemented a fully implicit finite-volume method in space and backward Euler scheme in time to solve the one-dimensional heat equation at each pixel in our 50 x 50 km DTM. Due to the non-linear dependencies of the parameters mentioned above, Newton's method is employed as the non-linear solver together with the Gauss-Seidel method as the iterative linear solver in each Newton iteration. The software is written in OpenCL and runs in parallel on the GPU cores, which allows for fast computation of large areas and long time scales.

Glaser, P.↗
NASA NTRS2004

Improved Ionospheric Electrodynamic Models and Application to Calculating Joule Heating Rates

Improved techniques have been developed for empirical modeling of the high-latitude electric potentials and magnetic field aligned currents (FAC) as a function of the solar wind parameters. The FAC model is constructed using scalar magnetic Euler potentials, and functions as a twin to the electric potential model. The improved models have more accurate field values as well as more accurate boundary locations. Non-linear saturation effects in the solar wind-magnetosphere coupling are also better reproduced. The models are constructed using a hybrid technique, which has spherical harmonic functions only within a small area at the pole. At lower latitudes the potentials are constructed from multiple Fourier series functions of longitude, at discrete latitudinal steps. It is shown that the two models can be used together in order to calculate the total Poynting flux and Joule heating in the ionosphere. An additional model of the ionospheric conductivity is not required in order to obtain the ionospheric currents and Joule heating, as the conductivity variations as a function of the solar inclination are implicitly contained within the FAC model's data. The models outputs are shown for various input conditions, as well as compared with satellite measurements. The calculations of the total Joule heating are compared with results obtained by the inversion of ground-based magnetometer measurements. Like their predecessors, these empirical models should continue to be a useful research and forecast tools.

Weimer, D. R.↗
NASA NTRS2006

AEROELASTIC SIMULATION TOOL FOR INFLATABLE BALLUTE AEROCAPTURE

A multidisciplinary analysis tool is under development for predicting the impact of aeroelastic effects on the functionality of inflatable ballute aeroassist vehicles in both the continuum and rarefied flow regimes. High-fidelity modules for continuum and rarefied aerodynamics, structural dynamics, heat transfer, and computational grid deformation are coupled in an integrated multi-physics, multi-disciplinary computing environment. This flexible and extensible approach allows the integration of state-of-the-art, stand-alone NASA and industry leading continuum and rarefied flow solvers and structural analysis codes into a computing environment in which the modules can run concurrently with synchronized data transfer. Coupled fluid-structure continuum flow demonstrations were conducted on a clamped ballute configuration. The feasibility of implementing a DSMC flow solver in the simulation framework was demonstrated, and loosely coupled rarefied flow aeroelastic demonstrations were performed. A NASA and industry technology survey identified CFD, DSMC and structural analysis codes capable of modeling non-linear shape and material response of thin-film inflated aeroshells. The simulation technology will find direct and immediate applications with NASA and industry in ongoing aerocapture technology development programs.

Liever, P. A.↗
NASA NTRS2017

Performance of an X-Ray Microcalorimeter with a 240 Micron Absorber and a 50 Micron TES Bilayer

We have been developing superconducting transition-edge sensor (TES) microcalorimeters for a variety of potential astrophysics missions, including Athena. The X-ray Integral Field Unit (X-IFU) instrument on this mission requires close-packed pixels on a 0.25 mm pitch, and high quantum efficiency between 0.2 and 12 keV. The traditional approach within our group has been to use square TES bilayers on molybdenum and gold that are between 100 and 140 microns in size, deposited on silicon nitride membranes to provide a weak thermal conductance to a 50 mK heat bath temperature. It has been shown that normal metal stripes on top of the bilayer are needed to keep the unexplained noise at a level consistent with the expected based upon estimates for the non-equilibrium non-linear Johnson noise.In this work we describe a new approach in which we use a square TES bilayer that is 50 microns in size. While the weak link effect is much stronger in this size of TES, we have found that excellent spectral performance can be achieved without the need for any normal metal strips on top of the TES. A spectral performance of 1.58 eV at 6 KeV has been achieved, the best resolution seen in any of our devices with this pixel size. The absence of normal metal stripes has led to more uniform transition shapes, and more reliable excellent spectral performance. The smaller TES size has meant that that the thermal conductance to the heat bath, determined by the perimeter length of the TES and the membrane thickness, is lower than on previous devices, and thus has a lower count rate capability. This is an advantage for low count-rate applications where the slower speed enables easier multiplexing in the read-out, thus potential higher multiplexing factors. In order to recover the higher count rate capabilities, a potential path exits using thicker silicon nitride membranes to increase the thermal conductance to the heat bath.

Miniussi, Antoine R.↗
NASA NTRS1999

Analysis of Unsteady Flow Phenomena: Shock - Vortex and Shock - Boundary Layer Interactions

The interaction of shock waves with vortices has received much attention in the past, mainly because shock-vortex interaction closely models the interaction of a shock wave with the coherent structures of a turbulent flow-field, and is a key feature in the broad-band shock noise for supersonic jets in off-project conditions. Chu and Kovasznay have shown that a weak disturbance in a viscous heat conducting fluid can be decomposed as the sum of three basic modes, namely acoustic, vortical and entropy mode; the interaction of any of these modes with a shock wave gives rise to all three disturbance modes downstream of the shock. The vortical mode is important since it constitutes the basis of the coherent structures that have been observed to dominate turbulence for low- to moderate-flow speed. Hollingsworth et al. have experimentally investigated the interaction of a cylindrical shock-induced starting vortex with a plane normal shock, and have shown that the interaction generates a cylindrical acoustic pulse that exhibits a quadrupolar structure consisting of four alternate compression and expansion regions centered around the transmitted vortex. The investigations of Hollingsworth and Richards have been extended by Dosanjh and Weeks that have analyzed the interaction of a columnar spiral vortex with a normal shock wave, thus obtaining quantitative measurements and confirming the generation of a progressive cylindrical wavefront of alternate compression-expansion nature. Naumann and Hermanns' have experimentally addressed the non-linear aspects of shock-vortex interaction, and have shown that the interaction causes both a diffraction and a reflection of the shock with a pattern consisting of either a regular-or a Mach-reflection depending on the shock and the vortex strengths. An attempt to theoretically explain the production of sound from the shock-vortex interaction was carried out by Ribner. Pao and Salas have numerically studied two-dimensional shock-vortex interactions.

Grasso, Francesco↗
NASA NTRS1993

The importance of properties in modeling

Casting and welding of superalloys, stainless steel and titanium alloys are processes which can be improved through modeling of heat flow, fluid flow, residual stress development, and microstructural evolution. These simulations require inputs of thermophysical data, some of which involves the partially or totally liquid state. In particular, these processes involve melting, flow in the liquid, and solidification. Modeling of such processes can lead to an improved understanding of defects such as shrinkage, inclusions, cracks, incomplete filling (or penetration), macrosegregation, improper grain structure, and deviations from dimensional specifications. Effective modeling can shorten process development time and improve quality. An approach to these problems is to develop efficient models; validate through correlations with thermal, distortion, and microstructural data; run parametric studies; extract knowledge based rules; and apply to adaptive closed loop control systems. With the appropriate pre- and post-processing, such analyses can be made 'user friendly'. This would include graphical user interfaces as well as realistic images and color maps. In such form, these models can be used for sensitivity analyses, which are useful in defining appropriate sensors and in the development of control strategies. Such modeling can be done at several levels, e.g., the MARO level, modeling large scale phenomena such as heat and fluid flow or material deformation; the MICRO level, modeling the development of dendrites, grains or precipitates; or at the NANO level, modeling point defects, dislocations, stacking faults, etc. There are many computational issues associated with these simulations, e.g. computational efficiency and accuracy. In addition, there are many materials issues, not the least of which is the availability of accurate high temperature thermophysical data for complex alloys. This would include latent heat of fusion, temperature dependent heat capacity and thermal conductivity (for liquid and solid), viscosity, surface tension, thermal expansion, mechanical properties, etc. Preliminary data is frequently gathered from the literature; however, this is often not available for modern alloys. If additional data are required, measurements can be used; however, these are costly, time consuming and can be erroneous due to a lack of testing standards or impure materials. Microstructural predictors can be extracted from thermal information, e.g. cooling rate and thermal gradient; the prediction of microstructure is dependent on solidus and liquidus temperature, mushy zone permeability, the solidification curve, volume changes, phase transformations, alloying effects (such as surface tension or viscosity), mold/metal reactions, metal/environment reactions, etc. Defect maps may be needed to predict the onset of shrinkage, hot cracking or 'freckling'. Constants may be needed for stress relaxation, dendrite coarsening, vaporization, etc. Visualization was used as a tool to better comprehend complex data sets associated with the analysis of directional solidification (including crystal growth) and welding. Examples include not only isotherms, but also cooling rate, growth rate and thermal gradient. The latter two are not single valued scalars, but rather time and space dependent vector fields. Efficient models were developed for both casting and welding to predict heat flow and the relationship to dendrite and grain growth. These codes include many of the non-linear effects, e.g. radiation, which dominate these processes. The home-built FDM code(s) were designed to be useful not only to the scientist, but also to the process engineer. Special output can be requested to compare directly to experimental data. Visualization procedures were developed to visualize critical results, e.g. fusion zone width at the surface opposite that where the arc is applied ('penetration'). Both elaborate and simplified distortion analyses were carried out. It is clear that extensive mechanical property data are critical in order to accurately predict residual stress patterns. A scheme is currently being developed to integrate these modeling tools into a set of control algorithms; however, the success of this approach is critically dependent on the availability of accurate high temperature thermophysical data.

Giamei, A. F.↗
NASA NTRS1995

Cloud Ozone Dust Imager (CODI): Investigation and Technical Plan - Volume 1

The Cloud Ozone Dust Imager (CODI) is proposed to investigate the current climatic balance of the Mars atmosphere, with particular emphasis on the important but poorly understood roles which dust and water ice aerosols play in this balance. The large atmospheric heating (20-50 K) resulting from global dust storms around Mars perihelion is well recognized. However, groundbased observations of Mars atmospheric temperatures, water vapor, and clouds since the Viking missions have identified a much colder, cloudier atmosphere around Mars aphelion that may prove as important as global dust storms in determining the interannual and long-term behavior of the Mars climate. The key climate issues CODI is designed to investigate are: 1) the degree to which non-linear interactions between atmospheric dust heating, water vapor saturation, and cloud nucleation influence the seasonal and interannual variability of the Mars atmosphere, and 2) whether the strong orbital forcing of atmospheric dust loading, temperatures and water vapor saturation determines the long-term balance of Mars water, as reflected in the north-south hemispheric asymmetries of atmospheric water vapor and polar water ice abundances. The CODI experiment will measure the daily, seasonal and (potentially) interannual variability of atmospheric dust and cloud opacities, and the key physical properties of these aerosols which determine their role in the climate cycles of Mars. CODI is a small (1.2 kg), fixed pointing camera, in which four wide-angle (+/- 70 deg) lenses illuminate fixed filters and CCD arrays. Simultaneous sky/surface imaging of Mars is obtained at an angular resolution of 0.28 deg/pixel for wavelengths of 255, 336, 502, and 673 nm (similar to Hubble Space Telescope filters). These wavelengths serve to measure atmospheric ozone (255 and 336 nm), discriminate ice and dust aerosols (336 and 673 nm), and construct color images (336, 502, and 673 nm). The CODI images are detected on four 512 x 512 pixel arrays, as partitioned on two 1024 x 1024 CCD's operated in frame transfer mode. The center of the CODI field-of-view is canted 40 deg from the zenith direction to obtain sky brightness measurements and a 20 deg surface field-of-view. Daily image observations will be conducted when the Sun is greater than or equal to 5 deg outside the edge of the CODI field-of-view, and twilight and nighttime imaging will obtained on a weekly basis. The 673 nm channel includes a polarizer wheel to obtain sky/surface polarimetry. A dust cover protects the entire lens assemblies of all four CODI channels. This opaque dust cover, which is normally opened for CODI imaging, includes a small fixed mirror and transparent window positioned above the 673 nm lens, to redirect the 673 nm field-of-view to the surface for descent imaging. Fixed pointing, internal data buffering, low operating power (2-4 W for less than or equal to 30 seconds), selective data transmission, and simple operational characteristics of the CODI experiment place minimum resource and operational demands on the Mars Surveyor 1998 lander. The CODI science goals are optimized for, but not restricted to, a low-latitude landing site (20 deg S-30 deg N). The primary CODI measurement objectives are the opacities, wave forms, particle properties (size, shape, and alignment), and heights of clouds; the opacities, particle properties, and vertical distribution of dust; and the opacity and vertical distribution of ozone. The variability of cloud, ozone, and dust opacities will be determined on diurnal, daily, and seasonal timescales. Wind velocities will be determined from cloud motions and wave characteristics; and the temporal variability of atmospheric water vapor, with limited altitude information, will be inferred from the CODI ozone observations. Secondary measurement objectives include limited descent imaging capability, surface uv-visible photometry and polarimetry, photochemistry, and meteorite infall rates.

Clancy, R. Todd↗
NASA NTRS1992

Bypass transition in compressible boundary layers

Transition to turbulence in aerospace applications usually occurs in a strongly disturbed environment. For instance, the effects of free-stream turbulence, roughness and obstacles in the boundary layer strongly influence transition. Proper understanding of the mechanisms leading to transition is crucial in the design of aircraft wings and gas turbine blades, because lift, drag and heat transfer strongly depend on the state of the boundary layer, laminar or turbulent. Unfortunately, most of the transition research, both theoretical and experimental, has focused on natural transition. Many practical flows, however, defy any theoretical analysis and are extremely difficult to measure. Morkovin introduced in his review paper the concept of bypass transition as those forms of transition which bypass the known mechanisms of linear and non-linear transition theories and are currently not understood by experiments. In an effort to better understand the mechanisms leading to transition in a disturbed environment, experiments are conducted studying simpler cases, viz. the effects of free stream turbulence on transition on a flat plate. It turns out that these experiments are very difficult to conduct, because generation of free stream turbulence with sufficiently high fluctuation levels and reasonable homogeneity is non trivial. For a discussion see Morkovin. Serious problems also appear due to the fact that at high Reynolds numbers the boundary layers are very thin, especially in the nose region of the plate where the transition occurs, which makes the use of very small probes necessary. The effects of free-stream turbulence on transition are the subject of this research and are especially important in a gas turbine environment, where turbulence intensities are measured between 5 and 20 percent, Wang et al. Due to the fact that the Reynolds number for turbine blades is considerably lower than for aircraft wings, generally a larger portion of the blade will be in a laminar transitional state. The effects of large free stream turbulence in compressible boundary layers at Mach numbers are examined both in the subsonic and transonic regime using direct numerical simulations. The flow is computed over a flat plate and curved surface. while many applications operate in the transonic regime. Due the nature of their numerical scheme, a non-conservation formulation of the Navier-Stokes equations, it is a non-trivial extension to compute flow fields in the transonic regime. This project aims at better understanding the effects of large free-stream turbulence in compressible boundary layers at mach number both in the subsonic and transonic regime using direct numerical simulations. The present project aims at computing the flow over a flat plate and curved surface. This research will provide data which can be used to clarify mechanisms leading to transition in an environment with high free stream turbulence. This information is useful for the development of turbulence models, which are of great importance for CFD applications, and are currently unreliable for more complex flows, such as transitional flows.

Vandervegt, J. J.↗
NASA NTRS1996

Non-linear resonance of fluids in a crystal growth cavity

In the microgravity environment, the effect of gravity on fluid motion is much reduced. Hence, secondary effects such as vibrations, jitters, surface tension, capillary effects, and electromagnetic forces become the dominant mechanism of fluid convection. Numerous studies have been conducted to investigate fluid behavior in microgravity with the ultimate goal of developing processes with minimal influence from convection. Industrial applications such as crystal growth from solidification of melt and protein growth for pharmatheutical application are just a few examples of the vast potential benefit that can be reaped from material processing in space. However, a space laboratory is not immune from all undesirable disturbances and it is imperative that such disturbances be well understood, quantifiable, and controlled. Non-uniform and transient accelerations such as vibrations, jitters, and impulsive accelerations exist as a result of crew activities, space vehicle maneuvering, and the operations of on-board equipment. Measurements conducted on-board a U.S. Spacelab showed the existence of vibrations in the frequency range of 1 to 100 Hz with a dominant mode of 17 Hz and harmonics of 54 Hz. The observed vibration is not limited to any coordinate plane but exists in all directions. Similar situation exists on-board the Russian MIR space station. Due to the large structure of its design, the future International Space Station will have its own characteristic vibration spectrum. It is well known that vibration can exert substantial influence on heat and mass transfer processes, thus hindering any attempts to achieve a diffusion-limited process. Experiments on vibration convection for a liquid-filled enclosure under one-g environment showed the existence of different flow regimes as vibration frequency and intensity changes. Results showed the existence of a resonant frequency, near which the enhancement is the strongest, and the existence of a high frequency asymptote. Numerical simulations of vibration convection have been conducted by Yurkov, Fu and Shieh, and by Wang. These analyses considered a two-dimensional air-filled cell under weightlessness condition and showed results similar to those of the experiments. It is not yet known whether resonance convection can be triggered by jitter alone or whether it requires the interaction of jitter with other convective forces in low gravity. An order of magnitude analysis, however, can be used to show the dependence of the resonance frequency on the fluid Prandtl number. Even though the onset of resonance convection may depend on other factors, results indicates that fluids with low Prandtl numbers are more susceptible to resonance than those with high Prandtl numbers. The current study is aimed at gaining additional insights to this problem using germanium as working fluid. Germanium was chosen for this analysis because of its common usage in solidification process and its relatively low Prandtl number (Pr = 0.02).

Wang, Francis C.↗
NASA NTRS2000

Combustion Response of AP Composite Propellants

The Cohen & Strand model for AP composite propellants is applied as boundary conditions, one for AP and one for binder, in solving the heat conduction equation in each to compute linear and non-linear combustion response properties for each and for the aggregate propellant. Iterations couple AP and binder through the quasi-steady flame processes. Illustrative results for linear response functions (pressure-coupled and velocity-coupled) are presented for a monomodal AP propellant varying AP size, pressure and crossflow speed, and examples of non-linear responses to arbitrary waveforms are shown. A quantitative comparison with response function data is limited to one well-characterized research formulation. Mechanistic implications are discussed, including recommendations for future modeling work.

Shusser, Michael↗
NASA NTRS2004

Single Droplet Combustion of Decane in Microgravity: Experiments and Numerical Modeling

This paper presents experimental data on single droplet combustion of decane in microgravity and compares the results to a numerical model. The primary independent experiment variables are the ambient pressure and oxygen mole fraction, pressure, droplet size (over a relatively small range) and ignition energy. The droplet history (D(sup 2) history) is non-linear with the burning rate constant increasing throughout the test. The average burning rate constant, consistent with classical theory, increased with increasing ambient oxygen mole fraction and was nearly independent of pressure, initial droplet size and ignition energy. The flame typically increased in size initially, and then decreased in size, in response to the shrinking droplet. The flame standoff increased linearly for the majority of the droplet lifetime. The flame surrounding the droplet extinguished at a finite droplet size at lower ambient pressures and an oxygen mole fraction of 0.15. The extinction droplet size increased with decreasing pressure. The model is transient and assumes spherical symmetry, constant thermo-physical properties (specific heat, thermal conductivity and species Lewis number) and single step chemistry. The model includes gas-phase radiative loss and a spherically symmetric, transient liquid phase. The model accurately predicts the droplet and flame histories of the experiments. Good agreement requires that the ignition in the experiment be reasonably approximated in the model and that the model accurately predict the pre-ignition vaporization of the droplet. The model does not accurately predict the dependence of extinction droplet diameter on pressure, a result of the simplified chemistry in the model. The transient flame behavior suggests the potential importance of fuel vapor accumulation. The model results, however, show that the fractional mass consumption rate of fuel in the flame relative to fuel vaporized is close to 1.0 for all but the lowest ambient oxygen mole fractions.

Dietrich, D. L.↗