Simulations Of Disks Research Articles

Probing the eccentricity in protostellar discs: Modelling kinematics and morphologies

Context. Protostellar discs are mostly modelled as circular structures of gas and dust orbiting a protostar. However, a number of physical mechanisms, for example, the presence of a (sub)stellar companion or initial axial asymmetry, can cause the gas and dust orbital motion to become eccentric. Theoretical studies have revealed that, when present, disc eccentricity is expected to occur with predictable profiles that can be long-lasting and potentially observable in protostellar systems. Aims. We construct an analytical model predicting the typical features of the kinematics and morphology of eccentric protostellar discs, with the final goal of characterising the observational appearance of eccentricity in discs. Methods. We validate the model using a numerical simulation of a circumbinary disc (where the binary makes the disc eccentric). We finally post-process the simulation with Monte Carlo radiative transfer to study how eccentric features would appear through the ‘eyes’ of ALMA. Results. Besides the motion of the material on eccentric Keplerian orbits in the disc orbital plane, the most characteristic eccentric feature emerging from the analytical model is strong vertical motion with a typical anti-symmetric pattern (with respect to the disc line of pericentres). A circumbinary disc with a ≈ 40 au eccentric cavity (ecav = 0.2), carved by an abin = 15 au binary, placed at a distance d = 130 pc, is expected to host in its upper emission surface vertical oscillations up to vz ~ 400 m s−1 close to the cavity edge, that is to say, well within ALMA spectral and spatial resolution capabilities. A residual spiral pattern in the vertical velocity Δvz ~ 150 m s−1 of the simulation cannot be captured by the theoretical model, we speculate it to be possibly linked to the presence of a companion in the system.

Modeling of Blockage and Wake Effect: Comparison with Field data

We apply engineering approaches for estimating the induction zone in front of a wind turbine and escalating to the global blockage at wind farm and cluster of farms scale as well accounting for wind farm wake. The methods include the vortex cylinder model and the self-similar model for blockage, and Jensen and Eddy-viscosity models for the wake estimation. To calibrate each model, Reynold Average Navier-Stokes with Actuator disc simulations are employed. Every model’s efficacy is assessed in comparison to the validation data of an operated wind farm in Scotland with more than 100 turbines of 10 MW each. We also examine the widely used process within the wind industry for transferring the application of calibrated wake and blockage models from one wind farm. According to our simulation analysis, blockage causes sideways acceleration and upstream slowdown. This also result in a gradient in the power produced by each turbine on the farm, which, if not properly accounted for, appears to be a wake model constraint. Furthermore, the response of different turbine type in another environment is not effectively predicted by the blockage and wake model, which is based on a particular turbine type and environment.

The interplay between forming planets and photoevaporating discs

Context. Disc winds and planet–disc interactions are two crucial mechanisms that define the structure, evolution, and dispersal of protoplanetary discs. While winds are capable of removing material from discs, eventually leading to their dispersal, massive planets can shape their disc by creating sub-structures such as gaps and spiral arms. Aims. We studied the interplay between an X-ray photoevaporative disc wind and the sub-structures generated due to planet–disc interactions to determine how their mutual interactions affect the disc’s and the planet’s evolution. Methods. We performed 3D hydrodynamic simulations of viscous discs (α = 6.9 × 10−4) that host a Jupiter-like planet and undergo X-ray photoevaporation. We traced the gas flows within the disc and wind and measured the rate of accretion onto the planet, as well as the gravitational torque that is acting on it. Results. Our results show that the planetary gap removes the wind’s pressure support, allowing wind material to fall back into the gap. This opens new pathways for material from the inner disc (and part of the outer disc) to be redistributed through the wind towards the gap. Consequently, the gap becomes shallower and the flow of mass across the gap in both directions is significantly increased, as is the planet’s mass-accretion rate (by factors of ≈5 and ≈2, respectively). Moreover, the wind-driven redistribution results in a denser inner disc and a less dense outer disc, which, combined with the recycling of a significant portion of the inner wind, leads to longer lifetimes for the inner disc, contrary to the expectation in a planet-induced photoevaporation scenario that has been proposed in the past.

Thermo-mechanical effects and microstructural evolution-coupled numerical simulation on the hot forming processes of superalloy turbine disk

Abstract Macroscopic deformation and microstructural evolution simultaneously exist in the hot forming processes of superalloy. In order to effectively and accurately study the hot forming processes of superalloy turbine disk with the numerical simulation method, a multi-scale finite element model of GH4065 superalloy turbine disk involving macroscopic and microscopic aspects was established by defining macro- and micromaterial model of superalloy, hot forming processing parameters, and boundary conditions. Via the numerical simulations of superalloy turbine disk, the macroscopic material flow and microstructural evolution behaviors in the hot forming processes of superalloy turbine disk were studied. Besides, the macroscopic deformation and microstructure distribution states after the hot forming processes were revealed and analyzed. A corresponding hot forming physical test of superalloy turbine disk was conducted to verify the results of the numerical simulation. Via the qualitative and quantitative analyses, it was concluded that the macroscopic deformation and microstructural evolution in the hot forming processes of superalloy turbine disk can be accurately predicted by the numerical simulation method.

Introducing two improved methods for approximating radiative cooling in hydrodynamical simulations of accretion discs

ABSTRACT The evolution of many astrophysical systems depends strongly on the balance between heating and cooling, in particular star formation in giant molecular clouds and the evolution of young protostellar systems. Protostellar discs are susceptible to the gravitational instability, which can play a key role in their evolution and in planet formation. The strength of the instability depends on the rate at which the system loses thermal energy. To study the evolution of these systems, we require radiative cooling approximations because full radiative transfer is generally too expensive to be coupled to hydrodynamical models. Here, we present two new approximate methods for computing radiative cooling that make use of the polytropic cooling approximation. This approach invokes the assumption that each parcel of gas is located within a spherical pseudo-cloud, which can then be used to approximate the optical depth. The first method combines the methods introduced by Stamatellos et al. and Lombardi et al. to overcome the limitations of each method at low and high optical depths, respectively. The second method, the ‘modified Lombardi’ method, is specifically tailored for self-gravitating discs. This modifies the scale height estimate from the method of Lombardi et al. using the analytical scale height for a self-gravitating disc. We show that the modified Lombardi method provides an excellent approximation for the column density in a fragmenting disc, a regime in which the existing methods fail to recover the clumps and spiral structures. We therefore recommend this improved radiative cooling method for more realistic simulations of self-gravitating discs.

Time-dependent, long-term hydrodynamic simulations of the inner protoplanetary disk. III. The influence of photoevaporation

The final stages of a protoplanetary disk are essential for our understanding of the formation and evolution of planets. Photoevaporation is an important mechanism that contributes to the dispersal of an accretion disk and has significant consequences for the disk's lifetime. However, the combined effects of photoevaporation and star-disk interaction have not been investigated in previous studies. A photoevaporative disk evolution model including a detailed formulation of the inner star-disk interaction region will improve the understanding of the final stages of disk evolution. We combined an implicit disk evolution model with a photoevaporative mass-loss profile. By including the innermost disk regions down to 0.01 AU, we could calculate the star-disk interaction, the stellar spin evolution, and the transition from an accreting disk to the propeller regime self-consistently. Starting from an early Class II star-disk system (with an age of 1 Myr), we calculated the long-term evolution of the system until the disk becomes almost completely dissolved. Photoevaporation has a significant effect on disk structure and evolution. The radial extent of the dead zone decreases, and the number of episodic accretion events (outbursts) is reduced by high stellar X-ray luminosities. Reasonable accretion rates ($> 10^ M_ / yr $) in combination with photoevaporative gaps are possible for a dead zone that is still massive enough to develop episodic accretion events. Furthermore, the stellar spin evolution during the Class II evolution is less affected by the star-disk interaction in the case of high X-ray luminosities. Our results suggest that the formation of planets, especially habitable planets, in the dead zone is strongly impaired in the case of strong X-ray luminosities. Additionally, the importance of the star-disk interaction during the Class II phase with respect to the stellar spin evolution is reduced.

Particle Bounce Stick Behavior in the Rotating Frame of Reference

Abstract Particle deposition is a major damage mechanism for gas turbine components, especially in the secondary air system. Predicting the transport and deposition of ingested atmospheric contaminants is of great interest in component-level simulations. Bounce stick models predict deposition upon collision with a wall during Lagrangian particle tracking; they consider a range of physical phenomena, including van der Waals forces and plastic deformation. The effect of the rotating frame of reference on particle collision physics has thus far been neglected in the literature despite the significant centrifugal (CF) loads experienced by many components of interest for deposition studies, for example, turbine blades. The collision physics of the rotating frame are discussed here using low-order models and Monte Carlo style simulations. The present work aims to provide a conceptual framework for the inclusion of CF forces in collision physics. No significant effect was found on the rebound velocities of particles experiencing “worst case” CF forces. However, differences were observed in the sticking probability between concave and convex rotating surfaces, where the CF forces act in opposing directions to either push the particle into the surface or detach it. A force-based analogy of the popular critical velocity model was developed to study the phenomenon. It was used in a Monte Carlo style simulation of a rotor disk, finding that the variation of critical velocity due to CF forces was likely to bias deposition toward concave surfaces. The collision physics of the rotating frame were found to be unintuitive, with complex implications for modeling deposition in gas turbine components. However, they were consequential in determining the distribution of deposition through the system, hence should be included in particle deposition simulations in computational fluid dynamics (CFD).

Insights into the Galactic Bulge Chemodynamical Properties from Gaia Data Release 3

We explore the chemodynamical properties of the Galaxy in the azimuthal velocity V ϕ and metallicity [Fe/H] space using red giant stars from Gaia Data Release 3. The row-normalized V ϕ –[Fe/H] maps form a coherent sequence from the bulge to the outer disk, clearly revealing the thin/thick disk and the Splash. The metal-rich stars display bar-like kinematics, while the metal-poor stars show dispersion-dominated kinematics. The intermediate-metallicity population (−1 < [Fe/H]< − 0.4) can be separated into two populations, one that is bar-like, i.e., dynamically cold ( σVR∼80 km s−1) and fast-rotating (V ϕ ≳ 100 km s−1), and the Splash, which is dynamically hot ( σVR∼110 km s−1) and slow-rotating (V ϕ ≲ 100 km s−1). We compare the observations in the bulge region with an Auriga simulation where the last major merger event occurred ∼10 Gyr ago: only stars born around the time of the merger reveal a Splash-like feature in the V ϕ –[Fe/H] space, suggesting that the Splash is likely merger-induced, predominantly made up of heated disk stars and the starburst associated with the last major merger. Since the Splash formed from the proto-disk, its lower metallicity limit coincides with that of the thick disk. The bar formed later from the dynamically hot disk with [Fe/H] > − 1 dex, with the Splash not participating in the bar formation and growth. Moreover, with a set of isolated evolving N-body disk simulations, we confirm that a nonrotating classical bulge can be spun up by the bar and develop cylindrical rotation, consistent with the observations for the metal-poor stars.

Properties of the Line-of-sight Velocity Field in the Hot and X-Ray-emitting Circumgalactic Medium of Nearby Simulated Disk Galaxies

The hot, X-ray-emitting phase of the circumgalactic medium of massive galaxies is believed to be the reservoir of baryons from which gas flows onto the central galaxy and into which feedback from active galactic nuclei and stars inject mass, momentum, energy, and metals. These effects shape the velocity fields of the hot gas, which can be observed via the Doppler shifting and broadening of emission lines by X-ray integral field units. In this work, we analyze the gas kinematics of the hot circumgalactic medium of Milky Way–mass disk galaxies from the TNG50 simulation with synthetic observations to determine how future instruments can probe this velocity structure. We find that the hot phase is often characterized by outflows from the disk driven by feedback processes, radial inflows near the galactic plane, and rotation, although in some systems the velocity field is more disorganized and turbulent. With a spectral resolution of ∼1 eV, fast and hot outflows (∼200–500 km s−1) can be measured, depending on the orientation of the galaxy on the sky. The rotation velocity of the hot phase (∼100–200 km s−1) can be measured using line shifts in edge-on galaxies, and is slower than that of colder gas phases but similar to stellar rotation velocities. By contrast, the slow inflows (∼50–100 km s−1) are difficult to measure in projection with these other components, but may be detected in multicomponent spectral fits. We find that the velocity measured is sensitive to which emission lines are used. Measuring these flows will constrain theories of how the gas in these galaxies evolves.

Plasmoid identification and statistics in two-dimensional Harris sheet and GRMHD simulations

ABSTRACT Magnetic reconnection is a ubiquitous phenomenon for magnetized plasma and leads to the rapid reconfiguration of magnetic field lines. During reconnection events, plasma is heated and accelerated until the magnetic field lines enclose and capture the plasma within a circular configuration. These so-called plasmoids could therefore observationally manifest themselves as hotspots, which are associated with flaring behaviour in supermassive black hole systems, such as Sagittarius A*. We have developed a novel algorithm for identifying plasmoid structures, which incorporates watershed and custom closed contouring steps. From the identified structures, we determine the plasmoids’ plasma characteristics and energetics in magnetohydrodynamical simulations. The algorithm’s performance is showcased for a high-resolution suite of axisymmetric ideal and resistive magnetohydrodynamical simulations of turbulent accretion discs surrounding a supermassive black hole. For validation purposes, we also evaluate several Harris current sheets that are well-investigated in the literature. We recover the characteristic power-law distribution of plasmoid sizes for both the black hole and Harris sheet simulations. This indicates that while the dynamics are vastly different, with different dominant plasma instabilities, the plasmoid creation characteristics are similar. Plasmoid formation rates for resistive general relativistic magnetohydrodynamical simulations are significantly higher than for their ideal counterpart. Moreover, the largest identified plasmoids are consistent with sizes typically assumed for semi-analytical interpretation of observations. We recover a positive correlation between the plasmoid formation rate and decreases in black-hole-horizon-penetrating magnetic flux, during which the accretion flow is temporarily halted. These results demonstrate the efficacy of the newly developed algorithm which has enabled an extensive quantitative analysis of plasmoid formation in black hole accretion simulations.

Growth and evolution of low-mass planets in pressure bumps

Observations of protoplanetary disks have revealed dust rings that are likely due to the presence of pressure bumps in the disk. Because these structures tend to trap drifting pebbles, it has been proposed that pressure bumps may play an important role in the planet formation process. In this paper, we investigate the orbital evolution of a 0.1 M⊕ protoplanet embedded in a pressure bump using 2D hydrodynamical simulations of protoplanetary disks consisting of gas and pebbles. We examine the role of thermal forces generated by the pebble accretion-induced heat release, taking into account the feedback between the luminosity and the eccentricity. We also study the effect of the pebble-scattered flow on the planet’s orbital evolution. Due to the accumulation of pebbles at the pressure bump, the planet’s accretion luminosity is high enough to induce significant eccentricity growth through thermal forces. Accretion luminosity is also responsible for vortex formation at the planet’s position through baroclinic effects, which cause the planet to escape from the dust ring if dust feedback on the gas is neglected. Including the effect of the dust feedback leads to weaker vortices, which enable the planet to remain close to the pressure maximum on an eccentric orbit. Simulations in which the planet mass is allowed to increase as a consequence of pebble accretion result in the formation of giant planet cores with masses in the range of 5–20 M⊕ over ~2 × 104 yr. This occurs for moderate values of the Stokes number, St ≈ 0.01, such that the pebble drift velocity is not too high and the dust ring mass not too small. Our results suggest that pressure bumps mays be preferred locations for the formation of giant planets, but this requires a moderate level of grain growth within the disk.

Galactic disc heating by density granulation in fuzzy dark matter simulations

Abstract Fuzzy dark matter (FDM), an attractive dark matter candidate comprising ultralight bosons (axions) with a particle mass ma ∼ 10−22 eV, is motivated by the small-scale challenges of cold dark matter and features a kpc-size de Broglie wavelength. Quantum wave interference inside an FDM halo gives rise to stochastically fluctuating density granulation; the resulting gravitational perturbations could drive significant disc thickening, providing a natural explanation for galactic thick discs. Here we present the first self-consistent simulations of FDM haloes and stellar discs, exploring ma = 0.2–1.2 × 10−22 eV and halo masses Mh = 0.7–2.8 × 1011 M⊙. Disc thickening is observed in all simulated systems. The disc heating rates are approximately constant in time and increase substantially with decreasing ma, reaching dh/dt ≃ 0.04 (0.4) kpc Gyr−1 and $d\sigma _z^2/dt \simeq 4$ (150) km2s−2Gyr−1 for ma = 1.2 (0.2) × 10−22 eV and $M_{\rm h}=7\times 10^{10} \, \rm {M}_{\odot }$, where h is the disc scale height and σz is the vertical velocity dispersion. These simulated heating rates agree within a factor of two with the theoretical estimates of Chiang et al., confirming that the rough estimate of Church et al. overpredicts the granulation-driven disc heating rate by two orders of magnitude. However, the simulation-inferred heating rates scale less steeply than the theoretically predicted relation $d\sigma ^2_z/dt \propto m_a^{-3}$. Finally, we examine the applicability of the Fokker–Planck approximation in FDM granulation modelling and the robustness of the ma exclusion bound derived from the Galactic disc kinematics.

Cloud properties across spatial scales in simulations of the interstellar medium

Molecular clouds (MCs) are structures of dense gas in the interstellar medium (ISM) that extend from ten to a few hundred parsecs and form the main gas reservoir available for star formation. Hydrodynamical simulations of a varying complexity are a promising way to investigate MCs evolution and their properties. However, each simulation typically has a limited range in resolution and different cloud extraction algorithms are used, which complicates the comparison between simulations. In this work, we aim to extract clouds from different simulations covering a wide range of spatial scales. We compare their properties, such as size, shape, mass, internal velocity dispersion, and virial state. We applied the Hop cloud detection algorithm on (M)HD numerical simulations of stratified ISM boxes and isolated galactic disk simulations that were produced using Flash Ramses and Arepo We find that the extracted clouds are complex in shape, ranging from round objects to complex filamentary networks in all setups. Despite the wide range of scales, resolution, and sub-grid physics, we observe surprisingly robust trends in the investigated metrics. The mass spectrum matches in the overlap between simulations without rescaling and with a high-mass power-law index of -1 for logarithmic bins of mass, in accordance with theoretical predictions. The internal velocity dispersion scales with the size of the cloud as $ $ for large clouds ($R pc $). For small clouds we find larger sigma compared to the power-law scaling, as seen in observations, which is due to supernova-driven turbulence. Almost all clouds are gravitationally unbound with the virial parameter scaling as $ vir $, which is slightly flatter compared to observed scaling but in agreement given the large scatter. We note that the cloud distribution towards the low-mass end is only complete if the more dilute gas is also refined, rather than only the collapsing regions.

Measuring the dynamical length of galactic bars

ABSTRACT We define a physically motivated measure for galactic bar length, called the dynamical length. The dynamical length of the bar corresponds to the radial extent of the trapped orbits that are the backbone supporting the bar feature. We propose a direct observational technique well suited to integral field unit spectroscopy to measure it. Identifying these orbits and using the dynamical length is a more faithful tracer of the secular evolution and influence of the bar. We demonstrate the success of our observational technique for recovering the maximal bar-parenting orbit in a range of simulations, and to show its promise we perform its measurement on a real galaxy. We also study the difference between traditionally used ellipse fit approaches to determine bar length and the dynamical length proposed here in a wide range of bar-forming N-body simulations of a stellar disc and dark matter halo. We find that ellipse fitting may severely overestimate measurements of the bar length by a factor of 1.5–2.5 relative to the extent of the orbits that are trapped and actually comprise the bar. This bias leads to overestimates of both bar mass and the ratio of corotation radius to bar length, i.e. the bar speed, affecting inferences about the evolution of bars in the real Universe.

The vertical structure of galactic discs: non-local gravity versus dark matter

ABSTRACT Recent isolated galactic simulations show that the morphology of galactic discs in modified gravity differs from that of the standard dark matter model. In this study, we focused on the vertical structure of galactic discs and compared the bending instability in the vertical direction for both paradigms. To achieve this, we utilized high-resolution N-body simulations to construct two models in a specific non-local gravity theory (NLG) and the standard dark matter model and compared their stability against the bending perturbations. Our numerical results demonstrate that the outer regions of the disc are more susceptible to the instability in NLG, whereas the disc embedded in the dark matter halo is more unstable in the central regions. We then interpret these results based on the dispersion relation of the bending waves. To do so, we presented an analytical study to derive the dispersion relation in NLG. Our numerical results align with the predictions of our analytical models. Consequently, we conclude that the analysis of bending instability in galactic discs offers an explanation for the distinct vertical structures observed in simulated galactic discs under these two theories. These findings represent a significant step towards distinguishing between the modified gravity and dark matter models.

The dynamics of accretion flows near to the innermost stable circular orbit

ABSTRACT Accretion flows are fundamentally turbulent systems, yet are classically modelled with viscous theories only valid on length scales significantly greater than the typical size of turbulent eddies in the flow. We demonstrate that, while this will be a reasonable bulk description of the flow at large radii, this must break down as the flow approaches absorbing boundaries, such as the innermost stable circular orbit (ISCO) of a black hole disc. This is because in a turbulent flow large velocity fluctuations can carry a fluid element over the ISCO from a finite distance away, from which it will not return, a process without analogy in conventional models. This introduces a non-zero directional bias into the velocity fluctuations in the near-ISCO disc. By studying reduced random walk problems, we derive a number of implications of the presence of an absorbing boundary in an accretion context. In particular, we show that the average velocity with which a typical fluid element crosses the ISCO is much larger than is assumed in traditional theories. This enhanced velocity modifies the thermodynamic properties of black hole accretion flows on both sides of the ISCO. In particular, thermodynamic quantities for larger ISCO stresses no longer display pronounced cusps at the ISCO in this new formalism, a result with relevance for a number of observational probes of the intra-ISCO region. Finally, we demonstrate that these extended models reproduce the trans-ISCO behaviour observed in GRMHD simulations of thin discs.

Cosmic-Ray Acceleration of Galactic Outflows in Multiphase Gas

We investigate the dynamical interaction between cosmic rays (CRs) and the multiphase interstellar medium (ISM) using numerical magnetohydrodynamic (MHD) simulations with a two-moment CR solver and TIGRESS simulations of star-forming galactic disks. We previously studied the transport of CRs within TIGRESS outputs using a “postprocessing” approach, and we now assess the effects of the MHD backreaction to CR pressure. We confirm our previous conclusion that there are three quite different regimes of CR transport in multiphase ISM gas, while also finding that simulations with “live MHD” predict a smoother CR pressure distribution. The CR pressure near the midplane is comparable to other pressure components in the gas, but the scale height of CRs is far larger. Next, with a goal of understanding the role of CRs in driving galactic outflows, we conduct a set of controlled simulations of the extraplanar region above z = 500 pc, with imposed boundary conditions flowing from the midplane into this region. We explore a range of thermal and kinematic properties for the injected thermal gas, encompassing both hot, fast-moving outflows, and cooler, slower-moving outflows. The boundary conditions for CR energy density and flux are scaled from the supernova rate in the underlying TIGRESS model. Our simulations reveal that CRs efficiently accelerate extraplanar material if the latter is mostly warm/warm-hot gas, in which CRs stream at the Alfvén speed, and the effective sound speed increases as density decreases. In contrast, CRs have very little effect on fast, hot outflows where the Alfvén speed is small, even when the injected CR momentum flux exceeds the injected MHD momentum flux.

Kinematic signatures of planet–disk interactions in vertical shear instability-turbulent protoplanetary disks

Context. Planets are thought to form inside weakly ionized regions of protoplanetary disks, where turbulence creates ideal conditions for solid growth. However, the nature of this turbulence is still uncertain. In fast cooling parts of this zone the vertical shear instability (VSI) can operate, inducing a low level of gas turbulence and large-scale gas motions. Resolving the kinematic signatures of active VSI could reveal the origin of turbulence in planet-forming disk regions. However, an exploration of kinematic signatures of the interplay between VSI and forming planets is needed for a correct interpretation of radio interferometric observations. A robust detection of VSI would lead the way to a deeper understanding of the impact of gas turbulence on planet formation. Aims. The objective of this study is to explore the effect of VSI on the disk substructures triggered by an embedded fairly massive planet. We focus on the impact of this interplay on CO kinematic observations with the ALMA interferometer. Methods. We conducted global 3D hydrodynamical simulations of VSI-unstable disks with and without embedded massive planets, exploring Saturn- and Jupiter-mass cases. We studied the effect of planets on the VSI gas dynamics, and made a comparison with viscous disks. Post-processing the simulations with a radiative transfer code, we examined the kinematic signatures expected in CO molecular line emission, varying disk inclination. Further, we simulated deep ALMA high-resolution observations of our synthetic images, to test the observability of VSI and planetary signatures. Results. The embedded planet produces a damping of the VSI along a radial region, most effective at the disk midplane. For the Saturn case, the VSI modes are distorted by the planet’s spirals producing mixed kinematic signatures. For the Jupiter case, the planet’s influence dominates the overall disk gas kinematics. Conclusions. The presence of massive planets embedded in the disk can weaken the VSI large-scale gas flows, limiting its observability in CO kinematic observations with ALMA.

Emergent Nucleosynthesis from a 1.2 s Long Simulation of a Black Hole Accretion Disk

We simulate a black hole accretion disk system with full-transport general relativistic neutrino radiation magnetohydrodynamics for 1.2 s. This system is likely to form after the merger of two compact objects and is thought to be a robust site of r-process nucleosynthesis. We consider the case of a black hole accretion disk arising from the merger of two neutron stars. Our simulation time coincides with the nucleosynthesis timescale of the r-process (∼1 s). Because these simulations are time-consuming, it is common practice to run for a “short” duration of approximately 0.1–0.3 s. We analyze the nucleosynthetic outflow from this system and compare the results of stopping at 0.12 and 1.2 s. We find that the addition of mass ejected in the longer simulation as well as more favorable thermodynamic conditions from emergent viscous ejecta greatly impacts the nucleosynthetic outcome. We quantify the error in nucleosynthetic outcomes between short and long cuts.

Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks

Context. The vertical shear instability (VSI) is a hydrodynamical instability predicted to produce turbulence in magnetically inactive regions of protoplanetary disks. The regions in which this instability can occur and the physical phenomena leading to its saturation are a current matter of research. Aims. We explore the secondary instabilities triggered by the nonlinear evolution of the VSI and their role in its saturation. We also expand on previous investigations on stability regions by considering temperature stratifications enforced by stellar irradiation and radiative cooling, and including the effects of dust-gas collisions and molecular line emission. Methods. We modeled the gas-dust mixture in a circumstellar disk around a T Tauri star by means of high-resolution axisymmetric radiation-hydrodynamical simulations including stellar irradiation with frequency-dependent opacities, considering different degrees of depletion of small dust grains. Results. The flow pattern produced by the interplay of the axisymmetric VSI modes and the baroclinic torque forms bands of nearly uniform specific angular momentum. In the high-shear regions in between these bands, the Kelvin–Helmholtz instability (KHI) is triggered. A third instability mechanism, consisting of an amplification of eddies by baroclinic torques, forms meridional vortices with Mach numbers up to ∼0.4. Our stability analysis suggests that protoplanetary disks can be VSI-unstable in surface layers up to tens of au for reasonably high gas emissivities. Conclusions. The significant transfer of kinetic energy to small-scale eddies produced by the KHI and possibly even the baroclinic acceleration of eddies limit the maximum energy of the VSI modes, likely leading to the saturation of the VSI. Depending on the gas molecular composition, the VSI can operate at surface layers even in regions where the midplane is stable. This picture is consistent with current observations of disks showing thin midplane millimeter-sized dust layers while appearing vertically extended in optical and near-infrared wavelengths.