BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Chao-Chin Yang

Chao-Chin Yang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.COastro-ph.GAastro-ph.SRastro-ph.IMMachine Learningphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

14works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2024arXiv

Particle clustering in turbulence: Prediction of spatial and statistical properties with deep learning

We investigate the utility of deep learning for modeling the clustering of particles that are aerodynamically coupled to turbulent fluids. Using a Lagrangian particle module within the Athena++ hydrodynamics code, we simulate the dynamics of particles in the Epstein drag regime within a periodic domain of isotropic forced hydrodynamic turbulence. This setup is an idealized model relevant to the collisional growth of micron to mm-sized dust particles in early stage planet formation. The simulation data are used to train a U-Net deep learning model to predict gridded three-dimensional representations of the particle density and velocity fields, given as input the corresponding fluid fields. The trained model qualitatively captures the filamentary structure of clustered particles in a highly non-linear regime. We assess model fidelity by calculating metrics of the density field (the radial distribution function) and of the velocity field (the relative velocity and the relative radial velocity between particles). Although trained only on the spatial fields, the model predicts these statistical quantities with errors that are typically <10%. Our results suggest that, given appropriately expanded training data, deep learning could complement direct numerical simulations in predicting particle clustering within turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Direct Formation of Planetary Embryos in Self-Gravitating Disks

Giant planets have been discovered at large separations from the central star. Moreover, a striking number of young circumstellar disks have gas and/or dust gaps at large orbital separations, potentially driven by embedded planetary objects. To form massive planets at large orbital separations through core accretion within disk lifetime, however, an early solid body to seed pebble and gas accretion is desirable. Young protoplanetary disks are likely self-gravitating, and these gravitoturbulent disks may efficiently concentrate solid material at the midplane driven by spiral waves. We run 3D local hydrodynamical simulations of gravitoturbulent disks with Lagrangian dust particles to determine whether particle and gas self-gravity can lead to the formation of dense solid bodies, seeding later planet formation. When self-gravity between dust particles is included, solids of size $\mathrm{St} = 0.1$ to $1$ concentrate within the gravitoturbulent spiral features and collapse under their own self-gravity into dense clumps up to several $M_{\oplus}$ in mass at wide orbits. Simulations with dust that drift most efficiently, $\mathrm{St}=1$, form the most massive clouds of particles, while simulations with smaller dust particles, $\mathrm{St}=0.1$, have clumps with masses an order of magnitude lower. When the effect of dust backreaction onto the gas is included, dust clumps become smaller by a factor of a few but more numerous. The existence of large solid bodies at an early stage of the disk can accelerate the planet formation process, particularly at wide orbital separations, and potentially explain planets distant from the central stars and young protoplanetary disks with substructures.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Formation of Dust Rings and Gaps in Non-ideal MHD Disks Through Meridional Gas Flows

Rings and gaps are commonly observed in the dust continuum emission of young stellar disks. Previous studies have shown that substructures naturally develop in the weakly ionized gas of magnetized, non-ideal MHD disks. The gas rings are expected to trap large mm/cm-sized grains through pressure gradient-induced radial dust-gas drift. Using 2D (axisymmetric) MHD simulations that include ambipolar diffusion and dust grains of three representative sizes (1~mm, 3.3~mm, and 1~cm), we show that the grains indeed tend to drift radially relative to the gas towards the centers of the gas rings, at speeds much higher than in a smooth disk because of steeper pressure gradients. However, their spatial distribution is primarily controlled by meridional gas motions, which are typically much faster than the dust-gas drift. In particular, the grains that have settled near the midplane are carried rapidly inwards by a fast accretion stream to the inner edges of the gas rings, where they are lifted up by the gas flows diverted away from the midplane by a strong poloidal magnetic field. The flow pattern in our simulation provides an attractive explanation for the meridional flows recently inferred in HD 163296 and other disks, including both "collapsing" regions where the gas near the disk surface converges towards the midplane and a disk wind. Our study highlights the prevalence of the potentially observable meridional flows associated with the gas substructure formation in non-ideal MHD disks and their crucial role in generating rings and gaps in dust.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Kozai-Lidov oscillations triggered by a tilt instability of detached circumplanetary discs

Circumplanetary discs can be linearly unstable to the growth of disc tilt in the tidal potential of the star-planet system. We use three-dimensional hydrodynamical simulations to characterize the disc conditions needed for instability, together with its long term evolution. Tilt growth occurs for disc aspect ratios, evaluated near the disc outer edge, of $H/r\gtrsim 0.05$, with a weak dependence on viscosity in the wave-like regime of warp propagation. Lower mass giant planets are more likely to have circumplanetary discs that satisfy the conditions for instability. We show that the tilt instability can excite the inclination to above the threshold where the circumplanetary disc becomes unstable to Kozai--Lidov (KL) oscillations. Dissipation in the Kozai--Lidov unstable regime caps further tilt growth, but the disc experiences large oscillations in both inclination and eccentricity. Planetary accretion occurs in episodic accretion events. We discuss implications of the joint tilt--KL instability for the detectability of circumplanetary discs, for the obliquity evolution of forming giant planets, and for the formation of satellite systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Formation of the polar debris disc around 99 Herculis

We investigate the formation mechanism for the observed nearly polar aligned (perpendicular to the binary orbital plane) debris ring around the eccentric orbit binary 99 Herculis. An initially inclined nonpolar debris ring or disc will not remain flat and will not evolve to a polar configuration, due to the effects of differential nodal precession that alter its flat structure. However, a gas disc with embedded well coupled solids around the eccentric binary may evolve to a polar configuration as a result of pressure forces that maintain the disc flatness and as a result of viscous dissipation that allows the disc to increase its tilt. Once the gas disc disperses, the debris disc is in a polar aligned state in which there is little precession. We use three-dimensional hydrodynamical simulations, linear theory, and particle dynamics to study the evolution of a misaligned circumbinary gas disc and explore the effects of the initial disc tilt, mass, and size. We find that for a wide range of parameter space, the polar alignment timescale is shorter than the lifetime of the gas disc. Using the observed level of alignment of 3 deg. from polar, we place an upper limit on the mass of the gas disc of about 0.014 M_sun at the time of gas dispersal. We conclude that the polar debris disc around 99 Her can be explained as the result of an initially moderately inclined gas disc with embedded solids. Such a disc may provide an environment for the formation of polar planets.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Morphological signatures induced by dust back reaction in discs with an embedded planet

Recent observations have revealed a gallery of substructures in the dust component of nearby protoplanetary discs, including rings, gaps, spiral arms, and lopsided concentrations. One interpretation of these substructures is the existence of embedded planets. Not until recently, however, most of the modelling effort to interpret these observations ignored the dust back reaction to the gas. In this work, we conduct local-shearing-sheet simulations for an isothermal, inviscid, non-self-gravitating, razor-thin dusty disc with a planet on a fixed circular orbit. We systematically examine the parameter space spanned by planet mass ($0.1M_\text{th} \le M_\text{p} \le 1M_\text{th}$, where $M_\text{th}$ is the thermal mass), dimensionless stopping time ($10^{-3} \le τ_\text{s} \le 1$), and solid abundance ($0 < Z \le 1$). We find that when the dust particles are tightly coupled to the gas ($τ_\text{s} < 0.1$), the spiral arms are less open and the gap driven by the planet becomes deeper with increasing $Z$, consistent with a reduced speed of sound in the approximation of a single dust-gas mixture. By contrast, when the dust particles are marginally coupled ($0.1 \lesssim τ_\text{s} \lesssim 1$), the spiral structure is insensitive to $Z$ and the gap structure in the gas can become significantly skewed and unidentifiable. When the latter occurs, the pressure maximum radially outside of the planet is weakened or even extinguished, and hence dust filtration by a low-mass ($M_\text{p} < M_\text{th}$) planet could be reduced or eliminated. Finally, we find that the gap edges where the dust particles are accumulated as well as the lopsided large-scale vortices driven by a massive planet, if any, are unstable, and they are broken into numerous small-scale dust-gas vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Initial mass function of planetesimals formed by the streaming instability

The streaming instability is a mechanism to concentrate solid particles into overdense filaments that undergo gravitational collapse and form planetesimals. However, it remains unclear how the initial mass function of these planetesimals depends on the box dimensions of numerical simulations. To resolve this, we perform simulations of planetesimal formation with the largest box dimensions to date, allowing planetesimals to form simultaneously in multiple filaments that can only emerge within such large simulation boxes. In our simulations, planetesimals with sizes between 80 km and several hundred kilometers form. We find that a power law with a rather shallow exponential cutoff at the high-mass end represents the cumulative birth mass function better than an integrated power law. The steepness of the exponential cutoff is largely independent of box dimensions and resolution, while the exponent of the power law is not constrained at the resolutions we employ. Moreover, we find that the characteristic mass scale of the exponential cutoff correlates with the mass budget in each filament. Together with previous studies of high-resolution simulations with small box domains, our results therefore imply that the cumulative birth mass function of planetesimals is consistent with an exponentially tapered power law with a power-law exponent of approximately -1.6 and a steepness of the exponential cutoff in the range of 0.3-0.4.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Integration of Particle-Gas Systems with Stiff Mutual Drag Interaction

Numerical simulation of numerous mm/cm-sized particles embedded in a gaseous disk has become an important tool in the study of planet formation and in understanding the dust distribution in observed protoplanetary disks. However, the mutual drag force between the gas and the particles can become so stiff, particularly because of small particles and/or strong local solid concentration, that an explicit integration of this system is computationally formidable. In this work, we consider the integration of the mutual drag force in a system of Eulerian gas and Lagrangian solid particles. Despite the entanglement between the gas and the particles under the particle-mesh construct, we are able to devise a numerical algorithm that effectively decomposes the globally coupled system of equations for the mutual drag force and makes it possible to integrate this system on a cell-by-cell basis, which considerably reduces the computational task required. We use an analytical solution for the temporal evolution of each cell to relieve the time-step constraint posed by the mutual drag force as well as to achieve the highest degree of accuracy. To validate our algorithm, we use an extensive suite of benchmarks with known solutions in one, two, and three dimensions, including the linear growth and the nonlinear saturation of the streaming instability. We demonstrate numerical convergence and satisfactory consistency in all cases. Our algorithm can for example be applied to model the evolution of the streaming instability with mm/cm-sized pebbles at high mass loading, which has important consequences for the formation scenarios of planetesimals.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Feeding Zone of Planetesimal Formation by the Streaming Instability

The streaming instability is a promising mechanism to overcome the barriers in direct dust growth and lead to the formation of planetesimals. Most previous studies of the streaming instability, however, were focused on a local region of a protoplanetary disk with a limited simulation domain such that only one filamentary concentration of solids has been observed. The characteristic separation between filaments is therefore not known. To address this, we conduct the largest-scale simulations of the streaming instability to date, with computational domains up to 1.6 gas scale heights both horizontally and vertically. The large dynamical range allows the effect of vertical gas stratification to become prominent. We observe more frequent merging and splitting of filaments in simulation boxes of high vertical extent. We find multiple filamentary concentrations of solids with an average separation of about 0.2 local gas scale heights, much higher than the most unstable wavelength from linear stability analysis. This measures the characteristic separation of planetesimal forming events driven by the streaming instability and thus the initial feeding zone of planetesimals.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Temperature Fluctuations driven by Magnetorotational Instability in Protoplanetary Disks

The magnetorotational instability (MRI) drives magnetized turbulence in sufficiently ionized regions of protoplanetary disks, leading to mass accretion. The dissipation of the potential energy associated with this accretion determines the thermal structure of accreting regions. Until recently, the heating from the turbulence has only been treated in an azimuthally averaged sense, neglecting local fluctuations. However, magnetized turbulence dissipates its energy intermittently in current sheet structures. We study this intermittent energy dissipation using high resolution numerical models including a treatment of radiative thermal diffusion in an optically thick regime. Our models predict that these turbulent current sheets drive order unity temperature variations even where the MRI is damped strongly by Ohmic resistivity. This implies that the current sheet structures where energy dissipation occurs must be well resolved to correctly capture the flow structure in numerical models. Higher resolutions are required to resolve energy dissipation than to resolve the magnetic field strength or accretion stresses. The temperature variations are large enough to have major consequences for mineral formation in disks, including melting chondrules, remelting calcium-aluminum rich inclusions, and annealing silicates; and may drive hysteresis: current sheets in MRI active regions could be significantly more conductive than the remainder of the disk.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hydrodynamical Simulations of the Barred Spiral Galaxy NGC 1097

NGC 1097 is a nearby barred spiral galaxy believed to be interacting with the elliptical galaxy NGC 1097A located to its northwest. It hosts a Seyfert 1 nucleus surrounded by a circumnuclear starburst ring. Two straight dust lanes connected to the ring extend almost continuously out to the bar. The other ends of the dust lanes attach to two main spiral arms. To provide a physical understanding of its structural and kinematical properties, two-dimensional hydrodynamical simulations have been carried out. Numerical calculations reveal that many features of the gas morphology and kinematics can be reproduced provided that the gas flow is governed by a gravitational potential associated with a slowly rotating strong bar. By including the self-gravity of the gas disk in our calculation, we have found the starburst ring to be gravitationally unstable which is consistent with the observation in \citet{hsieh11}. Our simulations show that the gas inflow rate is 0.17 M$_\sun$ yr$^{-1}$ into the region within the starburst ring even after its formation, leading to the coexistence of both a nuclear ring and a circumnuclear disk.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planetesimal and Protoplanet Dynamics in a Turbulent Protoplanetary Disk: Ideal Stratified Disks

Due to the gravitational influence of density fluctuations driven by magneto-rotational instability in the gas disk, planetesimals and protoplanets undergo diffusive radial migration as well as changes in other orbital properties. The magnitude of the effect on particle orbits can have important consequences for planet formation scenarios. We use the local-shearing-box approximation to simulate an ideal, isothermal, magnetized gas disk with vertical density stratification and simultaneously evolve numerous massless particles moving under the gravitational field of the gas and the host star. We measure the evolution of the particle orbital properties, including mean radius, eccentricity, inclination, and velocity dispersion, and its dependence on the disk properties and the particle initial conditions. Although the results converge with resolution for fixed box dimensions, we find the response of the particles to the gravity of the turbulent gas correlates with the horizontal box size, up to 16 disk scale heights. This correlation indicates that caution should be exercised when interpreting local-shearing-box models involving gravitational physics of magneto-rotational turbulence. Based on heuristic arguments, nevertheless, the criterion L_h / R ~ O(1), where L_h is the horizontal box size and R is the distance to the host star, is proposed to possibly circumvent this conundrum. If this criterion holds, we can still conclude that magneto-rotational turbulence seems likely to be ineffective at driving either diffusive migration or collisional erosion under most circumstances.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thermal-Instability-Driven Turbulent Mixing in Galactic Disks: I. Effective Mixing of Metals

Observations show that radial metallicity gradients in disk galaxies are relatively shallow, if not flat, especially at large galactocentric distances and for galaxies in the high-redshift universe. Given that star formation and metal production are centrally concentrated, this requires a mechanism to redistribute metals. However, the nature of this mechanism is poorly understood, let alone quantified. To address this problem, we conduct magnetohydrodynamical simulations of a local shearing sheet of a thin, thermally unstable, gaseous disk driven by a background stellar spiral potential, including metals modeled as passive scalar fields. Contrary to what a simple α prescription for the gas disk would suggest, we find that turbulence driven by thermal instability is very efficient at mixing metals, regardless of the presence or absence of stellar spiral potentials or magnetic fields. The timescale for homogenizing randomly distributed metals is comparable to or less than the local orbital time in the disk. This implies that turbulent mixing of metals is a significant process in the history of chemical evolution of disk galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Rayleigh Adjustment of Narrow Barriers in Protoplanetary Discs

Sharp density features in protoplanetary discs, for instance at the edge of a magnetically dead zone, have recently been proposed as effective barriers to slow down or even stop the problematically fast migration of planetary cores into their central star. Density features on a radial scale approaching the disc vertical scale height might not exist, however, since they could be Rayleigh (or more generally Solberg-Hoiland) unstable. Stability must be checked explicitly in one-dimensional viscous accretion disc models because these instabilities are artificially eliminated in the process of reducing the full set of axisymmetric equations. The disc thermodynamics, via the entropy stratification, and its vertical structure also influence stability when sharp density features are present. We propose the concept of Rayleigh adjustment for viscous disc models: any density feature that violates Rayleigh stability (or its generalization) should be diffused radially by hydrodynamical turbulence on a dynamical time-scale, approaching marginal stability in a quasi-continuous manner.

0 citations0 reviewsNo code linkedOpen access