freestream milne

freestream-milne (c) Derek Everett

Objetivo

freestream-milne resuelve numéricamente la ecuación de Boltzmann sin colisiones suponiendo que

1. las partículas no tienen masa
2. la función de distribución inicial es isotrópica en el espacio de momento azimutal $px,py$
3. la distribución inicial es gaussiana $(y - eta)$ en el caso 3D, o función delta en $(y-eta)$ en caso 2D.

Dado un perfil de densidad de energía inicial $T^{00}$ , calcula el tensor de energía de tensión. $T^{munu}$ en un momento posterior, y también una descomposición hidrodinámica viscosa.

OpenMP se utiliza para acelerar el cálculo.

Uso

Los parámetros y opciones se pueden configurar en freestream_input, que se lee en tiempo de ejecución. Se describen a continuación OUTPUTFORMAT (int) 1: disco, 2: escribir en vectores en la memoria

BARYON (int) 0: corriente libre la densidad de energía, 1: corriente libre la densidad bariónica y la densidad de energía

IC_ENERGY (int) 1: Gaussiano suave, 2: MC Guassiano, 3: leer archivo delimitado por tabulaciones para densidad de energía, 4: leer formato de bloque superMC, 5: leer desde vector en la memoria

IC_BARYON (int) 1: Gaussiano suave, 2: MC Guassiano, 3: leer archivo delimitado por tabulaciones para densidad bariónica

ETA_WIDTH (flotante) utilizado para crear un perfil 3D de densidad de energía a partir de 2D MC Glb (solo para el bloque superMC)

ETA_FLAT (flotador) utilizado para crear un perfil 3D de densidad de energía a partir de 2D MC Glb (solo para el bloque superMC)

SIGMA (flotador) estándar. desviación de la función de distribución de una partícula f ~ exp((y - eta)^2 / 2 SIGMA^2)

SIGMA_B (flotante) estándar. desviación de la función de distribución de una partícula de cargas bariónicas f_B ~ exp((y - eta)^2 / 2 SIGMA_B^2)

DIM_X (int) número de puntos de la cuadrícula en x

DIM_Y (int) número de puntos de la cuadrícula en y

DIM_ETA (int) número de puntos de la cuadrícula en eta

DIM_RAP (int) número de puntos de la cuadrícula en rapidez para la suma de Riemann

DIM_PHIP (int) número de puntos de la cuadrícula en phi_p para la suma de Riemann

Espaciado de red DX (flotante) en x [fm]

Espaciado de red DY (flotación) en y [fm]

Espaciado de celosía DETA (flotante) en eta

Espaciado de celosía DRAP (flotante) en rapidez (diferencial) DESPRECADO

Intervalo de tiempo adecuado de transmisión gratuita DTAU (flotante) [fm/c]

TAU0 (flotante) inicio de la transmisión gratuita en el momento adecuado [fm/c]

EOS_TYPE (int) 1: EoS conforme, 2: parametrización de Wuppertal-Bhudapest

E_FREEZE (flotante) [GeV/fm^3] usado para calcular la fracción de energía contenida dentro de la superficie de densidad de energía > E_FREEZE