Turbulence in particle laden midplane layers of planet forming disks

We examine the settled particle layers of planet forming disks in which the streaming instability (SI) is thought to be either weak or inactive. A suite of low-to-moderate resolution three-dimensional simulations in a $0.2H$ sized box, where $H$ is the pressure scale height, are performed using PENCIL for two Stokes numbers, \St$=0.04$ and $0.2$, at 1\% disk metallicity. We find a complex of Ekman-layer jet-flows emerge subject to three co-acting linearly growing processes: (1) the Kelvin-Helmholtz instability (KHI), (2) the planet-forming disk analog of the baroclinic Symmetric Instability (SymI), and (3) a later-time weakly acting secondary transition process, possibly a manifestation of the SI, producing a radially propagating pattern state. For \St$=0.2$, KHI is dominant and manifests as off-midplane axisymmetric rolls, while for \St$=0.04$ the axisymmetric SymI mainly drives turbulence. SymI is analytically developed in a model disk flow, predicting that it becomes strongly active when the Richardson number (Ri) of the particle-gas midplane layer transitions below 1, exhibiting growth rates \le\sqrt{2/\Ri - 2}\cdot\Omega, where $\Omega$ is local disk rotation rate. For fairly general situations absent external sources of turbulence it is conjectured that the SI, when and if initiated, emerges out of a turbulent state primarily driven and shaped by at least SymI and/or KHI. We also find that turbulence produced in $256^3$ resolution simulations are not statistically converged and that corresponding $512^3$ simulations may be converged for \St$=0.2$. Furthermore, we report that our numerical simulations significantly dissipate turbulent kinetic energy on scales less than 6-8 grid points.Comment: 55 pages, 27 figures, accepted for publication in Ap

Length and Velocity Scales in Protoplanetary Disk Turbulence

In the theory of protoplanetary disk turbulence, a widely adopted \emph{ansatz}, or assumption, is that the turnover frequency of the largest turbulent eddy, $\Omega_L$, is the local Keplerian frequency $\Omega_K$. In terms of the standard dimensionless Shakura-Sunyaev $\alpha$ parameter that quantifies turbulent viscosity or diffusivity, this assumption leads to characteristic length and velocity scales given respectively by $\sqrt{\alpha}H$ and $\sqrt{\alpha}c$, in which $H$ and $c$ are the local gas scale height and sound speed. However, this assumption is not applicable in cases when turbulence is forced numerically or driven by some natural processes such as Vertical Shear Instability. Here we explore the more general case where $\Omega_L\ge\Omega_K$ and show that under these conditions, the characteristic length and velocity scales are respectively $\sqrt{\alpha/R'}H$ and $\sqrt{\alpha R'}c$, where $R'\equiv \Omega_L/\Omega_K$ is twice the Rossby number. It follows that \alpha=\alphat/R', where \sqrt{\alphat} c is the root-mean-square average of the turbulent velocities. Properly allowing for this effect naturally explains the reduced particle scale heights produced in shearing box simulations of particles in forced turbulence, and may help with interpreting recent edge-on disk observations; more general implications for observations are also presented. For $R'>1$ the effective particle Stokes numbers are increased, which has implications for particle collision dynamics and growth, as well as for planetesimal formation.Comment: Accepted for publication in Ap

Back to the Roots: Predicting the Source Domain of Metaphors using Contrastive Learning

Metaphors frame a given target domain using concepts from another, usually more concrete, source domain. Previous research in NLP has focused on the identification of metaphors and the interpretation of their meaning. In contrast, this paper studies to what extent the source domain can be predicted computationally from a metaphorical text. Given a dataset with metaphorical texts from a finite set of source domains, we propose a contrastive learning approach that ranks source domains by their likelihood of being referred to in a metaphorical text. In experiments, it achieves reasonable performance even for rare source domains, clearly outperforming a classification baseline

Magnetically Driven Turbulence in the Inner Regions of Protoplanetary Disks

Given the important role turbulence plays in the settling and growth of dust grains in protoplanetary disks, it is crucial that we determine whether these disks are turbulent and to what extent. Protoplanetary disks are weakly ionized near the mid-plane, which has led to a paradigm in which largely laminar magnetic field structures prevail deeper in the disk, with angular momentum being transported via magnetically launched winds. Yet, there has been little exploration on the precise behavior of the gas within the bulk of the disk. We carry out 3D, local shearing box simulations that include all three low-ionization effects (Ohmic diffusion, ambipolar diffusion, and the Hall effect) to probe the nature of magnetically driven gas dynamics 1-30 AU from the central star. We find that gas turbulence can persist with a generous yet physically motivated ionization prescription (order unity Elsasser numbers). The gas velocity fluctuations range from 0.03-0.09 of the sound speed $c_s$ at the disk mid-plane to $\sim c_s$ near the disk surface, and are dependent on the initial magnetic field strength. However, the turbulent velocities do not appear to be strongly dependent on the field polarity, and thus appear to be insensitive to the Hall effect. The mid-plane turbulence has the potential to drive dust grains to collision velocities exceeding their fragmentation limit, and likely reduces the efficacy of particle clumping in the mid-plane, though it remains to be seen if this level of turbulence persists in disks with lower ionization levels.Comment: 23 pages, 19 figures, 2 tables, submitted to Ap

Planet formation: The case for large efforts on the computational side

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.Comment: White paper submitted to the Astro2020 decadal surve

Planet formation: The case for large efforts on the computational side

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade

Analysis of Enhanced OSPF for Routing Lightpaths in Optical Mesh Networks

We discuss enhancements to the OSPF protocol for routing and topology discovery in optical mesh networks. OSPF's opaque LSA mechanism is used to extend OSPF to disseminate optical resource related information through optical LSAs. Standard link-state database flooding mechanisms are used for distribution of optical LSAs. Each optical LSA carries optical resource information pertaining to a single optical link bundle between two adjacent OXCs, allowing for fine granularity changes in topology to be incorporated in path computation algorithms. OSPF packets are carried over a single IP control channel between adjacent OXCs. We analyze the performance of OSPF with optical extensions. Specifically, we compute control channel bandwidth used due to LSA updates

Knowing the Robot’s World

The term robotics is used in broader aspects now-a-days. One of its main implementation is in the game development. Robotics & artificial intelligence are used in this field almost interchangeably. In this demonstration project of Patnaik’s algorithm in the field of artificial intelligence we are trying to demonstrate how the intelligent robots traces path during a game. This is implemented using the most fundamental graphics tool known i.e. BGI graphics in C. Our demonstration project will also help in getting the idea of a robot’s path tracing algo