Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Unity”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

SN Ia Standardization on the Rise: Evidence for the Cosmological Importance of Pre-maximum Measurements

We introduce SALT2X, an extension of the SALT2 model for SN Ia light curves. SALT2X separates the light-curve-shape parameter x 1 into an ${x}_{1}^{r}$ and ${x}_{1}^{f}$ for the rise and fall portions of the light curve. Using the Joint Lightcurve Analysis SN sample, we analyze the importance of the rising and falling portions of the light curve for cosmological standardization using a modified version of the Unified Nonlinear Inference for Type Ia cosmologY (UNITY) framework. We find strong evidence that ${x}_{1}^{r}$ has a stronger correlation with peak magnitude than ${x}_{1}^{f}$. We observe evidence that standardizing on the rise affects the color standardization relation, and reduces the size of the host-galaxy standardization and the unexplained ("intrinsic") luminosity dispersion. Since SNe Ia typically rise more quickly than they decline, a faster observing cadence in future surveys will be necessary to maximize the gain from this work and to continue to explore the impacts of decoupling the rising and falling portions of SN Ia light curves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

A Measurement of the Cosmic Microwave Background Lensing Potential and Power Spectrum from 500 deg 2 of SPTpol Temperature and Polarization Data

In this work, we present a measurement of the cosmic microwave background (CMB) lensing potential using 500 deg$^2$ of 150 GHz data from the SPTpol receiver on the South Pole Telescope. The lensing potential is reconstructed with signal-to-noise per mode greater than unity at lensing multipoles $L \lesssim 250$, using a quadratic estimator on a combination of CMB temperature and polarization maps. We report measurements of the lensing potential power spectrum in the multipole range of $100< L < 2000$ from sets of temperature-only, polarization-only, and minimum-variance estimators. We measure the lensing amplitude by taking the ratio of the measured spectrum to the expected spectrum from the best-fit $\Lambda$CDM model to the $\textit{Planck}$ 2015 TT+lowP+lensing dataset. For the minimum-variance estimator, we find $A_{\rm{MV}} = 0.944 \pm 0.058{\rm (Stat.)}\pm0.025{\rm (Sys.)}$; restricting to only polarization data, we find $A_{\rm{POL}} = 0.906 \pm 0.090 {\rm (Stat.)} \pm 0.040 {\rm (Sys.)}$. Considering statistical uncertainties alone, this is the most precise polarization-only lensing amplitude constraint to date (10.1 $\sigma$), and is more precise than our temperature-only constraint. We perform null tests and consistency checks and find no evidence for significant contamination.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Anatomy of a Cooling Flow: The Feedback Response to Pure Cooling in the Core of the Phoenix Cluster

We present new, deep observations of the Phoenix cluster from the Chandra X-ray Observatory, the Hubble Space Telescope, and the Karl Jansky Very Large Array. These data provide an order of magnitude improvement in depth and/or angular resolution at X-ray, optical, and radio wavelengths, yielding an unprecedented view of the core of the Phoenix cluster. We find that the one-dimensional temperature and entropy profiles are consistent with expectations for pure-cooling hydrodynamic simulations and analytic descriptions of homogeneous, steady-state cooling flow models. In the inner ~10 kpc, the cooling time is shorter by an order of magnitude than any other known cluster, while the ratio of the cooling time to freefall time approaches unity, signaling that the ICM is unable to resist multiphase condensation on kpc scales. When we consider the thermodynamic profiles in two dimensions, we find that the cooling is highly asymmetric. The bulk of the cooling in the inner ~20 kpc is confined to a low-entropy filament extending northward from the central galaxy. We detect a substantial reservoir of cool (10^4 K) gas (as traced by the [OII] doublet), which is coincident with the low-entropy filament. The bulk of this cool gas is draped around and behind a pair of X-ray cavities, presumably bubbles that have been inflated by radio jets, which are detected for the first time on kpc scales. These data support a picture in which AGN feedback is promoting the formation of a multiphase medium via a combination of ordered buoyant uplift and locally enhanced turbulence. These processes ought to counteract the tendency for buoyancy to suppress condensation, leading to rapid cooling along the jet axis. The recent mechanical outburst has sufficient energy to offset cooling, and appears to be coupling to the ICM via a cocoon shock, raising the entropy in the direction orthogonal to the radio jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Galaxy–Halo Connection of DESI Luminous Red Galaxies with Subhalo Abundance Matching

Abstract We use subhalo abundance and age distribution matching to create magnitude-limited mock galaxy catalogs at z ∼ 0.43, 0.52, and 0.63 with z -band and 3.4 μ m W 1-band absolute magnitudes and r − z and r − W 1 colors. From these magnitude-limited mocks, we select mock luminous red galaxy (LRG) samples according to the ( r − z )-based (optical) and ( r − W 1)-based (infrared) selection criteria for the LRG sample of the Dark Energy Spectroscopic Instrument (DESI) survey. Our models reproduce the number densities, luminosity functions, color distributions, and projected clustering of the DESI Legacy Surveys that are the basis for DESI LRG target selection. We predict the halo occupation statistics of both optical and IR DESI LRGs at fixed cosmology and assess the differences between the two LRG samples. We find that IR-based SHAM modeling represents the differences between the optical and IR LRG populations better than using the z band and that age distribution matching overpredicts the clustering of LRGs, implying that galaxy color is uncorrelated with halo age in the LRG regime. Both the optical and IR DESI LRG target selections exclude some of the most luminous galaxies that would appear to be LRGs based on their position on the red sequence in optical color–magnitude space. Both selections also yield populations with a nontrivial LRG–halo connection that does not reach unity for the most massive halos. We find that the IR selection achieves greater completeness (≳90%) than the optical selection across all redshift bins studied.

Berti, Angela M. (ORCID:0000000335826649)↗
DOE OSTI2023

Constraining Inputs to Realistic Kilonova Simulations through Comparison to Observed r -process Abundances

Kilonovae, one source of electromagnetic emission associated with neutron star mergers, are powered by the decay of radioactive isotopes in the neutron-rich merger ejecta. Models for kilonova emission consistent with the electromagnetic counterpart to GW170817 predict characteristic abundance patterns, determined by the relative balance of different types of material in the outflow. Assuming that the observed source is prototypical, this inferred abundance pattern in turn must match r -process abundances deduced by other means, such as what is observed in the solar system. We report on analysis comparing the input mass-weighted elemental compositions adopted in our radiative transfer simulations to the mass fractions of elements in the Sun, as a practical prototype for the potentially universal abundance signature from neutron star mergers. We characterize the extent to which our parameter inference results depend on our assumed composition for the dynamical and wind ejecta and examine how the new results compare to previous work. We find that a dynamical ejecta composition calculated using the FRDM2012 nuclear mass and FRLDM fission models with extremely neutron-rich ejecta ($Y_e$ = 0.035) along with moderately neutron-rich ($Y_e$ = 0.27) wind ejecta composition yields a wind-to-dynamical mass ratio of $M_w$/$M_d$ = 0.47, which best matches the observed AT2017gfo kilonova light curves while also producing the bestmatching abundance of neutron capture elements in the solar system, though, allowing for systematics, the ratio may be as high as of order unity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Electron Heating in the Transrelativistic Perpendicular Shocks of Tilted Accretion Flows

Abstract General relativistic magnetohydrodynamic (GRMHD) simulations of black hole tilted disks—where the angular momentum of the accretion flow at large distances is misaligned with respect to the black hole spin—commonly display standing shocks within a few to tens of gravitational radii from the black hole. In GRMHD simulations of geometrically thick, optically thin accretion flows, applicable to low-luminosity sources like Sgr A* and M87*, the shocks have transrelativistic speed, moderate plasma beta (the ratio of ion thermal pressure to magnetic pressure is β pi1 ∼ 1–8), and low sonic Mach number (the ratio of shock speed to sound speed is M s ∼ 1–6). We study such shocks with 2D particle-in-cell simulations, and we quantify the efficiency and mechanisms of electron heating for the special case of preshock magnetic fields perpendicular to the shock direction of propagation. We find that the postshock electron temperature T e2 exceeds the adiabatic expectation T e2,ad by an amount T e 2 / T e 2 , ad − 1 ≃ 0.0016 M s 3.6 , nearly independent of the plasma beta and of the preshock electron-to-ion temperature ratio T e1 / T i1 , which we vary from 0.1 to unity. We investigate the heating physics for M s ∼ 5–6 and find that electron superadiabatic heating is governed by magnetic pumping at T e1 / T i1 = 1, whereas heating by B -parallel electric fields (i.e., parallel to the local magnetic field) dominates at T e1 / T i1 = 0.1. Our results provide physically motivated subgrid prescriptions for electron heating at the collisionless shocks seen in GRMHD simulations of black hole accretion flows.

Astronomy & Astrophysics↗
DOE OSTI2024

Multizone Modeling of Black Hole Accretion and Feedback in 3D GRMHD: Bridging Vast Spatial and Temporal Scales

Simulating accretion and feedback from the horizon scale of supermassive black holes (SMBHs) out to galactic scales is challenging because of the vast range of scales involved. Elaborating on H. Cho et al., we describe and test a "multizone" technique, which is designed to tackle this difficult problem in three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. While short-timescale variability should be interpreted with caution, the method is demonstrated to be well-suited for finding dynamical steady states over a wide dynamic range. We simulate accretion on a nonspinning SMBH ($a\ast$ = 0) using initial conditions and the external galactic potential from a large-scale galaxy simulation and achieve a steady state over eight decades in radius. As found in H. Cho et al., the density scales with radius as ρ ∝ r –1 inside the Bondi radius R B , which is located at R B = 2 × 10 5 r g (≈60 pc for M87), where r g is the gravitational radius of the SMBH; the plasma-β is ~ unity, indicating an extended magnetically arrested state; the mass accretion rate $\dot{M}$ is ≈1% of the analytical Bondi accretion rate ${\dot{M}}_{{\rm{B}}};$ and there is continuous energy feedback out to ≈100R B (or beyond > kpc) at a rate $\approx 0.02\dot{M}{c}^{2}$. Surprisingly, no ordered rotation in the external medium survives as the magnetized gas flows to smaller radii, and the final steady solution is very similar to when the exterior has no rotation. Using the multizone method, we simulate GRMHD accretion over a wide range of Bondi radii, R B ~ 10 2 –10 7 r g , and find that $\dot{M}/{\dot{M}}_{{B}}\approx {({R}_{{B}}/6\,{r}_{g})}^{-0.5}$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The Multiple Paths to Merger of Unequal-mass Black Hole Binaries in the Disks of Active Galactic Nuclei

The accretion disks that power active galactic nuclei (AGN) are thought to house populations of stars and compact objects; after forming binaries these compact objects may merge, begetting gravitational waves (GWs) such as those detected by LIGO and VIRGO. We present a comprehensive study of the early evolution of binaries within AGN disks as their orbits are influenced by the surrounding gas, focusing on eccentric and unequal-mass binaries. Nearly equal-mass binaries behave similarly to their equal-mass counterparts: Prograde binaries inspiral, albeit somewhat slowly, and have their eccentricities damped; retrograde binaries inspiral ∼2–3 times faster than their prograde counterparts, and those with near-equal masses are driven quickly toward near-unity eccentricities. However, the primaries in retrograde binaries with mass ratios of m 2 /m 1 ≲ 0.4 experience significantly weaker headwinds and retain substantial accretion disks that help damp binary eccentricities, slowing binary inspirals. Additionally, we find that while accretion drives prograde binaries toward equal masses thanks to the exchange of material between the primary and secondary accretion disks, retrograde binaries are driven slowly toward more extreme mass ratios. Prograde binaries, and generally those with low mass ratios, likely accrete for multiple e-folding timescales before merger. On the other hand, high-mass-ratio retrograde binaries may merge before accreting substantially, potentially approaching merger with detectable eccentricity. Future ground-based GW observatories, with their broader frequency coverage, should be particularly useful for studying these populations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Simulation to a Newborn Supernova Remnant from a Low-mass Iron Core Star

Supernova remnant observations show a high degree of asymmetry, mixing, and inhomogeneity. These asymmetries are seeded during the early seconds of the explosion and are further enhanced and modified as the shock and ejecta move through the stellar progenitor and into the circumstellar medium. We present simulations of a 9.6 M⊙ zero-metallicity progenitor initialized after shock revival and evolved for several years when the ejecta is in the circumstellar medium. A suite of 1D and 2D simulations examines the effects of neutron-star wind and radioactive decay heating. In 1D, decay heating forms a low-density bubble that suppresses the reverse shock. While in 2D, the heating is localized to metal-rich pockets, inflating them and compressing the surrounding material into dense shells. In 3D, the neutron-star wind and decay heating modify the plume morphology, producing more large-scale structures. The extended plume morphology leads to an asymmetrical shock breakout. After breakout, the leading plumes cannot keep up with the shock front, resulting in deceleration and fragmentation by the reverse shock while retaining the large-scale asymmetry. The projected ejecta morphology and velocities are strongly viewing angle dependent. The relatively uniform metal-rich distribution does not resemble the strongly inhomogeneous ejecta structure of Cas A. The 160-isotope decay network shows that 24.4% of the radioactive heating comes from decay chains other than the canonical 56Ni chain. The low explosion energy, low 56Ni yield, and Ni/Fe ratio greater than unity suggest an observational signature similar to an electron capture supernova.

Neopane, Sudarshan [University of Tennessee (UT)]↗
DOE OSTI2017

Mira Variable Stars from LAMOST DR4 Data: Emission Features, Temperature Types, and Candidate Selection

Based on an extensive spectral study of a photometrically confirmed sample of Mira variables, we find a relationship between the relative Balmer emission-line strengths and spectral temperatures of O-rich Mira stars. The F{sub Hδ}/F{sub Hγ} flux ratio increases from less than unity to five as stars cool down from M0 to M10, which is likely driven by increasing TiO absorption above the deepest shock-emitting regions. We also discuss the relationship between the equivalent widths of the Balmer emission lines and the photometric luminosity phase of our Mira sample stars. Using our 291 Mira spectra as templates for reference, 191 Mira candidates are newly identified from the LAMOST DR4 catalog. We summarize the criteria adopted to select Mira candidates based on emission-line indices and molecular absorption bands. This enlarged spectral sample of Mira variables has the potential to contribute significantly to our knowledge of the optical properties of Mira stars and will facilitate further studies of these late-type, long-period variables.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Analytical Models of Exoplanetary Atmospheres. IV. Improved Two-stream Radiative Transfer for the Treatment of Aerosols

We present a novel generalization of the two-stream method of radiative transfer, which allows for the accurate treatment of radiative transfer in the presence of strong infrared scattering by aerosols. We prove that this generalization involves only a simple modification of the coupling coefficients and transmission functions in the hemispheric two-stream method. This modification originates from allowing the ratio of the first Eddington coefficients to depart from unity. At the heart of the method is the fact that this ratio may be computed once and for all over the entire range of values of the single-scattering albedo and scattering asymmetry factor. We benchmark our improved two-stream method by calculating the fraction of flux reflected by a single atmospheric layer (the reflectivity) and comparing these calculations to those performed using a 32-stream discrete-ordinates method. We further compare our improved two-stream method to the two-stream source function (16 streams) and delta-Eddington methods, demonstrating that it is often more accurate at the order-of-magnitude level. Finally, we illustrate its accuracy using a toy model of the early Martian atmosphere hosting a cloud layer composed of carbon dioxide ice particles. The simplicity of implementation and accuracy of our improved two-stream method renders it suitable for implementation in three-dimensional general circulation models. In other words, our improved two-stream method has the ease of implementation of a standard two-stream method, but the accuracy of a 32-stream method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ALMA 0.1–0.2 arcsec resolution imaging of the NGC 1068 Nucleus: compact dense molecular gas emission at the putative AGN location

We present the results of our ALMA Cycle 2 high angular resolution (0.″1–0.″2) observations of the nuclear region of the nearby well-studied type-2 active galactic nucleus (AGN), NGC 1068, at HCN J = 3–2 and HCO{sup +} J = 3–2 emission lines. For the first time, due to a higher angular resolution than previous studies, we clearly detected dense molecular gas emission at the putative AGN location, identified as a ∼1.1 mm (∼266 GHz) continuum emission peak, by separating this emission from brighter emission located at 0.″5–2.″0 on the eastern and western sides of the AGN. The estimated intrinsic molecular emission size and dense molecular mass, which are thought to be associated with the putative dusty molecular torus around an AGN, were ∼10 pc and ∼several × 10{sup 5} M {sub ⊙}, respectively. HCN-to-HCO{sup +} J = 3–2 flux ratios substantially higher than unity were found throughout the nuclear region of NGC 1068. The continuum emission displayed an elongated morphology along the direction of the radio jet located at the northern side of the AGN, as well as a weak spatially-resolved component at ∼2.″0 on the southwestern side of the AGN. The latter component most likely originated from star formation, with the estimated luminosity more than one order of magnitude lower than the luminosity of the central AGN. No vibrationally excited ( v {sub 2} = 1f) J = 3–2 emission lines were detected for HCN and HCO{sup +} across the field of view.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

AN ULTRA-FAST X-RAY DISK WIND IN THE NEUTRON STAR BINARY GX 340+0

We present a spectral analysis of a brief Chandra /HETG observation of the neutron star low-mass X-ray binary GX 340+0. The high-resolution spectrum reveals evidence of ionized absorption in the Fe K band. The strongest feature, an absorption line at approximately 6.9 keV, is required at the 5 σ level of confidence via an F -test. Photoionization modeling with XSTAR grids suggests that the line is the most prominent part of a disk wind with an apparent outflow speed of v = 0.04 c . This interpretation is preferred at the 4 σ level over a scenario in which the line is H-like Fe xxvi at a modest redshift. The wind may achieve this speed owing to its relatively low ionization, enabling driving by radiation pressure on lines; in this sense, the wind in GX 340+0 may be the stellar-mass equivalent of the flows in broad absorption line quasars. If the gas has a unity volume filling factor, the mass ouflow rate in the wind is over 10{sup −5} M {sub ⊙} yr{sup −1}, and the kinetic power is nearly 10{sup 39} erg s{sup −1} (or, 5–6 times the radiative Eddington limit for a neutron star). However, geometrical considerations—including a small volume filling factor and low covering factor—likely greatly reduce these values.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

DAMPING OF THE MILKY WAY BAR BY MANIFOLD-DRIVEN SPIRALS

We describe a new phenomenon of “bar damping” that may have played an important role in shaping the Milky Way bar and bulge as well as its spiral structure. We use a collisionless N -body simulation of a Milky Way–like galaxy initially composed of a dark matter halo and an exponential disk with a Toomre parameter slightly above unity. In this configuration, dominated by the disk in the center, a bar forms relatively quickly, after 1 Gyr of evolution. This is immediately followed by the formation of two manifold-driven spiral arms and the outflow of stars that modifies the potential in the vicinity of the bar, apparently shifting the position of the L {sub 1}/ L {sub 2} Lagrange points. This modification leads to the shortening of the bar and the creation of a next generation of manifold-driven spiral arms at a smaller radius. The process repeats itself a few times over the next 0.5 Gyr resulting in further substantial weakening and shortening of the bar. The time when the damping comes to an end coincides with the first buckling episode in the bar that rebuilds the orbital structure so that no more new spiral arms are formed. The morphology of the bar and the spiral structure at this time show remarkable similarity to the present properties of the Milky Way. Later on, the bar starts to grow rather steadily again, weakened only by subsequent buckling episodes occurring at more distant parts of the disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

High radiation pressure on interstellar dust computed by light-scattering simulation on fluffy agglomerates of magnesium-silicate grains with metallic-iron inclusions

Recent space missions have provided information on the physical and chemical properties of interstellar grains such as the ratio β of radiation pressure to gravity acting on the grains in addition to the composition, structure, and size distribution of the grains. Numerical simulation on the trajectories of interstellar grains captured by Stardust and returned to Earth constrained the β ratio for the Stardust samples of interstellar origin. However, recent accurate calculations of radiation pressure cross-sections for model dust grains have given conflicting stories in the β ratio of interstellar grains. The β ratio for model dust grains of so-called “astronomical silicate” in the femto-kilogram range lies below unity, in conflict with β∼1 for the Stardust interstellar grains. Here, I tackle this conundrum by re-evaluating the β ratio of interstellar grains on the assumption that the grains are aggregated particles grown by coagulation and composed of amorphous MgSiO{sub 3} with the inclusion of metallic iron. My model is entirely consistent with the depletion and the correlation of major rock-forming elements in the Local Interstellar Cloud surrounding the Sun and the mineralogical identification of interstellar grains in the Stardust and Cassini missions. I find that my model dust particles fulfill the constraints on the β ratio derived from not only the Stardust mission but also the Ulysses and Cassini missions. My results suggest that iron is not incorporated into silicates but exists as metal, contrary to the majority of interstellar dust models available to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Retrieval Analysis of the Emission Spectrum of WASP-12b: Sensitivity of Outcomes to Prior Assumptions and Implications for Formation History

We analyze the emission spectrum of the hot Jupiter WASP-12b using our HELIOS-R retrieval code and HELIOS-K opacity calculator. When interpreting Hubble and Spitzer data, the retrieval outcomes are found to be prior-dominated. When the prior distributions of the molecular abundances are assumed to be log-uniform, the volume mixing ratio of HCN is found to be implausibly high. A VULCAN chemical kinetics model of WASP-12b suggests that chemical equilibrium is a reasonable assumption even when atmospheric mixing is implausibly rigorous. Guided by (exo)planet formation theory, we set Gaussian priors on the elemental abundances of carbon, oxygen, and nitrogen with the Gaussian peaks being centered on the measured C/H, O/H, and N/H values of the star. By enforcing chemical equilibrium, we find substellar O/H and stellar to slightly superstellar C/H for the dayside atmosphere of WASP-12b. The superstellar carbon-to-oxygen ratio is just above unity, regardless of whether clouds are included in the retrieval analysis, consistent with Madhusudhan et al. Furthermore, whether a temperature inversion exists in the atmosphere depends on one’s assumption for the Gaussian width of the priors. Our retrieved posterior distributions are consistent with the formation of WASP-12b in a solar-composition protoplanetary disk, beyond the water iceline, via gravitational instability or pebble accretion (without core erosion) and migration inward to its present orbital location via a disk-free mechanism, and are inconsistent with both in situ formation and core accretion with disk migration, as predicted by Madhusudhan et al. We predict that the interpretation of James Webb Space Telescope WASP-12b data will not be prior-dominated.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Planet-induced Vortices with Dust Coagulation in Protoplanetary Disks

In this work, we study how the dust coagulation/fragmentation will influence the evolution and observational appearances of vortices induced by a massive planet embedded in a low-viscosity disk by performing global 2D high-resolution hydrodynamical simulations. Within the vortex, due to its higher gas surface density and steeper pressure gradients, dust coagulation, fragmentation, and drift (to the vortex center) are all quite efficient, producing dust particles ranging from 1 μm to ~1.0 cm, as well as an overall high dust-to-gas ratio (above unity). In addition, the dust size distribution is quite nonuniform inside the vortex, with the mass-weighted average dust size at the vortex center (~4.0 mm) being a factor of ~10 larger than other vortex regions. Both large (~millimeter) and small (tens of microns) particles contribute strongly to affect the gas motion within the vortex. As such, we find that the inclusion of dust coagulation has a significant impact on the vortex lifetime and the typical vortex lifetime is about 1000 orbits. Furthermore, after the initial gaseous vortex is destroyed, the dust spreads into a ring with a few remaining smaller gaseous vortices with a high dust concentration and a large maximum size (~millimeter). At late time, the synthetic dust continuum images for the coagulation case show as a ring inlaid with several hot spots at the 1.33 mm band, while only distinct hot spots remain at 7.0 mm.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Scaling of Turbulent Viscosity and Resistivity: Extracting a Scale-dependent Turbulent Magnetic Prandtl Number

Turbulent viscosity νt and resistivity ηt are perhaps the simplest models for turbulent transport of angular momentum and magnetic fields, respectively. The associated turbulent magnetic Prandtl number Pr t ≡ ν t /η t has been well recognized to determine the final magnetic configuration of accretion disks. Here, we present an approach to determining these "effective transport" coefficients acting at different length scales using coarse-graining and recent results on decoupled kinetic and magnetic energy cascades. By analyzing the kinetic and magnetic energy cascades from a suite of high-resolution simulations, we show that our definitions of ν t , η t , and Pr t have power-law scalings in the "decoupled range." We observe that Pr t ≈ 1–2 at the smallest inertial-inductive scales, increasing to ≈5 at the largest scales. Furthermore, based on physical considerations, our analysis suggests that Pr t has to become scale independent and of order unity in the decoupled range at sufficiently high Reynolds numbers (or grid resolution) and that the power-law scaling exponents of velocity and magnetic spectra become equal. In addition to implications for astrophysical systems, the scale-dependent turbulent transport coefficients offer a guide for large-eddy simulation modeling.

79 ASTRONOMY AND ASTROPHYSICS↗