Описание программы. Структура программы довольно проста

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

 

,

,

,

,

,

,

.

 

Таблица 2.1. – Название таблицы

Параметр	Значение
M	352 грамм моль-1
r внеш	0,065 м
r внут	0,0001 м
Ω;	2π 1700 с-1
T 0	300 K
P	10665 Па
K	2π 10 - 2π 200 м-1
cp	385 Дж кг-1 К-1
μ;	1,83 10-5 Па с
J	2π 868 – 2π 17376 с-1
Λ;	0,0061 Вт м-1 К-1
R	8,314462 м2 кг с-2 К-1 Моль-1
G	1,067

 

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

 

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

,

где

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

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

 

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

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

,

где a и b – искомые параметры,

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

 

 

Рис.8.