/
Zlalex
/
PLASMA_SIM
Обзор
Документация
Войти
/
Zlalex
/
PLASMA_SIM
Код
Запросы
0
Задачи
Вики
Пакеты
0
Релизы
0
Аналитика
Безопасность
dev_univers
TempDestribution.py
22 строки
1 KB
Zlalex
update of velocity computing
01 мар 2025, 13:25
01 мар 2025, 13:25
0118752
Код
Авторство
О чём код?
import numpy as np import matplotlib.pyplot as plt from tqdm import tqdm from Config import * def fusion_rate(n_D:np.array, n_T:np.array, T_e:np.array): """Скорость реакций слияния (дейтерий-дейтерий и дейтерий-тритий).""" reaction_rate_DD = 2.8e-14 * (n_D**2) * np.exp(-1.0 / T_e) reaction_rate_DT = 2.8e-14 * (n_D*n_T) * np.exp(-1.0 / T_e) R_DD = reaction_rate_DD * n_D * n_D R_DT = reaction_rate_DT * n_D * n_T return (R_DD, R_DT) def P_br(T_e:np.array, n_D:np.array, n_T:np.array,n_He:np.array): Z_eff_2 = 1.7 return 8.5 * 1.34e-36 * Z_eff_2 * (n_D + n_T + 2 * n_He) * (n_D + n_T + 2 * n_He) *( Lx/2) **2 * np.pi * Ly * np.sqrt(T_e)/511 def UpdateTemp(T_e:np.array, n_D:np.array, n_T:np.array,n_He:np.array): R_DD, R_DT = fusion_rate(n_D, n_T, T_e) n_D, n_T, n_He = -1*(2* R_DD + R_DT), - R_DT, -R_DT + n_He dE = R_DT*17.6 * 1.6e-13 + R_DD * 7.3 *1.6e-13 - P_br(T_e, n_D, n_T, n_He) return (dE / k_B, n_D, n_T, n_He)