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 的地表內所包含的能量分數