Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed average tracking”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

On the accuracy of ERS-1 orbit predictions

Since the launch of ERS-1, the D-PAF (German Processing and Archiving Facility) provides regularly orbit predictions for the worldwide SLR (Satellite Laser Ranging) tracking network. The weekly distributed orbital elements are so called tuned IRV's and tuned SAO-elements. The tuning procedure, designed to improve the accuracy of the recovery of the orbit at the stations, is discussed based on numerical results. This shows that tuning of elements is essential for ERS-1 with the currently applied tracking procedures. The orbital elements are updated by daily distributed time bias functions. The generation of the time bias function is explained. Problems and numerical results are presented. The time bias function increases the prediction accuracy considerably. Finally, the quality assessment of ERS-1 orbit predictions is described. The accuracy is compiled for about 250 days since launch. The average accuracy lies in the range of 50-100 ms and has considerably improved.

Koenig, Rolf↗
NASA NTRS2011

Uncertainty Analysis in Space Radiation Protection

Space radiation is comprised of high energy and charge (HZE) nuclei, protons, and secondary radiation including neutrons. The uncertainties in estimating the health risks from galactic cosmic rays (GCR) are a major limitation to the length of space missions, the evaluation of potential risk mitigation approaches, and application of the As Low As Reasonably Achievable (ALARA) principle. For long duration space missio ns, risks may approach radiation exposure limits, therefore the uncertainties in risk projections become a major safety concern and methodologies used for ground-based works are not deemed to be sufficient. NASA limits astronaut exposures to a 3% risk of exposure induced death (REID) and protects against uncertainties in risks projections using an assessment of 95% confidence intervals in the projection model. We discuss NASA s approach to space radiation uncertainty assessments and applications for the International Space Station (ISS) program and design studies of future missions to Mars and other destinations. Several features of NASA s approach will be discussed. Radiation quality descriptions are based on the properties of radiation tracks rather than LET with probability distribution functions (PDF) for uncertainties derived from radiobiology experiments at particle accelerators. The application of age and gender specific models for individual astronauts is described. Because more than 90% of astronauts are never-smokers, an alternative risk calculation for never-smokers is used and will be compared to estimates for an average U.S. population. Because of the high energies of the GCR limits the benefits of shielding and the limited role expected for pharmaceutical countermeasures, uncertainty reduction continues to be the optimal approach to improve radiation safety for space missions.

Cucinotta, Francis A.↗
DOE OSTI2026

Characterization of runaway electron impact on instrumented sacrificial limiters on DIII-D

Instrumented sacrificial limiter heads, both domed (proud) and flat (flush) are used in DIII-D runaway electron (RE) wall strikes to study the wall impact dynamics with high spatial and time resolution. The approximate structure of the RE wetted area and heating depth on the domed limiter heads were predicted qualitatively using orbit-tracking simulations, although a strong left–right asymmetry (about the magnetic field direction) was not captured well by the simulations. It is hypothesized that this difference is perhaps due to the local 3D magnetic field perturbation of the dome limiter head. The average kinetic energy K and pitch angle θ of REs striking the limiter head were estimated from the spatial distribution of local HXR emission and were estimated to be roughly K ≈ 4 MeV and θ ≈ 0.2 rad. These values are roughly consistent with in-plasma values estimated before the loss event, indicating that RE kinetic energy and pitch angle are not drastically altered when transporting to the wall. Large shot–shot variations (1–10 kJ) in energy deposition into the limiter head were observed and were explained by shot–shot variations in locked magneto-hydrodynamics mode toroidal phase. For the largest deposited energies (10 kJ), graphite material failure and explosive dust release was observed, and the depth of material failure at higher energy deposition was successfully reproduced using modelling of volumetric energy deposition and brittle failure. The presence of energetic (keV) level ion impact during the RE wall strike was confirmed by three different surface analysis techniques. The ratio of energetic ion to RE flux appears to be larger on flat surfaces within the RE wetted area, although the energetic ion flux and total energy flux due to energetic ions have not yet been quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

Depths and Geologic Setting of Northern Hemisphere Gullies (and Comparison to Their Southern Counterparts)

Gullies are recognized as having an incised (often sinuous) channel, alcove, and apron morphology as defined. Previous detailed work has focused on the characteristics of southern hemisphere gullies; here we examine the distribution, morphology and depth of northern hemisphere gullies. Gully locations (both hemispheres) were compiled from the literature, not a global search, and thus the gullies studied here are the subset of the whole for which corresponding MOLA tracks are available. We find that the heads of the majority of N. hemisphere gullies commence in a specific cliff-forming layer, and average 250 m below the local surface. This observation suggests that the rock layer has a genetic relationship to gullies and is difficult at present to reconcile with gully formation from snowmelt, which predicts no correlation to rock layers. Gully locations and depths correspond to specific geologic units. We propose that in addition to ground ice and a melting mechanism, gully formation requires groundwater flow along or within specific rock units.

Gilmore, Martha S.↗
DOE OSTI2023

Influence of Methylammonium Chloride on Wide-Bandgap Halide Perovskites Films for Solar Cells

Wide-bandgap perovskites are of paramount importance as the photoactive layer of the top cell in high-efficiency tandem solar cells. Comparably high Br contents are required to widen the perovskite bandgap. However, the increase in Br content causes heterogeneous halide distribution and photoinstability. Here, in this study, the positive effect of the additive methylammonium chloride (MACl) on the optical and electronic properties of Br-rich perovskite, deposited using N-methyl-2-pyrrolidone (NMP) as co-solvent and the gas quenching method, is investigated. Simultaneous in situ grazing-incidence wide-angle X-ray scattering and photoluminescence spectroscopy are used to track the evolution of the structural and optoelectronic properties of the perovskites with different amounts of Br and MACl during the spin-coating and thermal annealing steps. The formation mechanism is elucidated in the presence of MACl. It is observed that chloride ions inhibit the intermediate phases, favoring the formation of a perovskite phase with higher crystallinity. Nano X-ray fluorescence mapping recognizes Br-richer and poorer nanometric domains, whose average sizes reduce for samples with MACl. In conclusion, it is demonstrated that adding MACl affects the formation of wide-bandgap perovskites via destabilization of the intermediate phases and acts on the homogenization of the halide distribution, leading to improved solar cell performances.

14 SOLAR ENERGY↗
DOE OSTI2021

First 3D vector tracking of helium recoils for fast neutron measurements at Super$\mathrm{KEKB}$

Here we present results from the first deployment of novel, high definition, compact gas Time Projection Chambers (TPCs) with pixel chip readout as part of the BEAST II beam background measurement project at SuperKEKB. The TPCs provide detailed 3D imaging of ionization from neutron-induced nuclear recoils in a helium and carbon dioxide target gas mixture at standard temperature and pressure. We present the TPC performance and the neutron backgrounds observed during the initial stage of collider commissioning. We find excellent electron background rejection, leading to background-free nuclear recoil measurements above 50 keV ee , despite the extreme high-background environment. We measure an angular resolution better than 20° for recoil tracks longer than 1.7 mm, corresponding to an average ionization energy of approximately 100 keV ee . We also obtain the full 3D vector direction of helium recoils by utilizing charge profile measurements along the recoil axis, with a correct head/tail assignment efficiency of approximately 80%. With this performance, we present comparisons between measured and simulated event rates, recoil energy spectra, and directional distributions originating from beam–gas and Touschek beam losses at SuperKEKB. We utilize head/tail recognition to distinguish neutron components traveling with positive radial velocity in the Belle II coordinate system from those traveling in the opposite direction. Finally, we present a novel method of discriminating beam–gas interactions from Touschek beam losses that can eliminate the need for dedicated accelerator runs for background measurements. This method is still statistics-limited. However, future studies should be able to verify this method, which in turn could lead to neutron background analysis runs symbiotic with normal Belle II operation. The capabilities demonstrated here also suggest that high definition recoil imaging in gas TPCs is applicable to low energy, low-background experiments, such as directional dark matter searches.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2004

Role of shielding in modulating the effects of solar particle events: Monte Carlo calculation of absorbed dose and DNA complex lesions in different organs

Distributions of absorbed dose and DNA clustered damage yields in various organs and tissues following the October 1989 solar particle event (SPE) were calculated by coupling the FLUKA Monte Carlo transport code with two anthropomorphic phantoms (a mathematical model and a voxel model), with the main aim of quantifying the role of the shielding features in modulating organ doses. The phantoms, which were assumed to be in deep space, were inserted into a shielding box of variable thickness and material and were irradiated with the proton spectra of the October 1989 event. Average numbers of DNA lesions per cell in different organs were calculated by adopting a technique already tested in previous works, consisting of integrating into "condensed-history" Monte Carlo transport codes--such as FLUKA--yields of radiobiological damage, either calculated with "event-by-event" track structure simulations, or taken from experimental works available in the literature. More specifically, the yields of "Complex Lesions" (or "CL", defined and calculated as a clustered DNA damage in a previous work) per unit dose and DNA mass (CL Gy-1 Da-1) due to the various beam components, including those derived from nuclear interactions with the shielding and the human body, were integrated in FLUKA. This provided spatial distributions of CL/cell yields in different organs, as well as distributions of absorbed doses. The contributions of primary protons and secondary hadrons were calculated separately, and the simulations were repeated for values of Al shielding thickness ranging between 1 and 20 g/cm2. Slight differences were found between the two phantom types. Skin and eye lenses were found to receive larger doses with respect to internal organs; however, shielding was more effective for skin and lenses. Secondary particles arising from nuclear interactions were found to have a minor role, although their relative contribution was found to be larger for the Complex Lesions than for the absorbed dose, due to their higher LET and thus higher biological effectiveness. c2004 COSPAR. Published by Elsevier Ltd. All rights reserved.

Non-NASA Center↗
NASA NTRS1992

Cosmic ray LET spectra and doses on board Cosmos-2044 biosatellite

Results of the experiments on board Cosmos-2044 (Biosatellite 9) are presented. Various nuclear track detectors (NTD) (dielectric, AgCl-based, nuclear emulsions) were used to obtain the LET spectra inside and outside the satellite. The spectra from the different NTDs have proved to be in general agreement. The results of LET spectra calculations using two different models are also presented. The resultant LET distributions are used to calculate the absorbed and equivalent doses and the orbit-averaged quality factors (QF) of the cosmic rays (CR). Absorbed dose rates inside (approximately 20 g cm-2 shielding) and outside (1 g cm-2) the spacecraft, omitting electrons, were found to be 4.8 and 8.6 mrad d-1, respectively, while the corresponding equivalent doses were 8.8 and 19.7 mrem d-1. The effects of the flight parameters on the total fluence of, and on the dose from, the CR particles are analyzed. Integral dose distributions of the detected particles are also determined. The LET values which separate absorbed and equivalent doses into 50% intervals are estimated. The CR-39 dielectric NTD is shown to detect 20-30% of the absorbed dose and 60-70% of the equivalent dose in the Cosmos-2044 orbit. The influence of solar activity phase on the magnitude of CR flux is discussed.

Non-NASA Center↗
NASA NTRS1999

Aircraft NO(x) had no Unique Fingerprint on Sonex; Lightning Dominated Fresh NO(x) Sources

Key questions to which SONEX was directed were the following: Can aircraft corridors be detected? Is there a unique tracer for aircraft NO(x)? Can a "background" NO(x) (or NO(y) be defined? What fraction of NO(x) measured during SONEX was from aircraft? How representative was SONEX of the North Atlantic in 1997 and how typical of other years? We attempt to answer these questions through species-species correlations, probability distribution functions (PDFs), and meteorological history. There is not a unique aircraft tracer, largely due to the high variability of air mass origins and tracer ratios, which render "average" quantities meaningless. The greatest NO and NO(y) signals were associated with lightning and convective NO sources. Well-defined background CO, NO(y) and NO(y)/ozone ratio appear in subsets of two cross-track flights with subtropical origins and five flights with predominantly mid-latitude air. Forty percent of the observations on these 7 flights showed NO(y)/ozone to be above background, evidently due to unreacted NO(x). This NO(x) is a combination of aircraft, lightning and surface pollution injected by convection. The strongly subtropical signatures in SONEX observations, confirmed by pv (potential vorticity) values along flight tracks, argues for most of the unreacted NO(x) originating from lightning. Potential vorticity statistics along SONEX flight tracks in 1992-1998, and for the North Atlantic as a whole, show the SONEX meteorological environment to be representative of the North Atlantic flight corridor in the October-November period.

Thompson, A.↗
NASA NTRS1995

Cosmic ray LET spectra and doses on board Cosmos-2044 biosatellite

Results of the experiments on board Cosmos-2044 (Biosatellite 9) are presented. Various nuclear track detectors (NTD) (dielectric, AgCl-based, nuclear emulsions) were used to obtain the Linear Energy Transfer (LET) spectra inside and outside the satellite. The spectra from the different NTDs have proved to be in general agreement. The results of LET spectra calculations using two different models are also presented. The resultant LET distributions are used to calculate the absorbed and equivalent doses and the orbit-averaged quality factors (QF) of the cosmic rays (CR). Absorbed dose rates inside (approximately 20 g cm (exp -2) shielding) and outside (1 g cm(exp -2) the spacecraft, omitting electrons, were found to be 4.8 and 8.6 mrad d (exp -1), respectively, while the corresponding equivalent doses were 8.8 and 19.7 mrem d(exp -1). The effects of the flight parameters on the total fluence of, and on the dose from the CR particles are analyzed. Integral dose distributions of the detected particles are also determined. The LET values which separate absorbed and equivalent doses into 50% intervals are estimated. The CR-39 dielectric NTD is shown to detect 20-30% of the absorbed dose and 60-70% of the equivalent dose in the Cosmos-2044 orbit. The influence of solar activity phase on the magnitude of CR flux is discussed.

Watts, J. W., Jr.↗
NASA NTRS1994

Independent pixel and Monte Carlo estimates of stratocumulus albedo

Monte Carlo radiative transfer methods are employed here to estimate the plane-parallel albedo bias for marine stratocumulus clouds. This is the bias in estimates of the mesoscale-average albedo, which arises from the assumption that cloud liquid water is uniformly distributed. The authors compare such estimates with those based on a more realistic distribution generated from a fractal model of marine stratocumulus clouds belonging to the class of 'bounded cascade' models. In this model the cloud top and base are fixed, so that all variations in cloud shape are ignored. The model generates random variations in liquid water along a single horizontal direction, forming fractal cloud streets while conserving the total liquid water in the cloud field. The model reproduces the mean, variance, and skewness of the vertically integrated cloud liquid water, as well as its observed wavenumber spectrum, which is approximately a power law. The Monte Carlo method keeps track of the three-dimensional paths solar photons take through the cloud field, using a vectorized implementation of a direct technique. The simplifications in the cloud field studied here allow the computations to be accelerated. The Monte Carlo results are compared to those of the independent pixel approximation, which neglects net horizontal photon transport. Differences between the Monte Carlo and independent pixel estimates of the mesoscale-average albedo are on the order of 1% for conservative scattering, while the plane-parallel bias itself is an order of magnitude larger. As cloud absorption increases, the independent pixel approximation agrees even more closely with the Monte Carlo estimates. This result holds for a wide range of sun angles and aspect ratios. Thus, horizontal photon transport can be safely neglected in estimates of the area-average flux for such cloud models. This result relies on the rapid falloff of the wavenumber spectrum of stratocumulus, which ensures that the smaller-scale variability, where the radiative transfer is more three-dimensional, contributes less to the plane-parallel albedo bias than the larger scales, which are more variable. The lack of significant three-dimensional effects also relies on the assumption of a relatively simple geometry. Even with these assumptions, the independent pixel approximation is accurate only for fluxes averaged over large horizontal areas, many photon mean free paths in diameter, and not for local radiance values, which depend strongly on the interaction between neighboring cloud elements.

Cahalan, Robert F.↗
NASA NTRS2010

Eye Tracking Metrics for Workload Estimation in Flight Deck Operation

Flight decks of the future are being enhanced through improved avionics that adapt to both aircraft and operator state. Eye tracking allows for non-invasive analysis of pilot eye movements, from which a set of metrics can be derived to effectively and reliably characterize workload. This research identifies eye tracking metrics that correlate to aircraft automation conditions, and identifies the correlation of pilot workload to the same automation conditions. Saccade length was used as an indirect index of pilot workload: Pilots in the fully automated condition were observed to have on average, larger saccadic movements in contrast to the guidance and manual flight conditions. The data set itself also provides a general model of human eye movement behavior and so ostensibly visual attention distribution in the cockpit for approach to land tasks with various levels of automation, by means of the same metrics used for workload algorithm development.

Ellis, Kyle↗
NASA NTRS1996

Comments on SEE: Comparative Advantages and and Experimental Consequences

The Satellite Energy Exchange experiment measures the periodic, near-miss encounters between a sheppard satellite and a small test body (satellite) in approximately the same orbit about a primary. Several important experimental requirements have been chosen to enhance capabilities: (a) The satellite be flown in a sun-synchronous orbit at an altitude of about 1350 Km, (b) Passive temperature system stabilized by spacecraft axial rotation with sunshade baffles at the end of the spacecraft, (c) Test bodies with different material composition be available for experiments, (d) The containment spacecraft fly about the sheppard mass in a zero-g environment whereas the test bodies, experience average zero-g environment over an orbital period, (e) Primary attitude and station-keeping uses magnetic field alignment plus micro-Newton thrusters such as Field Emission Electric Propulsion, and (f) Very low power (nW) laser tracking systems minimize impulse delivered to test bodies. With the above conditions, SEE has the capabilities: (1) Long duration (several years life-time) flight experiment (2) Long-term, active (with historical time record), self-calibration of satellite mass distribution (capsule geodesy) over lifetime of the spacecraft. (3) Novel passive thermal stabilization systems designed to attain cryogenic temperatures around 78K. (4) Novel spacecraft stabilization systems. (5) Ability to measure G to 1 part in 10(exp 6-7) depending on ultimate duration of experiment. (6) Ability to place limits on both temporal and spacial variations on G. (7) Ability to set experimental limits on the Post Newtonian parameters (PPN) alpha(2) and zeta(2). (8) Ability to measure (or place limits on) the non Einsteinian eccentricity of the Earth-Sun system (and the parameter alpha(1)) for long duration flight. (9) Ability to measure Delta((dot)-G)/G to 1 part in 10(exp 12-13). The MiniSTEP, competes in a limited way with Project SEE. It is designed to improve the measurement of the equivalence principle by seven orders of magnitude using active, low temperature (1.8 K) cooling for SQUID based, differential superconducting circuits. The experiment consists of a small cylinder concentrically located within a larger cylinder at its null gravitational point. The satellite is operated in zero-g mode using four differential accelerometers consisting to two test bodies of different material composition. The SQUIDS are needed to measure test body motion to precisions of 10(exp -18) over a four orbit period. The entire satellite moves in a very precise zero-g mode since the accelerometers are rigidly attached to the satellite. This limits the experiment to an approximately six month due to limitations on helium storage used in cryogenic cooling and thrust control to maintain the zero-g operation.

Smalley, Larry L.↗
DOE OSTI2025

Scaling Arctic landscape and permafrost features improves active layer depth modeling

Tundra ecosystems in the Arctic store up to 40% of global below-ground organic carbon but are exposed to the fastest climate warming on Earth. However, accurately monitoring landscape changes in the Arctic is challenging due to the complex interactions among permafrost, micro-topography, climate, vegetation, and disturbance. This complexity results in high spatiotemporal variability in permafrost distribution and active layer depth (ALD). Moreover, these key tundra processes interact at different scales, and an observational mismatch can limit our understanding of intrinsic connections and dynamics between above and below-ground processes. Consequently, this could limit our ability to model and anticipate how ALD will respond to climate change and disturbances across tundra ecosystems. In this paper, we studied the fine-scale heterogeneity of ALD and its connections with land surface characteristics across spatial and spectral scales using a combination of ground, unoccupied aerial system, airborne, and satellite observations. We showed that airborne sensors such as AVIRIS-NG and medium-resolution satellite Earth observation systems like Sentinel-2 can capture the average ALD at the landscape scale. We found that the best observational scale for ALD modeling is heavily influenced by the vegetation and landform patterns occurring on the landscape. Landscapes characterized by small-scale permafrost features such as polygon tussock tundra require high-resolution observations to capture the intrinsic connections between permafrost and small-scale land surface and disturbance patterns. Conversely, in landscapes dominated by water tracks and shrubs, permafrost features manifest at a larger scale and our model results indicate the best performance at medium resolution (5 m), outperforming both higher (0.4 m) and lower resolution (10 m) models. This transcends our study to show that permafrost response to climate change may vary across dominant ecosystem types, driven by different above- and below-ground connections and the scales at which these connections are happening. We thus recommend tailoring observational scales based on landforms and characteristics for modeling permafrost distribution, thereby mitigating the influences of spatial-scale mismatches and improving the understanding of vegetation and permafrost changes for the Arctic region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Tracking of Marangoni driven motion during laser powder bed fusion

Laser Powder Bed Fusion (LPBF) can be used to create Functionally Graded Materials (FGMs) by selectively distributing property augmenting particles throughout the build matrix. With this method parts can be created with nonhomogeneous properties tuned to spatially varied, pre-determined design constraints. To unlock this capability, the fluid motion that occurs in the molten liquid phase must be thoroughly understood. More importantly, the actual tracking of motion during the LPBF process is needed. Further, a numerical model of heat transfer and fluid flow in LPBF process was created. Using results from the model the motion of particles starting at many different locations was tracked. The average change in position for these particles based on starting location was calculated. The effects of non-dimensional parameters on melt pool size were examined. These results were validated experimentally and compared to what is available in the literature.

42 ENGINEERING↗
DOE OSTI2025

Long-Term Photovoltaic System Performance in Cold, Snowy Climates

As countries around the world transition towards renewable energy, there is increasing interest in using photovoltaic (PV) technologies to help decarbonize northern and alpine communities due to their scalability and affordability. However, a barrier to large-scale adoption of PV in cold climates is long-term performance uncertainty under snowfall, freeze-thaw cycles, low temperatures, and high winds. In this work, we provide a comprehensive review of published silicon degradation rates in cold Koppen-Geiger climate classifications of Dfb (humid continental), Dfc (subarctic), and ET (tundra). We first analyze the system degradation rates of three subarctic ground-mounted photovoltaic sites in North America using the RdTools year-on-year method: an Al-BSF double-axis tracking site in Fairbanks, Alaska (65degrees N); a PERC and silicon heterojunction bifacial vertical and south-tilted site in Fairbanks, Alaska; and a PERC south-facing fixed-tilt site in Fort Simpson, Northwest Territories (62degrees N). Degradation rates of these newly analyzed sites vary between -0.4%/year and -1.5%/year. Combining these data with previously reported cold climate degradation rates, we show that the distribution of cold climate degradation peaks at -0.1%/year to -0.2%/year but has a large tail with rates above -0.5%/year. The average reported cold climate degradation rate is -0.45%/year, whereas the median value is -0.33%/year. These results suggest that despite frequent freeze-thaw cycles and potential exposure to high wind and snow loads, PV systems in cold climates tend to degrade slower than PV systems in warmer climates. The limited sample size of reported degradation rates in cold climates (27) motivates the need for further data acquisition and monitoring efforts as new technologies are deployed.

14 SOLAR ENERGY↗
DOE OSTI2023

Gas-Liquid Flow Modeling for Renewable Fuels Production

Aerobic/anaerobic and gas fermentation pathways have emerged as promising new technologies for the generation of renewable fuels/chemicals from biomass derived sugars, and mixtures of greenhouse/energy rich gas streams (CO2/CH4/H2/CO) via microbial action. Example pathways include sugars-to-ethanol conversion, biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the scale-up, design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of at-scale (~500 m3) bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for gas uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM [1], to perform these simulations, which has been validated against small-scale reactors in our previous work [2]. There is currently a knowledge-gap regarding bubble-size distributions when using gas mixtures with vastly different properties, which can have a significant impact overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a population balance model (PBM) for bubble-size-distribution modeling that is validated against small-scale experiments in our solver with an uncertainty quantification study for bubble coalescence and break-up model parameters. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient are analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors. This work will also examine the performance of different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor mass-transfer coefficient, gas hold-up, and dissolved gas distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented. [1] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998. [2] Rahimi, M., Sitaraman, H., Humbird, D. and Stickel, J., Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution, Chemical Engineering Research and Design, 139: 283-295.

BIOMASS FUELS↗
NASA NTRS2003

Depth and Distribution of CO2 Snow on Mars

The dynamic role of volatiles on the surface of Mars has been a subject of longstanding interest. In the pre-Viking era, much of the debate was necessarily addressed by theoretical considerations. A particularly influential treatment by Leighton and Murray put forth a simple model relying on solar energy balance, and led to the conclusion that the most prominent volatile exchanging with the atmosphere over seasonal cycles is carbon dioxide. Their model suggested that due to this exchange, atmospheric CO2 partial pressure is regulated by polar ice. While current thinking attributes a larger role to H2O ice than did the occasional thin polar coating this model predicted, the CO2 cycle appears to be essentially correct. There are a number of observational constraints on the seasonal exchange of surface volatiles with the atmosphere. The growth and retreat of polar CO2 frost is visible from Earth-based telescopes and from spacecraft in Mars orbit, both at visible wavelengths and in thermal IR properties of the surface. Recently, variations in Gamma ray and neutron fluxes have also been used to infer integrated changes in CO2 mass on the surface. Measurements made by Viking's Mars Atmospheric Water Detector experiment were sensitive to atmospheric H2O vapor abundance. Surface condensates and their transient nature were detected by the Viking landers. The study here is motivated by recent data collected by the Mars Global Surveyor, affording the opportunity to not only detect the lateral distribution of volatiles, but also to constrain the variable volumes of the reservoirs. We elaborate on a technique first employed by Smith et al. By examining averages of a large number of topographic measurements collected by the Mars Orbiter Laser Altimeter (MOLA), that study showed that the zonal pattern of deposition and sublimation of CO2 can be determined. In their first approach, reference surfaces were fit to all measurements in narrow latitude annuli, and the time dependent variations about those mean surfaces were examined. In their second approach, height measurements from pairs of tracks that cross on the surface were interpolated and differenced, forming a set of crossover residuals. These residuals were then examined as a function of time and latitude. The initial studies averaged over longitude to maximize signal and minimize noise in order to isolate the expected small signal. In this follow-up study we now attempt to extract the elevation change pattern also as a function of longitude, and we focus on the crossover approach.

Aharonson, Oded↗