Search NASA⌕ Search

SEARCH · Search NASA

Results for “Step Height”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Validation of turbulence intensity as simulated by the Weather Research and Forecasting model off the US northeast coast

We report turbulence intensity (TI) is often used to quantify the strength of turbulence in wind energy applications and serves as the basis of standards in wind turbine design. Thus, accurately characterizing the spatiotemporal variability in TI should lead to improved predictions of power production. Nevertheless, turbulence measurements over the ocean are far less prevalent than over land due to challenges in instrumental deployment, maintenance, and operation. Atmospheric models such as mesoscale (weather prediction) and large-eddy simulation (LES) models are commonly used in the wind energy industry to assess the spatial variability of a given site. However, the TI derivation from atmospheric models has not been well examined. An algorithm is proposed in this study to realize online calculation of TI in the Weather Research and Forecasting (WRF) model. Simulated TI is divided into two components depending on scale, including sub-grid (parameterized based on turbulence kinetic energy (TKE)) and grid resolved. The sensitivity of sea surface temperature (SST) on simulated TI is also tested. An assessment is performed by using observations collected during a field campaign conducted from February to June 2020 near the Woods Hole Oceanographic Institution Martha's Vineyard Coastal Observatory. Results show that while simulated TKE is generally smaller than the lidar-observed value, wind speed bias is usually small. Overall, this leads to a slight underestimation in sub-grid-scale estimated TI. Improved SST representation subsequently reduces model biases in atmospheric stability as well as wind speed and sub-grid TI near the hub height. Large TI events in conjunction with mesoscale weather systems observed during the studied period pose a challenge to accurately estimating TI from models. Due to notable uncertainty in accurately simulating those events, this suggests summing up sub-grid and resolved TI may not be an ideal solution. Efforts in further improving skills in simulating mesoscale flow and cloud systems are necessary as the next steps.

17 WIND ENERGY↗
DOE OSTI2024

Implementation of Surface Tension on a Reacting Flow Solver, PeleLM: Preprint

In liquid rocket engines, the fuel is supplied to the combustion chamber in the liquid state though injectors. Such fuel undergoes atomization, vaporization, and combustion processes. To design reliable and efficient injectors, it is required to understand the full processes. This research is part of an effort to develop a full atomization-vaporization-combustion solver from first principles. As an initial step to tackle the atomization process, a multiphase flow solver is under development. For the development, a library of the volume of fluid scheme for multiphase, IRL is coupled with a reacting Navier-Stokes equation solver, PeleLM. Furthermore, as the surface tension has considerable effects on spray breakup. surface tension is implemented in the momentum equation using the continuum surface force model and the improved height function technique.

height function↗
NASA NTRS2003

A 2-D Piston Effect Solution for the Relaxation in the MISTE Cell

When a large density stratification is no longer a problem in a microgravity environment, one would like to increase the sample size in order to increase the signal-to-noise ratio for a specific heat measurement. To reduce the equilibration time associated with the large sample size, we designed a cylindrical cell containing a stack of plates that separate the bulk fluid into 60 equally thin layers. To understand the thermal behavior of the whole cell, we analyzed the thermal behavior of a 2-D composite system of a cylindrical near-critical fluid layer in contact with a cylindrical copper plate. In this 2-D analysis, the circumference boundary of the two cylindrical layers is subjected to a step temperature change. The solution of this 2-D composite system includes the piston effect that speeds up the equilibration in the near-critical fluid layer and the pure diffusion in the copper plate. The results of this analysis indicate that the characteristic length for the equilibration of the stacked cell is determined by an effective thickness of a single fluid layer instead of the total height of the cylindrical cell.

Weilert, Mark↗
DOE OSTI2020

The Development of a Generalized Riser Flow Regime Map Based Upon Higher Moment and Chaotic Statistics Using Electrical Capacitance Volume Tomography (ECVT)

Dynamic analyses have been applied to the temporal signals from an Electro Capacitance Volume Tomography instrument located near mid-height on the riser of an industrial-scale cold-flow circulating fluidized bed to characterize gas-solids flow behavior in the riser. Twelve capacitance electrodes surround the cylindrical riser over a height of 1.3 m. The instrument used a neural network deconvolution algorithm to determine the spatially resolved solids fraction recorded at 52 Hz. Experiments were carried out over a range of gas and solids flows in the transport regime using a Geldart Group B bed material, high density polyethylene with mean particle size of 880 μm. The radial solids distribution was found to vary from one-time step to the next between profiles typical of laminar and turbulent flow. The duration of time spent in each of these flow profiles depended upon the operating regime – dilute, core-annular, or fast fluidized bed. The chaotic structure of the temporal data was characterized using the three conventional approaches: the first 4 moments from the distribution of signal in time, system memory parameters from the autocorrelation function and the Hurst exponent, and analysis of the correlationentropy and correlation dimension of the attractor. These signal analysis techniques were used to clearly distinguish differences between different transport operating regimes. Specifically, it was experimentally observed that a riser transitions from core annular flow profile to dilute and dense regimes via increasing the frequency of short term transients to either dilute or dense flow profiles, respectively. A regime map was generated based upon these dynamics using solids flux and gas velocity axes. Fast fluidized, core annular, and dilute each exhibited different degree of dynamic characteristics typical of fluid dominated or particle compromising behavior. It should be noted that the magnitude for the different statistics was in the same range regardless of the regime, it was the radial profile for the statistic that changed and subsequently identified that there was a change in the regime. Finally, a reduced regime map was developed consisting of plotting the gas velocity normalized by the upper transport velocity versus the solids flux normalized by the saturation carrying capacity. The use of this reduced plot allowed the data from widely different conditions to be plotted and compared on the same<p>graph. Note that in many instances, some of the statistics identified the operating point as being in one regime while others indicated that it was in another indicating a transition region between dilute or core annular regimes and between the core annular and fast fluidization regimes. This now provides a tool that can be used to optimize process performance, identify changes in operating states, or replicate process dynamics during process scaling or changing operating parameters. </p>

Breault, Ronald↗
DOE OSTI2022

Characterization of a modified printed optical particle spectrometer for high-frequency and high-precision laboratory and field measurements

Abstract. The printed optical particle spectrometer (POPS) is a lightweight, low-cost instrument for measurements of aerosol number concentrations and size distributions. This work reports on modifications of the Handix Scientific commercial version of the POPS to facilitate its use in multi-instrument aerosol sampling systems. The flow system is modified by replacing the internal pump with a needle valve and a vacuum pump. The instrument is integrated into closed-flow systems by routing the sheath flow from filtered inlet air. A high-precision multichannel analyzer (MCA) card is added to sample the analog pulse signal. The MCA card is polled at 10 Hz frequency using an external data acquisition system and improves upon the count-rate limitation associated with the POPS internal data acquisition system. The times required to change the concentration between 90 % and 10 % and vice versa for a step change in concentration were measured to be 0.17 and 0.41 s at a flow rate of 5 cm3 s−1. This yields a sampling frequency of ∼ 1–2 Hz, below which the amplitude of measured fluctuations is captured with > 70 % efficiency. The modified POPS was integrated into the dual tandem differential mobility analyzer system to explore the coalescence of dimer particles. Results show that the pulse-height response increases upon dimer coalescence. The magnitude of the increase is broadly consistent with the change in light-scattering amplitude predicted by the T-matrix method. It is anticipated that this modified version of the POPS will extend the utilization of the technique for a range of field and laboratory applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Single-Crystal Polycationic Polymers Obtained by Single-Crystal-to-Single-Crystal Photopolymerization

The efficient preparation of single-crystalline ionic polymers and fundamental understanding of their structure–property relationships at the molecular level remains a challenge in chemistry and materials science. Herein, we observe the single-crystal structure of a highly ordered polycationic polymer (polyelectrolyte) and its proton conductivity. The polyelectrolyte single crystals can be prepared on a gram-scale in quantitative yield, by taking advantage of an ultraviolet/sunlight-induced topochemical polymerization, from a tricationic monomer—a self-complementary building block possessing a preorganized conformation. A single-crystal-to-single-crystal photopolymerization was revealed unambiguously by in situ single-crystal X-ray diffraction analysis, which was also employed to follow the progression of molecular structure from the monomer, to a partially polymerized intermediate, and, finally, to the polymer itself. Collinear polymer chains are held together tightly by multiple Coulombic interactions involving counterions to form two-dimensional lamellar sheets (1 nm in height) with sub-nanometer pores (5 Å). The polymer is extremely stable under 254 nm light irradiation and high temperature (above 500 K). The extraordinary mechanical strength and environmental stability—in combination with its impressive proton conductivity (~3 × 10 –4 S cm –1 )—endow the polymer with potential applications as a robust proton-conducting material. By marrying supramolecular chemistry with macromolecular science, the outcome represents a major step toward the controlled synthesis of single-crystalline polyelectrolyte materials with perfect tacticity.

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

Sap Velocity Data for Urban Trees in Chicago, Illinois (2024-2025)

This dataset contains uncorrected sap velocity measurements using the heat ratio method (HRM) collected using ICT International SFM1x sensors at five urban sites in Chicago, Illinois, as part of the DOE CROCUS project. The data includes continuous monitoring of sap velocity from various tree species, including Maples (Acer spp.): Sugar Maple (Acer saccharum), Silver Maple (Acer saccharinum), and Red Maple (Acer rubrum); Oaks (Quercus spp.): Swamp White Oak (Quercus bicolor); American Elm (Ulmus americana); Honey Locust (Gleditsia triacanthos); Cottonwood (Populus deltoides); and Tree of Heaven (Ailanthus altissima) across Chicago State University (CSU), Northeastern Illinois University (NEIU), Northwestern University (NU), University of Illinois Chicago (UIC), and West Woodlawn "Blacks in Green" (BIG). These include both street trees and those in urban park locations. Measurements were collected at 15-20 minute intervals, depending on the sensor, and transmitted via Long Range Wide Area Network (LoRaWAN) protocols. The wireless data was collected by Sage Network (https://sagecontinuum.org/) nodes. The dataset includes sensor ID, Global Positioning System (GPS) coordinates, tree species (common and scientific names), tree identification number, diameter at breast height (DBH in cm), uncorrected sap velocity measurements (cm/hr) from both inner and outer probes, and Sage Node identifiers so the data can be mapped to related variables such as air quality and wind speed that were collected on the Sage nodes. All timestamps are in local Chicago time (CDT/CST). Quality control flags are provided using a 3-bit binary system indicating physical range violations (< -10 or > 60 cm/hr), step spikes (absolute difference > 36 cm/hr), and stuck sensor conditions (> 10 consecutive identical values). These are raw data, not corrected for wood anatomy or species-specific characteristics. Data is provided in comma separated (CSV) format. This dataset is part of a larger collection of CROCUS environmental monitoring data, including linked datasets from Air Quality Transmitter (AQT) sensors, Weather Transmitter (WXT) sensors, and Multi-Function Research LoRaWAN (MFR) Nodes. DOIs for the supporting data are provided as part of this data package.

Chicago↗
NASA NTRS2024

Airborne lidar measurements of atmospheric CO 2 column concentrations to cloud tops made during the 2017 ASCENDS/ABoVE campaign

We measured the column-averaged atmospheric CO 2 mixing ratio (XCO 2 ) to a variety of cloud tops with an airborne pulsed multi-wavelength integrated path differential absorption (IPDA) lidar during NASA's 2017 ASCENDS/ABoVE airborne campaign. Measurements of height-resolved atmospheric backscatter profiles allow this lidar to retrieve XCO 2 to cloud tops, as well as to the ground, with accurate knowledge of the photon path length. We validated these measurements with those from an onboard in situ CO 2 sensor during spiral-down maneuvers. These lidar measurements were 2–3 times better than those from previous airborne campaigns due to our using a wavelength step-locked laser transmitter and a high-efficiency detector for this campaign. Precisions of 0.6 parts per million (ppm) were achieved for 10 s average measurements to mid-level clouds and 0.9 ppm to low-level clouds at the top of the planetary boundary layer. This study demonstrates the lidar's capability to fill in XCO 2 measurement gaps in cloudy regions and to help resolve the vertical and horizontal distributions of atmospheric CO 2 . Future airborne campaigns and spaceborne missions with this capability can be used to improve atmospheric transport modeling, flux estimation and carbon data assimilation.

Jianping Mao↗
NASA NTRS2002

Onboard Autonomous Rock Shape Analysis for Mars Rovers

In the quest to maximize the scientific return of future robotic missions, it is imperative that our rovers be capable of determining the importance of the science they collect so that they may prioritize the acquisition and relay of that data. As an important step in this process, we present an automated technique to allow a rover to classify the shape and other geologic characteristics of rocks from two-dimensional photographic images and three-dimensional stereographically produced data. Experiments were conducted in the Matlab environment using images returned by JPL's Mars Pathfinder mission. Our method begins by first segmenting the rocks from the background using a combination of image intensity and height data. Various metrics are then used to classify the region's sphericity, roundness, and other geometric properties. Preliminary experiments to determine the most useful metrics were conducted by characterizing the two-dimensional rock shape while the three-dimensional shape was later studied with metrics derived from these two-dimensional techniques. Seven measures were developed and implemented. The performance of each measure was characterized by analyzing imagesfrom the Pathfinder mission and ranking the rocks according to the measured properties. Combined, the measures would provide a tool by which an automated rover could discover a greater amount of information about the data it collects, leading to a more productive mission.

image segmentation↗
DOE OSTI2024

Identification of Climatological Representative Days in the Mid-Atlantic for High-Fidelity Offshore Wind Energy Modeling

The goal of reaching 30 GW of offshore wind energy by 2030 becomes more realistic with the continued approval of offshore wind energy areas by the Biden Administration. In the Mid-Atlantic, where wind energy projects are in the most advanced stages of development, there is increased research focus on the eventual interaction of these wind farms. These interactions, in the form of wakes and cluster wakes, or wakes from multiple wind farms, could have detrimental effects on power production and forecastability for downwind wind farms (Pryor et al. 2022, Golbazi et al. 2022, Rosencrans et al. 2023). To help alleviate these issues, numerical simulations in the form of numerical weather prediction (NWP) and large eddy simulations (LES) can provide insight into when cluster wake situations may occur, but running such simulations can be expensive and difficult to run for multiple years. In this study, we leverage and build upon existing techniques in the literature (Fischereit et al. 2022) to identify climatologically representative days for wind energy areas in the Mid-Atlantic where conditions would promote cluster wake situations. We select meteorological variables (wind speed, wind direction, atmospheric stability, boundary-layer height, TKE) critical to understanding wind energy production and wake propagation. We then consider two different NWP datasets of varying spatial and temporal resolution: ERA5 provides data at hourly intervals from 1940 to present at 0.25 deg (31 km) spatial resolution (Hersbach et al. 2020), and the NOW-23 dataset provides data at 5-minute resolution for 21 years at 2-km spatial resolution (Bodini et al. 2020). Our first step is to compare these two datasets for an overlapping 21-year time period. Initial results show that the required number of days to represent the long-term climate increases with each additional variable considered. In their study of the German Bight, Fischereit et al. (2022) found that they could represent the long-term wind and wave climate in a "near-perfect" way with -180 days, by reaching a Perkins Skill Score (PSS) of 0.9; our investigation of the mid-Atlantic wind resource region with ERA5 and NOW-23 data suggests that we will need -100 days to reach a PSS of 0.9. As we expand our parameter space to include multiple variables, the number of required days will likely grow. These results will ultimately be used to select case studies to best represent cluster wake conditions that apply to this region for the lifetime of likely wind farms in this mid-Atlantic region.

clusterwakes↗
NASA NTRS1993

Altitude profile of aerosols on Mars from measurements of its thermal radiation on limb

Measurements of the thermal (range 7 - 13 micrometers) radiation of Mars with the high space resolution (approximately 2 km) were made by the TERMOSKAN experiment of the Phobos mission. Some of the results were published earlier but only the surface radiation was analyzed in detail. However some part of these measurements was made near the limb of the planet. The atmosphere gives an important input here in the planetary thermal radiation. Beyond the limb the atmosphere is the only source. The task of this work is to estimate some characteristics of the atmosphere using brightness profiles of the thermal radiation near the limb. An appropriate model of the temperature profile T(h) is necessary for such an analysis. A set of T(h) models (nominal, maximal and minimal) was defined using various sources including MARSGRAM, Viking-1 lander data, its theoretical considerations and boundary layer models. On the next step the possible input of the atmospheric gaseous emissions (wing of CO2 15 micrometer band) was estimated. It was found that even for the maximal T(h) this input is no more than a few percents of the measured radiation beyond the limb. Consequently the aerosols are responsible for almost all measured emission. The analysis of the observed profile showed that these aerosols have two components: (1) exponential with the scale height about 10 km and (2) some layered structure (two layers with maxima about 23 and 33 km consisted probably of ice).

Moroz, V. I.↗
NASA NTRS1978

Effects of ultrathin oxides in conducting MIS structures on GaAs

Schottky barrier-type GaAs baseline devices (semiconductor surface etched and then immediately metalized) and GaAs conducting metal oxide-semiconductor devices are fabricated and characterized. The baseline surfaces (no purposeful oxide) are prepared by a basic or an acidic etch, while the surface for the MIS devices are prepared by oxidizing after the etch step. The metallizations used are thin-film Au, Ag, Pd, and Al. It is shown that the introduction of purposeful oxide into these Schottky barrier-type structures examined on n-type GaAs modifies the barrier formation, and that thin interfacial layers can modify barrier formation through trapping and perhaps chemical reactions. For Au- and Pd-devices, enhanced photovoltaic performance of the MIS configuration is due to increased barrier height.

Childs, R. B.↗
NASA NTRS1983

An interactive system for compositing digital radar and satellite data

This paper describes an approach for compositing digital radar data and GOES satellite data for meteorological analysis. The processing is performed on a user-oriented image processing system, and is designed to be used in the research mode. It has a capability to construct PPIs and three-dimensional CAPPIs using conventional as well as Doppler data, and to composite other types of data. In the remapping of radar data to satellite coordinates, two steps are necessary. First, PPI or CAPPI images are remapped onto a latitude-longitude projection. Then, the radar data are projected into satellite coordinates. The exact spherical trigonometric equations, and the approximations derived for simplifying the computations are given. The use of these approximations appears justified for most meteorological applications. The largest errors in the remapping procedure result from the satellite viewing angle parallax, which varies according to the cloud top height. The horizontal positional error due to this is of the order of the error in the assumed cloud height in mid-latitudes. Examples of PPI and CAPPI data composited with satellite data are given for Hurricane Frederic on 13 September 1979 and for a squall line on 2 May 1979 in Oklahoma.

Heymsfield, G. M.↗
NASA NTRS2013

Lightning NOx Statistics Derived by NASA Lightning Nitrogen Oxides Model (LNOM) Data Analyses

What is the LNOM? The NASA Marshall Space Flight Center (MSFC) Lightning Nitrogen Oxides Model (LNOM) [Koshak et al., 2009, 2010, 2011; Koshak and Peterson 2011, 2013] analyzes VHF Lightning Mapping Array (LMA) and National Lightning Detection Network(TradeMark) (NLDN) data to estimate the lightning nitrogen oxides (LNOx) produced by individual flashes. Figure 1 provides an overview of LNOM functionality. Benefits of LNOM: (1) Does away with unrealistic "vertical stick" lightning channel models for estimating LNOx; (2) Uses ground-based VHF data that maps out the true channel in space and time to < 100 m accuracy; (3) Therefore, true channel segment height (ambient air density) is used to compute LNOx; (4) True channel length is used! (typically tens of kilometers since channel has many branches and "wiggles"); (5) Distinction between ground and cloud flashes are made; (6) For ground flashes, actual peak current from NLDN used to compute NOx from lightning return stroke; (7) NOx computed for several other lightning discharge processes (based on Cooray et al., 2009 theory): (a) Hot core of stepped leaders and dart leaders, (b) Corona sheath of stepped leader, (c) K-change, (d) Continuing Currents, and (e) M-components; and (8) LNOM statistics (see later) can be used to parameterize LNOx production for regional air quality models (like CMAQ), and for global chemical transport models (like GEOS-Chem).

Koshak, William↗
NASA NTRSDate not supplied

Demonstration of Chromospheric Magnetic Mapping with CLASP2.1

Probing the magnetic nature of the Sun’s chromosphere requires measurement of the polarization profiles of relevant magnetically sensitive spectral lines, many of which are in the ultraviolet spectrum, necessitating observations above the absorbing terrestrial atmosphere. The CLASP series of sounding rocket missions were designed to develop and test a technique for observing the Sun in ultraviolet light, and for quantifying the polarization of that light. By demonstrating successful measurement and interpretation of the polarization in hydrogen Lyman-alpha and the Mg II h and k spectral lines, these missions are crucial steps towards routine quantitative characterization of the local thermal and magnetic conditions in the solar chromosphere. In the most recent observations, CLASP2.1, the spectrograph slit was scanned across an active region plage to acquire a two-dimensional map of Stokes V/I, to demonstrate the ability of UV spectropolarimetry to yield chromospheric magnetic fields over a large area. The technique yields a set of simultaneous line-of-sight magnetograms at multiple heights within the plage atmosphere. By combining the CLASP2.1 measurements with magnetograms from Hinode/SOT or SDO/HMI, a wide range of atmospheric heights are mapped, from the photosphere to the upper chromosphere.

David E McKenzie↗
DOE OSTI2024

Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations

There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Numerical simulations of the convective flame in white dwarfs

A first step toward better understanding of the mechanism driving convective flames in exploding white dwarfs is presented. The propagation of the convective flame is examined using a two-dimensional implicit hydrodynamical code. The large scales of the instability are captured by the grid while the scales that are smaller than the grid resolution are approximated by a mixing-length approximation. It is found that largescale perturbations (of order of the pressure scale height) do grow significantly during the expansion, leading to a very nonspherical burning front. The combustion rate is strongly enhanced (compared to the unperturbed case) during the first second, but later the expansion of the star suppresses the flame speed, leading to only partial incineration of the nuclear fuel. Our results imply that large-scale perturbations by themselves are not enough to explain the mechanism by which convective flames are driven, and a study of the whole spectrum of relevant perturbations is needed. The implications of these preliminary results on future simulations, in the context of current models for Type Ia supernovae, are discussed.

Livne, Eli↗
NASA NTRS2000

Candle Flames in Non-Buoyant Atmospheres

This paper addresses the behavior of a candle flame in a long-duration, quiescent microgravity environment both on the space Shuttle and the Mir Orbiting Station. On the Shuttle, the flames became dim blue after an initial transient where there was significant yellow (presumably soot) in the flame. The flame lifetimes were typically less than 60 seconds. The safety-mandated candlebox that contained the candle flame inhibited oxygen transport to the flame and thus limited the flame lifetime. The flames on the Mir were similar, except that the yellow luminosity persisted longer into the flame lifetime because of a higher initial oxygen concentration, The Mir flames burned for as long as 45 minutes. The difference in the flame lifetime between the Shuttle and Mir flames was primarily the redesigned candlebox that did not inhibit oxygen transport to the flame. In both environments, the flame intensity and the height-to-width ratio gradually decreased as the ambient oxygen content in the sealed chamber slowly decreased. Both sets of experiments showed spontaneous, axisymmetric flame oscillations just prior to extinction. The paper also presents a numerical model of a candle flame. The formulation is two-dimensional and time-dependent in the gas phase with constant specific heats, thermal conductivity and Lewis number (although different species can have different Lewis numbers), one-step finite-rate kinetics, and gas-phase radiative losses from CO2 and H2O. The treatment of the liquid/wick phase assumes that the, fuel evaporates from a constant diameter sphere connected to an inert cone. The model predicts a steady flame with a shape and size quantitatively similar to the Shuttle and Mir flames. The computation predicts that the flame size will increase slightly with increasing ambient oxygen mole fraction. The model also predicts pre-extinction flame oscillations if the rate of decrease in ambient oxygen is small enough, such as that which would occur for a flame burning in a sealed ambient.

Dietrich, D. L.↗