Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 505 records · Page 28

FORMULATION OF NON-STEADY-STATE DUST FORMATION PROCESS IN ASTROPHYSICAL ENVIRONMENTS

The non-steady-state formation of small clusters and the growth of grains accompanied by chemical reactions are formulated under the consideration that the collision of key gas species (key molecule) controls the kinetics of dust formation process. The formula allows us to evaluate the size distribution and condensation efficiency of dust formed in astrophysical environments. We apply the formulation to the formation of C and MgSiO{sub 3} grains in the ejecta of supernovae, as an example, to investigate how the non-steady effect influences the formation process, condensation efficiency f{sub con,{sub ∞}}, and average radius a{sub ave,{sub ∞}} of newly formed grains in comparison with the results calculated with the steady-state nucleation rate. We show that the steady-state nucleation rate is a good approximation if the collision timescale of key molecule τ{sub coll} is much smaller than the timescale τ{sub sat} with which the supersaturation ratio increases; otherwise the effect of the non-steady state becomes remarkable, leading to a lower f{sub con,{sub ∞}} and a larger a{sub ave,{sub ∞}}. Examining the results of calculations, we reveal that the steady-state nucleation rate is applicable if the cooling gas satisfies Λ ≡ τ{sub sat}/τ{sub coll} ∼> 30 during the formation of dust, and find that f{sub con,{sub ∞}} and a{sub ave,{sub ∞}} are uniquely determined by Λ{sub on} at the onset time t{sub on} of dust formation. The approximation formulae for f{sub con,{sub ∞}} and a{sub ave,{sub ∞}} as a function of Λ{sub on} could be useful in estimating the mass and typical size of newly formed grains from observed or model-predicted physical properties not only in supernova ejecta but also in mass-loss winds from evolved stars.

79 ASTRONOMY AND ASTROPHYSICS↗

MEASUREMENTS OF CYCLOTRON FEATURES AND PULSE PERIODS IN THE HIGH-MASS X-RAY BINARIES 4U 1538–522 AND 4U 1907+09 WITH THE INTERNATIONAL GAMMA-RAY ASTROPHYSICS LABORATORY

We present a spectral and timing analysis of International Gamma-Ray Astrophysics Laboratory (INTEGRAL) observations of two high-mass X-ray binaries, 4U 1538–522 and 4U 1907+09. Our timing measurements for 4U 1538–522 find the pulse period to have exhibited a spin-up trend until approximately 2009, after which there is evidence for a torque reversal, with the source beginning to spin down to the most recently measured period of 525.407 ± 0.001 s. The most recent INTEGRAL observations of 4U 1907+09 are not found to yield statistically significant pulse periods due to the significantly lower flux from the source compared with 4U 1538–522. A spectral model consisting of a power-law continuum with an exponential cutoff and modified by two cyclotron resonance scattering features is found to fit both sources well, with the cyclotron scattering features detected at ∼22 and ∼49 keV for 4U 1538–522 and at ∼18 and ∼36 keV for 4U 1907+09. The spectral parameters of 4U 1538–522 are generally not found to vary significantly with flux and there is little to no variation across the torque reversal. Examining our results in conjunction with previous work, we find no evidence for a correlation between cyclotron line energy and luminosity for 4U 1538–522. 4U 1907+09 shows evidence for a positive correlation between cyclotron line energy and luminosity, which would make it the fourth, and lowest luminosity, cyclotron line source to exhibit this relationship.

79 ASTRONOMY AND ASTROPHYSICS↗

Astrophysical tests of modified gravity: Constraints from distance indicators in the nearby universe

We use distance measurements in the nearby universe to carry out new tests of gravity, surpassing other astrophysical tests by over two orders of magnitude for chameleon theories. The three nearby distance indicators—cepheids, tip of the red giant branch (TRGB) stars, and water masers—operate in gravitational fields of widely different strengths. This enables tests of scalar-tensor gravity theories because they are screened from enhanced forces to different extents. Inferred distances from cepheids and TRGB stars are altered (in opposite directions) over a range of chameleon gravity theory parameters well below the sensitivity of cosmological probes. Using published data, we have compared cepheid and TRGB distances in a sample of unscreened dwarf galaxies within 10 Mpc. We use a comparable set of screened galaxies as a control sample. We find no evidence for the order unity force enhancements expected in these theories. Using a two-parameter description of the models (the coupling strength and background field value), we obtain constraints on both the chameleon and symmetron screening scenarios. In particular we show that f(R) models with background field values f {sub R0} above 5 × 10{sup –7} are ruled out at the 95% confidence level. We also compare TRGB and maser distances to the galaxy NGC 4258 as a second test for larger field values. While there are several approximations and caveats in our study, our analysis demonstrates the power of gravity tests in the local universe. We discuss the prospects for additional improved tests with future observations.

79 ASTRONOMY AND ASTROPHYSICS↗

New determination of the {sup 2}H(d,p){sup 3}H and {sup 2}H(d,n){sup 3}He reaction rates at astrophysical energies

The cross sections of the {sup 2}H(d,p){sup 3}H and {sup 2}H(d,n){sup 3}He reactions have been measured via the Trojan Horse method applied to the quasi-free {sup 2}H({sup 3}He,p {sup 3}H){sup 1}H and {sup 2}H({sup 3}He,n {sup 3}He){sup 1}H processes at 18 MeV off the proton in {sup 3}He. For the first time, the bare nucleus S(E) factors have been determined from 1.5 MeV, across the relevant region for standard Big Bang nucleosynthesis, down to the thermal energies of deuterium burning in the pre-main-sequence (PMS) phase of stellar evolution, as well as of future fusion reactors. Both the energy dependence and the absolute value of the S(E) factors deviate by more than 15% from the available direct data and existing fitting curves, with substantial variations in the electron screening by more than 50%. As a consequence, the reaction rates for astrophysics experience relevant changes, with a maximum increase of up to 20% at the temperatures of the PMS phase. From a recent primordial abundance sensitivity study, it turns out that the {sup 2}H(d,n){sup 3}He reaction is quite influential on {sup 7}Li, and the present change in the reaction rate leads to a decrease in its abundance by up to 10%. The present reaction rates have also been included in an updated version of the FRANEC evolutionary code to analyze their influence on the central deuterium abundance in PMS stars with different masses. The largest variation of about 10%-15% pertains to young stars (≤1 Myr) with masses ≥1 M {sub ☉}.

79 ASTRONOMY AND ASTROPHYSICS↗

Theoretical hot methane line lists up to T = 2000 K for astrophysical applications

The paper describes the construction of complete sets of hot methane lines based on accurate ab initio potential and dipole moment surfaces and extensive first-principle calculations. Four line lists spanning the [0-5000] cm{sup –1} infrared region were built at T = 500, 1000, 1500, and 2000 K. For each of these four temperatures, we have constructed two versions of line lists: a version for high-resolution applications containing strong and medium lines and a full version appropriate for low-resolution opacity calculations. A comparison with available empirical databases is discussed in detail for both cold and hot bands giving a very good agreement for line positions, typically <0.1-0.5 cm{sup –1} and ∼5% for intensities of strong lines. Together with numerical tests using various basis sets, this confirms the computational convergence of our results for the most important lines, which is the major issue for theoretical spectra predictions. We showed that transitions with lower state energies up to 14,000 cm{sup –1} could give significant contributions to the methane opacity and have to be systematically taken into account. Our list at 2000 K calculated up to J = 50 contains 11.5 billion transitions for I > 10{sup –29} cm mol{sup –1}. These new lists are expected to be quantitatively accurate with respect to the precision of available and currently planned observations of astrophysical objects with improved spectral resolution.

79 ASTRONOMY AND ASTROPHYSICS↗

The photochemistry of pyrimidine in realistic astrophysical ices and the production of nucleobases

Nucleobases, together with deoxyribose/ribose and phosphoric acid, are the building blocks of DNA and RNA for all known life. The presence of nucleobase-like compounds in carbonaceous chondrites delivered to the Earth raises the question of an extraterrestrial origin for the molecules that triggered life on our planet. Whether these molecules are formed in interstellar/protostellar environments, in small parent bodies in the solar system, or both, is currently unclear. Recent experiments show that the UV irradiation of pyrimidine (C{sub 4}H{sub 4}N{sub 2}) in H{sub 2}O-rich ice mixtures that contain NH{sub 3}, CH{sub 3}OH, or CH{sub 4} leads to the formation of the pyrimidine-based nucleobases uracil, cytosine, and thymine. In this work, we discuss the low-temperature UV irradiation of pyrimidine in realistic astrophysical ice mixtures containing H{sub 2}O, CH{sub 3}OH, and NH{sub 3}, with or without CH{sub 4}, to search for the production of nucleobases and other prebiotic compounds. These experiments show the presence of uracil, urea, glycerol, hexamethylenetetramine, small amino acids, and small carboxylic acids in all samples. Cytosine was only found in one sample produced from ices irradiated with a higher UV dose, while thymine was not found in any sample, even after irradiation with a higher UV dose. Results are discussed to evaluate the role of the photochemistry of pyrimidine in the inventory of organic molecules detected in meteorites and their astrophysical/astrobiological implications.

79 ASTRONOMY AND ASTROPHYSICS↗

Scaling of magneto-quantum-radiative hydrodynamic equations: from laser-produced plasmas to astrophysics

We introduce the equations of magneto-quantum-radiative hydrodynamics. By rewriting them in a dimensionless form, we obtain a set of parameters that describe scale-dependent ratios of characteristic hydrodynamic quantities. We discuss how these dimensionless parameters relate to the scaling between astrophysical observations and laboratory experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Atomic chemistry in turbulent astrophysical media. I. EFFECT OF Atomic cooling

We carry out direct numerical simulations of turbulent astrophysical media that explicitly track ionizations, recombinations, and species-by-species radiative cooling. The simulations assume solar composition and follows the evolution of hydrogen, helium, carbon, oxygen, sodium, and magnesium, but they do not include the presence of an ionizing background. In this case, the medium reaches a global steady state that is purely a function of the one-dimensional turbulent velocity dispersion, σ{sub 1D}, and the product of the mean density and the driving scale of turbulence, nL. Our simulations span a grid of models with σ{sub 1D} ranging from 6 to 58 km s{sup −1} and nL ranging from 10{sup 16} to 10{sup 20} cm{sup −2}, which correspond to turbulent Mach numbers from M = 0.2 to 10.6. The species abundances are well described by single-temperature estimates whenever M is small, but local equilibrium models can not accurately predict the global equilibrium abundances when M≳1. To allow future studies to account for nonequilibrium effects in turbulent media, we gather our results into a series of tables, which we will extend in the future to encompass a wider range of elements, compositions, and ionizing processes.

79 ASTRONOMY AND ASTROPHYSICS↗

ON THE ANISOTROPIC NATURE OF MRI-DRIVEN TURBULENCE IN ASTROPHYSICAL DISKS

The magnetorotational instability (MRI) is thought to play an important role in enabling accretion in sufficiently ionized astrophysical disks. The rate at which MRI-driven turbulence transports angular momentum is intimately related to both the strength of the amplitudes of the fluctuations on various scales and the degree of anisotropy of the underlying turbulence. This has motivated several studies to characterize the distribution of turbulent power in spectral space. In this paper we investigate the anisotropic nature of MRI-driven turbulence using a pseudo-spectral code and introduce novel ways for providing a robust characterization of the underlying turbulence. We study the growth of the MRI and the subsequent transition to turbulence via parasitic instabilities, identifying their potential signature in the late linear stage. We show that the general flow properties vary in a quasi-periodic way on timescales comparable to ∼10 inverse angular frequencies, motivating the temporal analysis of its anisotropy. We introduce a 3D tensor invariant analysis to quantify and classify the evolution of the anisotropy of the turbulent flow. This analysis shows a continuous high level of anisotropy, with brief sporadic transitions toward two- and three-component isotropic turbulent flow. This temporal-dependent anisotropy renders standard shell averaging especially when used simultaneously with long temporal averages, inadequate for characterizing MRI-driven turbulence. We propose an alternative way to extract spectral information from the turbulent magnetized flow, whose anisotropic character depends strongly on time. This consists of stacking 1D Fourier spectra along three orthogonal directions that exhibit maximum anisotropy in Fourier space. The resulting averaged spectra show that the power along each of the three independent directions differs by several orders of magnitude over most scales, except the largest ones. Our results suggest that a first-principles theory to describe fully developed MRI-driven turbulence will likely have to consider the anisotropic nature of the flow at a fundamental level.

79 ASTRONOMY AND ASTROPHYSICS↗

A HYDROCHEMICAL HYBRID CODE FOR ASTROPHYSICAL PROBLEMS. I. CODE VERIFICATION AND BENCHMARKS FOR A PHOTON-DOMINATED REGION (PDR)

A two-dimensional hydrochemical hybrid code, KM2, is constructed to deal with astrophysical problems that would require coupled hydrodynamical and chemical evolution. The code assumes axisymmetry in a cylindrical coordinate system and consists of two modules: a hydrodynamics module and a chemistry module. The hydrodynamics module solves hydrodynamics using a Godunov-type finite volume scheme and treats included chemical species as passively advected scalars. The chemistry module implicitly solves nonequilibrium chemistry and change of energy due to thermal processes with transfer of external ultraviolet radiation. Self-shielding effects on photodissociation of CO and H{sub 2} are included. In this introductory paper, the adopted numerical method is presented, along with code verifications using the hydrodynamics module and a benchmark on the chemistry module with reactions specific to a photon-dominated region (PDR). Finally, as an example of the expected capability, the hydrochemical evolution of a PDR is presented based on the PDR benchmark.

79 ASTRONOMY AND ASTROPHYSICS↗

SHOCK STRENGTH IN SPACE AND ASTROPHYSICAL PLASMAS

The Rankine–Hugoniot (R–H) jump conditions are the most important and frequently used equations in studies of the properties of space and astrophysical plasmas during their passage through shock discontinuities. This paper revisits the R–H conditions for shocks, develops the formulation of the compression ratio, and examines its range of values and properties. The analysis expresses the downstream thermodynamic variables and the compression ratio as functions of the upstream thermodynamic variables, either for equal or different polytropic indices upstream and downstream of the shock. In the general case of space plasmas with an oblique magnetic field, the compression ratio is given by a quartic polynomial, which is reduced to a cubic trinomial when the upstream/downstream polytropic indices are equal. The special cases of magnetic fields that are perpendicular or parallel to the shock normal are also examined. In any case, the compression ratio polynomial has one degree larger order, when the upstream/downstream polytropic indices are different. Emphasis is placed on the maximum value of the compression ratio, which is known to be ∼4 for adiabatic polytropic index ∼5/3. However, the compression ratio can be much larger if the upstream/downstream polytropic indices are not equal to each other and less than one. Several other issues are investigated: (i) the entropic condition, showing that statistical mechanics and thermodynamics lead to the same relation of entropy variation; (ii) the effect of kappa distributions on jump conditions; and (iii) the upper limit of the upstream temperature for a shock to exist.

79 ASTRONOMY AND ASTROPHYSICS↗

Vector interaction enhanced bag model for astrophysical applications

For quark matter studies in astrophysics the thermodynamic bag model (tdBAG) has been widely used. Despite its success it fails to account for various phenomena expected from QCD. We suggest a straightforward extension of tdBAG in order to take the dynamical breaking of chiral symmetry and the influence of vector interactions explicitly into account. As for tdBAG the model mimics confinement in a phenomenological approach. It is based on an analysis of the Nambu–Jona-Lasinio (NJL) model at finite density. Furthermore, we demonstrate how NJL and bag models in this regime follow from the more general and QCD-based framework of the Dyson–Schwinger equations in a medium by assuming simple gluon contact interaction. Based on our simple and novel model, we construct quark hadron hybrid equations of state and systematically study chiral and deconfinement phase transitions, the appearance of s-quarks, and the role of vector interaction. We further study these aspects for matter in β-equilibrium at zero temperature, with particular focus on the current ∼2 M{sub ⊙} maximum mass constraint for neutron stars. Our approach indicates that the currently only theoretical evidence for the hypothesis of stable strange matter is an artifact of tdBAG and results from neglecting the dynamical breaking of chiral symmetry.

79 ASTRONOMY AND ASTROPHYSICS↗

SHOCK WAVE STRUCTURE IN ASTROPHYSICAL FLOWS WITH AN ACCOUNT OF PHOTON TRANSFER

For an accurate treatment of the shock wave propagation in high-energy astrophysical phenomena, such as supernova shock breakouts, gamma-ray bursts and accretion disks, knowledge of radiative transfer plays a crucial role. In this paper we consider one-dimensional (1D) special relativistic radiation hydrodynamics by solving the Boltzmann equation for radiative transfer. The structure of a radiative shock is calculated for a number of shock tube problems, including strong shock waves, and relativistic- and radiation-dominated cases. Calculations are performed using an iterative technique that consistently solves the equations of relativistic hydrodynamics and relativistic comoving radiative transfer. A comparison of radiative transfer solutions with the Eddington approximation and the M1 closure is made. A qualitative analysis of moment equations for radiation is performed and the conditions for the existence of jump discontinuity for non-relativistic cases are investigated numerically.

79 ASTRONOMY AND ASTROPHYSICS↗

Thresholded power law size distributions of instabilities in astrophysics

Power-law-like size distributions are ubiquitous in astrophysical instabilities. There are at least four natural effects that cause deviations from ideal power law size distributions, which we model here in a generalized way: (1) a physical threshold of an instability; (2) incomplete sampling of the smallest events below a threshold x{sub 0}; (3) contamination by an event-unrelated background x{sub b}; and (4) truncation effects at the largest events due to a finite system size. These effects can be modeled in the simplest terms with a “thresholded power law” distribution function (also called generalized Pareto [type II] or Lomax distribution), N(x)dx∝(x+x{sub 0}){sup −a}dx, where x{sub 0} > 0 is positive for a threshold effect, while x{sub 0} < 0 is negative for background contamination. We analytically derive the functional shape of this thresholded power law distribution function from an exponential growth evolution model, which produces avalanches only when a disturbance exceeds a critical threshold x{sub 0}. We apply the thresholded power law distribution function to terrestrial, solar (HXRBS, BATSE, RHESSI), and stellar flare (Kepler) data sets. We find that the thresholded power law model provides an adequate fit to most of the observed data. Major advantages of this model are the automated choice of the power law fitting range, diagnostics of background contamination, physical instability thresholds, instrumental detection thresholds, and finite system size limits. When testing self-organized criticality models that predict ideal power laws, we suggest including these natural truncation effects.

79 ASTRONOMY AND ASTROPHYSICS↗

BOW SHOCK FRAGMENTATION DRIVEN BY A THERMAL INSTABILITY IN LABORATORY ASTROPHYSICS EXPERIMENTS

The role of radiative cooling during the evolution of a bow shock was studied in laboratory-astrophysics experiments that are scalable to bow shocks present in jets from young stellar objects. The laboratory bow shock is formed during the collision of two counterstreaming, supersonic plasma jets produced by an opposing pair of radial foil Z-pinches driven by the current pulse from the MAGPIE pulsed-power generator. The jets have different flow velocities in the laboratory frame, and the experiments are driven over many times the characteristic cooling timescale. The initially smooth bow shock rapidly develops small-scale nonuniformities over temporal and spatial scales that are consistent with a thermal instability triggered by strong radiative cooling in the shock. The growth of these perturbations eventually results in a global fragmentation of the bow shock front. The formation of a thermal instability is supported by analysis of the plasma cooling function calculated for the experimental conditions with the radiative packages ABAKO/RAPCAL.

79 ASTRONOMY AND ASTROPHYSICS↗

ELECTRON IMPACT DISSOCIATION X {sup 1}{Sigma}{sup +} {sub g} {yields} b {sup 3}{Sigma} {sub u} {sup +} AND EXCITATIONS X {sup 1}{Sigma}{sup +} {sub g} {yields} a {sup 3}{Sigma} {sub g} {sup +} AND X {sup 1}{Sigma}{sup +} {sub g} {yields} B {sup 1}{Sigma} {sub u} {sup +} OF MOLECULAR HYDROGEN IN NONTHERMAL ASTROPHYSICAL PLASMAS

We investigate the electronic transitions X {sup 1}{Sigma}{sup +} {sub g} {yields} b {sup 3}{Sigma} {sub u} {sup +}, X {sup 1}{Sigma}{sup +} {sub g} {yields} a {sup 3}{Sigma} {sub g} {sup +}, and X {sup 1}{Sigma}{sup +} {sub g} {yields} B {sup 1}{Sigma} {sub u} {sup +} of molecular hydrogen by studying electron impacts in astrophysical Lorentzian plasmas. Useful fitting formulae for the X {sup 1}{Sigma}{sup +} {sub g} {yields} b {sup 3}{Sigma} {sub u} {sup +}, X {sup 1}{Sigma}{sup +} {sub g} {yields} a {sup 3}{Sigma} {sub g} {sup +}, and X {sup 1}{Sigma}{sup +} {sub g} {yields} B {sup 1}{Sigma} {sub u} {sup +} excitation cross sections are employed in order to obtain the electronic excitation rate coefficients of H{sub 2} as functions of the spectral index and temperature. In low-temperature regions, it is found that the excitation rate coefficients R{sub b{sup 3}{Sigma}{sub u{sup {sub +}}}}, R{sub a{sup 3}{Sigma}{sub g{sup {sub +}}}}, and R{sub B{sub {sup 1}{Sigma}{sub u{sup {sub +}}}}} of H{sub 2} in non-Maxwellian plasmas are smaller than those in Maxwellian plasmas. However, in high-temperature regions, the excitation rate coefficients of H{sub 2} in non-Maxwellian plasmas are greater than those in Maxwellian plasmas. It is also shown that the X {sup 1}{Sigma}{sup +} {sub g} {yields} b {sup 3}{Sigma} {sub u} {sup +} excitation rate coefficient is the main contributor in low-temperature regions. In contrast, it is found that the X {sup 1}{Sigma}{sup +} {sub g} {yields} B {sup 1}{Sigma} {sub u} {sup +} electronic excitation is dominant in high-temperature regions.

79 ASTRONOMY AND ASTROPHYSICS↗

STARLIB: A NEXT-GENERATION REACTION-RATE LIBRARY FOR NUCLEAR ASTROPHYSICS

STARLIB is a next-generation, all-purpose nuclear reaction-rate library. For the first time, this library provides the rate probability density at all temperature grid points for convenient implementation in models of stellar phenomena. The recommended rate and its associated uncertainties are also included. Currently, uncertainties are absent from all other rate libraries, and, although estimates have been attempted in previous evaluations and compilations, these are generally not based on rigorous statistical definitions. A common standard for deriving uncertainties is clearly warranted. STARLIB represents a first step in addressing this deficiency by providing a tabular, up-to-date database that supplies not only the rate and its uncertainty but also its distribution. Because a majority of rates are lognormally distributed, this allows the construction of rate probability densities from the columns of STARLIB. This structure is based on a recently suggested Monte Carlo method to calculate reaction rates, where uncertainties are rigorously defined. In STARLIB, experimental rates are supplemented with: (1) theoretical TALYS rates for reactions for which no experimental input is available, and (2) laboratory and theoretical weak rates. STARLIB includes all types of reactions of astrophysical interest to Z = 83, such as (p, {gamma}), (p, {alpha}), ({alpha}, n), and corresponding reverse rates. Strong rates account for thermal target excitations. Here, we summarize our Monte Carlo formalism, introduce the library, compare methods of correcting rates for stellar environments, and discuss how to implement our library in Monte Carlo nucleosynthesis studies. We also present a method for accessing STARLIB on the Internet and outline updated Monte Carlo-based rates.

79 ASTRONOMY AND ASTROPHYSICS↗

MPI-AMRVAC FOR SOLAR AND ASTROPHYSICS

In this paper, we present an update to the open source MPI-AMRVAC simulation toolkit where we focus on solar and non-relativistic astrophysical magnetofluid dynamics. We highlight recent developments in terms of physics modules, such as hydrodynamics with dust coupling and the conservative implementation of Hall magnetohydrodynamics. A simple conservative high-order finite difference scheme that works in combination with all available physics modules is introduced and demonstrated with the example of monotonicity-preserving fifth-order reconstruction. Strong stability-preserving high-order Runge-Kutta time steppers are used to obtain stable evolutions in multi-dimensional applications, realizing up to fourth-order accuracy in space and time. With the new distinction between active and passive grid cells, MPI-AMRVAC is ideally suited to simulate evolutions where parts of the solution are controlled analytically or have a tendency to progress into or out of a stationary state. Typical test problems and representative applications are discussed with an outlook toward follow-up research. Finally, we discuss the parallel scaling of the code and demonstrate excellent weak scaling up to 30, 000 processors, allowing us to exploit modern peta-scale infrastructure.

79 ASTRONOMY AND ASTROPHYSICS↗