freestream milne

freestream-milne (c) Derek Everett

Zweck

freestream-milne löst die kollisionsfreie Boltzmann-Gleichung numerisch unter der Annahme, dass

1. die Teilchen sind masselos
2. Die anfängliche Verteilungsfunktion ist im azimutalen Impulsraum isotrop $px, py$
3. Die anfängliche Verteilung ist eine Gaußsche Verteilung $(y - eta)$ im 3D-Fall oder Delta-Funktion in $(y-eta)$ im 2D-Fall.

Gegeben sei ein anfängliches Energiedichteprofil $T^{00}$ , berechnet es den Spannungsenergietensor $T^{munu}$ zu einem späteren Zeitpunkt, und auch eine viskose hydrodynamische Zersetzung.

Um die Berechnung zu beschleunigen, wird OpenMP verwendet.

Verwendung

Parameter und Optionen können in freestream_input eingestellt werden, das zur Laufzeit eingelesen wird. Sie werden unten beschrieben. OUTPUTFORMAT (int) 1: Festplatte, 2: Schreiben in Vektoren im Speicher

BARYON (int) 0 : Freestream die Energiedichte, 1 : Freestream die Baryonendichte und Energiedichte

IC_ENERGY (int) 1: Glatter Gaussian, 2: MC Guassian, 3: tabulatorgetrennte Datei für Energiedichte lesen, 4: SuperMC-Blockformat lesen, 5: Aus Vektor im Speicher lesen

IC_BARYON (int) 1: Smooth Gaußian, 2: MC Guassian, 3: tabulatorgetrennte Datei für Baryonendichte lesen

ETA_WIDTH (float) wird verwendet, um ein 3D-Profil der Energiedichte aus 2D MC Glb zu erstellen (nur für SuperMC-Block)

ETA_FLAT (Float) wird verwendet, um ein 3D-Profil der Energiedichte aus 2D MC Glb zu erstellen (nur für SuperMC-Block)

SIGMA (Float) std. Abweichung einer Partikelverteilungsfunktion f ~ exp((y - eta)^2 / 2 SIGMA^2)

SIGMA_B (Float) std. Abweichung einer Teilchenverteilungsfunktion der Baryonenladungen f_B ~ exp((y - eta)^2 / 2 SIGMA_B^2)

DIM_X (int) Anzahl der Gitterpunkte in x

DIM_Y (int) Anzahl der Gitterpunkte in y

DIM_ETA (int) Anzahl der Gitterpunkte in Eta

DIM_RAP (int) Anzahl der Gitterpunkte in Schnelligkeit für die Riemann-Summe

DIM_PHIP (int) Anzahl der Gitterpunkte in phi_p für Riemann-Summe

DX (float) Gitterabstand in x [fm]

DY (float) Gitterabstand in y [fm]

DETA (Float)-Gitterabstand in Eta

DRAP (Float)-Gitterabstand in Schnelligkeit (Differential) VERALTET

DTAU (Float) kostenloses Streaming im richtigen Zeitintervall [fm/c]

TAU0 (float) Beginn des kostenlosen Streamings zur richtigen Zeit [fm/c]

EOS_TYPE (int) 1: konformes EoS, 2: Wuppertal-Bhudapest-Parametrisierung

E_FREEZE (float) [GeV/fm^3] wird verwendet, um den Anteil der in der Oberfläche enthaltenen Energie mit einer Energiedichte > E_FREEZE zu berechnen