Source author record

Alexei G. Kritsuk

Alexei G. Kritsuk appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

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

18 published item(s)

preprint2020arXiv

Conformal invariance of weakly compressible two-dimensional turbulence

We study conformal invariance of vorticity clusters in weakly compressible two-dimensional turbulence at low Mach numbers. On the basis of very high resolution direct numerical simulation we demonstrate the scaling invariance of the inverse cascade with scaling close to Kolmogorov prediction. In this range of scales, the statistics of zero-vorticity isolines are found to be compatible with those of critical percolation, thus generalizing the results obtained in incompressible Navier-Stokes turbulence.

preprint2020arXiv

Physical Processes in Star Formation

Star formation is a complex multi-scale phenomenon that is of significant importance for astrophysics in general. Stars and star formation are key pillars in observational astronomy from local star forming regions in the Milky Way up to high-redshift galaxies. From a theoretical perspective, star formation and feedback processes (radiation, winds, and supernovae) play a pivotal role in advancing our understanding of the physical processes at work, both individually and of their interactions. In this review we will give an overview of the main processes that are important for the understanding of star formation. We start with an observationally motivated view on star formation from a global perspective and outline the general paradigm of the life-cycle of molecular clouds, in which star formation is the key process to close the cycle. After that we focus on the thermal and chemical aspects in star forming regions, discuss turbulence and magnetic fields as well as gravitational forces. Finally, we review the most important stellar feedback mechanisms.

preprint2013arXiv

A supersonic turbulence origin of Larson's laws

We revisit the origin of Larson's scaling laws describing the structure and kinematics of molecular clouds. Our analysis is based on recent observational measurements and data from a suite of six simulations of the interstellar medium, including effects of self-gravity, turbulence, magnetic field, and multiphase thermodynamics. Simulations of isothermal supersonic turbulence reproduce observed slopes in linewidth-size and mass-size relations. Whether or not self-gravity is included, the linewidth-size relation remains the same. The mass-size relation, instead, substantially flattens below the sonic scale, as prestellar cores start to form. Our multiphase models with magnetic field and domain size 200 pc reproduce both scaling and normalization of the first Larson law. The simulations support a turbulent interpretation of Larson's relations. This interpretation implies that: (i) the slopes of linewidth-size and mass-size correlations are determined by the inertial cascade; (ii) none of the three Larson laws is fundamental; (iii) instead, if one is known, the other two follow from scale invariance of the kinetic energy transfer rate. It does not imply that gravity is dynamically unimportant. The self-similarity of structure established by the turbulence breaks in star-forming clouds due to the development of gravitational instability in the vicinity of the sonic scale. The instability leads to the formation of prestellar cores with the characteristic mass set by the sonic scale. The high-end slope of the core mass function predicted by the scaling relations is consistent with the Salpeter power-law index.

preprint2013arXiv

Energy cascade and scaling in supersonic isothermal turbulence

Supersonic turbulence plays an important role in a number of extreme astrophysical and terrestrial environments, yet its understanding remains rudimentary. We use data from a three-dimensional simulation of supersonic isothermal turbulence to reconstruct an exact fourth-order relation derived analytically from the Navier-Stokes equations (Galtier and Banerjee, Phys. Rev. Lett., vol. 107, 2011, p. 134501). Our analysis supports a Kolmogorov-like inertial energy cascade in supersonic turbulence previously discussed on a phenomenological level. We show that two compressible analogues of the four-fifths law exist describing fifth- and fourth-order correlations, but only the fourth-order relation remains `universal' in a wide range of Mach numbers from incompressible to highly compressible regimes. A new approximate relation valid in the strongly supersonic regime is derived and verified. We also briefly discuss the origin of bottleneck bumps in simulations of compressible turbulence.

preprint2013arXiv

Enzo: An Adaptive Mesh Refinement Code for Astrophysics

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

preprint2012arXiv

Flux Correlations in Supersonic Isothermal Turbulence

Using data from a large-scale three-dimensional simulation of supersonic isothermal turbulence, we have tested the validity of an exact flux relation derived analytically from the Navier--Stokes equation by Falkovich, Fouxon and Oz [2010 New relations for correlation functions in Navier--Stokes turbulence. J. Fluid Mech. 644, 465]. That relation, for compressible barotropic fluids, was derived assuming turbulence generated by a large-scale force. However, compressible turbulence in simulations is usually initialized and maintained by a large-scale acceleration, as in gravity-driven astrophysical flows. We present a new approximate flux relation for isothermal turbulence driven by a large-scale acceleration, and find it in reasonable agreement with the simulation results.

preprint2012arXiv

The Two States of Star Forming Clouds

We examine the effects of self-gravity and magnetic fields on supersonic turbulence in isothermal molecular clouds with high resolution simulations and adaptive mesh refinement. These simulations use large root grids (512^3) to capture turbulence and four levels of refinement to capture high density, for an effective resolution of 8,196^3. Three Mach 9 simulations are performed, two super-Alfvénic and one trans-Alfvénic. We find that gravity splits the clouds into two populations, one low density turbulent state and one high density collapsing state. The low density state exhibits properties similar to non-self-gravitating in this regime, and we examine the effects of varied magnetic field strength on statistical properties: the density probability distribution function is approximately lognormal; velocity power spectral slopes decrease with field strength; alignment between velocity and magnetic field increases with field; the magnetic field probability distribution can be fit to a stretched exponential. The high density state is characterized by self-similar spheres; the density PDF is a power-law; collapse rate decreases with increasing mean field; density power spectra have positive slopes, P(ρ,k) \propto k; thermal-to-magnetic pressure ratios are unity for all simulations; dynamic-to-magnetic pressure ratios are larger than unity for all simulations; magnetic field distribution is a power-law. The high Alfvén Mach numbers in collapsing regions explain recent observations of magnetic influence decreasing with density. We also find that the high density state is found in filaments formed by converging flows, consistent with recent Herschel observations. Possible modifications to existing star formation theories are explored.

preprint2012arXiv

What Shapes the Structure of Molecular Clouds: Turbulence or Gravity?

We revisit the origin of Larson's scaling relations, which describe the structure and kinematics of molecular clouds, based on recent observations and large-scale simulations of supersonic turbulence. Using dimensional analysis, we first show that both linewidth-size and mass-size correlations observed on scales 0.1-50 pc can be explained by a simple conceptual theory of compressible turbulence without resorting to the often assumed virial equilibrium or detailed energy balance condition. The scaling laws can be consistently interpreted as a signature of supersonic turbulence with no need to invoke gravity. We then show how self-similarity of structure established by the turbulence breaks in star-forming clouds through development of gravitational instabilities in the vicinity of the sonic scale, l_s~0.1 pc, leading to the formation of prestellar cores.

preprint2011arXiv

Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

preprint2011arXiv

On the Density Distribution in Star-forming Interstellar Clouds

We use deep adaptive mesh refinement simulations of isothermal self-gravitating supersonic turbulence to study the imprints of gravity on the mass density distribution in molecular clouds. The simulations show that the density distribution in self-gravitating clouds develops an extended power-law tail at high densities on top of the usual lognormal. We associate the origin of the tail with self-similar collapse solutions and predict the power index values in the range from -7/4 to -3/2 that agree with both simulations and observations of star-forming molecular clouds.

preprint2011arXiv

Turbulent Clustering of Protoplanetary Dust and Planetesimal Formation

We study clustering of inertial particles in turbulent flows and discuss its applications to dust particles in protoplanetary disks. Using numerical simulations, we compute the radial distribution function (RDF), which measures the probability of finding particle pairs at given distances, and the probability density function of the particle concentration. The clustering statistics depend on the Stokes number, $St$, defined as the ratio of the particle friction timescale, $τ_{\rm p} $, to the Kolmogorov timescale in the flow. In the dissipation range, the clustering intensity strongly peaks at $St \simeq 1$, and the RDF for $St \sim 1$ shows a fast power-law increase toward small scales, suggesting that turbulent clustering may considerably enhance the particle collision rate. Clustering at inertial-range scales is of particular interest to the problem of planetesimal formation. At these scales, the strongest clustering is from particles with $τ_{\rm p}$ in the inertial range. Clustering of these particles occurs primarily around a scale where the eddy turnover time is $\simτ_{\rm p}$. Particles of different sizes tend to cluster at different locations, leading to flat RDFs between different particles at small scales. In the presence of multiple particle sizes, the overall clustering strength decreases as the particle size distribution broadens. We discuss particle clustering in recent models for planetesimal formation. We point out that, in the model based on turbulent clustering of chondrule-size particles, the probability of finding strong clusters that can seed planetesimals may have been significantly overestimated. We discuss various clustering mechanisms in simulations of planetesimal formation by gravitational collapse of dense clumps of meter-size particles, in particular the contribution from turbulent clustering due to the limited numerical resolution.

preprint2010arXiv

Interstellar Turbulence and Star Formation

We provide a brief overview of recent advances and outstanding issues in simulations of interstellar turbulence, including isothermal models for interior structure of molecular clouds and larger-scale multiphase models designed to simulate the formation of molecular clouds. We show how self-organization in highly compressible magnetized turbulence in the multiphase ISM can be exploited in simple numerical models to generate realistic initial conditions for star formation.

preprint2009arXiv

Piecewise Parabolic Method on a Local Stencil for Magnetized Supersonic Turbulence Simulation

Stable, accurate, divergence-free simulation of magnetized supersonic turbulence is a severe test of numerical MHD schemes and has been surprisingly difficult to achieve due to the range of flow conditions present. Here we present a new, higher order-accurate, low dissipation numerical method which requires no additional dissipation or local "fixes" for stable execution. We describe PPML, a local stencil variant of the popular PPM algorithm for solving the equations of compressible ideal magnetohydrodynamics. The principal difference between PPML and PPM is that cell interface states are evolved rather that reconstructed at every timestep, resulting in a compact stencil. Interface states are evolved using Riemann invariants containing all transverse derivative information. The conservation laws are updated in an unsplit fashion, making the scheme fully multidimensional. Divergence-free evolution of the magnetic field is maintained using the higher order-accurate constrained transport technique of Gardiner and Stone. The accuracy and stability of the scheme is documented against a bank of standard test problems drawn from the literature. The method is applied to numerical simulation of supersonic MHD turbulence, which is important for many problems in astrophysics, including star formation in dark molecular clouds. PPML accurately reproduces in three-dimensions a transition to turbulence in highly compressible isothermal gas in a molecular cloud model. The low dissipation and wide spectral bandwidth of this method make it an ideal candidate for direct turbulence simulations.

preprint2009arXiv

Self-organization in Turbulent Molecular Clouds: Compressional versus Solenoidal Modes

We use three-dimensional numerical simulations to study self-organization in supersonic turbulence in molecular clouds. Our numerical experiments describe decaying and driven turbulent flows with an isothermal equation of state, sonic Mach numbers from 2 to 10, and various degrees of magnetization. We focus on properties of the velocity field and, specifically, on the level of its potential (dilatational) component as a function of turbulent Mach number, magnetic field strength, and scale. We show how extreme choices of either purely solenoidal or purely potential forcing can reduce the extent of the inertial range in the context of periodic box models for molecular cloud turbulence. We suggest an optimized forcing to maximize the effective Reynolds number in numerical models.

preprint2009arXiv

The Power Spectrum of Turbulence in NGC 1333: Outflows or Large-Scale Driving?

Is the turbulence in cluster-forming regions internally driven by stellar outflows or the consequence of a large-scale turbulent cascade? We address this question by studying the turbulent energy spectrum in NGC 1333. Using synthetic 13CO maps computed with a snapshot of a supersonic turbulence simulation, we show that the VCS method of Lazarian and Pogosyan provides an accurate estimate of the turbulent energy spectrum. We then apply this method to the 13CO map of NGC 1333 from the COMPLETE database. We find the turbulent energy spectrum is a power law, E(k) k^-beta, in the range of scales 0.06 pc < ell < 1.5 pc, with slope beta=1.85\pm 0.04. The estimated energy injection scale of stellar outflows in NGC 1333 is ell_inj 0.3 pc, well resolved by the observations. There is no evidence of the flattening of the energy spectrum above the scale ell_inj predicted by outflow-driven simulations and analytical models. The power spectrum of integrated intensity is also a nearly perfect power law in the range of scales 0.16 pc < ell < 7.9 pc, with no feature above ell_inj. We conclude that the observed turbulence in NGC 1333 does not appear to be driven primarily by stellar outflows.

preprint2007arXiv

Two regimes of Turbulent Fragmentation and the stellar IMF from Primordial to Present Day Star Formation

The Padoan and Nordlund model of the stellar initial mass function (IMF) is derived from low order statistics of supersonic turbulence, neglecting gravity (e.g. gravitational fragmentation, accretion and merging). In this work the predictions of that model are tested using the largest numerical experiments of supersonic hydrodynamic (HD) and magneto-hydrodynamic (MHD) turbulence to date (~1000^3 computational zones) and three different codes (Enzo, Zeus and the Stagger Code). The model predicts a power law distribution for large masses, related to the turbulence energy power spectrum slope, and the shock jump conditions. This power law mass distribution is confirmed by the numerical experiments. The model also predicts a sharp difference between the HD and MHD regimes, which is recovered in the experiments as well, implying that the magnetic field, even below energy equipartition on the large scale, is a crucial component of the process of turbulent fragmentation. These results suggest that the stellar IMF of primordial stars may differ from that in later epochs of star formation, due to differences in both gas temperature and magnetic field strength. In particular, we find that the IMF of primordial stars born in turbulent clouds may be narrowly peaked around a mass of order 10 solar masses, as long as the column density of such clouds is not much in excess of 10^22 cm^-2.

preprint2005arXiv

Adaptive Mesh Refinement for Supersonic Molecular Cloud Turbulence

We performed a series of three-dimensional numerical simulations of supersonic homogeneous Euler turbulence with adaptive mesh refinement (AMR) and effective grid resolution up to 1024^3 zones. Our experiments describe non-magnetized driven supersonic turbulent flows with an isothermal equation of state. Mesh refinement on shocks and shear is implemented to cover dynamically important structures with the highest resolution subgrids and calibrated to match the turbulence statistics obtained from the equivalent uniform grid simulations. We found that at a level of resolution slightly below 512^3, when a sufficient integral/dissipation scale separation is first achieved, the fraction of the box volume covered by the AMR subgrids first becomes smaller than unity. At the higher AMR levels subgrids start covering smaller and smaller fractions of the whole volume, which scale with the Reynolds number as Re^{-1/4}. We demonstrate the consistency of this scaling with a hypothesis that the most dynamically important structures in intermittent supersonic turbulence are strong shocks with a fractal dimension of two. We show that turbulence statistics derived from AMR simulations and simulations performed on uniform grids agree surprisingly well, even though only a fraction of the volume is covered by AMR subgrids. Based on these results, we discuss the signature of dissipative structures in the statistical properties of supersonic turbulence and their role in overall flow dynamics.

preprint1996arXiv

On the Fundamental Properties of Dynamically Hot Galaxies

A two-component isothermal equilibrium model is applied to reproduce basic structural properties of dynamically hot stellar systems immersed in their massive dark haloes. The origin of the fundamental plane relation for giant ellipticals is naturally explained as a consequence of dynamical equilibrium in the context of the model. The existence of two galactic families displaying different behaviour in the luminosity--surface-brightness diagram is shown to be a result of a smooth transition from dwarfs, dominated by dark matter near the centre, to giants dominated by the luminous stellar component. The comparison of empirical scaling relations with model predictions suggests that probably a unique dissipative process was operating during the violent stage of development of stellar systems in the dark haloes, and the depth of the potential well controlled the observed luminosity of the resulting galaxies. The interpretation also provides some restrictions on the properties of dark haloes implied by the fundamental scaling laws.