Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accretion”

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 883 records · Page 49

Broadband X-Ray Constraints on the Accreting Black Hole in Quasar 4C 74.26

X-ray data for quasar 4C 74.26 have previously been modeled with a broad Fe Kα emission line and reflection continuum originating in the inner part of the accretion disk around the central supermassive black hole (SMBH), i.e., the strong-gravity regime. We modeled broadband X-ray spectra from Suzaku and NuSTAR with MYTORUS, self-consistently accounting for Fe Kα line emission, as well as direct and reflected continuum emission, from matter with a finite column density. A narrow Fe Kα emission line originating in an X-ray reprocessor with solar Fe abundance far from the central SMBH is sufficient to produce excellent fits for all spectra. For the first time, we are able to measure the global column density, out of the line of sight, to be in the range ∼1.5 to ∼2.9 ´ 10(exp 24) per sq.cm, i.e., in the Compton-thick regime, while the line-of-sight column density is Compton-thin in all observations. The Fe Kα emission line is unresolved in all observations but one. The Compton-scattered continuum from distant matter removes the need for relativistic broadening of the Fe Kα emission line, which is required for SMBH spin measurements. The resolved line observation can alternatively be modeled with a relativistic model but we do not find evidence for a truncated accretion disk model. We conclude that the X-ray emission in these 4C 74.26 data is unlikely to originate in the inner accretion disk region and thus cannot be used to measure SMBH spin.

Tzanavaris, P.↗

Three Dimensional Surface Redefinition Method for Computational Ice Accretion Solvers

Computational tools have been increasing in maturity and are thus commonly used in the engineering design process. Advancements in NASA’s current state of the art computational ice accretion tool are required to tackle the icing challenges of tomorrow. GlennICE, a next generation ice accretion solver, is under development at NASA to tackle these challenges. One of the elements of this migration is transitioning from a quasi-three dimensional strip theory based ice accretion methodology to a fully three dimensional methodology. This requires construction of a method to redefine or extrude a discretized or tessellated surface geometry based on the predicted volumetric ice growth for each tessellated surface triangle. This paper describes the methodology that is employed in the GlennICE software, and assesses the performance of the method in replicating a well defined analytical test case.

Computational↗

Effect of Geometric Fidelity on the Aerodynamics of a Swept Wing with Glaze Ice Accretion

Aerodynamic assessment of icing effects on swept wings is an important component of a larger effort to improve three-dimensional icing simulation capabilities. An understanding of ice-shape geometric fidelity on iced-wing aerodynamics and the associated flowfield features are needed to guide the development and validation of ice-accretion simulation tools. To this end, wind-tunnel testing was carried out for 8.9% and 13.3% scale semispan wing models based upon the Common Research Model airplane configuration. Various levels of geometric fidelity of an artificial ice shape representing a glaze-ice accretion on a swept wing were investigated. The highest fidelity artificial ice shape reproduced all of the three-dimensional features associated with the glaze ice accretion. The lowest fidelity artificial ice shapes were simple, spanwise-varying horn ice geometries intended to represent the maximum ice thickness on the wing upper surface. The results presented in this paper show that the addition of grit roughness to some lower-fidelity artificial ice shapes resulted in favorable lift and pitching moment comparisons to the wing with the highest fidelity artificial ice shape. In the range of 4.3 to 7.4 deg. angle of attack, surface oil flow visualization and pressure data show that the wing with the two lower fidelity simulations clearly demonstrated a leading edge vortex dominated flowfield, referred to as type I. For the wing with the high fidelity ice shape, the flowfield at lower angles of attack was characterized by streamwise-running, counter-rotating vortical flow referred to as type II. Between 6.4 and 7.4 deg. angle of attack, the effect of the type II flow structures was significantly altered and gave way to the type I leading edge vortex. This means that for angles of attack 7.4 deg. and higher, the wing with all three configurations exhibited the same type of flowfield. This helps to explain why there is reasonably good agreement in the lift and pitching moment coefficients among these configurations.

Aircraft icing↗

Three Dimensional Surface Redefinition Method for Computational Ice Accretion Solvers

Computational tools have been increasing in maturity and are thus commonly used in the engineering design process. Advancements in NASA’s current state of the art computational ice accretion tool are required to tackle the icing challenges of tomorrow. GlennICE, a next generation ice accretion solver, is under development at NASA to tackle these challenges. One of the elements of this migration is transitioning from a quasi-three dimensional strip theory based ice accretion methodology to a fully three dimensional methodology. This requires construction of a method to redefine or extrude a discretized or tessellated surface geometry based on the predicted volumetric ice growth for each tessellated surface triangle. This paper describes the methodology that is employed in the GlennICE software, and assesses the performance of the method in replicating a well definied analytical test case.

Computational↗

Timing the Pulsations of The Accreting Millisecond Pulsar SAX J1808.4-3658 During Its 2019 Outburst

In this paper we present a coherent timing analysis of the 401 Hz pulsations of the accreting millisecond X-ray pulsar SAX J1808.4-3658 during its 2019 outburst. Using observations collected with the Neutron Star Interior Composition Explorer (NICER), we establish the pulsar spin frequency and orbital phase during its latest epoch. We find that the 2019 outburst shows a pronounced evolution in pulse phase over the course of the outburst. These phase shifts are found to correlate with the source flux, and are interpreted in terms of hot-spot drift on the stellar surface, driven by changes in the mass accretion rate. Additionally, we find that the long-term evolution of the pulsar spin frequency shows evidence for a modulation at the Earth's orbital period, enabling pulsar timing based astrometry of this accreting millisecond pulsar.

Peter Bult↗

NuSTAR Observation of GROJ1744−28 at Low Mass Accretion Rate

Context.Neutron stars in low-mass X-ray binaries (LMXBs) are important systems to study the physics of accretion onto compact ob-jects. The system GRO J1744−28 is particularly interesting, as it usually shows clear pulsations as well as X-ray bursts. Additionally,there are claims for a magnetic field of 5×1011G through the detection of a Cyclotron Resonant Scattering Feature (CRSF).Aims.We present the spectral analysis of GRO J1744−28 using∼29 ks of NuSTAR data taken in 2017 February at a low luminosityof 3.2×1036erg s−1(3–50 keV). Our goal is to study the variability of the source spectrum with pulse phase and to search for the claimed CRSF. Methods.The continuum spectrum is modeled with an absorbed power-law with exponential cut-off, and an additional iron line component. We find no obvious indications for a CRSF and therefore perform a detailed cyclotron line search using statistical methods.We perform this search on pulse phase-averaged as well as phase-resolved spectra. Results.GRO J1744−28 was observed at one of the lowest luminosities studied so far. The previously detected Type II X-ray bursts re absent. Clear pulsations at a period of 2.141124(9) Hz are detected. The pulse profile shows an indication of a secondary peak,which was not seen at higher flux. The 4σupper limit for the strength of a CRSF in the 3–20 keV band is 0.07 keV, lower than the strength of the line found at higher luminosity. Conclusions.The detection of pulsations shows that the source did not enter the “propeller” regime, even though the source flux of4.15×10−10erg cm−2s−1was almost one order of magnitude below the threshold for the propeller regime claimed in previous studies on this source. The transition into the propeller regime in GRO J1744−28 must therefore be below a luminosity of 3.2×1036erg s−1(3–50 keV), which implies a surface magnetic field.2.9×1011G and mass accretion rate.1.7×1016g s−1. A change of the CRSF depth as function of luminosity is not unexpected and has been observed in other sources. This result possibly implies a change in emission geometry as function of mass accretion rate to reduce the depth of the line below our detection limit.

Ole König↗

Fitting strategies of accretion column models and application to the broadband spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT dataset of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies.

X-rays: binaries↗

The First Swift Intensive AGN Accretion Disk Reverberation Mapping Survey

Swift intensive accretion disk reverberation mapping of four AGN yielded light curves sampled ∼200–350 times in 0.3–10 keV X-ray and six UV/optical bands. Uniform reduction and cross-correlation analysis of these data sets yields three main results: (1) The X-ray/UV correlations are much weaker than those within the UV/optical, posing severe problems for the lamp-post reprocessing model in which variations in a central X-ray corona drive and power those in the surrounding accretion disk. (2) The UV/optical interband lags are generally consistent with t μ l4 3 as predicted by the centrally illuminated thin accretion disk model. While the average interband lags are somewhat larger than predicted, these results alone are not inconsistent with the thin disk model given the large systematic uncertainties involved. (3) The one exception is the U band lags, which are on average a factor of ∼2.2 larger than predicted from the surrounding band data and fits. This excess appears to be due to diffuse continuum emission from the broad-line region (BLR). The precise mixing of disk and BLR components cannot be determined from these data alone. The lags in different AGN appear to scale with mass or luminosity. We also find that there are systematic differences between the uncertainties derived by JAVELIN versus more standard lag measurement techniques, with JAVELIN reporting smaller uncertainties by a factor of 2.5 on average. In order to be conservative only standard techniques were used in the analyses reported herein.

R. Edelson↗

Mini-disk accretion onto spinning black hole binaries: quasi-periodicities and outflows

We perform a full 3D general relativistic magnetohydrodynamical (GRMHD) simulation of an equal-mass, spinning, binary black hole approaching merger, surrounded by a circumbinary disk and with mini-disks around each black hole. For this purpose, we evolve the ideal GRMHD equations on top of an approximated spacetime for the binary that is valid in every position of space, including the black hole horizons, during the inspiral regime. We use relaxed initial data for the circumbinary disk from a previous long-term simulation, where the accretion is dominated by an m = 1 overdensity called the lump. We compare our new spinning simulation with a previous non-spinning run, studying how spin influences the mini-disk properties. We analyze the accretion from the inner edge of the lump to the black hole, focusing on the angular momentum budget of the fluid around the mini-disks. We find that mini-disks in the spinning case have more mass over a cycle than the non-spinning case. However, in both cases we find most of the mass received by the black holes is delivered by direct plunging of material from the lump. We also analyze the morphology and variability of the electromagnetic fluxes and we find they share the same periodicities of the accretion rate. In the spinning case, we find that the outflows are stronger than the non-spinning case. Our results will be useful to understand and produce realistic synthetic light curves and spectra, which can be used in future observations.

Astrophysics↗

Fitting Strategies of Accretion Column Models and Application to the Broadband Spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT data-set of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies

Philipp Thalhammer↗

Testing SgrA{sup *} with the spectrum of its accretion structure

SgrA{sup *} is the supermassive black hole candidate at the center of the Galaxy and an ideal laboratory to test general relativity. Following previous work by other authors, we use the Polish doughnut model to describe an optically thin and constant angular momentum ion torus in hydrodynamical equilibrium and model the accretion structure around SgrA{sup *}. The radiation mechanisms are bremsstrahlung, synchrotron emission, and inverse Compton scattering. We compute the spectrum as seen by a distant observer in Kerr and non-Kerr spacetimes and we study how an accurate measurement can constrain possible deviations form the Kerr solution. As in the case of emission from a thin accretion disk, we find a substantial degeneracy between the determination of the spin and of possible deviations from the Kerr geometry, even when the parameters of the ion torus are fixed. This means that this technique cannot independently test the nature of SgrA{sup *} even in the presence of good data and with the systematics under control. However, it might do it in combination with other measurements (black hole shadow, radio pulsar, etc.)

79 ASTRONOMY AND ASTROPHYSICS↗

Formation and evolution of protostellar accretion discs – I. Angular-momentum budget, gravitational self-regulation, and numerical convergence

Here, we investigate the formation and early evolution of a protostellar disc from a magnetized prestellar core using non-ideal magnetohydrodynamic (MHD) simulations including ambipolar diffusion and Ohmic dissipation. The dynamical contraction of the prestellar core ultimately leads to the formation of a first hydrostatic core, after ambipolar diffusion decouples the magnetic field from the predominantly neutral gas. The hydrostatic core accumulates angular momentum from the infalling material, evolving into a rotationally supported torus; this ‘first hydrostatic torus’ then forms an accreting protostar and a rotationally supported disc. The disc spreads out by gravitational instability, reaching ~30 au in diameter at ~3 kyr after protostar formation. The total mass and angular momentum of the protostar–disc system are determined mainly by accretion of gas from an infalling pseudo-disc, which has low specific angular momentum because of magnetic braking; their removal from the protostar–disc system by outflow and disc magnetic braking are negligible, in part because the magnetic field is poorly coupled there. The redistribution of angular momentum within the protostar–disc system is facilitated mainly by gravitational instability; this allows formation of relatively large discs even when the specific angular momentum of infalling material is low. We argue that such discs should remain marginally unstable as they grow (with Toomre Q ~ 1–2), an idea that is broadly consistent with recent observational estimates for Class 0/I discs. We discuss the numerical convergence of our results, and show that properly treating the inner boundary condition is crucial for achieving convergence at an acceptable computational cost.

79 ASTRONOMY AND ASTROPHYSICS↗

Disc tearing leads to low and high frequency quasi-periodic oscillations in a GRMHD simulation of a thin accretion disc

Black hole X-ray binaries (BHXRBs) display a wide range of variability phenomena, from long duration spectral state changes to short-term broad-band variability and quasi-periodic oscillations (QPOs). A particularly puzzling aspect is the production of QPOs, which – if properly understood – could be used as a powerful diagnostic tool of black hole accretion and evolution. In this work, we analyse a high-resolution 3D general relativistic magnetohydrodynamic simulation of a geometrically thin accretion disc, which is tilted by 65° with respect to the black hole spin axis. We find that the Lense–Thirring torque from the rapidly spinning 10 M ⊙ black hole causes several sub-discs to tear off within ~10–20 gravitational radii. Tearing occurs in cycles on time-scales of seconds. During each tearing cycle, the inner sub-disc precesses for 1–5 periods before it falls into the black hole. We find a precession frequency of ~3 Hz, consistent with observed low-frequency QPOs. In addition, we find a high frequency QPO (HFQPO) with centroid frequency of ~55 Hz in the power spectra of the mass-weighted radius of the inner disc and the radial mass flux. This signal is caused by radial epicyclic oscillations of a dense ring of gas at the tearing radius, which suggests a corresponding modulation of the X-ray light curve and may thus explain some of the observed HFQPOs.

79 ASTRONOMY AND ASTROPHYSICS↗

Probing plasma physics with spectral index maps of accreting black holes on event horizon scales

The Event Horizon Telescope (EHT) collaboration has produced the first resolved images of the supermassive black holes at the centre of our galaxy and at the centre of the elliptical galaxy M87. As both technology and analysis pipelines improve, it will soon become possible to produce spectral index maps of black hole accretion flows on event horizon scales. Here, in this work, we predict spectral index maps of both M87* and Sgr A* by applying the general relativistic radiative transfer (GRRT) code IPOLE to a suite of general relativistic magnetohydrodynamic (GRMHD) simulations. We analytically show that the spectral index increases with increasing magnetic field strength, electron temperature, and optical depth. Consequently, spectral index maps grow more negative with increasing radius in almost all models, since all of these quantities tend to be maximized near the event horizon. Additionally, photon ring geodesics exhibit more positive spectral indices, since they sample the innermost regions of the accretion flow with the most extreme plasma conditions. Spectral index maps are sensitive to highly uncertain plasma heating prescriptions (the electron temperature and distribution function). However, if our understanding of these aspects of plasma physics can be tightened, even the spatially unresolved spectral index around 230 GHz can be used to discriminate between models. In particular, Standard and Normal Evolution (SANE) flows tend to exhibit more negative spectral indices than Magnetically Arrested Disc (MAD) flows due to differences in the characteristic magnetic field strength and temperature of emitting plasma.

(magnetohydrodynamics) MHD↗

Local models of two-temperature accretion disc coronae – II. Ion thermal conduction and the absence of disc evaporation

We use local stratified shearing-box simulations with magnetic field-aligned thermal conduction to study an idealized model of the coupling between a cold, radiatively efficient accretion disc, and an overlying, hot, two-temperature corona. Evaporation of a cold disc by conduction from the hot corona has been proposed as a means of mediating the soft-to-hard state transitions observed in X-ray binary systems. We model the coronal plasma in our local disc patch as an MHD fluid subject to both free-streaming ion conduction and a parametrized cooling function that captures the collisional transfer of energy from hot ions to colder, rapidly cooling leptons. In all of our models, independent of the initial net vertical magnetic flux (NF) threading the disc, we find no evidence of disc evaporation. The ion heat flux into the disc is radiated away before conduction can heat the disc’s surface layers. When an initial NF is present, steady-state temperature, density, and outflow velocities in our model coronae are unaffected by conduction. Instead of facilitating disc evaporation, thermal conduction is more likely to feed the disc with plasma condensing out of the corona, particularly in flows without NF. Our work indicates that uncertainties in the amount of NF threading the disc hold far greater influence over whether or not the disc will evaporate into a radiatively inefficient accretion flow compared to thermal conduction. We speculate that a change in net flux mediates disc truncation/evaporation.

79 ASTRONOMY AND ASTROPHYSICS↗

Nuclear-dominated accretion flows in two dimensions – II. Ejecta dynamics and nucleosynthesis for CO and ONe white dwarfs

We study mass ejection from accretion discs formed in the merger of a white dwarf with a neutron star or black hole. These discs are mostly radiatively inefficient and support nuclear fusion reactions, with ensuing outflows and electromagnetic transients. Here we perform time-dependent, axisymmetric hydrodynamic simulations of these discs including a physical equation of state, viscous angular momentum transport, a coupled 19-isotope nuclear network, and self-gravity. We find no detonations in any of the configurations studied. Our global models extend from the central object to radii much larger than the disc. We evolve these global models for several orbits, as well as alternate versions with an excised inner boundary to much longer times. We obtain robust outflows, with a broad velocity distribution in the range 10 2 –10 4 km s -1 . The outflow composition is mostly that of the initial white dwarf, with burning products mixed in at the |${\lesssim } 10\rm {-}30{{\ \rm per\ cent}}$| level by mass, including up to ~10 -2 M ⊙ of 56 Ni. These heavier elements (plus 4 He) are ejected within ≲ 40° of the rotation axis, and should have higher average velocities than the lighter elements that make up the white dwarf. These results are in broad agreement with previous one- and two-dimensional studies, and point to these systems as progenitors of rapidly rising (~ few day) transients. If accretion on to the central BH/NS powers a relativistic jet, these events could be accompanied by high-energy transients with peak luminosities ~10 47 –10 50 erg s -1 and peak durations of up to several minutes, possibly accounting for events like CDF-S XT2.

79 ASTRONOMY AND ASTROPHYSICS↗

The Properties of Reconnection Current Sheets in GRMHD Simulations of Radiatively Inefficient Accretion Flows

We report that non-ideal magnetohydrodynamic (MHD) effects may play a significant role in determining the dynamics, thermal properties, and observational signatures of radiatively inefficient accretion flows onto black holes. In particular, particle acceleration during magnetic reconnection events may influence black hole spectra and flaring properties. We use representative general relativistic magnetohydrodynamic (GRMHD) simulations of black hole accretion flows to identify and explore the structures and properties of current sheets as potential sites of magnetic reconnection. In the case of standard and normal evolution (SANE) disks, we find that in the reconnection sites, the plasma beta ranges from 0.1 to 1000, the magnetization ranges from 10 -4 to 1, and the guide fields are weak compared with the reconnecting fields. In magnetically arrested (MAD) disks, we find typical values for plasma beta from 10 -2 to 10 3 , magnetizations from 10 -3 to 10, and typically stronger guide fields, with strengths comparable to or greater than the reconnecting fields. These are critical parameters that govern the electron energy distribution resulting from magnetic reconnection and can be used in the context of plasma simulations to provide microphysics inputs to global simulations. Lastly, we also find that ample magnetic energy is available in the reconnection regions to power the fluence of bright X-ray flares observed from the black hole in the center of the Milky Way.

43 PARTICLE ACCELERATORS↗

Continuum Reverberation Mapping of the Accretion Disks in Two Seyfert 1 Galaxies

We present optical continuum lags for two Seyfert 1 galaxies, MCG+08-11-011 and NGC 2617, using monitoring data from a reverberation mapping campaign carried out in 2014. Our light curves span the ugriz filters over four months, with median cadences of 1.0 and 0.6 days for MCG+08-11-011 and NGC 2617, respectively, combined with roughly daily X-ray and near-UV data from Swift for NGC 2617. We find lags consistent with geometrically thin accretion-disk models that predict a lag-wavelength relation of τ ∝ λ 4/3 . However, the observed lags are larger than predictions based on standard thin-disk theory by factors of 3.3 for MCG+08-11-011 and 2.3 for NGC 2617. Furthermore, these differences can be explained if the mass accretion rates are larger than inferred from the optical luminosity by a factor of 4.3 in MCG+08-11-011 and a factor of 1.3 in NGC 2617, although uncertainty in the SMBH masses determines the significance of this result. While the X-ray variability in NGC 2617 precedes the UV/optical variability, the long (2.6 day) lag is problematic for coronal reprocessing models.

79 ASTRONOMY AND ASTROPHYSICS↗