Simulations Of Disks Research Articles

How do the successive buckling events affect a galaxy bar and stellar disc? Potential observable signatures for spotting the buckling action – I

ABSTRACT Until now, observations have caught up only a handful of galaxies in ongoing buckling action. Interestingly, N-body simulations over the past several decades show that almost every bar buckles or vertically thickens as soon as it reaches its peak strength during its evolution and leads to box/peanut/x (BPX) shapes. In order to understand the effect of multiple buckling events on the observable properties of galactic bar and disc, we perform an N-body simulation of a Milky Way-type disc. The axisymmetric galaxy disc forms a bar within a Gyr of its evolution and the bar undergoes two successive buckling events. We report that the time-spans of these two buckling events are 220 Myr and 1 Gyr, which have almost similar strengths of the bending modes. As a result of these two buckling events, the full lengths of BPX shapes are around 5.8 and 8.6 kpc, which are around two-thirds of the full bar length at the end of each buckling event. We find that the first buckling occurs at a smaller scale (radius $\lt $3 kpc) with a shorter time-span affecting the larger length-scales of the disc, which is quantified in terms of changes in $m=$2 and $m=$ 4 Fourier modes. While the second buckling occurs at larger scales (radius $\approx$6 kpc) affecting the inner disc the most. Finally, we provide observable kinematic signatures (i.e. quadrupolar patterns of the line-of-sight velocities), which can potentially differentiate the successive buckling events.

Shock-induced Partial Alignment in Geometrically Thick Tilted Accretion Disks Around Black Holes

We carry out idealized three-dimensional general-relativistic magnetohydrodynamic simulations of prograde, weakly magnetized, and geometrically thick accretion flows where the gas distribution is misaligned from the black hole (BH) spin axis. We evolve the disk for three BH spins: a = 0.5, 0.75, and 0.9375, and we contrast them with a standard aligned disk simulation with a = 0.9375. The tilted disks achieve a warped and twisted steady-state structure, with the outer disk misaligning further away from the BH and surpassing the initial 24° misalignment. However, closer to the BH, there is evidence of partial alignment, as the inclination angle decreases with radius in this regime. Standing shocks also emerged in proximity to the BH, roughly at ∼6 gravitational radii. We show that these shocks act to partially align the inner disk with the BH spin. The rate of alignment increases with increasing BH spin magnitude, but in all cases is insufficient to fully align the gas before it accretes. Additionally, we present a toy model of orbit crowding that can predict the location of the shocks in moderate-to-fast rotating BHs, illustrating a potential physical origin for the behavior seen in simulations—with possible applications in determining the positions of shocks in real misaligned astrophysical systems.

The origin of the H α line profiles in simulated disc galaxies

ABSTRACT Observations of ionized H $\alpha$ gas in high-redshift disc galaxies have ubiquitously found significant line broadening, $\sigma _{\rm H\,\alpha }\sim 10{\!-\!}100\, {\rm km\, s^{-1}}$. To understand whether this broadening reflects gas turbulence within the interstellar medium (ISM) of galactic discs, or arises from out-of-plane emission in mass-loaded outflows, we perform radiation hydrodynamic simulations of isolated Milky Way-mass disc galaxies in a gas-poor (low-redshift) and gas rich (high-redshift) condition and create mock H $\alpha$ emission line profiles. We find that the majority of the total (integrated) ${\rm H\,\alpha }$ emission is confined within the ISM, with extraplanar gas contributing ${\sim} 45~{{\ \rm per\ cent}}$ of the extended profile wings ($v_z\ge 200$${\, \rm {km\, s^{-1}} }$) in the gas-rich galaxy. This substantiates using the ${\rm H\,\alpha }$ emission line as a tracer of mid-plane disc dynamics. We investigate the relative contribution of diffuse and dense ${\rm H\,\alpha }$ emitting gas, corresponding to diffuse ionized gas (DIG; $\rho \lesssim 0.1\, {\rm cm^{-3}}$, $T\sim 8\, 000$ K) and H ii regions ($\rho \gtrsim 10\, {\rm cm^{-3}}$, $T\sim 10\, 000$ K), respectively, and find that DIG contributes $f_{\rm DIG}\lesssim 10~{{\ \rm per\ cent}}$ of the total L$_{\rm H\alpha }$. However, the DIG can reach upwards of $\sigma _{\rm H\,\alpha } \sim 60{\!-\!}80\, {\rm km\, s^{-1}}$ while the H ii regions are much less turbulent $\sigma _{\rm H\,\alpha }\sim 10{\!-\!}40\, {\rm km\, s^{-1}}$. This implies that the $\sigma _{\rm H\,\alpha }$ observed using the full ${\rm H\,\alpha }$ emission line is dependent on the relative ${\rm H\,\alpha }$ contribution from DIG/H ii regions and a larger $f_{\rm DIG}$ would shift $\sigma _{\rm H\,\alpha }$ to higher values. Finally, we show that $\sigma _{\rm H\,\alpha }$ evolves, in both the DIG and H ii regions, with the galaxy gas fraction. Our high-redshift equivalent galaxy is roughly twice as turbulent, except for in the DIG which has a more shallow evolution.

The active role of co-evolving haloes in stellar bar formation

ABSTRACT We use idealized N-body simulations of equilibrium discs in live and static haloes to study how dark matter co-evolution impacts the assembly of stellar particles into a bar and the halo response. Initial conditions correspond to a marginally unstable disc according to commonly used disc stability criteria, and are evolved for the equivalent of about 150 disc dynamical times (10 Gyr). An extensive convergence study ensures accurate modelling of the bar formation process. Live haloes lead to the formation of a strong bar, but the same disc remains unbarred when evolved in a static halo. Neither seeded disc instabilities nor longer (60 Gyr) simulations result in the formation of a bar when the halo is static. When the live halo is replaced with a static analogue at later times, the previously robust bar slowly dissipates, suggesting that (1) the co-evolution of the disc and halo is critical for the assembly and long-term survival of bars in marginally unstable discs and (2) global disc stability criteria must be modified for discs in the presence of live haloes. In our live halo runs, a ‘dark bar’ grows synchronously with the stellar bar. Processes that inhibit the transfer of angular momentum between the halo and disc may stabilize a galaxy against bar formation, and can lead to the dissolution of the bar itself. This raises further questions about the puzzling stability of observed discs that are marginally unstable, but unbarred.

Two-dimensional simulations of disks in close binaries

Context. Previous simulations of cataclysmic variables studied either the quiescence, or the outburst state in multiple dimensions or they simulated complete outburst cycles in one dimension using simplified models for the gravitational torques. Aims. We self-consistently simulate complete outburst cycles of normal and superoutbursts in cataclysmic variable systems in two dimensions. We study the effect of different α viscosity parameters, mass transfer rates, and binary mass ratios on the disk luminosities, outburst occurrence rates, and superhumps. Methods. We simulate non-isothermal, viscous accretion disks in cataclysmic variable systems using a modified version of the FARGO code with an updated equation of state and a cooling function designed to reproduce s-curve behavior. Results. Our simulations can model complete outburst cycles using the thermal tidal instability model. We find higher superhump amplitudes and stronger gravitational torques than previous studies, resulting in better agreement with observations.

Angular momentum transfer in cosmological simulations of Milky Way-mass discs

Angular momentum transfer in cosmological simulations of Milky Way-mass discs

Resolution analysis of magnetically arrested disc simulations

ABSTRACT Polarization measurements by the Event Horizon Telescope from M87* and Sgr A* suggest that there is a dynamically strong, ordered magnetic field, typical of what is expected of a magnetically arrested accretion disc (MAD). In such discs, the strong poloidal magnetic field can suppress the accretion flow and cause episodic flux eruptions. Recent work shows that general relativistic magnetohydrodynamic (GRMHD) MAD simulations feature dynamics of turbulence and mixing instabilities that are becoming resolved at higher resolutions. We perform a convergence study of MAD states exceeding the status quo by an order of magnitude in resolution. We use existing 3D simulations performed with the H-AMR code, up to a resolution of 5376 × 2304 × 2304 in a logarithmic spherical-polar grid. We find consistent time-averaged disc properties across all resolutions. However, higher resolutions reveal signs of inward angular momentum transport attributed to turbulent convection, particularly evident when mixing instabilities occur at the surfaces of flux tubes during flux eruptions. Additionally, we see wave-like features in the jet sheath, which become more prominent at higher resolutions, that may induce mixing between the jet and disc. At higher resolutions, we observe the sheath to be thinner, resulting in increased temperature, reduced magnetization, and greater variability. Those differences could affect the dissipation of energy, which would eventually result in distinct observable radiative emission from high-resolution simulations. With higher resolutions, we can delve into crucial questions about horizon-scale physics and its impact on the dynamics and emission properties of larger-scale jets.

Warps and breaks in circumbinary discs

ABSTRACT Disc warping, and possibly disc breaking, has been observed in protoplanetary discs around both single and multiple stars. Large warps can break the disc, producing multiple observational signatures. In this work, we use comparisons of disc time-scales to derive updated formulae for disc breaking, with better predictions as to when and where a disc is expected to break and how many breaks could occur. Disc breaking is more likely for discs with small inner cavities, cooler temperatures, and steeper power-law profiles, such that thin, polar-aligning discs are more likely to break. We test our analytical formulae using three-dimensional grid-based simulations of protoplanetary discs warped by the gravitational torque of an inner binary. We reproduce the expected warp behaviours in different viscosity regimes and observe disc breaking at locations in agreement with our derived equations. As our simulations only show disc breaking when disc viscosity is low, we also consider a viscous criterion for disc breaking, where rapid alignment to the precession vector can prevent a break by reducing the maximum misalignment between neighbouring rings. We apply these results to the GW Orionis circumtriple disc and find that the precession induced from the central stars can break the disc if it is relatively thin. We expect repeated or multiple disc breaking to occur for discs with sufficiently steep power-law profiles. We simulate a polar-aligning disc around an eccentric binary with steep power-law profiles and observe two separate breaking events at locations in rough agreement with our analytical predictions.

Exploring Viscosity Effects in Accretion Disk Simulations: A Code Comparison between FLASH and Athena

This project investigates the impact of viscosity on the evolution of black hole accretion disks through a comparative analysis between FLASH and Athena simulations. Our focus lies in understanding the role of viscosity in overall disk evolution. We present mass accretion rate plots that reveal near identical patterns in disk evolution when varying viscosity. We emphasize the need for additional code-comparison research, including exploration of a wider variety of codes and physics included to enhance our understanding and confidence in simulated models.

Inner dusty regions of protoplanetary discs – III. The role of non-radial radiation pressure in dust dynamics

ABSTRACT We explore the dynamical behaviour of dust particles that populate the surface of inner optically thick protoplanetary discs. This is a disc region with the hottest dust and is of a great importance for planet formation and dust evolution, but we still struggle to understand all the forces that shape this environment. In our approach, we combine results from two separate numerical studies, one is the wind velocity and density distributions obtained from magnetohydrodynamical simulations of accretion discs, and the other is a high-resolution multigrain dust radiation transfer. In our previous paper in the series, we described the methodology for utilizing these results as an environmental input for the integration of dust trajectories driven by gravity, gas drag, and radiation pressure. Now we have two improvements, we incorporate time changes in the wind density and velocity, and we implement the non-radial radiation pressure force. We applied our analysis on the Herbig Ae and T Tau stars. We confirm that the radiation pressure force can lead to dust outflow, especially in the case of more luminous stars. Additionally, it opposes dust accretion at the inner disc edge and reduces dust settling. These effects are enhanced by the disc wind, especially in the zone where the stellar and the disc magnetic fields meet. Our results suggest that dust grains can stay in the hottest disc region for an extended period and then end up ejected into the outer disc regions.

Rapid Protoplanet Formation in Vortices: Three-dimensional Local Simulations with Self-gravity

Disk vortices, seen in numerical simulations of protoplanetary disks and found observationally in Atacama Large Millimeter/submillimeter Array and Very Large Array images of these objects, are promising sites for planet formation given their pebble trapping abilities. Previous works have shown a strong concentration of pebbles in vortices, but gravitational collapse has only been shown in low-resolution, two-dimensional, global models. In this Letter, we aim to study the pebble concentration and gravitational collapse of pebble clouds in vortices via high-resolution, three-dimensional, local models. We performed simulations of the dynamics of gas and solids in a local shearing box where the gas is subject to convective overstability, generating a persistent giant vortex. We find that the vortex produces objects of Moon and Mars mass, with a mass function of power-law dlnN/dlnM=−1.6±0.3 . The protoplanets grow rapidly, doubling in mass in about five orbits, following pebble accretion rates. The mass range and mass doubling rate are in broad agreement with previous low-resolution global models. We conclude that Mars-mass planetary embryos are the natural outcome of planet formation inside the disk vortices seen in millimeter and radio images of protoplanetary disks.

Galaxies with grains: unraveling dust evolution and extinction curves with hydrodynamical simulations

Dust in galaxies is an important tracer of galaxy properties and their evolution over time. The physical origin of the grain size distribution, the dust chemical composition, and, hence, the associated ultraviolet-to-optical extinctions in diverse galaxies remains elusive. To address this issue, we introduce a model for dust evolution in the RAMSES code for simulations of galaxies with a resolved multiphase interstellar medium. Dust is modelled as a fluid transported with the gas component, and is decomposed into two sizes, 5 nm and 0.1 μm, and two chemical compositions for carbonaceous and silicate grains. This dust model includes the growth of dust by accretion of elements from the gas phase and by the release of dust in stellar ejecta, the destruction by thermal sputtering, supernovae, and astration, and the exchange of dust mass between the two main populations of grain sizes by coagulation and shattering. Using a suite of isolated disc simulations with different masses and metallicities, the simulations can explore the role of these processes in shaping the key properties of dust in galaxies. The simulated Milky Way analogue reproduces the dust-to-metal mass ratio, depletion factors, size distribution and extinction curves of the Milky Way. Galaxies with lower metallicities reproduce the observed decrease in the dust-to-metal mass ratio with metallicity at around a few 0.1 Z⊙. This break in the dust-to-metal ratio corresponds to a galactic gas metallicity threshold that marks the transition from an ejecta-dominated to an accretion-dominated grain growth, and that is different for silicate and carbonaceous grains, with ≃0.1 Z⊙ and ≃0.5 Z⊙ respectively. This leads to more Magellanic Cloud-like extinction curves, i.e. with steeper slopes in the ultraviolet and a weaker bump feature at 2175 Å, in galaxies with lower masses and lower metallicities. Steeper slopes in these galaxies are caused by the combination of the higher efficiency of gas accretion by silicate relative to carbonaceous grains and by the low rates of coagulation that preserves the amount of small silicate grains. Weak bumps are due to the overall inefficient accretion growth of carbonaceous dust at low metallicity, whose growth is mostly supported by the release of large grains in SN ejecta. We also show that the formation of CO molecules is a key component to limit the ability of carbonaceous dust to grow, in particular in low-metallicity gas-rich galaxies.

Magnetorotational dynamo can generate large-scale vertical magnetic fields in 3D GRMHD simulations of accreting black holes

ABSTRACT Jetted astrophysical phenomena with black hole engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed three-dimensional global general relativistic magnetohydrodynamical simulations of accretion discs to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the non-linear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their non-linear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disc. This leads to a larger characteristic size of the, now advected, magnetic field when compared to the local disc height. We characterize the complete dynamo mechanism with two time-scales: one for the local magnetic field generation, $t_{\rm gen}$, and one for the large-scale scale advection, $t_{\rm adv}$. Whereas the dynamo we describe is non-linear, we explore the potential of linear mean field models to replicate its core features. Our findings indicate that traditional $\alpha$-dynamo models, often computed in stratified shearing box simulations, are inadequate and that the effective large-scale dynamics is better described by the shear current effects or stochastic $\alpha$-dynamos.

Vertical shear instability with dust evolution and consistent cooling times

Context. Gas in protoplanetary disks mostly cools via thermal accommodation with dust particles. Thermal relaxation is thus highly sensitive to the local dust size distributions and the spatial distribution of the grains. So far, the interplay between thermal relaxation and gas turbulence has not been dynamically modeled in hydrodynamic simulations of protoplanetary disks with dust. Aims. We aim to study the effects of the vertical shear instability (VSI) on the thermal relaxation times, and vice versa. We are particularly interested in the influence of the initial dust grain size on the VSI and whether the emerging turbulence is sustained over long timescales. Methods. We ran three axisymmetric hydrodynamic simulations of a protoplanetary disk including four dust fluids that initially resemble MRN size distributions of different initial grain sizes. From the local dust densities, we calculated the thermal accommodation timescale of dust and gas and used the result as the thermal relaxation time of the gas in our simulation. We included the effect of dust growth by applying the monodisperse dust growth rate and the typical growth limits. Results. We find that the emergence of the VSI is strongly dependent on the initial dust grain size. Coagulation also counteracts the emergence of hydrodynamic turbulence in our simulations, as shown by others before. Starting a simulation with larger grains (100 μm) generally leads to a less turbulent outcome. While the inner disk regions (within ∼70 au) develop turbulence in all three simulations, we find that the simulations with larger particles do not develop VSI in the outer disk. Conclusions. Our simulations with dynamically calculated thermal accommodation times based on the drifting and settling dust distribution show that the VSI, once developed in a disk, can be sustained over long timescales, even if grain growth is occurring. The VSI corrugates the dust layer and even diffuses the smaller grains into the upper atmosphere, where they can cool the gas. Whether the instability can emerge for a specific stratification depends on the initial dust grain sizes and the initial dust scale height. If the grains are initially ≳100 μm and if the level of turbulence is initially assumed to be low, we find no VSI turbulence in the outer disk regions.

The backreaction of stellar wobbling on accretion discs of massive protostars

black In recent years, it has been demonstrated that massive stars see their infant circumstellar medium shaped into a large irradiated, gravitationally unstable accretion disc during their early formation phase. Such discs constitute the gas reservoir from which nascent high-mass stars gain a substantial fraction of their mass by episodic accretion of dense gaseous circumstellar clumps, simultaneously undergoing accretion-driven bursts and producing close-orbit spectroscopic companions of the young high-mass stellar object. We aim to evaluate the effects of stellar motion caused by the disc non-axisymmetric gravitational field on the disc evolution and its spatial morphology. In particular, we analyse the disc's propensity to gravitational instability and fragmentation and the disc's appearance in synthetic millimetre band images pertinent to the ALMA facility. We employed three-dimensional radiation-hydrodynamical simulations of the surroundings of a young massive star in the non-inertial spherical coordinate system, adopting the highest spatial resolution to date and including the indirect star-disc gravitational potential caused by the asymmetries in the circumstellar disc. The resulting disc configurations were post-processed with the radiation transfer tool RADMC-3D and CASA software to obtain synthetic images of the disc. We confirm that the early evolution of the accretion disc is notably different when stellar wobbling is taken into account. The redistribution of angular momentum in the system makes the disc smaller and rounder, reduces the number of circumstellar gaseous clumps formed via disc gravitational fragmentation, and prevents the ejection of gaseous clumps from the disc. The synthetic predictive images at millimetre wavelengths of the accretion disc that includes stellar wobbling are in better agreement with the observations of the surroundings of massive young stellar objects, namely black AFGL 4176 mm1, G17.64+0.16, and G353.273 than our simulations of numerical hydrodynamics that omit this physical mechanism. Our work confirms that stellar wobbling is an essential ingredient to account for in numerical simulations of accretion discs of massive protostars.

Dynamics near the inner dead-zone edges in a proprotoplanetary disk

Abstract We perform three-dimensional global non-ideal magnetohydrodynamic simulations of a protoplanetary disk containing the inner dead-zone edge. We take into account realistic diffusion coefficients of the Ohmic resistivity and ambipolar diffusion based on detailed chemical reactions with single-size dust grains. We found that the conventional dead zone identified by the Elsässer numbers of the Ohmic resistivity and ambipolar diffusion is divided into two regions: “the transition zone” and “the coherent zone.” The coherent zone has the same properties as the conventional dead zone, and extends outside of the transition zone in the radial direction. Between the active and coherent zones, we discover the transition zone, the inner edge of which is identical to that of the conventional dead zone. The transition zone extends out over the regions where thermal ionization determines diffusion coefficients. The transition zone has completely different physical properties than the conventional dead zone, the so-called undead zone, and the zombie zone. The combination of amplification of the radial magnetic field owing to the ambipolar diffusion and a steep radial gradient of the Ohmic diffusivity causes the efficient evacuation of the net vertical magnetic flux from the transition zone within several rotations. Surface gas accretion occurs in the coherent zone but not in the transition zone. The presence of the transition zone prohibits mass and magnetic flux transport from the coherent zone to the active zone. Mass accumulation occurs at both edges of the transition zone as a result of mass supply from the active and coherent zones.

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.