The propagation of chaos for a rarefied gas of hard spheres

Size: px
Start display at page:

Download "The propagation of chaos for a rarefied gas of hard spheres"

Transcription

1 The propagation of chaos for a rarefied gas of hard spheres Ryan Denlinger 1 1 University of Texas at Austin 35th Annual Western States Mathematical Physics Meeting Caltech February 13, 2017 Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 1 / 30

2 Outline 1 Introduction & Motivation 2 Lanford s theorem 3 Dispersive inequalities for hard spheres 4 Ideas behind the proof 5 References Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 2 / 30

3 Introduction & Motivation D. Hilbert, J.C. Maxwell and L. Boltzmann (Wikipedia) Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 3 / 30

4 Introduction & Motivation Hilbert s Sixth Problem To treat in the same manner, by means of axioms, those physical sciences in which already today mathematics plays an important part; in the first rank are the theory of probabilities and mechanics. (Wikipedia) Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 4 / 30

5 Introduction & Motivation Hilbert s Sixth Problem To treat in the same manner, by means of axioms, those physical sciences in which already today mathematics plays an important part; in the first rank are the theory of probabilities and mechanics. (Wikipedia) Certain areas of continuum mechanics require dissipation: compressible fluid dynamics (shock formation), incompressible avier-stokes (viscosity), etc. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 4 / 30

6 Introduction & Motivation Hilbert s Sixth Problem To treat in the same manner, by means of axioms, those physical sciences in which already today mathematics plays an important part; in the first rank are the theory of probabilities and mechanics. (Wikipedia) Certain areas of continuum mechanics require dissipation: compressible fluid dynamics (shock formation), incompressible avier-stokes (viscosity), etc. In any case kinetic energy is converted into thermal energy in an irreversible manner until some steady state is achieved (hence turning back the clock would violate the 2nd law of thermodynamics). Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 4 / 30

7 Introduction & Motivation Hilbert s Sixth Problem To treat in the same manner, by means of axioms, those physical sciences in which already today mathematics plays an important part; in the first rank are the theory of probabilities and mechanics. (Wikipedia) Certain areas of continuum mechanics require dissipation: compressible fluid dynamics (shock formation), incompressible avier-stokes (viscosity), etc. In any case kinetic energy is converted into thermal energy in an irreversible manner until some steady state is achieved (hence turning back the clock would violate the 2nd law of thermodynamics). The laws of physics are time-reversible! How should we resolve this apparent paradox? Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 4 / 30

8 Introduction & Motivation Physical systems contain many degrees of freedom, e.g. = Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 5 / 30

9 Introduction & Motivation Physical systems contain many degrees of freedom, e.g. = We can observe O(1) degrees of freedom. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 5 / 30

10 Introduction & Motivation Physical systems contain many degrees of freedom, e.g. = We can observe O(1) degrees of freedom. Boltzmann claims that - for a macroscopic system in equilibrium - all microscopic configurations that are consistent with observations are equally likely. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 5 / 30

11 Introduction & Motivation Physical systems contain many degrees of freedom, e.g. = We can observe O(1) degrees of freedom. Boltzmann claims that - for a macroscopic system in equilibrium - all microscopic configurations that are consistent with observations are equally likely. Corollary: dynamical processes take us from less probable states to more probable states. This is quantified using entropy, i.e. the log of the number of microstates consistent with our observations. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 5 / 30

12 Introduction & Motivation Physical systems contain many degrees of freedom, e.g. = We can observe O(1) degrees of freedom. Boltzmann claims that - for a macroscopic system in equilibrium - all microscopic configurations that are consistent with observations are equally likely. Corollary: dynamical processes take us from less probable states to more probable states. This is quantified using entropy, i.e. the log of the number of microstates consistent with our observations. Heuristically, this idea explains irreversibility (toy example: free expansion of a perfect gas; the reverse process is never observed). Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 5 / 30

13 Introduction & Motivation Physical systems contain many degrees of freedom, e.g. = We can observe O(1) degrees of freedom. Boltzmann claims that - for a macroscopic system in equilibrium - all microscopic configurations that are consistent with observations are equally likely. Corollary: dynamical processes take us from less probable states to more probable states. This is quantified using entropy, i.e. the log of the number of microstates consistent with our observations. Heuristically, this idea explains irreversibility (toy example: free expansion of a perfect gas; the reverse process is never observed). But how do you justify this conclusion from ewton s laws? ot obvious! Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 5 / 30

14 Introduction & Motivation Physical systems contain many degrees of freedom, e.g. = We can observe O(1) degrees of freedom. Boltzmann claims that - for a macroscopic system in equilibrium - all microscopic configurations that are consistent with observations are equally likely. Corollary: dynamical processes take us from less probable states to more probable states. This is quantified using entropy, i.e. the log of the number of microstates consistent with our observations. Heuristically, this idea explains irreversibility (toy example: free expansion of a perfect gas; the reverse process is never observed). But how do you justify this conclusion from ewton s laws? ot obvious! Ergodicity (if it holds) gives some insight but does not necessarily help us explain, e.g., compressible fluid flows. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 5 / 30

15 Lanford s theorem Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 6 / 30

16 Lanford s theorem Consider identical hard spheres in R d, of diameter ε > 0 (small), with the Boltzmann-Grad scaling ε d 1 = l 1 (l > 0 fixed). The spheres undergo free transport, interrupted by elastic collisions precisely when x i x j = ε for some 1 i < j. Elastic collision of hard spheres (ε = 1 in the figure) Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 7 / 30

17 Lanford s theorem Assume that the spheres are nearly independent at t = 0; that is, the joint probability f (0) factorizes to the maximum extent possible: f (0, x 1, v 1,..., x, v ) = Z 1 ( ) 1 xi x j >ε f 0 (x i, v i ) 1 i<j Here Z is a normalization factor ( partition function ), and f 0 (x, v) is nice (say, Schwartz class). i=1 Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 8 / 30

18 Lanford s theorem Assume that the spheres are nearly independent at t = 0; that is, the joint probability f (0) factorizes to the maximum extent possible: f (0, x 1, v 1,..., x, v ) = Z 1 ( ) 1 xi x j >ε f 0 (x i, v i ) 1 i<j Here Z is a normalization factor ( partition function ), and f 0 (x, v) is nice (say, Schwartz class). Define f (t) as the pushforward of f (0) through the dynamics, and also introduce the marginals f (s) (t) for 1 s : f (s) (t, Z s) = f (t, Z )dz s+1... dz R d( s) R d( s) i=1 Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 8 / 30

19 Lanford s theorem Recall f (0, x 1, v 1,..., x, v ) = Z 1 1 i<j 1 xi x j >ε Then we have (see (GSRT2014) for a quantitative version): ( ) f 0 (x i, v i ) Theorem (Lanford 1975; Gallagher, Saint-Raymond & Texier 2014) Consider the Boltzmann-Grad limit of hard spheres, ε d 1 = l 1, and let the particles be distributed initially like the function f (0) defined above, where f 0 is a function in the Schwartz class. Then for 0 t T (l, f 0 ), and any s, we have f (s) s (t) ft as pointwise a.e.; here f t (x, v) solves Boltzmann s kinetic equation, ( t + v x ) f t (x, v) = l 1 Q(f t, f t )(x, v) i=1 Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 9 / 30

20 Lanford s theorem Further Remarks Lanford s proof only holds on a short time interval T (l, f 0 ) = O(l). o hydrodynamic limit! Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 10 / 30

21 Lanford s theorem Further Remarks Lanford s proof only holds on a short time interval T (l, f 0 ) = O(l). o hydrodynamic limit! The time has been improved for a tagged particle in an equilibrium background (Lebowitz & Spohn 1982). This result has been further improved by Bodineau, Gallagher & Saint-Raymond (2015), in order to derive the Brownian motion of the tagged particle. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 10 / 30

22 Lanford s theorem Further Remarks Lanford s proof only holds on a short time interval T (l, f 0 ) = O(l). o hydrodynamic limit! The time has been improved for a tagged particle in an equilibrium background (Lebowitz & Spohn 1982). This result has been further improved by Bodineau, Gallagher & Saint-Raymond (2015), in order to derive the Brownian motion of the tagged particle. The total volume occupied by particles is negligible, so Lanford s theorem only applies to perfect gases (hence, no phase transitions). Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 10 / 30

23 Lanford s theorem Further Remarks Lanford s proof only holds on a short time interval T (l, f 0 ) = O(l). o hydrodynamic limit! The time has been improved for a tagged particle in an equilibrium background (Lebowitz & Spohn 1982). This result has been further improved by Bodineau, Gallagher & Saint-Raymond (2015), in order to derive the Brownian motion of the tagged particle. The total volume occupied by particles is negligible, so Lanford s theorem only applies to perfect gases (hence, no phase transitions). Lanford s proof is extremely sensitive to the specific form of f (0), and that ansatz is not (!!) propagated forwards in time. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 10 / 30

24 Lanford s theorem Further Remarks Lanford s proof only holds on a short time interval T (l, f 0 ) = O(l). o hydrodynamic limit! The time has been improved for a tagged particle in an equilibrium background (Lebowitz & Spohn 1982). This result has been further improved by Bodineau, Gallagher & Saint-Raymond (2015), in order to derive the Brownian motion of the tagged particle. The total volume occupied by particles is negligible, so Lanford s theorem only applies to perfect gases (hence, no phase transitions). Lanford s proof is extremely sensitive to the specific form of f (0), and that ansatz is not (!!) propagated forwards in time. Improving the time = Difficult open problem! Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 10 / 30

25 Lanford s theorem Problem: The convergence at time t > 0 in Lanford s theorem is almost everywhere, whereas at t = 0 we have uniform convergence. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 11 / 30

26 Lanford s theorem Problem: The convergence at time t > 0 in Lanford s theorem is almost everywhere, whereas at t = 0 we have uniform convergence. In fact it is possible to prove that uniform convergence fails at positive times. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 11 / 30

27 Lanford s theorem Problem: The convergence at time t > 0 in Lanford s theorem is almost everywhere, whereas at t = 0 we have uniform convergence. In fact it is possible to prove that uniform convergence fails at positive times. Hence we cannot iterate the theorem taking the evolved state as initial data. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 11 / 30

28 Lanford s theorem Problem: The convergence at time t > 0 in Lanford s theorem is almost everywhere, whereas at t = 0 we have uniform convergence. In fact it is possible to prove that uniform convergence fails at positive times. Hence we cannot iterate the theorem taking the evolved state as initial data. The inability to propagate the initial strong factorization means we cannot use the full structure of Boltzmann s equation (e.g., H-theorem, etc.) This may present a significant obstacle in any attempt to improve Lanford s theorem. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 11 / 30

29 Lanford s theorem Problem: The convergence at time t > 0 in Lanford s theorem is almost everywhere, whereas at t = 0 we have uniform convergence. In fact it is possible to prove that uniform convergence fails at positive times. Hence we cannot iterate the theorem taking the evolved state as initial data. The inability to propagate the initial strong factorization means we cannot use the full structure of Boltzmann s equation (e.g., H-theorem, etc.) This may present a significant obstacle in any attempt to improve Lanford s theorem. (Then again there may exist more probabilistic approaches which avoid using Boltzmann s equation explicitly - we don t know!) Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 11 / 30

30 Lanford s theorem A partial remedy is the notion of strong chaos, which first appears in F. King s PhD thesis (1975), and also in Appendix A of (vblls1980). Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 12 / 30

31 Lanford s theorem A partial remedy is the notion of strong chaos, which first appears in F. King s PhD thesis (1975), and also in Appendix A of (vblls1980). Define the sets } B s = {(X s, V s ) R 2ds 1 i < j s : x i x j (1) } B s = {(X s, V s ) R 2ds t 0, (X s V s t, V s ) B s (2) Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 12 / 30

32 Lanford s theorem A partial remedy is the notion of strong chaos, which first appears in F. King s PhD thesis (1975), and also in Appendix A of (vblls1980). Define the sets } B s = {(X s, V s ) R 2ds 1 i < j s : x i x j (1) } B s = {(X s, V s ) R 2ds t 0, (X s V s t, V s ) B s (2) Theorem (King (1975); van Beijeren et al (1980)) Consider the Boltzmann-Grad limit of hard spheres, ε d 1 = l 1, with initial distribution f (0). Assume that s the marginal f (s) (0) converges as uniformly on compact subsets of B s to f0 s, and f (s) (0) are nice enough. Then for 0 t T (l, f 0 ), the evolved marginals f (s) (t) also converge uniformly on compact subsets of B s, to a function f t, which solves Boltzmann s equation. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 12 / 30

33 Lanford s theorem The theorem on the previous slide lets us iterate as long as no singularities appear in the dynamics. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 13 / 30

34 Lanford s theorem The theorem on the previous slide lets us iterate as long as no singularities appear in the dynamics. It does not give us any information at points where x i x j ε, which is a significant limitation because the physical interaction is confined to the set x i x j = ε. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 13 / 30

35 Lanford s theorem The theorem on the previous slide lets us iterate as long as no singularities appear in the dynamics. It does not give us any information at points where x i x j ε, which is a significant limitation because the physical interaction is confined to the set x i x j = ε. In (BGSRS2016), quantitative estimates are employed to obtain a uniform error at distance scales x i x j ε log 1 ε, and their notion of chaos is propagated and the theorem may be iterated (absent singularities). Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 13 / 30

36 Lanford s theorem The theorem on the previous slide lets us iterate as long as no singularities appear in the dynamics. It does not give us any information at points where x i x j ε, which is a significant limitation because the physical interaction is confined to the set x i x j = ε. In (BGSRS2016), quantitative estimates are employed to obtain a uniform error at distance scales x i x j ε log 1 ε, and their notion of chaos is propagated and the theorem may be iterated (absent singularities). Our contribution gives uniform estimates at all scales x i x j > ε, after deletion of a small set of velocities (absent singularities). Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 13 / 30

37 Lanford s theorem Let ψs t be the flow of s hard spheres of diameter ε. Introduce the sets { } D s = Z s = (X s, V s ) R 2ds 1 i j < s, x i x j > ε (3) K s = { Z s = (X s, V s ) D s t > 0, ψ t s Z s = (X s V s t, V s ) } (4) { } Us η = Z s = (X s, V s ) D s inf v i v j > η (5) 1 i<j s Definition { Let F (0) = f (s) }1 s (0). The sequence {F (0)} is nonuniformly f 0 -chaotic if there exists κ (0, 1) such that, for all R > 0, lim were η(ε) = ε κ. ( f (s) ) s (0) f 0 1 Ks U s η(ε) 1 L s i=1 v i 2 2R 2 = 0 (6) Zs Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 14 / 30

38 Theorem (D. 2016) Consider the Boltzmann-Grad limit of hard spheres, ε d 1 = l 1, with initial distribution f (0). Assume that the initial marginals f (s) (0) are nice enough. Then if the sequence {F (0)} is nonuniformly f 0 -chaotic, then for 0 t T (l, f 0 ), the sequence {F (t)} is nonuniformly f t -chaotic, where f t solves Boltzmann s equation. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 15 / 30

39 Lanford s theorem Further Remarks Here nice enough always means bounded by Gaussian functions in v, uniformly in. This is an extremely strong bound; most experts expect these bounds to fail on large time intervals for large data. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 16 / 30

40 Lanford s theorem Further Remarks Here nice enough always means bounded by Gaussian functions in v, uniformly in. This is an extremely strong bound; most experts expect these bounds to fail on large time intervals for large data. It is possible to propagate (weak) chaos with weaker bounds, e.g. a recent result of Bodineau, Gallagher & Saint-Raymond uses an L 2 framework to derive linear hydrodynamic equations. Can we prove strong chaos in L 2? (Open question!) Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 16 / 30

41 Lanford s theorem Further Remarks Here nice enough always means bounded by Gaussian functions in v, uniformly in. This is an extremely strong bound; most experts expect these bounds to fail on large time intervals for large data. It is possible to propagate (weak) chaos with weaker bounds, e.g. a recent result of Bodineau, Gallagher & Saint-Raymond uses an L 2 framework to derive linear hydrodynamic equations. Can we prove strong chaos in L 2? (Open question!) Even better: strong chaos at the level of energy and entropy. This seems far out of reach, but if true, it would resolve Lanford s theorem once and for all. (Part of the problem is that uniqueness for Boltzmann s equation is unknown at this regularity!) Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 16 / 30

42 Dispersive inequalities for hard spheres Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 17 / 30

43 Dispersive inequalities for hard spheres Consider indentical hard spheres in R d, each with diameter ε > 0. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 18 / 30

44 Dispersive inequalities for hard spheres Consider indentical hard spheres in R d, each with diameter ε > 0. Suppose Z = (x 1, v 1,..., x, v ) is the initial configuration, and denote ψ t Z = Z (t) = (X (t), V (t)) the configuration at time t R. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 18 / 30

45 Dispersive inequalities for hard spheres Assume t > 0 and consider the following virial function: r Z (t) = ( xi (t) v i (t) v i (t) 2 t ) i=1 Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 19 / 30

46 Dispersive inequalities for hard spheres Assume t > 0 and consider the following virial function: r Z (t) = ( xi (t) v i (t) v i (t) 2 t ) i=1 Illner (1989) observed that this function is monotonic in time for almost every fixed data Z = (X, V ). Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 19 / 30

47 Dispersive inequalities for hard spheres Assume t > 0 and consider the following virial function: r Z (t) = ( xi (t) v i (t) v i (t) 2 t ) i=1 Illner (1989) observed that this function is monotonic in time for almost every fixed data Z = (X, V ). One can prove the following identity: r Z (T ) r Z (0) = ε k ω k (v ik (t k ) v j k (t k )) 0 where the sum indexes collisions along {Z (t)} 0 t T. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 19 / 30

48 Dispersive inequalities for hard spheres Illner s identity immediately implies that r Z (t) r Z (0) for a.e. Z and all t 0. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 20 / 30

49 Dispersive inequalities for hard spheres Illner s identity immediately implies that r Z (t) r Z (0) for a.e. Z and all t 0. This inequality is the only information we need to prove Lanford s theorem globally in time for small data near vacuum. (IP1989) Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 20 / 30

50 Dispersive inequalities for hard spheres Illner s identity immediately implies that r Z (t) r Z (0) for a.e. Z and all t 0. This inequality is the only information we need to prove Lanford s theorem globally in time for small data near vacuum. (IP1989) Indeed the inequality r Z (t) r Z (0) (t 0) allows us to compare the particle flow to free transport; this essentially forces the particles to undergo only finitely many collisions on average. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 20 / 30

51 Dispersive inequalities for hard spheres Recall f (0, x 1, v 1,..., x, v ) = Z 1 Then we have: 1 i<j 1 xi x j >ε ( ) f 0 (x i, v i ) Theorem (Illner & Pulvirenti 1986; Illner & Pulvirenti 1989) Consider the Boltzmann-Grad limit of hard spheres, ε d 1 = l 1, and let the particles be distributed initially like the function f (0) defined above, where f 0 is a function in the Schwartz class. Then if l is sufficiently large (depending on f 0 ) then for t 0, and any s, we have f (s) s (t) ft as pointwise a.e.; here f t (x, v) solves Boltzmann s kinetic equation, ( t + v x ) f t (x, v) = l 1 Q(f t, f t )(x, v) i=1 Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 21 / 30

52 Ideas behind the proof Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 22 / 30

53 Ideas behind the proof The proof follows the strategy of Lanford; in particular, it relies on the BBGKY hierarchy. The letters stand for Bogoliubov-Born-Green-Kirkwood-Yvon. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 23 / 30

54 Ideas behind the proof The proof follows the strategy of Lanford; in particular, it relies on the BBGKY hierarchy. The letters stand for Bogoliubov-Born-Green-Kirkwood-Yvon. The BBGKY hierarchy is the evolution equation satisfied by the marginals (t) of f (t). f (s) Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 23 / 30

55 Ideas behind the proof The proof follows the strategy of Lanford; in particular, it relies on the BBGKY hierarchy. The letters stand for Bogoliubov-Born-Green-Kirkwood-Yvon. The BBGKY hierarchy is the evolution equation satisfied by the marginals f (s) (t) of f (t). Lanford solved this equation on a small time interval to prove his theorem. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 23 / 30

56 Ideas behind the proof The proof follows the strategy of Lanford; in particular, it relies on the BBGKY hierarchy. The letters stand for Bogoliubov-Born-Green-Kirkwood-Yvon. The BBGKY hierarchy is the evolution equation satisfied by the marginals f (s) (t) of f (t). Lanford solved this equation on a small time interval to prove his theorem. We prove uniform bounds on the marginals f (s) (t) in L ; these bounds are global in time if the mean free path is large compared to f 0. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 23 / 30

57 Ideas behind the proof The proof follows the strategy of Lanford; in particular, it relies on the BBGKY hierarchy. The letters stand for Bogoliubov-Born-Green-Kirkwood-Yvon. The BBGKY hierarchy is the evolution equation satisfied by the marginals f (s) (t) of f (t). Lanford solved this equation on a small time interval to prove his theorem. We prove uniform bounds on the marginals f (s) (t) in L ; these bounds are global in time if the mean free path is large compared to f 0. Finally, following the classical strategy, we use the uniform bounds to control the probability of pathological collision trees. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 23 / 30

58 Ideas behind the proof { Let F (t) = f (s) }1 s (t) sphere BBGKY hierarchy: ( t + V s Xs ) f (s) C s+1 f (s+1) (t) = ( s)ε d 1 be a sequence of densities solving the hard (t) = l 1 C s+1 f (s+1) (t) s i=1 R d S d 1 ω (v s+1 v i ) f (s+1) (t,..., x i, v i,..., x i + εω, v s+1 )dωdv s+1 Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 24 / 30

59 Ideas behind the proof { Let F (t) = f (s) }1 s (t) sphere BBGKY hierarchy: ( t + V s Xs ) f (s) C s+1 f (s+1) (t) = ( s)ε d 1 be a sequence of densities solving the hard (t) = l 1 C s+1 f (s+1) (t) s i=1 R d S d 1 ω (v s+1 v i ) f (s+1) (t,..., x i, v i,..., x i + εω, v s+1 )dωdv s+1 Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 24 / 30

60 Ideas behind the proof Lanford has shown that if the following bounds hold at t = 0, for some β > 0 and µ R: sup sup sup e µs e 1 2 β s i=1 v i 2 (s) f (0, Z s) < (7) 1 s Z s D s then for some small T we have the following bounds for some β 1 < β and µ 1 < µ: sup sup sup sup e µ 1s e 1 2 β s 1 i=1 v i 2 (s) f (t, Z s) < (8) 0 t T 1 s Z s D s More generally we can assume that these bounds hold up to some T > 0 which may not be small. Such bounds have been justified under restrictive assumptions (e.g. large mean free path or tagged particle). Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 25 / 30

61 Ideas behind the proof Lanford s convergence proof is based on a series solution for the BBGKY hierarchy; this series counts all possible interactions, without any consideration for symmetries or cancellations. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 26 / 30

62 Ideas behind the proof Lanford s convergence proof is based on a series solution for the BBGKY hierarchy; this series counts all possible interactions, without any consideration for symmetries or cancellations. ( t + V s Xs ) f (s) (t) = l 1 C s+1 f (s+1) (t) Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 26 / 30

63 Ideas behind the proof Lanford s convergence proof is based on a series solution for the BBGKY hierarchy; this series counts all possible interactions, without any consideration for symmetries or cancellations. ( t + V s Xs ) f (s) (t) = l 1 C s+1 f (s+1) (t) f (s) t (t) = T s(t)f (s) (0) + l 1 T s (t t 1 )C s+1 f (s+1) (t 1 )dt 1 0 Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 26 / 30

64 Ideas behind the proof Lanford s convergence proof is based on a series solution for the BBGKY hierarchy; this series counts all possible interactions, without any consideration for symmetries or cancellations. ( t + V s Xs ) f (s) (t) = l 1 C s+1 f (s+1) (t) f (s) f (s) t (t) = T s(t)f (s) (0) + l 1 T s (t t 1 )C s+1 f (s+1) (t 1 )dt 1 0 s (t) = t l j j=0 0 tj 1 T s (t t 1 )C s+1... T s+j (t j )f (s+j) (0)dt j... dt 1 0 Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 26 / 30

65 Ideas behind the proof Lanford s convergence proof is based on a series solution for the BBGKY hierarchy; this series counts all possible interactions, without any consideration for symmetries or cancellations. ( t + V s Xs ) f (s) (t) = l 1 C s+1 f (s+1) (t) f (s) f (s) t (t) = T s(t)f (s) (0) + l 1 T s (t t 1 )C s+1 f (s+1) (t 1 )dt 1 0 s (t) = t l j j=0 0 tj 1 T s (t t 1 )C s+1... T s+j (t j )f (s+j) (0)dt j... dt 1 0 Problem: The series is only known to converge on a short time interval; this is why Lanford s proof breaks down after a short time. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 26 / 30

66 Ideas behind the proof Similar to (BGSR2015), we use the uniform estimate on the marginals, f (s) (t), to break the interval [0, T ] into manageable pieces. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 27 / 30

67 Ideas behind the proof Similar to (BGSR2015), we use the uniform estimate on the marginals, f (s) (t), to break the interval [0, T ] into manageable pieces. Let h > 0 be a small time interval. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 27 / 30

68 Ideas behind the proof Similar to (BGSR2015), we use the uniform estimate on the marginals, f (s) (t), to break the interval [0, T ] into manageable pieces. Let h > 0 be a small time interval. The BBGKY series expansion defines (informally) a branching process running backwards in time. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 27 / 30

69 Ideas behind the proof Similar to (BGSR2015), we use the uniform estimate on the marginals, f (s) (t), to break the interval [0, T ] into manageable pieces. Let h > 0 be a small time interval. The BBGKY series expansion defines (informally) a branching process running backwards in time. Good branching histories are those growing sub-exponentially at scale h: The branch points refer to fictitious collisions with fictitious particles! Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 27 / 30

70 Ideas behind the proof Following the scheme set forth by Lanford, we aim to show that typical branching histories converge (as ) to some limiting branching process governed by Boltzmann s equation. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 28 / 30

71 Ideas behind the proof Following the scheme set forth by Lanford, we aim to show that typical branching histories converge (as ) to some limiting branching process governed by Boltzmann s equation. The missing link is the estimation of pathological collision trees. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 28 / 30

72 Ideas behind the proof Following the scheme set forth by Lanford, we aim to show that typical branching histories converge (as ) to some limiting branching process governed by Boltzmann s equation. The missing link is the estimation of pathological collision trees. These are the trees which exhibit super-exponential growth of collisions, or contain recollisions (i.e. collisions which are not associated with a collision operator C s+k ). Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 28 / 30

73 Ideas behind the proof Following the scheme set forth by Lanford, we aim to show that typical branching histories converge (as ) to some limiting branching process governed by Boltzmann s equation. The missing link is the estimation of pathological collision trees. These are the trees which exhibit super-exponential growth of collisions, or contain recollisions (i.e. collisions which are not associated with a collision operator C s+k ). We estimate all the pathological trees using the uniform bounds, and show that the remaining trees are close to those generated by the Boltzmann process. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 28 / 30

74 Ideas behind the proof Following the scheme set forth by Lanford, we aim to show that typical branching histories converge (as ) to some limiting branching process governed by Boltzmann s equation. The missing link is the estimation of pathological collision trees. These are the trees which exhibit super-exponential growth of collisions, or contain recollisions (i.e. collisions which are not associated with a collision operator C s+k ). We estimate all the pathological trees using the uniform bounds, and show that the remaining trees are close to those generated by the Boltzmann process. This concludes the proof. Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 28 / 30

75 References T. Bodineau, I. Gallagher & L. Saint-Raymond. arxiv: (2015). (BGSR2015) T. Bodineau, I. Gallagher & L. Saint-Raymond. Invent. math. (2015). (BGSRS2016) T. Bodineau, I. Gallagher, L. Saint-Raymond, & S. Simonella. arxiv: (2016) R. Denlinger. arxiv: (2016). (GSRT2014) I. Gallagher, L. Saint-Raymond & B. Texier. Zurich Lec. Adv. Math (2014). R. Illner. Transport Theory and Stat. Phys. 18(1), (1989). (IP1989) R. Illner & M. Pulvirenti. Comm. Math. Phys. 121(1), (1989). F. King, Ph.D. thesis, O.E. Lanford. In: J. Moser (ed.), Dynamical Systems, Theory and Applications, Lec. otes in Physics, vol. 38, pp (1975) (vblls1980) H. van Beijeren, O. E. Lanford, J. L. Lebowitz, and H. Spohn. J. Stat. Phys. 22 no. 2, (1980). Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 29 / 30

76 Thank you! Ryan Denlinger (UT Austin) Propagation of chaos for hard spheres Caltech 30 / 30

Entropy and irreversibility in gas dynamics. Joint work with T. Bodineau, I. Gallagher and S. Simonella

Entropy and irreversibility in gas dynamics. Joint work with T. Bodineau, I. Gallagher and S. Simonella Entropy and irreversibility in gas dynamics Joint work with T. Bodineau, I. Gallagher and S. Simonella Kinetic description for a gas of hard spheres Hard sphere dynamics The system evolves under the combined

More information

VALIDITY OF THE BOLTZMANN EQUATION

VALIDITY OF THE BOLTZMANN EQUATION VALIDITY OF THE BOLTZMANN EQUATION BEYOND HARD SPHERES based on joint work with M. Pulvirenti and C. Saffirio Sergio Simonella Technische Universität München Sergio Simonella - TU München Academia Sinica

More information

From Newton s law to the linear Boltzmann equation without cut-off

From Newton s law to the linear Boltzmann equation without cut-off From Newton s law to the linear Boltzmann equation without cut-off Nathalie Ayi 1,2 1 Laboratoire J.A. Dieudonné, Université de Nice Sophia-Antipolis 2 Project COFFEE, INRIA Sophia Antipolis Méditerranée

More information

Dynamique d un gaz de sphères dures et équation de Boltzmann

Dynamique d un gaz de sphères dures et équation de Boltzmann Dynamique d un gaz de sphères dures et équation de Boltzmann Thierry Bodineau Travaux avec Isabelle Gallagher, Laure Saint-Raymond Outline. Linear Boltzmann equation & Brownian motion Linearized Boltzmann

More information

From Newton s law to the linear Boltzmann equation without cut-off

From Newton s law to the linear Boltzmann equation without cut-off From Newton s law to the linear Boltzmann equation without cut-off Nathalie Ayi (1) (1) Université Pierre et Marie Curie 27 Octobre 217 Nathalie AYI Séminaire LJLL 27 Octobre 217 1 / 35 Organisation of

More information

Diffusion d une particule marquée dans un gaz dilué de sphères dures

Diffusion d une particule marquée dans un gaz dilué de sphères dures Diffusion d une particule marquée dans un gaz dilué de sphères dures Thierry Bodineau Isabelle Gallagher, Laure Saint-Raymond Journées MAS 2014 Outline. Particle Diffusion Introduction Diluted Gas of hard

More information

Hilbert Sixth Problem

Hilbert Sixth Problem Academia Sinica, Taiwan Stanford University Newton Institute, September 28, 2010 : Mathematical Treatment of the Axioms of Physics. The investigations on the foundations of geometry suggest the problem:

More information

The existence of Burnett coefficients in the periodic Lorentz gas

The 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 information

Kinetic Transport Theory

Kinetic Transport Theory Lecture notes on Kinetic Transport Theory Christian Schmeiser 1 Preface Kinetic transport equations are mathematical descriptions of the dynamics of large particle ensembles in terms of a phase space (i.e.

More information

From molecular dynamics to kinetic theory and hydrodynamics

From molecular dynamics to kinetic theory and hydrodynamics From molecular dynamics to kinetic theory and hydrodynamics Thierry Bodineau, Isabelle Gallagher, and Laure Saint-Raymond Abstract. In these notes we present the main ingredients of the proof of the convergence

More information

Hydrodynamic Limits for the Boltzmann Equation

Hydrodynamic Limits for the Boltzmann Equation Hydrodynamic Limits for the Boltzmann Equation François Golse Université Paris 7 & Laboratoire J.-L. Lions golse@math.jussieu.fr Academia Sinica, Taipei, December 2004 LECTURE 2 FORMAL INCOMPRESSIBLE HYDRODYNAMIC

More information

CHEM-UA 652: Thermodynamics and Kinetics

CHEM-UA 652: Thermodynamics and Kinetics 1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 2 I. THE IDEAL GAS LAW In the last lecture, we discussed the Maxwell-Boltzmann velocity and speed distribution functions for an ideal gas. Remember

More information

From Hamiltonian particle systems to Kinetic equations

From Hamiltonian particle systems to Kinetic equations From Hamiltonian particle systems to Kinetic equations Università di Roma, La Sapienza WIAS, Berlin, February 2012 Plan of the lectures 1 Particle systems and BBKGY hierarchy (the paradigm of the Kinetic

More information

Isabelle Gallagher. Laure Saint-Raymond. Benjamin Texier FROM NEWTON TO BOLTZMANN: HARD SPHERES AND SHORT-RANGE POTENTIALS

Isabelle Gallagher. Laure Saint-Raymond. Benjamin Texier FROM NEWTON TO BOLTZMANN: HARD SPHERES AND SHORT-RANGE POTENTIALS Isabelle Gallagher Laure Saint-Raymond Benjamin Texier FROM EWTO TO BOLTZMA: HARD SPHERES AD SHORT-RAGE POTETIALS I. Gallagher Institut de Mathématiques de Jussieu UMR CRS 7586, Université Paris-Diderot

More information

Title of communication, titles not fitting in one line will break automatically

Title of communication, titles not fitting in one line will break automatically Title of communication titles not fitting in one line will break automatically First Author Second Author 2 Department University City Country 2 Other Institute City Country Abstract If you want to add

More information

The Boltzmann Equation and Its Applications

The Boltzmann Equation and Its Applications Carlo Cercignani The Boltzmann Equation and Its Applications With 42 Illustrations Springer-Verlag New York Berlin Heidelberg London Paris Tokyo CONTENTS PREFACE vii I. BASIC PRINCIPLES OF THE KINETIC

More information

The Hopf equation. The Hopf equation A toy model of fluid mechanics

The Hopf equation. The Hopf equation A toy model of fluid mechanics The Hopf equation A toy model of fluid mechanics 1. Main physical features Mathematical description of a continuous medium At the microscopic level, a fluid is a collection of interacting particles (Van

More information

From microscopic dynamics to kinetic equations

From microscopic dynamics to kinetic equations Facoltà di Scienze Matematiche, Fisiche e Naturali PhD Thesis in Mathematics, Sapienza Università di Roma XXV ciclo From microscopic dynamics to kinetic equations Chiara Sa rio Advisor: Prof. Mario Pulvirenti

More information

REGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS

REGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS 1 REGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS HENNING STRUCHTRUP ETH Zürich, Department of Materials, Polymer Physics, CH-8093 Zürich, Switzerland (on leave from University of Victoria,

More information

G : Statistical Mechanics Notes for Lecture 3 I. MICROCANONICAL ENSEMBLE: CONDITIONS FOR THERMAL EQUILIBRIUM Consider bringing two systems into

G : Statistical Mechanics Notes for Lecture 3 I. MICROCANONICAL ENSEMBLE: CONDITIONS FOR THERMAL EQUILIBRIUM Consider bringing two systems into G25.2651: Statistical Mechanics Notes for Lecture 3 I. MICROCANONICAL ENSEMBLE: CONDITIONS FOR THERMAL EQUILIBRIUM Consider bringing two systems into thermal contact. By thermal contact, we mean that the

More information

THE DERIVATION OF THE LINEAR BOLTZMANN EQUATION FROM A RAYLEIGH GAS PARTICLE MODEL

THE DERIVATION OF THE LINEAR BOLTZMANN EQUATION FROM A RAYLEIGH GAS PARTICLE MODEL THE DERIVATION OF THE LINEAR BOLTZMANN EQUATION FROM A RAYLEIGH GAS PARTICLE MODEL KARSTEN MATTHIES, GEORGE STONE, AND FLORIAN THEIL Abstract. A linear Boltzmann equation is derived in the Boltzmann-Grad

More information

Exponential tail behavior for solutions to the homogeneous Boltzmann equation

Exponential tail behavior for solutions to the homogeneous Boltzmann equation Exponential tail behavior for solutions to the homogeneous Boltzmann equation Maja Tasković The University of Texas at Austin Young Researchers Workshop: Kinetic theory with applications in physical sciences

More information

Derivation of Kinetic Equations

Derivation of Kinetic Equations CHAPTER 2 Derivation of Kinetic Equations As we said, the mathematical object that we consider in Kinetic Theory is the distribution function 0 apple f(t, x, v). We will now be a bit more precise about

More information

Goldstone modes in Lyapunov spectra of hard sphere systems

Goldstone modes in Lyapunov spectra of hard sphere systems PHYSICAL REVIEW E 70, 016207 (2004) Goldstone modes in Lyapunov spectra of hard sphere systems Astrid S. de Wijn* and Henk van Beijeren Institute for Theoretical Physics, Utrecht University, Leuvenlaan

More information

Physical Modeling of Multiphase flow. Boltzmann method

Physical Modeling of Multiphase flow. Boltzmann method with lattice Boltzmann method Exa Corp., Burlington, MA, USA Feburary, 2011 Scope Scope Re-examine the non-ideal gas model in [Shan & Chen, Phys. Rev. E, (1993)] from the perspective of kinetic theory

More information

Non-equilibrium phenomena and fluctuation relations

Non-equilibrium phenomena and fluctuation relations Non-equilibrium phenomena and fluctuation relations Lamberto Rondoni Politecnico di Torino Beijing 16 March 2012 http://www.rarenoise.lnl.infn.it/ Outline 1 Background: Local Thermodyamic Equilibrium 2

More information

G : Statistical Mechanics

G : Statistical Mechanics G25.2651: Statistical Mechanics Notes for Lecture 1 Defining statistical mechanics: Statistical Mechanics provies the connection between microscopic motion of individual atoms of matter and macroscopically

More information

1. Thermodynamics 1.1. A macroscopic view of matter

1. Thermodynamics 1.1. A macroscopic view of matter 1. Thermodynamics 1.1. A macroscopic view of matter Intensive: independent of the amount of substance, e.g. temperature,pressure. Extensive: depends on the amount of substance, e.g. internal energy, enthalpy.

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

Fluctuation theorems. Proseminar in theoretical physics Vincent Beaud ETH Zürich May 11th 2009

Fluctuation theorems. Proseminar in theoretical physics Vincent Beaud ETH Zürich May 11th 2009 Fluctuation theorems Proseminar in theoretical physics Vincent Beaud ETH Zürich May 11th 2009 Outline Introduction Equilibrium systems Theoretical background Non-equilibrium systems Fluctuations and small

More information

Particles approximation for Vlasov equation with singular interaction

Particles approximation for Vlasov equation with singular interaction Particles approximation for Vlasov equation with singular interaction M. Hauray, in collaboration with P.-E. Jabin. Université d Aix-Marseille GRAVASCO Worshop, IHP, November 2013 M. Hauray (UAM) Particles

More information

Hypocoercivity for kinetic equations with linear relaxation terms

Hypocoercivity for kinetic equations with linear relaxation terms Hypocoercivity for kinetic equations with linear relaxation terms Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (A JOINT

More information

Verschränkung versus Stosszahlansatz: The second law rests on low correlation levels in our cosmic neighborhood

Verschränkung versus Stosszahlansatz: The second law rests on low correlation levels in our cosmic neighborhood Verschränkung versus Stosszahlansatz: The second law rests on low correlation levels in our cosmic neighborhood Talk presented at the 2012 Anacapa Society Meeting Hamline University, Saint Paul, Minnesota

More information

Kinetic models of Maxwell type. A brief history Part I

Kinetic models of Maxwell type. A brief history Part I . A brief history Part I Department of Mathematics University of Pavia, Italy Porto Ercole, June 8-10 2008 Summer School METHODS AND MODELS OF KINETIC THEORY Outline 1 Introduction Wild result The central

More information

at time t, in dimension d. The index i varies in a countable set I. We call configuration the family, denoted generically by Φ: U (x i (t) x j (t))

at time t, in dimension d. The index i varies in a countable set I. We call configuration the family, denoted generically by Φ: U (x i (t) x j (t)) Notations In this chapter we investigate infinite systems of interacting particles subject to Newtonian dynamics Each particle is characterized by its position an velocity x i t, v i t R d R d at time

More information

Complex systems: Self-organization vs chaos assumption

Complex systems: Self-organization vs chaos assumption 1 Complex systems: Self-organization vs chaos assumption P. Degond Institut de Mathématiques de Toulouse CNRS and Université Paul Sabatier pierre.degond@math.univ-toulouse.fr (see http://sites.google.com/site/degond/)

More information

From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray

From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park

More information

Potential Descending Principle, Dynamic Law of Physical Motion and Statistical Theory of Heat

Potential Descending Principle, Dynamic Law of Physical Motion and Statistical Theory of Heat Potential Descending Principle, Dynamic Law of Physical Motion and Statistical Theory of Heat Tian Ma and Shouhong Wang Supported in part by NSF, ONR and Chinese NSF http://www.indiana.edu/ fluid Outline

More information

Entropy A measure of molecular disorder

Entropy A measure of molecular disorder Entropy A measure of molecular disorder Second Law uses Entropy, S, to identify spontaneous change. Restatement of Second Law: The entropy of the universe tends always towards a maximum (S universe > 0

More information

Continuum Mechanics Fundamentals

Continuum Mechanics Fundamentals Continuum Mechanics Fundamentals James R. Rice, notes for ES 220, 12 November 2009; corrections 9 December 2009 Conserved Quantities Let a conseved quantity have amount F per unit volume. Examples are

More information

Complex Analysis MATH 6300 Fall 2013 Homework 4

Complex Analysis MATH 6300 Fall 2013 Homework 4 Complex Analysis MATH 6300 Fall 2013 Homework 4 Due Wednesday, December 11 at 5 PM Note that to get full credit on any problem in this class, you must solve the problems in an efficient and elegant manner,

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

FUNDAMENTALS OF CHEMISTRY Vol. II - Irreversible Processes: Phenomenological and Statistical Approach - Carlo Cercignani

FUNDAMENTALS OF CHEMISTRY Vol. II - Irreversible Processes: Phenomenological and Statistical Approach - Carlo Cercignani IRREVERSIBLE PROCESSES: PHENOMENOLOGICAL AND STATISTICAL APPROACH Carlo Dipartimento di Matematica, Politecnico di Milano, Milano, Italy Keywords: Kinetic theory, thermodynamics, Boltzmann equation, Macroscopic

More information

Non-equilibrium mixtures of gases: modeling and computation

Non-equilibrium mixtures of gases: modeling and computation Non-equilibrium mixtures of gases: modeling and computation Lecture 1: and kinetic theory of gases Srboljub Simić University of Novi Sad, Serbia Aim and outline of the course Aim of the course To present

More information

MACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM. Contents AND BOLTZMANN ENTROPY. 1 Macroscopic Variables 3. 2 Local quantities and Hydrodynamics fields 4

MACROSCOPIC 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 information

5 Measure theory II. (or. lim. Prove the proposition. 5. For fixed F A and φ M define the restriction of φ on F by writing.

5 Measure theory II. (or. lim. Prove the proposition. 5. For fixed F A and φ M define the restriction of φ on F by writing. 5 Measure theory II 1. Charges (signed measures). Let (Ω, A) be a σ -algebra. A map φ: A R is called a charge, (or signed measure or σ -additive set function) if φ = φ(a j ) (5.1) A j for any disjoint

More information

Fluid Dynamics from Kinetic Equations

Fluid Dynamics from Kinetic Equations Fluid Dynamics from Kinetic Equations François Golse Université Paris 7 & IUF, Laboratoire J.-L. Lions golse@math.jussieu.fr & C. David Levermore University of Maryland, Dept. of Mathematics & IPST lvrmr@math.umd.edu

More information

The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method

The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method Alexis Vasseur, and Yi Wang Department of Mathematics, University of Texas

More information

YAN GUO, JUHI JANG, AND NING JIANG

YAN GUO, JUHI JANG, AND NING JIANG LOCAL HILBERT EXPANSION FOR THE BOLTZMANN EQUATION YAN GUO, JUHI JANG, AND NING JIANG Abstract. We revisit the classical ork of Caflisch [C] for compressible Euler limit of the Boltzmann equation. By using

More information

Lecture 8. The Second Law of Thermodynamics; Energy Exchange

Lecture 8. The Second Law of Thermodynamics; Energy Exchange Lecture 8 The Second Law of Thermodynamics; Energy Exchange The second law of thermodynamics Statistics of energy exchange General definition of temperature Why heat flows from hot to cold Reading for

More information

arxiv:cond-mat/ v2 [cond-mat.stat-mech] 25 Sep 2000

arxiv:cond-mat/ v2 [cond-mat.stat-mech] 25 Sep 2000 technical note, cond-mat/0009244 arxiv:cond-mat/0009244v2 [cond-mat.stat-mech] 25 Sep 2000 Jarzynski Relations for Quantum Systems and Some Applications Hal Tasaki 1 1 Introduction In a series of papers

More information

Une approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck

Une approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck Une approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (EN

More information

Ψ ν ), if ν is ω-normal, Ψω Ψ ν

Ψ ν ), if ν is ω-normal, Ψω Ψ ν Entropy production VOJKAN JAŠIĆ 1, CLAUDE-ALAIN PILLET 2 1 Department of Mathematics and Statistics McGill University 85 Sherbrooke Street West Montreal, QC, H3A 2K6, Canada jaksic@math.mcgill.ca 2 CPT-CNRS,

More information

7 The Navier-Stokes Equations

7 The Navier-Stokes Equations 18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and

More information

Relaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni

Relaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)

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

n v molecules will pass per unit time through the area from left to

n v molecules will pass per unit time through the area from left to 3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

NON-EQUILIBRIUM THERMODYNAMICS

NON-EQUILIBRIUM THERMODYNAMICS NON-EQUILIBRIUM THERMODYNAMICS S. R. DE GROOT Professor of Theoretical Physics University of Amsterdam, The Netherlands E MAZUR Professor of Theoretical Physics University of Leiden, The Netherlands DOVER

More information

Lecture 8. The Second Law of Thermodynamics; Energy Exchange

Lecture 8. The Second Law of Thermodynamics; Energy Exchange Lecture 8 The Second Law of Thermodynamics; Energy Exchange The second law of thermodynamics Statistics of energy exchange General definition of temperature Why heat flows from hot to cold Reading for

More information

6.730 Physics for Solid State Applications

6.730 Physics for Solid State Applications 6.730 Physics for Solid State Applications Lecture 25: Chemical Potential and Equilibrium Outline Microstates and Counting System and Reservoir Microstates Constants in Equilibrium Temperature & Chemical

More information

Particles approximation for Vlasov equation with singular interaction

Particles approximation for Vlasov equation with singular interaction Particles approximation for Vlasov equation with singular interaction M. Hauray, in collaboration with P.-E. Jabin. Université d Aix-Marseille Oberwolfach Worshop, December 2013 M. Hauray (UAM) Particles

More information

Spring_#8. Thermodynamics. Youngsuk Nam

Spring_#8. Thermodynamics. Youngsuk Nam Spring_#8 Thermodynamics Youngsuk Nam ysnam1@khu.ac.krac kr Ch.7: Entropy Apply the second law of thermodynamics to processes. Define a new property called entropy to quantify the secondlaw effects. Establish

More information

On the Boltzmann equation: global solutions in one spatial dimension

On the Boltzmann equation: global solutions in one spatial dimension On the Boltzmann equation: global solutions in one spatial dimension Department of Mathematics & Statistics Colloque de mathématiques de Montréal Centre de Recherches Mathématiques November 11, 2005 Collaborators

More information

Statistical Mechanics

Statistical Mechanics 42 My God, He Plays Dice! Statistical Mechanics Statistical Mechanics 43 Statistical Mechanics Statistical mechanics and thermodynamics are nineteenthcentury classical physics, but they contain the seeds

More information

DSMC simulations of thermal escape

DSMC simulations of thermal escape DSMC simulations of thermal escape Alexey N. Volkov, R.E. Johnson, O.J. Tucker, J.T. Erwin Materials Science & Engineering University of Virginia, USA Financial support is provided by NASA through Planetary

More information

Onsager theory: overview

Onsager theory: overview Onsager theory: overview Pearu Peterson December 18, 2006 1 Introduction Our aim is to study matter that consists of large number of molecules. A complete mechanical description of such a system is practically

More information

The general reason for approach to thermal equilibrium of macroscopic quantu

The general reason for approach to thermal equilibrium of macroscopic quantu The general reason for approach to thermal equilibrium of macroscopic quantum systems 10 May 2011 Joint work with S. Goldstein, J. L. Lebowitz, C. Mastrodonato, and N. Zanghì. Claim macroscopic quantum

More information

in Bounded Domains Ariane Trescases CMLA, ENS Cachan

in Bounded Domains Ariane Trescases CMLA, ENS Cachan CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline

More information

Computational Methods of Physical Problems with Mixed Scales

Computational Methods of Physical Problems with Mixed Scales Computational Methods of Physical Problems with Mixed Scales Shi Jin University of Wisconsin-Madison IPAM Tutorial, March 10-12, 2008 Scales and Physical Laws Figure from E & Engquist, AMS Notice Connections

More information

Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I

Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I J. A. Carrillo collaborators: T. Goudon (Lille), P. Lafitte (Lille) and F. Vecil (UAB) (CPDE 2005), (JCP, 2008), (JSC, 2008) ICREA

More information

Entropy generation and transport

Entropy 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 information

Analisis para la ecuacion de Boltzmann Soluciones y Approximaciones

Analisis para la ecuacion de Boltzmann Soluciones y Approximaciones Analisis para la ecuacion de Boltzmann Soluciones y Approximaciones Irene M. Gamba Department of Mathematics and ICES The University of Texas at Austin Mar del PLata, Julio 2012 Collaborators: R. Alonso,

More information

Lecture II: The Lorentz gas

Lecture II: The Lorentz gas Lecture II: The Lorentz gas Applications of measure rigidity: from number theory to statistical mechanics Simons Lectures, Stony Brook October 203 Jens Marklof University of Bristol http://www.maths.bristol.ac.uk

More information

Simplified Hyperbolic Moment Equations

Simplified Hyperbolic Moment Equations Simplified Hyperbolic Moment Equations Julian Koellermeier and Manuel Torrilhon Abstract Hyperbolicity is a necessary property of model equations for the solution of the BGK equation to achieve stable

More information

Mathematical Structures of Statistical Mechanics: from equilibrium to nonequilibrium and beyond Hao Ge

Mathematical Structures of Statistical Mechanics: from equilibrium to nonequilibrium and beyond Hao Ge Mathematical Structures of Statistical Mechanics: from equilibrium to nonequilibrium and beyond Hao Ge Beijing International Center for Mathematical Research and Biodynamic Optical Imaging Center Peking

More information

fiziks Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES

fiziks Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES Content-Thermodynamics & Statistical Mechanics 1. Kinetic theory of gases..(1-13) 1.1 Basic assumption of kinetic theory 1.1.1 Pressure exerted by a gas 1.2 Gas Law for Ideal gases: 1.2.1 Boyle s Law 1.2.2

More information

1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer 11(2) (1989),

1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer 11(2) (1989), Real Analysis 2, Math 651, Spring 2005 April 26, 2005 1 Real Analysis 2, Math 651, Spring 2005 Krzysztof Chris Ciesielski 1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer

More information

J10M.1 - Rod on a Rail (M93M.2)

J10M.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 information

Direct Modeling for Computational Fluid Dynamics

Direct Modeling for Computational Fluid Dynamics Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.

More information

Random Walks A&T and F&S 3.1.2

Random Walks A&T and F&S 3.1.2 Random Walks A&T 110-123 and F&S 3.1.2 As we explained last time, it is very difficult to sample directly a general probability distribution. - If we sample from another distribution, the overlap will

More information

Chapter 19. Chemical Thermodynamics

Chapter 19. Chemical Thermodynamics Chapter 19. Chemical Thermodynamics 19.1 Spontaneous Processes Chemical thermodynamics is concerned with energy relationships in chemical reactions. We consider enthalpy and we also consider entropy in

More information

Information in Biology

Information in Biology Lecture 3: Information in Biology Tsvi Tlusty, tsvi@unist.ac.kr Living information is carried by molecular channels Living systems I. Self-replicating information processors Environment II. III. Evolve

More information

Periodic Sinks and Observable Chaos

Periodic Sinks and Observable Chaos Periodic Sinks and Observable Chaos Systems of Study: Let M = S 1 R. T a,b,l : M M is a three-parameter family of maps defined by where θ S 1, r R. θ 1 = a+θ +Lsin2πθ +r r 1 = br +blsin2πθ Outline of Contents:

More information

Continuum Probability and Sets of Measure Zero

Continuum Probability and Sets of Measure Zero Chapter 3 Continuum Probability and Sets of Measure Zero In this chapter, we provide a motivation for using measure theory as a foundation for probability. It uses the example of random coin tossing to

More information

18.13 Review & Summary

18.13 Review & Summary 5/2/10 10:04 PM Print this page 18.13 Review & Summary Temperature; Thermometers Temperature is an SI base quantity related to our sense of hot and cold. It is measured with a thermometer, which contains

More information

Modeling and Analysis of Dynamic Systems

Modeling and Analysis of Dynamic Systems Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 34 Outline 1 Lecture 7: Recall on Thermodynamics

More information

The physics of information: from Maxwell s demon to Landauer. Eric Lutz University of Erlangen-Nürnberg

The physics of information: from Maxwell s demon to Landauer. Eric Lutz University of Erlangen-Nürnberg The physics of information: from Maxwell s demon to Landauer Eric Lutz University of Erlangen-Nürnberg Outline 1 Information and physics Information gain: Maxwell and Szilard Information erasure: Landauer

More information

84 My God, He Plays Dice! Chapter 12. Irreversibility. This chapter on the web informationphilosopher.com/problems/reversibility

84 My God, He Plays Dice! Chapter 12. Irreversibility. This chapter on the web informationphilosopher.com/problems/reversibility 84 My God, He Plays Dice! This chapter on the web informationphilosopher.com/problems/reversibility Microscopic In the 1870 s, Ludwig Boltzmann developed his transport equation and his dynamical H-theorem

More information

From billiards to thermodynamics. Renato Feres Washington University Department of Mathematics

From billiards to thermodynamics. Renato Feres Washington University Department of Mathematics From billiards to thermodynamics Renato Feres Washington University Department of Mathematics 1 / 21 Billiards Engraving from Charles Cotton s 1674 book, The Complete Gamester 2 / 21 Math billiards - different

More information

Particle-Simulation Methods for Fluid Dynamics

Particle-Simulation Methods for Fluid Dynamics Particle-Simulation Methods for Fluid Dynamics X. Y. Hu and Marco Ellero E-mail: Xiangyu.Hu and Marco.Ellero at mw.tum.de, WS 2012/2013: Lectures for Mechanical Engineering Institute of Aerodynamics Technical

More information

Piecewise Smooth Solutions to the Burgers-Hilbert Equation

Piecewise Smooth Solutions to the Burgers-Hilbert Equation Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang

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

Mean Field Limit for Stochastic Particle Systems

Mean Field Limit for Stochastic Particle Systems Mean Field Limit for Stochastic Particle Systems Pierre-Emmanuel Jabin and Zhenfu Wang Abstract We review some classical and more recent results for the derivation of mean field equations from systems

More information

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods Using Quadrature-Based Projection Methods Julian Köllermeier, Manuel Torrilhon 29th International Symposium on Rarefied Gas Dynamics, Xi an July 14th, 2014 Julian Köllermeier, Manuel Torrilhon 1 / 12 Outline

More information

Master 2 Macro I. Lecture notes #12 : Solving Dynamic Rational Expectations Models

Master 2 Macro I. Lecture notes #12 : Solving Dynamic Rational Expectations Models 2012-2013 Master 2 Macro I Lecture notes #12 : Solving Dynamic Rational Expectations Models Franck Portier (based on Gilles Saint-Paul lecture notes) franck.portier@tse-fr.eu Toulouse School of Economics

More information

Introduction. Statistical physics: microscopic foundation of thermodynamics degrees of freedom 2 3 state variables!

Introduction. Statistical physics: microscopic foundation of thermodynamics degrees of freedom 2 3 state variables! Introduction Thermodynamics: phenomenological description of equilibrium bulk properties of matter in terms of only a few state variables and thermodynamical laws. Statistical physics: microscopic foundation

More information

Chapter 4: Going from microcanonical to canonical ensemble, from energy to temperature.

Chapter 4: Going from microcanonical to canonical ensemble, from energy to temperature. Chapter 4: Going from microcanonical to canonical ensemble, from energy to temperature. All calculations in statistical mechanics can be done in the microcanonical ensemble, where all copies of the system

More information

Unique equilibrium states for geodesic flows in nonpositive curvature

Unique equilibrium states for geodesic flows in nonpositive curvature Unique equilibrium states for geodesic flows in nonpositive curvature Todd Fisher Department of Mathematics Brigham Young University Fractal Geometry, Hyperbolic Dynamics and Thermodynamical Formalism

More information

On the Interior Boundary-Value Problem for the Stationary Povzner Equation with Hard and Soft Interactions

On the Interior Boundary-Value Problem for the Stationary Povzner Equation with Hard and Soft Interactions On the Interior Boundary-Value Problem for the Stationary Povzner Equation with Hard and Soft Interactions Vladislav A. Panferov Department of Mathematics, Chalmers University of Technology and Göteborg

More information