Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational astrophysics”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Feasibility study of an orbiting laboratory for testing CSI technology

A concept for an orbiting laboratory for testing Controls-Structures Integration (CSI) technology is described. The CSI-Star concept reflects a lower cost, higher risk approach. The concept supports demonstration and validation testing for critical CSI technologies at a cost of $20M to $26M with a 1-year reliability of approximately 0.9. The Ball Aerospace QuickStar bus is the carrier for the CSI test article. QuickStar is launched as a secondary payload on the McDonnell Douglas Delta 2. The QuickStar/Delta 2 approach is flight proven. The CSI test article is a 20 foot, 1 Hz, truss beam which is deployed from the QuickStar bus. The test article is well instrumented for quality system identification. The laboratory provides three layers of active control consisting of global vibration suppression along the truss beam, vibration isolation between the beam and instrument platforms, and vibration compensation through the use of gimbaled platforms which point lasers relative to optical sensor targets. The configuration simulates the dynamics of multi-instrument science platforms such as those of the Earth Observation System (EOS) while maintaining strong ties to astrophysics missions such as the Optical Interferometer. Uplink/downlink services and a reprogrammable computer provide flexibility for long-term investigations by members of the CSI community (NASA, DoD, academia, and industry). CSI-Star fills the gap between short-term experiments, which have been conducted primarily on the Shuttle, and future science missions which require the technology. The on-orbit maturity of CSI technology must be established to obtain acceptance by project managers and to promote injection of the technology into future science missions.

Bicos, Andrew S.↗
NASA NTRS2000

Photoionization and Recombination

Theoretically self-consistent calculations for photoionization and (e + ion) recombination are described. The same eigenfunction expansion for the ion is employed in coupled channel calculations for both processes, thus ensuring consistency between cross sections and rates. The theoretical treatment of (e + ion) recombination subsumes both the non-resonant recombination ("radiative recombination"), and the resonant recombination ("di-electronic recombination") processes in a unified scheme. In addition to the total, unified recombination rates, level-specific recombination rates and photoionization cross sections are obtained for a large number of atomic levels. Both relativistic Breit-Pauli, and non-relativistic LS coupling, calculations are carried out in the close coupling approximation using the R-matrix method. Although the calculations are computationally intensive, they yield nearly all photoionization and recombination parameters needed for astrophysical photoionization models with higher precision than hitherto possible, estimated at about 10-20% from comparison with experimentally available data (including experimentally derived DR rates). Results are electronically available for over 40 atoms and ions. Photoionization and recombination of He-, and Li-like C and Fe are described for X-ray modeling. The unified method yields total and complete (e+ion) recombination rate coefficients, that can not otherwise be obtained theoretically or experimentally.

Nahar, Sultana N.↗
NASA NTRS2004

Spectroscopic Determination of the Physical Conditions in Hot Optically Thin Sources

Astrophysical Plasma Emission Database - APED consists of atomic data, primarily theoretical, needed to calculate X-ray through optical emission spectra of hot thermal plasmas. These data are supplemented by experimental values, such as laboratory plasma wave- lengths, and are validated by experiment only rarely. Thus the comparison of predicted spectra to astrophysical spectra serves as invaluable feedback on the quality of the models. Following up on unfavorable comparisons, we have communicated atomic data needs to atomic theorists and experimental plasma physicists so that improvements can be made to the database. The database is now reasonably complete from approx. 1 to 30 A, and from - 90 to 180 A, with significant gaps remaining in the range 30 to 90 A. Astrophysical Plasma Emission Code - APEC uses the APED data as inputs and computes level populations and line emissivities using a rate matrix solver. Tabulated models for a fine grid of temperatures and densities can be applied to observational data. The current public version of APEC output uses an input ionization balance model, assuming CIE; however, we have a working version to calculate the ionization state of the gas. We will begin testing non-equilibrium ionization (NEI) models over the next year - we are now finalizing the atomic database needed to perform the ionization state calculations. At that time, the code will be made public. We have developed optimized error codes, additional testing protocols and substantial documentation, in preparation for this public release. Secondarily, we have begun to investigate the use of APED and APEC for X-ray photoionized plasma, beginning with opacity modeling (e.g. warm absorbers). Application to bright AGN Chandra grating spectra suggests that the accurate wavelengths in APED are of great importance to fitting the data. We are exploring additional subroutines to APEC to extend its usefulness for X-ray photoionized plasma, e.g. to include collisional or photo-excitation of metastable levels for absorption.

Oliversen, Ronald J.↗
NASA NTRS1995

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

Turbulent transport in the solar nebula

It is likely that turbulence played a major role in the evolution of the solar nebula, which is the flattened disk of dust and gas out of which our solar system formed. Relevant turbulent processes include the transport of angular momentum, mass, and heat, which were critically important to the formation of the solar system. This research will break ground in the modeling of compressible turbulence and its effects on the evolution of the solar nebula. The computational techniques which were developed should be of interest to researchers studying other astrophysical disk systems (e.g., active galactic nuclei), as well as turbulence modelers outside the astrophysics community.

Thompson, Kevin W.↗
NASA NTRS1995

Experimental Evaluation and Workload Characterization for High-Performance Computer Architectures

This research is conducted in the context of the Joint NSF/NASA Initiative on Evaluation (JNNIE). JNNIE is an inter-agency research program that goes beyond typical.bencbking to provide and in-depth evaluations and understanding of the factors that limit the scalability of high-performance computing systems. Many NSF and NASA centers have participated in the effort. Our research effort was an integral part of implementing JNNIE in the NASA ESS grand challenge applications context. Our research work under this program was composed of three distinct, but related activities. They include the evaluation of NASA ESS high- performance computing testbeds using the wavelet decomposition application; evaluation of NASA ESS testbeds using astrophysical simulation applications; and developing an experimental model for workload characterization for understanding workload requirements. In this report, we provide a summary of findings that covers all three parts, a list of the publications that resulted from this effort, and three appendices with the details of each of the studies using a key publication developed under the respective work.

El-Ghazawi, Tarek A.↗
NASA NTRS2010

Modeling of Photoionized Plasmas

In this paper I review the motivation and current status of modeling of plasmas exposed to strong radiation fields, as it applies to the study of cosmic X-ray sources. This includes some of the astrophysical issues which can be addressed, the ingredients for the models, the current computational tools, the limitations imposed by currently available atomic data, and the validity of some of the standard assumptions. I will also discuss ideas for the future: challenges associated with future missions, opportunities presented by improved computers, and goals for atomic data collection.

Kallman, Timothy R.↗
NASA NTRS1994

Gas dynamics for accretion disk simulations

The behavior of accretion disks can largely be understood in terms of the basic physical processes of mass, energy, and momentum conservation. Despite this, detailed modeling of these systems using modern computational techniques is challenging and controversial. Disturbing differences exist between methods used widely in astrophysics, namely Eulerian finite-difference techniques and particle codes such as SPH. Therefore neither technique is fully satisfactory for accretion disk simulations. This paper describes a new fully Lagrangian method designed to resolve these difficulties.

Whitehurst, R.↗
NASA NTRS1994

Power spectrum, correlation function, and tests for luminosity bias in the CfA redshift survey

We describe and apply a method for directly computing the power spectrum for the galaxy distribution in the extension of the Center for Astrophysics Redshift Survey. Tests show that our technique accurately reproduces the true power spectrum for k greater than 0.03 h Mpc(exp -1). The dense sampling and large spatial coverage of this survey allow accurate measurement of the redshift-space power spectrum on scales from 5 to approximately 200 h(exp -1) Mpc. The power spectrum has slope n approximately equal -2.1 on small scales (lambda less than or equal 25 h(exp -1) Mpc) and n approximately -1.1 on scales 30 less than lambda less than 120 h(exp -1) Mpc. On larger scales the power spectrum flattens somewhat, but we do not detect a turnover. Comparison with N-body simulations of cosmological models shows that an unbiased, open universe CDM model (OMEGA h = 0.2) and a nonzero cosmological constant (CDM) model (OMEGA h = 0.24, lambda(sub zero) = 0.6, b = 1.3) match the CfA power spectrum over the wavelength range we explore. The standard biased CDM model (OMEGA h = 0.5, b = 1.5) fails (99% significance level) because it has insufficient power on scales lambda greater than 30 h(exp -1) Mpc. Biased CDM with a normalization that matches the Cosmic Microwave Background (CMB) anisotropy (OMEGA h = 0.5, b = 1.4, sigma(sub 8) (mass) = 1) has too much power on small scales to match the observed galaxy power spectrum. This model with b = 1 matches both Cosmic Background Explorer Satellite (COBE) and the small-scale power spect rum but has insufficient power on scales lambda approximately 100 h(exp -1) Mpc. We derive a formula for the effect of small-scale peculiar velocities on the power spectrum and combine this formula with the linear-regime amplification described by Kaiser to compute an estimate of the real-space power spectrum. Two tests reveal luminosity bias in the galaxy distribution: First, the amplitude of the pwer spectrum is approximately 40% larger for the brightest 50% of galaxies in volume-limited samples that have M(sub lim) greater than M*. This bias in the power spectrum is independent of scale, consistent with the peaks-bias paradigm for galaxy formation. Second, the distribution of local density around galaxies shows that regions of moderate and high density contain both very bright (M less than M* = -19.2 + 5 log h) and fainter galaxies, but that voids preferentially harbor fainter galaxies (approximately 2 sigma significance level).

Park, Changbom↗
NASA NTRS2016

Computing Highly Accurate Spectroscopic Line Lists for Characterization of Planetary Atmospheres: CO2 and SO2 Line Lists Needed for Modeling Venus

Over the last decade, it has become apparent that the most effective approach for determining highly accurate rotational and rovibrational line lists for molecules of interest in planetary atmospheres and other astrophysical environments is through a combination of highresolution laboratory experiments coupled with state-of-the art ab initio quantum chemistry methods. The approach involves computing the most accurate potential energy surface (PES) possible using state-of-the art electronic structure methods, followed by computing rotational and rovibrational energy levels using an exact variational method to solve the nuclear Schrödinger equation. Then, reliable experimental data from high-resolution experiments is used to refine the ab initio PES in order to improve the accuracy of the computed energy levels and transition energies. From the refinement step, we have been able to achieve an accuracy of approximately 0.015 cm-1 for rovibrational transition energies, and even better for purely rotational transitions. This combined "experiment / theory" approach allows for determination of essentially a complete line list, with hundreds of millions of transitions, and having the transition energies and intensities be highly accurate. Our group has successfully applied this approach to determine highly accurate line lists for NH3, CO2 and isotopologues, and SO2 and isotopologues. Here I will report our latest results for CO2 and SO2 including all isotopologues. Comparisons to the available data in HITRAN2012 and other available databases will be shown, though we note that our line lists for SO2 are significantly more complete than any other databases. Since it is important to span a large temperature range in order to model the spectral signature of Venus as well as exoplanets, we will demonstrate how the spectra change on going from low temperatures (100 K) to higher temperatures (500 K to 1500 K).

potential energy surface (PES)5405 PLANETARY SCIEN↗
NASA NTRS1982

Highly versatile computer-controlled television detector system

A description is presented of a television detector system which has been designed to accommodate a wide range of applications. It is currently being developed for use in X-ray diffraction, X-ray astrophysics, and electron microscopy, but it is also well suited for astronomical observations. The image can be integrated in a large, high-speed memory system, in the memory of a computer system, or the target of the TV tube or CCD array. The detector system consists of a continuously scanned, intensified SIT vidicon with scan and processing electronics which generate a digital image that is integrated in the detector memory. Attention is given to details regarding the camera system, scan control and image processing electronics, the memory system, and aspects of detector performance.

Kalata, K.↗
NASA NTRS2003

Research in Computational Astrobiology

We report on several projects in the field of computational astrobiology, which is devoted to advancing our understanding of the origin, evolution and distribution of life in the Universe using theoretical and computational tools. Research projects included modifying existing computer simulation codes to use efficient, multiple time step algorithms, statistical methods for analysis of astrophysical data via optimal partitioning methods, electronic structure calculations on water-nuclei acid complexes, incorporation of structural information into genomic sequence analysis methods and calculations of shock-induced formation of polycylic aromatic hydrocarbon compounds.

Chaban, Galina↗
NASA NTRS2021

Modeling of the Tau and Muon Neutrino-induced Optical Cherenkov Signals from Upward-moving Extensive Air Showers

We present a detailed modeling and computation methodology to determine the optical Cherenkovsignals produced by upward-moving extensive air showers (EASs) induced byτ-leptons and muons,sourced from the interaction of high-energy astrophysical neutrinos interacting in the Earth. Fol-lowing and extending the physics modeling and Cherenkov signal simulations performed in [1], thisscheme encompasses a new, state-of-the-art computation of the muon neutrino propagation insidethe Earth and the contribution to theτ-lepton muon decay channel. The modeling takes into ac-count all possibleτ-lepton decay and muon energy loss channels that feed the optical Cherenkovemission, produced by both tau and muon initiated EASs. The EAS modeling uses the electronenergy, angular, and lateral distributions in the EAS and their evolution as well as the wavelengthdependence of the Cherenkov emission and its atmospheric attenuation. The results presented hereare focused on the detection capabilities of sub-orbital (balloon-borne) and orbital (satellite) basedinstruments. The latter case was calculated for POEMMAato compare to that presented in [1],specifically including the muon-decay channel ofτ-leptons and the muonic EAS Cherenkov signalfrom muon neutrino interactions in the Earth. By detailing all these individual contributions tothe optical Cherenkov emission and detection, we show how the ensemble that includes muonicchannels provides a large detection capability for space-based, high-energy cosmic neutrino detec-tion. Specifically, we show that for neutrino energies∼<10 PeV, the upward-EAS sensitivity due tomuon neutrino interactions in the Earth begin to dominate over that for tau neutrino interactions,effectively extending the neutrino sensitivity to lower energies.

A L Cummings↗
NASA NTRS1993

Visualization techniques to aid in the analysis of multispectral astrophysical data sets

The goal of this project was to support the scientific analysis of multi-spectral astrophysical data by means of scientific visualization. Scientific visualization offers its greatest value if it is not used as a method separate or alternative to other data analysis methods but rather in addition to these methods. Together with quantitative analysis of data, such as offered by statistical analysis, image or signal processing, visualization attempts to explore all information inherent in astrophysical data in the most effective way. Data visualization is one aspect of data analysis. Our taxonomy as developed in Section 2 includes identification and access to existing information, preprocessing and quantitative analysis of data, visual representation and the user interface as major components to the software environment of astrophysical data analysis. In pursuing our goal to provide methods and tools for scientific visualization of multi-spectral astrophysical data, we therefore looked at scientific data analysis as one whole process, adding visualization tools to an already existing environment and integrating the various components that define a scientific data analysis environment. As long as the software development process of each component is separate from all other components, users of data analysis software are constantly interrupted in their scientific work in order to convert from one data format to another, or to move from one storage medium to another, or to switch from one user interface to another. We also took an in-depth look at scientific visualization and its underlying concepts, current visualization systems, their contributions and their shortcomings. The role of data visualization is to stimulate mental processes different from quantitative data analysis, such as the perception of spatial relationships or the discovery of patterns or anomalies while browsing through large data sets. Visualization often leads to an intuitive understanding of the meaning of data values and their relationships by sacrificing accuracy in interpreting the data values. In order to be accurate in the interpretation, data values need to be measured, computed on, and compared to theoretical or empirical models (quantitative analysis). If visualization software hampers quantitative analysis (which happens with some commercial visualization products), its use is greatly diminished for astrophysical data analysis. The software system STAR (Scientific Toolkit for Astrophysical Research) was developed as a prototype during the course of the project to better understand the pragmatic concerns raised in the project. STAR led to a better understanding on the importance of collaboration between astrophysicists and computer scientists. Twenty-one examples of the use of visualization for astrophysical data are included with this report. Sixteen publications related to efforts performed during or initiated through work on this project are listed at the end of this report.

Brugel, E. W.↗
NASA NTRS2020

An Efficient Focal Plane Alignment Methodology with Application to the ASTERIA Nanosatellite Space Telescope

The process of aligning an instrument focal plane toprecursor optics such as a telescope can be a difficult and timeconsuming process. The final alignment must achieve sufficientperformance at multiple locations across the field of view. Inthe case of instruments with infinite conjugates, such as a spacetelescope, there are practical challenges in illuminating the system at multiple off-axis field points simultaneously. Iterationduring the alignment process can cause a lengthy search for localminima within the performance cost function.This paper presents a focal plane alignment methodology thatavoids the above difficulties and obtains the desired level ofoptical performance in an efficient manner. The approach isbased on systematically measuring spot size across a set of fieldpoints and defocus positions, and then using a least squaresfit to determine the optimal location for the focal plane. Thisapproach has the additional benefit of revealing the amount offield curvature present in the instrument. It also provides theanalyst with a direct calculation of the shims needed to place thefocal plane at the computed best-fit location.This approach was successfully applied to the Arcsecond SpaceTelescope Enabling Research In Astrophysics (ASTERIA) payload. ASTERIA is a cubesat mission that was deployed intolow-Earth orbit in November 2017 and operated until December2019. In addition to presenting the theoretical basis for themethodology, this paper will present laboratory measurementsobtained during the payload alignment campaign. This approach is applicable to future space-borne optical instrumentsthat require an efficient methodology for focal plane alignmentwith limited cost or schedule resources.

Smith, Matthew W↗
NASA NTRS1987

Numerical tests of a modified full implicit continuous Eulerian (FICE) scheme with projected normal characteristic boundary conditions for MHD flows

A numerical method has been developed based on a modified full implicit continuous Eulerian (FICE) scheme and projected normal characteristic boundary conditions for simulating MHD flows which undergo a long process of evolution. An astrophysical flow is chosen for illustration of this procedure, and numerical tests are made to verify the computational stability and physically realistic solution. Three computational tests have been accomplished; they are tests of solving methods, characteristic boundary condition, and time steps. The tests show that the program from the modified FICE scheme with proper boundary conditions and time steps can be made numerically stable for a time long enough to obtain physically plausible solutions.

Wu, S. T.↗