freestream milne

Фристрим-Милн (с) Дерек Эверетт

Цель

freestream-milne численно решает бесстолкновительное уравнение Больцмана, предполагая, что

1. частицы не имеют массы
2. начальная функция распределения изотропна в пространстве азимутальных импульсов $px, py$
3. начальное распределение является гауссовским $(y - eta)$ в 3D случае или дельта-функция в $(y-eta)$ в 2D случае.

Учитывая начальный профиль плотности энергии $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 (float) используется для создания 3D-профиля плотности энергии из 2D MC Glb (только для блока superMC)

ETA_FLAT (float) используется для создания 3D-профиля плотности энергии из 2D MC Glb (только для блока superMC)

СИГМА (плавающее) стандарт. отклонение функции распределения одной частицы 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) количество точек сетки в эта

DIM_RAP (int) количество точек сетки по быстроте для суммы Римана

DIM_PHIP (int) количество точек сетки в phi_p для суммы Римана

Шаг решетки DX (плавающее) в x [фм]

DY (float) шаг решетки в y [фм]

DETA (float) шаг решетки в эта

DRAP (float) шаг решетки по быстроте (дифференциал) УСТАРЕЛО

DTAU (float) соответствующий временной интервал бесплатной потоковой передачи [fm/c]

TAU0 (float) время начала бесплатной потоковой передачи [fm/c]

EOS_TYPE (int) 1: конформный EoS, 2: параметризация Вупперталь-Бхудапешт

E_FREEZE (плавающее число) [ГэВ/фм^3] используется для расчета доли энергии, содержащейся на поверхности с плотностью энергии > E_FREEZE