The Janus Face of Turbulent Pressure

Size: px
Start display at page:

Download "The Janus Face of Turbulent Pressure"

Transcription

1 CRC 963 Astrophysical Turbulence and Flow Instabilities with thanks to Christoph Federrath, Monash University Patrick Hennebelle, CEA/Saclay Alexei Kritsuk, UCSD yt-project.org Seminar über Astrophysik, Hamburg, November 2012

2 Overview Chandrasekhar (1951), etc.: Compressible turbulence increases the effective pressure of the gas: c 2 eff = c2 s u2 But is it really that simple?

3 Overview Chandrasekhar (1951), etc.: Compressible turbulence increases the effective pressure of the gas: c 2 eff = c2 s u2 But is it really that simple? Different routes: 1. Subgrid-scale turbulent pressure 2. Virial turbulent pressure 3.

4 Scale Separation in Turbulence Simulations

5 Decomposition of the Equations of Fluid Dynamics Compressible Navier-Stokes equation for momentum: ρu + (ρu u) = P + σ + F t Smoothing over length scale by a spatial filter yields: ρu + ρu u = P + σ +F t }{{} 0 if L/Re 3/4 Define ρ = ρ, ū = ρu / ρ, etc., then t ρū + ρu u = P + F Looks almost like the compressible Euler equation for ρū...

6 Turbulent Stresses... if we substitute the turbulence stress tensor τ = ρu u + ρū ū for the non-linear interactions across so that t ρū + ρū ū = P + τ + F Trace of τ defines the turbulent pressure: τ ii ρū 2 ρu 2 = 2ρK = 3P sgs ρk is the energy density of unresolved (l ) turbulent velocity fluctuations Effective pressure term in the Euler equation is given by P ( ( ) P + P sgs = P + 2 ) 3 ρk

7 How to Compute the Turbulent Pressure From the decomposition follows an additional PDE for ρk (e. g., Germano 1992, Schmidt et al. 2006): [ ] ρk + ( ρūk) = C t }{{} κ ρ K 1/2 K }{{} advection SGS transport (diffusion) + τ ij S ij 2 3 ρk d }{{} turbulent cascade C ɛ ρ K 3/2 }{{} dissipation Production rate Σ = τ ij Sij, where (Woodward et al. 2006, Schmidt & Federrath 2011) τ ij = C 1 ρ K 1/2 S ij }{{} linear eddy-viscosity part 2C 2 ρk 2ū i,kū j,k ū 2 } {{ } non-linear part Supersonic turbulence: C and C (1 C 2) ρkδ ij

8 The subgrid-scale turbulent pressure LES of Supersonic Turbulence Denstrophy 12 (ρ1/2 v) 2 in a 5123 LES (Enzo), Mrms 5.5 (Schmidt & Federrath 2011) solenoidal forcing compressive forcing

9 The subgrid-scale turbulent pressure LES of Supersonic Turbulence SGS turbulence energy ρk in a 5123 LES (Enzo), Mrms 5.5 (Schmidt & Federrath 2011) solenoidal forcing compressive forcing

10 SGS Turbulence Energy Statistics K sgs Ρ 0 V LES with grid resolutions changing from 64 3 to : Mean SGS turbulence energy follows power-laws with slope β 0.8 (steeper than Kolmogorov scaling) K sgs Ρ 0 V 2 K sgs Ρ 0 V

11 Effective Pressure Microscopic equation of state for isothermal gas: P ρ Effective pressure on the grid scale: P eff = P+ 2 3 ρk ρc2 0 ( ) 3 M2 sgs LES, solenoidal forcing, M rms Locally, the turbulent pressure can largely exceed the thermal pressure 1 cdf sgs log 10 P eff P log 10 Ρ Ρ 0

12 The General Virial Theorem The Euler equation including gravitational forces can be contracted with r and integrated over the volume V By defining I = ρr 2 d 3 x, E int = V E int = V φ E grav = ρx i d 3 x, E kin = V E kin = x i V it follows that 1 2Ï = 2(E int + E kin ) + E grav + V V V 3 2 P d3 x 1 2 ρu2 d 3 x ( Pδij + u 2) x i n j ds

13 Virial Equilibrium Under the assumptions equlibrium (Ï = 0) surfaces terms negligible no external gravitational potential we obtain the well-known relation 2(E int + E kin ) + E grav 0 The first two terms can be expressed as where 2(E int + E kin ) = 3 (P + P t ) P = P d 3 x and P t = 1 V 3 V ρu 2 d 3 x Applicable to molecular clouds, even to substructure in molecular clouds (clumps, cores)?

14 Numerical Clump Mass Distributions Schmidt, Kern, Federrath, and Klessen (2010): Isothermal turbulence simulations with assumed mass scale given by L = 0.6 pc and n = 10 4 cm 3 (λ 0 J /L 0.2) Simply connected regions exceeding their Bonnor-Ebert mass are identified by a clump finding algorithm (Padoan & Nordlund 2007) 2E int = 3P E grav: 3(P + P t) E grav:

15 Hennebelle-Chabrier Theory of the Clump Mass Function Hennebelle & Chabrier (2008): Based on the Press-Schechter statistical formalism for hierarchical objects, the mass distribution is related to the PDF P(δ) of δ = log(ρ/ ρ): N (M) ρ ( M dm dr ) 1 dδ dr exp(δ)p(δ) Clump scale R is connected to its mass M By assuming the Jeans- or Bonnor-Ebert-mass for gravitationally unstable cores, we have M = R = exp( δ/2), Can now calculate N (M) for given P(δ) if R L where M = M R MJ 0, R = λ 0 J

16 The Clump Mass Function with Turbulent Support Define effective Jeans mass M J,eff through the effective sound speed (in virial equilbr.) ceff 2 = c δu2 (R) [ c ( ) η ] R 3 M rms L New parameter: turbulent Mach number scaled down to Jeans length at the mean density 3 M = M rms ( λ 0 J /L ) η Can compute N (M) for given P(δ), M, and η for HD turbulence with soln./compr. forcing (Federrath et al. 2010, Schmidt et al. 2010) With these definitions: ( ( M = R 1 + M 2 R 2η) 1 + M 2, δ = log R 2 R 2η )

17 Turbulent Clouds Isolated system: virial equilibrium, turbulent pressure calculated by Bonazolla (1992) Substructure: boundary terms, turbulent pressure? white

18 Caveats Bonazzola et al. (1992) showed that effective sound speed is only valid on( length scales integral length scale CMF theory and numerical analysis are based on the assumption of local virial equilibrium for the critical mass Questioned by Ballesteros-Paredes (2006) Dib et al. (2007) demonstrate with numerical data that cores/clumps are non-equilibrium structures: significant contributions from boundary terms undergoing supersonic compression while not necessarily gravitationally bound Flat slope can also be explained by magnetohydrodynamic turbulence without turbulent pressure support

19 The Rate of Compression Equation for the rate of compression (d = u): Dd Dt = 2 φ Λ turb Λ therm = 4πG(ρ ρ 0 ) }{{} gravity 1 ( ω 2 S 2) 1 ( 2 P + 1ρ ) } 2 {{} ρ ρ P }{{} turbulence thermal pressure hydrostatic equilibrium for d = ω = S = 0 Jeans criterion follows in the limit δ = ρ/ρ free-fall collapse with τff (Gρ) 1/2 for 4πG(ρ ρ 0 )/ Λ 0 Schmidt, Collins, and Kritsuk (2012): Compute statistics of positive/negative support for self-gravitating turbulence simulations: Λ + = { Λ if Λ 0 0 otherwise Λ = { Λ if Λ 0 0 otherwise

20 Self-Gravitating Turbulence Kritsuk et al. (2011): α 0.25 and M rms 6 Deep AMR simulation with 5 levels of refinement, refinement factor 4

21 Thermal Support Mean thermal support is positive Asymptotic relation: 2 ρλ therm P for P 100 Intermittency: strong fluctuations produced by filaments/cores undergoing strong compressions Strong support around δ 10 7 corresponds to flattening of the density PDF

22 Turbulent Support However, mean turbulent support is negative Asymptotic relation: 2 ρλ turb ρω 2 for ρω 2 1 Negative pressure caused by shear and supersonic compression overcompensates positive pressure of turbulent eddies

23 Power Spectra 500 P Ω k P turb k P S k k 2 Fourier-transform of root-grid data (smoothing factor up to 1024) Burgers scaling ω(l) l 1/2 implies a flat power spectrum: ( ) 1 P turb (k) = 2 ˆω i ˆω i ŜijŜij k 2 dω k = 1 2 [P ω(k) P S (k)] k =k

24 Total Support For 100 δ 10 4, net support has negative floor with Λ 4πGρ 0 δ: supersonic gas compression triggers gravitational collapse this corresponds to the tail of the density PDF with slope to the 1.7 (compared to 1.5 for pressure-free collapse solution by Penston 1969)

25 Resume In the LES picture, something like Chandrasekhar s effective pressure follows rigorously from scale separation Applies approximately to substructure (clumps/cores) in turbulent star-forming clouds if local virial equilibrium can be assumed and boundary terms can be neglected Turbulent support against gravity appears to be negative in the non-linear non-equilibrium regime, so the nature of turbulent pressure is profoundly different from the thermal and magnetic pressures Are these conclusions somehow in contradiction? No, because Turbulent pressure is just the diagonal part of the full turbulence stress tensor Support is locally given by the pressure Laplacian, not the pressure itself

Subgrid Scale Physics in Galaxy Simulations

Subgrid Scale Physics in Galaxy Simulations CRC 963 Astrophysical Turbulence and Flow Instabilities with thanks to Harald Braun, Jan Frederic Engels, Jens Niemeyer, IAG Ann Almgren and John Bell, LBNL Christoph Federrath, Monash University yt-project.org

More information

Local support against gravity in magnetoturbulent fluids

Local support against gravity in magnetoturbulent fluids MNRAS 431, 3196 3215 (2013) Advance Access publication 2013 April 3 doi:10.1093/mnras/stt399 Local support against gravity in magnetoturbulent fluids W. Schmidt, 1 D. C. Collins 2,3 anda.g.kritsuk 3 1

More information

Alexei Kritsuk. CalSpace IGPP/UCR Conference Palm Springs, March 27-30, 2006 c Alexei Kritsuk

Alexei Kritsuk.   CalSpace IGPP/UCR Conference Palm Springs, March 27-30, 2006 c Alexei Kritsuk Supersonic Turbulence in Molecular Clouds 1 Alexei Kritsuk University of California, San Diego In collaboration with Rick Wagner, Mike Norman, & Paolo Padoan, UCSD http://akpc.ucsd.edu/isothermal CalSpace

More information

Turbulence simulations with ENZO and FLASH3

Turbulence simulations with ENZO and FLASH3 Turbulence simulations with ENZO and FLASH3... not yet Adaptive-mesh simulations with FLASH Christoph Federrath Institute for Theoretical Astrophysics Heidelberg Collaborators: Ralf Klessen, Robi Banerjee,

More information

Mixing in Highly Compressible Turbulence

Mixing in Highly Compressible Turbulence Mixing in Highly Compressible Turbulence Liubin Pan Evan Scannapieco School of Earth and Space Exploration Arizona State University Astrophysical motivation: Heavy elements from supernova and stellar winds

More information

Computational Fluid Dynamics 2

Computational Fluid Dynamics 2 Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2

More information

Theory of star formation

Theory of star formation Theory of star formation Monday 8th 17.15 18.00 Molecular clouds and star formation: Introduction Tuesday 9th 13.15 14.00 Molecular clouds: structure, physics, and chemistry 16.00 16.45 Cloud cores: statistics

More information

A supersonic turbulence origin of Larson s laws

A supersonic turbulence origin of Larson s laws MNRAS 436, 3247 3261 (2013) Advance Access publication 2013 October 17 doi:10.1093/mnras/stt1805 A supersonic turbulence origin of Larson s laws Alexei G. Kritsuk, 1 Christoph T. Lee 1 and Michael L. Norman

More information

2. Conservation Equations for Turbulent Flows

2. Conservation Equations for Turbulent Flows 2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging

More information

The Initial Mass Function Elisa Chisari

The Initial Mass Function Elisa Chisari The Initial Mass Function AST 541 Dec 4 2012 Outline The form of the IMF Summary of observations Ingredients of a complete model A Press Schechter model Numerical simulations Conclusions The form of the

More information

Modelling of turbulent flows: RANS and LES

Modelling of turbulent flows: RANS and LES Modelling of turbulent flows: RANS and LES Turbulenzmodelle in der Strömungsmechanik: RANS und LES Markus Uhlmann Institut für Hydromechanik Karlsruher Institut für Technologie www.ifh.kit.edu SS 2012

More information

A supersonic turbulence origin of Larson s laws

A supersonic turbulence origin of Larson s laws Mon. Not. R. Astron. Soc., 1 17 (213) Printed 9 October 213 (MN LATEX style file v2.2) A supersonic turbulence origin of Larson s laws Alexei G. Kritsuk, 1 Christoph T. Lee 1 and Michael L. Norman 1,2

More information

TWO REGIMES OF TURBULENT FRAGMENTATION AND THE STELLAR INITIAL MASS FUNCTION FROM PRIMORDIAL TO PRESENT-DAY STAR FORMATION

TWO REGIMES OF TURBULENT FRAGMENTATION AND THE STELLAR INITIAL MASS FUNCTION FROM PRIMORDIAL TO PRESENT-DAY STAR FORMATION The Astrophysical Journal, 661:972Y981, 2007 June 1 # 2007. The American Astronomical Society. All rights reserved. Printed in U.S.A. A TWO REGIMES OF TURBULENT FRAGMENTATION AND THE STELLAR INITIAL MASS

More information

Lecture 22 Stability of Molecular Clouds

Lecture 22 Stability of Molecular Clouds Lecture 22 Stability of Molecular Clouds 1. Stability of Cloud Cores 2. Collapse and Fragmentation of Clouds 3. Applying the Virial Theorem References Myers, Physical Conditions in Molecular Clouds in

More information

dt 2 = 0, we find that: K = 1 2 Ω (2)

dt 2 = 0, we find that: K = 1 2 Ω (2) 1 1. irial Theorem Last semester, we derived the irial theorem from essentially considering a series of particles which attract each other through gravitation. The result was that d = K + Ω (1) dt where

More information

Simulations of Supersonic Turbulence in Molecular Clouds: Evidence for a New Universality

Simulations of Supersonic Turbulence in Molecular Clouds: Evidence for a New Universality NUMERICAL MODELING OF SPACE PLASMA FLOWS: ASTRONUM-28 ASP Conference Series, Vol. 46, 29 Nikolai V. Pogorelov, Edouard Audit, Phillip Colella, and Gary P. Zank, eds. Simulations of Supersonic Turbulence

More information

Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace

Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and

More information

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

Fundamentals of Fluid Dynamics: Elementary Viscous Flow Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research

More information

arxiv: v1 [astro-ph.ga] 16 Jul 2009

arxiv: v1 [astro-ph.ga] 16 Jul 2009 Analytical theory for the initial mass function: II. Properties of the flow. Patrick Hennebelle arxiv:0907.2765v1 [astro-ph.ga] 16 Jul 2009 Laboratoire de radioastronomie, UMR CNRS 8112, École normale

More information

Quick Recapitulation of Fluid Mechanics

Quick Recapitulation of Fluid Mechanics Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change

More information

The nature of the velocity field in molecular clouds I. The non-magnetic case

The nature of the velocity field in molecular clouds I. The non-magnetic case Mon. Not. R. Astron. Soc. 390, 769 780 (2008) doi:10.1111/j.1365-2966.2008.13778.x The nature of the velocity field in molecular clouds I. The non-magnetic case Enrique Vázquez-Semadeni, 1 Ricardo F. González,

More information

Enrique Vázquez-Semadeni. Centro de Radioastronomía y Astrofísica, UNAM, México

Enrique Vázquez-Semadeni. Centro de Radioastronomía y Astrofísica, UNAM, México Enrique Vázquez-Semadeni Centro de Radioastronomía y Astrofísica, UNAM, México 1 Javier Ballesteros-Paredes Centro de Radioastronomía y Astrofísica, UNAM, México 2 Collaborators: Javier Ballesteros-Paredes

More information

Incompressible MHD simulations

Incompressible MHD simulations Incompressible MHD simulations Felix Spanier 1 Lehrstuhl für Astronomie Universität Würzburg Simulation methods in astrophysics Felix Spanier (Uni Würzburg) Simulation methods in astrophysics 1 / 20 Outline

More information

Turbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i

Turbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)

More information

A localised subgrid scale model for fluid dynamical simulations in astrophysics. I. Theory and numerical tests ABSTRACT. N cr N 4/3 Re 3.

A localised subgrid scale model for fluid dynamical simulations in astrophysics. I. Theory and numerical tests ABSTRACT. N cr N 4/3 Re 3. A&A 450, 265 281 (2006) DOI: 10.1051/0004-6361:20053617 c ESO 2006 Astronomy & Astrophysics A localised subgrid scale model for fluid dynamical simulations in astrophysics I. Theory and numerical tests

More information

Regularization modeling of turbulent mixing; sweeping the scales

Regularization modeling of turbulent mixing; sweeping the scales Regularization modeling of turbulent mixing; sweeping the scales Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) D 2 HFest, July 22-28, 2007 Turbulence

More information

The Johns Hopkins Turbulence Databases (JHTDB)

The Johns Hopkins Turbulence Databases (JHTDB) The Johns Hopkins Turbulence Databases (JHTDB) HOMOGENEOUS BUOYANCY DRIVEN TURBULENCE DATA SET Data provenance: D. Livescu 1 Database Ingest and Web Services: C. Canada 1, K. Kalin 2, R. Burns 2 & IDIES

More information

Gravity or Turbulence?

Gravity or Turbulence? Gravity or Turbulence? On the dynamics of Molecular Clouds Javier Ballesteros-Paredes On Sabbatical at Institut für Theoretische Astrophysik, University of Heidelberg Instituto de Radioastronomía y Astrofísica,

More information

Lecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.

Lecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. Lecture 14 Turbulent Combustion 1 We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. In a fluid flow, turbulence is characterized by fluctuations of

More information

Mapping the core mass function to the initial mass function

Mapping the core mass function to the initial mass function doi:10.1093/mnras/stv872 Mapping the core mass function to the initial mass function Dávid Guszejnov and Philip F. Hopkins TAPIR, Mailcode 350-17, California Institute of Technology, Pasadena, CA 91125,

More information

formation of the cosmic large-scale structure

formation of the cosmic large-scale structure formation of the cosmic large-scale structure Heraeus summer school on cosmology, Heidelberg 2013 Centre for Astronomy Fakultät für Physik und Astronomie, Universität Heidelberg August 23, 2013 outline

More information

Collapse of magnetized dense cores. Is there a fragmentation crisis?

Collapse of magnetized dense cores. Is there a fragmentation crisis? Collapse of magnetized dense cores Is there a fragmentation crisis? Patrick Hennebelle (ENS-Observatoire de Paris) Collaborators: Benoît Commerçon, Andréa Ciardi, Sébastien Fromang, Romain Teyssier, Philippe

More information

Introduction to Turbulence and Turbulence Modeling

Introduction to Turbulence and Turbulence Modeling Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly

More information

Multiscale Computation of Isotropic Homogeneous Turbulent Flow

Multiscale Computation of Isotropic Homogeneous Turbulent Flow Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible

More information

Dissipation Scales & Small Scale Structure

Dissipation Scales & Small Scale Structure Dissipation Scales & Small Scale Structure Ellen Zweibel zweibel@astro.wisc.edu Departments of Astronomy & Physics University of Wisconsin, Madison and Center for Magnetic Self-Organization in Laboratory

More information

Enrique Vázquez-Semadeni. Centro de Radioastronomía y Astrofísica, UNAM, México

Enrique Vázquez-Semadeni. Centro de Radioastronomía y Astrofísica, UNAM, México Enrique Vázquez-Semadeni Centro de Radioastronomía y Astrofísica, UNAM, México 1 Collaborators: CRyA UNAM: Abroad: Javier Ballesteros-Paredes Pedro Colín Gilberto Gómez Recent PhDs: Alejandro González

More information

arxiv:astro-ph/ v1 26 Sep 2003

arxiv:astro-ph/ v1 26 Sep 2003 Star Formation at High Angular Resolution ASP Conference Series, Vol. S-221, 2003 M.G. Burton, R. Jayawardhana & T.L. Bourke The Turbulent Star Formation Model. Outline and Tests arxiv:astro-ph/0309717v1

More information

Fundamental Issues in Star Formation

Fundamental Issues in Star Formation Fundamental Issues in Star Formation - Formation and statistical properties of dense molecular cloud cores (mass function of cores, scaling relations, gravitational boundedness, rotational properties)

More information

3 Hydrostatic Equilibrium

3 Hydrostatic Equilibrium 3 Hydrostatic Equilibrium Reading: Shu, ch 5, ch 8 31 Timescales and Quasi-Hydrostatic Equilibrium Consider a gas obeying the Euler equations: Dρ Dt = ρ u, D u Dt = g 1 ρ P, Dɛ Dt = P ρ u + Γ Λ ρ Suppose

More information

Chapter 7. Basic Turbulence

Chapter 7. Basic Turbulence Chapter 7 Basic Turbulence The universe is a highly turbulent place, and we must understand turbulence if we want to understand a lot of what s going on. Interstellar turbulence causes the twinkling of

More information

Mass Transfer in Turbulent Flow

Mass Transfer in Turbulent Flow Mass Transfer in Turbulent Flow ChEn 6603 References: S.. Pope. Turbulent Flows. Cambridge University Press, New York, 2000. D. C. Wilcox. Turbulence Modeling for CFD. DCW Industries, La Caada CA, 2000.

More information

arxiv: v1 [astro-ph] 5 Jun 2008

arxiv: v1 [astro-ph] 5 Jun 2008 Draft version December 15, 2017 Preprint typeset using L A TEX style emulateapj v. 08/13/06 DRIVEN AND DECAYING TURBULENCE SIMULATIONS OF LOW-MASS STAR FORMATION: FROM CLUMPS TO CORES TO PROTOSTARS Stella

More information

Hydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010

Hydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 Hydrodynamics Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 What is Hydrodynamics? Describes the evolution of physical systems (classical or quantum particles, fluids or fields) close to thermal

More information

Numerical Study of Compressible Isothermal Magnetohydrodynamic Turbulence

Numerical Study of Compressible Isothermal Magnetohydrodynamic Turbulence Numerical Study of Compressible Isothermal Magnetohydrodynamic Turbulence Junseong Park, Dongsu Ryu Dept. of Physics, Ulsan National Institute of Science and Technology. Ulsan, Korea 2016 KNAG meeting

More information

Physics of turbulent flow

Physics of turbulent flow ECL-MOD 3A & MSc. Physics of turbulent flow Christophe Bailly Université de Lyon, Ecole Centrale de Lyon & LMFA - UMR CNRS 5509 http://acoustique.ec-lyon.fr Outline of the course A short introduction to

More information

The Polymers Tug Back

The Polymers Tug Back Tugging at Polymers in Turbulent Flow The Polymers Tug Back Jean-Luc Thiffeault http://plasma.ap.columbia.edu/ jeanluc Department of Applied Physics and Applied Mathematics Columbia University Tugging

More information

Fluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14

Fluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14 Fluid Dynamics p.1/14 Fluid Dynamics Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/14 The equations of fluid dynamics are coupled PDEs that form an IVP (hyperbolic). Use

More information

Differential relations for fluid flow

Differential relations for fluid flow Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow

More information

Turbulence: Basic Physics and Engineering Modeling

Turbulence: Basic Physics and Engineering Modeling DEPARTMENT OF ENERGETICS Turbulence: Basic Physics and Engineering Modeling Numerical Heat Transfer Pietro Asinari, PhD Spring 2007, TOP UIC Program: The Master of Science Degree of the University of Illinois

More information

Frédérique Motte and Nicola Schneider (AIM Paris-Saclay, Obs. Bordeaux) Coordinated by Frédérique Motte, Annie Zavagno, and Sylvain Bontemps

Frédérique Motte and Nicola Schneider (AIM Paris-Saclay, Obs. Bordeaux) Coordinated by Frédérique Motte, Annie Zavagno, and Sylvain Bontemps Cloud structure and high-mass star formation in HOBYS, the Herschel imaging survey of OB Young Stellar objects Frédérique Motte and Nicola Schneider (AIM Paris-Saclay, Obs. Bordeaux) http://hobys-herschel.cea.fr

More information

AST242 LECTURE NOTES PART 5

AST242 LECTURE NOTES PART 5 AST242 LECTURE NOTES PART 5 Contents 1. Waves and instabilities 1 1.1. Sound waves compressive waves in 1D 1 2. Jeans Instability 5 3. Stratified Fluid Flows Waves or Instabilities on a Fluid Boundary

More information

CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics

CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS Multimedia Course on Continuum Mechanics Overview Introduction Fluid Mechanics What is a Fluid? Pressure and Pascal s Law Constitutive Equations in Fluids Fluid Models

More information

Predicting natural transition using large eddy simulation

Predicting natural transition using large eddy simulation Center for Turbulence Research Annual Research Briefs 2011 97 Predicting natural transition using large eddy simulation By T. Sayadi AND P. Moin 1. Motivation and objectives Transition has a big impact

More information

! what determines the initial mass function of stars (IMF)? dn /dm

! what determines the initial mass function of stars (IMF)? dn /dm ! what determines the initial mass function of stars (IMF)? dn /dm M ! what determines the initial mass function of stars (IMF)?! What determines the total mass of stars that can form in the cloud? dn

More information

Before we consider two canonical turbulent flows we need a general description of turbulence.

Before we consider two canonical turbulent flows we need a general description of turbulence. Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent

More information

1 Introduction to Governing Equations 2 1a Methodology... 2

1 Introduction to Governing Equations 2 1a Methodology... 2 Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................

More information

fluctuations based on the resolved mean flow

fluctuations based on the resolved mean flow Temperature Fluctuation Scaling in Reacting Boundary Layers M. Pino Martín CTR/NASA Ames, Moffett Field, CA 94035 Graham V. Candler Aerospace Engineering and Mechanics University of Minnesota, Minneapolis,

More information

Amplification of magnetic fields in core collapse

Amplification of magnetic fields in core collapse Amplification of magnetic fields in core collapse Miguel Àngel Aloy Torás, Pablo Cerdá-Durán, Thomas Janka, Ewald Müller, Martin Obergaulinger, Tomasz Rembiasz Universitat de València; Max-Planck-Institut

More information

Turbulence models and excitation of solar oscillation modes

Turbulence models and excitation of solar oscillation modes Center for Turbulence Research Annual Research Briefs Turbulence models and excitation of solar oscillation modes By L. Jacoutot, A. Wray, A. G. Kosovichev AND N. N. Mansour. Motivation and objectives

More information

The Truth about diffusion (in liquids)

The Truth about diffusion (in liquids) The Truth about diffusion (in liquids) Aleksandar Donev Courant Institute, New York University & Eric Vanden-Eijnden, Courant In honor of Berni Julian Alder LLNL, August 20th 2015 A. Donev (CIMS) Diffusion

More information

Magnetic Fields (and Turbulence) in Galaxy Clusters

Magnetic Fields (and Turbulence) in Galaxy Clusters Magnetic Fields (and Turbulence) in Galaxy Clusters Dongsu Ryu (UNIST, Ulsan National Institute of Science and Technology, Korea) with Hyesung Kang (Pusan Nat. U, Korea), Jungyeon Cho (Chungnam Nat. U.),

More information

Accretion Disks I. High Energy Astrophysics: Accretion Disks I 1/60

Accretion Disks I. High Energy Astrophysics: Accretion Disks I 1/60 Accretion Disks I References: Accretion Power in Astrophysics, J. Frank, A. King and D. Raine. High Energy Astrophysics, Vol. 2, M.S. Longair, Cambridge University Press Active Galactic Nuclei, J.H. Krolik,

More information

SAM GEEN (ITA/ZAH HEIDELBERG)

SAM GEEN (ITA/ZAH HEIDELBERG) SAM GEEN (ITA/ZAH HEIDELBERG) WITH PATRICK HENNEBELLE JUAN SOLER AND RALF KLESSEN Credit: Lost Valley Observatory Star formation is self regulating HII HII regions, regions, supernovae supernovae Molecular

More information

Structure formation and scaling relations in the ISM (large scale)

Structure formation and scaling relations in the ISM (large scale) Structure formation and scaling relations in the ISM (large scale) Bruce G. Elmegreen IBM T.J. Watson Research Center Yorktown Heights, NY 10598 USA bge@us.ibm.com February 2017 Gas velocity difference

More information

An evaluation of a conservative fourth order DNS code in turbulent channel flow

An evaluation of a conservative fourth order DNS code in turbulent channel flow Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical

More information

Energy transfer in compressible MHD turbulence

Energy transfer in compressible MHD turbulence Energy transfer in compressible MHD turbulence Philipp Grete Michigan State University in collaboration with Brian O Shea, Kris Beckwith, Wolfram Schmidt and Andrew Christlieb MIPSE Seminar University

More information

AY202a Galaxies & Dynamics Lecture 7: Jeans Law, Virial Theorem Structure of E Galaxies

AY202a Galaxies & Dynamics Lecture 7: Jeans Law, Virial Theorem Structure of E Galaxies AY202a Galaxies & Dynamics Lecture 7: Jeans Law, Virial Theorem Structure of E Galaxies Jean s Law Star/Galaxy Formation is most simply defined as the process of going from hydrostatic equilibrium to gravitational

More information

Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions

Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,

More information

Formation of massive stars: a review

Formation of massive stars: a review Formation of massive stars: a review Patrick Hennebelle Thanks to: Benoît Commerçon, Marc Joos, Andrea Ciardi Gilles Chabrier, Romain Teyssier How massive stars form? -Can we form massive stars in spite

More information

Anisotropic turbulence in rotating magnetoconvection

Anisotropic turbulence in rotating magnetoconvection Anisotropic turbulence in rotating magnetoconvection André Giesecke Astrophysikalisches Institut Potsdam An der Sternwarte 16 14482 Potsdam MHD-Group seminar, 2006 André Giesecke (AIP) Anisotropic turbulence

More information

arxiv:astro-ph/ v2 23 Oct 2001

arxiv:astro-ph/ v2 23 Oct 2001 Kolmogorov Burgers Model for Star Forming Turbulence Stanislav Boldyrev 1 Institute for Theoretical Physics, Santa Barbara, CA 93106 arxiv:astro-ph/0108300v2 23 Oct 2001 ABSTRACT The process of star formation

More information

LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS

LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,

More information

Introduction to Fluid Mechanics

Introduction to Fluid Mechanics Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the

More information

Turbulence Solutions

Turbulence Solutions School of Mechanical, Aerospace & Civil Engineering 3rd Year/MSc Fluids Turbulence Solutions Question 1. Decomposing into mean and fluctuating parts, we write M = M + m and Ũ i = U i + u i a. The transport

More information

A Finite-Element based Navier-Stokes Solver for LES

A Finite-Element based Navier-Stokes Solver for LES A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,

More information

ASTR 610 Theory of Galaxy Formation Lecture 16: Star Formation

ASTR 610 Theory of Galaxy Formation Lecture 16: Star Formation ASTR 610 Theory of Galaxy Formation Lecture 16: Star Formation Frank van den Bosch Yale University, spring 2017 Star Formation In this lecture we discuss the formation of stars. After describing the structure

More information

2. FLUID-FLOW EQUATIONS SPRING 2019

2. FLUID-FLOW EQUATIONS SPRING 2019 2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid

More information

Collapse of Massive Cloud Cores

Collapse of Massive Cloud Cores Collapse of Massive Cloud Cores Robi Banerjee ITA, University of Heidelberg Based on 3D MHD, AMR* Simulations * Adaptive Mesh Refinement Collapse of Hydrostatic Cores Cloud formation by TI Slowly rotating

More information

CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION

CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,

More information

10. Buoyancy-driven flow

10. Buoyancy-driven flow 10. Buoyancy-driven flow For such flows to occur, need: Gravity field Variation of density (note: not the same as variable density!) Simplest case: Viscous flow, incompressible fluid, density-variation

More information

The Effects of Radiative Transfer on Low-Mass Star Formation

The Effects of Radiative Transfer on Low-Mass Star Formation The Effects of Radiative Transfer on Low-Mass Star Formation Stella Offner NSF Fellow, ITC Dense Cores in Dark Clouds Oct 23 2009 Collaborators: Chris McKee (UC Berkeley), Richard Klein (UC Berkeley; LLNL),

More information

Subgrid-Scale Models for Compressible Large-Eddy Simulations

Subgrid-Scale Models for Compressible Large-Eddy Simulations Theoret. Comput. Fluid Dynamics (000 13: 361 376 Theoretical and Computational Fluid Dynamics Springer-Verlag 000 Subgrid-Scale Models for Compressible Large-Eddy Simulations M. Pino Martín Department

More information

FLASH Code Tutorial. part III sink particles & feedback. Robi Banerjee Hamburger Sternwarte

FLASH Code Tutorial. part III sink particles & feedback. Robi Banerjee Hamburger Sternwarte FLASH Code Tutorial part III sink particles & feedback Robi Banerjee Hamburger Sternwarte banerjee@hs.uni-hamburg.de Motivation: modelling of dense regions in collapse simulations, e.g. star formation

More information

Colloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.

Colloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p. Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš

More information

A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries

A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation

More information

An Introduction to Theories of Turbulence. James Glimm Stony Brook University

An Introduction to Theories of Turbulence. James Glimm Stony Brook University An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing

More information

ASP2062 Introduction to Astrophysics

ASP2062 Introduction to Astrophysics School of Physics and Astronomy ASP2062 2015 ASP2062 Introduction to Astrophysics Star formation II Daniel Price Key revision points 1. Star formation is a competition between gravity and pressure 2. Timescale

More information

Max Planck Institut für Plasmaphysik

Max Planck Institut für Plasmaphysik ASDEX Upgrade Max Planck Institut für Plasmaphysik 2D Fluid Turbulence Florian Merz Seminar on Turbulence, 08.09.05 2D turbulence? strictly speaking, there are no two-dimensional flows in nature approximately

More information

Turbulence, kinematics & galaxy structure in star formation in dwarfs. Mordecai-Mark Mac Low Department of Astrophysics

Turbulence, kinematics & galaxy structure in star formation in dwarfs. Mordecai-Mark Mac Low Department of Astrophysics Turbulence, kinematics & galaxy structure in star formation in dwarfs Mordecai-Mark Mac Low Department of Astrophysics Outline Turbulence inhibits star formation, but slowly Interplay between turbulence

More information

The density variance Mach number relation in supersonic turbulence I. Isothermal, magnetized gas

The density variance Mach number relation in supersonic turbulence I. Isothermal, magnetized gas Mon. Not. R. Astron. Soc. 423, 2680 2689 2012 doi:10.1111/j.1365-2966.2012.21075.x The density variance Mach number relation in supersonic turbulence I. Isothermal, magnetized gas F. Z. Molina, 1 S. C.

More information

Diffusive Transport Enhanced by Thermal Velocity Fluctuations

Diffusive Transport Enhanced by Thermal Velocity Fluctuations Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley

More information

Computational Astrophysics

Computational Astrophysics Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture

More information

Getting started: CFD notation

Getting started: CFD notation PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =

More information

Chapter 5. The Differential Forms of the Fundamental Laws

Chapter 5. The Differential Forms of the Fundamental Laws Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations

More information

Turbulent Magnetic Helicity Transport and the Rapid Growth of Large Scale Magnetic Fields

Turbulent Magnetic Helicity Transport and the Rapid Growth of Large Scale Magnetic Fields Turbulent Magnetic Helicity Transport and the Rapid Growth of Large Scale Magnetic Fields Jungyeon Cho Dmitry Shapovalov MWMF Madison, Wisconsin April 2012 The Large Scale Dynamo The accumulation of magnetic

More information

Direct Evidence for Two Fluid Effects in Molecular Clouds. Dinshaw Balsara & David Tilley University of Notre Dame

Direct Evidence for Two Fluid Effects in Molecular Clouds. Dinshaw Balsara & David Tilley University of Notre Dame Direct Evidence for Two Fluid Effects in Molecular Clouds Dinshaw Balsara & David Tilley University of Notre Dame 1 Outline Introduction earliest stages of star formation Theoretical background Magnetically

More information

CHAPTER 4. Basics of Fluid Dynamics

CHAPTER 4. Basics of Fluid Dynamics CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,

More information

Dynamic k-equation Model for Large Eddy Simulation of Compressible Flows. Xiaochuan Chai and Krishnan Mahesh

Dynamic k-equation Model for Large Eddy Simulation of Compressible Flows. Xiaochuan Chai and Krishnan Mahesh 40th Fluid Dynamics Conference and Exhibit 8 June - July 00, Chicago, Illinois AIAA 00-506 Dynamic k-equation Model for Large Eddy Simulation of Compressible Flows Xiaochuan Chai and Krishnan Mahesh University

More information

12.1 Viscous potential flow (VPF)

12.1 Viscous potential flow (VPF) 1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow

More information

An evaluation of LES for jet noise prediction

An evaluation of LES for jet noise prediction Center for Turbulence Research Proceedings of the Summer Program 2002 5 An evaluation of LES for jet noise prediction By B. Rembold, J. B. Freund AND M. Wang Large-eddy simulation (LES) is an attractive

More information