Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient Limiters”

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 793 records · Page 44

Dendritic solidification. III - Some further refinements to the model for dendritic growth under an imposed thermal gradient

Some further refinements to a simple model for dendritic solidification in a binary alloy melt under an imposed positive thermal gradient are presented. Two new expressions for the dendrite tip undercooling have been obtained and shown to yield a limiting value of Delta T sub 0 and very small growth rates. Here Delta T sub 0 is the equilibrium solidification range of the alloy. At very large growth rates, all three tip undercooling expressions reach the same limiting value depending on the value of a dimensionless parameter lambda which is related to the effective diffusion distance ahead of the dendrite tip. The predicted tip undercoolings are, however, somewhat lower at intermediate growth rates. An improved calculation for the solute buildup at the dendrite tip due to curvature effects is also included.

Laxmanan, V.↗

Design and Navier-Stokes analysis of hypersonic wind tunnel nozzles

Four hypersonic wind tunnel nozzles ranging in Mach number from 6 to 17 are designed with the method of characteristics and boundary layer approach (MOC/BL) and analyzed with a Navier-Stokes solver. Limitations of the MOC/BL approach when applied to thick high speed boundary layers with non-zero normal pressure gradients are investigated. Working gases include ideal air, thermally perfect nitrogen and virial CF4. Agreement between the design conditions and Navier-Stokes solutions for ideal air at Mach 6 is good. Thermally perfect nitrogen showed poor agreement at Mach 13.5 and Mach 17. Navier-Stokes solutions for CF4 are not obtained, but comparison of the effects of low gamma to those of high Mach number suggests that the Navier-Stokes solution would not compare well with design.

Benton, James R.↗

ADRPM-VII applied to the long-range acoustic detection problem

An acoustic detection range prediction model (ADRPM-VII) has been written for IBM PC/AT machines running on the MS-DOS operating system. The software allows the user to predict detection distances of ground combat vehicles and their associated targets when they are involved in quasi-military settings. The program can also calculate individual attenuation losses due to spherical spreading, atmospheric absorption, ground reflection and atmospheric refraction due to temperature and wind gradients while varying parameters effecting the source-receiver problem. The purpose here is to examine the strengths and limitations of ADRPM-VII by modeling the losses due to atmospheric refraction and ground absorption, commonly known as excess attenuation, when applied to the long range detection problem for distances greater than 3 kilometers.

Shalis, Edward↗

Low-frequency ionization-driven instability of a discrete auroral arc model

The low-frequency (time scales of tens of seconds) dynamics of the auroral ionosphere differs from that of the nonauroral ionosphere by the presence of strong, unstable space- and time-dependent ionospheric ionization produced by precipitating auroral electrons. If recombination is relatively unimportant (as at high ionospheric heights), it is shown that, in general, transport processes cannot remove this ionization as fast as it is created, and no equilibrium is possible. These nonequilibrium phenomena are investigated in the context of a nonlinear adiabatic auroral model, which has previously been studied in static situations. A linearized local perturbation analysis is given of what amounts to a current-driven E x B gradient-drift instability with an ionization source, as well as some exact nonlinear solutions valid in a finite but limited spatial region. These solutions show continuing motion of auroral potential and plasma density, as the aurora tries to shift its ionization problems from one place to another. The analysis gives clues to the possibility of generation of chaos and of fine-scale spatial structure.

Cornwall, John M.↗

Hypersonic flows as related to the National Aerospace plane

Experimental data for a series of 2-D and 3-D shock wave/boundary layer interaction flows at Mach 8.2 are presented. The test bodies, composed of simple geometric shapes fastened to a flat plate test bed, were designed to generate flows with varying degrees of pressure gradient, boundary layer separation, and turning angle. The data include surface pressure and heat transfer distributions as well as limited mean flowfield surveys both in the undisturbed and interaction regimes. The data are presented in a convenient form to be used to validate existing or future computational models of these hypersonic flows.

Kussoy, Marvin↗

Study of FES/CAST/HGS

The microgravity materials processing program has been instrumental in providing the crystal growth community with an experimental environment to better understand the phenomena associated with the growing of crystals. In many applications one may pursue the growth of large single crystals which cannot be grown on earth due to convective driven flows. A microgravity environment is characterized by neither convection of buoyancy. Consequently superior crystals are able to be grown in space. On the other hand, since neither convection nor buoyancy dominates the fluid flow in a microgravity environment, then lesser dominating phenomena can affect crystal growth, such as surface driven flows or diffusion limited solidification. In the case of experiments that are to be flown in space using the Fluid Experiments System (FES), diffusion limited growth should be the dominating phenomenon. The use of holographic and Schlieren optical techniques for studying the concentration gradients in solidification processes has been used by several investigators over the years. The Holographic Ground System (HGS) facility at MSFC has been a primary resource in researching this capability. Consequently scientific personnel have been able to utilize these techniques in both ground based research and in space experiments. An important event in the scientific utilization of the HGS facilities was the TGS (triglycine sulfate) Crystal Growth and the Casting and Solidification Technology (CAST) experiments that were flown on the International Microgravity Lab (IML) mission in March of this year. The preparation and processing of these space observations are the primary experiments reported in this work. This project provides some ground-based studies to optimize on the holographic techniques used to acquire information about the crystal growth processes flown on IML. Since the ground-based studies will be compared with the space-based experimental results, it is necessary to conduct sufficient ground based studies to best determine how the experiment in space worked. The current capabilities in computer based systems for image processing and numerical computation have certainly assisted in those efforts. As anticipated, this study has certainly shown that these advanced computing capabilities are helpful in the data analysis of such experiments.

Workman, Gary L.↗

An application of a two-equation model of turbulence to three-dimensional chemically reacting flows

A numerical study of three dimensional chemically reacting and non-reacting flowfields is conducted using a two-equation model of turbulence. A generalized flow solver using an implicit Lower-Upper (LU) diagonal decomposition numerical technique and finite-rate chemistry has been coupled with a low-Reynolds number two-equation model of turbulence. This flow solver is then used to study chemically reacting turbulent supersonic flows inside combustors with synergetic fuel injectors. The reacting and non-reacting turbulent combustor solutions obtained are compared with zero-equation turbulence model solutions and with available experimental data. The hydrogen-air chemistry is modeled using a nine-species/eighteen reaction model. A low-Reynolds number k-epsilon model was used to model the effect of turbulence because, in general, the low-Reynolds number k-epsilon models are easier to implement numerically and are far more general than algebraic models. However, low-Reynolds number k-epsilon models require a much finer near-wall grid resolution than high-Reynolds number models to resolve accurately the near-wall physics. This is especially true in complex flowfields, where the stiff nature of the near-wall turbulence must be resolved. Therefore, the limitations imposed by the near-wall characteristics and compressible model corrections need to be evaluated further. The gradient-diffusion hypothesis is used to model the effects of turbulence on the mass diffusion process. The influence of this low-Reynolds number turbulence model on the reacting flowfield predictions was studied parametrically.

Lee, J.↗

Elements of H-mode pedestal structure

Abstract This paper reviews current understanding of key physics elements that control the H-mode pedestal structure, which exists at the boundary of magnetically confined plasmas. The structure of interest is the width, height and gradient of temperature, density and pressure profiles in the pedestal. Emphasis is placed on understanding obtained from combined experimental, theoretical and simulation work and on results observed on multiple machines. Pedestal profiles are determined by the self-consistent interaction of sources, transport and magnetohydrodynamic limits. The heat source is primarily from heat deposited in the core and flowing to the pedestal. This source is computed from modeling of experimental data and is generally well understood. Neutrals at the periphery of the plasma provide the dominant particle source in current machines. This source has a complex spatial structure, is very difficult to measure and is poorly understood. For typical H-mode operation, the achievable pedestal pressure is limited by repetitive, transient magnetohydrodynamic instabilities. First principles models of peeling–ballooning modes are generally able to explain the observed limits. In some regimes, instability occurs below the predicted limits and these remain unexplained. Several mechanisms have been identified as plausible sources of heat transport. These include neoclassical processes for ion heat transport and several turbulent processes, driven by the steep pedestal gradients, as sources of electron and ion heat transport. Reduced models have successfully predicted the pedestal or density at the pedestal top. Firming up understanding of heat and particle transport remains a primary challenge for developing more complete predictive pedestal models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Laser Forming of Sheet Metal for In-Space Manufacturing Applications

As the in-space economy matures, the migration of typically Earth-bound manufacturing methods, such as metal forming and joining, up to space will be required to meet the increasing demand for refueling and supply depots, habitats, laboratories, lunar surface structures, and sites for deep-space exploration shipbuilding. Only in this way can the magnitude of individual structures begin to exceed what can be launched, deployed, docked, or inflated and allow true economies of scale to be unlocked in space. However, typical terrestrial methods of forming and shaping sheet metal require massive machinery with the ability to deflect reaction forces to a large mass in a fixed reference frame (Earth). To both limit up-mass required for manufacture and to mitigate the issues of reaction forces in the relative reference frames of space, a non-contact, laser-based method of forming metal sheet is proposed. NASA Marshall Space Flight Center (MSFC), in cooperation with DARPA and the University of Florida (UF), have begun investigation of non-contact methods of laser forming sheet metal in thermal vacuum (TVAC) conditions to great success. To date, several grades of aluminum, stainless steel, and titanium coupons have been formed via laser-induced internal stresses in these metals at both ambient atmosphere and vacuum conditions. Temperatures in these tests ranged from -120 to 60 C, and though the total energy required to induce strain does vary with workpiece temperature and environmental pressure, the method remains predictable and controllable for the materials investigated. It is believed that this method will revolutionize manufacturing in-space as it not only continues to build upon the capabilities of the laser (sensing, marking, cutting, drilling, joining, and now forming) as an all-in-one tool, but also reduces the need to counteract forming forces as would be required in mechanical bending. In theory, all required stress to induce bending is generated through thermal gradient and strain within the workpiece, meaning that the need for reaction control and fuel expense on-orbit will be limited. Archimedes once stated, “Give me a lever long enough and a fulcrum on which to place it, and I shall move the world”. With a laser as our lever and our workpiece its own fulcrum, we can move the world’s manufacturing into space.

Laser↗

Integral Reactor Containment Condensation Model and Experimental Validation

This NEUP funded project, NEUP 12-3630, is for experimental, numerical and analytical studies on high-pressure steam condensation phenomena in a steel containment vessel connected to a water cooling tank, carried out at Oregon State University (OrSU) and the University of Wisconsin at Madison (UW-Madison). In the three years of investigation duration, following the original proposal, the planned tasks have been completed: (1) Performed a scaling study for the full pressure test facility applicable to the reference design for the condensation heat transfer process during design basis accidents (DBAs), modified the existing test facility to route the steady-state secondary steam flow into the high pressure containment for controllable condensation tests, and extended the operations at negative gage pressure conditions (OrSU). (2) Conducted a series of DBA and quasi-steady experiments using the full pressure test facility to provide a reliable high pressure condensation database (OrSU). (3) Analyzed experimental data and evaluated condensation model for the experimental conditions, and predicted the prototypic containment performance under accidental conditions (UW-Madison). A film flow model was developed for the scaling analysis, and the results suggest that the 1/3 scaled test facility covers large portion of laminar film flow, leading to a lower average heat transfer coefficient comparing to the prototypic value. Although it is conservative in reactor safety analysis, the significant reduction of heat transfer coefficient (50%) could under estimate the prototypic condensation heat transfer rate, resulting in inaccurate prediction of the decay heat removal capability. Further investigation is thus needed to quantify the scaling distortion for safety analysis code validation. Experimental investigations were performed in the existing MASLWR test facility at OrST with minor modifications. A total of 13 containment condensation tests were conducted for pressure ranging from 4 to 21 bar with three different static inventories of non-condensable gas. Condensation and heat transfer rates were evaluated employing several methods, notably from measured temperature gradients in the HTP as well as measured condensate formation rates. A detailed mass and energy accounting was used to assess the various measurement methods and to support simplifying assumptions required for the analysis. Condensation heat fluxes and heat transfer coefficients are calculated and presented as a function of pressure to satisfy the objectives of this investigation. The major conclusions for those tests are summarized below: (1) In the steam blow-down tests, the initial condensation heat transfer process involves the heating-up of the containment heat transfer plate. An inverse heat conduction model was developed to capture the rapid transient transfer characteristics, and the analysis method is applicable to SMR safety analysis. (2) The average condensation heat transfer coefficients for different pressure conditions and non-condensable gas mass fractions were obtained from the integral test facility, through the measurements of the heat conduction rate across the containment heat transfer plate, and from the water condensation rates measurement based on the total energy balance equation. 15 (3) The test results using the measured HTP wall temperatures are considerably lower than popular condensation models would predict mainly due to the side wall conduction effects in the existing MASLWR integral test facility. The data revealed the detailed heat transfer characteristics of the model containment, important to the SMR safety analysis and the validation of associated evaluation model. However this approach, unlike separate effect tests, cannot isolate the condensation heat transfer coefficient over the containment wall, and therefore is not suitable for the assessment of the condensation heat transfer coefficient against system pressure and noncondensable gas mass fraction. (4) The average condensation heat transfer coefficients measured from the water condensation rates through energy balance analysis are appropriate, however, with considerable uncertainties due to the heat loss and temperature distribution on the containment wall. With the consideration of the side wall conduction effects, the results indicate that the measured heat transfer coefficients in the tests is about 20% lower than the prediction of Dehbi’s correlation, mainly due to the side wall conduction effects. The investigation also indicates an increase in the condensation heat transfer coefficient at high containment pressure conditions, but the uncertainties invoked with this method appear to be substantial. (5) Non-condensable gas in the tests has little effects on the condensation heat transfer at high elevation measurement ports. It does affect the bottom measurements near the water level position. The results suggest that the heavier non-condensable gas is accumulated in the lower portion of the containment due to stratification in the narrow containment space. The overall effects of the non-condensable gas on the heat transfer process should thus be negligible for tall containments of narrow condensation spaces in most SMR designs. Therefore, the previous correlations with noncondensable gas effects are not appropriate to those small SMR containments due to the very poor mixing of steam and non-condensable gas. The MELCOR simulation results agree with the experimental data reasonably well. However, it is observed that the MELCOR overpredicts the heat flux for all analyzed tests. The MELCOR predicts that the heat fluxes for CCT’s approximately range from 30 to 45 kW/m2 whereas the experimental data (averaged) ranges from about 25 to 40 kW/m2. This may be due to the limited availability of liquid film models included in MELCOR. Also, it is believed that due to complex test geometry, measured temperature gradients across the heat transfer plate may have been underestimated and thus the heat flux had been underestimated. The MELCOR model predicts a film thickness on the order of 100 microns, which agrees very well with film flow model developed in this study for scaling analysis. However, the expected differences in film thicknesses for near vacuum and near atmospheric test conditions are not significant. Further study on the behavior of condensate film is expected to refine the simulation results. Possible refinements include but are not limited to, the followings: CFD simulation focusing on the liquid film behavior and benchmarking with experimental analyses for simpler geometries. 16 1 INTRODUCTION This NEUP funded project, NEUP 12-3630, is for experimental, numerical and analytical studies on high-pressure steam condensation phenomena in a steel containment vessel connected to a water cooling tank, carried out at Oregon State University (OrSU) and the University of Wisconsin at Madison (UW-Madison). The experimental results are employed to validate the containment condensation model in reactor containment system safety analysis code for integral SMRs. Such a containment condensation model is important to demonstrate the adequate cooling. In the three years of investigation, following the original proposal, the following planned tasks have been completed: (1) Performed a scaling study for the full pressure test facility applicable to the reference design for the condensation heat transfer process during design basis accidents (DBAs), modified the existing test facility to route the steady-state secondary steam flow into the high pressure containment for controllable condensation tests, and extended the operations at negative gage pressure conditions (OrSU). (2) Conducted a series of DBA and quasi-steady experiments using the full pressure test facility to provide a reliable high pressure condensation database (OrSU). (3) Analyzed experimental data and evaluated condensation model for the experimental conditions, and predicted the prototypic containment performance under accidental conditions (UW-Madison). The results are applicable to integral Small Modular Reactor (SMR) designs, including NuScale, mPower, Westinghouse SMR, Holtec-160 and other integral reactors with small containments of relatively high pressures under accidental conditions. Testing has been conducted at the OrSU laboratory in the existing MASLWR (Multi-Application Small Light Water Reactor) integral test facility sponsored by the US Department of Energy. Its highpressure stainless steel containment model (~2 MPa) is scaled to the NuScale SMR currently under development at NuScale Power, Inc.. Minor modifications to the model containment have been made to control the non-condensable gas fraction and to utilize the secondary loop stable steam flow for condensation testing. UW-Madison has developed a containment condensation model, which leveraged previous validated containment heat transfer work carried out at UW-Madison, and extended the range of applicability of the model to integral SMR designs that utilize containment vessels of high heat transfer efficiencies. In this final report, the research background and literature survey are presented in Chapter 2 and 3, respectively. The test facility description and modifications are summarized in Chapter 4, and the scaling analysis is introduced in Chapter 5. The tests description, procedures, and data analysis are presented in Chapter 6, while the numerical modeling is presented in Chapter 7, followed by a conclusion section in Chapter 8.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Inferring morphology and strength of magnetic fields from proton radiographs

Proton radiography is an important diagnostic method for laser plasma experiments and is particularly important in the analysis of magnetized plasmas. The theory of radiographic image analysis has heretofore only permitted somewhat limited analysis of the radiographs of such plasmas. We furnish here a theory that remedies this deficiency. We show that to linear order in magnetic field gradients, proton radiographs are projection images of the MHD current along the proton trajectories. We demonstrate that in the linear regime (i.e., the small image contrast regime), the full structure of the projected perpendicular magnetic field can be reconstructed by solving a steady-state inhomogeneous 2-dimensional diffusion equation sourced by the radiograph fluence contrast data. Here, we explore the validity of the scheme with increasing image contrast, as well as limitations of the inversion method due to the Poisson noise, discretization errors, radiograph edge effects, and obstruction by laser target structures. We also provide a separate analysis that is suited to the inference of isotropic-homogeneous magnetic turbulence spectra. Finally, we discuss extension of these results to the nonlinear regime (i.e., the order unity image contrast regime).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dispersive transport and symmetry of the dispersion tensor in porous media

The macroscopic laws controlling the advection and diffusion of solute at the scale of the porous continuum are derived in a general manner that does not place limitations on the geometry and time evolution of the pore space. Special focus is given to the definition and symmetry of the dispersion tensor that is controlling how a solute plume spreads out. Here, we show that the dispersion tensor is not symmetric and that the asymmetry derives from the advective derivative in the pore-scale advection-diffusion equation. When flow is spatially variable across a voxel, such as in the presence of a permeability gradient, the amount of asymmetry can be large. As first shown by Auriault [J.-L. Auriault Transp. Porous Med. 85, 771 (2010)TPMEEI0169-391310.1007/s11242-010-9591-y] in the limit of low Péclet number, we show that at any Péclet number, the dispersion tensor D ij satisfies the flow-reversal symmetry D ij (+q)=D ji (-q) where q is the mean flow in the voxel under analysis; however, Reynold's number must be sufficiently small that the flow is reversible when the force driving the flow changes sign. We also demonstrate these symmetries using lattice-Boltzmann simulations and discuss some subtle aspects of how to measure the dispersion tensor numerically. Particularly, the numerical experiments demonstrate that the off-diagonal components of the dispersion tensor are antisymmetric which is consistent with the analytical dependence on the average flow gradients that we propose for these off-diagonal components.

36 MATERIALS SCIENCE↗

Experimental high gradient testing of a 17.1 GHz photonic band-gap accelerator structure

We report the design, fabrication and high gradient testing of a 17.1 GHz Photonic Band Gap accelerator structure. Photonic band-gap (PBG) structures are promising candidates for electron accelerators capable of high-gradient operation because they have the inherent damping of high order modes required to avoid beam breakup instabilities. The 17.1 GHz PBG structure tested was a single cell structure composed of a triangular array of round copper rods of radius 1.45 mm spaced by 8.05 mm. The test assembly consisted of the test PBG cell located between conventional (pillbox) input and output cells, with input power of up to 4 MW from a klystron supplied via a TM01 mode launcher. Breakdown at high gradient was observed by diagnostics including reflected power, downstream and upstream current monitors and visible light emission. The testing procedure was first benchmarked with a conventional disc-loaded waveguide structure, which reached a gradient of 87 MV/m at a breakdown probability of 1.19 × 10^-1 per pulse per meter. The PBG structure was tested with 100 ns pulses at power levels of less than 90 MV/m in order to limit the surface temperature rise to 120 K. The PBG structure reached up to 89 MV/m at a breakdown probability of 1.09×10^-1 per pulse per meter. These test results show that a PBG structure can simultaneously operate at high gradients and low breakdown probability, while also providing wakefield damping.

43 PARTICLE ACCELERATORS↗

Declines in ice cover are accompanied by light limitation responses and community change in freshwater diatoms

The rediscovery of diatom blooms embedded within and beneath the Lake Erie ice cover (2007–2012) ignited interest in psychrophilic adaptations and winter limnology. Subsequent studies determined the vital role ice plays in winter diatom ecophysiology as diatoms partition to the underside of ice, thereby fixing their location within the photic zone. Yet, climate change has led to widespread ice decline across the Great Lakes, with Lake Erie presenting a nearly “ice-free” state in several recent winters. It has been hypothesized that the resultant turbid, isothermal water column induces light limitation amongst winter diatoms and thus serves as a competitive disadvantage. To investigate this hypothesis, we conducted a physiochemical and metatranscriptomic survey that spanned spatial, temporal, and climatic gradients of the winter Lake Erie water column (2019–2020). Our results suggest that ice-free conditions decreased planktonic diatom bloom magnitude and altered diatom community composition. Diatoms increased their expression of various photosynthetic genes and iron transporters, which suggests that the diatoms are attempting to increase their quantity of photosystems and light-harvesting components (a well-defined indicator of light limitation). We identified two gene families which serve to increase diatom fitness in the turbid ice-free water column: proton-pumping rhodopsins (a potential second means of light-driven energy acquisition) and fasciclins (a means to “raft” together to increase buoyancy and co-locate to the surface to optimize light acquisition). With large-scale climatic changes already underway, our observations provide insight into how diatoms respond to the dynamic ice conditions of today and shed light on how they will fare in a climatically altered tomorrow.

54 ENVIRONMENTAL SCIENCES↗

Improved integrability and boundedness of solutions to some high-order variational problems

We give a series of results on the improved integrability and boundedness of solutions to several high-order variational problems with strengthened coercivity. In particular, we consider the homogeneous Dirichlet problem on the minimum of integral functionals and study variational inequalities with unilateral and bilateral obstacles and with integral and gradient constraints.

97 MATHEMATICS AND COMPUTING↗

Assembly and Succession of Iron Oxide Microbial Mat Communities in Acidic Geothermal Springs

Biomineralized ferric oxide microbial mats are ubiquitous features on Earth, are common in hot springs of Yellowstone National Park (YNP, WY, USA), and form due to direct interaction between microbial and physicochemical processes. The overall goal of this study was to determine the contribution of different community members to the assembly and succession of acidic high-temperature Fe(III)-oxide mat ecosystems. Spatial and temporal changes in Fe(III)-oxide accretion and the abundance of relevant community members were monitored over 70 days using sterile glass microscope slides incubated in the outflow channels of two acidic geothermal springs (pH = 3-3.5; temperature = 68-75°C) in YNP. Hydrogenobaculum spp. were the most abundant taxon identified during early successional stages (4-40 days), and have been shown to oxidize arsenite, sulfide, and hydrogen coupled to oxygen reduction. Iron-oxidizing populations of Metallosphaera yellowstonensis were detected within 4 days, and reached steady-state levels within 14-30 days, corresponding to visible Fe(III)-oxide accretion. Heterotrophic archaea colonized near 30 days, and emerged as the dominant functional guild after 70 days and in mature Fe(III)-oxide mats (1-2 cm thick). First-order rate constants of Fe(III)-oxide accretion ranged from 0.046 to 0.05 day -1 , and in situ microelectrode measurements showed that the oxidation of Fe(II) is limited by the diffusion of O2 into the Fe(III)-oxide mat. The formation of microterracettes also implicated O2 as a major variable controlling microbial growth and subsequent mat morphology. The assembly and succession of Fe(III)-oxide mat communities follows a repeatable pattern of colonization by lithoautotrophic organisms, and the subsequent growth of diverse organoheterotrophs. The unique geochemical signatures and micromorphology of extant biomineralized Fe(III)-oxide mats are also useful for understanding other Fe(II)-oxidizing systems.

58 GEOSCIENCES↗

Prediction of the Thermal Runaway Limit and Optimal Operation of Heat Transfer-Limited, Fixed-Bed Reactor Systems

Here, we derive a new prediction for thermal runaway starting from the alpha model for fixed-bed reactor systems. This method accounts for thermal resistance internal to the reactor tube and the radial temperature gradients that result. To showcase our method, we compare its predictions to other common criteria for thermal runaway using o-xylene oxidation as the example chemistry. Even in systems where internal heat transfer is negligible, the empirical practical design criterion for thermal runaway is inaccurate. For cases where internal heat transfer is relevant, our runaway limit is more stringent than limits derived from simpler 1-D models. To augment our work, we optimize the product yield with the thermal runaway constraint using orthogonal collocation. Using the alpha model, the results illustrate that the thermal runaway limit can be accurately determined using either numerical or analytical methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Energetic ion distributions on both sides of the Earth's magnetopause

The Active Magnetospheric Particle Tracer Explorer/Charge Composition Explorer (AMPTE/CCE) spacecraft, with an apogee of approximately 8.8R(sub E) and an inclination of approximately 4.3 deg, sampled the outer dayside equatorial magnetosphere for extended time periods and often crossed into the magnetosheath whenever the solar wind pressure was sufficiently high to compress the magnetopause to less than 8.8R(sub E). We have analyzed ion distributions on both sides of the magnetopause in order to investigate any local time and energy dependence, giving information about physical processes at the magnetopause and the bow shock. Particle measurements are from the Charge Energy Mass (CHEM)(1.5 to 300 keV/e) and Medium-Energy Particle Analyzer (MEPA)(10 keV to 2 MeV) instruments. The wide total energy range (1.5 keV to 2 MeV) covered describes the magnetospheric distributions quite well, and for the purpose of this study, adequately describes the high-energy part of the shocked solar wind. Thus both solar wind and magnetospheric components can be recognized in a mixed particle distribution. Case studies of representative magnetopause crossings at dawn, noon, and dusk, as well as a survey of several other crossings, indicated: (1) a local time and energy dependence of magnetosheath spectra at energies greater than or equal to 50 KeV; spectra were harder at the duskside than at the dawnside and also correlated with magnetospheric activity, (2) constantly much higher intensities in the magnetosphere than in the magnetosheath at energies greater than 10 KeV and an earthward gradient in the subsolar magnetosheath. In addition to the steady state magnetosheath population there exists a burst-type component indicative of a magnetospheric source, and most of the time this is recognized as a flux transfer event. Overall, the results about the origin of the greater than or equal to 50 KeV magnetosheath ions are consistent with the continuous leakage of magnetospheric particles across a tangential discontinuity magnetopause, locally distributed according to magnetospheric drift paths. Magnetic reconnection, although present, should not be a dominant source on average, because it is not continuous in time. Fermi acceleration should not be dominant because it predicts the opposite local time asymmetry, and shock drift acceleration should be a minor contributor at E greater than or equal to 50 keV because of upper-energy cutoff limitations. Our observations also indicate a significant magnetospheric contribution to energies as low as approximately 10 KeV, where the magnetosphere-magnetosheath intensity gradient reverses. However, in orderto examine the relative strength and local time distribution of all possible sources at these energies, a detailed analysis is required.

Paschalidis, N. P↗