III.H Zeroth Order Hydrodynamics
|
|
- Justina Murphy
- 5 years ago
- Views:
Transcription
1 III.H Zeroth Order Hydrodynamics As a first approximation, we shall assume that in local equilibrium, the density f 1 at each point in space can be represented as in eq.(iii.56), i.e. [ ] p m q, t)) f n( q, t) 1 ( p, q, t) = exp ( u(. (III.9) / (πmk B T( q, t)) mk B T( q, t) The choice of parameters clearly enforces d p f 1 = n, and p/m = u, as required. Average values are easily calculated from this Gaussian weight; in particular c α c β = k B T δαβ, m (III.94) leading to P and ε = αβ = nk B Tδ αβ, k B T. (III.95) Since the density f 1 is even in c, all odd expectation values vanish, and in particular h =. (III.96) The conservation laws in this approximation take the simple forms D t n = n α u α 1 md t u α = F α α (nk B T) n. (III.97) D t T = T α u α In the above expression, we have introduced the material derivative D t [ t + u β β ], (III.98) which measures the time variations of any quantity as it moves along the stream-lines set up by the average velocity field u. By combining the first and third equations, it is easy to get D t ln nt / =. (III.99) ( The quantity ln ) nt / is like a local entropy for the gas (see eq.(iii.67)), which according to the above equation is not changed along stream-lines. The zeroth order hydrodynamics thus predicts that the gas flow is adiabatic. This prevents the local equilibrium solution 68
2 of eq.(iii.9) from reaching a true global equilibrium form which necessitates an increase in entropy. To demonstrate that eqs.(iii.97) do not describe a satisfactory approach to equilibrium, examine the evolution of small deformations about a stationary ( u = ) state, in a uniform box (F = ), by setting { n( q, t) =n + ν( q, t). (III.1) T( q, t) =T + θ( q, t) We shall next expand eqs.(iii.97) to first order in the deviations (ν, θ, u). Note that to lowest order, D t = t +O(u), leading to the linearized zeroth order hydrodynamic equations t ν = n α u α k B T m t u α = α ν k B α θ n. (III.11) t θ = T α u α Normal modes of the system are obtained by Fourier transformations, [ ] A k, ω = d q dt exp i k q ωt A ( q, t), (III.1) where A stands for any of the three fields (ν, θ, u). The natural vibration frequencies are solutions to the matrix equation ν nk β ν ω u α = k B T k δ αβ k β B δ αβ k β u β. mn m θ Tk β θ (III.1) It is easy to check that this equation has the following modes, the first three with zero frequency: (a) Two modes describe shear flows in a uniform (n = n) and isothermal (T = T) fluid, in which the velocity varies along a direction normal to its orientation (e.g. u = f(x, t)ŷ). In terms of Fourier modes k u T ( k ) =, indicating transverse flows that are not relaxed in this zeroth order approximation. (b) A third zero frequency mode describes a stationary fluid with uniform pressure P = nk B T. While n and T may vary across space, their product is constant, insuring 69
3 that the fluid will not start moving due to pressure variations. The corresponding eigenvector of eq.(iii.1) is n v e =. (III.14) T (c) Finally, the longitudinal velocity ( u l k) combines with density and temperature variations in eigenmodes of the form n k v l = ω( k), with ω( k) = ±v l k, (III.15) T k where v l = 5 k B T, m (III.16) is the longitudinal sound velocity. Note that the density and temperature variations ( in this mode are adiabatic, i.e. the local entropy (proportional to ln ) nt / ) is left unchanged. We thus find that none of the conserved quantities relaxes to equilibrium in the zeroth order approximation. Shear flow and entropy modes persist forever, while the two sound modes have undamped oscillations. This is a deficiency of the zeroth order approximation which is removed by finding a better solution to the Boltzmann equation. III.I First Order Hydrodynamics While f 1 ( p, q, t) of eq.(iii.9) does set the right hand side of the Boltzmann equation to zero, it is not a full solution, as the left hand side causes its form to vary. The left hand side is a linear differential operator, which using the various notations introduced in the previous sections, can be written as [ ] [ ] p α F α L[f] t + α + F α f = D t + c α α + f. (III.17) m p α m c α It is simpler to examine the effect of L on lnf 1. which can be written as ln f 1 = ln nt / mc ln(πmkb ). (III.18) k B T 7
4 Using the relation (c /) = c β c β = c β u β, we get [ L ] ln f1 =D t ln nt / mc m + D t T + cα D t u α k B T k B T (III.19) α n α T mc m F α c α +c α + c α α T + c α c β α u β. n T k B T k B T k B T If the fields n, T, and u α, satisfy the zeroth order hydrodynamic eqs.(iii.97), we can simplify the above equation to [ L ln ] mc [ ] Fα α n α T α n α T F α f1 = α u α + c α + k B T k B T n T n T k B T mc m + c α α T + c α c β u αβ k B T k B T m δ αβ mc 5 cα = c α c β c u αβ + α T. k B T k B T T and linearize the collision operator as C [f 1, f 1 ] = d p d b v1 v f 1 (p 1 )f 1 (p ) [g(p 1 ) + g(p ) g(p 1 ) g(p )] (III.11) The characteristic time scale τ U for L is extrinsic, and can be made much larger than τ. The zeroth order result is thus exact in the limit (τ /τ U ) ; and corrections can be constructed in a perturbation series in (τ /τ U ). To this purpose, we set f 1 = f 1 (1 + g), f 1 (p 1 )C L [g]. (III.111) While linear, the above integral operator is still difficult to manipulate in general. As a first approximation, and noting its characteristic magnitude, we set C L [g] g. τ (III.11) This is known as the single collision time approximation, and from the linearized Boltzmann equation L[f 1 ] = f 1 C L [g], we obtain 1 g = τ f L [f [ 1] τ L ] ln f1, (III.11) 1 where we have kept only the leading term. Thus the first order solution is given by (using eq.(iii.11)) [ ] f 1 1 (p, q, t) = f τ µ m δ αβ mc 5 cα 1 ( p, q, t) 1 c α c β c u αβ τ K α T, k B T k B T T (III.114) 71
5 where τ µ = τ K = τ in the single collision time approximation. However, in writing the above equation, we have anticipated the possibility of τ µ = τ K which arises in more sophisticated treatments (although both times are still of order of τ ). 1 It is easy to check that d pf 1 = d pf 1 = n, and thus various local expectation values are calculated to first order as O 1 = 1 d p Of 1 (1 + g) = O + go. (III.115) n The calculation of averages over products of c α s, distributed according to the Gaussian weight of f 1, is greatly simplified by the use of Wick s theorem, which states that expectation value of the product is the sum over all possible products of paired expectation values, for example k B T c α c β c γ c δ = (δ αβ δ γδ + δ αγ δ βδ + δ αδ δ βγ ). (III.116) m (Expectation values involving a product of an odd number of c α s are zero by symmetry.) Using this result, it is easy to verify that pα 1 β T mc 5 = u α τ K c α c β = u α. (III.117) m T k B T The pressure tensor at first order is given by [ ] 1 = nm c α c β τ µ m δ µν αβ =nm c α c β cα c β c µ c ν c u µν k B T δ αβ =nk B Tδ αβ nk B Tτ µ u αβ u γγ. P 1 (III.118) (Using the above result, we can further verify that ε 1 = Finally, the heat flux is given by h 1 mc 1 nmτ K β T mc 5 α =n c α = c α c β c T k B T 5 nkb Tτ K = α T. m mc 1 / = kb T/, as before.) (III.119) At this order, we find that spatial variations in temperature generate a heat flow that tends to smooth them out, while shear flows are opposed by the off-diagonal terms in the pressure tensor. These effects are sufficient to cause relaxation to equilibrium, as can be seen by examining the modified behavior of the modes discussed previously. 7
6 (a) The pressure tensor now has an off diagonal term 1 P α= β = nk B Tτ µ u αβ µ ( α u β + β u α ), (III.1) where µ nk B Tτ µ is the viscosity coefficient. A shearing of the fluid (e.g. described by a velocity u y (x, t)) now leads to a a viscous force that opposes it (proportional to µ x u y ), causing its diffusive relaxation as discussed below. (b) Similarly, a temperature gradient leads to a heat flux h = K T, (III.11) where K = (5nkB Tτ K)/(m) is the coefficient of thermal conductivity of the gas. If the gas is at rest ( u =, and uniform P = nk B T), variations in temperature now satisfy n t ε = nk B t T = α ( K α T), t T = K T. (III.1) nk B This is the Fourier equation and shows that temperature variations relax by diffusion. We can discuss the behavior of all the modes by linearizing the equations of motion. The first order contribution to D t u α t u α is δ 1 1 ( t u α ) β δ 1 µ 1 P αβ α β + δ αβ γ γ u β, (III.1) mn mn where µ nk B Tτ µ. Similarly, the correction for D t T t θ, is given by δ 1 K ( t θ) α h α α α θ, k B n k B n (III.14) with K = (5nkB Tτ K )/(m). After Fourier transformation, the matrix equation (III.1) is modified to nδ αβ k β ν ν ω u α k = B T µ δ αβ k β i k k δ αβ + α k β k B δ αβ k mn mn m β u β. (III.15) θ Tδ αβ k β i θ Kk k B n We can ask how the normal mode frequencies calculated in the zeroth order approximation are modified at this order. It is simple to verify that the transverse (shear) normal models ( k u T = ) now have a frequency ω T = i µ k. mn (III.16) 7
7 The imaginary frequency implies that these modes are damped over a characteristic time τ T (k) 1/ ω T (λ) /(τ µ v ), where λ is the corresponding wavelength, and v k B T/m is a typical gas particle velocity. We see that the characteristic time scales grow as the square of the wavelength, which is characteristic of diffusive processes. In the remaining normal modes the velocity is parallel to k, and eq.(iii.15) reduces to ν nk ν ω u l = k B T u β mn k i 4µk k B mn m k. (III.17) θ θ Tk ikk k B n The determinant of the dynamical matrix is the product of the three eigen-frequencies, and to lowest order is given by Kk kbtk det(m) = i nk + O(τ ). (III.18) k B n mn At zeroth order the two sound modes have ω ± (k) = ±v l k, and hence the frequency of the isobaric mode is ω 1 det(m) Kk e(k) = i + O(τ ). (III.19) v k 5k B n l At first order, the longitudinal sound modes also turn into damped oscillations with frequencies ω± 1 (k) = ±v l k iγ. The simplest way to obtain the decay rates is to note that the trace of the dynamical matrix is equal to the sum of the eigenvalues, and hence 1 µ K ω ± (k) = ±v l k ik + + O(τ ). (III.1) mn 15k B n The damping of all normal modes guarantees the, albeit slow, approach of the gas to its final uniform and stationary equilibrium state. 74
8 MIT OpenCourseWare 8. Statistical Mechanics I: Statistical Mechanics of Particles Fall 1 For information about citing these materials or our Terms of Use, visit:
III.H Zeroth Order Hydrodynamics
III.H Zeroth Order Hydrodynaics As a first approxiation, we shall assue that in local equilibriu, the density f 1 at each point in space can be represented as in eq.iii.56, i.e. f 0 1 p, q, t = n q, t
More information8.333: Statistical Mechanics I Problem Set # 3 Due: 10/18/13. Kinetic Theory
8.333: Statistical Mechanics I Problem Set # 3 Due: 10/18/13 1. Poisson Brackets: Kinetic Theory (a) Show that for observable O(p(µ),q(µ)), do/dt = {O,H}, along the time trajectory of any micro state µ,
More information8.333: Statistical Mechanics I Fall 2007 Test 2 Review Problems
8.333: Statistical Mechanics I Fall 007 Test Review Problems The second in-class test will take place on Wednesday 10/4/07 from :30 to 4:00 pm. There will be a recitation with test review on Monday 10//07.
More informationAccurate representation of velocity space using truncated Hermite expansions.
Accurate representation of velocity space using truncated Hermite expansions. Joseph Parker Oxford Centre for Collaborative Applied Mathematics Mathematical Institute, University of Oxford Wolfgang Pauli
More informationI. Collective Behavior, From Particles to Fields
I. Collective Behavior, From Particles to Fields I.A Introduction The object of the first part of this course was to introduce the principles of statistical mechanics which provide a bridge between the
More informationChapter 6: The Rouse Model. The Bead (friction factor) and Spring (Gaussian entropy) Molecular Model:
G. R. Strobl, Chapter 6 "The Physics of Polymers, 2'nd Ed." Springer, NY, (1997). R. B. Bird, R. C. Armstrong, O. Hassager, "Dynamics of Polymeric Liquids", Vol. 2, John Wiley and Sons (1977). M. Doi,
More informationVIII.B Equilibrium Dynamics of a Field
VIII.B Equilibrium Dynamics of a Field The next step is to generalize the Langevin formalism to a collection of degrees of freedom, most conveniently described by a continuous field. Let us consider the
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More informationLattice Boltzmann Method
3 Lattice Boltzmann Method 3.1 Introduction The lattice Boltzmann method is a discrete computational method based upon the lattice gas automata - a simplified, fictitious molecular model. It consists of
More informationPolymer Dynamics and Rheology
Polymer Dynamics and Rheology 1 Polymer Dynamics and Rheology Brownian motion Harmonic Oscillator Damped harmonic oscillator Elastic dumbbell model Boltzmann superposition principle Rubber elasticity and
More informationMagnetohydrodynamic waves in a plasma
Department of Physics Seminar 1b Magnetohydrodynamic waves in a plasma Author: Janez Kokalj Advisor: prof. dr. Tomaž Gyergyek Petelinje, April 2016 Abstract Plasma can sustain different wave phenomena.
More information8.333: Statistical Mechanics I Fall 2005 Mid-term Quiz Review Problems
8.333: Statistical Mechanics I Fall 005 Mid-term Quiz 0 Review Problems The Mid-term quiz will take place on Monday 10/4/05 in room 1-190 from :30 to 4:00 pm. There will be a recitation with quiz review
More information1.3 Molecular Level Presentation
1.3.1 Introduction A molecule is the smallest chemical unit of a substance that is capable of stable, independent existence. Not all substances are composed of molecules. Some substances are composed of
More informationarxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
More informationFundamentals of wave kinetic theory
Fundamentals of wave kinetic theory Introduction to the subject Perturbation theory of electrostatic fluctuations Landau damping - mathematics Physics of Landau damping Unmagnetized plasma waves The plasma
More informationHydrodynamics. 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 informationJ10M.1 - Rod on a Rail (M93M.2)
Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight
More information8.334: Statistical Mechanics II Spring 2014 Test 2 Review Problems
8.334: Statistical Mechanics II Spring 014 Test Review Problems The test is closed book, but if you wish you may bring a one-sided sheet of formulas. The intent of this sheet is as a reminder of important
More informationLecture 2 Supplementary Notes: Derivation of the Phase Equation
Lecture 2 Supplementary Notes: Derivation of the Phase Equation Michael Cross, 25 Derivation from Amplitude Equation Near threshold the phase reduces to the phase of the complex amplitude, and the phase
More information14. Energy transport.
Phys780: Plasma Physics Lecture 14. Energy transport. 1 14. Energy transport. Chapman-Enskog theory. ([8], p.51-75) We derive macroscopic properties of plasma by calculating moments of the kinetic equation
More informationSolution set #7 Physics 571 Tuesday 3/17/2014. p 1. p 2. Figure 1: Muon production (e + e µ + µ ) γ ν /p 2
Solution set #7 Physics 571 Tuesday 3/17/2014 μ 1. The amplitude is Figure 1: Muon production ( e µ + µ ) it = ie2 s (v 2γ µ u 1 )(u 1 γ µ v 2 ), (1) so the spin averaged squared amplitude is T 2 = e4
More informationPhonons and lattice dynamics
Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationIn this chapter we will discuss the effect of shear viscosity on evolution of fluid, p T
Chapter 3 Shear viscous evolution In this chapter we will discuss the effect of shear viscosity on evolution of fluid, p T spectra, and elliptic flow (v ) of pions using a +1D relativistic viscous hydrodynamics
More informationMIT (Spring 2014)
18.311 MIT (Spring 014) Rodolfo R. Rosales May 6, 014. Problem Set # 08. Due: Last day of lectures. IMPORTANT: Turn in the regular and the special problems stapled in two SEPARATE packages. Print your
More informationTheory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus. K. C. Shaing
Theory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus K. C. Shaing Plasma and Space Science Center, and ISAPS, National Cheng Kung University, Tainan, Taiwan 70101, Republic
More informationFlow and Transport. c(s, t)s ds,
Flow and Transport 1. The Transport Equation We shall describe the transport of a dissolved chemical by water that is traveling with uniform velocity ν through a long thin tube G with uniform cross section
More informationLecture 20: Modes in a crystal and continuum
Physics 16a, Caltech 11 December, 218 The material in this lecture (which we had no time for) will NOT be on the Final exam. Lecture 2: Modes in a crystal and continuum The vibrational modes of the periodic
More informationIV. Classical Statistical Mechanics
IV. Classical Statistical Mechanics IV.A General Definitions Statistical Mechanics is a probabilistic approach to equilibrium macroscopic properties of large numbers of degrees of freedom. As discussed
More informationMODERN OPTICS. P47 Optics: Unit 9
MODERN OPTICS P47 Optics: Unit 9 Course Outline Unit 1: Electromagnetic Waves Unit 2: Interaction with Matter Unit 3: Geometric Optics Unit 4: Superposition of Waves Unit 5: Polarization Unit 6: Interference
More informationarxiv: v1 [nucl-th] 9 Jun 2008
Dissipative effects from transport and viscous hydrodynamics arxiv:0806.1367v1 [nucl-th] 9 Jun 2008 1. Introduction Denes Molnar 1,2 and Pasi Huovinen 1 1 Purdue University, Physics Department, 525 Northwestern
More informationI. INTRODUCTION AND HISTORICAL PERSPECTIVE
I. INTRODUCTION AND HISTORICAL PERSPECTIVE A. Failures of Classical Physics At the end of the 19th century, physics was described via two main approaches. Matter was described by Newton s laws while radiation
More information20. Alfven waves. ([3], p ; [1], p ; Chen, Sec.4.18, p ) We have considered two types of waves in plasma:
Phys780: Plasma Physics Lecture 20. Alfven Waves. 1 20. Alfven waves ([3], p.233-239; [1], p.202-237; Chen, Sec.4.18, p.136-144) We have considered two types of waves in plasma: 1. electrostatic Langmuir
More informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationThe Interaction of Light and Matter: α and n
The Interaction of Light and Matter: α and n The interaction of light and matter is what makes life interesting. Everything we see is the result of this interaction. Why is light absorbed or transmitted
More informationisocurvature modes Since there are two degrees of freedom in
isocurvature modes Since there are two degrees of freedom in the matter-radiation perturbation, there must be a second independent perturbation mode to complement the adiabatic solution. This clearly must
More information8 Ecosystem stability
8 Ecosystem stability References: May [47], Strogatz [48]. In these lectures we consider models of populations, with an emphasis on the conditions for stability and instability. 8.1 Dynamics of a single
More informationAdvanced Partial Differential Equations with Applications
MIT OpenCourseWare http://ocw.mit.edu 8.36 Advanced Partial Differential Equations with Applications Fall 9 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationTowards new relativistic hydrodynamcis from AdS/CFT
Towards new relativistic hydrodynamcis from AdS/CFT Michael Lublinsky Stony Brook with Edward Shuryak QGP is Deconfined QGP is strongly coupled (sqgp) behaves almost like a perfect liquid (Navier-Stokes
More informationCollective behavior, from particles to fields
978-0-51-87341-3 - Statistical Physics of Fields 1 Collective behavior, from particles to fields 1.1 Introduction One of the most successful aspects of physics in the twentieth century was revealing the
More informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.306 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Monday March 12, 2018. Turn it in (by 3PM) at the Math.
More informationIntroduction to Relaxation Theory James Keeler
EUROMAR Zürich, 24 Introduction to Relaxation Theory James Keeler University of Cambridge Department of Chemistry What is relaxation? Why might it be interesting? relaxation is the process which drives
More informationFrequency spectra at large wavenumbers in two-dimensional Hasegawa-Wakatani turbulence
Frequency spectra at large wavenumbers in two-dimensional Hasegawa-Wakatani turbulence Juhyung Kim and P. W. Terry Department of Physics, University of Wisconsin-Madison October 30th, 2012 2012 APS-DPP
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationarxiv: v1 [nucl-th] 18 Apr 2019
Prepared for submission to JHEP Emergence of hydrodynamical behavior in expanding quark-gluon plasmas arxiv:1904.08677v1 [nucl-th] 18 Apr 2019 Jean-Paul Blaizot a and Li Yan b a Institut de Physique Théorique,
More informationLecture 7. 1 Wavepackets and Uncertainty 1. 2 Wavepacket Shape Changes 4. 3 Time evolution of a free wave packet 6. 1 Φ(k)e ikx dk. (1.
Lecture 7 B. Zwiebach February 8, 06 Contents Wavepackets and Uncertainty Wavepacket Shape Changes 4 3 Time evolution of a free wave packet 6 Wavepackets and Uncertainty A wavepacket is a superposition
More informationFundamentals 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 informationD I S P E R S I O N AND D A M P I N G OF SECOND SOUND IN ) Abstract
Physics Vol. 3, No. 4. pp. 221-229, 1967. Physics Publishing Co. Printed in Great Britain. D I S P E R S I O N AND D A M P I N G OF SECOND SOUND IN N O N - I S O T R O P I C SOLIDS P H I L I P C. K WOK
More informationFrom the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1)
Lecture 9 9. Administration None. 9. Continuation of energy equation From the last time, we ended with an expression for the energy equation ρ D (e + ) u = ρg u + (τ u) q (9.) Where ρg u changes in potential
More informationI.I Stability Conditions
I.I Stability Conditions The conditions derived in section I.G are similar to the well known requirements for mechanical stability. A particle moving in an eternal potential U settles to a stable equilibrium
More information6 Inner Product and Hilbert Spaces
6 Inner Product and Hilbert Spaces 6. Motivation Of the different p-norms on R n, p = 2 is special. This is because the 2-norm (λ, λ 2,..., λ n ) 2 = λ 2 + λ2 2 + + λ2 n comes from the an inner product
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More informationNucleons in Nuclei: Interactions, Geometry, Symmetries
Nucleons in Nuclei: Interactions, Geometry, Symmetries Jerzy DUDEK Department of Subatomic Research, CNRS/IN 2 P 3 and University of Strasbourg, F-67037 Strasbourg, FRANCE September 28, 2010 Mathematial
More informationDamping of MHD waves in the solar partially ionized plasmas
Damping of MHD waves in the solar partially ionized plasmas M. L. Khodachenko Space Research Institute, Austrian Academy of Sciences, Graz, Austria MHD waves on the Sun Magnetic field plays the key role
More informationCHAPTER 8 ENTROPY GENERATION AND TRANSPORT
CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.1 CONVECTIVE FORM OF THE GIBBS EQUATION In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationLINEAR RESPONSE THEORY
MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior
More informationNon-Equilibrium Kinetics and Transport Processes in a Hypersonic Flow of CO 2 /CO/O 2 /C/O Mixture
Non-Equilibrium Kinetics and Transport Processes in a Hypersonic Flow of CO 2 /CO/O 2 /C/O Mixture E.V. Kustova, E.A. Nagnibeda, Yu.D. Shevelev and N.G. Syzranova Department of Mathematics and Mechanics,
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More informationPhysics 451/551 Theoretical Mechanics. G. A. Krafft Old Dominion University Jefferson Lab Lecture 18
Physics 451/551 Theoretical Mechanics G. A. Krafft Old Dominion University Jefferson Lab Lecture 18 Theoretical Mechanics Fall 18 Properties of Sound Sound Waves Requires medium for propagation Mainly
More informationNote the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size
L Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size 0.01L or smaller are subject to substantial viscous
More informationGravitational radiation
Lecture 28: Gravitational radiation Gravitational radiation Reading: Ohanian and Ruffini, Gravitation and Spacetime, 2nd ed., Ch. 5. Gravitational equations in empty space The linearized field equations
More informationConstraints on Fluid Dynamics From Equilibrium Partition Func
Constraints on Fluid Dynamics From Equilibrium Partition Function Nabamita Banerjee Nikhef, Amsterdam LPTENS, Paris 1203.3544, 1206.6499 J. Bhattacharya, S. Jain, S. Bhattacharyya, S. Minwalla, T. Sharma,
More informationSummary: Thermodynamic Potentials and Conditions of Equilibrium
Summary: Thermodynamic Potentials and Conditions of Equilibrium Isolated system: E, V, {N} controlled Entropy, S(E,V{N}) = maximum Thermal contact: T, V, {N} controlled Helmholtz free energy, F(T,V,{N})
More information6.730 Physics for Solid State Applications
6.730 Physics for Solid State Applications Lecture 5: Specific Heat of Lattice Waves Outline Review Lecture 4 3-D Elastic Continuum 3-D Lattice Waves Lattice Density of Modes Specific Heat of Lattice Specific
More informationQuasi-Harmonic Theory of Thermal Expansion
Chapter 5 Quasi-Harmonic Theory of Thermal Expansion 5.1 Introduction The quasi-harmonic approximation is a computationally efficient method for evaluating thermal properties of materials. Planes and Manosa
More informationFrom laser cooling to BEC First experiments of superfluid hydrodynamics
From laser cooling to BEC First experiments of superfluid hydrodynamics Alice Sinatra Quantum Fluids course - Complement 1 2013-2014 Plan 1 COOLING AND TRAPPING 2 CONDENSATION 3 NON-LINEAR PHYSICS AND
More informationF r (t) = 0, (4.2) F (t) = r= r a (t)
Chapter 4 Stochastic Equations 4.1 Langevin equation To explain the idea of introduction of the Langevin equation, let us consider the concrete example taken from surface physics: motion of an atom (adatom)
More informationWave Equation in One Dimension: Vibrating Strings and Pressure Waves
BENG 1: Mathematical Methods in Bioengineering Lecture 19 Wave Equation in One Dimension: Vibrating Strings and Pressure Waves References Haberman APDE, Ch. 4 and Ch. 1. http://en.wikipedia.org/wiki/wave_equation
More informationAnswers to Problem Set # 01, MIT (Winter-Spring 2018)
Answers to Problem Set # 01, 18.306 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Contents 1 Nonlinear solvable ODEs 2 1.1 Statement:
More informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
More informationPAPER 52 GENERAL RELATIVITY
MATHEMATICAL TRIPOS Part III Monday, 1 June, 2015 9:00 am to 12:00 pm PAPER 52 GENERAL RELATIVITY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight.
More informationLongitudinal Beam Dynamics
Longitudinal Beam Dynamics Shahin Sanaye Hajari School of Particles and Accelerators, Institute For Research in Fundamental Science (IPM), Tehran, Iran IPM Linac workshop, Bahman 28-30, 1396 Contents 1.
More informationAnisotropic fluid dynamics. Thomas Schaefer, North Carolina State University
Anisotropic fluid dynamics Thomas Schaefer, North Carolina State University Outline We wish to extract the properties of nearly perfect (low viscosity) fluids from experiments with trapped gases, colliding
More informationarxiv: v1 [cond-mat.stat-mech] 26 Jan 2008
Long Time Tail of the Velocity Autocorrelation Function in a Two-Dimensional Moderately Dense Hard Disk Fluid Masaharu Isobe Graduate School of Engineering, Nagoya Institute of Technology, Nagoya 466-8555,
More informationTutorial 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 informationExercise 1: Vertical structure of the lower troposphere
EARTH SCIENCES SCIENTIFIC BACKGROUND ASSESSMENT Exercise 1: Vertical structure of the lower troposphere In this exercise we will describe the vertical thermal structure of the Earth atmosphere in its lower
More informationThe existence of Burnett coefficients in the periodic Lorentz gas
The existence of Burnett coefficients in the periodic Lorentz gas N. I. Chernov and C. P. Dettmann September 14, 2006 Abstract The linear super-burnett coefficient gives corrections to the diffusion equation
More informationKinetic Models for Granular Flow
Journal of Statistical Physics, Vol. 97, Nos. 12, 1999 Kinetic Models for Granular Flow J. Javier Brey, 1 James W. Dufty, 2 and Andre s Santos 3 Received September 28, 1998; final March 29, 1999 The generalization
More informationLattice Boltzmann approach to liquid - vapour separation
Lattice Boltzmann approach to liquid - vapour separation T.Biciușcă 1,2, A.Cristea 1, A.Lamura 3, G.Gonnella 4, V.Sofonea 1 1 Centre for Fundamental and Advanced Technical Research, Romanian Academy Bd.
More informationPhysics 607 Exam 2. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2
Physics 607 Exam Please be well-organized, and show all significant steps clearly in all problems. You are graded on your work, so please do not just write down answers with no explanation! Do all your
More informationStructure formation. Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München
Structure formation Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München Structure formation... Random density fluctuations, grow via gravitational instability galaxies, clusters, etc. Initial perturbations
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II External Flows 1 Introduction In this chapter we will consider several fundamental flows, namely: the flat plate, the cylinder, the sphere, several other body shapes, and banks
More information(1) Transition from one to another laminar flow. (a) Thermal instability: Bernard Problem
Professor Fred Stern Fall 2014 1 Chapter 6: Viscous Flow in Ducts 6.2 Stability and Transition Stability: can a physical state withstand a disturbance and still return to its original state. In fluid mechanics,
More information2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasma physics Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas Transverse waves
More informationRelaxation. Ravinder Reddy
Relaxation Ravinder Reddy Relaxation What is nuclear spin relaxation? What causes it? Effect on spectral line width Field dependence Mechanisms Thermal equilibrium ~10-6 spins leads to NMR signal! T1 Spin-lattice
More information8.21 The Physics of Energy Fall 2009
MIT OpenCourseWare http://ocw.mit.edu 8.21 The Physics of Energy Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 8.21 Lecture 9 Heat Engines
More informationMACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM. Contents AND BOLTZMANN ENTROPY. 1 Macroscopic Variables 3. 2 Local quantities and Hydrodynamics fields 4
MACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM AND BOLTZMANN ENTROPY Contents 1 Macroscopic Variables 3 2 Local quantities and Hydrodynamics fields 4 3 Coarse-graining 6 4 Thermal equilibrium 9 5 Two systems
More informationFluid Dynamics and Balance Equations for Reacting Flows
Fluid Dynamics and Balance Equations for Reacting Flows Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Balance Equations Basics: equations of continuum mechanics balance equations for mass and
More informationA Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit
A Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit David Vener Department of Mathematics, MIT May 5, 3 Introduction In 8.366, we discussed the relationship
More informationMath 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations
Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The
More informationBulk and shear viscosities for the Gribov-Zwanziger plasma
EPJ Web of Conferences will be set by the publisher DOI: will be set by the publisher c Owned by the authors, published by EDP Sciences, 215 Bulk and shear viscosities for the Gribov-Zwanziger plasma Wojciech
More informationFundamental equations of relativistic fluid dynamics
CHAPTER VI Fundamental equations of relativistic fluid dynamics When the energy density becomes large as may happen for instance in compact astrophysical objects, in the early Universe, or in high-energy
More informationThermal and Statistical Physics Department Exam Last updated November 4, L π
Thermal and Statistical Physics Department Exam Last updated November 4, 013 1. a. Define the chemical potential µ. Show that two systems are in diffusive equilibrium if µ 1 =µ. You may start with F =
More information(Super) Fluid Dynamics. Thomas Schaefer, North Carolina State University
(Super) Fluid Dynamics Thomas Schaefer, North Carolina State University Hydrodynamics Hydrodynamics (undergraduate version): Newton s law for continuous, deformable media. Fluids: Gases, liquids, plasmas,...
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More information