SEARCH · Engineering Papers

Results for “parallel codes”

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.

SourceResearch connectionFrom yearThrough year

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

At least 973 records · Page 54

All sources NASA NTRS DOE OSTI

NASA NTRS1998

Molecular genetics of root gravitropism and waving in Arabidopsis thaliana

When Arabidopsis thaliana seedlings grow embedded in an agar-based medium, their roots grow vertically downward. This reflects their ability to sense the gravity vector and to position their tip parallel to it (gravitropism). We have isolated a number of mutations affecting root gravitropism in Arabidopsis thaliana. One of these mutations, named arg1, affects root and hypocotyl gravitropism without promoting defects in starch content or in the ability of seedlings' organs to respond to plant hormones. The ARG1 gene was cloned and shown to code for a protein with a J domain at its amino terminus and a second sequence motif found in several cytoskeleton binding proteins. Mutations in the AGR1 locus promote a strong defect in root gravitropism. Some alleles also confer an increased root resistance to exogenous ethylene and an increased sensitivity to auxin. AGR1 was cloned and found to encode a putative transmembrane protein which might be involved in polar auxin transport, or in regulating the differential growth response to gravistimulation. When Arabidopsis seedlings grow on the surface of agar-based media tilted backward, their roots wave. That wavy pattern of root growth derives from a combined response to gravity, touch and other surface-derived stimuli. It is accompanied by a reversible rotation of the root tip about its axis. A number of mutations affect the presence or the shape of root waves on tilted agar-based surfaces. One of them, wvc1, promotes the formation of compressed root waves under these conditions. The physiological and molecular analyses of this mutant suggest that a tryptophan-derived molecule other than IAA might be an important regulator of the curvature responsible for root waving.

NASA Discipline Plant Biology↗

NASA NTRS2005

High Resolution Aerospace Applications using the NASA Columbia Supercomputer

This paper focuses on the parallel performance of two high-performance aerodynamic simulation packages on the newly installed NASA Columbia supercomputer. These packages include both a high-fidelity, unstructured, Reynolds-averaged Navier-Stokes solver, and a fully-automated inviscid flow package for cut-cell Cartesian grids. The complementary combination of these two simulation codes enables high-fidelity characterization of aerospace vehicle design performance over the entire flight envelope through extensive parametric analysis and detailed simulation of critical regions of the flight envelope. Both packages. are industrial-level codes designed for complex geometry and incorpor.ats. CuStomized multigrid solution algorithms. The performance of these codes on Columbia is examined using both MPI and OpenMP and using both the NUMAlink and InfiniBand interconnect fabrics. Numerical results demonstrate good scalability on up to 2016 CPUs using the NUMAIink4 interconnect, with measured computational rates in the vicinity of 3 TFLOP/s, while InfiniBand showed some performance degradation at high CPU counts, particularly with multigrid. Nonetheless, the results are encouraging enough to indicate that larger test cases using combined MPI/OpenMP communication should scale well on even more processors.

Mavriplis, Dimitri J.↗

NASA NTRS2013

Absorption Efficiencies of Forsterite. I: DDA Explorations in Grain Shape and Size

We compute the absorption efficiency (Q(sub abs)) of forsterite using the discrete dipole approximation (DDA) in order to identify and describe what characteristics of crystal grain shape and size are important to the shape, peak location, and relative strength of spectral features in the 8 - 40 micron wavelength range. Using the DDSCAT code, we compute Q(sub abs) for non-spherical polyhedral grain shapes with a(sub eff) = 0.1 micron. The shape characteristics identified are: 1) elongation/reduction along one of three crystallographic axes; 2) asymmetry, such that all three crystallographic axes are of different lengths; and 3) the presence of crystalline faces that are not parallel to a specific crystallographic axis, e.g., non-rectangular prisms and (di)pyramids. Elongation/reduction dominates the locations and shapes of spectral features near 10, 11, 16, 23.5, 27, and 33.5 micron, while asymmetry and tips are secondary shape effects. Increasing grain sizes (0.1 - 1.0 micron) shifts the 10, 11 micron features systematically towards longer wavelengths and relative to the 11 micron feature increases the strengths and slightly broadens the longer wavelength features. Seven spectral shape classes are established for crystallographic a-, b-, and c-axes and include columnar and platelet shapes plus non-elongated or equant grain shapes. The spectral shape classes and the effects of grain size have practical application in identifying or excluding columnar, platelet or equant forsterite grain shapes in astrophysical environs. Identification of the shape characteristics of forsterite from 8 - 40 micron spectra provides a potential means to probe the temperatures at which forsterite formed.

circumstellar matter↗

NASA NTRS1991

An experimental and computational study of 2-D blade-vortex interaction

An experimental and computational study is carried out to investigate the dominant physical factors of 2D parallel blade-vortex interaction (BVI) and its noise generation. A shock tube was used to generate a starting vortex which interacted with a target airfoil. Double-exposed holographic interferometry and airfoil surface pressure measurements were employed to obtain quantitative data during the BVI. As a numerical approach, thin-layer Navier-Stokes code, with a multizonal grid, was also used to resolve the phenomena occuring in the BVI, especially in the head-on collision case.

Lee, Soogab↗

NASA NTRS1996

Distributed Large Data-Object Environments: End-to-End Performance Analysis of High Speed Distributed Storage Systems in Wide Area ATM Networks

We have developed and deployed a distributed-parallel storage system (DPSS) in several high speed asynchronous transfer mode (ATM) wide area networks (WAN) testbeds to support several different types of data-intensive applications. Architecturally, the DPSS is a network striped disk array, but is fairly unique in that its implementation allows applications complete freedom to determine optimal data layout, replication and/or coding redundancy strategy, security policy, and dynamic reconfiguration. In conjunction with the DPSS, we have developed a 'top-to-bottom, end-to-end' performance monitoring and analysis methodology that has allowed us to characterize all aspects of the DPSS operating in high speed ATM networks. In particular, we have run a variety of performance monitoring experiments involving the DPSS in the MAGIC testbed, which is a large scale, high speed, ATM network and we describe our experience using the monitoring methodology to identify and correct problems that limit the performance of high speed distributed applications. Finally, the DPSS is part of an overall architecture for using high speed, WAN's for enabling the routine, location independent use of large data-objects. Since this is part of the motivation for a distributed storage system, we describe this architecture.

Johnston, William↗

NASA NTRS2003

Particle Acceleration in Relativistic Jets due to Weibel Instability

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The jitter radiation (Medvedev 2000) from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.↗

NASA NTRS2003

Particle Acceleration and Emission in Relativistic Jets

Shock wave acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. While some Fermi acceleration may occur at the jet front, the majority of electron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The 'jitter' radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.↗

NASA NTRS2003

Particle Acceleration and Radiation associated with Magnetic Field Generation from Relativistic Collisionless Shocks

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. While some Fermi acceleration may occur at the jet front, the majority of electron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.↗

NASA NTRS2004

Particle Acceleration in Relativistic Jets Due to Weibel Instability

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a three-dimensional relativistic electromagnetic particle code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. While some Fermi acceleration may occur at the jet front, the majority of electron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.-I.↗

NASA NTRS2010

Adaptive Numerical Algorithms in Space Weather Modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.

Toth, Gabor↗

NASA NTRS2022

Numerical Results for Polarized Light Scattering in a Spherical Atmosphere

We report numerical results for polarized light reflection from the top of a Rayleigh scattering spherical atmosphere with height-dependent single scattering albedo over a dark surface. Michael Mishchenko considered this scenario back in the 1990’s, for a plane-parallel atmosphere of unit optical thickness (OT), for which radiance errors arising from neglecting polarization reaches their highest values. To further extend Mishchenko’s results, we consider a value of OT = 0.25, for which the effect of atmospheric curvature is pronounced. New results are generated using three state-of-the art radiative transfer (RT) codes. These are: the MYSTIC and MCSSA models, which simulate light scattering in a true- spherical atmosphere using Monte Carlo methods; and the discrete ordinate code VLIDORT, operating with a new multiple-scatter spherical correction designed to deliver reasonable approximations to spherical-medium scattering. In this work, we report results for both single and multiple scattering; this will help to support the validation of existing and future polarized spherical RT codes, especially those using approximative methods to deal with sphericity.

S Korkin↗

NASA NTRS1999

Point Spread Function and Transmittance Analyses for Conical and Hexapod Secondary Mirror Support Towers for the Next Generation Space Telescope (NGST)

Different secondary mirror support towers were modeled on the CODE V optical design/analysis program for the NGST Optical Telescope Assembly (OTA) B. The vertices of the NGST OTA B primary and secondary mirrors were separated by close to 9.0 m. One type of tower consisted of a hollow cone 6.0 m long, 2.00 m in diameter at the base, and 0.704 m in diameter at its top. The base of the cone was considered attached to the primary's reaction structure through a hole in the primary. Extending up parallel to the optical axis from the top of this cylinder were eight blades (pyramidal struts) 3.0 m long. A cross section of each these long blades was an isosceles triangle with a base of 0.010 m and a height of 0.100 m with the sharpest part of each triangle pointing inward. The eight struts occurred every 45 deg. The other type of tower was purely a hexapod arrangement and had no blades or cones. The hexapod consisted simply of six, very thin, circular struts, leaving in pairs at 12:00, 4:00, and 8:00 at the primary and traversing to the outer edge of the back of the secondary mount. At this mount, two struts arrived at each of 10:00, 2:00, and 6:00. The struts were attached to the primary mirror in a ring 3.5 m in diameter. They reached the back of the secondary mount, a circle 0.704 m in diameter. Transmittance analyses at two levels were performed on the secondary mirror support towers. Detailed transmittances were accomplished by the use of the CODE V optical design/analysis program and were compared to transmittance calculations that were almost back-of-the-envelope. Point spread function (PSF) calculations, including both diffraction and aberration effects, were performed on CODE V. As one goes out from the center of the blur (for a point source), the two types of support towers showed little difference between their PSF intensities until one reaches about the 3% level. Contours can be delineated on CODE V down to about 10 (exp -8) times the peak intensity, fine detailing the differences in the wings of the PSF. Between the first type of tower and the hexapod tower, below 3% of the peak intensity, the PSF contour maps well delineated significant differences in the outer portions of the PSFS. Such differences in the wings of the PSF significantly impact the resolution of brown dwarfs or extra solar planets from their star. The hexapod tower provided the best PSF.

Wilkerson, Gary W.↗

NASA NTRS1989

Rarefied gas flow through two-dimensional nozzles

A kinetic theory analysis is made of the flow of a rarefied gas from one reservoir to another through two-dimensional nozzles with arbitrary curvature. The Boltzmann equation simplified by a model collision integral is solved by means of finite-difference approximations with the discrete ordinate method. The physical space is transformed by a general grid generation technique and the velocity space is transformed to a polar coordinate system. A numerical code is developed which can be applied to any two-dimensional passage of complicated geometry for the flow regimes from free-molecular to slip. Numerical values of flow quantities can be calculated for the entire physical space including both inside the nozzle and in the outside plume. Predictions are made for the case of parallel slots and compared with existing literature data. Also, results for the cases of convergent or divergent slots and two-dimensional nozzles with arbitrary curvature at arbitrary knudsen number are presented.

De Witt, Kenneth J.↗

NASA NTRS1992

Numerical analysis of hypersonic turbulent film cooling flows

As a building block, numerical capabilities for predicting heat flux and turbulent flowfields of hypersonic vehicles require extensive model validations. Computational procedures for calculating turbulent flows and heat fluxes for supersonic film cooling with parallel slot injections are described in this study. Two injectant mass flow rates with matched and unmatched pressure conditions using the database of Holden et al. (1990) are considered. To avoid uncertainties associated with the boundary conditions in testing turbulence models, detailed three-dimensional flowfields of the injection nozzle were calculated. Two computational fluid dynamics codes, GASP and FDNS, with the algebraic Baldwin-Lomax and k-epsilon models with compressibility corrections were used. It was found that the B-L model which resolves near-wall viscous sublayer is very sensitive to the inlet boundary conditions at the nozzle exit face. The k-epsilon models with improved wall functions are less sensitive to the inlet boundary conditions. The testings show that compressibility corrections are necessary for the k-epsilon model to realistically predict the heat fluxes of the hypersonic film cooling problems.

Chen, Y. S.↗

NASA NTRS1999

Portable Parallel Programming for the Dynamic Load Balancing of Unstructured Grid Applications

The ability to dynamically adapt an unstructured -rid (or mesh) is a powerful tool for solving computational problems with evolving physical features; however, an efficient parallel implementation is rather difficult, particularly from the view point of portability on various multiprocessor platforms We address this problem by developing PLUM, tin automatic anti architecture-independent framework for adaptive numerical computations in a message-passing environment. Portability is demonstrated by comparing performance on an SP2, an Origin2000, and a T3E, without any code modifications. We also present a general-purpose load balancer that utilizes symmetric broadcast networks (SBN) as the underlying communication pattern, with a goal to providing a global view of system loads across processors. Experiments on, an SP2 and an Origin2000 demonstrate the portability of our approach which achieves superb load balance at the cost of minimal extra overhead.

Biswas, Rupak↗

NASA NTRS1989

Run-time scheduling and execution of loops on message passing machines

Sparse system solvers and general purpose codes for solving partial differential equations are examples of the many types of problems whose irregularity can result in poor performance on distributed memory machines. Often, the data structures used in these problems are very flexible. Crucial details concerning loop dependences are encoded in these structures rather than being explicitly represented in the program. Good methods for parallelizing and partitioning these types of problems require assignment of computations in rather arbitrary ways. Naive implementations of programs on distributed memory machines requiring general loop partitions can be extremely inefficient. Instead, the scheduling mechanism needs to capture the data reference patterns of the loops in order to partition the problem. First, the indices assigned to each processor must be locally numbered. Next, it is necessary to precompute what information is needed by each processor at various points in the computation. The precomputed information is then used to generate an execution template designed to carry out the computation, communication, and partitioning of data, in an optimized manner. The design is presented for a general preprocessor and schedule executer, the structures of which do not vary, even though the details of the computation and of the type of information are problem dependent.

Crowley, Kay↗

NASA NTRS1990

Run-time scheduling and execution of loops on message passing machines

Saltz, Joel↗

NASA NTRS1996

A New Approach to Parallel Interference Cancellation for CDMA

This paper introduces an improved nonlinear parallel interference cancellation scheme that significantly reduces the degrading effect of user interference with implementation complexity linear in the number of users. The scheme operates on the fact that parallel processing simultaneously removes from each user a part of the interference produced by the remaining users accessing the channel the amount being proportional to their reliability. The parallel processing can be done in multiple stages. Simulation results are given for a multitude of different situations, in particular those cases for which the analysis is too complex.

interference cancellation multiuser communication ↗