Results for “Cluster Dynamics Simulations” — Page 6

NASA NTRSDate not supplied

Regional and Model-Specific Response Types in A Global Gridded Crop Model Ensemble

Crop models are often employed to project crop yields under changing conditions such as global warming and associated management change for adaptation. Multi-model ensembles are promoted to enhance the robustness of projections, but questions remain on what causes often large differences between projections of individual models. Global Gridded Crop Models (GGCMs) are especially exposed to this question when applied for assessing climate change impacts, adaptation, environmental impacts of agricultural production, because their results are used in downstream analyses, such as in integrated assessment or economic modeling for projecting future land-use change. Even though global gridded crop models are often based on detailed field-scale models or have implemented similar modeling principles in other ecosystem models, global-scale models are subject to substantial uncertainties from both model structure and parametrization as well as from calibration and input data quality. AgMIP’s Global Gridded Crop Model Intercomparison (GGCMI) has thus set out to intercompare GGCMs in order to evaluate model performance, describe model uncertainties, identify inconsistencies within the ensemble and underlying reasons, and to ultimately improve models and modeling capacities. In phase 2 of the GGCMI activities, 12 modeling groups followed a modeling protocol that asked for up to 1404 31-year global simulations at 0.5 arc-degree spatial resolution to assess models’ sensitivities to changes in carbon dioxide (C; 4 different levels) temperature (T; 7 different offset levels), water supply (W; 9 levels), and nitrogen (N; 3 levels), the so-called CTWN experiment (Franke et al. 2020; http://dx.doi.org/10.5194/gmd-13-2315-2020). We here present analyses of model response types using impact response surfaces along the C, T, W, and N dimensions, respectively and collectively. Doing so, we can understand differences in simulated responses per driver rather than aggregated changes in yields. We find that models’ sensitivities to the individual driver dimensions are substantially different and often more different across models than across regions. A cluster analysis finds regional and model-specific patterns. There is some agreement across models with respect to the spatial patterns of response types but strong differences in the distribution of response type clusters across models suggests that models need to undergo further scrutiny. We suggest establishing standards in model process evaluation not only against historical dynamics but also against dedicated experiments across the CTWN dimensions.

crop models↗

NASA NTRS1978

Space Shuttle Plume Simulation Effect on Aerodynamics

Technology for simulating plumes in wind tunnel tests was not adequate to provide the required confidence in test data where plume induced aerodynamic effects might be significant. A broad research program was undertaken to correct the deficiency. Four tasks within the program are reported. Three of these tasks involve conducting experiments, related to three different aspects of the plume simulation problem: (1) base pressures; (2) lateral jet pressures; and (3) plume parameters. The fourth task involves collecting all of the base pressure test data generated during the program. Base pressures were measured on a classic cone ogive cylinder body as affected by the coaxial, high temperature exhaust plumes of a variety of solid propellant rockets. Valid data were obtained at supersonic freestream conditions but not at transonic. Pressure data related to lateral (separation) jets at M infinity = 4.5, for multiple clustered nozzles canted to the freestream and operating at high dynamic pressure ratios. All program goals were met although the model hardware was found to be large relative to the wind tunnel size so that operation was limited for some nozzle configurations.

Hair, L. M.↗

NASA NTRS2016

North American Megadroughts in the Common Era: Reconstructions and Simulations

During the Medieval Climate Anomaly (MCA), Western North America experienced episodes of intense aridity that persisted for multiple decades or longer. These megadroughts are well documented in many proxy records, but the causal mechanisms are poorly understood. General circulation models (GCMs) simulate megadroughts, but do not reproduce the temporal clustering of events during the MCA, suggesting they are not caused by the time history of volcanic or solar forcing. Instead, GCMs generate megadroughts through (1) internal atmospheric variability, (2) sea-surface temperatures, and (3) land surface and dust aerosol feedbacks. While no hypothesis has been definitively rejected, and no GCM has accurately reproduced all features (e.g., timing, duration, and extent) of any specific megadrought, their persistence suggests a role for processes that impart memory to the climate system (land surface and ocean dynamics). Over the 21st century, GCMs project an increase in the risk of megadrought occurrence through greenhouse gas forced reductions in precipitation and increases in evaporative demand. This drying is robust across models and multiple drought indicators, but major uncertainties still need to be resolved. These include the potential moderation of vegetation evaporative losses at higher atmospheric [CO2], variations in land surface model complexity, and decadal to multidecadal modes of natural climate variability that could delay or advance onset of aridification over the the next several decades. Because future droughts will arise from both natural variability and greenhouse gas forced trends in hydroclimate, improving our understanding of the natural drivers of persistent multidecadal megadroughts should be a major research priority.

atmospheric temperature↗

NASA NTRS1983

The evolution of clusters of galaxies. I - Very rich clusters

A multiple one-body Monte Carlo code is used to investigate the evolution of galaxies in a steady cluster potential under the influence of dynamical friction, two-body relaxation, tidal stripping, and galaxy mergers. The basic assumptions, estimated time scales, method and computer program, and effects of uncertainties in galaxy encounter physics are addressed. Numerical experiments in which the mass initially carried by galaxies and the mass function are varied are reported. It is shown that the formation of a very massive cluster galaxy depends critically on the number of galaxies, the initial division of cluster material between galaxies and a smooth intracluster medium, the mass spectrum of the galaxies, and chance. By the end of all simulations, less than half of the mass in the central regions of the cluster is bound to galaxies. It appears possible to produce any of the Bautz-Morgan classes from very similar initial conditions.

Richstone, D. O.↗

NASA NTRS2008

Three-Dimensional Signatures of Intermittent Magnetic Reconnection in Global Simulations of Dayside Magnetosphere Dynamics

We performed high resolution global MHD simulations of THEMIS dayside crossings events in May -June 2007. We found that magnetopause surface is not in steady-state even during extended periods of steady solar wind conditions. The so-called tilted reconnection lines become unstable due to formation of pressure bubbles, strong core field flux tubes, vortices, and traveling magnetic field cavities. The topology of FTEs differ from that in two dimension cartoons representing obliquely oriented quasi-2D flux rope. The structure of FTE is changing at spatial scales of 1 -2 Re. Closely located space probes can observe completely different signatures. Branches of bent flux rope can move in opposite directions. THEMIS and Cluster observations are consistent with signatures predicted by simulations.

Kuznetsova, M.M.↗

NASA NTRS2002

On the Formation of Massive Stars

We calculate numerically the collapse of slowly rotating, nonmagnetic, massive molecular clumps of masses 30,60, and 120 Stellar Mass, which conceivably could lead to the formation of massive stars. Because radiative acceleration on dust grains plays a critical role in the clump's dynamical evolution, we have improved the module for continuum radiation transfer in an existing two-dimensional (axial symmetry assumed) radiation hydrodynamic code. In particular, rather than using "gray" dust opacities and "gray" radiation transfer, we calculate the dust's wavelength-dependent absorption and emission simultaneously with the radiation density at each wavelength and the equilibrium temperatures of three grain components: amorphous carbon particles. silicates, and " dirty ice " -coated silicates. Because our simulations cannot spatially resolve the innermost regions of the molecular clump, however, we cannot distinguish between the formation of a dense central cluster or a single massive object. Furthermore, we cannot exclude significant mass loss from the central object(s) that may interact with the inflow into the central grid cell. Thus, with our basic assumption that all material in the innermost grid cell accretes onto a single object. we are able to provide only an upper limit to the mass of stars that could possibly be formed. We introduce a semianalytical scheme for augmenting existing evolutionary tracks of pre-main-sequence protostars by including the effects of accretion. By considering an open outermost boundary, an arbitrary amount of material could, in principal, be accreted onto this central star. However, for the three cases considered (30, 60, and 120 Stellar Mass originally within the computation grid), radiation acceleration limited the final masses to 3 1.6, 33.6, and 42.9 Stellar Mass, respectively, for wavelength-dependent radiation transfer and to 19.1, 20.1, and 22.9 Stellar Mass. for the corresponding simulations with gray radiation transfer. Our calculations demonstrate that massive stars can in principle be formed via accretion through a disk. The accretion rate onto the central source increases rapidly after one initial free-fall time and decreases monotonically afterward. By enhancing the nonisotropic character of the radiation field, the accretion disk reduces the effects of radiative acceleration in the radial direction - a process we call the "flashlight effect." The flashlight effect is further amplified in our case by including the effects of frequency-dependent radiation transfer. We conclude with the warning that a careful treatment of radiation transfer is a mandatory requirement for realistic simulations of the formation of massive stars.

Yorke, Harold W.↗

NASA NTRS2003

StarBooster Demonstrator Cluster Configuration Analysis/Verification Program

In order to study the flight dynamics of the cluster configuration of two first stage boosters and upper-stage, flight-testing of subsonic sub-scale models has been undertaken using two glideback boosters launched on a center upper-stage. Three high power rockets clustered together were built and flown to demonstrate vertical launch, separation and horizontal recovery of the boosters. Although the boosters fly to conventional aircraft landing, the centerstage comes down separately under its own parachute. The goal of the project has been to collect data during separation and flight for comparison with a six degree of freedom simulation. The configuration for the delta wing canard boosters comes from a design by Starcraft Boosters, Inc. The subscale rockets were constructed of foam covered in carbon or fiberglass and were launched with commercially available solid rocket motors. The first set of boosters built were 3-ft tall with a 4-ft tall centerstage, and two additional sets of boosters were made that were each over 5-ft tall with a 7.5 ft centerstage. The rocket cluster is launched vertically, then after motor bum out the boosters are separated and flown to a horizontal landing under radio-control. An on-board data acquisition system recorded data during both the launch and glide phases of flight.

DeTurris, Dianne J.↗

NASA NTRS2008

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗

NASA NTRS2010

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗

NASA NTRS1993

Using x ray images to detect substructure in a sample of 40 Abell clusters

Using a method for constraining the dynamical state of a galaxy cluster by examining the moments of its x-ray surface brightness distribution, we determine the statistics of cluster substructure for a sample of 40 Abell clusters. Using x-ray observations from the Einstein Observatory Imaging Proportional Counter (IPC), we measure the first moment M1(r), the ellipsoidal orientation angle theta2(r), and the axial ratio eta(r) at several different radii in the cluster. We determine the effects of systematics such as x-ray point source emission, telescope vignetting, Poisson noise, and characteristics of the IPC by measuring the same parameters on an ensemble of simulated cluster images. Due to the small band-pass of the IPC, the ICM emissivity is nearly independent of temperature so the intensity at each point in the IPC images is simply proportional to the emission measure calculated along the line of sight through the cluster (e.g. Fabricant et al. 1980). Therefore, barring a change superposition of two x-ray emitting clusters, a significant variation in the image centroid M1(r) as a function of radius indicates that the center of mass of the intra-cluster medium (ICM) varies with radius. We argue that such a configuration (essentially an m = 1 component in the ICM density distribution) is a non-equilibrium component; it results from an off-center subclump or a recent merger in the ICM.

Mohr, J. J.↗

NASA NTRS1990

Dynamics of groups around interacting double ellipticals: Measuring dark matter haloes

Binary galaxies, as binary stars, are important to measure masses, as suggested by Page (1952). Because three orbit parameters are measurable for galaxies at one instant of time, severe uncertainties remain in the orbit and mass determinations. These uncertainties can partly be overcome by statistical studies of selected samples and/or n-body simulations. Close double galaxies (and isolated galaxies) could also be useful to estimate dynamical masses if we can find test particles around them. Interacting elliptical pairs or dumb-bell galaxies are found with a large range, between 0-1200 km s(exp -1), of relative radial velocities. Standard 2-body orbit calculations, highly uncertain due to projection factors, suggest for the largest velocity differences very large galaxy masses, if the systems are bound and stationary. However, recent n-body simulations model these binaries as galaxies captured from hyperbolic orbits, requiring masses of order a few times 10(exp 11) solar maximum (Borne et al. 1988), but producing systems that are short lived. A different picture appears when we study observationally the dynamical mass of interacting double ellipticals using faint satellite galaxies. These satellites contribute little luminosity and, presumably, little mass to the system. The authors present results of two such groups, basically forming systems of test particles, around the dumb-bells NGC 4782/3 and IC 5049. They also briefly discuss the satellite group around the central dumb-bell in the cluster Sersic 40/6. Apparently, they detect large quantities of dark matter in the vicinity of these dumb-bell galaxies, because the system masses of approx. 4.5 times 10(exp 13) solar mass and 8 times 10(exp 13) solar mass for NGC 4782/3 and IC 5049, respectively, are quite high. Likewise, the mass of the Sersic 40/6 inner core is 7 times 10(exp 13) solar mass. The possibility that a common massive dark matter halo increases the merging times of these types of galaxies is suggested. Here, the authors assume H sub o = 100 km s (exp -1) Mpc(exp -1).

Quintana, H.↗

NASA NTRS2004

Chemical Reaction and Flow Modeling in Fullerene and Nanotube Production

The development of processes to produce fullerenes and carbon nanotubes has largely been empirical. Fullerenes were first discovered in the soot produced by laser ablation of graphite [1]and then in the soot of electric arc evaporated carbon. Techniques and conditions for producing larger and larger quantities of fullerenes depended mainly on trial and error empirical variations of these processes, with attempts to scale them up by using larger electrodes and targets and higher power. Various concepts of how fullerenes and carbon nanotubes were formed were put forth, but very little was done based on chemical kinetics of the reactions. This was mainly due to the complex mixture of species and complex nature of conditions in the reactors. Temperatures in the reactors varied from several thousand degrees Kelvin down to near room temperature. There are hundreds of species possible, ranging from atomic carbon to large clusters of carbonaceous soot, and metallic catalyst atoms to metal clusters, to complexes of metals and carbon. Most of the chemical kinetics of the reactions and the thermodynamic properties of clusters and complexes have only been approximated. In addition, flow conditions in the reactors are transient or unsteady, and three dimensional, with steep spatial gradients of temperature and species concentrations. All these factors make computational simulations of reactors very complex and challenging. This article addresses the development of the chemical reaction involved in fullerene production and extends this to production of carbon nanotubes by the laser ablation/oven process and by the electric arc evaporation process. In addition, the high-pressure carbon monoxide (HiPco) process is discussed. The article is in several parts. The first one addresses the thermochemical aspects of modeling; and considers the development of chemical rate equations, estimates of reaction rates, and thermodynamic properties where they are available. The second part addresses modeling of the arc process for fullerene and carbon nanotube production using O-D, 1-D and 2-D fluid flow models. The third part addresses simulations of the pulsed laser ablation process using time-dependent techniques in 2-D, and a steady state 2-D simulation of a continuous laser ablation process. The fourth part addresses steady state modeling in O-D and 2-D of the HiPco process. In each of the simulations, there is a variety of simplifications that are made that enable one to concentrate on one aspect or another of the process. There are simplifications that can be made to the chemical reaction models , e.g. reduction in number of species by lumping some of them together in a representative species. Other simulations are carried out by eliminating the chemistry altogether in order to concentrate on the fluid dynamics. When solving problems with a large number of species in more than one spatial dimension, it is almost imperative that the problem be decoupled by solving for the fluid dynamics to find the fluid motion and temperature history of "particles" of fluid moving through a reactor. Then one can solve the chemical rate equations with complex chemistry following the temperature and pressure history. One difficulty is that often mixing with an ambient gas is involved. Therefore, one needs to take dilution and mixing into account. This changes the ratio of carbon species to background gas. Commercially available codes may have no provision for including dilution as part of the input. One must the write special solvers for including dilution in decoupled problems. The article addresses both ful1erene production and single-walled carbon nanotube (SWNT) production. There are at least two schemes or concepts of SWNT growth. This article will only address growth in the gas phase by carbon and catalyst cluster growth and SW T formation by the addition of carbon. There are other models that conceive of SWNT growth as a phase separation process from clusters me up carbon and metal catalyst, with the carbon precipitating from the cluster as it cools. We will not deal with that concept in this article. Further research is needed to determine the rates at which these composite clusters form, evaporate, and segregate.

Scott, Carl D.↗

NASA NTRS1993

Testing higher-order Lagrangian perturbation theory against numerical simulation. 1: Pancake models

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasi-linear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the Zel'dovich approximation (hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of sigma is approximately 2). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-spectrum with power-index n = -1) using cross-correlation statistics employed in previous work. We found that for all statistical methods used the higher-order corrections improve the results obtained for the first-order solution up to the stage when sigma (linear theory) is approximately 1. While this improvement can be seen for all spatial scales, later stages retain this feature only above a certain scale which is increasing with time. However, third-order is not much improvement over second-order at any stage. The total breakdown of the perturbation approach is observed at the stage, where sigma (linear theory) is approximately 2, which corresponds to the onset of hierarchical clustering. This success is found at a considerable higher non-linearity than is usual for perturbation theory. Whether a truncation of the initial power-spectrum in hierarchical models retains this improvement will be analyzed in a forthcoming work.

Buchert, T.↗

NASA NTRS1976

Ultrahigh-accuracy body-pointing system for the Large Space Telescope

The Large Space Telescope (LST) is a 3-m diffraction-limited telescope. Pointing stability requirements necessary to assure diffraction-limited images are plus or minus 0.005 arc-sec, over possible experiment observation times of several hours. In order to determine whether these stringent pointing requirements could be met, a complex simulation model was defined which consisted of detailed dynamic representations of control moment gyros (CMGs) and reaction wheels (RWs), including their noise characteristics, dynamic sensor representations with sensor noise, shock mounts for the CMG actuators, a detailed representation of an image motion compensation (IMC) system, and a detailed flexible body structural model with all significant vehicle and solar panel bending modes. On the basis of both stability and performance studies utilizing this model, it was determined that a body-pointing system will meet LST requirements in the presence of CMG vibrational disturbances and sensor noise. The recommended system consists of three orthogonally mounted RWs for primary short-term control, and a cluster of CMG actuators for continuous RW desaturation and vehicle maneuvering.

Rybak, S. C.↗

NASA NTRS2010

Toward GEOS-6, A Global Cloud System Resolving Atmospheric Model

NASA is committed to observing and understanding the weather and climate of our home planet through the use of multi-scale modeling systems and space-based observations. Global climate models have evolved to take advantage of the influx of multi- and many-core computing technologies and the availability of large clusters of multi-core microprocessors. GEOS-6 is a next-generation cloud system resolving atmospheric model that will place NASA at the forefront of scientific exploration of our atmosphere and climate. Model simulations with GEOS-6 will produce a realistic representation of our atmosphere on the scale of typical satellite observations, bringing a visual comprehension of model results to a new level among the climate enthusiasts. In preparation for GEOS-6, the agency's flagship Earth System Modeling Framework [JDl] has been enhanced to support cutting-edge high-resolution global climate and weather simulations. Improvements include a cubed-sphere grid that exposes parallelism; a non-hydrostatic finite volume dynamical core, and algorithm designed for co-processor technologies, among others. GEOS-6 represents a fundamental advancement in the capability of global Earth system models. The ability to directly compare global simulations at the resolution of spaceborne satellite images will lead to algorithm improvements and better utilization of space-based observations within the GOES data assimilation system

Putman, William M.↗

NASA NTRS2011

Graphics Processing Unit (GPU) Acceleration of the Goddard Earth Observing System Atmospheric Model

The Goddard Earth Observing System 5 (GEOS-5) is the atmospheric model used by the Global Modeling and Assimilation Office (GMAO) for a variety of applications, from long-term climate prediction at relatively coarse resolution, to data assimilation and numerical weather prediction, to very high-resolution cloud-resolving simulations. GEOS-5 is being ported to a graphics processing unit (GPU) cluster at the NASA Center for Climate Simulation (NCCS). By utilizing GPU co-processor technology, we expect to increase the throughput of GEOS-5 by at least an order of magnitude, and accelerate the process of scientific exploration across all scales of global modeling, including: The large-scale, high-end application of non-hydrostatic, global, cloud-resolving modeling at 10- to I-kilometer (km) global resolutions Intermediate-resolution seasonal climate and weather prediction at 50- to 25-km on small clusters of GPUs Long-range, coarse-resolution climate modeling, enabled on a small box of GPUs for the individual researcher After being ported to the GPU cluster, the primary physics components and the dynamical core of GEOS-5 have demonstrated a potential speedup of 15-40 times over conventional processor cores. Performance improvements of this magnitude reduce the required scalability of 1-km, global, cloud-resolving models from an unfathomable 6 million cores to an attainable 200,000 GPU-enabled cores.

Putnam, Williama↗

NASA NTRS2000

Shape-dependent control of cell growth, differentiation, and apoptosis: switching between attractors in cell regulatory networks

Development of characteristic tissue patterns requires that individual cells be switched locally between different phenotypes or "fates;" while one cell may proliferate, its neighbors may differentiate or die. Recent studies have revealed that local switching between these different gene programs is controlled through interplay between soluble growth factors, insoluble extracellular matrix molecules, and mechanical forces which produce cell shape distortion. Although the precise molecular basis remains unknown, shape-dependent control of cell growth and function appears to be mediated by tension-dependent changes in the actin cytoskeleton. However, the question remains: how can a generalized physical stimulus, such as cell distortion, activate the same set of genes and signaling proteins that are triggered by molecules which bind to specific cell surface receptors. In this article, we use computer simulations based on dynamic Boolean networks to show that the different cell fates that a particular cell can exhibit may represent a preprogrammed set of common end programs or "attractors" which self-organize within the cell's regulatory networks. In this type of dynamic network model of information processing, generalized stimuli (e.g., mechanical forces) and specific molecular cues elicit signals which follow different trajectories, but eventually converge onto one of a small set of common end programs (growth, quiescence, differentiation, apoptosis, etc.). In other words, if cells use this type of information processing system, then control of cell function would involve selection of preexisting (latent) behavioral modes of the cell, rather than instruction by specific binding molecules. Importantly, the results of the computer simulation closely mimic experimental data obtained with living endothelial cells. The major implication of this finding is that current methods used for analysis of cell function that rely on characterization of linear signaling pathways or clusters of genes with common activity profiles may overlook the most critical features of cellular information processing which normally determine how signal specificity is established and maintained in living cells. Copyright 2000 Academic Press.

Review↗

NASA NTRS2003

The Ionization History of The Intergalactic Medium:

The funded project seeked a unified description of the ionization, physical structure, and evolution of the intergalactic medium (IGM) and quasar intervening absorption systems. We proposed to conduct theoretical studies of the IGM and QSO absorbers in the context of current theories of galaxy formation, developing and using numerical and analytical techniques aimed at a detailed modeling of cosmological radiative transfer, gas dynamics, and thermal and ionization evolution. The ionization history of the IGM has important implications for the metagalactic UV background, intergalactic helium absorption 21-cm tomography, metal absorption systems, fluctuations in the microwave background, and the cosmic rate of structure and star formation. All the original objectives of our program have been achieved, and the results widely used and quoted by the community. Indeed, they remain relevant as the level and complexity of research in this area has increased substantially since our proposal was submitted, due to new discoveries on galaxy formation and evolution, a flood of high-quality data on the distant universe, new theoretical ideas and direct numerical simulations of structure formation in hierarchical clustering theories.

Madau, Piero↗