Описание программы

Структура программы довольно проста. Первым делом производится расчет давления и плотности газа в роторе с помощью распределений (2.6) и (2.7), после чего решается система однородных дифференциальных уравнений (2.1)-(2.5), которая после линеаризации [14] принимает вид:

 

,

,

,

,

,

,

.

с граничными условиями (2.8) или (2.9).

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

,

где

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

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

 

 

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

 

Далее, предполагая вид кривой, полученный в прошлой главе:

,

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

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

Все вышеизложенные вычисления проводятся с заранее заданными шагом по частоте и параметрами, приведенными в таблице2.1.

 

Таблица 2.1. Расчетные параметры центрифуги

Параметр	Значение
M	352 грамм моль-1
rвнеш	0,065 м
rвнут	0,0001 м
Ω	2π 1700 с-1
T0	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 – показатель адиабаты.