Research Paper Notes on the Boltzmann Equation and its Hydrodynamic Limit

Page last edited 1,070 days ago
From LaTeX Notes Series
Jump to navigation Jump to search

Research Paper Notes on the Boltzmann Equation and its Hydrodynamic Limit

本文档除了包括推导,疑惑外,做读书重点的记录 <math></math>

参考文献

The Chapman-Enskog method and method of moments

  • On the Kinetic Theory of Rarefied Gases, Commun. Pure App. Math.2 (1949) 331-407, by Harold Grad
  • Transient Relativistic Thermodynamics and Kinetic Theory, Ann. Phys. 118 (1979) 341-372, W. Isreal and J. M. Stewart
  • Boltzmann Equation and Moment Equations in Curvilinear Coordinates, Fluid Dyn 2 (1967) 107–109, by E. M. Shakhov
  • Derivation of transient relativistic fluid dynamics from the Boltzmann equation, arXiv:1202.4551, by Gabriel S. Denicol, H. Niemi, E. Molnar, D.H. Rischke

Exact solution of the Boltzmann equation

  • A new exact solution of the relativistic Boltzmann equation and its hydrodynamic limit, arXiv:1408.5646, by Gabriel S. Denicol, Ulrich Heinz, Mauricio Martinez, Jorge Noronha, and Michael Strickland
  • Studying the validity of relativistic hydrodynamics with a new exact solution of the Boltzmann equation, arXiv:1408.7048, by Gabriel Denicol, Ulrich Heinz, Mauricio Martinez, Jorge Noronha, Michael Strickland

Anisotropic hydrodynamics, hydrodynamic attractor

  • The anisotropic non-equilibrium hydrodynamic attractor, arXiv:1709.06644, by Michael Strickland, Jorge Noronha, and Gabriel S. Denicol
  • Analytical attractor and the divergence of the slow-roll expansion in relativistic hydrodynamics, arXiv:1711.01657, by Gabriel S. Denicol and Jorge Noronha


On the Kinetic Theory of Rarefied Gases, Commun. Pure App. Math.2 (1949) 331-407, by Harold Grad

这是第一篇将玻尔兹曼方程按张量展开(14动量基展开),从而得到对应的流体力学方程的工作.

本文档除了包括推导,疑惑外,做读书重点的记录 <math></math>

Transient Relativistic Thermodynamics and Kinetic Theory, Ann. Phys. 118 (1979) 341-372, W. Isreal and J. M. Stewart

本文即著名的连接粘滞流体的Müller-Israel-Stewart模型. 基于一阶粘滞流体理论的因果性问题,这个工作有两个出发点. 第一是非平衡态下推广的热力学,它是基于熵流应该包含两阶项以引入延迟效应避免因果困难的现象学论述. 第二是利用moments展开由输运模型对流体力学方程的推导,它是基于这个展开可以导致双曲线型的偏微分方程避免因果性困难的数学角度论述. Israel正是广义相对论中曲面Israel条件的作者,本文不易读懂.

(1.1-1.4)

这里(1.1)是平衡态广延体系结果的推广,(1.2)同样是准静态过程热力学第一定律的推广,但被推广到偏离平衡态的情况. 如果只考虑准静态过程,则(1.3)是(1.1)的变分与(1.2)的差.在(1.2)被推广到非平衡态情况后,(1.3)也同样适用于平衡态到非平衡态的偏离. 考察从平衡态(1.1)出发到非平衡态的偏离(1.2),两者之和即(1.4). 又参见(2.13-15)和(2.17-20)附近的讨论

(1.5)

将(1.4)求散度,考虑守恒方程(2.1),以及(1.3)即得(1.5). 又参见(2.27)附近的讨论. 注意到这里的一个细节.以粒子数为例,在热力学关系中粒子数是可变<math>dn</math>,而粒子数守恒是指考虑粒子到在热力学过程中粒子没有湮灭,粒子数流守恒<math>\left(N^\mu\right)_{,\mu}=0</math>.更具体的,如通过选取参照系定义流体速度<math>u^\mu</math>的前提下,我们可以有<math>N^\mu =nu^\mu</math>.

本文档除了包括推导,疑惑外,做读书重点的记录 <math></math>

Boltzmann Equation and Moment Equations in Curvilinear Coordinates, Fluid Dyn 2 (1967) 107–109, by E. M. Shakhov

(1.9)

这是玻尔兹曼方程在曲线坐标系中的一般形式,在很多书籍和文献中都有直接应用.本地下载链接

这个推导中最重要的一点是,当以坐标变换的方式把动量分量也表达为曲线坐标系的形式,玻尔兹曼方程会牵涉到对曲线坐标系中动量分量的偏导,尽管最初的玻尔兹曼方程(1.1)不含外力,对动量的偏导无关.

这是因为在曲线坐标系中的动量分量不仅仅是在笛卡尔坐标系中动量分量的函数,还是曲线坐标系中坐标分量的函数.即由(1.4),

<math>q_i=q_i(x,y,z)</math>
<math>\xi_i=\xi_i(x,y,z,\xi_x,\xi_y,\xi_z)</math>

反之亦然,即

<math>\xi_x=\xi_x(q_1,q_2,q_3,\xi_1,\xi_2,\xi_3)</math>

这样原来的对坐标的偏导不再具有之前的形式.具体的

<math>\left.\frac{\partial f}{\partial q_i}\right|_{q_{j\ne i},\xi_1,\xi_2,\xi_3}\ne \left.\frac{\partial f}{\partial q_i}\right|_{q_{j\ne i},\xi_x,\xi_y,\xi_z}</math>.

所以方程右边的新增的项就是为了补偿这个区别,而把<math>\xi_x</math>等三个变量看做中间变量来得到需要的偏导

<math>\left.\frac{\partial f}{\partial q_i}\right|_{q_{j\ne i},\xi_1,\xi_2,\xi_3}=\left.\frac{\partial f}{\partial q_i}\right|_{q_{j\ne i},\xi_x,\xi_y,\xi_z}+\left.\frac{\partial f}{\partial \xi_x}\right|_{q_1,q_2,q_3,\xi_y,\xi_z}\left.\frac{\partial \xi_x}{\partial q_i}\right|_{q_{j\ne i},\xi_1,\xi_2,\xi_3}+\cdots</math>

这个结果实际上并不是直观的,很容易搞错.注意到经典力学中的广义坐标和动量的变换牵涉到完全类似的问题.

本文档除了包括推导,疑惑外,做读书重点的记录 <math></math>

Derivation of transient relativistic fluid dynamics from the Boltzmann equation, arXiv:1202.4551, by Gabriel S. Denicol, H. Niemi, E. Molnar, D.H. Rischke

本文为Gabriel经典之作,必须学习,理解,融通,超越.

(3)波色子对应<math>(1+f_k)</math>因子的来源.

甚至于文中引文中教科书都没有给出证明.因为对费米子<math>(1-f_k)</math>因子直接来源于泡利不相容原理.对玻色子同样的推演则无法实现.在网上搜索苦久.终于找到了一个推导.其实考虑在体系达到动态平衡的情况下,分布函数<math>f_k</math>对时间全微分为零,换言之,方程右边碰撞项的贡献等于零,被积函数为零意味着

<math>\begin{align}

f_pf_k(1\pm f_{p'})(1\pm f_{k'})-f_{p'}f_{k'}(1\pm f_p)(1\pm f_k)=0 \end{align}</math> 换言之

<math>\begin{align}

\frac{f_p}{(1\pm f_p)}\frac{f_k}{(1\pm f_k)}=\frac{f_{p'}}{(1\pm f_{p'})}\frac{f_{k'}}{(1\pm f_{k'})} \end{align}</math> 两边取对数后,是和守恒的形式.按Landau统计力学的思路,每一项对应两体碰撞过程的守恒量的函数,即

<math>\begin{align}

\ln\frac{f_p}{(1\pm f_p)}=-\beta{(\epsilon_p-\mu)} \end{align}</math> 而对应的<math>f_p</math>分布函数必然对应费米狄拉克或者波色爱因斯坦统计.

<math>\begin{align}

{f_p}=\frac{1}{exp(-\beta{\epsilon_p-\mu})\pm 1} \end{align}</math> 从而反过来解释了玻色子的<math>(1+f_k)</math>来源.

(10) 流体力学部分的定义

重要的定义,需要一一验证.自由度数目,对称的4X4能动张量,10=4+1+5

(11) 因为积分元必须是洛伦兹标量,所以在质壳上的积分导致分母上有<math>k^0</math>因子.又见(3)下面第二行对应表达式.故而从量纲上,可以理解<math>\langle E_k\rangle_0=n</math>.它的洛伦兹变换方式对应一个四矢量的零分量.

同理,能动张量的零分量,在随动系中除掉<math>k^0</math>因子后对应于<math>\langle E_k^2\rangle_0=nm</math>.它的洛伦兹变换方式对应一个两阶张量的零零分量.又参见比如A First Course in General Relativity by Bernard Schutz的(4.12)附近的讨论.

(13) 这里的<math>\tilde{f}</math>的定义在(3)下面.虽然不清楚为什么要这样定义,但是<math>\phi</math>还是承载了对平衡态偏离的信息,且刨去了<math>f\tilde{f}</math>因子.

(17) 这里对<math>k^\mu</math>的展开,由于(7),包括对速度平行部分<math>E_k</math>和垂直部分<math>k^{<\mu>}</math>的展开.这里先按后者展开,而展开系数为前者的函数.

对<math>k^{<\mu>}</math>采用的基,对应洛伦兹群不可约表示的,满足正交归一关系(16)的基(14)来展开.

(21) 为了满足归一条件(20),这里归一条件中又反过来和<math>k</math>有关.

(26) 这是本文对分布函数展开的一般表达式.注意到因为展开系数其实是<math>\rho</math>,要求解的其实是<math>\rho</math>对时间的依赖关系,比如(29-31).而他们是与宏观物理量直接有关的,见(41).

(42-43) 这里展开是关于两个小量进行的.第一个是关于Knudsen数,是指微观线度(平均自由程)远小于宏观线度(压强发现明显变化的尺度).第二个在文中被称为Reynolds数的倒数,是说系统偏离平衡态不远,即粘滞项远小于能量密度和压强.

(51-52) 这里由对矩阵对角化并且保留最小本征值的迟豫方程,对于其他本征值直接用渐进值代入.这个近似和之前他们的文章arXiv:1102.4780相关,其实,这里的一些表明上相对复杂的计算的物理实质在文章arXiv:1102.4780给出了更加清晰的分析.

(68) 这里在形式上把碰撞项(通过矩阵本征值)和粘滞系数联系起来了.

(69-70) 这里考虑固定散射截面的特殊情况,具体计算碰撞项相关的(45-46)中涉及的积分.

本文档除了包括推导,疑惑外,做读书重点的记录 <math></math>

A new exact solution of the relativistic Boltzmann equation and its hydrodynamic limit, arXiv:1408.5646, by Gabriel S. Denicol, Ulrich Heinz, Mauricio Martinez, Jorge Noronha, and Michael Strickland

本文档除了包括推导,疑惑外,做读书重点的记录 <math></math>

Studying the validity of relativistic hydrodynamics with a new exact solution of the Boltzmann equation, arXiv:1408.7048, by Gabriel Denicol, Ulrich Heinz, Mauricio Martinez, Jorge Noronha, Michael Strickland

这两篇文章讨论了把流体力学满足Gubser(共形)对称性的流体力学方程的解析解推广到Boltzmann方程的解.

由于这个解满足Bjorken对称性且包含对称的横向流,同时因为这个解是解析的,所以对于Boltzmann方程和很多流体力学方程的数值模型有引导作用.

(5) Gubser解涉及到两个坐标系之间的映射关系.在一个坐标系中对称性的表现很直观,而在另一个平直空间的坐标系对应物理问题的实际坐标空间.

本文档除了包括推导,疑惑外,做读书重点的记录 <math></math>

The anisotropic non-equilibrium hydrodynamic attractor, arXiv:1709.06644, by Michael Strickland, Jorge Noronha and Gabriel S. Denicol

这个工作是讨论了Boltzmann方程对应的各向异性流体力学aHydro以及粘滞流体力学的吸引子解.可视为是arXiv:1503.07514的更多的例子和数值验证.

首先,这个工作是从弛豫时间近似的Boltzmann方程出发,推导出对应的aHydro方程,然后和已经的流体力学方程的解比较MIS,DNMR,NS.比较对象还包括原Boltzmann方程的严格解,但是他的差别和aHydro的在慢滚近似下的解差别很小.最后的结果包括两个方面,第一是数值上展示了aHydro和DNMR模型的确存在吸引子解;第二是比较各模型的吸引子解和严格解的差别.

(4-5)是从Boltzmann方程出发推导出相应的粘滞流体力学方程.这里的解考虑了Bjorken和弓星对称性.

(15)是从Boltzmann方程出发推导出相应的aHydro的运动方程.值得指出,粘滞流体力学模型MIS或者DNMR其实就是各向异性的,所以(20)就是把两者统一以<math>\Pi</math>为自变量进行比较.由于(16)和(17),文中指出,这个方程实际上牵涉到无限阶的倒雷诺数展开.

(32-33)对MIS和DNMR模型,通过换元(27-29),两个运动方程只有一个与模型有关.通过进一步的变量变换最终得到了运动方程.在慢滚近似下,吸引子解可以解析的得到.慢滚近似相当于假设解对时间参数的导数足够小,所以解可以按其对时间参数的导数展开,并引入<math>\delta</math>作为展开参数.最低阶的展开相当于说在方程中忽略解对时间参数的导数,由于方程本身依赖于时间,所以这样得到的吸引子解仍然是时间参数的函数.数值结果说明,这个近似非常精确.

(39)接着讨论了aHydro的吸引子解,这在慢滚近似下没有简单的解析解.

(41)这是Boltzmann方程的严格解.

本文档除了包括推导,疑惑外,做读书重点的记录 <math></math>

Analytical attractor and the divergence of the slow-roll expansion in relativistic hydrodynamics, arXiv:1711.01657, by Gabriel S. Denicol and Jorge Noronha

这篇文章可以看做这两个作者之前文章arXiv:1608.07869 [nucl-th]在PRL不中后的加强版续作.

文章的主要内容是继续使用他们提出的和trans-series相关的展开办法得到的解和attrator解析解相比较,证明其实他们的解是唯一不发散的解.

另外,继承Jorge文章的一般风格,这篇文章清楚的讨论了得到attractor解得几种做法.比如大时间数值时的渐进解;初始时刻边界条件数值解;慢滚展开解等.清晰的说明了这些数值或者解析解法的具体过程.

本文中有些不清楚的地方.流体力学的因为和牛顿第二定律等价所以对应对时间的两阶方程.这样的方程对某一平衡态(吸引子)的微扰方程仍然是两阶的,所以趋向吸引子模式是准正规模是可以理解的.但是相关文章中通过变量变换,把方程写为两个一阶方程,其中一个退偶,导致趋向吸引子的方程为一阶,故准正规模的频率为纯虚.作为比较,简谐振动方程可以写为速度和坐标的一阶方程,但是两个方程始终互相耦合,微扰的运动方程仍然为两阶,可能出现准正模式.作者提出的新的展开方法,得到的是解存在初值和解本身的偏差,但是不清楚为什么相对误差却有准正模式的形式.

本文档除了包括推导,疑惑外,做读书重点的记录 <math></math>