Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive burn model”

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 73 records · Page 4

MULTI-DIMENSIONAL MODELING OF THE CFR ENGINE FOR THE INVESTIGATION OF SI NATURAL GAS COMBUSTION AND CONTROLLED END-GAS AUTOIGNITION

Engine knock and misfire are barriers to pathways leading to high-efficiency Spark-Ignited (SI) Natural Gas engines. The general tendency to knock is highly dependent on engine operating conditions and the fuel reactivity. The problem is further complicated by low emission limits and the wide range of chemical reactivity in pipeline quality natural gas. Depending on the region and the source of the natural gas, its reactivity, described by its methane number (analogous to the octane number for liquid SI fuels) can span from 65 - 95. In order to realize diesel-like efficiencies, SI natural gas engines must be designed to operate at high BMEP near knock limits over a wide range of fuel reactivity. This requires a deep understanding regarding the combustion-engine interactions pertaining to flame propagation and end-gas autoignition (EGAI). However, EGAI, if controlled, provides an opportunity to increase SI natural gas engine efficiency by increasing combustion rate and the total burned fuel, mitigating the effects of the slow flame speeds of natural gas fuels which generally reduce BMEP and increase unburned hydrocarbon emissions. For this reason, in order to study EGAI phenomenon, the present work highlights multi-dimensional computational fluid dynamics (CFD) models of the Cooperative Fuel Research (CFR) engine. The CFR engine models are used to investigate fuel-engine interactions that lead to EGAI with natural gas, including effects of fuel reactivity, engine operating parameters, and exhaust gas recirculation (EGR). A Three-Pressure Analysis, performed with GT-Power, was used to estimate initial and boundary conditions for the three-dimensional CFD model. CONVERGE CFD v2.4 was used for the three-dimensional CFD modeling where the level set G-Equation model and SAGE detailed chemical kinetics solver were used. An assessment of the different modeling approaches is also provided to evaluate their limitations, advantages and disadvantages, and for which situations they are most applicable. Model validation was performed with experimental data taken with a CFR engine over varying compression ratio, CA50, EGR fraction, and IMEP and shows good agreement in Peak Cylinder Pressure (PCP), PCP crank angle, and the location of the 10%, 50%, and 90% mass fraction burned (CA10, CA50, and CA90, respectively). The models can predict the onset crank angle and pressure rise rate for light, medium, and heavy EGAI under a variety of fuel reactivities and engine operating conditions.

Bestel, Diego↗

Revised (Mixed-Effects) Estimation for Forest Burning Emissions of Gases and Smoke, Fire/Emission Factor Typology, and Potential Remote Sensing Classification of Types for Ozone and Black-Carbon Simulation

We summarize recent progress (a) in correcting biomass burning emissions factors deduced from airborne sampling of forest fire plumes, (b) in understanding the variability in reactivity of the fresh plumes sampled in ARCTAS (2008), DC3 (2012), and SEAC4RS (2013) airborne missions, and (c) in a consequent search for remotely sensed quantities that help classify forest-fire plumes. Particle properties, chemical speciation, and smoke radiative properties are related and mutually informative, as pictures below suggest (slopes of lines of same color are similar). (a) Mixed-effects (random-effects) statistical modeling provides estimates of both emission factors and a reasonable description of carbon-burned simultaneously. Different fire plumes will have very different contributions to volatile organic carbon reactivity; this may help explain differences of free NOx(both gas- and particle-phase), and also of ozone production, that have been noted for forest-fire plumes in California. Our evaluations check or correct emission factors based on sequential measurements (e.g., the Normalized Ratio Enhancement and similar methods). We stress the dangers of methods relying on emission-ratios to CO. (b) This work confirms and extends many reports of great situational variability in emissions factors. VOCs vary in OH reactivity and NOx-binding. Reasons for variability are not only fuel composition, fuel condition, etc., but are confused somewhat by rapid transformation and mixing of emissions. We use "unmixing" (distinct from mixed-effects) statistics and compare briefly to approaches like neural nets. We focus on one particularly intense fire the notorious Yosemite Rim Fire of 2013. In some samples, NOx activity was not so suppressed by binding into nitrates as in other fires. While our fire-typing is evolving and subject to debate, the carbon-burned delta(CO2+CO) estimates that arise from mixed effects models, free of confusion by background-CO2 variation, should provide a solid base for discussion. (c) We report progress using promising links we find between emissions-related "fire types" and promising features deducible from remote observations of plumes, e.g., single scatter albedo, Angstrom exponent of scattering, Angstrom exponent of absorption, (CO column density)/(aerosol optical depth).

Remote sensing↗

Deflagration to Detonation Transition Update: XDDT Code Modularization

A legacy FORTRAN 77 implementation of the Baer–Nunziato two-phase mixture theory for deflagration-to-detonation transition (DDT) in reactive granular materials—hereafter the XDDT (eXplosive DDT) code—has been modularized to Fortran 90 with modular structure, external input files, and adaptive mesh capability. During validation, two code defects were identified and corrected: an inconsistency in the nodal solid pressure evaluation and a nonphysical burn-front tracking criterion. The ignition criterion was also corrected to use the granular surface temperature from the interface heat transfer model, matching the original Baer implementation. An initial attempt to validate against Figure 3 of the original Baer and Nunziato (1986) paper revealed that the code’s detonation velocity on a 201-node mesh (5.5 km/s) was approximately 21% below the expected Chapman–Jouguet value for 70% TMD HMX (∼7 km/s). Validation was redirected to the piston-driven DDT experiments of McAfee et al. (1989), Shot B-9036, for which well-characterized ionization-pin data are available. With the compaction-burn coefficient calibrated to 𝐶 𝛼 = 75, the XDDT code reproduces the DDT transition time to within 0.4% and produces a steady-state detonation velocity within 4% of the McAfee experimental value of 6.36 km/s. The burn model was generalized to support pressure-dependent exponents, enabling application to nitrocellulose-based ball propellants (TS3659) with a cube-root pressure dependence. Validation against the Sandusky/Baer PDC82 piston-impact experiment yielded a reactive wave velocity of 2.3–2.8 km/s, in good agreement with the experimental value of ∼2.2 km/s, and wave coalescence within 5% of the experimental timing. The mathematical model, input parameter requirements, and a roadmap for extending XDDT to PETN with an autocatalytic burn model are presented.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Numerical modeling of plasma assisted deflagration to detonation transition in a microscale channel

Here, this work numerically studies the plasma assisted deflagration to detonation transition (DDT) of H 2 /O 2 mixtures in a microscale channel with detailed chemistry and transport. The results show that the DDT onset time is non-monotonically dependent on the discharge pulse number. The DDT is accelerated with small pulse numbers, whereas retarded with large ones. Two different DDT regimes, respectively at a small and large plasma discharge number, via acoustic choking of the burned gas and plasma-enhanced reactivity gradient without acoustic choking, are observed. Without plasma discharge, pronounced pressure and temperature gradients in front of the flame are generated by acoustic compression after the choking of the burned gas, triggering DDT via autoignition. With small plasma pulse numbers, the plasma-generated species enhance the ignition kinetics and lead to an increased reactivity in the boundary layer. After the choking of the burned gas, the plasma-enhanced reactivity advances the sequence of autoignition near the wall, strengthens ignition-shock wave coupling, and accelerates DDT. However, with a large discharge pulse number, a direct autoignition initiating DDT can occur without the acoustic choking of the burned gas due to the strongly accelerated reactivity and elevated temperature. In this case, DDT onset is retarded because the elevated temperature increases sonic velocity and the increased reactivity accelerates fuel oxidation in front of the flame, decelerating the formation of a leading shock and subsequent pressure buildup ahead of the flame. The present modeling reveals that no matter with or without plasma discharge, DDT is initiated by autoignition in thermal, pressure, and reactivity gradient fields via the Zel'dovich gradient mechanism. The acoustic choking of the burned gas may not be the necessary condition of DDT with strong plasma-enhanced reactivity gradient. This work provides an answer to the experimentally observed non-monotonic DDT onset time by plasma, which provides guidance to control DDT in advanced detonation engines and fire safety of hydrogen-fueled catalytic reactors in microchannels by non-equilibrium plasma discharge.

33 ADVANCED PROPULSION SYSTEMS↗

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion

The classical Landau/Levich models of liquid propellant combustion, which serve as seminal examples of hydrodynamic instability in reactive systems, have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and/or temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity, surface tension and viscosity on the hydrodynamic instability of the propagating liquid/gas interface. In particular, a composite asymptotic expression, spanning three distinguished wavenumber regimes, is derived for both cellular and pulsating hydrodynamic neutral stability boundaries A(sub p)(k), where A(sub p) is the pressure sensitivity of the burning rate and k is the disturbance wavenumber. For the case of cellular (Landau) instability, the results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limiting case of weak gravity, it is shown that cellular hydrodynamic instability in this context is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(l) wavenumber disturbances. It is also demonstrated that, in the large wavenumber regime, surface tension and both liquid and gas viscosity all produce comparable stabilizing effects in the large-wavenumber regime, thereby providing significant modifications to previous analyses of Landau instability in which one or more of these effects were neglected. In contrast, the pulsating hydrodynamic stability boundary is found to be insensitive to gravitational and surface-tension effects, but is more sensitive to the effects of liquid viscosity, which is a significant stabilizing effect for O(l) and higher wavenumbers. Liquid-propellant combustion is predicted to be stable (i.e., steady and planar) only for a range of negative pressure sensitivities that lie between the two types of hydrodynamic stability boundaries.

Margolis, Stephen B.↗

Type Ia supernovae deflagration-to-detonation transition explosions powered by the Zel’dovich reactivity gradient mechanism

ABSTRACT Our aim in this work is to identify and explain the necessary conditions required for an energetic explosion of a Chandrasekhar-mass white dwarf. We construct and analyse weakly compressible turbulence models with nuclear burning effects for carbon/oxygen plasma at a density expected for the deflagration-to-detonation transition (DDT) to occur. We observe the formation of carbon deflagrations and transient carbon detonations at early times. As turbulence becomes increasingly inhomogeneous, sustained carbon detonations are initiated by the Zel’dovich reactivity gradient mechanism. The fuel is suitably preconditioned by the action of compressive turbulent modes with wavelength comparable to the size of resolved turbulent eddies; no acoustic wave is involved in this process. Oxygen detonations are initiated, aided either by reactivity gradients or by collisions of carbon detonations. The observed evolutionary time-scales are found to be sufficiently short for the above process to occur in the expanding, centrally ignited massive white dwarf. The inhomogeneous conditions produced prior to the DDT might be of consequence for the chemical composition of the outer ejecta regions of Type Ia supernovae from the single degenerate channel, and offer the potential for validation of the proposed model.

Brooker, E.↗

Agricultural Fires in the Southeastern U.S. During SEAC4RS: Emissions of Trace Gases and Particles and Evolution of Ozone, Reactive Nitrogen, and Organic Aerosol

Emissions from 15 agricultural fires in the southeastern U.S. were measured from the NASA DC-8 research aircraft during the summer 2013 Studies of Emissions and Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys (SEAC4RS) campaign. This study reports a detailed set of emission factors (EFs) for 25 trace gases and 6 fine particle species. The chemical evolution of the primary emissions in seven plumes was examined in detail for ~1.2 h. A Lagrangian plume cross-section model was used to simulate the evolution of ozone (O3), reactive nitrogen species, and organic aerosol (OA). Observed EFs are generally consistent with previous measurements of crop residue burning, but the fires studied here emitted high amounts of SO2 and fine particles, especially primary OA and chloride. Filter-based measurements of aerosol light absorption implied that brown carbon (BrC) was ubiquitous in the plumes. In aged plumes, rapid production of O3, peroxyacetyl nitrate (PAN), and nitrate was observed with (Delta)O3/(Delta)CO, (Delta)PAN/(Delta)NOy, and (Delta)nitrate/(Delta)NOy reaching approx. 0.1, approx. 0.3, and approx.0.3. For five selected cases, the model reasonably simulated O3 formation but underestimated PAN formation. No significant evolution of OA mass or BrC absorption was observed. However, a consistent increase in oxygen-to-carbon (O/C) ratios of OA indicated that OA oxidation in the agricultural fire plumes was much faster than in urban and forest fire plumes. Finally, total annual SO2, NOx, and CO emissions from agricultural fires in Arkansas, Louisiana, Mississippi, and Missouri were estimated (within a factor of approx. 2) to be equivalent to approx. 2% SO2 from coal combustion and approx. 1% NOx and approx. 9% CO from mobile sources.

OA mass or BrC absorption↗

Modeling Volatility-Based Aerosol Phase State Predictions in the Amazon Rainforest

Organic aerosol (OA) is a complex matrix of various constituents—fresh (primary organic aerosols—POA) and aged via oxidation (secondary organic aerosols—SOA), generated from biogenic, anthropogenic, and biomass burning sources. The viscosity of OA can be critical in influencing new particle formation, reactive uptake processes that impact evaporation-growth kinetics, and the lifetime of particles in the atmosphere. We report this work utilizes a well-defined relationship between volatility and viscosity for pure compounds, which we incorporated within the Weather Research and Forecasting Model coupled to chemistry (WRF-Chem) to simulate the phase state and viscosity of bulk OA during the dry-to-wet transition season (September–October) in the Amazon rainforest during 2014. Our simulations indicate spatial and temporal heterogeneity in aerosol phase state often not captured by global-scale models. We show the strong role of water associated with organic aerosol ($w_s$) as the dominant factor that can be used to quantitatively estimate OA viscosity. Analysis of WRF-Chem simulations across the entire atmospheric column indicates a strong inverse log-linear relationship between $w_s$ and OA viscosity with a correlation coefficient approaching 1, in the background and biomass burning-influenced conditions. At high altitudes where relative humidity (RH) and temperatures are low, our simulations indicate that OA exists in a semisolid-/solid-like phase state, consistent with previous studies. OA hygroscopicity is strongly correlated (ca. –0.8) with OA viscosity at RH ca. 30–50%, but this RH range is found mostly at low OA concentrations and the middle troposphere (ca. 6–10 km altitudes) in our simulated domain. OA hygroscopicity is uncorrelated with viscosity at higher-RH (near surface) and lower-RH (upper troposphere) regimes. At the urban site near surface, where day–night differences in RH are significant, RH is found to drive the phase state. At the background forested site near surface, where day–night RH differences are small, biomass burning-influenced OA is semisolid and a significant OA associated with background conditions is liquid-like. Simulations indicate a long tail of OA viscosity frequency distributions extending in the semisolid/solid regimes over background biogenic-influenced conditions due to the role of low-volatility OA components such as monoterpene oxidation products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Control-oriented core-SOL-divertor model to address integrated burn and divertor control challenges in ITER

The real-time regulation of a burning plasma’s temperature and density, or burn control, will be necessary to produce high fusion power in future tokamaks like ITER. This is made more challenging due to the plasma’s nonlinear characteristics and the interdependence between the core-plasma and edge-plasma regions. For example, a raising plasma temperature leads to increasing reactivity and therefore to more alpha-particle heating, which further increases temperature. Furthermore, a raise of the fusion power increases the heat flow through the scrape-off-layer (SOL), which can compromise the integrity of the divertor without proper safeguards. For control design, a model-based approach is attractive because it can directly incorporate the nonlinear, coupled, burning-plasma dynamics into the design. To facilitate this design approach, a control-oriented core-SOL-divertor (CSD) model is presented in this work. In this CSD model, a core-plasma model captures the nonlinear dynamics of the core’s density and temperature, and a SOL-divertor model defines the plasma conditions at the separatrix and divertor including the heat load on the target plates. The core-plasma and SOL-divertor models are coupled through the exchange of various variables. In particular, the SOL-divertor model yields the separatrix temperature and the influx of recycled particles into the core-plasma. Further, these variables influence the power and particle balances captured by the core-plasma model. In return, the core-plasma model determines the intensity of the heat and particles fluxes across the separatrix, and this outflow strongly impacts the SOL-divertor model. Therefore, the power and density of the core-plasma, which can be readily modulated through external heating systems and pellet injection, can be viewed as control knobs for the SOL-divertor region in addition to the gas puffing. In simulations of the CSD model, it is demonstrated how external actuation can be utilized to meet burn control and divertor control objectives simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion at Normal and Reduced Gravity

The classical Landau/Levich models of liquid-propellant combustion, despite their relative simplicity, serve as seminal examples that correctly describe the onset of hydrodynamic instability in reactive systems. Recently, these two separate models have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed either numerically or analytically in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity and other parameters on the hydrodynamic instability of the propagating liquid/gas interface. In particular, an analytical expression is derived for the neutral stability boundary A(p)(k), where A(p) is the pressure sensitivity of the burning rate and k is the wavenumber of the disturbance. The results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity (both liquid and gas) and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for negative values of A(p). In the limiting case of weak gravity, it is shown that hydrodynamic instability in liquid-propellant combustion is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumber disturbances. it is also demonstrated that, in general, surface tension and the viscosity of both the liquid and gas phases each produce comparable stabilizing effects in the long-wavenumber regime, thereby providing important modifications to previous analyses in which one or more of these effects were neglected.

Margolis, S. B.↗

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion at Normal and Reduced Gravity

The classical Landau/Levich models of liquid-propellant combustion, despite their relative simplicity, serve as seminal examples that correctly describe the onset of hydrodynamic instability in reactive systems. Recently, these two separate models have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed either numerically or analytically in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity and other parameters on the hydrodynamic instability of the propagating liquid/gas interface. In particular, an analytical expression is derived for the neutral stability boundary A(sub p)(k), where A(sub p) is the pressure sensitivity of the burning rate and k is the wavenumber of the disturbance. The results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity (both liquid and gas) and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limiting case of weak gravity, it is shown that hydrodynamic instability in liquid-propellant combustion is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumber disturbances. It is also demonstrated that, in general, surface tension and the viscosity of both the liquid and gas phases each produce comparable stabilizing effects in the large-wavenumber regime, thereby providing important modifications to previous analyses in which one or more of these effects were neglected.

Margolis, Stephen B.↗

Development of an Effective Finite-rate Oxidation Model for NuSil-coated Charred Carbon Preform Ablators

Recently, a detailed effective finite-rate surface chemistry model was developed for the oxidation of FiberForm [1] using the molecular beam experimental data of Poovathingal et al., [2]. FiberForm is the major building block of the thermal protection system (TPS) material Phenolic Impregnated Carbon Ablator (PICA), commonly used by NASA. The surface chemistry model consists of detailed surface reaction mechanisms such as adsorption, desorption, and several types of Langmuir-Hinshelwood (LH) reactions to characterize the oxygen-carbon interactions at the surface. This model provides excellent agreement with the experimental data for oxidation product compositions and corresponding translational energy distributions. Further, an effective oxidation model was constructed that captures the equivalent interaction of oxygen inside the microstructure (including multiple surface collisions) via a modified reactivity for a smooth wall boundary condition. This enables the use of this model directly into Computational Fluid Dynamics (CFD) codes and Material Response (MR) codes to accurately simulate the gas–surface interactions within FiberForm without using the detailed micro-structure. However, this effective model is valid only for the virgin FiberForm. As the TPS material undergoes ablation, the preform carbon burns and turns into char. In addition, PICA is also coated with a protective silicon coating called NuSil for the purpose of mitigating the spread of phenolic dust, and limit contamination during clean room operations. Fig. 1 shows a X-ray microtomography image of a charred NuSil-coated PICA. The NuSil layer introduces species containing silicon into the product mixture. Further, the reactivity of the carbon within the char layer is different from the virgin FiberForm as shown in Fig. 2. Recently a new set of molecular beam experiments were performed on this NuSil-coated TPS material [3] with the same type of oxygen beam used in the previous experiments. Using the latest experimental data, the previously developed effective model will be extended to account for the charred carbon as well as the NuSil coating. Finally, this new effective model with three phases – preform carbon, char, and NuSil; will be compared and validated against the experimental product compositions.

K Swaminathan Gopalan↗

Acetone in theGlobal Troposphere: Its Possible Role as a Global Source of PAN

Oxygenated hydrocarbons are thought to be important components of the atmosphere but, with the exception of formaldehyde, very little about their distribution and fate is known. Aircraft measurements of acetone (CH3COCH3), PAN (CH3CO3NO2) and other organic species (e. g. acetaldehyde, methanol and ethanol) have been performed over the Pacific, the southern Atlantic, and the subarctic atmospheres. Sampled areas extended from 0 to 12 km altitude over latitudes of 70 deg N to 40 deg S. All measurements are based on real time in-situ analysis of cryogenically preconcentrated air samples. Substantial concentrations of these oxygenated species (10-2000 ppt) have been observed at all altitudes and geographical locations in the troposphere. Important sources include, emissions from biomass burning, plant and vegetation, secondary oxidation of primary non-methane hydrocarbons, and man-made emissions. Direct measurements within smoke plumes have been used to estimate the biomass burning source. Photochemistry studies are used to suggest that acetone could provide a major source of peroxyacetyl radicals in the atmosphere and play an important role in sequestering reactive nitrogen. Model calculations show that acetone photolysis contributes significantly to PAN formation in the middle and upper troposphere.

Singh, H. B.↗

Rationalizing Burned Carbon with Carbon Monoxide Exported from South America

We present several estimates cross-checking the fluxes of carbon to the atmosphere from burning, comparing models that are based on simple land-surface parameterizations and atmospheric transport dynamics. Both estimates made by NASA Ames and USP modeling techniques are quite high compared to some detailed satellite/land-use studies of emissions. The flux of carbon liberated to the atmosphere via biomass burning is important for several reasons. This flux is a fundamental statistic for the parameterization of the large-scale flux of gases controlling the reactive greenhouse gases methane and ozone. Similarly, it is central to the estimation of the translocation of nitrogen and pyrodenitrification in the tropics. Thirdly, CO2 emitted from rainforest clearing contributes directly to carbon lost from the rainforest system as it contributes to greenhouse gas forcing. While CO2 from pasturage, agriculture, etc, is considered to be reabsorbed seasonally, and so "off budget" for the carbon cycle, it must also be accounted. CO2 anomalies related to daily weather and interannual climatic variation are strong enough to perturb our scientific perception of long-term carbon storage trends. We compare fluxes deduced from land-use statistics (originally, W.M. Hao) and from satellite hot pixels (A. Setzer) with atmospheric fluxes determined by the mesoscale/continental scale models RAMS and MM5, and point to some new work with highly resolved global models (the NASA Data Assimilation Office's GEOS4). Our simulations are tied to events, so that measured tracers like CO tie the models directly to the burning and meteorology of a specific period. We point out a particular sensitivity in estimates based on CO, and indicate how analysis of CO2 along with other biomass-burning tracers may lead to an improved multi-species estimator of carbon burned.

Chatfield, R.↗

SCALE Modeling of the Sodium Cooled Fast-Spectrum Advanced Burner Test Reactor

This report documents the modeling and simulation of a sodium-cooled fast reactor (SFR) as part of a U.S. Nuclear Regulatory Commission–sponsored project to assess the modeling and simulation capabilities for accident progression, source term, and consequence analysis for advanced reactor technologies with the Oak Ridge National Laboratory code SCALE and the Sandia National Laboratories (SNL) code MELCOR. Based on publicly available benchmark specifications, a fully heterogeneous 3D SCALE model of the 250 MWth Advanced Burner Test Reactor (ABTR) was developed to demonstrate SCALE’s capabilities for full-core reactivity analysis, fuel inventory prediction, and decay heat analysis of an SFR. The benchmark specifications contain modeling details for the ABTR core at the beginning of equilibrium cycle (BOEC) at operating conditions; they were derived from a 2006 preconceptual design report produced by Argonne National Laboratory. The ABTR was designed to demonstrate reactor-based transmutation of transuranics, that is, to “burn” transuranics recovered from light-water reactor (LWR) spent fuel. The ABTR’s fuel is designed to operate in 4 month cycles using uranium/transuranic (U/TRU) metallic fuel, with a TRU content of approximately 20%, at a conversion ratio of approximately 0.6. Various reactivity calculations were performed with SCALE for the ABTR and, where possible, compared with results available in the open literature. Additionally, SCALE was used to perform a full-core depletion calculation over the 4 month cycle to obtain the nuclide inventory at the end of equilibrium cycle (EOEC). These nuclide inventories, decay heat, power profiles, and reactivity feedback coefficients at EOEC represent the initial conditions for analyzing severe accident scenarios with MELCOR.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Numerical calculation of shock-induced initiation of detonations

Results of numerical calculations of the impact of steel cylinders and spheres on the plastic bonded high explosive PBX 9501 are described. The calculations were carried out by a reactive, multicomponent, two dimensional, Eulerian hydrodynamic computer code, 2DE. The 2DE computer code is a finite difference code that uses the donor acceptor cell method to compute mixed cell fluxes. The parameters in the Forest Fire burn model are developed from experiments where the induced shock approximates a plane wave and are applied, in this case, to a situation where the induced shock is a divergent wave with curvature that depends on the size and shape of the projectile. The calculated results are compared with results from experiments involving instrumented mock and live high explosives, with projectiles of varying size, shapes, and velocities.

Cort, G. E.↗

Experimental and Kinetic Modeling Study of 3-Methyl-2-butenol (Prenol) Oxidation

Longer chain alcohols with 4–5 carbon atoms are attractive alternative fuels as they can be derived from biological sources and since their combustion leads to lower exhaust gas levels of NO x and soot compared to commercial fossil fuels. The auto-ignition behavior of fuels that contain both a hydroxyl group and a C=C double bond in their molecular structure is not well established in the literature. Understanding the influence of these functional groups on the ignition behavior of fuels is critical in the development of tailor-made fuels for advanced combustion engines. In this study, ignition delay times of an unsaturated alcohol, 3-methyl-2-butenol (prenol), are measured using a high-pressure shock tube and a rapid compression machine at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 in “air” in the temperature range 600–1400 K. Furthermore, a detailed kinetic model is developed and validated using the new experimental data in this study. In addition, speciation data in a jet-stirred reactor, ignition delay times, and laminar burning velocities available in the literature were also used to validate the new kinetic model. Fuel flux and sensitivity analyses are performed using this new model to determine the important fuel consumption pathways and critical reactions that affect the reactivity of prenol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗