Описание программы. 
Вязкое затухание звуковых волн в сильных центробежных полях

Структура программы довольно проста. Она имеет цикл по частоте возмущений, который начинается с того, что считываются входные параметры, приведённые в таблице 2, далее по формулам (28) и (29) вычисляется распределение плотности и давления по радиусу, после чего решается система однородных дифференциальных уравнений (27), которая после линеаризации [16] принимает вид:

,.

,

,

,

,

,

.

Таблица 2.1.

Параметр	Значение
M	352 грамм моль-1
rвнеш	0,065 м
rвнут	0,0001 м
?	2р 1700 с-1
T0	300 K
P	10 665 Па
K	2р 10 — 2р 200 м-1
cp	385 Дж кг-1 К-1
м	1,83 10-5 Па с
J	2р 868 — 2р 17 376 с-1
Л	0,0061 Втм-1 К-1
R	8,314 462 м2 кг с-2 К-1 Моль-1
G	1,067

где M — молярная масса гексофторида урана (UF₆), r_внеш -внешняя граница расчетов, r_внут — внутренняя граница расчетов, щ — круговая частота вращения ротора, T₀ — температура на внешней границе, p — давление на внешней границе, k — волновое число задаваемых возмущений, c_p — удельная теплоемкость гексофторида урана UF₆ при постоянном давлении, м -динамическая вязкость гексофторида урана UF₆, j — частота задаваемых возмущений, л — теплопроводность гексофторида урана UF₆, R — универсальная газовая постоянная, g — показатель адиабаты.

Далее по формуле:

.

где.

считается суммарная энергия волны.

По полученным значениям строится график зависимости энергии волны от частоты вращения ротора.

Далее, предполагая вид кривой:

,

где a и b — искомые параметры, проводится интерполяция полученных значений к данному виду.

Рис. 8.