Researcher profile

Orkan M. Umurhan

Orkan M. Umurhan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
12works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 published item(s)

preprint2022arXiv

Global Modeling of Nebulae With Particle Growth, Drift, and Evaporation Fronts. II. The Influence of Porosity on Solids Evolution

Incremental particle growth in turbulent protoplanetary nebulae is limited by a combination of barriers that can slow or stall growth. Moreover, particles that grow massive enough to decouple from the gas are subject to inward radial drift which could lead to the depletion of most disk solids before planetesimals can form. Compact particle growth is probably not realistic. Rather, it is more likely that grains grow as fractal aggregates which may overcome this so-called radial drift barrier because they remain more coupled to the gas than compact particles of equal mass. We model fractal aggregate growth and compaction in a viscously evolving solar-like nebula for a range of turbulent intensities $α_{\rm{t}} = 10^{-5}-10^{-2}$. We do find that radial drift is less influential for porous aggregates over much of their growth phase; however, outside the water snowline fractal aggregates can grow to much larger masses with larger Stokes numbers more quickly than compact particles, leading to rapid inward radial drift. As a result, disk solids outside the snowline out to $\sim 10-20$ AU are depleted earlier than in compact growth models, but outside $\sim 20$ AU material is retained much longer because aggregate Stokes numbers there remain lower initially. Nevertheless, we conclude even fractal models will lose most disk solids without the intervention of some leap-frog planetesimal forming mechanism such as the Streaming Instability (SI), though conditions for the SI are generally never satisfied, except for a brief period %for a brief stage around $\sim 0.2$ Myr at the snowline for $α_{\rm{t}}=10^{-5}$.

preprint2022arXiv

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Ω$, where $Ω$ 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.

preprint2021arXiv

A minimal model for vertical shear instability in protoplanetary accretion disks

The Vertical Shear Instability is an axisymmetric effect suggested to drive turbulence in the magnetically inactive zones of protoplanetary accretion disks. Here we examine its physical mechanism in analytically tractable ``minimal models" in three settings that include a uniform density fluid, a stratified atmosphere, and a shearing-box section of a protoplanetary disk. Each of these analyses show that the vertical shear instability's essence is similar to the slantwise convective symmetric instability in the mid-latitude Earth atmosphere, in the presence of vertical shear of the baroclinic jet stream, as well as mixing in the top layers of the Gulf Stream. We show that in order to obtain instability the fluid parcels' slope should exceed the slope of the mean absolute momentum in the disk radial-vertical plane. We provide a detailed and mutually self-consistent physical explanation from three perspectives: in terms of angular momentum conservation, as a dynamical interplay between a fluid's radial and azimuthal vorticity components, and from an energy perspective involving a generalized Solberg-Høiland Rayleigh condition. Furthermore, we explain why anelastic dynamics yield oscillatory unstable modes and isolate the oscillation mechanism from the instability one.

preprint2016arXiv

Critical layers and protoplanetary disk turbulence

A linear analysis of the zombie vortex instability is performed in a stratified shearing sheet setting for three model barotropic shear flows. The linear analysis is done by utilizing a Green's function formulation to resolve the critical layers of the associated normal-mode problem. The instability is the result of a resonant interaction between a Rossby wave and a gravity wave which we refer to as Z-modes. The associated critical layer is the location where the Doppler shifted frequency of a distant Rossby wave equals the local Brunt-Vaisala frequency. The minimum required Rossby number for instability, Ro= 0.2, is confirmed for parameter values reported in the literature. It is also found that the shear layer supports the instability in the limit where stratification vanishes. The zombie vortex instability is examined in a jet model, finding that the instability can occur for Ro= 0.05. Nonlinear vorticity forcing due to unstable Z-modes is shown to result in the creation of a jet flow at the critical layer emerging as the result of the competition between the vertical lifting of perturbation radial vorticity and the radial transport of perturbation vertical vorticity. We find that the picture of this instability leading to a form of nonlinearly driven self-replicating pattern of creation and destruction is warranted: a parent jet spawns a growing child jet at associated critical layers. A mature child jet creates a next generation of child jets at associated critical layers while simultaneously destroying its parent jet via the Rossby wave instability.

preprint2016arXiv

On the mechanism of self gravitating Rossby interfacial waves in proto-stellar accretion discs

The dynamical response of edge waves under the influence of self-gravity is examined in an idealized two-dimensional model of a proto-stellar disc, characterized in steady state as a rotating vertically infinite cylinder of fluid with constant density except for a single density interface at some radius r0. The fluid in basic state is prescribed to rotate with a Keplerian profile $Ω_k(r)\sim r^{-3/2}$ modified by some additional azimuthal sheared flow. A linear analysis shows that there are two azimuthally propagating edge waves, kin to the familiar Rossby waves and surface gravity waves in terrestrial studies, which move opposite to one another with respect to the local basic state rotation rate at the interface. Instability only occurs if the radial pressure gradient is opposite to that of the density jump (unstably stratified) where self-gravity acts as a wave stabilizer irrespective of the stratification of the system. The propagation properties of the waves are discussed in detail in the language of vorticity edge waves. The roles of both Boussinesq and non-Boussinesq effects upon the stability and propagation of these waves with and without the inclusion of self-gravity are then quantified. The dynamics involved with self-gravity non- Boussinesq effect is shown to be a source of vorticity production where there is a jump in the basic state density, in addition, self-gravity also alters the dynamics via the radial main pressure gradient, which is a Boussinesq effect . Further applications of these mechanical insights are presented in the conclusion including the ways in which multiple density jumps or gaps may or may not be stable.

preprint2016arXiv

The Geology of Pluto and Charon Through the Eyes of New Horizons

NASA's New Horizons spacecraft has revealed the complex geology of Pluto and Charon. Pluto's encounter hemisphere shows ongoing surface geological activity centered on a vast basin containing a thick layer of volatile ices that appears to be involved in convection and advection, with a crater retention age no greater than $\approx$10 Ma. Surrounding terrains show active glacial flow, apparent transport and rotation of large buoyant water-ice crustal blocks, and pitting, likely by sublimation erosion and/or collapse. More enigmatic features include tall mounds with central depressions that are conceivably cryovolcanic, and ridges with complex bladed textures. Pluto also has ancient cratered terrains up to ~4 Ga old that are extensionally fractured and extensively mantled and perhaps eroded by glacial or other processes. Charon does not appear to be currently active, but experienced major extensional tectonism and resurfacing (probably cryovolcanic) nearly 4 billion years ago. Impact crater populations on Pluto and Charon are not consistent with the steepest proposed impactor size-frequency distributions proposed for the Kuiper belt.

preprint2016arXiv

Vortex formation in protoplanetary discs induced by the vertical shear instability

We present the results of 2D and 3D hydrodynamic simulations of idealized protoplanetary discs that examine the formation and evolution of vortices by the vertical shear instability (VSI). In agreement with recent work, we find that discs with radially decreasing temperature profiles and short thermal relaxation time-scales, are subject to the axisymmetric VSI. In three dimensions, the resulting velocity perturbations give rise to quasi-axisymmetric potential vorticity perturbations that break-up into discrete vortices, in a manner that is reminiscent of the Rossby wave instability. Discs with very short thermal evolution time-scales (i.e. τ<0.1 local orbit periods) develop strong vorticity perturbations that roll up into vortices that have small aspect ratios (χ<2) and short lifetimes (~ a few orbits). Longer thermal time-scales give rise to vortices with larger aspect ratios (6<χ<10), and lifetimes that depend on the entropy gradient. A steeply decreasing entropy profile leads to vortex lifetimes that exceed the simulation run times of hundreds of orbital periods. Vortex lifetimes in discs with positive or weakly decreasing entropy profiles are much shorter, being 10s of orbits at most, suggesting that the subcritical baroclinic instability plays an important role in sustaining vortices against destruction through the elliptical instability. Applied to the outer regions of protoplanetary discs, where the VSI is most likely to occur, our results suggest that vortices formed by the VSI are likely to be short lived structures.

preprint2014arXiv

Interacting vorticity waves as an instability mechanism for MHD shear instabilities

The interacting vorticity wave formalism for shear flow instabilities is extended here to the magnetohydrodynamic (MHD) setting, to provide a mechanistic description for the stabilising and destabilising of shear instabilities by the presence of a background magnetic field. The interpretation relies on local vorticity anomalies inducing a non-local velocity field, resulting in action-at-a-distance. It is shown here that the waves supported by the system are able to propagate vorticity via the Lorentz force, and waves may interact; existence of instability then rests upon whether the choice of basic state allows for phase-locking and constructive interference of the vorticity waves via mutual interaction. To substantiate this claim, we solve the instability problem of two representative basic states, one where a background magnetic field stabilises an unstable flow and the other where the field destabilises a stable flow, and perform relevant analyses to show how this mechanism operates in MHD.

preprint2013arXiv

Intercomparison of General Circulation Models for Hot Extrasolar Planets

We compare five general circulation models (GCMs) which have been recently used to study hot extrasolar planet atmospheres (BOB, CAM, IGCM, MITgcm, and PEQMOD), under three test cases useful for assessing model convergence and accuracy. Such a broad, detailed intercomparison has not been performed thus far for extrasolar planets study. The models considered all solve the traditional primitive equations, but employ different numerical algorithms or grids (e.g., pseudospectral and finite volume, with the latter separately in longitude-latitude and `cubed-sphere&#39; grids). The test cases are chosen to cleanly address specific aspects of the behaviors typically reported in hot extrasolar planet simulations: 1) steady-state, 2) nonlinearly evolving baroclinic wave, and 3) response to fast timescale thermal relaxation. When initialized with a steady jet, all models maintain the steadiness, as they should -- except MITgcm in cubed-sphere grid. A very good agreement is obtained for a baroclinic wave evolving from an initial instability in pseudospectral models (only). However, exact numerical convergence is still not achieved across the pseudospectral models: amplitudes and phases are observably different. When subject to a typical `hot-Jupiter&#39;-like forcing, all five models show quantitatively different behavior -- although qualitatively similar, time-variable, quadrupole-dominated flows are produced. Hence, as have been advocated in several past studies, specific quantitative predictions (such as the location of large vortices and hot regions) by GCMs should be viewed with caution. Overall, in the tests considered here, pseudospectral models in pressure coordinate (PEBOB and PEQMOD) perform the best and MITgcm in cubed-sphere grid performs the worst.

preprint2012arXiv

Linear and nonlinear evolution of the vertical shear instability in accretion discs

(Abridged) We analyse the stability and evolution of power-law accretion disc models. These have midplane densities that follow radial power-laws, and have either temperature or entropy distributions that are power-law functions of cylindrical radius. We employ two different hydrodynamic codes to perform 2D-axisymmetric and 3D simulations that examine the long-term evolution of the disc models as a function of the power-law indices of the temperature or entropy, the thermal relaxation time of the fluid, and the viscosity. We present a stability analysis of the problem that we use to interpret the simulation results. We find that disc models whose temperature or entropy profiles cause the equilibrium angular velocity to vary with height are unstable to the growth of modes with wavenumber ratios |k_R/k_Z| >> 1 when the thermodynamic response of the fluid is isothermal, or the thermal evolution time is comparable to or shorter than the local dynamical time scale. These discs are subject to the Goldreich-Schubert-Fricke (GSF) or `vertical shear&#39; linear instability. Development of the instability involves excitation of vertical breathing and corrugation modes in the disc, with the corrugation modes in particular being a feature of the nonlinear saturated state. Instability operates when the dimensionless disc kinematic viscosity nu < 10^{-6} (Reynolds numbers Re>H c_s/nu > 2500). In 3D the instability generates a quasi-turbulent flow, and the Reynolds stress produces a fluctuating effective viscosity coefficient whose mean value reaches alpha ~ 6 x 10^{-4} by the end of the simulation. The vertical shear instability in disc models which include realistic thermal physics has yet to be examined. Should it occur, however, our results suggest that it will have significant consequences for their internal dynamics, transport properties, and observational appearance.

preprint2012arXiv

Non-Dissipative Saturation of the Magnetorotational Instability in Thin Disks

A new non-dissipative mechanism is proposed for the saturation of the axisymmetric magnetorotational (MRI) instability in thin Keplerian disks that are subject to an axial magnetic field. That mechanism relies on the energy transfer from the MRI to stable magnetosonic (MS) waves. Such mode interaction is enabled due to the vertical stratification of the disk that results in the discretization of its MRI spectrum, as well as by applying the appropriate boundary conditions. A second order Duffing-like amplitude equation for the initially unstable MRI modes is derived. The solutions of that equation exhibit bursty nonlinear oscillations with a constant amplitude that signifies the saturation level of the MRI. Those results are verified by a direct numerical solution of the full nonlinear reduced set of thin disk magnetohydrodynamics equations.

preprint2009arXiv

Global transient dynamics of three-dimensional hydrodynamical disturbances in a thin viscous accretion disk

Thin viscous Keplerian accretion disks are considered asymptotically stable, even though they can show significant dynamic activity on short timescales. In this paper the dynamics of non-axisymmetric hydrodynamical disturbances of disks are investigated analytically building upon the steady state three-dimensional structure and evolution of axisymmetric perturbations explored in previous work. Assuming a polytropic equation of state solutions are found by means of an asymptotic expansion in the small parameter measuring the ratio of the disk thickness to characteristic radius. In-depth analysis shows that every perturbation that disturbs the radial velocity induces significant transient growth in the (acoustic) energy of the evolving disturbance. This effect is most evident in the density and vertical velocity. The transient growth observed is tied to the non-separable nature of the solutions where, in particular, pattern evolution is controlled by a similarity variable composed of the radial coordinate and time. This leads to growing winding perturbations that display successive radial peaks and troughs. We argue that these transient non-axisymmetric structures may precipitate secondary instabilities which, consequently, may be a critical element for a new alternative picture of turbulence arousal in non-magnetized astrophysical disks.