/
Zlalex
/
PLASMA_SIM
Обзор
Документация
Войти
/
Zlalex
/
PLASMA_SIM
Код
Запросы
0
Задачи
Вики
Пакеты
0
Релизы
0
Аналитика
Безопасность
main
TempDestribution.py
18 строк
760 B
Zlalex
Added correct coding-decoding files for solver
13 дек 2024, 12:37
13 дек 2024, 12:37
6cc8235
Код
Авторство
О чём код?
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 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 return (R_DT*MevToK(17.6) + R_DD * MevToK(7.3), n_D, n_T, n_He)