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

Структура программы довольно проста. Она имеет цикл по частоте возмущений, который начинается с того, что считываются входные параметры, приведённые в таблице 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 — показатель адиабаты.

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

.

где.

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

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

Рис. 6. Пример расчетного графика зависимости энергии волны от её частоты

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

,

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

Рис. 7.