Search NASA⌕ Search

SEARCH · Search NASA

Results for “conditional probability distribution”

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.

323 records · Page 18

Ventifact Formation in the Mojave Desert: Field Analogs for Martian Processes

Two field studies in the Mojave Desert, California, shed light on processes of ventifact formation. The field sites are located on a ridge at Little Cowhole Mountain, which lies approximately 12 km south of Baker, and on an unnamed ridge situated along the northern boundary of the Mojave River Sink (Rasor Road site). The rocks at Little Cowhole Mountain are a blue-grey marble/dolomite, whereas those at Rasor Road are Miocene volcanic rocks (basalt). At both sites the abrasive agent is a fine-grained aeolian sand which was probably derived largely from the Mojave River. There are minimal modem inputs of sand to either site: abrasion occurs as a result of unique climatic and topographic conditions which allow pre-existing sand to be recycled from one aspect of the ridge to the other. Climatic conditions are well suited for ventifact formation. Owing, to the dry climate (marked by low average relative humidity, infrequent dew, and low annual rainfall), rates of chemical weathering are low. Where resurfacing of the rocks by sand abrasion proceeds at a rate greater than weathering, the ventifacts are considered "active." Active ventifacts are found atop and straddling the ridge crests, in the zone of maximum wind velocity and sediment supply. Inactive ventifacts occur where modem weathering Processes exceed abrasion rates; principally on the basal two-thirds of the hillslope, where wind velocity and sediment supply are lower. At intermediate locations between the slope base and crest, ventifacts are either active or inactive, depending on local conditions. The presence of relict ventifacts at the study sites, as well as elsewhere in the eastern Mojave Desert, suggests that the conditions for venti fact formation must have been more intense and extensive in the past. Together, the studies illustrate that the processes that interact to form ventifacts are highly complex, and must be studied at many scales. Small-scale effects, such as local topography, plant cover, or even the spatial distribution of boulders, strongly influence the formation of each individual ventifact. Mesoscale effects (hillslope form, total availability of sediment, seasonality of winds, etc.) determine the distribution of overall erosional energy and the location of active and relict forms.

Laity, J. E.↗

The galactic globular cluster system

We explore correlations between various properties of Galactic globular clusters, using a database on 143 objects. Our goal is identify correlations and trends which can be used to test and constrain theoretical models of cluster formation and evolution. We use a set of 13 cluster parameters, 9 of which are independently measured. Several arguments suggest that the number of clusters still missing in the obscured regions of the Galaxy is of the order of 10, and thus the selection effects are probably not severe for our sample. Known clusters follow a power-law density distribution with a slope approximately -3.5 to -4, and an apparent core with a core radius approximately 1 kpc. Clusters show a large dynamical range in many of their properties, more so for the core parameters (which are presumably more affected by dynamical evolution) than for the half-light parameters. There are no good correlations with luminosity, although more luminous clusters tend to be more concentrated. When data are binned in luminosity, several trends emerge: more luminous clusters tend to have smaller and denser cores. We interpret this as a differential survival effect, with more massive clusters surviving longer and reaching more evolved dynamical states. Cluster core parameters and concentrations also correlate with the position in the Galaxy, with clusters closer to the Galactic center or plane being more concentrated and having smaller and denser cores. These trends are more pronounced for the fainter (less massive) clusters. This is in agreement with a picture where tidal shocks form disk or bulge passages accelerate dynamical evolution of clusters. Cluster metallicities do not correlate with any other parameter, including luminosity and velocity dispersion; the only detectable trend is with the position in the Galaxy, probably reflecting Zinn's disk-halo dichotomy. This suggests that globular clusters were not self-enriched systems. Velocity dispersions show excellent correlations with luminosity and surface brightness. Their origin is not well understood, but they may well reflect initial conditions of cluster formation, and perhaps even be used to probe the initial density perturbation spectrum on a approximately 10(exp 6) solar mass scale. Core radii and concentrations play a role of a 'second parameter' in these correlations. While a global manifold of cluster properties has a high statistical dimensionality (D greater than 4), a subset of structural, photometric, and dynamical parameters forms a statistically three-dimensional family, as expected from objects following King models; we propose to call this set of quantities the King Manifold. Some of the observed correlations may be usable as distance indicator relations for globular clusters.

Djorgovski, S.↗

Utilitarian models of the solar nebula

Models of the primitive solar nebula based on a combination of theory, observations of T Tauri stars, and global conservation laws are presented. The models describe the motions of nebular gas, mixing of interstellar material during the formation of the nebula, and evolution of thermal structure in terms of several characteristic parameters. The parameters describe key aspects of the protosolar cloud (its rotation rate and collapse rate) and the nebula (its mass relative to the Sun, decay time, and density distribution). For most applications, the models are heuristic rather than predicted. Their purpose is to provide a realistic context for the interpretation of solar system data, and to distinquish those nebular characteristics that can be specified with confidence, independently of the assumtions of particular models, form those that are poorly constrained. It is demonstrated that nebular gas typically experienced large radial excursions during the evolution of the nebula and that both inward and outward mean radial velocities on the order of meters per second occured in the terrestrial planet region, with inward velocities predominant for most ofthe evolution. However, the time history of disk size, surface density, and radial velocities are sensitive to the total angular momentun of the protosolar cloud, which cannot be constrained by purely theoretical considerations.It is shown that a certain amount of 'formational' mixing of interstellar material was an inevitable consequenc of nebular mass and angular momentum transport during protostellar collapse, regardless of the specific transport mechanisms invloved. Even if the protosolar cloud was initially homogeneous, this mixing was important because it had the effect of mingling presolar material that had experienced different degrees of thermal processing during collapse and passage through the accertion shock. Nebular thermal structure is less sensitive to poorly constrained parameters than is dynamical history. A simple criterion is derived for the condition that silicate grains are evaporated at midplane, and it is argued that this condition was probably fulfilled early in nebular history. Cooling of a hot nebula due tocoagulation of dust and consequent local reduction of optical depth is examined, and it is shown how such a process leads naturally to an enrichment of rock-forming elements in the gas phase.

Cassen, Patrick↗

Probability of Decompression Sickness and Venous Gas Emboli from 49 NASA Hypobaric Chamber Tests with Reference to Exploration Atmosphere

Introduction: Decompression sickness (DCS) is a complex biophysical event; it combines human perception of pain, for instance, and the presence of a gas phase in the tissues. Living tissues are complex and dynamic. Micronuclei and later bubbles may or may not form given what appears to be the same conditions. Even when bubbles grow, symptoms may or may not develop under what appears to be the same conditions. Therefore, at this time it is appropriate to consider DCS as a probabilistic rather than a deterministic event. Methods: Probabilistic models about hypobaric DCS and venous gas emboli (VGE) require a large amount of quality research data, a definition of decompression dose using physical and physiologic variables, and a flexible analytical approach that can quantify the association between each outcome and all covariates of interest (assuming independence between DCS and VGE) and then ultimately be extended to acknowledge dependencies between DCS and VGE. Our DCS and VGE data are from 1,031 hypobaric decompressions from 1983 to 2016. A total of 577 humans participated in 49 hypobaric chamber tests to evaluate denitrogenation procedures used by astronauts in the Space Shuttle and International Space Station programs. We defined decompression dose as the ratio of computed nitrogen tension in a theoretical 360-minute half-time compartment to ambient pressure, which accounts for denitrogenation and exposure pressure as well as explanatory variables such as age, sex, body mass index, and the presence or absence of ambulation as part of exercise at the exposure pressure. A parametric survival model, using a log-logistic distribution, was used to quantify the time to development of DCS, VGE, and Grade IV VGE. Results: Our survival estimates are applicable to simple hypobaric decompressions, such as depressurizations in 5 to 30 minutes to exposure pressures between 4 and 10 pounds per square inch absolute (psia) and after minutes to hours of denitrogenation, either under resting or exercise conditions to accelerated denitrogenation. The regressions are applicable to exposures between 2 to 6 hours and under conditions of ambulation or no ambulation as part of exercise at the test pressure. We estimate that an exposure to 4.3 psia with simulated extravehicular activity (EVA) that includes ambulation after equilibration to the exploration atmosphere at 8.2 psia with a 34% oxygen atmosphere will result in 3.1% DCS (1.8% to 5.2 95% confidence interval), 23.2% VGE (16.7 to 31.2%), and 8.5% Grade IV VGE (4.7 to 14.7%) in equal samples of men and women exposed for 6 hours. Discussion: Probabilistic models for DCS, VGE, and Grade IV VGE can be used to inform those that plan future EVAs. Their applications are useful to quantify the risk of DCS and VGE in astronauts that perform EVAs in low-pressure space suits while in space or while exploring the surfaces of the moon or Mars.

Johnny Conkin↗

Transient Reliability Analysis Capability Developed for CARES/Life

The CARES/Life software developed at the NASA Glenn Research Center provides a general-purpose design tool that predicts the probability of the failure of a ceramic component as a function of its time in service. This award-winning software has been widely used by U.S. industry to establish the reliability and life of a brittle material (e.g., ceramic, intermetallic, and graphite) structures in a wide variety of 21st century applications.Present capabilities of the NASA CARES/Life code include probabilistic life prediction of ceramic components subjected to fast fracture, slow crack growth (stress corrosion), and cyclic fatigue failure modes. Currently, this code can compute the time-dependent reliability of ceramic structures subjected to simple time-dependent loading. For example, in slow crack growth failure conditions CARES/Life can handle sustained and linearly increasing time-dependent loads, whereas in cyclic fatigue applications various types of repetitive constant-amplitude loads can be accounted for. However, in real applications applied loads are rarely that simple but vary with time in more complex ways such as engine startup, shutdown, and dynamic and vibrational loads. In addition, when a given component is subjected to transient environmental and or thermal conditions, the material properties also vary with time. A methodology has now been developed to allow the CARES/Life computer code to perform reliability analysis of ceramic components undergoing transient thermal and mechanical loading. This means that CARES/Life will be able to analyze finite element models of ceramic components that simulate dynamic engine operating conditions. The methodology developed is generalized to account for material property variation (on strength distribution and fatigue) as a function of temperature. This allows CARES/Life to analyze components undergoing rapid temperature change in other words, components undergoing thermal shock. In addition, the capability has been developed to perform reliability analysis for components that undergo proof testing involving transient loads. This methodology was developed for environmentally assisted crack growth (crack growth as a function of time and loading), but it will be extended to account for cyclic fatigue (crack growth as a function of load cycles) as well.

Nemeth, Noel N.↗

Critical review of electrical conductivity measurements and charge distribution analysis of magnesium oxide

The electrical conductivity sigma of MgO single crystals shows a sharp increase at 500-800 C, in particular of sigma surface, generally attributed to surface contamination. Charge Distribution Analysis (CDA), a new technique providing information on fundamental properties that was previously unavailable, allows for the determination of surface charges, their sign and associated internal electric field. Data on 99.99% purity, arc-fusion grown MgO crystals show that mobile charge carriers start to appear in the bulk of the MgO crystals between 200 and 400 C when sigma (measured by conventional techniques) is in t he 10(exp -14) to 10(exp -16) /omega/cm range. Above 500 C, as sigma increases to 10(exp -6) to 10(exp -7)/omega/cm, more charges appear giving rise to a strong positive surface charge supported by a strong internal field. This indicates that charges are generated in the bulk and diffuse to the surface by an internally controlled process. On the basis of their positive sign they are identified as holes (defect electrons). Because of the low cation content of these very pure MgO crystals, theses holes cannnot be associated with transition metal impurties. Instead, they are associated with the O(2-) sublattice, e.g. consist of O(-) states or positive holes. This conclusion is supported by magnetic susceptibility data showing the appearance of 1000 +/- 500 ppm paramagnetic species between 200-500 C. The magnetic data are consistent with strongly coupled, diamagnetic O(-) pairs below 200-500 C, chemically equivalent to peroxy anions, O2(2-), and probably associated with cation vacancies in the MgO matrix. The formation of O2(2-) in arc-fusion grown MgO crystals is very unexpected because of the highly reducing growth conditions. Their presence implies an internal redox reaction involving dissolved 'water' by which OH(-) pairs convert to O2(2-) plus H2 molecules. This redox conversion is supported by mass spectroscopic measurements of the H2 release from highly OH(-)-doped, finely divided MgO and by wet-chemical analysis of its oxidant concentration.

Freund, Friedemann↗

An Overview of the NCC Spray/Monte-Carlo-PDF Computations

This paper advances the state-of-the-art in spray computations with some of our recent contributions involving scalar Monte Carlo PDF (Probability Density Function), unstructured grids and parallel computing. It provides a complete overview of the scalar Monte Carlo PDF and Lagrangian spray computer codes developed for application with unstructured grids and parallel computing. Detailed comparisons for the case of a reacting non-swirling spray clearly highlight the important role that chemistry/turbulence interactions play in the modeling of reacting sprays. The results from the PDF and non-PDF methods were found to be markedly different and the PDF solution is closer to the reported experimental data. The PDF computations predict that some of the combustion occurs in a predominantly premixed-flame environment and the rest in a predominantly diffusion-flame environment. However, the non-PDF solution predicts wrongly for the combustion to occur in a vaporization-controlled regime. Near the premixed flame, the Monte Carlo particle temperature distribution shows two distinct peaks: one centered around the flame temperature and the other around the surrounding-gas temperature. Near the diffusion flame, the Monte Carlo particle temperature distribution shows a single peak. In both cases, the computed PDF's shape and strength are found to vary substantially depending upon the proximity to the flame surface. The results bring to the fore some of the deficiencies associated with the use of assumed-shape PDF methods in spray computations. Finally, we end the paper by demonstrating the computational viability of the present solution procedure for its use in 3D combustor calculations by summarizing the results of a 3D test case with periodic boundary conditions. For the 3D case, the parallel performance of all the three solvers (CFD, PDF, and spray) has been found to be good when the computations were performed on a 24-processor SGI Origin work-station.

Raju, M. S.↗

HST eclipse mapping of dwarf nova OY Carinae in quiescence: An 'Fe II curtain' with Mach approx. = 6 velocity dispersion veils the white dwarf

Hubble Space Telescope (HST) observations of the eclipsing dwarf nova OY Car in its quiescent state are used to isolate the ultraviolet spectrum (1150-2500 A at 9.2 A Full Width at Half Maximum (FWHM) resolution) of the white dwarf, the accretion disk, and the bright spot. The white dwarf spectrum has a Stark-broadened photospheric L(alpha) absorption, but is veiled by a forest of blended Fe II features that we attribute to absorption by intervening disk material. A fit gives T(sub w) approx. = 16.5 x 10(exp 3) K for the white dwarf with a solar-abundance, log g = 8 model atmosphere, and T approx. = 10(exp 4) K, n(sub e) approx. = 10(exp 13)/cu cm, N(sub H) approx. = 10(exp 22) sq cm, and velocity dispersion delta V approx. = 60 km/s for the veil of homogeneous solar-abundance local thermodynamic equilibrium (LTE) gas. The veil parameters probably measure characteristic physical conditions in the quiescent accretion disk or its chromosphere. The large velocity dispersion is essential for a good fit; it lowers (chi square)/778 from 22 to 4. Keplerian shear can produce the velocity dispersion if the veiling gas is located at R approx. = 5 R(sub W) with (delta R)/R approx. = 0.3, but this model leaves an unobscured view to the upper hemisphere of the white dwarf, incompatible with absorptions that are up to 80% deep. The veiling gas may be in the upper atmosphere of the disk near its outer rim, but we then require supersonic (Mach approx. = 6) but sub-Keplerian (delta V/V(sub Kep) approx. = 0.07) velocity disturbances in this region to produce both the observed radial velocity dispersion and vertical motions sufficient to elevate the gas to z/R = cos i = 0.12. Such motions might be driven by the gas stream, since it may take several Kepler periods to reestablish the disk's vertical hydrostatic equilibrium. The temperature and column density of the gas we see as Fe II absorption in the ultraviolet are similar to what is required to produce the strong Balmer jump and line emissions seen in optical spectra of OY Car and similar quiescent dwarf novae. The outer accretion disk is detected at mid-eclipse with a spectrum that rises from 0.05 to 0.3 mJy between 2000 and 2500 A, consistent with combinations of cool blackbodies, blended Fe II emission lines, and Balmer continuum emission. The total disk flux density is 0.5 mJy at 2500 A, and this shallow disk eclipse implies a roughly flat surface brightness distribution. The bright spot, somewhat bluer than the disk, has a flux density rising from 0.05 to 0.15 mJy between 1600 and 2500 A. The C IV emission line has a broad shallow eclipse, but the radial velocity variations observed during the eclipse do not clearly distinguish between a disk or wind origin. The only possible indications of boundary layer emission are fast UV flares that appear to arise from near the central object -- not from the bright spot.

Horne, Keith↗

SAO and Kelvin Waves in the EuroGRIPS GCMS and the UK Meteorological Offices Analyses

This work is an intercomparison of four tropospheric-stratospheric climate models, the Unified Model (UM) of the U.K. Meteorological Office (UKMO), the model of the Free University in Berlin (FUB). the ARPEGE-climat model of the National Center for Meteorological Research (CNRM), and the Extended UGAMP GCM (EUGCM) of the Center for Global Atmospheric Modelling (CGAM), against the UKMO analyses. This comparison has been made in the framework of the "GSM-Reality Intercomparison Project for SPARC" (GRIPS). SPARC (Stratospheric Processes and their Role in Climate) aims are to investigate the effects of the middle atmosphere on climate and the GRIPS purpose is to organized a comprehensive assessment of current Middle Atmosphere-Climate Models (MACMs). The models integrations were made without identical contraints e.g. boundary conditions, incoming solar radiation). All models are able to represent the dominant features of the extratropical circulation. In this paper, the structure of the tropical winds and the strengths of the Kelvin waves are examined. Explanations for the differences exhibited. between the models. as well as between models and analyses, are also proposed. In the analyses a rich spectrum of waves (eastward and westward) is present and contributes to drive the SAO (SemiAnnual Oscillation) and the QBO (Quasi-Biennal Oscillation). The amplitude of the Kelvin waves is close to the one observed in UARS (Upper Atmosphere Research Satellite) data. In agreement with observations, the Kelvin waves generated in the models propagate into the middle atmosphere as wave packets which underlines convective forcing origin. In most models, slow Kelvin waves propagate too high and are hence overestimated in the upper stratosphere and in the mesosphere, except for the UM which is more diffusive. These waves are not sufficient to force realistic westerlies of the QBO or SAO westerly phases. If the SAO is represented by all models only two of them are able to generate westerlies between 10 hPa and 50 hPa. The importance of the role played by subgrided gravity waves is more and more recognized. Actually, the EUGCM which includes a parametrization of gravity waves with a non-zero phase speed is able to simulate. with however some unrealistic features, clear easterly to westerly transitions as well as westerlies downward propagations. Thermal damping is also important in the westerlies forcing in the stratosphere. The model ARPEGE-climat shows more westerlies in the stratosphere than tile other three models probably due to the use of a simplified scheme to predict the ozone distribution in the middle atmosphere.

Amodei, M.↗

A Report on Stochastic Fairness Queueing (SFQ) Experiments

SRI International (SRI) has developed an improved queueing algorithm, known as Stochastic Fairness Queueing (SFQ), for best-effort traffic (i.e., traffic that does not require any guaranteed service). SFQ is a probablistic variant of strict fair queueing where instead of a single queue being allocated per flow, a fixed number of queues are used and a hash function maps the IP source and destination to a particular queue. A seed to the hash function is also perturbed occasionally to help distribute the flows amongst different queues when more than one flow maps to the same queue during the lifetime of the flow. SFQ provides 'fair' access by trying to ensure that each flow from source to destination host obtains equal access to the available bandwidth. This report covers a series of experiments performed on DARTnet evaluating the behavior and performance of SFQ against a FIFO queueing discipline. These experiments were designed to show SFQ's advantages and performance, and include tests demonstrating: Fair utilization of available resources; Starvation prevention; Graceful degradation under overload conditions; and Resource usage. In general, the experiments do show that SFQ is better than FIFO queueing at allocating bandwidth equally among a set of flows. SFQ also prevents a stream from dominating the available bandwidth, which seems to be a tendency with FIFO queueing (i.e., if a flow demands more than its share of the available bandwidth, with FIFO queueing that stream receives a disproportionate amount when compared to flows demanding less than their share). Furthermore, SFQ seems to reward 'nice' users of the network by providing a lower variance in delay and more throughput when their resource demand is less than their available share. Both SFQ and FIFO queueing seem to degrade fairly well as the network becomes saturated and to recover well as the network becomes less congested. Not unexpectedly, FIFO queueing is a little more efficient than SFQ-the delays are less and the throughput slightly higher because SFQ requires more processing. However, the performance difference between the two queueing disciplines is relatively small. However, the experiments do point out some interesting behavior. FIFO queueing can behave better than SFQ with seed perturbation. We recommend further evaluation of the hash function and the seed perturbation technique. There are probably weaknesses in their current selection that cause this unexpected behavior. SFQ also seems to possess good scaling properties. To verify this, more experiments with a larger number of streams from more hosts need to be executed and examined, including the staggered introduction of streams. Staggering the streams may prove important, because graphs in the degradation experiment revealed some unexpected increases and decreases in throughput, which should be examined. This may again be due to the interaction of the hash function with the seed perturbation but it may also be related to some other unknown problem.

Denny, Barbara A.↗

Convective instabilities of electromagnetic ion cyclotron waves in the outer magnetosphere

The path-integrated linear growth of electromagnetic ion cyclotron waves in the outer (L is greater than or equal to 7) magnetosphere is investigated using a realistic thermal plasma distribution with an additional anisotropic energetic ring current H(+) to provide free energy for instability. The results provide a realistic simulation of the recent Active Magneto- spheric Particle Tracer Explorers (AMPTE) observations. For conditions typical of the dayside magnetosphere, high plasma beta effects reduce the group velocity and significantly increase the spatial growth rates for left-handed polarized instabilities just below the helium gyrofrequency Omega(sub He(+)), and on the guided mode above Omega(sub He(+)) but below the cross over frequency omega(sub cr). Relatively high densities, typical of the afternoon local time sector, favor these low group velocity effects for predominantly field-aligned waves. Lower densities, typical of those found in the early morning local time sector, increase the group velocity but allow strong convective instabilities at high normalized frequencies well above Omega(sub He(+)). These waves are reflected in the magnetosphere and can exist for several equatorial transits without significant damping. They are left-handed polarized only on the first equatorial crossing and become linearly polarized for the remainder of the ray path. Consequently, these waves should be observed with basically linear polarization at all frequencies and all latitudes in the early morning local time sector. Wave growth below Omega(sub He(+)) is severely limited owing to the narrow bandwidth for instability and the small resonant path lengths. In the afternoon sector, where plasma densities can exceed 10(exp 7)/cu m, intense convective amplification is possible both above and below Omega(sub He(+)). Waves below Omega(sub He(+)) are not subject to reflection when the O(+) concentration is small and therefore should be observed with left-handed polarization near the equator and essentially linear polarization at higher latitudes. Since the He(+) concentration is usually large in the afternoon sector, guided mode waves above Omega(sub He(+)) reflect to form a background distribution with basically linear polarization. We suggest that the strong left-handed polarized emissions observed by AMPTE in the afternoon sector near the equator are probably due to strongly growing low group velocity waves at frequencies just below Omega(sub He(+)), and on the guided mode above Omega(sub He(+)).

Horne, Richard B.↗

Interactive Soil Dust Aerosol Model in the GISS GCM: Sensitivity of the Soil Dust Cycle to Radiative Properties of Soil Dust Aerosols - Part 1

The sensitivity of the soil dust aerosol cycle to the radiative forcing by soil dust aerosols is studied. Four experiments with the NASA/GISS atmospheric general circulation model, which includes a soil dust aerosol model, are compared, all using a prescribed climatological sea surface temperature as lower boundary condition. In one experiment, dust is included as dynamic tracer only (without interacting with radiation), whereas dust interacts with radiation in the other simulations. Although the single scattering albedo of dust particles is prescribed to be globally uniform in the experiments with radiatively active dust, a different single scattering albedo is used in those experiments to estimate whether regional variations in dust optical properties, corresponding to variations in mineralogical composition among different source regions, are important for the soil dust cycle and the climate state. On a global scale, the radiative forcing by dust generally causes a reduction in the atmospheric dust load corresponding to a decreased dust source flux. That is, there is a negative feedback in the climate system due to the radiative effect of dust. The dust source flux and its changes were analyzed in more detail for the main dust source regions. This analysis shows that the reduction varies both with the season and with the single scattering albedo of the dust particles. By examining the correlation with the surface wind, it was found that the dust emission from the Saharan/Sahelian source region and from the Arabian peninsula, along with the sensitivity of the emission to the single scattering albedo of dust particles, are related to large scale circulation patterns, in particular to the trade winds during Northern Hemisphere winter and to the Indian monsoon circulation during summer. In the other regions, such relations to the large scale circulation were not found. There, the dependence of dust deflation to radiative forcing by dust particles is probably dominated by physical processes with short time scales. The experiments show that dust radiative forcing can lead to significant changes both in the soil dust cycle and in the climate state. To estimate dust concentration and radiative forcing by dust more accurately, dust size distributions and dust single scattering albedo in the model should be a function of the source region, because dust concentration and climate response to dust radiative forcing are sensitive to dust radiative parameters.

Perlwitz, Jan↗

CFD Study of Turbo-Ramjet Interactions in Hypersonic Airbreathing Propulsion System

Advanced airbreathing propulsion systems used in Mach 4-6 mission scenarios, usually involve turbo-ramjet configurations. As the engines transition from turbojet to ramjet, there is an operational envelope where both engines operate simultaneously. In the first phase of our study, an over/under nozzle configuration was analyzed. The two plumes from the turbojet and ramjet interact at the end of a common 2-D cowl, where they both reach an approximate Mach 3.0 condition and then jointly expand to Mach 3.6 at the common nozzle exit plane. For the problem analyzed, the turbojet engine operates at a higher nozzle pressure ratio than the ramjet, causes the turbojet plume overpowers the ramjet plume, deflecting it approximately 12 degrees downward and in turn the turbojet plume is deflected 6 degrees upward. In the process, shocks were formed at the deflections and a shear layer formed at the confluence of the two jets. This particular case was experimentally tested and the data were used to compare with a computational fluid dynamics (CFD) study using the PARC2D code. The CFD results were in good agreement with both static pressure distributions on the cowl separator and on nozzle walls. The thrust coefficients were also in reasonable agreement. In addition, inviscid relationships were developed around the confluence point, where the two exhaust jets meet, and these results compared favorably with the CFD results. In the second phase of our study, a 3-D CFD solution was generated to compare with the 2-D solution. The major difference between the 2-D and 3-D solutions was the interaction of the shock waves, generated by the plume interactions, on the sidewall. When a shock wave interacts with a sidewall and sidewall boundary layer, it is called a glancing shock sidewall interaction. These interactions entrain boundary layer flow down the shockline into a vortical flow pattern. The 3-D plots show the streamlines being entrained down the shockline. The pressure of the flow also decreases slightly as the sidewall is approached. Other difference between the 2-D and 3-D solutions were a lowering of the nozzle thrust coefficient value from 0.9850 (2-D) to 0.9807 (3-D), where the experimental value was 0.9790. In the third phase of our study, a different turbo-ramjet configuration was analyzed. The confluence of a supersonic turbojet and a subsonic ramjet in the turbine based combined-cycle (TBCC) propulsion system was studied by a 2-D CFD code. In the analysis, Mach 1.4 primary turbojet was mixed with the subsonic ramjet secondary flow in an ejector mode operation. Reasonable agreements were obtained with the supplied I-D TBCC solutions. For low downstream backpressure, the Fabri choke condition (Break-Point condition) was observed in the secondary flow within mixing zone. For sufficient high downstream backpressure, the Fabri choke no longer exist, the ramjet flow was reduced and the ejector flow became backpressure dependent. Highly non-uniform flow at ejector exit were observed, indicated that for smooth downstream combustion, the mixing of the two streams probably required some physical devices.

Chang, Ing↗

Radiation-Driven Warping: Nonisothermal Disks - 2

Recent work by Pringle and by Maloney, Begelman, & Pringle has shown that geometrically thin, optically thick, accretion disks are unstable to warping driven by radiation torque from the central source. This work was confined to isothermal (i.e., surface density Sigma varies as R(sup -3/2) disks. In this paper we generalize the study of radiation-driven warping to include general power-law surface density distributions, Sigma varies as R(sup -delta).We consider the range from Delta = 3/2 (the isothermal case) to Delta = -3/2, which corresponds to a radiation-pressure-supported disk; this spans the range of surface density distributions likely to be found in real astrophysical disks. In all cases there are an infinite number of zero-crossing solutions (i.e., solutions that cross the equator), which are the physically relevant modes if the outer boundary of the disk is required to lie in a specified plane. However, unlike the isothermal disk, which is the degenerate case, the frequency eigenvalues for Delta does not equal 3/2 are all distinct. In all cases the location of the zero moves outward from the steady state (pure precession) value with increasing growth rate; thus, there is a critical minimum size for unstable disks. Modes with zeros at smaller radii are damped. The critical radius and the steady state precession rate depend only weakly on Delta. An additional analytic solution has been found for Delta = 1. The case Delta = 1 divides the solutions into two qualitatively different regimes. For Delta greater than or equal to 1, the fastest growing modes have maximum warp amplitude, close to the disk outer edge, and the ratio of Beta(sub max) to the warp amplitude at the disk inner edge, Beta(sub o), is much greater than 1. For Delta less than 1, Beta(sub max/Beta(sub o) approximately equals 1, and the warp maximum steadily approaches the origin as Delta decreases. This implies that nonlinear effects must be important if the warp extends to the disk inner edge for Delta greater than or equal to 1, but for Delta less than 1 nonlinearity will be important only if the warp amplitude is large at the origin. Because of this qualitative difference in the shapes of the warps, the effects of shadowing of the central source by the warp will also be very different in the two regimes of Delta. This has important implications for radiation-driven warping in X-ray binaries, for which the value of Delta characterizing the disk is likely to be less than unity. In real accretion disks the outer boundary condition is likely to be different from the zero-crossing condition that we have assumed. In accretion disks around massive black holes in active galactic nuclei, the disk will probably become optically thin before the outer disk boundary is reached, whereas in X-ray binaries there will be an outer disk region (outside the circularization radius) in which the inflow velocity is zero but angular momentum is still transported. We show that in both these cases the solutions are similar to the zero-crossing eigenfunctions.

Maloney, Philip R.↗

Robotic and Human-Tended Collaborative Drilling Automation for Subsurface Exploration

Future in-situ lunar/martian resource utilization and characterization, as well as the scientific search for life on Mars, will require access to the subsurface and hence drilling. Drilling on Earth is hard - an art form more than an engineering discipline. Human operators listen and feel drill string vibrations coming from kilometers underground. Abundant mass and energy make it possible for terrestrial drilling to employ brute-force approaches to failure recovery and system performance issues. Space drilling will require intelligent and autonomous systems for robotic exploration and to support human exploration. Eventual in-situ resource utilization will require deep drilling with probable human-tended operation of large-bore drills, but initial lunar subsurface exploration and near-term ISRU will be accomplished with lightweight, rover-deployable or standalone drills capable of penetrating a few tens of meters in depth. These lightweight exploration drills have a direct counterpart in terrestrial prospecting and ore-body location, and will be designed to operate either human-tended or automated. NASA and industry now are acquiring experience in developing and building low-mass automated planetary prototype drills to design and build a pre-flight lunar prototype targeted for 2011-12 flight opportunities. A successful system will include development of drilling hardware, and automated control software to operate it safely and effectively. This includes control of the drilling hardware, state estimation of both the hardware and the lithography being drilled and state of the hole, and potentially planning and scheduling software suitable for uncertain situations such as drilling. Given that Humans on the Moon or Mars are unlikely to be able to spend protracted EVA periods at a drill site, both human-tended and robotic access to planetary subsurfaces will require some degree of standalone, autonomous drilling capability. Human-robotic coordination will be important, either between a robotic drill and humans on Earth, or a human-tended drill and its visiting crew. The Mars Analog Rio Tinto Experiment (MARTE) is a current project that studies and simulates the remote science operations between an automated drill in Spain and a distant, distributed human science team. The Drilling Automation for Mars Exploration (DAME) project, by contrast: is developing and testing standalone automation at a lunar/martian impact crater analog site in Arctic Canada. The drill hardware in both projects is a hardened, evolved version of the Advanced Deep Drill (ADD) developed by Honeybee Robotics for the Mars Subsurface Program. The current ADD is capable of 20m, and the DAME project is developing diagnostic and executive software for hands-off surface operations of the evolved version of this drill. The current drill automation architecture being developed by NASA and tested in 2004-06 at analog sites in the Arctic and Spain will add downhole diagnosis of different strata, bit wear detection, and dynamic replanning capabilities when unexpected failures or drilling conditions are discovered in conjunction with simulated mission operations and remote science planning. The most important determinant of future 1unar and martian drilling automation and staffing requirements will be the actual performance of automated prototype drilling hardware systems in field trials in simulated mission operations. It is difficult to accurately predict the level of automation and human interaction that will be needed for a lunar-deployed drill without first having extensive experience with the robotic control of prototype drill systems under realistic analog field conditions. Drill-specific failure modes and software design flaws will become most apparent at this stage. DAME will develop and test drill automation software and hardware under stressful operating conditions during several planned field campaigns. Initial results from summer 2004 tests show seven identifi distinct failure modes of the drill: cuttings-removal issues with low-power drilling into permafrost, and successful steps at executive control and initial automation.

Glass, Brian↗

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Space Policy Directive-1 has led to NASA partnerships with commercial entities on procurement which includes the development of the Human Landing System (HLS) [1]. With the goal of delivering human crew to the lunar surface by 2024, system uncertainties become an important obstacle to the maturation of multiple new, driving technologies and mission concepts of the HLS program. As unmitigated uncertainties have previously led to failed development programs, these risks and their impacts must be understood and handled to ensure program success [2]. Sources of uncertainty include novel engine designs and configurations, increased reliance on cryogenic fluid management(CFM), and refueling technologies—which propagate as high-level performance metrics such as overall propellant mass and engine performance. Also, the occurrence of operational uncertainties—e.g. launch conditions or need to abort during the mission—can cause cascading effects on the rest of the mission that are difficult to definitively quantify, and are outside the scope of control. These concrete examples and other occurrences can be categorized as either epistemic or aleatory uncertainties.Epistemic uncertainty arises due to a lack of knowledge and can be alleviated with design and program maturation. Aleatory uncertainty is due to the inherent randomness of the system and cannot be directly reduced, unlike epistemic uncertainty. Robust design and probabilistic methods can compensate for aleatory effects. A taxonomy of uncertainty is referred to for this work [3]. In this paper, a probabilistic methodology to handle uncertainties has been demonstrated on a three-element HLS concept [1, 4], which allows tracking of current best estimates of the concept and assessment of concept design robustness against uncertainties. A sample case has been completed for this abstract, and an expansion on the methodology will be included in the final paper. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. This abstract will cover the basics of the approach used, and further details and justifications will be in the final paper.The mission profile associated with this three-element concept (Fig 1) was modeled as a set of mission events that facilitated mass changes, idles, or spacecraft maneuvers. The mission profile scope starts with each element’s NRHO orbit insertion and aggregation and ends at post-sortie rendezvous with Orion. More detail on the mission profile will be in the final paper. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as the physics framework to model the HLS architecture for applying the probabilistic methodology [5, 6]. Specifically, a parametric representation of the lander, ascent, and transfer elements and the mission profile of each element was established, with vehicle and mission parameters available as inputs to allow for a dynamic model. Each vehicle stage was modeled with high-level performance metrics, using Isp and propellant mass fraction (PMF) to remain parametric. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. These parameters were frozen at their nominal values for the purposes of baselining architecture performance and sizing the vehicle appropriately based on reference documentation [1]. Range estimating—a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities—is traditionally used with Mass Equipment Lists (MELs), but has been adapted with operational parameters as well as vehicle parameters in theDYREQT model to capture mission uncertainty alongside vehicle uncertainty [7, 3]. This method was selected due to its application and insight on a system from a bottom-up perspective, independence from historical rules of thumb, and ability to generate sensitivities based on design decisions and uncertainties. As a sample case for the abstract, the boiloff rates of the vehicle elements and the loiter times during the mission (simulating launch time variations and changing window of opportunities) were used with range estimating to provide preliminary results. To perform the range estimation portion of this methodology (depicted in Fig. 3, further details in final paper), the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the sample set of uncertainty parameters; 5,000 cases via Latin Hypercube Sampling were computed on the DYREQT architecture model. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Range estimating was applied to the surrogates instead of the actual models, which saves computational expense due to the bulk of cases needed for the Monte Carlo simulation as part of range estimating. Uncertainty parameters were sampled independently from triangular distributions using the DoE ranges as ‘min’ and ‘max’, and the nominal value as ‘most likely’. Based engineering intuition, some uncertainty parameters are correlated—e.g. if the main propellant has a high boil-off rate, the oxidizer should follow suit as both are related to CFM technology.While a Monte Carlo simulation samples all inputs as independent, the results would show model correlations; thus, it is efficient to sample the inputs as correlated. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo. A table for the DoE ranges and probability distribution parameters is shown in Table 1, and more details on Correlated Monte Carlo Simulations will be discussed in the final paper. The model’s resulting DoE showed that multivariate polynomial equations fit via least squares method captured its behavior accurately for the sample case. For the Correlated Monte Carlo Simulation, a positive correlation between fuel and oxidizer boiloff rates was used as a demonstration. 10,000 cases were computed with the surrogates and the launched masses for each vehicle element was collated. The results can be displayed in a probability density function (PDF), showing the impact of the uncertainty parameters chosen. Integrating the PDFs will yield a cumulative distribution function (CDF) that shows the cumulative probability of a given value on the x-axis. For the sample case, the elements’ launch mass margin was calculated and represented in as CDFs, as a demonstrated representation of figures of merit for the HLS concept. For the lander and ascent elements, the NRHO mass insertion limit is 16t; the transfer element has a limit of 30t [1]. It can be seen with Figure 2 that this probabilistic methodology can provide insight into mass margin with respect to the uncertainties being modeled. Currently, the results show that the lander (descent) vehicle element has the most restrictive design space; it is the only element to show a 10% probability of negative margin. Further analysis on the Monte Carlo results will show sensitivities for driving constraints and parameters for architecture feasibility, which can lead to establishing potential mission rules.The combination of range estimating with a parametric architecture model for HLS demonstrated the capability of this probabilistic methodology in a sample case. As the HLS development progresses, this methodology has the potential for keeping current best estimates of architecture performance for awarded concepts due to the flexibility in DYREQT’s modeling framework and its parametric nature. Concept maturation and increased epistemic knowledge can be injected into the model probabilistic modeling, and thus continue to track probability of mission success.

Stephanie Y Zhu↗