freestream milne

自由流米尔恩 (c) 德里克·埃弗雷特 (Derek Everett)

目的

freestream-milne 数值求解无碰撞玻尔兹曼方程，假设：

1. 粒子没有质量
2. 初始分布函数在方位角动量空间中是各向同性的 $px, py$
3. 初始分布是高斯分布 $(y - eta)$在 3D 情况下，或 delta 函数 $(y-eta)$在二维情况下。

给定初始能量密度分布 $T^{00}$ ，它计算应力能量张量 $T^{munu}$随后，还会发生粘性流体动力分解。

OpenMP用于加速计算。

用法

参数和选项可以在 freestream_input 中设置，该输入在运行时读入。它们的描述如下 OUTPUTFORMAT (int) 1 : 磁盘 , 2 : 写入内存中的向量

BARYON (int) 0：自由流能量密度，1：自由流重子密度和能量密度

IC_ENERGY (int) 1：平滑高斯，2：MC 高斯，3：读取制表符分隔文件的能量密度，4：读取 superMC 块格式，5：从内存中的向量读取

IC_BARYON (int) 1：平滑高斯分布，2：MC 高斯分布，3：读取制表符分隔文件的重子密度

ETA_WIDTH（浮点）用于从 2D MC Glb 创建能量密度的 3D 剖面（仅适用于 superMC 块）

ETA_FLAT（浮点）用于从 2D MC Glb 创建能量密度的 3D 剖面（仅适用于 superMC 块）

SIGMA（浮点）标准。单个粒子分布函数的偏差 f ~ exp((y - eta)^2 / 2 SIGMA^2)

SIGMA_B（浮点）标准。重子电荷一粒子分布函数的偏差 f_B ~ exp((y - eta)^2 / 2 SIGMA_B^2)

DIM_X (int) x 中的网格点数

DIM_Y (int) y 中的网格点数

DIM_ETA (int) eta 中的网格点数

DIM_RAP (int) 黎曼和的快速网格点数

DIM_PHIP (int) 黎曼和的 phi_p 中的网格点数

DX（浮点）晶格间距，单位为 x [fm]

DY（浮点）y 中的晶格间距 [fm]

DETA（浮点）晶格间距（以 eta 为单位）

DRAP（浮动）快速晶格间距（微分）已弃用

DTAU（浮点）自由流适当时间间隔 [fm/c]

TAU0（浮点）免费流媒体开始适当时间 [fm/c]

EOS_TYPE (int) 1：保形 EoS，2：Wuppertal-Bhudapest 参数化

E_FREEZE (float) [GeV/fm^3] 用于计算能量密度> E_FREEZE 的表面内包含的能量分数