Exponential moments for the homogeneous Kac equation

Similar documents
Exponential tail behavior for solutions to the homogeneous Boltzmann equation

arxiv: v2 [math.ap] 3 Mar 2017

Entropic structure of the Landau equation. Coulomb interaction

Classical inequalities for the Boltzmann collision operator.

On the Boltzmann equation: global solutions in one spatial dimension

Analisis para la ecuacion de Boltzmann Soluciones y Approximaciones

Analisis para la ecuacion de Boltzmann Soluciones y Approximaciones

Maja Tasković. Program Associate Mathematical Sciences Research Institute, six weeks in Fall 2015

A conservative spectral method for the Boltzmann equation with anisotropic scattering and the grazing collisions limit

Copyright by Ricardo Jose Alonso 2008

Maja Tasković. address:

Global Weak Solution to the Boltzmann-Enskog equation

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

The Boltzmann Equation and Its Applications

Maja Tasković. Address: Department of Mathematics University of Pennsylvania David Rittenhouse Lab 209 South 33rd Street Philadelphia PA 19104

Kinetic models of Maxwell type. A brief history Part I

ALMOST EXPONENTIAL DECAY NEAR MAXWELLIAN

A model of alignment interaction for oriented particles with phase transition

Kinetic theory of gases

arxiv: v2 [math.ap] 13 Oct 2010

SOME APPLICATIONS OF THE METHOD OF MOMENTS FOR THE HOMOGENEOUS BOLTZMANN AND KAC EQUATIONS

Modelling and numerical methods for the diffusion of impurities in a gas

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

Alignment processes on the sphere

Existence of global solutions to the Cauchy problem for the inelastic Boltzmann equation with near-vacuum data

Hypocoercivity for kinetic equations with linear relaxation terms

Slowing down the neutrons

Semigroup factorization and relaxation rates of kinetic equations

Fluid Dynamics from Kinetic Equations

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

Different types of phase transitions for a simple model of alignment of oriented particles

Hypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th

A model of alignment interaction for oriented particles with phase transition

in Bounded Domains Ariane Trescases CMLA, ENS Cachan

A model of alignment interaction for oriented particles with phase transition

Entropy-dissipation methods I: Fokker-Planck equations

Exponential methods for kinetic equations

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

Different types of phase transitions for a simple model of alignment of oriented particles

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

A Selective Modern History of the Boltzmann and Related Eq

Non-equilibrium mixtures of gases: modeling and computation

YAN GUO, JUHI JANG, AND NING JIANG

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

Fourier Law and Non-Isothermal Boundary in the Boltzmann Theory

Hölder regularity estimation by Hart Smith and Curvelet transforms

Tanaka Theorem for Inelastic Maxwell Models

Spectral - Lagrangian methods for Collisional Models of Non - Equilibrium Statistical States

A quick introduction to regularity structures

Monte Carlo method with negative particles

The Equipartition Theorem

Anomalous transport of particles in Plasma physics

1 Lyapunov theory of stability

BOLTZMANN KINETIC THEORY FOR INELASTIC MAXWELL MIXTURES

Moment Inequalities and High-Energy Tails for Boltzmann Equations with Inelastic Interactions

ON A NEW CLASS OF WEAK SOLUTIONS TO THE SPATIALLY HOMOGENEOUS BOLTZMANN AND LANDAU EQUATIONS

ABOUT THE REGULARIZING PROPERTIES OF THE NON CUT{OFF KAC EQUATION. Laurent Desvillettes. 45, Rue d'ulm Paris Cedex 05. April 19, 2002.

Quiz 3 for Physics 176: Answers. Professor Greenside

Transport coefficients in plasmas spanning weak to strong correlation

Monte Carlo methods for kinetic equations

hal , version 1-22 Nov 2009

Review of last lecture I

VALIDITY OF THE BOLTZMANN EQUATION

The Cauchy problem and BEC stability for the quantum Boltzmann-Condensation system for bosons at very low temperature

Uncertainty Quantification and hypocoercivity based sensitivity analysis for multiscale kinetic equations with random inputs.

Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations

Lecture Note 1: Background

On a New Diagram Notation for the Derivation of Hyperbolic Moment Models

Divergence of the gradient expansion and the applicability of fluid dynamics Gabriel S. Denicol (IF-UFF)

Fluid Approximations from the Boltzmann Equation for Domains with Boundary

ON THE SELF-SIMILAR ASYMPTOTICS FOR GENERALIZED NON-LINEAR KINETIC MAXWELL MODELS

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution

Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma

Physics Dec The Maxwell Velocity Distribution

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution

Physics 607 Final Exam

Physics 6303 Lecture 11 September 24, LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation

Mathematics 104 Fall Term 2006 Solutions to Final Exam. sin(ln t) dt = e x sin(x) dx.

Chaotic motion. Phys 420/580 Lecture 10

Numerical methods for kinetic equations

Lattice Boltzmann Method for Moving Boundaries

THE VLASOV-POISSON-BOLTZMANN SYSTEM FOR SOFT POTENTIALS

Math 32B Discussion Session Session 3 Notes August 14, 2018

A comparative study of no-time-counter and majorant collision frequency numerical schemes in DSMC

Topics for the Qualifying Examination

The propagation of chaos for a rarefied gas of hard spheres

Classical Fourier Analysis

Uncertainty Quantification for multiscale kinetic equations with random inputs. Shi Jin. University of Wisconsin-Madison, USA

Physics 342 Lecture 22. The Hydrogen Atom. Lecture 22. Physics 342 Quantum Mechanics I

From Hamiltonian particle systems to Kinetic equations

Stability of Stochastic Differential Equations

Rate of Heating and Cooling

Convergence to self-similarity for the Boltzmann equation for strongly inelastic Maxwell molecules

On Weak Solutions to the Linear Boltzmann Equation with Inelastic Coulomb Collisions

Chaotic motion. Phys 750 Lecture 9

Physical Modeling of Multiphase flow. Boltzmann method

arxiv: v1 [math.ap] 6 Apr 2016

Qualifying Exam. Aug Part II. Please use blank paper for your work do not write on problems sheets!

Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions and Power Like Tails

Recent advances in kinetic theory for mixtures of polyatomic gases

Transcription:

Exponential moments for the homogeneous Kac equation Maja Tasković University of Pennsylvania Young Researchers Workshop: Stochastic and deterministic methods in kinetic theory, Duke University, November 29, 2016. In collaboration with: Milana Pavić-Čolić (University of Novi Sad, Serbia) Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 1 / 23

Outline 1 Introduction 2 Previous results 3 Our contribution Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 2 / 23

Introduction Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 3 / 23

Similarities and differences b/w Kac and Boltzmann equations Similarities between the Kac equation and the Boltzmann equations: Both model the evolution of rarefied gases via a distribution function f (t, x, v). Binary interactions between particles. Quadratic collision operator. Our assumption: spatial homogeneity, i.e. we assume the distribution function f (t, v) does not depend on x. Differences between the Kac equation and the Boltzmann equation: The dimension of spatial and velocity spaces is d 2 for the Boltzmann equation, while for the Kac equation d = 1. Conservation laws: collisions modeled by the Boltzmann equation conserve mass, momentum and energy, while those for the Kac equation conserve only mass and energy (but not the momentum). Pre-post collisional velocity laws are be different. Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 4 / 23

The Kac equation The Kac equation (1959, Kac) models a 1-dimensional gas in which collisions conserve the mass and the energy, but not the momentum. The spatially homogeneous Kac equation reads tf (t, v) = R π π ( f f f f ) bk ( θ ) dθ dv, with the standard abbreviations f := f (t, v ), f := f (t, v ), f := f (t, v ). Conservation of energy v 2 + v 2 = v 2 + v 2 implies that the pre and post-collisional velocities v, v and v, v can be related by introducing a parameter θ [ π, π] as: v = v cos θ v sin θ, v = v sin θ + v cos θ. By a change of variables, one can also use the following relation: v = v 2 + v 2 cos θ, v = v 2 + v 2 sin θ. Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 5 / 23

The Boltzmann equation The Boltzmann equation (late 1860s and 1870s, Maxwell and Boltzmann) models evolution of a d-dimensional gas (d 2) in which particles interact via binary collisions. The spatially homogeneous Boltzmann equation reads: ( tf = f f ) ( ) f f v v γ b v v B σ R d S d 1 v v dσ dv, Conservation of momentum and energy implies the following relation between pre and post-collisional velocities (v, v ) and (v, v ): v v = v + v 2 = v + v 2 v v σ + 2 v v σ, σ S d 1. 2 Notation: u = v v, u = v v, γ (0, 1]: hard potentials, γ = 0: Maxwell molecules, γ ( d, 0): soft potentials. v v u σ θ û u v v O Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 6 / 23

The angular cross-section Angular kernels b K ( θ ) and b B (cos θ) often have a non-integrable singularity at θ = 0. However, such singularity is often cut-off ( b K ( θ ) L 1 ([0, π]) or b B (û σ) L 1 (S d 1 )). Cutoff simplifies the analysis of the collision operator by enabling its splitting into the gain and loss terms Q(f, f ) = Q + (f, f ) Q (f, f ). It was believed that this removal of the non-integrability does not influence the equation significantly. However, it has been observed that singularity carries regularizing effect. This motivates further study of the non-cutoff regime: Non-cutoff (Kac): π π π π b K ( θ ) dθ =, b K ( θ ) sin β θ dθ <, β (0, 2] Non-cutoff (Boltzmann): π 0 π 0 b B (cos θ) sin d 2 θ dθ =, b b (cos θ) sin β θ sin d 2 θ dθ <, β (0, 2] Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 7 / 23

Summary: Kac equation and Boltzmann equation for Maxwell molecules Kac equation: tf = π R π ( f f f f ) bk ( θ ) dθ dv, π π π π b K ( θ ) dθ =, b K ( θ ) sin β θ dθ <, β (0, 2] Boltzmann equation for Maxwell molecules: ( tf = f f ) f f bb R d S d 1 ( v v ) σ v v dσ dv, π 0 π 0 b B (cos θ) sin d 2 θ dθ =, b b (cos θ) sin β θ sin d 2 θ dθ <, β (0, 2] Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 8 / 23

What do we study? We say that f has an exponential tail in the L 1 if, for some rate α > 0 and some order s > 0, the following norm is finite: f L 1 (t) := f (t, v) e α s v dv <. exp(α v s ) R d Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 9 / 23

What do we study? We say that f has an exponential tail in the L 1 if, for some rate α > 0 and some order s > 0, the following norm is finite: f L 1 (t) := f (t, v) e α s v dv <. exp(α v s ) R d We study whether these norms propagate in time: propagation of exponentially weighted L 1 norms means: s v dv < C 0, for some α 0, s > 0 R d f 0 (v) e α 0 α, C > 0, t 0 : R d f (t, v) e α s v dv < C. Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 9 / 23

What do we study? We say that f has an exponential tail in the L 1 if, for some rate α > 0 and some order s > 0, the following norm is finite: f L 1 (t) := f (t, v) e α s v dv <. exp(α v s ) R d We study whether these norms propagate in time: propagation of exponentially weighted L 1 norms means: s v dv < C 0, for some α 0, s > 0 R d f 0 (v) e α 0 α, C > 0, t 0 : R d f (t, v) e α s v dv < C. generation of exponentially weighted L 1 norms means f 0 (v) v q dv < C 0, for some q N R d s α, s, C > 0, t > 0 : f (t, v) e α v dv < C. R d Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 9 / 23

Why exponential moments? They provide information about tail behavior of solutions (exponential decay of f (t, v) for large velocities). First step towards pointwise exponential bounds L exp. In fact, in 1972 Arkeryd asked if the following is true for the Boltzmann equation: f 0 (v) c e α 0 v 2 f (t, v) c 1 e α 1 v 2. First results in this direction were in the L 1 setting by Bobylev in 1984 and 1997. Original question was addressed by Gamba, Panferov and Villani in 2009 in the cutoff case, by relating L exp norm with L 1 exp norms. Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 10 / 23

Previous results Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 11 / 23

Previous results 1984 Bobylev: Boltzmann equation for Maxwell molecules, s = 2 1993 Desvillettes: Kac equation, s = 1, s = 2 1997 Bobylev: Hard sphere, s = 2. This paper opened the doors for many extensions: 2004 Bobylev-Gamba-Panferov 2009 Gamba-Panferov-Villani 2006 Mouhot 2006 2013 Alonso-Canizo-amba-Mouhot 2012 Lu-Mouhot 2012 2016 T.-Alonso-Gamba-Pavlović Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 12 / 23

Previous results 1984 Bobylev: Boltzmann equation for Maxwell molecules, s = 2 1993 Desvillettes: Kac equation, s = 1, s = 2 1997 Bobylev: Hard sphere, s = 2. This paper opened the doors for many extensions: 2004 Bobylev-Gamba-Panferov 2009 Gamba-Panferov-Villani 2006 Mouhot 2006 2013 Alonso-Canizo-amba-Mouhot 2012 Lu-Mouhot 2012 2016 T.-Alonso-Gamba-Pavlović Cutoff case for γ > 0: generation propagation 0 Non-cutoff case for γ > 0: generation propagation γ 1 2 s = order of the exponential e α v s 0 γ 1 2 s = order of the exponential e α v s Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 12 / 23

Desvillettes 1993 For the Kac equation with constant angular kernel and for orders s = 1 and s = 2, Desvillettes 1993 proved propagation of exponential moments of integer orders s = 1 and s = 2. His elegant proof works directly with exponential moments. For example, for s = 2, let the exponential moment of order 2 be denoted by M(t, α) = f (t, v)e α v 2 dv. Multiply the Kac equation with the exponential weight and integrate in v: tm(t, α) = K (f, f )e α v 2 dv R π ( = f f e α v 2 e α v 2) dθ R R π 2π dvdv π ( = f f e α( v 2 + v 2 ) cos 2 θ e α v 2) dθ 2π dvdv = R R π π π R M 2 (t, α cos 2 θ) M(t, α)m(0, 0) dθ 2π Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 13 / 23

Bobylev 1997 The first groundbreaking work on exponential moments for hard potentials and cutoff, was done by Bobylev in 1997. The key idea is to Taylor expand exponential function to reformulate (formally) the question to showing summability of polynomial moments renormalized by Gamma functions: f (t, v) e α s v dv = R d R d f (t, v) v qs α q dv = q! q=0 q=0 m qs(t) α q Γ(q + 1). Because of the required summability, constants need to be sharper compared to results concerning polynomial moments. Estimates are developed by looking for ordinary differential inequalities for m qs, and then deriving: term-by-term estimates partial sum estimates Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 14 / 23

Classical technique Multiply the Boltzmann equation with v 2q : tf v 2q = Q(f, f ) v 2q. Integrate in velocity. The weak form of the collision operator yields m 2q(t)= Q(f, f ) v 2q dv R d = 1 ( ( f f v 2q + v 2q v 2q v 2q) ) u γ b(û σ) dσ dvdv 2 R 2d S d 1 Then, look for a bound of the right-hand-side in terms of polynomial moments. For the bounds to be good enough it is important that the highest-order moment comes with a negative sign. m 2q K 1 m 2q+γ + K 2 m 2q q+1 2 ( ) q ) + K 3 ε q (m2k+γ m 2(q k) + m 2k m 2(q k)+γ k k=1 Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 15 / 23

Our contribution Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 16 / 23

Our main result Theorem (Pavic-Colic, T. 2016) Suppose initial datum f 0 0 has finite mass, energy and entropy. Let f (t, v) be an associated weak solution to the Kac equation with the angular kernel satisfying: If π π π π b K ( θ ) dθ =, b K ( θ ) sin β θ dθ <, β (0, 2] s 4 2 + β, then for every α 0 > 0 there exists 0 < α α 0 and a constant C > 0 (depending only on the initial data and κ) so that if f 0 (v)e α 0 v s dv M 0 <, R then f (t, v)e α v s dv C, t 0. R Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 17 / 23

Interpretation 2 s 1 4 2+β 0.5 1 1.5 2 β Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 18 / 23

Strategy We use recent techniques of T.-Alonso-Gamba-Pavlovic 2016, which consists of studying generalized exponential moments, so-called Mittag-Leffler moments. They are L 1 norms weighted with Mittag-Leffler functions: E a(x) := q=0 x q Γ(aq + 1) Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 19 / 23

Strategy We use recent techniques of T.-Alonso-Gamba-Pavlovic 2016, which consists of studying generalized exponential moments, so-called Mittag-Leffler moments. They are L 1 norms weighted with Mittag-Leffler functions: E a(x) := q=0 x q Γ(aq + 1) They generalize exponentials and are known to asymptotically behave like exp: E a(x) e x1/a, for x 1. Because of this asymptotic behavior, finiteness of exponential moment is equivalent to the finiteness of the corresponding Mittag-Leffler moment. So, we work with: Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 19 / 23

Strategy We use recent techniques of T.-Alonso-Gamba-Pavlovic 2016, which consists of studying generalized exponential moments, so-called Mittag-Leffler moments. They are L 1 norms weighted with Mittag-Leffler functions: E a(x) := q=0 x q Γ(aq + 1) They generalize exponentials and are known to asymptotically behave like exp: E a(x) e x1/a, for x 1. Because of this asymptotic behavior, finiteness of exponential moment is equivalent to the finiteness of the corresponding Mittag-Leffler moment. So, we work with: Definition (Mittag-Leffler moment) The Mittag-Leffler moment of f of order s and rate α > 0 is introduced via: ML (α,s) = q=0 m 2q α 2q/s Γ( 2 s q + 1). Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 19 / 23

Why Mittag-Leffler moments? part 1 We develop the partial sum technique for the non-cutoff case. The story begins as before m 2q(t)= K (f, f ) v 2q dv R = 1 π ( f f v 2q + v 2q v 2q v 2q) b K ( θ ) dθ dv dv. 2 R R π First challenge is angular singularity. Need to exploit certain cancellation properties that become visible after an application of the Taylor expansion to the test functions. This leads to an ODI for polynomial moments with another challenge: compared to the Grad s cutoff, there are two extra powers of q in the last term of the inequality: m 2q K 1 m 2q + K 2 m 2q 2 + K 3 ε β,q q (q 1) q+1 2 k=1 ( ) q 2 m 2k m 2(q k) k 1 Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 20 / 23

Why Mittag-Leffler moments? part 2 Sequence ε q depends on the angular kernel b(cos θ). It has the following property: β [0, 2] q 1 β 2 εβ,q 0, as q, β = 0 (Grad s cutoff) q ε 0,q 0, as q, β = 2 (Full non-cutoff) ε 2,q 0, as q. Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 21 / 23

Why Mittag-Leffler moments? part 2 Sequence ε q depends on the angular kernel b(cos θ). It has the following property: β [0, 2] q 1 β 2 εβ,q 0, as q, β = 0 (Grad s cutoff) q ε 0,q 0, as q, β = 2 (Full non-cutoff) ε 2,q 0, as q. Need to reduce the quadratic power of q(q 1) in m 2q... + ε β,q q (q 1) q+1 2 k=1 ( ) q 2 m 2k m 2(q k) k 1 Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 21 / 23

Why Mittag-Leffler moments? part 2 Sequence ε q depends on the angular kernel b(cos θ). It has the following property: β [0, 2] q 1 β 2 εβ,q 0, as q, β = 0 (Grad s cutoff) q ε 0,q 0, as q, β = 2 (Full non-cutoff) ε 2,q 0, as q. Need to reduce the quadratic power of q(q 1) in m 2q... + ε β,q q (q 1) q+1 2 k=1 ( ) q 2 m 2k m 2(q k) k 1 Going towards partial sums, divide all moments by appropriate Gamma functions: m 2q εβ,q q (q 1)... + Γ(aq + 1) Γ(aq + 1) q+1 2 k=1 ( q 2 k 1 ) m 2k Γ(ak + 1) m 2(q k) Γ(a(q k) + 1) Γ(aq + 2)B(ak + 1, a(q k) + 1) Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 21 / 23

Why Mittag-Leffler moments? part 3 Use q+1 2 k=1 ( q 2 k 1 ) B(ak + 1, a(q k) + 1) C a 1 q a+1 Last inequality holds only if a > 1, and hence the need for Mittag-Leffler moments. This yields ε β,q q 2 a in an ODE for partial sums. For this sequence to converge to zero as q, we need: 2 a = 2 2 s 1 β 2 Hence the order of the exponential moment depends on the signularity rate of the angular kernel as: s 4 2 + β. Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 22 / 23

Possible further questions 1 L exp bounds for the Kac equation 2 Exponential moments for soft potentials γ < 0? 3 Lower bounds? 4 Relation with the rate of convergence to the equilibrium? Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 23 / 23

Possible further questions 1 L exp bounds for the Kac equation 2 Exponential moments for soft potentials γ < 0? 3 Lower bounds? 4 Relation with the rate of convergence to the equilibrium? Thank you!!! Maja Tasković (UPenn) Exponential moments for the homogeneous Kac equation 23 / 23