Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Source Term Analysis”

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 1,153 records · Page 64

Long-Term Simulation of Dust Distribution with the GOCART Model: Correlation with the North Atlantic Oscillation

Global distribution of aeolian dust is simulated from 1981 to 1996 with the Goddard Ozone Chemistry Aerosol Radiation and Transport (GOCART) model. The results are assessed with in-situ measurements and the Total Ozone Mapping Spectrometer (TOMS) aerosol products. The annual budget over the different continents and oceans are analyzed. It is found that there is a maximum of 25% difference of global annual emission from the minimum in 1996 to the maximum in 1988. There is a downward trend of dust emission over Africa and East Asia, of 6 and 2 Tg/yr, respectively. The inter-annual variability of dust distribution is analyzed over the North Atlantic and Africa. It is found that in winter most of the North Atlantic and Africa dust loading is correlated with the North Atlantic Oscillation. The GOCART model indicates that a controlling factor of such correlation can be attributed to dust emission from the Sahel. The Bodele depression is the major dust source in winter and its inter-annual variability is highly correlated with the NAO. However, it is not possible to conclude without further analysis that the North Atlantic Oscillation is forcing the inter-annual variability of dust emission and in-turn dust concentration over the North Atlantic.

Ginoux, P.↗

Hydraulic characterization of a fault zone from fracture distribution

A quantitative assessment of how faults control the migration of geofluids is critical in many areas of geosciences. Here we integrated geological fieldwork, quantitative analysis of the fractures distribution and numerical modeling to build a geometrical representation of a fault zone and to characterize its hydraulic properties. Our target is a fault located in the Majella Mountain (Italy). We collected 21 scan lines across the fault profile in order to characterize its architecture. The numerical modeling of the fracture network of the damage zones and their hydraulic parameters was performed using both commercial (Move ® ) and open source software (dfnWorks and PFLOTRAN). Move ® was used to build a representative model of the fault zone using fracture spacing as a proxy, and to model the hydraulic parameters of the different fault domains. dfnWorks and PFLOTRAN were employed to infer the hydraulic parameters of the damage zones of the fault and then upscale these properties to an equivalent continuum domain, suitable for fluid flow simulations through the whole fault zone. Our findings show how even in a relatively small area it is possible to describe changes in terms of hydraulic properties of a fault zone and to build models capable to represent these variations.

58 GEOSCIENCES↗

Combined Meteorological and Hydrologic Uncertainties Shape Projections of Future Soil Moisture in the Eastern United States

Physical hazards pose risks to many critical systems. Designing adaptive measures to mitigate these risks is challenging due to large uncertainties in modeling future hazards and the associated sectoral responses. Here, we help address this challenge in a hydrologic context by examining the combined role of meteorological forcing and hydrologic parameter uncertainties in shaping projections of future soil moisture. By encoding a simple conceptual water balance model in a differentiable programming framework, we facilitate fast runtimes and an efficient calibration, enabling an improved uncertainty analysis. We characterize uncertainty in model parameters by calibrating against different target data sets and by using several loss functions. We then convolve the resulting parameter ensemble with a set of Earth system model projections to produce a large ensemble (2,340 members) of daily soil moisture simulations. Focusing on the eastern United States, we find that most ensemble members project a drying of soils across the region, although some simulate wetter conditions throughout this century. Our ensemble shows an increase in the frequency and intensity of dry extremes while there is less agreement for wet extremes. We conduct sensitivity analyses on several soil moisture signatures to measure the relative influence of meteorological and hydrologic uncertainties across space and time. Both meteorological and hydrologic factors contribute consistently to uncertainty surrounding long-term trends, while changes to both wet and dry soil extremes are typically more sensitive to hydrologic parameter uncertainty. Our results underscore the need to account for varied sources of uncertainty when developing long-term hydrometeorological projections.

Lafferty, David C. [University of Illinois Urbana‐↗

Interannual Oscillations in Earth Rotation

At seasonal and shorter time scales, polar motion and length‐of‐day (LOD) change are dominated by mass redistribution in the Earth climate system, including the atmosphere, oceans, and hydrosphere. Long‐term polar motion and LOD change are believed to be forced by solid Earth geophysical process, such as glacial isostatic adjustment. Decadal LOD variations are also connected to core‐mantle coupling. This study provides a comprehensive analysis of interannual oscillations in all three components of Earth rotation using accurately measured polar motion and LOD time series, and geophysical excitations computed from atmospheric, oceanic, and hydrological models over the period 1962 to 2018. Strong interannual oscillations remain in all three components after detrending and band‐pass filtering to retain variations at periods between 2 and 8 years. Polar motion Y is dominated by a 3.65‐year oscillation, plus others at periods of 2.5 and 5.9 years. Polar motion X shows a strong oscillation near 5.9 years, plus smaller variations at shorter periods. A 5.9‐year LOD variation has been recognized in earlier studies, but oscillations at periods of 2.36, 3.65, and 4.9 years are also observed. Atmosphere, ocean, and hydrosphere (AOH) sources largely account for 2.5‐ and 3.65‐year components in polar motion Y, but not for the 5.9‐year component in polar motion X and Y. The 2.36‐, 3.65‐, and 4.9‐year components in LOD can largely be accounted for by AOH sources. The 5.9‐year LOD component becomes more prominent after AOH sources have been subtracted. The presence of a 5.9‐year variation in all three components of Earth rotation suggests a source in Earth's interior, probably from the core.

Jianli Chen↗

Baseline Vector Repeatability at the Sub-Millimeter Level Enabled by Radio Interferometer Phase Delays of Intra-Site Baselines

We report the results of position ties for short baselines at eight geodetic sites based on phase delays that are extracted from global geodetic very-long-baseline interferometry (VLBI) observations rather than dedicated short-baseline experiments. An analysis of phase delay observables at X band from two antennas at the Geodetic Observatory Wettzell, Germany, extracted from 107 global 24-hr VLBI sessions since 2019 yields weighted root-mean-square scatters about the mean baseline vector of 0.3, 0.3, and 0.8 mm in the east, north, and up directions, respectively. Position ties are also obtained for other short baselines between legacy antennas and nearby, newly built antennas. They are critical for maintaining a consistent continuation of the realization of the terrestrial reference frame, especially when including the new VGOS network. The phase delays of the baseline WETTZ13N–WETTZELL enable an investigation of sources of error at the sub-millimeter level. We found that a systematic variation of larger than 1 mm can be introduced to the Up estimates of this baseline vector when atmospheric delays were estimated. Although the sub-millimeter repeatability has been achieved for the baseline vector WETTZ13N–WETTZELL, we conclude that long term monitoring should be conducted for more short baselines to assess the instrumental effects, in particular the systematic differences between phase delays and group delays, and to find common solutions for reducing them. This will be an important step toward the goal of global geodesy at the 1 mm level.

geodetic VLBI↗

Towards a Robust Sampling Approach: A Computational Review and Design

As pointed out in several other works, the estimation of the reliability of the electrical grid can not be conducted without the estimation of the stochastic phenomena of electicity demand and electricity production by variable renewable sources. Therefore, sampling procedures have become integral in the design of engineering structures and analysis. Commonly referred to as Monte Carlo uncertainty analysis or integration, the general objective of these procedures is to establish specifics about the uncertainty of an output characteristic of such an engineering system, given uncertainty about its input characteristics. These sampling procedures are applied in a context in which establishing such specifics cannot be performed through other means. Variance-reduction techniques are designed to lessen the variability among estimators to estimate statistics of those output uncertainties. Importance-sampling techniques, on the other hand, are designed to reduce the number of samples needed to estimate a particular statistic—e.g., a tail probability. The combination of these approaches can reduce the computational burden considerably for a particular estimator and statistic. Importance sampling—geared and designed as it is toward improving a particular estimator—suffers, unfortunately, from the unintended consequence of reducing the performance of other estimators in terms of their variance. The objective of this paper is to offer an alternative sampling procedure where this variance does not grow unacceptably large for a suite of estimators. Moreover, it is anticipated that, with additional knowledge of how an engineering output characteristic responds to its input characteristics, tuning parameters of the input’s sampling procedure can be set to improve the output characteristic’s estimation. This work proves the effectiveness of the suggested alternative approach. Such positive outcome will lead to a decrease of the computational burden of performing stochastic optimization of integrated energy systems (e.g., components dispatch, and portfolio composition). In particular, capturing the contribution to the overall system cost of rare and unlikely events and patterns of the electricity demand and production will become less computationally expensive. This is due to the fact that the approach demonstrated here will allow the sampling of those rare occurrences more frequently without misrepresenting their probabilistic impacts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Modeling of the Synfuel Production Process: Process models of Fischer-Tropsch production with electricity and hydrogen provided by various scales of nuclear plants

Synthetic fuels (synfuels), also known as electro-fuels (E-fuels), are hydrocarbon fuels produced from waste CO2 streams and water electrolysis, with electricity as the primary source of energy. To achieve substantial reductions in greenhouse gas (GHG) emissions, electricity sources must release zero carbon or near-zero carbon, as is the case with solar, wind, hydro, and nuclear power. Nuclear power is one of the largest and steadiest domestic sources of clean energy in the United States. Moreover, nuclear power has the potential to produce hydrogen economically for less than $2/kg, reaching the DOE near-term target price. Thus, using nuclear power to produce synfuels has the unique potential to significantly reduce the GHG emissions of hydrocarbon fuels production and end-use applications. Fisher-Tropsch or FT fuel (a mixture of naphtha, jet fuel, and diesel) is of great interest because it is a drop-in fuel that can be blended with conventional petroleum counterparts and is compatible with existing infrastructure. By using the ASPEN Plus model, this report develops FT fuel production models on three scales, corresponding to nuclear plants with capacities of 1000 MWe, 437 Mwe, and 100 MWe, respectively. The FT model case with energy from a 437-MWe nuclear plant is used as a baseline case. This report summarizes the baseline ASPEN Plus model results with a detailed mass and energy analysis. Our modeled facility produces 507 MT/day (185,000 gal/day) of FT fuel by converting 255 MT/day of hydrogen and 1,580 MT/day of CO2. The FT fuel production energy efficiency from hydrogen and electricity energy inputs is 70% (lowerheating-value or LHV-based). Including the high-temperature electrolyzer in the system boundary, the FT fuel production LHV efficiency from electricity and thermal energy inputs is 51%, considering 39.8 kWh/kg of electricity and 6.86 kWh/kg of thermal energy use from a nuclear plant for hydrogen production. The FT production efficiency can potentially be increased by further integrating the heat exchange between nuclear plant and FT process, and this study is underway. The carbon conversion ratio in the baseline case is 99%, with process CO2 capture and recirculation and oxy-combustion using the oxygen by-product from water electrolysis. The hydrogen consumption is 1.38 kg/gal-FT fuel and the CO2 consumption is 8.56 kg/gal-FT fuel in the baseline case. With different FT production scales determined by the nuclear plant capacity, the FT model was scaled using the same operating parameters, which led to the same conversion efficiency regardless of scale. However, the different FT plant scales will impact the economics of FT fuel production; this impact will be examined in the next phase of this study.

Zang, Guiyan↗

SoDaH: the SOils DAta Harmonization database, an open-source synthesis of soil data from research networks, version 1.0

Data collected from research networks present opportunities to test theories and develop models about factors responsible for the long-term persistence and vulnerability of soil organic matter (SOM). Synthesizing datasets collected by different research networks presents opportunities to expand the ecological gradients and scientific breadth of information available for inquiry. Synthesizing these data is challenging, especially considering the legacy of soil data that have already been collected and an expansion of new network science initiatives. To facilitate this effort, here we present the SOils DAta Harmonization database (SoDaH; https://lter.github.io/som-website, last access: 22 December 2020), a flexible database designed to harmonize diverse SOM datasets from multiple research networks. SoDaH is built on several network science efforts in the United States, but the tools built for SoDaH aim to provide an open-access resource to facilitate synthesis of soil carbon data. Moreover, SoDaH allows for individual locations to contribute results from experimental manipulations, repeated measurements from long-term studies, and local- to regional-scale gradients across ecosystems or landscapes. Finally, we also provide data visualization and analysis tools that can be used to query and analyze the aggregated database. The SoDaH v1.0 dataset is archived and available at https://doi.org/10.6073/pasta/9733f6b6d2ffd12bf126dc36a763e0b4 (Wieder et al., 2020).

54 ENVIRONMENTAL SCIENCES↗

HUMAN-ARTIFICIAL INTELLIGENCE TEAMING FOR THE U.S. NAVY: DEVELOPING A HOLISTIC RESEARCH ROADMAP

With the ever-increasing deluge of data and demand for warfighters to make decisions upon its analysis, U.S. defense strategy has prioritized the development of artificial intelligence (AI)/machine learning (ML) systems that can analyze multi-source data streams and suggest courses of action. However, there is a history of systems that have failed to be adopted by the warfighter due not only to unsolved technical challenges, but also a lack of usability or contributions to mission effectiveness, perceived or otherwise. To avoid this, research into human-AI teaming shows promise for developing AI systems that work with frontline operators. A recent National Academies of Sciences, Engineering, and Medicine report (NASEM, 2022) presented 57 research objectives in this area; however, the U.S. Navy requires a more-focused set of priorities, as it is impossible to tackle every priority. A workshop involving 23 human factors scientists, computer scientists, and active-duty sailors was organized at the Naval Information Warfare Center Pacific, resulting in a set of five research priorities spanning near-, mid-, and far-term time frames. This panel will summarize the results of this workshop, with a focus on the big questions both going into this workshop and coming out of it. The panel participants come from government, academia, and industry, providing perspective from the different kinds of organizations required to accomplish these research goals.

Wong, Jason↗

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗

Machine learning from RANS and LES to inform coarse grid simulations

Nuclear system thermal hydraulic analysis has historically relied on computationally inexpensive 1D codes. However, such tools are unable to capture multiscale multidimensional effects in large nuclear reactor enclosures. On the other hand, simulations with higher fidelity can be too expensive for such purposes. One of the ways to reduce computational cost is to perform simulations on a coarse grid, which, unfortunately, introduces large discretization errors. In this paper, two high-to-low data-driven approaches are investigated: (1) a coarse grid turbulence model to predict eddy viscosity and (2) correction of errors in coarse grid velocity fields. The approaches aim to reduce grid- and turbulence model-induced errors in coarse grid Reynolds-averaged Navier–Stokes (RANS) simulations. Two sources of high-fidelity data, RANS and large eddy simulations (LES), are explored. To extract the eddy viscosity from the LES data, an inverse optimization problem is solved. However, the LES eddy viscosity is shown to be comparable to the RANS eddy viscosity in terms of error reduction. Therefore, the directly available RANS eddy viscosity was used to develop a coarse grid data-driven turbulence model. Additionally, error correction in velocity is used to reduce the remaining uncertainties and bring the results closer to reality. In conclusion, the performance of the frameworks is demonstrated for a scaled upper plenum of a gas-cooled reactor facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Bandwidth extension and conversion efficiency improvements beyond phase matching limitations using cavity-enhanced OPCPA

The conversion efficiency and phase matching bandwidth of ultrafast optical parametric amplification (OPA) are constrained by the dispersion and nonlinear coefficient of the employed crystal as well as pulse shaping effects. In our work we show that an enhancement cavity resonant with the pump seeded at the full repetition rate of the pump laser can automatically reshape the small-signal gain in optical parametric chirped-pulse amplification (OPCPA) to achieve close-to-optimal operation. This new method termed cavity-enhanced OPCPA or C-OPCPA significantly increases both the gain bandwidth and the conversion efficiency, in addition to boosting gain for high-repetition-rate amplification. The goal in C-OPCPA is to arrive at a condition of impedance matching at all temporal coordinates, such that, in the absence of linear losses, all the incident pump power is dissipated in the nonlinear loss element, i.e., converted to signal and idler. The use of a low finesse enhancement cavity resonant with a low average power (<1W) and a high repetition rate (78MHz) pump source is shown to achieve more than 50% conversion efficiency into signal and idler from the coupled pump in an optical parametric process, whereas an equivalent amount of pump power in a single-pass configuration leads to negligible conversion. Additionally, the gain bandwidth is extended by a factor of 3-4 beyond the phase-matching limit. Our empirical observations are corroborated by a numerical analysis of depletion optimizing the single-pass case, which assesses the underlying impedance matching that is responsible for the observed performance improvements.

47 OTHER INSTRUMENTATION↗

Fluidized Bed Gasification For Conversion of Biomass and Waste Materials to Renewable Hydrogen

This project was undertaken in order to study the potential for hydrogen production, at low cost, from mixtures of biomass and municipal solid waste (MSW). This approach allows for the production of hydrogen with a very low fossil carbon burden, while taking advantage of tipping fees (associated with MSW) to improved process economics. The team sourced three primary feedstocks (wood, MSW, and waste plastics) and characterized them comprehensively using established techniques with a long track record in the field of gasification. All three primary feedstocks were highly reactive and lost most of their mass during initial devolatilization. The production of tars, including heavy tars, was quite high, and was most problematic in the case of the MSW and Waste Plastics feedstocks. Little practical difference was identified between the MSW and Waste Plastics materials, and the addition of bed-forming materials (dolomite and brown alumina) to the feedstocks was found to reduce production of tars during devolatilization under thermogravimetric analysis and/or Fischer Assay conditions. A series of four tests in a lab-scale bubbling-fluidized-bed gasifier, at 50 psig of pressure and about 825 C, confirmed these findings. Pellet feedstocks, broken into fragments, were used for these tests, and pellets comprised of 50% MSW and 50% biomass were found to be the best option in terms of fossil carbon burden, economic potential, and gasification characteristics. Tests were then undertaken in a pilot-scale gasifier facility based on the GTI U-Gas fluidized-bed gasification technology. The feedstock handling train of the 20 TPD U-Gas pilot-scale gasifier located in Des Plaines, IL, was operated under simulated gasification conditions, and the 50/50 pellets were found to be very robust and unproblematic. An improved design for the forward end of the feedstock injection screw of the gasifier was developed and installed. The design approach was based on improved passive cooling of the front-most shroud at the end of the screw, since this approach was found in comprehensive modeling studies to be more than sufficient to accomplish the project objectives associated with this phase of the work, while also avoiding thermal gradients that could have caused heat-stress-induced damage to the refractory around the feedstock inlet if an active cooling approach had been applied. Careful technoeconomic analysis (TEA) of two possible 1000 TPD facilities was carried out. The TEA of conversion of MSW with corn stover in one case, and MSW with woody feedstock in the other case, showed that both had the potential to provide hydrogen at about $1/kg (minimum selling price, 2018 dollar basis). Of the two TEA cases, the one that was based on the conversion of wood in the southeastern USA was found to have slightly better economic potential. The other TEA case was based on a real location in Nebraska and called for corn stover feedstock conversion along with MSW. An underserved communities outreach program plan was developed in cooperation with personnel from the Nebraska Public Power District.

08 HYDROGEN↗

Machine Learning Downscaling of SoilMERGE in the United States Southern Great Plains

SoilMERGE (SMERGE) is a root-zone soil moisture (RZSM) product that covers the entire continental United States and spans 1978 to 2019. Machine learning techniques, Random Forest (RF), eXtreme Gradient Boosting (XGBoost), and Gradient Boost (GBoost) downscaled SMERGE to spatial resolutions straddling the field scale domain (100 to 3000 m). Study area was northern Oklahoma and southern Kansas. The coarse resolution of SMERGE (0.125 degree) limits this product’s utility. To validate downscaled results in situ data from four sources were used that included: United States Department of Energy Atmospheric Radiation Measurement (ARM) observatory, United States Climate Reference Network (USCRN), Soil Climate Analysis Network (SCAN), and Soil moisture Sensing Controller and oPtimal Estimator (SoilSCAPE). In addition, RZSM retrievals from NASA’s Airborne Microwave Observatory of Subcanopy and Surface (AirMOSS) campaign provided a nearly spatially continuous comparison. Three periods were examined: era 1 (2016 to 2019), era 2 (2012 to 2015), and era 3 (2003 to 2007). During eras 1 and 2, RF outperformed XGBoost and GBoost, whereas during era 3 no model dominated. Performance was better during eras 1 and 2 as opposed to the pre-L band era 3. Improvements across all eras, regions, and models realized from downscaling included an increase in correlation from 0.03 to 0.42 and a decrease in ub RMSE from -0.0005 to -0.0118 m 3 /m 3 . This study demonstrates the feasibility of SMERGE downscaling opening the prospect for the development of a long-term RZSM dataset at a more desirable field-scale resolution with the potential to support diverse hydrometeorological and agricultural applications.

54 ENVIRONMENTAL SCIENCES↗

The Analysis of Using Satellite Soil Moisture Observations for Flood Detection, Evaluating Over the Thailand’s Great Flood of 2011

A flood monitoring and warning system provides critical information that can protect property and save lives. A basin-scale flood monitoring system requires an effective observation platform that offers extensive ground coverage of flood conditions, low latency, and high spatiotemporal resolution. While satellite imagery offers substantial spatial flood extent in detail due to its high spatial resolution, the coarse temporal resolution and cloud obstruction limit its near real-time application. Daily soil moisture data derived from satellite sensors at a scale of a few km can be used to monitor extreme wet surface conditions arising in flood occurrences. This study analyses the flood detection capabilities of several sources of soil moisture information, including the Soil Moisture and Ocean Salinity mission (SMOS), the Advanced Microwave Scanning Radiometer on EOS (AMSR-E), the Advanced SCATterometer (ASCAT) on MetOp, the Global Land Data Assimilation System (GLDAS), and the WaterGAP Global Hydrology Model (WGHM). In addition to soil moisture, the analysis includes measurements of surface reflectance from the Moderate Resolution Imaging Spectroradiometer (MODIS), precipitation measurements from the Tropical Rainfall Measuring Mission (TRMM), and terrestrial water storage estimates from the Gravity Recovery And Climate Experiment (GRACE) as proxies for flood inundations. The analysis was conducted over the Chao Phraya River Basin (CPB) in Thailand, where the Great Flood of 2011 led to one of the most significant economic losses in the country's history. Satellite-derived soil moisture exhibits a stronger correlation with the flood inundations than the precipitation, model-derived soil moisture, and terrestrial water storage data. SMOS soil moisture observation agrees best with the MODIS-derived flood extent/occurrence, both in terms of spatial distribution and timing, and providing approximated flood lead-time of one week or longer. A neural network constructed from SMOS and MODIS data is used to predict flood intensity/occurrence (given soil moisture input) with a predicted time window from eight days to thirty-two days. The short-term prediction (e.g., eight days) achieves the highest accuracy with an averaged recovery rate of approximately 60% (correlation coefficient). This study's results suggest a potential application of satellite soil moisture data in assisting flood monitoring and warning systems.

Satellite soil moisture↗

Thermodynamics of the internal combustion engine

The ubiquitous automobile has played both the role of hero and villain during the course of its history. At the present time it is being cast as the villain with increasing frequency because it has been indicted as a major source of air pollution. It is just this characterization that is responsible for our renewed interest in the automobile. When we speak of the automobile as a source of air pollution, we really mean the heart of the automobile, the internal combustion engine. The prototype of the modern internal combustion engine was first successfully operated by Nicholas August Otto in 1876. In spite of its long history, the internal combustion engine is still poorly understood in terms of the details of the physical processes that take place during its operations. The internal combustion engine operates on what is known as the Otto cycle, and Otto cycle calculations have been with us for many years. The emphasis in earlier calculations was usually on the performance rather than on the pollution aspects. Performance characteristics can often be estimated reasonably accurately even with a considerable degree of simplification. We are interested in the pollution aspects of the internal combustion engine, and therefore the calculations must be performed with more care. Our extensive experience with equilibrium calculations places us in a position to carry out such careful calculations. We have just recently begun such an analysis, and its objective is the prediction of exhaust gas properties. I will describe our progress by first reviewing the Otto cycle, and then I will illustrate the effect of various assumptions on the calculated results.

Frank J. Zeleznik↗

Broad-Band Measurements of Cen X-3 With XTE and CGRO

Centaurus X-3 has played a key role in the development of our understanding of galactic x-ray binary sources. Timing analysis of the UHURU x-ray observations for the luminous Cen X-3 source (L approximately 10(exp 38) erg/s) revealed the first evidence for coherent x-ray pulsations from an object in a binary system (Giaconni 1971; Schreier 1972). It was quickly understood that the luminous pulsed x-ray emission could be generated by the accretion of matter from a companion star onto a rotating neutron star and led to the adoption of binary star models as the fundamental model for galactic x-ray sources (e.g. Pringle and Rees 1972; Lamb 1973). Based on modeling and refined observations since the original measurements, we now believe that Cen X-3 is a high mass x-ray binary system that contains a disk-fed pulsar with a period of 4.84 seconds that is in a 2.087 day orbit around an O-star companion. Since the pulsar discovery, its period has been intermittently monitored and those studies show a long term spin-up of the pulsar punctuated by short intervals of spin-down (e.g. Finger 1994). The implied torques are thought to originate from the interaction of an accretion disk with the magnetic field of a neutron star (Ghosh and Lamb 1979).

Vestrand, W. Thomas↗

Stochastic Analysis of Orbital Lifetimes of Spacecraft

A document discusses (1) a Monte-Carlo-based methodology for probabilistic prediction and analysis of orbital lifetimes of spacecraft and (2) Orbital Lifetime Monte Carlo (OLMC)--a Fortran computer program, consisting of a previously developed long-term orbit-propagator integrated with a Monte Carlo engine. OLMC enables modeling of variances of key physical parameters that affect orbital lifetimes through the use of probability distributions. These parameters include altitude, speed, and flight-path angle at insertion into orbit; solar flux; and launch delays. The products of OLMC are predicted lifetimes (durations above specified minimum altitudes) for the number of user-specified cases. Histograms generated from such predictions can be used to determine the probabilities that spacecraft will satisfy lifetime requirements. The document discusses uncertainties that affect modeling of orbital lifetimes. Issues of repeatability, smoothness of distributions, and code run time are considered for the purpose of establishing values of code-specific parameters and number of Monte Carlo runs. Results from test cases are interpreted as demonstrating that solar-flux predictions are primary sources of variations in predicted lifetimes. Therefore, it is concluded, multiple sets of predictions should be utilized to fully characterize the lifetime range of a spacecraft.

Sasamoto, Washito↗