Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel cluster”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

Sub-problem Optimization With Regression and Neural Network Approximators

Design optimization of large systems can be attempted through a sub-problem strategy. In this strategy, the original problem is divided into a number of smaller problems that are clustered together to obtain a sequence of sub-problems. Solution to the large problem is attempted iteratively through repeated solutions to the modest sub-problems. This strategy is applicable to structures and to multidisciplinary systems. For structures, clustering the substructures generates the sequence of sub-problems. For a multidisciplinary system, individual disciplines, accounting for coupling, can be considered as sub-problems. A sub-problem, if required, can be further broken down to accommodate sub-disciplines. The sub-problem strategy is being implemented into the NASA design optimization test bed, referred to as "CometBoards." Neural network and regression approximators are employed for reanalysis and sensitivity analysis calculations at the sub-problem level. The strategy has been implemented in sequential as well as parallel computational environments. This strategy, which attempts to alleviate algorithmic and reanalysis deficiencies, has the potential to become a powerful design tool. However, several issues have to be addressed before its full potential can be harnessed. This paper illustrates the strategy and addresses some issues.

Guptill, James D.↗

Hot Chips and Hot Interconnects for High End Computing Systems

I will discuss several processors: 1. The Cray proprietary processor used in the Cray X1; 2. The IBM Power 3 and Power 4 used in an IBM SP 3 and IBM SP 4 systems; 3. The Intel Itanium and Xeon, used in the SGI Altix systems and clusters respectively; 4. IBM System-on-a-Chip used in IBM BlueGene/L; 5. HP Alpha EV68 processor used in DOE ASCI Q cluster; 6. SPARC64 V processor, which is used in the Fujitsu PRIMEPOWER HPC2500; 7. An NEC proprietary processor, which is used in NEC SX-6/7; 8. Power 4+ processor, which is used in Hitachi SR11000; 9. NEC proprietary processor, which is used in Earth Simulator. The IBM POWER5 and Red Storm Computing Systems will also be discussed. The architectures of these processors will first be presented, followed by interconnection networks and a description of high-end computer systems based on these processors and networks. The performance of various hardware/programming model combinations will then be compared, based on latest NAS Parallel Benchmark results (MPI, OpenMP/HPF and hybrid (MPI + OpenMP). The tutorial will conclude with a discussion of general trends in the field of high performance computing, (quantum computing, DNA computing, cellular engineering, and neural networks).

Saini, Subhash↗

Parallelization of a Six Degree of Freedom Entry Vehicle Trajectory Simulation Using OpenMP and OpenACC

The art and science of writing parallelized software, using methods such as Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC), is dominated by computer scientists. Engineers and non-computer scientists looking to apply these techniques to their project applications face a steep learning curve, especially when looking to adapt their original single threaded software to run multi-threaded on graphics processing units (GPUs). There are significant changes in mindset that must occur; such as how to manage memory, the organization of instructions, and the use of if statements (also known as branching). The purpose of this work is twofold: 1) to demonstrate the applicability of parallelized coding methodologies, OpenMP and OpenACC, to tasks outside of the typical large scale matrix mathematics; and 2) to discuss, from an engineer’s perspective, the lessons learned from parallelizing software using these computer science techniques. This work applies OpenMP, on both multi-core central processing units (CPUs) and Intel® Xeon Phi™ 7210, and OpenACC on GPUs. These parallelization techniques are used to tackle the simulation of thousands of entry vehicle trajectories through the integration of six degree of freedom (DoF) equations of motion (EoM). The forces and moments acting on the entry vehicle, and used by the EoM, are estimated using multiple models of varying levels of complexity. Several benchmark comparisons are made on the execution of six DoF trajectory simulation: single thread Intel® Xeon® E5-2670 CPU, multi-thread CPU using OpenMP, multi-thread Xeon Phi™ 7210 using OpenMP, and multi-thread NVIDIA® Tesla® K40 GPU using OpenACC. These benchmarks are run on the Pleiades Supercomputer Cluster at the National Aeronautics and Space Administration (NASA) Ames Research Center (ARC), and a Xeon Phi™ 7210 node at NASA Langley Research Center (LaRC).

Green, Justin S.↗

An Extended Lagrangian Method

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method,' is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. The present method and the Arbitrary Lagrangian-Eulerian (ALE) method have a similarity in spirit-eliminating the cross-streamline numerical diffusion. For this purpose, we suggest a simple grid constraint condition and utilize an accurate discretization procedure. This grid constraint is only applied to the transverse cell face parallel to the local stream velocity, and hence our method for the steady state problems naturally reduces to the streamline-curvature method, without explicitly solving the steady stream-coordinate equations formulated a priori. Unlike the Lagrangian method proposed by Loh and Hui which is valid only for steady supersonic flows, the present method is general and capable of treating subsonic flows and supersonic flows as well as unsteady flows, simply by invoking in the same code an appropriate grid constraint suggested in this paper. The approach is found to be robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.

Liou, Meng-Sing↗

New Insight into Short-Wavelength Solar Wind Fluctuations from Vlasov Theory

The nature of solar wind (SW) turbulence below the proton gyroscale is a topic that is being investigated extensively nowadays, both theoretically and observationally. Although recent observations gave evidence of the dominance of kinetic Alfven waves (KAWs) at sub-ion scales with omega < omega(sub ci), other studies suggest that the KAW mode cannot carry the turbulence cascade down to electron scales and that the whistler mode (i.e., omega > omega (sub ci)) is more relevant. Here, we study key properties of the short-wavelength plasma modes under limited, but realistic, SW conditions, Typically Beta(sub i) approx. > Beta (sub e) ∼ 1 and for high oblique angles of propagation 80 deg <= Theta (sub kB) < 90 deg as observed from the Cluster spacecraft data. The linear properties of the plasma modes under these conditions are poorly known, which contrasts with the well-documented cold plasma limit and/or moderate oblique angles of propagation (Theta (sub kB) < 80 deg). Based on linear solutions of the Vlasov kinetic theory, we discuss the relevance of each plasma mode (fast, Bernstein, KAW, whistler) in carrying the energy cascade down to electron scales. We show, in particular, that the shear Alfven mode (known in the magnetohydrodynamic limit) extends at scales kappa rho (sub i) approx. > 1 to frequencies either larger or smaller than omega (sub ci), depending on the anisotropy kappa (parallel )/ kappa(perpendicular). This extension into small scales is more readily called whistler (omega > omega (sub ci)) or KAW (omega < omega (sub ci)) although the mode is essentially the same. This contrasts with the well-accepted idea that the whistler branch always develops as a continuation at high frequencies of the fast magnetosonic mode. We show, furthermore, that the whistler branch is more damped than the KAW one, which makes the latter the more relevant candidate to carry the energy cascade down to electron scales. We discuss how these new findings may facilitate resolution of the controversy concerning the nature of the small-scale turbulence, and we discuss the implications for present and future spacecraft wave measurements in the SW.

waves↗

Peripheral sensory processing in mammalian gravity receptors - Observations of ciliary tuft configurations

Scanning electron microscopy was used to study dynamic polarizations of clustered cells of the anterior part of rat saccular macula and to shed light on the possible roles of two types of hair cells integrated into the same neural circuitry: those with short stereocilia and long kinocilium (ss/lk), and those with long stereocilia and still longer kinocilium (ls/lk). It was found that the ss/lk-type cells could be further subdivided into two types, whereas the ls/lk cells consisted of four major kinds. It was also found that the kinocilium was most often fixed in a recovery stroke position (curved basally, and the upper portion projected back over the tuft) and that the kinocilia were not aligned in parallel in any given part of a macula, even though each cilium pointed in the proper direction relative to the striola line. The possibility of a relationship between the ciliary tuft morphology and the function of the hair cell of which it is a part is discussed.

Ross, Muriel D.↗

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.↗

The Analysis of the Patterns of Radiation-Induced DNA Damage Foci by a Stochastic Monte Carlo Model of DNA Double Strand Breaks Induction by Heavy Ions and Image Segmentation Software

To create a generalized mechanistic model of DNA damage in human cells that will generate analytical and image data corresponding to experimentally observed DNA damage foci and will help to improve the experimental foci yields by simulating spatial foci patterns and resolving problems with quantitative image analysis. Material and Methods: The analysis of patterns of RIFs (radiation-induced foci) produced by low- and high-LET (linear energy transfer) radiation was conducted by using a Monte Carlo model that combines the heavy ion track structure with characteristics of the human genome on the level of chromosomes. The foci patterns were also simulated in the maximum projection plane for flat nuclei. Some data analysis was done with the help of image segmentation software that identifies individual classes of RIFs and colocolized RIFs, which is of importance to some experimental assays that assign DNA damage a dual phosphorescent signal. Results: The model predicts the spatial and genomic distributions of DNA DSBs (double strand breaks) and associated RIFs in a human cell nucleus for a particular dose of either low- or high-LET radiation. We used the model to do analyses for different irradiation scenarios. In the beam-parallel-to-the-disk-of-a-flattened-nucleus scenario we found that the foci appeared to be merged due to their high density, while, in the perpendicular-beam scenario, the foci appeared as one bright spot per hit. The statistics and spatial distribution of regions of densely arranged foci, termed DNA foci chains, were predicted numerically using this model. Another analysis was done to evaluate the number of ion hits per nucleus, which were visible from streaks of closely located foci. In another analysis, our image segmentaiton software determined foci yields directly from images with single-class or colocolized foci. Conclusions: We showed that DSB clustering needs to be taken into account to determine the true DNA damage foci yield, which helps to determine the DSB yield. Using the model analysis, a researcher can refine the DSB yield per nucleus per particle. We showed that purely geometric artifacts, present in the experimental images, can be analytically resolved with the model, and that the quantization of track hits and DSB yields can be provided to the experimentalists who use enumeration of radiation-induced foci in immunofluorescence experiments using proteins that detect DNA damage. An automated image segmentaiton software can prove useful in a faster and more precise object counting for colocolized foci images.

Ponomarev, Artem↗

Large-Scale NASA Science Applications on the Columbia Supercluster

Columbia, NASA's newest 61 teraflops supercomputer that became operational late last year, is a highly integrated Altix cluster of 10,240 processors, and was named to honor the crew of the Space Shuttle lost in early 2003. Constructed in just four months, Columbia increased NASA's computing capability ten-fold, and revitalized the Agency's high-end computing efforts. Significant cutting-edge science and engineering simulations in the areas of space and Earth sciences, as well as aeronautics and space operations, are already occurring on this largest operational Linux supercomputer, demonstrating its capacity and capability to accelerate NASA's space exploration vision. The presentation will describe how an integrated environment consisting not only of next-generation systems, but also modeling and simulation, high-speed networking, parallel performance optimization, and advanced data analysis and visualization, is being used to reduce design cycle time, accelerate scientific discovery, conduct parametric analysis of multiple scenarios, and enhance safety during the life cycle of NASA missions. The talk will conclude by discussing how NAS partnered with various NASA centers, other government agencies, computer industry, and academia, to create a national resource in large-scale modeling and simulation.

Brooks, Walter↗

Fast Multipole Methods for Three-Dimensional N-body Problems

We are developing computational tools for the simulations of three-dimensional flows past bodies undergoing arbitrary motions. High resolution viscous vortex methods have been developed that allow for extended simulations of two-dimensional configurations such as vortex generators. Our objective is to extend this methodology to three dimensions and develop a robust computational scheme for the simulation of such flows. A fundamental issue in the use of vortex methods is the ability of employing efficiently large numbers of computational elements to resolve the large range of scales that exist in complex flows. The traditional cost of the method scales as Omicron (N(sup 2)) as the N computational elements/particles induce velocities at each other, making the method unacceptable for simulations involving more than a few tens of thousands of particles. In the last decade fast methods have been developed that have operation counts of Omicron (N log N) or Omicron (N) (referred to as BH and GR respectively) depending on the details of the algorithm. These methods are based on the observation that the effect of a cluster of particles at a certain distance may be approximated by a finite series expansion. In order to exploit this observation we need to decompose the element population spatially into clusters of particles and build a hierarchy of clusters (a tree data structure) - smaller neighboring clusters combine to form a cluster of the next size up in the hierarchy and so on. This hierarchy of clusters allows one to determine efficiently when the approximation is valid. This algorithm is an N-body solver that appears in many fields of engineering and science. Some examples of its diverse use are in astrophysics, molecular dynamics, micro-magnetics, boundary element simulations of electromagnetic problems, and computer animation. More recently these N-body solvers have been implemented and applied in simulations involving vortex methods. Koumoutsakos and Leonard (1995) implemented the GR scheme in two dimensions for vector computer architectures allowing for simulations of bluff body flows using millions of particles. Winckelmans presented three-dimensional, viscous simulations of interacting vortex rings, using vortons and an implementation of a BH scheme for parallel computer architectures. Bhatt presented a vortex filament method to perform inviscid vortex ring interactions, with an alternative implementation of a BH scheme for a Connection Machine parallel computer architecture.

Koumoutsakos, P.↗

Numerical Propulsion System Simulation: An Overview

The cost of implementing new technology in aerospace propulsion systems is becoming prohibitively expensive and time consuming. One of the main contributors to the high cost and lengthy time is the need to perform many large-scale hardware tests and the inability to integrate all appropriate subsystems early in the design process. The NASA Glenn Research Center is developing the technologies required to enable simulations of full aerospace propulsion systems in sufficient detail to resolve critical design issues early in the design process before hardware is built. This concept, called the Numerical Propulsion System Simulation (NPSS), is focused on the integration of multiple disciplines such as aerodynamics, structures and heat transfer with computing and communication technologies to capture complex physical processes in a timely and cost-effective manner. The vision for NPSS, as illustrated, is to be a "numerical test cell" that enables full engine simulation overnight on cost-effective computing platforms. There are several key elements within NPSS that are required to achieve this capability: 1) clear data interfaces through the development and/or use of data exchange standards, 2) modular and flexible program construction through the use of object-oriented programming, 3) integrated multiple fidelity analysis (zooming) techniques that capture the appropriate physics at the appropriate fidelity for the engine systems, 4) multidisciplinary coupling techniques and finally 5) high performance parallel and distributed computing. The current state of development in these five area focuses on air breathing gas turbine engines and is reported in this paper. However, many of the technologies are generic and can be readily applied to rocket based systems and combined cycles currently being considered for low-cost access-to-space applications. Recent accomplishments include: (1) the development of an industry-standard engine cycle analysis program and plug 'n play architecture, called NPSS Version 1, (2) A full engine simulation that combines a 3D low-pressure subsystem with a 0D high pressure core simulation. This demonstrates the ability to integrate analyses at different levels of detail and to aerodynamically couple components, the fan/booster and low-pressure turbine, through a 3D computational fluid dynamics simulation. (3) Simulation of all of the turbomachinery in a modern turbofan engine on parallel computing platform for rapid and cost-effective execution. This capability can also be used to generate full compressor map, requiring both design and off-design simulation. (4) Three levels of coupling characterize the multidisciplinary analysis under NPSS: loosely coupled, process coupled and tightly coupled. The loosely coupled and process coupled approaches require a common geometry definition to link CAD to analysis tools. The tightly coupled approach is currently validating the use of arbitrary Lagrangian/Eulerian formulation for rotating turbomachinery. The validation includes both centrifugal and axial compression systems. The results of the validation will be reported in the paper. (5) The demonstration of significant computing cost/performance reduction for turbine engine applications using PC clusters. The NPSS Project is supported under the NASA High Performance Computing and Communications Program.

Lytle, John K.↗

The Omics of Stem Cell Mediated Regeneration: A Pilot Single Cell RNA-Seq Study of Mechanotransduction

Mechanical forces are potent modulators of stem cell based tissue regenerative mechanisms, inducing cell fate decisions and tissue specific commitment. A unique platform for investigating mechanotransduction is spaceflight, where microgravity and altered fluid mechanics provide a loading-null experimental condition. Seminal investigations of regenerative capacity in a wholly regenerative species, the newt model, and in a variety of totipotent and adult stem cell populations have demonstrated the detrimental effects of unloading on maintenance of stem cell based regeneration. Of particular interest is the observation that unloading interferes with the transition of stem cell pools from proliferative state to differentiation commitment. In this work we sought to test the hypothesis that gravity mechanotransduction regulates stem cell tissue regenerative processes by modulating stem cell proliferation and differentiation fates at specific cell cycle stages. To do this, clonally-derived ESCs were plated on a collagen matrix and expanded for 36 hours before re-plating on a non-adherent culture dish in the absence of leukemia inhibitory factor (LIF) to form spheroid aggregate EBs. After formation, the EBs were transferred to a collagen matrix coated culture dishes and given 4 days to allow implantation and outgrowth. In parallel, totipotent ESCs were plated 24 hours before mechanical stimulation on collagen matrix culture dishes in the presence of LIF to maintain totipotency and serve as un-differentiation committed controls. The EBs and ESCs were then subjected to either a 60 minute pulse of gravity (static loading) or 60 minutes of cyclic stretch (dynamic loading) mechanotransduction. Six hours post-stimulation, we used a 10X Genomics Single Cell controller to generate bar-coded single cell Illumina libraries and sequenced expressomes for 5,000 static loaded cells, representative of a change in gravity mechanotransduction, 5,000 dynamic loaded cells, representative of tissue loading associate with physiologic function, and 5,000 unstimulated 1g control cells. The comparison of these 3 libraries by cluster assignment based on like gene expression patterns show substantial alteration in cluster geometry due to mechanical loading. Specifically the mechanically loaded EB outgrowth cells to retain potency markers (PAX6, SOX2, CD34) and suppress early commitment markers (Dhh, VCAN, Igf1). Whereas the EBs cultured under the non-stimulated conditions display clear departure from the ESC expressome with lineage commitment markers upregulated and several tissue specific markers being expressed (BMP "early musculoskeletal development, Mesp1" early cardiovascular cell lineage). These markers are not seen in the mechano-stimulated cultures or the totipotent ESC cultures. Comparison of like clusters between our experimental conditions revealed an array of regenerative and stem cell genes are significantly mechano-regulated. Of particular importance CDKN1a/p21, a gene shown by previous investigation of our research team to be significantly upregulated in unloading, was suppressed in the static and dynamic loaded EBS. In addition to CDKN1a/p21 many genes related to cell cycle and transitory differentiation markers had elevated expression in the mechano-stimulated EBs, but surprisingly these trends were not observed in the ESC cultures. This study is the first of its kind investigating for mechano-signaling and mechano-regulated pathways, and has alre

Single cell RNA-sequecing↗

Investigating the Feedback Path in a Jet-Surface Resonant Interaction

A resonant interaction between an 8:1 aspect ratio rectangular jet and flat-plates, placed parallel to the jet, is studied experimentally. For certain locations of the plate relative to the jet, the resonance takes place with a loud accompanying tone. The sound pressure level spectra are often marked by multiple peaks. The frequencies of the spectral peaks are studied as a function of the streamwise length of the plate, its relative location to the jet as well as the jet Mach number. It is demonstrated that the tones are not due to a simple feedback between the plate's trailing edge and the nozzle's exit; the leading edge of the plate also comes into play in the frequency selection. With parametric variation, it is found that there is an order in the most energetic spectral peaks; their frequencies cluster in distinct bands. The 'fundamental', i.e., the lowest frequency band is explained by an acoustic feedback involving diffraction at the plate's leading edge.

Resonance↗

Crater Flux Transfer Events: Highroad to the X Line?

We examine Cluster observations of a so-called magnetosphere crater FTE, employing data from five instruments (FGM, CIS, EDI, EFW, and WHISPER), some at the highest resolution. The aim of doing this is to deepen our understanding of the reconnection nature of these events by applying recent advances in the theory of collisionless reconnection and in detailed observational work. Our data support the hypothesis of a stratified structure with regions which we show to be spatial structures. We support the bulge-like topology of the core region (R3) made up of plasma jetting transverse to reconnected field lines. We document encounters with a magnetic separatrix as a thin layer embedded in the region (R2) just outside the bulge, where the speed of the protons flowing approximately parallel to the field maximizes: (1) short (fraction of a sec) bursts of enhanced electric field strengths (up to approximately 30 mV/m) and (2) electrons flowing against the field toward the X line at approximately the same time as the bursts of intense electric fields. R2 also contains a density decrease concomitant with an enhanced magnetic field strength. At its interface with the core region, R3, electric field activity ceases abruptly. The accelerated plasma flow profile has a catenary shape consisting of beams parallel to the field in R2 close to the R2/R3 boundary and slower jets moving across the magnetic field within the bulge region. We detail commonalities our observations of crater FTEs have with reconnection structures in other scenarios. We suggest that in view of these properties and their frequency of occurrence, crater FTEs are ideal places to study processes at the separatrices, key regions in magnetic reconnection. This is a good preparation for the MMS mission.

Farrugia, C. J.↗

Stereo-Based Region-Growing using String Matching

We present a novel stereo algorithm based on a coarse texture segmentation preprocessing phase. Matching is performed using a string comparison. Matching sub-strings correspond to matching sequences of textures. Inter-scanline clustering of matching sub-strings yields regions of matching texture. The shape of these regions yield information concerning object's height, width and azimuthal position relative to the camera pair. Hence, rather than the standard dense depth map, the output of this algorithm is a segmentation of objects in the scene. Such a format is useful for the integration of stereo with other sensor modalities on a mobile robotic platform. It is also useful for localization; the height and width of a detected object may be used for landmark recognition, while depth and relative azimuthal location determine pose. The algorithm does not rely on the monotonicity of order of image primitives. Occlusions, exposures, and foreshortening effects are not problematic. The algorithm can deal with certain types of transparencies. It is computationally efficient, and very amenable to parallel implementation. Further, the epipolar constraints may be relaxed to some small but significant degree. A version of the algorithm has been implemented and tested on various types of images. It performs best on random dot stereograms, on images with easily filtered backgrounds (as in synthetic images), and on real scenes with uncontrived backgrounds.

Mandelbaum, Robert↗

Africa's Megafans and Their Tectonic Setting

Megafans are a really extensive continental sediment bodies, fluvially derived, and fan-shaped in planform. Only those >80 km long were included in this study. Africa's megafans were mapped for purposes of both comprehensive geomorphic description and as a method of mapping by remote sensing large probable fluvial sediment bodies (we exclude sediment bodies deposited in well defined, modern floodplains and coastal deltas). Our criteria included a length dimension of >80 km and maximum width >40 km, partial cone morphology, and a radial drainage pattern. Visible and especially IR imagery were used to identify the features, combined with topographic SRTM data. We identified 99 megafans most of which are unstudied thus far. Their feeder rivers responsible for depositing megafan sediments rise on, and are consequent drainages oriented down the slopes of the swells that have dominated African landscapes since approximately 34 Ma (the high points in Africa's so-called basin-and-swell topography [1]). Most megafans (66%) have developed along these consequent rivers relatively near the swell cores, oriented radially away from the swells. The vast basins between the swells provide accommodation for megafan sediment wedges. Although clearly visible remotely, most megafans are inactive as a result of incision by the feeder river (which then no longer operates on the fan surface). Two tectonic settings control the location of Africa's megafans, 66% on swell flanks, and 33% related to rifts. (i) Swell flanks Most megafans are apexed relatively near the core of the parent swell, and are often clustered in groups: e.g., six on the west and north flanks of the Hoggar Swell (Algeria), seven on the north and south flanks of the Tibesti Swell (Libya-Chad borderlands), twelve on the west flank of the Ethiopian Swell, four on the east flank of the East African Swell (Kenya), Africa's largest, and eight around Angola's Bié Swell (western Zambia, northern Namibia). A cluster of possible fans lies on the western margin of the Congo Basin (Mayombe Swell), and on the coastal slopes of the Namibia Swell. Sheer size may have militated aginst the recognition of many megafans: the largest in the Sahara are the Teghahart (378 km, Hoggar Swell, Algeria), and the Wadi Albalata (340 km, Uweinat Swell, Egypt). In southern Africa the largest are the Cubango (320 km, Bié Swell, Angola/ Namibia), and the Limpopo (230 km, Mozambique). (ii) Rift zones (a) Steer's horns basins-wide depressions centered on rifts. The largest contiguous group (n=14) developed in a steer's-horns basin occupies the wide Muglad depression (200-350 km, South Sudan). Four rift-related megafans lie SE of Lake Chad (Chad). Nine megafans occupy the complex Anza Rift in Kenya/South Somalia. The Salamat megafan (Chad), is unusual because it oriented parallel with the linked Salamat, Doseo and Doba rift axes, and is consequently one of the longest in Africa (465 km). (b) Rift depressions sensu stricto. Most rifts are too narrow to provide a transverse dimension large enough to accommodate megafans. Although well-known, the Okavango Rift (NW Botswana, NE Namibia) is unique in Africa in hosting three megafans within identifiable faulted margins. The Nile megafan is Africa's largest (476 km) and comprises the vast Sudd wetland (South Sudan). An explanation for its remarkable size may be its location in a depression at the junction of two conducive tectonic zones, the East African Swell margin and the Muglad steer's-horns depression. Discharge of the River Nile, the largest in the region, has allowed the Nile megafan to outcompete neighboring megafans for space.

Wilkinson, M. J.↗

Physiological Maturation of Regenerating Hair Cells

The bullfrog saccule, a sensor of gravity and substrate-borne vibration, is a model system for hair cell transduction. Saccular hair cells also increase in number throughout adult life and rapidly recover after hair cell damage, making this organ an ideal system for studying hair cell development, repair, and regeneration. We have used of hair cell and supporting cell immunocytochemical markers to identify damaged hair cells and hair cell precursors in organotypic cultures of the bullfrog saccule. We then used an innovative combination of confocal, electron, and time-lapse microscopy to study the fate of damaged hair cells and the origin of new hair cells after gentamicin ototoxicity in normal and mitotically blocked saccular cultures. These studies have shown that gentamicin ototoxicity produces both lethal and sublethal hair cell damage. They have also shown that hair cell recovery in this organ takes place by both the repair of sublethally damaged hair cells and by the replacement of lost hair cells by mitotic regeneration. In parallel studies, we have used biophysical and molecular biological techniques to study the differentiation and innervation of developing, repairing, and regenerating hair cells. More specifically, we have used RT-PCR to obtain the bullfrog homologues of L-type voltage- gated calcium (L-VGCC) and large-conductance Ca(2+)-activated potassium (BK) channel genes. We have then obtained probes for these genes and, using in situ hybridization, begun to examine their expression in the bullfrog saccule and amphibian papilla. We have also used fluorescent-labeled channel toxins and channel toxin derivatives to determine the time of appearance of L-type voltage-gated calcium (L-VGCC) and Ca(2+)-activated potassium (BK) channels and to study dynamic changes in the number, distribution, and co-localization of these proteins in developing, repairing, and regenerating hair cells. Using time-lapse microscopy, we are also studying the dynamic relationship between ion channel clustering and synaptic formation in hair cells and afferent neurons. In future studies, we will determine when hair cell precursors acquire electrical tuning, and, using whole-cell patch-clamp techniques, identify and characterize their L-VGCC and BK currents. We will also use biophysical techniques to determine the number of L-VGCC and BK channels and the size and gating kinetics of their underlying L-VGCC and BK conductances, correlating these variables with the amplitude and frequency of membrane oscillations produced by intracellular current steps. We expect these studies to determine how hair cells regulate ion channel expression to achieve specific physiological responses.

Baird, Richard A.↗

CFD Simulation of Liquid Rocket Engine Injectors

Detailed design issues associated with liquid rocket engine injectors and combustion chamber operation require CFD methodology which simulates highly three-dimensional, turbulent, vaporizing, and combusting flows. The primary utility of such simulations involves predicting multi-dimensional effects caused by specific injector configurations. SECA, Inc. and Engineering Sciences, Inc. have been developing appropriate computational methodology for NASA/MSFC for the past decade. CFD tools and computers have improved dramatically during this time period; however, the physical submodels used in these analyses must still remain relatively simple in order to produce useful results. Simulations of clustered coaxial and impinger injector elements for hydrogen and hydrocarbon fuels, which account for real fluid properties, is the immediate goal of this research. The spray combustion codes are based on the FDNS CFD code' and are structured to represent homogeneous and heterogeneous spray combustion. The homogeneous spray model treats the flow as a continuum of multi-phase, multicomponent fluids which move without thermal or velocity lags between the phases. Two heterogeneous models were developed: (1) a volume-of-fluid (VOF) model which represents the liquid core of coaxial or impinger jets and their atomization and vaporization, and (2) a Blob model which represents the injected streams as a cloud of droplets the size of the injector orifice which subsequently exhibit particle interaction, vaporization, and combustion. All of these spray models are computationally intensive, but this is unavoidable to accurately account for the complex physics and combustion which is to be predicted, Work is currently in progress to parallelize these codes to improve their computational efficiency. These spray combustion codes were used to simulate the three test cases which are the subject of the 2nd International Workshop on-Rocket Combustion Modeling. Such test cases are considered by these investigators to be very valuable for code validation because combustion kinetics, turbulence models and atomization models based on low pressure experiments of hydrogen air combustion do not adequately verify analytical or CFD submodels which are necessary to simulate rocket engine combustion. We wish to emphasize that the simulations which we prepared for this meeting are meant to test the accuracy of the approximations used in our general purpose spray combustion models, rather than represent a definitive analysis of each of the experiments which were conducted. Our goal is to accurately predict local temperatures and mixture ratios in rocket engines; hence predicting individual experiments is used only for code validation. To replace the conventional JANNAF standard axisymmetric finite-rate (TDK) computer code 2 for performance prediction with CFD cases, such codes must posses two features. Firstly, they must be as easy to use and of comparable run times for conventional performance predictions. Secondly, they must provide more detailed predictions of the flowfields near the injector face. Specifically, they must accurately predict the convective mixing of injected liquid propellants in terms of the injector element configurations.

Farmer, Richard↗