Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Athena”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

LOCAL SIMULATIONS OF INSTABILITIES DRIVEN BY COMPOSITION GRADIENTS IN THE ICM

The distribution of helium in the intracluster medium (ICM) permeating galaxy clusters is not well constrained due to the very high plasma temperature. Therefore, the plasma is often assumed to be homogeneous. A nonuniform helium distribution can, however, lead to biases when measuring key cluster parameters. This has motivated one-dimensional models that evolve the ICM composition assuming that the effects of magnetic fields can be parameterized or ignored. Such models for nonisothermal clusters show that helium can sediment in the cluster core, leading to a peak in concentration offset from the cluster center. The resulting profiles have recently been shown to be linearly unstable when the weakly collisional character of the magnetized plasma is considered. In this paper, we present a modified version of the MHD code Athena, which makes it possible to evolve a weakly collisional plasma subject to a gravitational field and stratified in both temperature and composition. We thoroughly test our implementation and confirm excellent agreement against several analytical results. In order to isolate the effects of composition, in this initial study we focus our attention on isothermal plasmas. We show that plasma instabilities, feeding off gradients in composition, can induce turbulent mixing and saturate by rearranging magnetic field lines and alleviating the composition gradient. Composition profiles that increase with radius lead to instabilities that saturate by driving the average magnetic field inclination to roughly 45°. We speculate that this effect may alleviate the core insulation observed in homogeneous settings, with potential consequences for the associated cooling flow problem.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

IRON OPACITY BUMP CHANGES THE STABILITY AND STRUCTURE OF ACCRETION DISKS IN ACTIVE GALACTIC NUCLEI

Accretion disks around supermassive black holes have regions where the Rosseland mean opacity can be larger than the electron scattering opacity due to the large number of bound–bound transitions in iron. We study the effects of this iron opacity “bump” on the thermal stability and vertical structure of radiation-pressure-dominated accretion disks, utilizing three-dimensional radiation magnetohydrodynamic (MHD) simulations in the local shearing box approximation. The simulations self-consistently calculate the heating due to MHD turbulence caused by magneto-rotational instability and radiative cooling by using the radiative transfer module based on a variable Eddington tensor in Athena. For a 5 × 10{sup 8} solar mass black hole with ∼3% of the Eddington luminosity, a model including the iron opacity bump maintains its structure for more than 10 thermal times without showing significant signs of thermal runaway. In contrast, if only electron scattering and free–free opacity are included as in the standard thin disk model, the disk collapses on the thermal timescale. The difference is caused by a combination of (1) an anti-correlation between the total optical depth and the midplane pressure, and (2) enhanced vertical advective energy transport. These results suggest that the iron opacity bump may have a strong impact on the stability and structure of active galactic nucleus (AGN) accretion disks, and may contribute to a dependence of AGN properties on metallicity. Since this opacity is relevant primarily in UV emitting regions of the flow, it may help to explain discrepancies between observation and theory that are unique to AGNs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Novel Machine Learning Approach to Disentangle Multitemperature Regions in Galaxy Clusters

The hot intracluster medium (ICM) surrounding the heart of galaxy clusters is a complex medium that comprises various emitting components. Although previous studies of nearby galaxy clusters, such as the Perseus, the Coma, or the Virgo cluster, have demonstrated the need for multiple thermal components when spectroscopically fitting the ICM’s X-ray emission, no systematic methodology for calculating the number of underlying components currently exists. In turn, underestimating or overestimating the number of components can cause systematic errors in the emission parameter estimations. In this paper, we present a novel approach to determining the number of components using an amalgam of machine learning techniques. Synthetic spectra containing a various number of underlying thermal components were created using well-established tools available from the Chandra X-ray Observatory. The dimensions of the training set was initially reduced using principal component analysis and then categorized based on the number of underlying components using a random forest classifier. Our trained and tested algorithm was subsequently applied to Chandra X-ray observations of the Perseus cluster. Our results demonstrate that machine learning techniques can efficiently and reliably estimate the number of underlying thermal components in the spectra of galaxy clusters, regardless of the thermal model (MEKAL versus APEC). We also confirm that the core of the Perseus cluster contains a mix of differing underlying thermal components. We emphasize that although this methodology was trained and applied on Chandra X-ray observations, it is readily portable to other current (e.g., XMM-Newton, eROSITA) and upcoming (e.g., Athena, Lynx, XRISM) X-ray telescopes. The code is publicly available at https://github.com/XtraAstronomy/Pumpkin.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ON HELIUM MIXING IN QUASI-GLOBAL SIMULATIONS OF THE INTRACLUSTER MEDIUM

The assumption of a spatially uniform helium distribution in the intracluster medium (ICM) can lead to biases in the estimates of key cluster parameters if composition gradients are present. The helium concentration profile in galaxy clusters is unfortunately not directly observable. Current models addressing the putative sedimentation are one-dimensional and parametrize the presence of magnetic fields in a crude way, ignoring the weakly collisional, magnetized nature of the medium. When these effects are considered, a wide variety of instabilities can play an important role in the plasma dynamics. In a series of recent papers, we have developed the local, linear theory of these instabilities and addressed their nonlinear development with a modified version of Athena. Here, we extend our study by developing a quasi-global approach that we use to simulate the mixing of helium as induced by generalizations of the heat-flux-driven buoyancy instability (HBI) and the magnetothermal instability, which feed off thermal and composition gradients. In the inner region of the ICM, mixing can occur over a few gigayears, after which the average magnetic field inclination angle is ∼30°–50°, resulting in an averaged Spitzer parameter higher by about 20% than the value obtained in homogeneous simulations. In the cluster outskirts the instabilities are rather inefficient, due to the shallow gradients. This suggests that composition gradients in cluster cores might be shallower than one-dimensional models predict. More quantitative statements demand more refined models that can incorporate the physics driving the sedimentation process and simultaneously account for the weakly collisional nature of the plasma.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Multiscale Pressure-Balanced Structures in Three-dimensional Magnetohydrodynamic Turbulence

Observations of solar wind turbulence indicate the existence of multiscale pressure-balanced structures (PBSs) in the solar wind. In this work, we conduct a numerical simulation to investigate multiscale PBSs and in particular their formation in compressive magnetohydrodynamic turbulence. By the use of the higher-order Godunov code Athena, a driven compressible turbulence with an imposed uniform guide field is simulated. The simulation results show that both the magnetic pressure and the thermal pressure exhibit a turbulent spectrum with a Kolmogorov-like power law, and that in many regions of the simulation domain they are anticorrelated. The computed wavelet cross-coherence spectra of the magnetic pressure and the thermal pressure, as well as their space series, indicate the existence of multiscale PBSs, with the small PBSs being embedded in the large ones. These multiscale PBSs are likely to be related to the highly oblique-propagating slow-mode waves, as the traced multiscale PBS is found to be traveling in a certain direction at a speed consistent with that predicted theoretically for a slow-mode wave propagating in the same direction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Three-phase Interstellar Medium in Galaxies Resolving Evolution with Star Formation and Supernova Feedback (TIGRESS): Algorithms, Fiducial Model, and Convergence

We introduce TIGRESS, a novel framework for multi-physics numerical simulations of the star-forming interstellar medium (ISM) implemented in the Athena MHD code. The algorithms of TIGRESS are designed to spatially and temporally resolve key physical features, including: (1) the gravitational collapse and ongoing accretion of gas that leads to star formation in clusters; (2) the explosions of supernovae (SNe), both near their progenitor birth sites and from runaway OB stars, with time delays relative to star formation determined by population synthesis; (3) explicit evolution of SN remnants prior to the onset of cooling, which leads to the creation of the hot ISM; (4) photoelectric heating of the warm and cold phases of the ISM that tracks the time-dependent ambient FUV field from the young cluster population; (5) large-scale galactic differential rotation, which leads to epicyclic motion and shears out overdense structures, limiting large-scale gravitational collapse; (6) accurate evolution of magnetic fields, which can be important for vertical support of the ISM disk as well as angular momentum transport. We present tests of the newly implemented physics modules, and demonstrate application of TIGRESS in a fiducial model representing the solar neighborhood environment. We use a resolution study to demonstrate convergence and evaluate the minimum resolution Δx required to correctly recover several ISM properties, including the star formation rate, wind mass-loss rate, disk scale height, turbulent and Alfvénic velocity dispersions, and volume fractions of warm and hot phases. For the solar neighborhood model, all these ISM properties are converged at Δx⩽8 pc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Modeling UV Radiation Feedback from Massive Stars. I. Implementation of Adaptive Ray-tracing Method and Tests [Modeling UV Radiation Feedback from Massive Stars: I. Implementation of Adaptive Ray Tracing Algorithm]

Here, we present an implementation of an adaptive ray-tracing (ART) module in the Athena hydrodynamics code that accurately and efficiently handles the radiative transfer involving multiple point sources on a three-dimensional Cartesian grid. We adopt a recently proposed parallel algorithm that uses nonblocking, asynchronous MPI communications to accelerate transport of rays across the computational domain. We validate our implementation through several standard test problems, including the propagation of radiation in vacuum and the expansions of various types of H II regions. Additionally, scaling tests show that the cost of a full ray trace per source remains comparable to that of the hydrodynamics update on up to $\sim {10}^{3}$ processors. To demonstrate application of our ART implementation, we perform a simulation of star cluster formation in a marginally bound, turbulent cloud, finding that its star formation efficiency is 12% when both radiation pressure forces and photoionization by UV radiation are treated. We directly compare the radiation forces computed from the ART scheme with those from the M 1 closure relation. Although the ART and M 1 schemes yield similar results on large scales, the latter is unable to resolve the radiation field accurately near individual point sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Non-isobaric Thermal Instability

Multiphase media have very complex structure and evolution. Accurate numerical simulations are necessary to make advances in our understanding of this rich physics. Because simulations can capture both the linear and nonlinear evolution of perturbations with a relatively wide range of sizes, it is important to thoroughly understand the stability of condensation and acoustic modes between the two extreme wavelength limits of isobaric and isochoric instability as identified by Field. Partially motivated by a recent suggestion that large non-isobaric clouds can "shatter" into tiny cloudlets, we revisit the linear theory to survey all possible regimes of thermal instability. We uncover seven regimes in total, one of which allows three unstable condensation modes. Using the code Athena++, we determine the numerical requirements to properly evolve small amplitude perturbations of the entropy mode into the nonlinear regime. Here, our 1D numerical simulations demonstrate that for a typical AGN cooling function, the nonlinear evolution of a single eigenmode in an isobarically unstable plasma involves increasingly larger amplitude oscillations in cloud size, temperature, and density as the wavelength increases. Such oscillations are the hallmark behavior of non-isobaric multiphase gas dynamics and may be observable as correlations between changes in brightness and the associated periodic redshifts and blueshifts in systems that can be spatially resolved. Intriguingly, we discuss regimes and derive characteristic cloud sizes for which the saturation process giving rise to these oscillations can be so energetic that the cloud may indeed break apart. However, we dub this process "splattering" instead of "shattering," as it is a different fragmentation mechanism that is triggered when the cloud suddenly "lands" on the stable cold branch of the equilibrium curve.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Catalogue for Astrophysical Turbulence Simulations (CATS)

Turbulence is a key process in many fields of astrophysics. Advances in numerical simulations of fluids over the last several decades have revolutionized our understanding of turbulence and related processes such as star formation and cosmic ray propagation. However, data from numerical simulations of astrophysical turbulence are often not made public. We introduce a new simulation-oriented database for the astronomical community: the Catalogue for Astrophysical Turbulence Simulations (CATS), located at www.mhdturbulence.com. CATS includes magnetohydrodynamic (MHD) turbulent box simulation data products generated by the public codes athena++, arepo, enzo, and flash. CATS also includes several synthetic observational data sets, such as turbulent HI data cubes. We also include measured power spectra and three-point correlation functions from some of these data. We discuss the importance of open-source statistical and visualization tools for the analysis of turbulence simulations such as those found in CATS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

X-Ray Performance of Critical-angle Transmission Grating Prototypes for the Arcus Mission

Arcus is a proposed Explorer Class soft X-ray grating spectrometer. It aims to explore cosmic feedback by mapping hot gases within and between galaxies and galaxy clusters and characterizing jets and winds from supermassive black holes and to investigate the dynamics of protoplanetary disks and stellar accretion. Arcus features 12 m focal-length grazing-incidence silicon pore optics (SPO) developed for the Athena mission. Critical-angle transmission (CAT) gratings efficiently disperse high diffraction orders onto CCDs. We report new and improved X-ray performance results for Arcus-like CAT gratings, including a record resolving power for two coaligned CAT gratings. Multiple Arcus prototype grating facets were illuminated by an SPO at the PANTER facility. The facets consist of 32 × 32.5 mm 2 patterned silicon membranes, bonded to metal frames. The bonding angle is adjusted according to the measured average tilt angle of the grating bars in the membrane. Two simultaneously illuminated facets show a minor broadening of the Al-K α doublet in the 18th and 21st orders with the best-fit record effective resolving power of ${R}_{G}\approx {1.3}_{-0.5}^{+\infty }\times {10}^{4}$ (3σ), about three to four times the Arcus requirement. We measured the diffraction efficiency of quasi-fully illuminated gratings at O-K wavelengths in orders 4–7 in an Arcus-like configuration and compare results with synchrotron spot measurements. After corrections for geometrical effects and bremsstrahlung continuum, we find agreement between full and spot illumination at the two different facilities, as well as with the models used for Arcus effective area predictions. We find that these flight-like gratings meet the diffraction efficiency and greatly exceed the resolving power Arcus requires.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Simulation of Binary-single Interactions in AGN Disk. I. Gas-enhanced Binary Orbital Hardening

Stellar-mass binary black hole (BBH) mergers within the accretion disks of active galactic nuclei may contribute to gravitational wave (GW) events detected by ground-based GW detectors. In particular, the interaction between a BBH and a single stellar-mass black hole (sBH), known as the binary-single interaction (BSI) process, can potentially lead to GW events with detectable nonzero eccentricity. Previous studies of the BSI process, which neglected the effects of gas, showed that BSIs contribute non-negligibly to GW events in a coplanar disk environment. In this work, we conduct a series of two-dimensional hydrodynamical and N-body simulations to explore the BSI in a gas environment by coupling REBOUND with Athena++. We perform 360 simulation runs, spanning parameters in disk surface density Σ 0 and impact parameter b. We find that the gas-induced energy dissipation within the three-body system becomes significant if the encounter velocity between the sBHs is sufficiently large (≫c s ). Our simulation results indicate that approximately half of the end states of the BSI are changed by gas. Furthermore, at higher gas density, the number of close encounters during the BSI process will increase, and the end-state BBHs tend to be more compact. Consequently, the presence of gas may shorten the GW merger timescale for end-state BBHs and increase the three-body merger rate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project

Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope, and high- precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). In this work we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD, and KORAL. Agreement among the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. Thus, we decide that the community of GRMHD codes is mature, capable, and consistent on these test problems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Diagnosing the Kinematics of the Tori in Active Galactic Nuclei with the Velocity-resolved Reverberation Mapping of the Narrow Iron K α Line

The properties of the dusty tori in active galactic nuclei (AGNs) have been investigated in detail, mainly focusing on the geometry and components; however, the kinematics of the torus are still not clear. The narrow iron K α line at 6.4 keV is thought to be produced by the X-ray reflection from the torus. Thus, the velocity-resolved reverberation mapping of it is able to constrain the kinematics of the torus. Such effort is limited by the spectral resolution of current charged coupled device (CCD) detectors and should be possible with the microcalorimeter on the next generation X-ray satellite. In this paper, we first construct the response functions of the torus under a uniform inflow, a Keplerian rotation, and a uniform outflow. Then the energy-dependent light curve of the narrow iron K α line is simulated according to the performance of the X-ray Integral Field Unit in Athena. Finally, the energy-dependent cross-correlation function is calculated to reveal the kinematic signal. According to our results, 100 observations with 5 ks exposure of each are sufficient to distinguish the above three velocity fields. Although the real geometry and velocity field of the torus could be more complex than we assumed, the present result proves the feasibility of the velocity-resolved reverberation mapping of the narrow iron K α line. The combination of the dynamics of the torus with those of the broad-line region and the host galaxy is instructive for the understanding of the feeding and feedback process of AGNs.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1984

CONC/11: A computer program for calculating the performance of dish-type solar thermal collectors and power systems

The CONC/11 computer program designed for calculating the performance of dish-type solar thermal collectors and power systems is discussed. This program is intended to aid the system or collector designer in evaluating the performance to be expected with possible design alternatives. From design or test data on the characteristics of the various subsystems, CONC/11 calculates the efficiencies of the collector and the overall power system as functions of the receiver temperature for a specified insolation. If desired, CONC/11 will also determine the receiver aperture and the receiver temperature that will provide the highest efficiencies at a given insolation. The program handles both simple and compound concentrators. The CONC/11 is written in Athena Extended FORTRAN (similar to FORTRAN 77) to operate primarily in an interactive mode on a Sperry 1100/81 computer. It could also be used on many small computers. A user's manual is also provided for this program.

Jaffe, L. D.↗
NASA NTRS1987

Evaluating Performances of Solar-Energy Systems

CONC11 computer program calculates performances of dish-type solar thermal collectors and power systems. Solar thermal power system consists of one or more collectors, power-conversion subsystems, and powerprocessing subsystems. CONC11 intended to aid system designer in comparing performance of various design alternatives. Written in Athena FORTRAN and Assembler.

Jaffe, L. D.↗
NASA NTRS1991

Analyzing Large Reflector Antenna Structures

Iterative Design of Antenna Structures (IDEAS) computer program performs finite-element analysis and optimization of design, with special features for analysis and design of microwave antennas and associated structural components. Primary feature to perform root-mean-square-best-fitting path-length-error and boresight-pointing-error analyses for Cassegrain microwave antenna structure loaded by gravity, wind, or other environmental effects. Flexible and easy to use. Written in ATHENA FORTRAN and ASSEMBLER.

Levy, Roy↗
NASA NTRS1990

Continuation of research into software for space operations support, volume 1

A prototype workstation executive called the Hardware Independent Software Development Environment (HISDE) was developed. Software technologies relevant to workstation executives were researched and evaluated and HISDE was used as a test bed for prototyping efforts. New X Windows software concepts and technology were introduced into workstation executives and related applications. The four research efforts performed included: (1) Research into the usability and efficiency of Motif (an X Windows based graphic user interface) which consisted of converting the existing Athena widget based HISDE user interface to Motif demonstrating the usability of Motif and providing insight into the level of effort required to translate an application from widget to another; (2) Prototype a real time data display widget which consisted of research methods for and prototyping the selected method of displaying textual values in an efficient manner; (3) X Windows performance evaluation which consisted of a series of performance measurements which demonstrated the ability of low level X Windows to display textural information; (4) Convert the Display Manager to X Window/Motif which is the application used by NASA for data display during operational mode.

Collier, Mark D.↗
NASA NTRS1992

Update of GRASP/Ada reverse engineering tools for Ada

The GRASP/Ada project (Graphical Representations of Algorithms, Structures, and Processes for Ada) has successfully created and prototyped a new algorithmic level graphical representation of Ada software, the Control Structure Diagram (CSD). The primary impetus for creation of the CSD was to improve the comprehension efficiency of Ada software and, as a result, improve reliability and reduce costs. The emphasis was on the automatic generation of the CSD from Ada PDL or source code to support reverse engineering and maintenance. The CSD has the potential to replace traditional prettyprinted Ada source code. In Phase 1 of the GRASP/Ada project, the CSD graphical constructs were created and applied manually to several small Ada programs. A prototype (Version 1) was designed and implemented using FLEX and BISON running under VMS on a VAS 11-780. In Phase 2, the prototype was improved and ported to the Sun 4 platform under UNIX. A user interface was designed and partially implemented using the HP widget toolkit and the X Windows System. In Phase 3, the user interface was extensively reworked using the Athena widget toolkit and X Windows. The prototype was applied successfully to numerous Ada programs ranging in size from several hundred to several thousand lines of source code. Following Phase 3, the prototype was evaluated by software engineering students at Auburn University and then updated with significant enhancements to the user interface including editing capabilities. Version 3.2 of the prototype was prepared for limited distribution to facilitate further evaluation. The current prototype provides the capability for the user to generate CSD's from Ada PDL or source code in a reverse engineering as well as forward engineering mode with a level of flexibility suitable for practical application.

Cross, James H., II↗