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 109 records · Page 6

Pentium Pro inside: A treecode at 430 Gigaflops on ASCI Red - 1

As an entry for the 1997 Gordon Bell performance prize, we present results from two methods of solving the gravitational N-body problem on the Intel Teraflops system at Sandia National Laboratory (ASCI Red). The first method, an O(N2) algorithm, obtained 635 Gigaflops for a 1 million particle problem on 6800 Pentium Pro processors. The second solution method, a tree-code which scales as O(N log N), sustained 170 Gigaflops over a continuous 9.4 hour period on 4096 processors, integrating the motion of 322 million mutually interacting particles in a cosmology simulation, while saving over 100 Gigabytes of raw data. Additionally, the tree-code sustained 430 Gigaflops on 6800 processors for the first 5 time-steps of that simulation. This tree-code solution is approximately 105 times more efficient than the O(N2) algorithm for this problem. As an entry for the 1997 Gordon Bell price/performance prize, we present two calculations from the disciplines of astrophysics and fluid dynamics. The simulations were performed on two 16 Pentium Pro processor Beowulf-class computers (Loki and Hyglac) constructed entirely from commodity personal computer technology, at a cost of roughly $50k each in September, 1996. The price of an equivalent system in August 1997 is less than $30. At Los Alamos, Loki performed a gravitational tree-code N-body simulation of galaxy formation using 9.75 million particles, which sustained an average of 879 Mflops over a ten day period, and produced roughly 10 Gbytes of raw data.

Warren, M. S.↗

Computing Rotational, Rovibrational, and Vibrational Spectra for Astronomical Observations: High Accuracy Line Lists for High Temperatures, Limited Line Lists for Biosignature Molecules, and PAH Emission Spectra

Over the last several years, our group has been involved in developing approaches to compute highly accurate spectroscopic constants and vibrational frequencies for small transient molecules that may be used in the interpretation and assignment of high-resolution laboratory experiments as well as high-resolution astronomical spectra. Additionally, we have used the computed spectroscopic constants to simulate purely rotational and rovibrational spectra so that these may be compared directly with high- resolution astronomical observations, and we have worked on developing approaches that can be applied to much larger molecules, such as PAHs, where we can explicitly determine an harmonic corrections to vibrational modes as well as take into account intensity sharing due to resonances. Another part of our work in spectroscopic signatures involves computing highly accurate line lists for common molecules, such as CO2, SO2, and NH3, which occur in many astrophysical environments, including the atmospheres of exoplanets, and often need to have their lines identified in high- resolution observations in order to determine which lines are due to other molecules. In order to characterize the atmospheres of hot exoplanets, these line lists need to be very accurate and extend to very high energies. I will discuss our latest work in these areas of astrochemical spectroscopy research.

Lee, Timothy J.↗

Computing rovibrational, vibrational, and cascade emission spectra for comparison to astronomical observations

Over the last several years, our group has been involved in developing approaches to compute highly accurate spectroscopic constants and vibrational frequencies for small transient molecules that may be used in the interpretation and assignment of high-resolution laboratory experiments as well as high-resolution astronomical spectra. Additionally, we have used the computed spectroscopic constants to simulate purely rotational and rovibrational spectra so that these may be compared directly with high-resolution astronomical observations, and we have worked on developing approaches that can be applied to much larger molecules, such as polycyclic aromatic hydrocarbon (PAH) molecules, where we can explicitly determine anharmonic corrections to vibrational frequencies as well as take into account intensity sharing due to resonances. In addition, for PAH molecules we have used the anharmonic data to compute a library of temperature-dependent anharmonic vibrational spectra, with proper inclusion of polyad resonances, which can then be used to model the cascade emission spectra of PAH molecules – the type of spectra directly observed by astronomers. Another part of our work in spectroscopic signatures involves computing highly accurate line lists for common molecules, such as CO2, SO2, and NH3, which occur in many astrophysical environments, including the atmospheres of exoplanents, and often need to have their lines identified in high-resolution observations in order to determine which lines are due to other molecules. In order to characterize the atmospheres of hot exoplanets, these line lists need to be very accurate and extend to very high energies. I will discuss our latest work in these areas of astrochemical spectroscopy research.

Timothy J Lee↗

Visualization techniques to aid in the analysis of multi-spectral astrophysical data sets

This report describes our project activities for the period Sep. 1991 - Oct. 1992. Our activities included stabilizing the software system STAR, porting STAR to IDL/widgets (improved user interface), targeting new visualization techniques for multi-dimensional data visualization (emphasizing 3D visualization), and exploring leading-edge 3D interface devices. During the past project year we emphasized high-end visualization techniques, by exploring new tools offered by state-of-the-art visualization software (such as AVS3 and IDL4/widgets), by experimenting with tools still under research at the Department of Computer Science (e.g., use of glyphs for multidimensional data visualization), and by researching current 3D input/output devices as they could be used to explore 3D astrophysical data. As always, any project activity is driven by the need to interpret astrophysical data more effectively.

Domik, Gitta↗

Heavy ion beams in extended materials - Computational methods and experiment

The transport of heavy ion beams in extended materials is a problem of interest in accelerator and space shielding, radiation therapy, and astrophysical and radiobiological studies. The beam particles change their energy and direction of motion through atomic/molecular collisions and undergo occasional radical transformation in nuclear collision. In health physics applications, a heavy ion beam of initially well defined radiation quality is transformed into a complex mixture of diverse quality components after passing through a modest amount of material. This transformation of radiation quality must be understood to adequately explain the biological response of tissue to heavy ion radiation. A theoretical/experimental program to define an ion beam and its products in extended matter is described.

Wilson, J. W.↗

The statistics of gravitational lenses. III - Astrophysical consequences of quasar lensing

The method of Schmidt and Green (1983) for calculating the luminosity function of quasars is combined with gravitational-lensing theory to compute expected properties of lensed systems. Multiple quasar images produced by galaxies are of order 0.001 of the observed quasars, with the numbers over the whole sky calculated to be (0.86, 120, 1600) to limiting B magnitudes of (16, 19, 22). The amount of 'false evolution' is small except for an interesting subset of apparently bright, large-redshift objects for which minilensing by starlike objects may be important. Some of the BL Lac objects may be in this category, with the galaxy identified as the parent object really a foreground object within which stars have lensed a background optically violent variable quasar.

Ostriker, J. P.↗

Origins of interstellar and solar system: Carbonaceous materials

Carbon is a crucial atom in cosmochemistry. It is well-established that carbon is synthesized in stellar interiors after the main sequence, is ejected by red giants as small carbonaceous grains during their 'carbon star' phase, resides in the interstellar medium, and was later incorporated into the solar system. The mechanisms of carbon grain formation and later chemical processing are complex because, with only small thermodynamic differences, carbon can take on a bewildering variety of forms: diamond; oxides; carbides; graphite; aliphatic hydrocarbons; polycyclic aromatic hydrocarbons (PAH's); fullerenes; amorphous carbon; and other compounds. These are evidence for many of the forms of carbon found in astronomical observations. We seek to understand the possible astrophysical sites and conditions of the origins of different forms of carbon by combining state-of-the-art capabilities of carbon chemistry with astrophysical modeling. The work is a collaboration between Prof. Frenklach, a leading carbon materials scientist with both laboratory and computer modeling expertise and Prof. Feigelson, an astrophysicist with interests in star formation. The largest effort under this grant was devoted to developing this concept into a comprehensive quantitative model. In addition to explaining the astronomical properties of red giants producing carbonaceous grains, our model also can incorporate recent meteoritic findings. Finally, our induced nucleation grain formation model provides a natural explanation for the widespread presence of PAH emission bands in the Galactic interstellar medium. A brief synopsis of other activities sponsored under this grant and a list of publications from this grant is included.

Feigelson, Eric D.↗

Data Analysis Challenges for Multi-Messenger Astrophysics

Recent multi-messenger observations of gravitational-wave and high-energy neutrino sources together with electromagnetic signatures have opened new ways of observing the Universe. These promise a future in which physics and astronomy will be advanced by combining observations and data from across the electromagnetic spectrum with gravitational waves and neutrinos. We consider the challenges the field is facing in fully utilizing data for multi-messenger astrophysics. Such data come from heterogeneous detector networks and standards, and their analysis is often time-critical to guide further observations. In this area, science capabilities depend on the interplay among observation, theory and computational/modeling work. Advances in data science and computing present additional opportunities and considerations in analyzing such data. We invited ADASS participants to a Birds of a Feather session to engage in discussion on the challenges and opportunities in data analysis for multimessenger astrophysics.

Peter S Shawhan↗

Stochastic Template Bank for Gravitational Wave Searches for Precessing Neutron Star-Black Hole Coalescence Events

Gravitational wave searches to date have largely focused on non-precessing systems. Including precession effects greatly increases the number of templates to be searched over. This leads to a corresponding increase in the computational cost and can increase the false alarm rate of a realistic search. On the other hand, there might be astrophysical systems that are entirely missed by non-precessing searches. In this paper we consider the problem of constructing a template bank using stochastic methods for neutron star-black hole binaries allowing for precession, but with the restrictions that the total angular momentum of the binary is pointing toward the detector and that the neutron star spin is negligible relative to that of the black hole. We quantify the number of templates required for the search, and we explicitly construct the template bank. We show that despite the large number of templates, stochastic methods can be adapted to solve the problem. We quantify the parameter space region over which the non-precessing search might miss signals.

gravitational wave↗

Plasma-environment effects on K lines of astrophysical interest III. IPs, K thresholds, radiative rates, and Auger widths in Fe ix – Fe xvi

Aims. In the context of black-hole accretion disks, we aim to compute the plasma-environment e ects on the atomic parameters used to model the decay of K-vacancy states in moderately charged iron ions, namely Fe ix – Fe xvi. Methods. We used the fully relativistic multiconfiguration Dirac–Fock method approximating the plasma electron–nucleus and electron–electron screenings with a time-averaged Debye–Hückel potential. Results. We report modified ionization potentials, K-threshold energies, wavelengths, radiative emission rates, and Auger widths for plasmas characterized by electron temperatures and densities in the ranges 10^(5)-10^(7) K and 10^(18)-10^(22)/cu. cm. Conclusions. This study confirms that the high-resolution X-ray spectrometers onboard the future XRISM and Athena space missions will be capable of detecting the lowering of the K edges of these ions due to the extreme plasma conditions occurring in accretion disks around compact objects.

J. Deprince↗

Theoretical study of the X-ray emission from astrophysical shock waves

Theoretical X-ray emission spectra are needed to interpret the X-ray emission observed by many low and moderate resolution X-ray instruments, and to provide diagnosis of physical conditions for high resolution spectra. Over the past decade, a set of model codes which compute the X-ray and XUV emission for a wide set of physical conditions, including high or low densities, photoionized gas, and time-dependent ionization balance was developed. In the past year, the atomic rate coefficients in the code was improved. Further capabilities were added, and applied to several astrophysical problems.

Raymond, J.↗

The MSFC/J70 orbital atmosphere model and the data bases for the MSFC solar activity prediction technique

The MSFC/J70 Orbital Atmospheric Density Model, a modified version of the Smithsonian Astrophysical Observatory Jacchia 1970 model is explained. The algorithms describing the MSFC/J70 model are included as well as listing of the computer program. The 13-month smoothed values of solar flux (F sub 10.7) and geomagnetic index (S sub p), which are required as inputs for the MSFC/J70 model, are also included and discussed.

Johnson, D. L.↗

Tests and comparisons of satellite-derived geoids with Skylab altimeter data

During the Skylab 4 mission, the S-193 radar altimeter was operated nearly continuously for a revolution around the world on Jan. 31, 1974. This direct measurement to the sea surface has provided an independent basis for the evaluation of the precision of global geoids computed from satellite-derived earth gravity models. This paper presents comparisons between the Skylab data and several recent gravity models published by Goddard Space Flight Center, the Smithsonian Astrophysical Observatory, and the National Oceanic and Atmospheric Administration. The differences between the altimeter geoid and the satellite geoids were as large as 20 m, rms values ranging from 8 to 10 m. These differences also indicated a systematic long-wavelength variation (about 100 deg) not related to error in the Skylab orbits. Truncation of the models to degree and order 8 did not eliminate the long-wavelength variation, but in every case the rms agreement between the satellite geoids and the altimeter geoid was slightly improved. Orbits computed with the truncated models were found to be inferior to those computed with the complete models.

Marsh, J. G.↗

Visualization techniques to aid in the analysis of multi-spectral 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.

Brugel, Edward W.↗

The Prospective Search for Highly Ionized Technetium in Hot (pre-) White Dwarfs

The discovery of technetium (Tc) in the atmospheres of red giants by Merrill (1952, ApJ, 116, 21) constituted convincing proof that s-process nucleosynthesis is indeed occurring in evolved stars. In principle, Tc should still be present in the atmospheres of hot post-AGB stars and (pre-) white dwarfs although, due to radioactive decay, it should be present in decreasing quantities along post-AGB evolution. The recent discovery of a large number of trans-iron group elements in hot white dwarfs with atomic numbers in the range A = 30-56 (Zn to Ba) raises the prospect that Tc (A = 43) may also be detected. However, this is currently not feasible because no atomic data exist for ionization stages beyond Tc ii. As an initial step, we calculated atomic energy levels and oscillator strengths of Tc iv-vi and used these data to compute non-local thermodynamic equilibrium (NLTE) model atmospheres to estimate at which minimum abundance level Tc could be detected. We show that Tc lines can be found in ultraviolet spectra of hot white dwarfs provided Tc is as abundant as other detected trans-Fe elements. We find that radiative levitation can keep Tc in large, easily detectable quantities in the atmosphere. A direct identification of Tc lines is still not feasible because wavelength positions cannot be computed with necessary precision. Laboratory measurements are necessary to overcome this problem. Our results suggest that such efforts are beneficial to the astrophysical community.

prospective↗

Magnetic Field Effects on Plasma Plumes

Here, we will discuss our numerical studies of plasma jets and loops, of basic interest for plasma propulsion and plasma astrophysics. Space plasma propulsion systems require strong guiding magnetic fields known as magnetic nozzles to control plasma flow and produce thrust. Propulsion methods currently being developed that require magnetic nozzles include the VAriable Specific Impulse Magnetoplasma Rocket (VASIMR) [1] and magnetoplasmadynamic thrusters. Magnetic nozzles are functionally similar to de Laval nozzles, but are inherently more complex due to electromagnetic field interactions. The two crucial physical phenomenon are thrust production and plasma detachment. Thrust production encompasses the energy conversion within the nozzle and momentum transfer to a spacecraft. Plasma detachment through magnetic reconnection addresses the problem of the fluid separating efficiently from the magnetic field lines to produce maximum thrust. Plasma jets similar to those of VASIMR will be studied with particular interest in dual jet configurations, which begin as a plasma loops between two nozzles. This research strives to fulfill a need for computational study of these systems and should culminate with a greater understanding of the crucial physics of magnetic nozzles with dual jet plasma thrusters, as well as astrophysics problems such as magnetic reconnection and dynamics of coronal loops.[2] To study this problem a novel, hybrid kinetic theory and single fluid magnetohydrodynamic (MHD) solver known as the Magneto-Gas Kinetic Method is used.[3] The solver is comprised of a "hydrodynamic" portion based on the Gas Kinetic Method and a "magnetic" portion that accounts for the electromagnetic behaviour of the fluid through source terms based on the resistive MHD equations. This method is being further developed to include additional physics such as the Hall effect. Here, we will discuss the current level of code development, as well as numerical simulation results

Ebersohn, F.↗

Information technologies for astrophysics circa 2001

It is easy to extrapolate current trends to see where technologies relating to information systems in astrophysics and other disciplines will be by the end of the decade. These technologies include mineaturization, multiprocessing, software technology, networking, databases, graphics, pattern computation, and interdisciplinary studies. It is easy to see what limits our current paradigms place on our thinking about technologies that will allow us to understand the laws governing very large systems about which we have large datasets. Three limiting paradigms are saving all the bits collected by instruments or generated by supercomputers; obtaining technology for information compression, storage and retrieval off the shelf; and the linear mode of innovation. We must extend these paradigms to meet our goals for information technology at the end of the decade.

Denning, Peter J.↗

Information technologies for astrophysics circa 2001

It is easy to extrapolate current trends to see where technologies relating to information systems in astrophysics and other disciplines will be by the end of the decade. These technologies include miniaturization, multiprocessing, software technology, networking, databases, graphics, pattern computation, and interdisciplinary studies. It is less easy to see what limits our current paradigms place on our thinking about technologies that will allow us to understand the laws governing very large systems about which we have large data sets. Three limiting paradigms are as follows: saving all the bits collected by instruments or generated by supercomputers; obtaining technology for information compression, storage, and retrieval off the shelf; and the linear model of innovation. We must extend these paradigms to meet our goals for information technology at the end of the decade.

Denning, Peter J.↗