SEARCH · Engineering Papers
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.
A parallel algorithm for 3D FFTs is implemented as a series of local 1D FFTs combined with data transposes. This allows the use of vendor supplied (often fully optimized) sequential 1D FFTs. The FFTs are carried out in-place by using an in-place data transpose across the processors.
We discuss attempts to explain the steepening of magnetohydrodynamic waves upstream of the Earth's bow shock and the formation of short-wavelength 'shocklets' in terms of simple analytic models. The effort involves use of the Derivative Nonlinear Schroedinger Equation (DNLS) to describe the wave evolution. We review shortcomings of previous attempts to model foreshock wave phenomena in terms of wave packet evolution according to the DNLS. It is pointed out that: (1) the oblique propagation of the waves and (2) their growth by unstable particle distribution are conducive to the processes of wave steepening and shocklet formation. Numerical solutions of the DNLS with initial conditions corresponding to envelope-modulated, obliquely propagating wave packets yield evolved wave packets with many properties similar to the observed waves.
In this talk, we will discuss how plasma distribution functions from the Magnetospheric Multiscale (MMS) mission transformed the picture on how the solar wind is shocked. For quasi-parallel shocks with typical Alfven Mach numbers, we will examine how ion deceleration, reflection, and thermalization take place, and illustrate how the new picture expands on the framework envisioned by Schwartz and Burgess in 1991. For quasi-perpendicular shocks with the range of Alfven Mach numbers (~1-60) captured by MMS, we will examine examples showing the processes of ion and electron thermalization go beyond the existing paradigm and call for future investigation. The knowledge advanced by MMS measurements provides a ground to understand planetary shocks in our solar system as well as astrophysical shocks arising from stellar wind interaction with exoplanets and supernova remnants.
In 1993, a proposal at the Japan.-U.S. Cooperation in Space Program Workshop lead to a subsequent series of satellite communications experiments and demonstrations, under the title of Trans-Pacific High Data Rate Satellite Communications Experiments. The first of which is a joint collaboration between government and industry teams in the United States and Japan that successfully demonstrated distributed high definition video (HDV) post-production on a global scale using a combination of high data rate satellites and terrestrial fiber optic asynchronous transfer mode (ATM) networks. The HDV experiment is the first GIBN experiment to establish a dual-hop broadband satellite link for the transmission of digital HDV over ATM. This paper describes the team's effort in using the NASA Advanced Communications Technology Satellite (ACTS) at rates up to OC-3 (155 Mbps) between Los Angeles and Honolulu, and using Intelsat at rates up to DS-3 (45 Mbps) between Kapolei and Tokyo, with which HDV source material was transmitted between Sony Pictures High Definition Center (SPHDC) in Los Angeles and Sony Visual Communication Center (VCC) in Shinagawa, Tokyo. The global-scale connection also used terrestrial networks in Japan, the States of Hawaii and California. The 1.2 Gbps digital HDV stream was compressed down to 22.5 Mbps using a proprietary Mitsubishi MPEG-2 codec that was ATM AAL-5 compatible. The codec: employed four-way parallel processing. Improved versions of the codec are now commercially available. The successful post-production activity performed in Tokyo with a HDV clip transmitted from Los Angeles was predicated on the seamless interoperation of all the equipment between the sites, and was an exciting example in deploying a global-scale information infrastructure involving a combination of broadband satellites and terrestrial fiber optic networks. Correlation of atmospheric effects with cell loss, codec drop-out, and picture quality were made. Current efforts in the Trans-Pacific series plan to examine the use of Internet Protocol (IP)-related technologies over such an infrastructure. The use of IP allows the general public to be an integral part of the exciting activities, helps to examine issues in constructing the solar-system internet, and affords an opportunity to tap the research results from the (reliable) multicast and distributed systems communities. The current Trans- Pacific projects, including remote astronomy and digital library (visible human) are briefly described.
Growing demand for visual, user-friendly representation of information inspires search for the new methods of image transmission. Currently used in-series (sequential) methods of information processing are inherently slow and are designed mainly for transmission of one or two dimensional arrays of data. Conventional transmission of data by fibers requires many fibers with array of laser diodes and photodetectors. In practice, fiber bundles are also used for transmission of images. Image is formed on the fiber-optic bundle entrance surface and each fiber transmits the incident image to the exit surface. Since the fibers do not preserve phase, only 2D intensity distribution can be transmitted in this way. Each single mode fiber transmit only one pixel of an image. Multimode fibers may be also used, so that each mode represent different pixel element. Direct transmission of image through multimode fiber is hindered by the mode scrambling and phase randomization. To overcome these obstacles wavelength and time-division multiplexing have been used, with each pixel transmitted on a separate wavelength or time interval. Phase-conjugate techniques also was tested in, but only in the unpractical scheme when reconstructed image return back to the fiber input end. Another method of three-dimensional imaging over single mode fibers was demonstrated in, using laser light of reduced spatial coherence. Coherence encoding, needed for a transmission of images by this methods, was realized with grating interferometer or with the help of an acousto-optic deflector. We suggest simple practical holographic method of image transmission over single multimode fiber or over fiber bundle with coherent light using filtering by holographic optical elements. Originally this method was successfully tested for the single multimode fiber. In this research we have modified holographic method for transmission of laser illuminated images over commercially available fiber bundle (fiber endoscope, or fiberscope).
Kusunose, Mineshige & Yamada (1996; hereafter KMY) extended the model of Kusunose & Mineshige (1995) to the Galactic black hole candidates by considering nonthermal electron injection with gamma(EQ\0(,\s\up2(less than),\s\do-l(_))) 10. The effects of pair escape and advection on the disk structure and general relativistic effects on the emission spectrum were also examined. They found that the energy spectral index (alpha)(sub x) of the power law X-rays is about-0.8 and-2.0 when 1(sub soft)/1 = 0.2 and 2, respectively, where 1(sub soft)/1 is the ratio of the compactness of the injected soft photons to that of the gravitational energy. The power law index was found to be nearly independent of the mass accretion which is consistent with the observed luminosity independence. The model with small 1(sub soft)/1 (less than 1) shows promise for explaining the low state observed in Galactic black hole candidates. Model fits were provided for GX339-4 and Cyg X- 1 data from COMPTEL and OSSE on the Compton Gamma Ray Observatory. The difference in emission spectra between thermal disks and the model of KMY appears only in the energy range greater than 100 keV. Li, Kusunose and Liang (1996) studied stochastic particle acceleration to produce nonthermal particle distributions which then were used in the model of Kusunose & Mineshige (1995) to model the spectrum above 1 Mev from GBHC's. Under certain conditions, stochastic electron acceleration overcame Coulomb and Compton losses resulting in a suprathermal electron population. Good fits were obtained by COMPTEL and OSSE observations of Cyg X-1 and GRO J0422+23. Kusunose & Mineshige (1996a) examined the role of electron-positron pairs in advection-dominated disks. They found that the results for advection-dominated disks without pairs are not qualitatively changed by including pairs. Summaries of work sponsored by this grant are given in Wheeler, Kim, Moscoso, Kusunose & Mineshige (1996) and Kusunose (1996) Work was also done on developing a model for an e(+-) pair wind from the inner disk region of a black hole. The model consists of three zones: a pair production/annihilation zone at the base of the wind, a pair annihilation zone slightly further out from the disk and a wind zone extending to infinity where no annihilation occurs. The model assumes an input X-ray / gamma-ray power-law spectrum as a function of photon energy and angular distribution. Pairs are created in the pair production/annihilation zone via photon-photon, photonparticle and particle-particle collisions. The bulk velocity of the pairs is obtained from the momentum component perpendicular to the disk taking into account the radiation pressure acceleration of the pairs. Energy balance in the pair production/annihilation zone is calculated by assuming that the momentum component parallel to the disk is thermalized and by taking into account heating/cooling via anisotropic thermal Compton scattering along with the cooling processes of bremsstrahlung and synchrotron radiation. The pair density is also calculated. The outer boundary of the pair production/annihilation zone is where the optical depth for photon-photon pair produciton is unity. Beyond this point only pair annihilation will occur in the pair annihilation zone. When the pair density becomes small, pair annihilation will be negligible and the pairs will flow freely to infinity. In this model we have found that the X-ray power-law spectral index, alpha(sub x) is the primary parameter which determines the density, temperature and velocity of the pair production/annihilation zone (Moscoso, Kusunose & Wheeler 1996).
Observations from several experiments on board the Dynamics Explorer 1 and 2 (DE 1 and 2) spacecraft and ground-based radar measurements from the Chatanika radar are combined in order to examine the details of ionospheric/magnetospheric coupling in the local evening sector. DE 1 and DE 2 were in coplanar polar orbits that provided measurements almost simultaneously in time and magnetically coincident with the Chatanika radar from L = 3 to L = 17. The coupling processes are inferred from the density, temperature, composition, and angular distributions of the low-energy plasma observed from the E region of the ionosphere to magnetospheric altitudes of 2.5 earth radii. Plasma characteristics of the plasmasphere, main trough, auroral zone, and polar cap can be studied in this data set. The observations imply that as L increases, the dominant coupling mechanism between the ionosphere and magnetosphere in the measured energy range changes from equilibrium diffusion to perpendicular acceleration and finally to parallel acceleration.
The initiation of the breakdown process for axially symmetric compressible flows is investigated using a numerical solution of the conservation equations for mass, momentum, and energy. The vortex is isolated, with its axis parallel to the direction of the main stream, and the core radius is small compared to the breakdown length. Computations for several flowfields indicate that the breakdown of the solution is shifted further downstream with increasing Mach number until breakdown is no longer observed. In the subsonic case, the influence of the initial temperature distribution on the breakdown length of the solution is more pronounced than in the supersonic case, with heating of the core enhancing breakdown, and cooling delaying it. The breakdown of the solution is seen to always occur for nonvanishing axial velocity components.
A portable computer code called BLASIM is developed at NASA LeRC to assess the ice impact damage on aircraft engine blades. In addition to the ice impact analyses, the code is also capable of carrying out static, dynamic, resonance margin and flutter analyses. The blade can be solid, hollow, superhybrid or composite material. An optional preprocessor (input generator) is also developed to generate input to the code through interactive process. The blade geometry can be defined either by a series of airfoils at discrete input stations or by a finite element grid. The code employs a coarse fixed finite element mesh with triangular plate finite elements and has quick turnaround time. The ice piece is modeled as an equivalent spherical object and has the velocity opposite to that of the aircraft with direction parallel to the engine axis. For the local impact damage assessment, the impact force is considered as a distributed load acting over a region around the impact point and the average radial strain of the finite elements along the leading edge is taken as a measure of the local damage. To estimate the damage at the blade root, the impact is considered to be an impulse and a combined stress failure criteria is employed. Parametric studies for local and root ice impact damage, and post-impact dynamics are discussed for solid and composite blades.
This thesis studies a forward recovery strategy using checkpointing and optimistic execution in parallel and distributed systems. The approach uses replicated tasks executing on different processors for forwared recovery and checkpoint comparison for error detection. To reduce overall redundancy, this approach employs a lower static redundancy in the common error-free situation to detect error than the standard N Module Redundancy scheme (NMR) does to mask off errors. For the rare occurrence of an error, this approach uses some extra redundancy for recovery. To reduce the run-time recovery overhead, look-ahead processes are used to advance computation speculatively and a rollback process is used to produce a diagnosis for correct look-ahead processes without rollback of the whole system. Both analytical and experimental evaluation have shown that this strategy can provide a nearly error-free execution time even under faults with a lower average redundancy than NMR.
This paper presents our experiences of parallelizing the sequential implementation of NAS benchmarks using compiler directives on SGI Origin2000 distributed shared memory (DSM) system. Porting existing applications to new high performance parallel and distributed computing platforms is a challenging task. Ideally, a user develops a sequential version of the application, leaving the task of porting to new generations of high performance computing systems to parallelization tools and compilers. Due to the simplicity of programming shared-memory multiprocessors, compiler developers have provided various facilities to allow the users to exploit parallelism. Native compilers on SGI Origin2000 support multiprocessing directives to allow users to exploit loop-level parallelism in their programs. Additionally, supporting tools can accomplish this process automatically and present the results of parallelization to the users. We experimented with these compiler directives and supporting tools by parallelizing sequential implementation of NAS benchmarks. Results reported in this paper indicate that with minimal effort, the performance gain is comparable with the hand-parallelized, carefully optimized, message-passing implementations of the same benchmarks.
Multidisciplinary design, analysis and optimization involves modeling the interactions of complex systems across a variety of disciplines. The optimization of such systems can be a computationally expensive exercise with multiple levels of nested nonlinear solvers running under an optimizer.The application of optimal control in project development often involves performing trajectory optimization for fixed vehicle designs or parametric sweeps across some key vehicle properties.This information is then relayed to the subsystem design teams who update their designs and relay some bulk characteristics back to the trajectory optimization procedure.This iteration is then repeated until the design closes.However, with increasing interest in more tightly coupled systems, such as electric and hybrid-electric aircraft propulsion and boundary layer ingestion, this process is prone to ignore subtle coupling between vehicle subsystem designs and vehicle operation on a given mission.Integrating trajectory optimization into a tightly coupled multidisciplinary design procedure can be computationally prohibitive, depending on the complexity of the subsystem analyses and the optimal control technique applied.To address these issues a new optimal control software tool, Dymos, has been developed.Dymos is built upon NASA's OpenMDAO software and can leverage its capabilities to efficiently compute gradients for the optimization and optimize complex models in parallel on distributed memory systems.This report provides some explanation into the numerical methods employed in Dymos and provides several use cases that demonstrate its performance on traditional optimal control problems and improvements ino techniques have been used extensively in recent decades to solve a variety of optimal control problems, typically in the form of aerospace vehicle trajectory optimization.
Johnson Space Center's Mission Control Center is a space vehicle, space program agnostic facility. The current operational design is essentially identical to the original facility architecture that was developed and deployed in the mid-90's. In an effort to streamline the support costs of the mission critical facility, the Mission Operations Division (MOD) of Johnson Space Center (JSC) has sponsored an exploratory project to evaluate and inject current state-of-the-practice Information Technology (IT) tools, processes and technology into legacy operations. The general push in the IT industry has been trending towards a data-centric computer infrastructure for the past several years. Organizations facing challenges with facility operations costs are turning to creative solutions combining hardware consolidation, virtualization and remote access to meet and exceed performance, security, and availability requirements. The Operations Technology Facility (OTF) organization at the Johnson Space Center has been chartered to build and evaluate a parallel Mission Control infrastructure, replacing the existing, thick-client distributed computing model and network architecture with a data center model utilizing virtualization to provide the MCC Infrastructure as a Service. The OTF will design a replacement architecture for the Mission Control Facility, leveraging hardware consolidation through the use of blade servers, increasing utilization rates for compute platforms through virtualization while expanding connectivity options through the deployment of secure remote access. The architecture demonstrates the maturity of the technologies generally available in industry today and the ability to successfully abstract the tightly coupled relationship between thick-client software and legacy hardware into a hardware agnostic "Infrastructure as a Service" capability that can scale to meet future requirements of new space programs and spacecraft. This paper discusses the benefits and difficulties that a migration to cloud-based computing philosophies has uncovered when compared to the legacy Mission Control Center architecture. The team consists of system and software engineers with extensive experience with the MCC infrastructure and software currently used to support the International Space Station (ISS) and Space Shuttle program (SSP).
We describe how high spatial resolution imaging of circumstellar dust at a wavelength of about 10 micron, combined with knowledge of the source spectral energy distribution, can yield useful information about the sizes of the individual dust grains responsible for the infrared emission. Much can be learned even when only upper limits to source size are available. In parallel with high-resolution single-telescope imaging that may resolve the more extended mid-infrared sources, we plan to apply these less direct techniques to interpretation of future observations from two-element optical interferometers, where quite general arguments may be made despite only crude imaging capability. Results to date indicate a tendency for circumstellar grain sizes to be rather large compared to the Mathis-Rumpl-Nordsieck size distribution traditionally thought to characterize dust in the general interstellar medium. This may mean that processing of grains after their initial formation and ejection from circumstellar atmospheres adjusts their size distribution to the ISM curve; further mid-infrared observations of grains in various environments would help to confirm this conjecture.
The Medical Data Architecture (MDA) project supports the Exploration Medical Capability (ExMC) risk to minimize or reduce the risk of adverse health outcomes and decrements in performance due to in-flight medical capabilities on human exploration missions. To mitigate this risk, the ExMC MDA project addresses the technical limitations identified in ExMC Gap Med 07: We do not have the capability to comprehensively process medically-relevant information to support medical operations during exploration missions. This gap identifies that the current in-flight medical data management includes a combination of data collection and distribution methods that are minimally integrated with on-board medical devices and systems. Furthermore, there are a variety of data sources and methods of data collection. For an exploration mission, the seamless management of such data will enable a more medically autonomous crew than the current paradigm. The medical system requirements are being developed in parallel with the exploration mission architecture and vehicle design. ExMC has recognized that in order to make informed decisions about a medical data architecture framework, current methods for medical data management must not only be understood, but an architecture must also be identified that provides the crew with actionable insight to medical conditions. This medical data architecture will provide the necessary functionality to address the challenges of executing a self-contained medical system that approaches crew health care delivery without assistance from ground support. Hence, the products supported by current prototype development will directly inform exploration medical system requirements.
The Medical Data Architecture (MDA) project supports the Exploration Medical Capability (ExMC) risk to minimize or reduce the risk of adverse health outcomes and decrements in performance due to in-flight medical capabilities on human exploration missions. To mitigate this risk, the ExMC MDA project addresses the technical limitations identified in ExMC Gap Med 07: We do not have the capability to comprehensively process medically-relevant information to support medical operations during exploration missions. This gap identifies that the current in-flight medical data management includes a combination of data collection and distribution methods that are minimally integrated with on-board medical devices and systems. Furthermore, there are a variety of data sources and methods of data collection. For an exploration mission, the seamless management of such data will enable a more medically autonomous crew than the current paradigm. The medical system requirements are being developed in parallel with the exploration mission architecture and vehicle design. ExMC has recognized that in order to make informed decisions about a medical data architecture framework, current methods for medical data management must not only be understood, but an architecture must also be identified that provides the crew with actionable insight to medical conditions. This medical data architecture will provide the necessary functionality to address the challenges of executing a self-contained medical system that approaches crew health care delivery without assistance from ground support. Hence, the products supported by current prototype development will directly inform exploration medical system requirements.
Recent developments in the field of nonequilibrium thermodynamics associated with viscous flows are examined and related to developments to the understanding of specific phenomena in aerodynamics and aeroacoustics. A key element of the nonequilibrium theory is the principle of minimum entropy production rate for steady dissipative processes near equilibrium, and variational calculus is used to apply this principle to several examples of viscous flow. A review of nonequilibrium thermodynamics and its role in fluid motion are presented. Several formulations are presented of the local entropy production rate and the local energy dissipation rate, two quantities that are of central importance to the theory. These expressions and the principle of minimum entropy production rate for steady viscous flows are used to identify parallel-wall channel flow and irrotational flow as having minimally dissipative velocity distributions. Features of irrotational, steady, viscous flow near an airfoil, such as the effect of trailing-edge radius on circulation, are also found to be compatible with the minimum principle. Finally, the minimum principle is used to interpret the stability of infinitesimal and finite amplitude disturbances in an initially laminar, parallel shear flow, with results that are consistent with experiment and linearized hydrodynamic stability theory. These results suggest that a thermodynamic approach may be useful in unifying the understanding of many diverse phenomena in aerodynamics and aeroacoustics.
The National Aeronautics and Space Administration (NASA) Data Assimilation Office (DAO) at the Goddard Space Flight Center (GSFC) has developed the GEOS DAS, a data assimilation system that provides production support for NASA missions and will support NASA's Earth Observing System (EOS) in the coming years. The GEOS DAS will be used to provide background fields of meteorological quantities to EOS satellite instrument teams for use in their data algorithms as well as providing assimilated data sets for climate studies on decadal time scales. The DAO has been involved in prototyping parallel implementations of the GEOS DAS for a number of years and is now embarking on an effort to convert the production version from shared-memory parallelism to distributed-memory parallelism using the portable Message-Passing Interface (MPI). The GEOS DAS consists of two main components, an atmospheric General Circulation Model (GCM) and a Physical-space Statistical Analysis System (PSAS). The GCM operates on data that are stored on a regular grid while PSAS works with observational data that are scattered irregularly throughout the atmosphere. As a result, the two components have different data decompositions. The GCM is decomposed horizontally as a checkerboard with all vertical levels of each box existing on the same processing element(PE). The dynamical core of the GCM can also operate on a rotated grid, which requires communication-intensive grid transformations during GCM integration. PSAS groups observations on PEs in a more irregular and dynamic fashion.