Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical performance”

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 127 records · Page 7

Extended hot-gas halos around starburst galaxies

A reanalysis of Einstein Imaging Proportional Counter (IPC) data and new observations from the Ginga Large Area (Proportional) counters (LAC) indicate the presence of extended X-ray emission (10-50 kpc) around the starburst galaxy M82. Here, we discuss our model of this emission, which was obtained by performing numerical hydrodynamic simulations of the starburst event to much later times and larger scales than were previously considered. For our models, we adopted a supernova rate of 0.1/yr, and an extended low-density static halo that is bound to the galaxy. There are three stages to the evolution of the wind-blown bubble and the propagation of the shock front: the bubble expands in an almost uniform density disk gas, with a deceleration of the shock front (t is less than or approximately 3.6 Myr); breakout from the disk and the upward acceleration of the shock front (3.6 Myr is less than or approximately t is less than or approximately 18 Myr); propagation into the halo, leading to a more spherical system and shock deceleration (18 Myr is less than or approximately t). For a halo density of 10(exp -3)/cu cm, the outflow reaches a distance of 40-50 kpc from the center of the starburst galaxy in 50 Myr. We calculated the time evolution of the X-ray luminosity and found that the extended starburst emits 3 x 10(exp 39) erg/s to 10(exp 40) in the Ginga LAC band and approximately 10(exp 41) erg/s in the Einstein band. The degree of the ionization equilibrium in the outflow and its effect on the iron K alpha line emission are discussed.

Tomisaka, Kohji↗

Parallel-vector out-of-core equation solver for computational mechanics

A parallel/vector out-of-core equation solver is developed for shared-memory computers, such as the Cray Y-MP machine. The input/ output (I/O) time is reduced by using the a synchronous BUFFER IN and BUFFER OUT, which can be executed simultaneously with the CPU instructions. The parallel and vector capability provided by the supercomputers is also exploited to enhance the performance. Numerical applications in large-scale structural analysis are given to demonstrate the efficiency of the present out-of-core solver.

Qin, J.↗

Fourier/Chebyshev methods for the incompressible Navier-Stokes equations in finite domains

A fully spectral numerical scheme for the incompressible Navier-Stokes equations in domains which are infinite or semi-infinite in one dimension. The domain is not mapped, and standard Fourier or Chebyshev expansions can be used. The handling of the infinite domain does not introduce any significant overhead. The scheme assumes that the vorticity in the flow is essentially concentrated in a finite region, which is represented numerically by standard spectral collocation methods. To accomodate the slow exponential decay of the velocities at infinity, extra expansion functions are introduced, which are handled analytically. A detailed error analysis is presented, and two applications to Direct Numerical Simulation of turbulent flows are discussed in relation with the numerical performance of the scheme.

Corral, Roque↗

Steady-State and Transient Boundary Element Methods for Coupled Heat Conduction

Boundary element algorithms for the solution of steady-state and transient heat conduction are presented. The algorithms are designed for efficient coupling with computational fluid dynamic discretizations and feature piecewise linear elements with offset nodal points. The steady-state algorithm employs the fundamental solution approach; the integration kernels are computed analytically based on linear shape functions, linear elements, and variably offset nodal points. The analytic expressions for both singular and nonsingular integrands are presented. The transient algorithm employs the transient fundamental solution; the temporal integration is performed analytically and the nonsingular spatial integration is performed numerically using Gaussian quadrature. A series solution to the integration is derived for the instance of a singular integrand. The boundary-only character of the algorithm is maintained by integrating the influence coefficients from initial time. Numerical results are compared to analytical solutions to verify the current boundary element algorithms. The steady-state and transient algorithms are numerically shown to be second-order accurate in space and time, respectively.

Kontinos, Dean A.↗

Sensitivity of Regional Climate to Deforestation in the Amazon Basin

The deforestation results in several adverse effect on the natural environment. The focus of this paper is on the effects of deforestation on land-surface processes and regional climate of the Amazon basin. In general, the effect of deforestation on climate are likely to depend on the scale of the defrosted area. In this study, we are interested in the effects due to deforestation of areas with a scale of about 250 km. Hence, a meso-scale climate model is used in performing numerical experiments on the sensitivity of regional climate to deforestation of areas with that size. It is found that deforestation results in less net surface radiation, less evaporation, less rainfall, and warmer surface temperature. The magnitude of the of the change in temperature is of the order 0.5 C, the magnitudes of the changes in the other variables are of the order of IO%. In order to verify some of he results of the numerical experiments, the model simulations of net surface radiation are compared to recent observations of net radiation over cleared and undisturbed forest in the Amazon. The results of the model and the observations agree in the following conclusion: the difference in net surface radiation between cleared and undisturbed forest is, almost, equally partioned between net solar radiation and net long-wave radiation. This finding contributes to our understanding of the basic physics in the deforestation problem.

Eltahir, Elfatih A. B.↗

New Parallel Algorithms for Structural Analysis and Design of Aerospace Structures

Subspace and Lanczos iterations have been developed, well documented, and widely accepted as efficient methods for obtaining p-lowest eigen-pair solutions of large-scale, practical engineering problems. The focus of this paper is to incorporate recent developments in vectorized sparse technologies in conjunction with Subspace and Lanczos iterative algorithms for computational enhancements. Numerical performance, in terms of accuracy and efficiency of the proposed sparse strategies for Subspace and Lanczos algorithm, is demonstrated by solving for the lowest frequencies and mode shapes of structural problems on the IBM-R6000/590 and SunSparc 20 workstations.

Nguyen, Duc T.↗

Global Artificial Boundary Conditions for Computation of External Flow Problems with Propulsive Jets

We propose new global artificial boundary conditions (ABC's) for computation of flows with propulsive jets. The algorithm is based on application of the difference potentials method (DPM). Previously, similar boundary conditions have been implemented for calculation of external compressible viscous flows around finite bodies. The proposed modification substantially extends the applicability range of the DPM-based algorithm. In the paper, we present the general formulation of the problem, describe our numerical methodology, and discuss the corresponding computational results. The particular configuration that we analyze is a slender three-dimensional body with boat-tail geometry and supersonic jet exhaust in a subsonic external flow under zero angle of attack. Similarly to the results obtained earlier for the flows around airfoils and wings, current results for the jet flow case corroborate the superiority of the DPM-based ABC's over standard local methodologies from the standpoints of accuracy, overall numerical performance, and robustness.

Tsynkov, Semyon↗

Interpreting Broad Double-Peaked Emission Lines in Active Galactic Nuclei

The principal objectives of this project were to probe the inner regions of active galactic nuclei and to test general relativity in the strong-field limit. The approach takes advantage of broad atomic line emission observed from material deep in the potential well of an active galactic nucleus which contains key information as to the physics of the system. Line profiles in a wide range of wavebands from optical to X-ray have provided compelling evidence of the existence of a relativistic accretion disk around a supermassive black hole in a number of galaxies. The simplest model posits a geometrically thin disk in Keplerian orbit, with general relativistic effects in evidence. This model is the point of departure for the proposed work. We developed a high-performance numerical code to calculate photon trajectories in a Schwarzschild or Kerr metric and implemented it on parallel supercomputers. This code includes a general purpose ray tracer that calculates line profiles, light curves, and other observable quantities for a wide variety of emitter configurations. The versatility comes from the fact that the ray tracing algorithm does not depend on any symmetries regarding emitter locations. The speed comes from parallel implementation which enables us to sample hitherto unattainable volumes of disk model parameter space. During the period 1 March 1997 through 28 February 1998, two papers, supported in whole or in part by this grant, were published in refereed journals. They are reproduced in their entirety in the next two sections of this report.

Halpern, Jules↗

Nonlinear Insolation Forcing: A Physical Mechanism for Climate Change

This paper focuses on recent advances in the understanding of nonlinear insolation forcing for climate change. The amplitude-frequency resonances in the insolation variations induced by the Earth's changing obliquity are emergent and may provide a physical mechanism to drive the glaciation cycles. To establish the criterion that nonlinear insolation forcing is responsible for major climate changes, the cooperative phenomena between the frequency and amplitude of the insolation are defined as insolation pulsation. Coupling of the insolation frequency and amplitude variations has established an especially new and interesting series of insolation pulses. These pulses would modulate the insolation in such a way that the mode of insolation variations could be locked to generate the 100-kyr ice age cycle which is a long-time geophysical puzzle. The nonlinear behavior of insolation forcing is tested by energy balance and ice sheet climate models and the physical mechanism behind this forcing is explained in terms of pulse duration in the incoming solar radiation. Calculations of the solar energy flux at the top of the atmosphere show that the duration of the negative and positive insolation pulses is about 2 thousand years which is long enough to prolong glaciation into deep ice ages and cause rapid melting of large ice sheets in the high latitudes of the northern hemisphere. We have performed numerical simulations of climate response to nonlinear insolation forcing for the past 2 million years. Our calculated results of temperature fluctuations are in good agreement with the climate cycles as seen in the terrestrial biogenic silica (BDP-96-2) data as well as in the marine oxygen isotope (delta(sup 18)O) records.

Liu, H. S.↗

ER-2 Observations of Precipitation Systems During TRMM-LBA

The NASA ER-2 performed numerous flights over precipitation systems in Rondonia, Brazil. The ER-2 carried a payload including the ER-2 Doppler Radar (EDOP), the Advanced Microwave Precipitation Radiometer (AMPR), the Lightning Instrument Package, and other instruments. This presentation will overview the types of data sets collected during TRMM-LBA (Tropical Rainfall Measuring Mission Satellite-Large Scale Biosphere Atmosphere Experiment in Amazonia) with particular emphasis on EDOP measurements. Numerous cases of convection ranging from weak to very intense, were overflown by the ER-2. Two TRMM overpasses were coincident with ER-2 flights which allowed for intercomparisons between the Precipitation Radiometer (PR), EDOP, and the S-POL (S-band Polarimetric Radar) and TOGA (Tropical Oceans and Global Atmosphere) ground-based radars. Preliminary results from this comparison will be presented as well as initial selection of case studies and efforts involving vertical motions in convection.

Heymsfield, Gerald↗

Theory and Modeling in Support of Tether

This final report summarizes the work performed by SAIC's Applied Physics Operation on the modeling and support of Tethered Satellite System missions (TSS-1 and TSS-1R). The SAIC team, known to be Theory and Modeling in Support of Tether (TMST) investigation, was one of the original twelve teams selected in July, 1985 for the first TSS mission. The accomplishments described in this report cover the period December 19, 1985 to September 31, 1999 and are the result of a continuous effort aimed at supporting the TSS missions in the following major areas. During the contract period, the SAIC's TMST investigation acted to: Participate in the planning and the execution on both of the TSS missions; Provide scientific understanding on the issues involved in the electrodynamic tether system operation prior to the TSS missions; Predict ionospheric conditions encountered during the re-flight mission (TSS-lR) based on realtime global ionosounde data; Perform post mission analyses to enhance our understanding on the TSS results. Specifically, we have 1) constructed and improved current collection models and enhanced our understanding on the current-voltage data; 2) investigated the effects of neutral gas in the current collection processes; 3) conducted laboratory experiments to study the discharge phenomena during and after tether-break; and 4) perform numerical simulations to understand data collected by plasma instruments SPES onboard the TSS satellite; Design and produce multi-media CD that highlights TSS mission achievements and convey the knowledge of the tether technology to the general public. Along with discussions of this work, a list of publications and presentations derived from the TMST investigation spanning the reporting period is compiled.

Chang, C. L.↗

Flight Dynamics Analysis Branch End of Fiscal Year 1999 Report

This document summarizes the major activities and accomplishments carried out by the Goddard Space Flight Center (GSFC)'s Flight Dynamics Analysis Branch (FDAB), Code 572, in support of flight projects and technology development initiatives in Fiscal Year (FY) 1999. The document is intended to serve as both an introduction to the type of support carried out by the FDAB (Flight Dynamics Analysis Branch), as well as a concise reference summarizing key analysis results and mission experience derived from the various mission support roles assumed over the past year. The major accomplishments in the FDAB in FY99 were: 1) Provided flight dynamics support to the Lunar Prospector and TRIANA missions among a variety of spacecraft missions; 2) Sponsored the Flight Mechanics Symposium; 3) Supported the Consultative Committee for Space Data Systems (CCSDS) workshops; 4) Performed numerous analyses and studies for future missions; 5) Started the Flight Dynamics Analysis Branch Lab for in-house mission analysis and support; and 6) Complied with all requirements in support of GSFC IS09000 certification.

Stengle, Thomas↗

ER-2 Observations of Precipitation Systems During TRMM-LBA

The NASA ER-2 performed numerous flights over precipitation systems in Rondonia, Brazil. The ER-2 carried a payload including the ER-2 Doppler Radar (EDOP), the Advanced Microwave Precipitation Radiometer (AMPR), the Lightning Instrument Package, and other instruments. This presentation will overview the types of data sets collected during TRMM-LBA with particular emphasis on EDOP measurements. Numerous cases of convection ranging from weak to very intense, were overflown by the ER-2. Two TRMM overpasses were coincident with ER-2 flights which allowed for intercomparisons between the Precipitation Radiometer (PR), EDOP, and the S-POL and TOGA ground-based radars. Preliminary results from this comparison will be presented as well as initial selection of case studies and efforts involving vertical motions in convection.

Heymsfield, Gerald↗

Optical Control of Intersubband Absorption in a Multiple Quantum Well-Embedded Semiconductor Microcravity

Optical intersubband response of a multiple quantum well (MQW)-embedded microcavity driven by a coherent pump field is studied theoretically. The n-type doped MQW structure with three subbands in the conduction band is sandwiched between a semi-infinite medium and a distributed Bragg reflector (DBR). A strong pump field couples the two upper subbands and a weak field probes the two lower subbands. To describe the optical response of the MQW-embedded microcavity, we adopt a semi-classical nonlocal response theory. Taking into account the pump-probe interaction, we derive the probe-induced current density associated with intersubband transitions from the single-particle density-matrix formalism. By incorporating the current density into the Maxwell equation, we solve the probe local field exactly by means of Green's function technique and the transfer-matrix method. We obtain an exact expression for the probe absorption coefficient of the microcavity. For a GaAs/Al(sub x)Ga(sub 1-x)As MQW structure sandwiched between a GaAs/AlAs DBR and vacuum, we performed numerical calculations of the probe absorption spectra for different parameters such as pump intensity, pump detuning, and cavity length. We find that the probe spectrum is strongly dependent on these parameters. In particular, we find that the combination of the cavity effect and the Autler-Townes effect results in a triplet in the optical spectrum of the MQW system. The optical absorption peak value and its location can be feasibly controlled by varying the pump intensity and detuning.

Liu, Ansheng↗

Comparison of Implicit Collocation Methods for the Heat Equation

We combine a high-order compact finite difference scheme to approximate spatial derivatives arid collocation techniques for the time component to numerically solve the two dimensional heat equation. We use two approaches to implement the collocation methods. The first one is based on an explicit computation of the coefficients of polynomials and the second one relies on differential quadrature. We compare them by studying their merits and analyzing their numerical performance. All our computations, based on parallel algorithms, are carried out on the CRAY SV1.

Kouatchou, Jules↗

The Response of Disks to Oscillatory Modes in Galaxies

Recent studies suggest that galaxies can oscillate in normal modes with essentially no damping over a Hubble time. These modes may play an important role in the structure and evolution of disk/halo systems. Motivated by the possibility that normal mode oscillations exist in real galaxies, we are investigating the response of galactic disks to halo oscillations. The goal of these investigations is to ascertain whether or not observational signatures exist for such oscillations. Our approach is to perform numerical experiments on the response of a self-gravitating disk to a time-varying halo potential. We assume that a significant fraction of the mass in a galaxy is in a dark halo. The halo oscillates and the luminous disk material responds to these oscillations. Preliminary results are reported for disks embedded in a radially oscillating gravitational potential. The equilibrium initial disk is represented by an exponential density profile. Considerable care was taken to build an initial disk model that was "stable" over long time scales. A control experiment was run with the disk in a static halo potential. The disk responds to the time-varying potential by developing a ring structure, which forms and disappears during each halo oscillation cycle. The density of stars becomes depressed in an annular region at the radius where the disk epicycle frequency is equal to the halo oscillation period. This pattern of response persists over time periods approaching a Hubble time. In the oscillating potential, a bar develops in the inner disk. This bar is absent when the halo remains static. Specific targets of this study include the implications for large-scale disk structure, the gas dynamical response of the interstellar medium in such systems, and the inflow of material into the central regions of the galaxy.

Smith, Bruce F.↗

The Fate of Ejecta from Hyperion

Ejecta from Saturn's moon Hyperion is subject to powerful perturbations from nearby Titan, which control its ultimate fate. We have performed numerical integrations to simulate a simplified system consisting of Saturn (including oblateness), Tethys, Dione, Titan, Hyperion, Iapetus, and the Sun (treated simply as a massive satellite). In addition, 1050 massless particles were ejected from Hyperion at five different points in its orbit. These particles started more or less evenly distributed over latitude and longitude, 1 km above Hyperion's mean radius, and were ejected radially outward at speeds 10\% faster than its escape speed. Only about 4\% of the particles survived for the 100,000-year course of the integration, while $\sim$8/% escaped from the Saturnian system. Titan accreted $\sim$77\% of all the particles, while Hyperion reaccreted only $\sim$5\%. This may help to account for Hyperion's rugged shape. Three particles hit Rhea and 2 hit Dione, but $\sim$5\% of the particles were removed when they penetrated within 150,000 km of Saturn. Most removals occurred within the first few thousand years. In general, ejecta from Hyperion are much more widely scattered than previously thought, and cross the orbits of all of the other classical satellites.

Lissauer, Jack J.↗

Building Muscles, Keeping Muscles: Protein Turnover During Space Flight

As we age we lose muscle mass and strength. The problem is a matter of use it or lose it and more - a fact to which any active senior can attest. An imbalance in the natural cycle of protein turnover may be a contributing factor to decreased muscle mass. But the answer is not so simple, since aging is associated with changes in hormones, activity levels, nutrition, and often, disease. The human body constantly uses amino acids to build muscle protein, which then breaks down and must be replaced. When protein turnover gets out of balance, so that more protein breaks down than the body can replace, the result is muscle loss. This is not just the bane of aging, however. Severely burned people may have difficulty building new muscle long after the burned skin has been repaired. Answers to why we lose muscle mass and strength - and how doctors can fix it - may come from space. Astronauts usually eat a well-balanced diet and maintain an exercise routine to stay in top health. During long-duration flight, they exercise regularly to reduce the muscle loss that results from being in a near-weightless environment. Despite these precautions, astronauts lose muscle mass and strength during most missions. They quickly recover after returning to Earth - this is a temporary condition in an otherwise healthy population. Members of the STS-107 crew are participating in a study of the effects of space flight, hormone levels, and stress on protein turnover. When we are under stress, the body responds with a change in hormone levels. Researchers hypothesize that this stress-induced change in hormones along with the near-weightlessness might result in the body synthesizing less muscle protein, causing muscles to lose their strength and size. Astronauts, who must perform numerous duties in a confined and unusual environment, experience some stress during their flight, making them excellent candidates for testing the researchers' hypothesis.

Arny Ferrando↗