Файл: Динамика и управление ядерным ракетным двигателем [Текст] 1974. - 253 с.pdf

Переход к одномерному описанию может быть осуществ­ лен либо феноменологически (исходя из понимания сути происходящих явлений сразу можно дать описание по средним величинам), либо формально. В том и другом слу­ чаях конечный результат получается одинаковым, но формальный путь имеет то преимущество, что позволяет проследить связь вводимых средних величин с локальными величинами, входящими в уравнения исходной системы.

Процедура формального перехода к одномерным урав­ нениям и некоторые дополнительные соображения изложе­ ны в приложении 1. Ниже приведена одномерная система уравнений, полученная таким путем из системы (2.13):

S {dp! dt) -[- (dGIdz) =- 0;

dG/dt + d(GV)ldz= —S {dp/dz)~l1GV;

(2.14)

S TpTcT (dTJdf) -- knpT S v— Па (Тт— T);

G = pSV-, p = pRT;

R — R (T, p); ср = ср{Т,р)- а — а (T, G),

где G — расход теплоносителя; S — общая площадь сече­ ния каналов для теплоносителя; S T — площадь сечения твэла, — приведенный коэффициент гидравлического

сопротивления • ( | г = т/ (4 ^5 ), где £ — коэффициент гидравлического сопротивления), П — параметр; а — коэф­ фициент теплопередачи. Индекс осреднения в системе

(2.14) опущен.

Одномерная система уравнений, соответствующая си­ стеме (2.12), получается из системы (2.14) отбрасыванием всех частных производных по времени в уравнениях тепло­ носителя.

Иногда систему (2.14) целесообразно использовать в не­ сколько иной форме, введя понятие импульса. Для этого, пользуясь известным выражением для импульса

преобразуем уравнение движения в системе (2.14). Диф­ ференцируя (2.15) по г и складывая полученное выражение с уравнением движения, имеем

Заменяя в системе (2.14) уравнение движения уравне­ нием (2.16) и добавляя для полноты системы соотношение (2.15), получаем вместо (2.14) следующую систему урав­ нений:

5 (ЗрIdt) + 3G/dz = 0;

dGldt + d£fldz = — ^GV;

^ r ( f + v f ) = n « ( T , - T y ,

(2.17)

S r pTcT (dTyldt) = knpTST— Па (Тт— T); p = pRT\ G = pSV; Cf = GV + Sp-,

R = R (T, p); cp = cp (T, p); а = а (7\ G).

Особенно существенно система (2.17) упрощается, когда оправдано квазистационарное описание теплоносителя и последний можно считать совершенным газом. В этом слу­ чае, используя известные соотношения [6]:

Вспоминая затем [6], что 3 ~ р* Sf (X), где р* — давление торможения; f (X) — известная газодинамическая функция, и, учитывая, что при малых Я можно принять / (X) = 1, вместо (2.18) получаем уравнение

Газодинамическая функция X/z (X), входящая в А, при малых X может быть аппроксимирована выражением

где q (X) — газодинамическая функция, которую можно выразить через параметры торможения, пользуясь изве­ стным соотношением

34

где т — величина, характеризующая теплоноситель и за­ висящая от k и R. Учитывая соотношения (2.20) и (2.21), из (2.18) получаем следующую форму уравнения движения:

циент. Заметим, что в уравнении (2.22) скорость оказа­ лась исключенной. Учитывая далее, что при малых к мож­

STрт ст (dTT/dt)=knpTST— Па (Гт— Т*).

Вследствие малости к параметры торможения р*, Т* могут

быть заменены в системе (2.23) статическими параметрами

Р, т.

В заключение отметим, что системы уравнений (2.14), (2.17), (2.23) могут быть с некоторыми изменениями при­ менены для описания процессов в других элементах кон­ струкции реактора (замедлитель, корпус твэла и т. д.).

§ 3. Динамика системы твэл — теплоноситель вблизи номинального режима

жима будем изучать, линеаризуя одномерные уравнения, полученные в § 2, и составляя выражения для частотных характеристик реакций выходных координат твэла на входные воздействия. При линеаризации будем исходить из уравнений (2.17), внося в них следующие дополнитель­ ные упрощения:

1) пренебрежем членами с производными по t в урав­ нениях динамики теплоносителя; это соответствует предпо­ ложению о быстром в сравнении с тепловыми процессами затухании газодинамических процессов (некоторые асим­ птотические свойства этих процессов изучаются в при­ ложении 2);

2) пренебрежем кинетической энергией теплоносителя в сравнении с его энтальпией;

3)зависимость коэффициента теплоотдачи а от расхода

итемпературы теплоносителя примем (согласно работе [7]) имеющей вид

(ср, К и р считаем зависящими лишь от температуры тепло­ носителя);

4) будем пренебрегать изменением в динамике геомет­ рических размеров и плотности твэла, а также коэффициен­ та kn в выражении для тепловыделения: q = £n«pTS T; 5) скорость потока теплоносителя на номинальном ре­ жиме будем считать меньшей изотермической скорости

звука [8].

Передаточные функции теплоносителя. Упрощения 1 и 2 позволяют решать уравнения нагрева твэла и теплоносите­ ля независимо от остальных уравнений, входящих в (2.17). Линеаризация этих уравнений дает

(2.25)

где Ф = Д777’б; # т = А 7\/Т б — отклонения температур соответственно теплоносителя и твэла, отнесенные к не­ которому постоянному значению Тб; ѵ = Аq/q0\ у = АG'G0— относительные отклонения соответственно тепловыделения и расхода теплоносителя (согласно предположениям 1 и 4, они не зависят от х)\ Ѳ— перепад температур между твэлом и теплоносителем; индекс «О» означает «на номинальном режиме».

При получении (2.25) использованы соотношения для

36

Переходя в (2.25) к преобразованию Лапласа при нуле­ вых начальных условиях и исключая изображение От, получаем обыкновенное дифференциальное уравнение 1-го порядка относительно изображения О:

[{рІѵ0) + А{х,р)]Ъ= Ы {х,р)Ъ -В{х,р) у, (2.27)

где черта сверху означает изображение соответствующей величины по Лапласу;

({>; = ATj/T6 — отнесенное к Тб отклонение температуры теплоносителя на входе в реактор) дает выражение для за­ висимости О (х, р) от изображений входных воздействий Ог, ѵ и у:

s

37

X = х/1 (здесь Гб—температура теплоносителя на выходе твэла). Представленные характеристики дают основание для аппроксимации динамических свойств теплоносителя

Рис. 2.6. Частотная характеристика реакции температуры теплоносителя на изменение расхода.

характеристиками простейших звеньев — не выше 2-го порядка.

Усредняя эти функции по х с весом Фо (х) [2] (где Фо (х) — интенсивность потока тепловых нейтронов в но­ минальном режиме), получаем передаточные функции

Wri (р), Wjn (р) и Wtg (р) соответственно от ifit ѵ и у

39