ТЕОРЕТИЧЕСКАЯ ОЦЕНКА ВРЕМЕНИ ЗАХОЛАЖИВАНИЯ БАКА ЖИДКОГО ВОДОРОДА ПРИ ИСПЫТАНИИ НА ПРОЧНОСТЬ Текст научной статьи по специальности «Строительство и архитектура»

УДК 62-621.2 doi: 10.18698/0536-1044-2021-12-83-89

Теоретическая оценка времени захолаживания бака жидкого водорода при испытании на прочность

Н.О. Борщев, О.А. Юранев

АО «ЦНИИМаш»

Theoretical Estimate of Cooling Time of a Liquid Hydrogen Tank during Structural Tests

N.O. Borschev, O.A. Yuranev

JSC "TsNIIMash"

В России продолжаются работы по созданию ракетно-космических транспортных средств на базе криогенных компонентов топлива (жидких водорода, кислорода и метана). При этом важное значение приобретают вопросы отработки температурной прочности топливных баков. При прочностных испытаниях имитируются эксплуатационные температуры испытуемого объекта, так как температурное состояние влияет на прочность и жесткость конструкции. Исходя из этого при экспериментальной наземной отработке водородных баков их необходимо захолаживать до температуры кипения водорода (20 К). В АО «ЦНИИМаш» ведутся работы по созданию гелиевой системы, способной захолодить крупногабаритные конструкции до температуры 20 К. Так как температура кипения гелия (4 К) ниже, чем у водорода, его можно использовать в газообразном состоянии для захолаживания конструкции. До настоящего времени охлаждение баков происходило только заливкой жидкого азота, поэтому моделирование их температурного состояния при испытаниях проведено только для этого случая. Для определения целесообразности применения подобного метода захолажи-вания выполнена оценка времени охлаждения до 20 К крупногабаритных емкостей с помощью газообразного гелия на примере захолаживания водородного бака, по своим размерам характерного для второй ступени перспективной ракеты-носителя среднего класса.

Ключевые слова: лучисто-конвективный теплообмен, газообразный гелий, алюминиевый бак, конвективная теплоотдача, водородный бак, уравнение теплопроводности

Russian enterprises continue developing rocket and space vehicles based on cryogenic propel-lants, i.e. liquid hydrogen, oxygen, and methane. Hence, the issues of fuel tanks' thermal strength are increasingly important. During structural tests, the operating temperatures of the test object should be simulated, since the temperature condition affects the strength and rigidity of the structure. Consequently, during ground-based experimental tests, hydrogen tanks must be cooled down to 20 K, the boiling point of hydrogen. JSC TsNIIMash is developing a helium system capable of cooling large-sized structures to a temperature of 20 K. Helium can be used in a gaseous state to cool down the structure, since the boiling point of helium, 4 K, is lower than the boiling point of hydrogen. Until now, the tanks were cooled only by filling with liquid nitrogen, therefore the temperature state of the tanks during the tests was simulated only for this case. In order to determine the applicability of the method developed, the cooling time of large-sized containers was estimated by cooling a hydrogen tank, which by its dimensions is typical for an advanced medium-class second stage launcher, to 20 K by gaseous helium.

Keywords: radiant-convective heat exchange, gaseous helium, aluminum tank, convective heat transfer, hydrogen tank, thermal conductivity equation

Отработка температурной прочности топливных баков — важный этап разработки ракетно-космических транспортных средств на базе криогенных компонентов топлива: жидких водорода, кислорода и метана [1]. Специалисты АО «ЦНИИМаш» занимаются созданием гелиевой системы, способной захолодить крупногабаритные баки до температуры 20 К [2-5].

Суть охлаждения бака газообразным гелием заключается в следующем. Через бак прокачивают холодный газообразный гелий (рис. 1), температуру которого можно понизить до 5...6 К. Начальное давление в баке принимают равным атмосферному, считая, что воздушная среда в баке заменена газообразным гелием.

Процесс охлаждения можно осуществить двумя способами: циркуляцией гелия с температурой не выше 20 К и с температурой, изменяемой по ступенчатому закону. Для расчетов выбран первый способ.

Сверху и снизу бака установлены силовые кольца для передачи на него статических усилий. Так как силовые кольца могут соприкасаться со шпангоутами бака, они способны изменить внешний приток теплоты, который необходимо учитывать в расчетах. Обычно такие кольца изготавливают из недорогой и технологичной стали 20 и имеют массу около 5 т.

Для расчетов выбран модельный бак вместимостью примерно 100 м3, выполненный из алюминиевого сплава и состоящий из цилиндрической обечайки диаметром й = 3,8 м и высотой к = 8 м и двух сферических днищ с радиусом кривизны 3 м. Масса бака принята равной 1200 кг. Предполагаем, что он соответствует размерам водородного бака, который может войти в состав второй ступени перспективной ракеты космического назначения среднего класса, изготовленной по тандемной схеме.

Принимаем, что на баке смонтирован (методом напыления) слой теплоизоляции типа ППУ толщиной 0,15 м с постоянным коэффициентом теплопроводности, равным 0,03 Вт/(м-К). Считаем, что теплоизоляция имеет постоянную плотность 90 кг/м3, а теплоемкость, зависящая от температуры, определяется методом интерполяции.

Цель работы — определение времени захо-лаживания водородного топливного ракетного бака при испытании на прочность.

Тепловая физико-математическая модель объекта исследования. Анализ работ [6-11]

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

Такая тепловая модель хорошо зарекомендовала себя при оценке теплового состояния конструкций [6-11], поэтому ее взяли за основу. Тепловую модель доработали в части учета неконтролируемых тепловых потерь через теплоизоляцию бака, которая подвержена сложному лучисто-конвективному нагреву от внешней среды в свободном объеме и зависимости всех теплофизических характеристик от температуры.

Указанный подход позволяет получить более точное распределение температуры по стенке бака, что критически важно при его захолажи-вании газообразным гелием, так как в первые моменты времени за счет большой разницы температур стенки и теплоносителя происходит большой теплосъем от конструкции бака.

Предполагаем, что в бак поступает газообразный гелий с постоянной температурой на входе Твх = 20 К. Главным фактором захолажи-вания стенок бака является конвективная теплоотдача к ядру газообразного гелия, а также

Рис. 1. Схема охлаждения до 20 К крупногабаритного топливного водородного бака 1 с помощью газообразного гелия Не(г)

его расход на входе и выходе через дренажное отверстие.

Схема охлаждения до 20 К крупногабаритного топливного бака с помощью газообразного гелия Не(г), поступающего от гелиевого рефрижератора (ГР), приведена на рис. 1.

Чтобы решить поставленную задачу, используем уравнение теплопроводности для бака и его теплозащитного покрытия (ТЗП) [12].

Проинтегрируем уравнение теплопроводности в декартово — евклидовом пространстве по трем координатам и запишем подынтегральные выражения для каждого оператора дифференцирования.

Для алюминиевой оболочки считаем, что перетекание тепловой энергии происходит только в двух направлениях в силу малости толщины алюминиевой стенки. Таким образом, происходит кондуктивный теплоприток со стороны ТЗП и перетекание тепла вдоль бака:

I-

0дх

Л(ТА1 >

а(тА1)

дТ(М, т)

дТ(М, т)

дх

йх =

дх

-АМт А1 )дт (М,т) дх ^ '

дх

= -ак (ТА1 - ТНе) +

Атзп (ттзп)

->ТЗП

(Т ТЗП — Т А1 ).

I -

1 ду

= —А(ТА1) ду 1 ;

А(Т) дТ(М, т)

дТ(М, т)

ду

ду

--ЧТ ")

йу =

дТ(М, т)

ду

I-

• дх

ЦТТЗП)

= —А(ТТЗП ) дх

— дЦ(ТТЗП)

дх

дТ(М, т)

дх

дТ(М, т)

дх дТ(М, т)

йх =

дх

Атзп (ттзп)

->ТЗП

(ТТЗП — ТА1 ) + етЗпа(ТТЗП4 — Т4р ) —

— акв (ТТЗП, ТСр)(ТТЗП — ТСр),

где еТЗП — излучательно-поглощательная способность ТЗП; а — постоянная Стефана — Больцмана, а = 5,67 -10-8; Тср — температура внешней среды; акв — коэффициент свободной конвективной теплоотдачи от алюминиевой стенки к ядру газообразного гелия.

Для стального кольца считаем, что перетекание тепловой энергии происходит только в одном направлении из-за его сравнительно малой толщины. Со стороны алюминиевой стенки происходит кондуктивная отдача тепловой энергии, с другой стороны — лучисто-конвективный теплообмен с внешней средой:

дх

ЧТ}1)

дТ(М, т)

дх

4

дх дх

Л(Г)

йх = дТ(М, т)

дх

где А и АТЗП — коэффициенты теплопроводности алюминиевой оболочки и теплозащитного покрытия; Т(М, т) — температура бака в произвольной точке декартового пространства; ТА1, ТНе и ТТЗП — температура алюминиевой оболочки, газообразного гелия и ТЗП соответственно; 1х и 1у — длины стенок бака в соответствующих направлениях; ак — коэффициент конвективной теплоотдачи от алюминиевой стенки к ядру газообразного гелия; 8ТЗП — толщина ТЗП.

Для теплозащитного покрытия принимаем, что перетекание тепловой энергии происходит только в одном направлении в силу малого изменения температуры гелия. Со стороны алюминиевой стенки происходит кондуктивная отдача тепловой энергии, с другой стороны — лучисто-конвективный теплообмен с внешней средой:

= ^^ (( — ТГ)+■^^ (Тст — ТА1)+

8А1 8А1

+ е,а[(Тст )4 — Тс4р ] — «кв (Т^т —Тср),

где Тст — температура стального кольца; 8А1 — толщина алюминиевой оболочки; E¡ — излуча-тельно-поглощательная способность материла;

Таким образом, получаем систему обыкновенных дифференциальных уравнений (ОДУ), описывающих тепловое состояние конструкции [13].

ОДУ для изотермических узлов алюминиевой конструкции имеет вид

дТ. А1 (т) с (ТА1 )Р(ТА1 =

—акын (ТА1 — ТНе ) + А(ТА-)(ТА1 —ТА )+

н

А ТЗП (ТТЗП)

А1

Нтзп ТА1 (0 ) = 20К;

(ТТЗП — ТА1 );

а вши _

к — "

Nud X

d^ He

d

акв =

Nu X (T) l '

где С — теплоемкость материала; р — плотность материала; Дт — шаг интегрирования по времени; кА1 и кТзп — высота алюминиевого бака и ТЗП; ТА1 (0) — значение температуры алюминиевой оболочки в начальный момент времени; акын — коэффициент вынужденной конвективной теплоотдачи; Ыий — число Нуссельта для диаметра бака й; АНе — коэффициент теплопроводности газообразного гелия.

Число Нуссельта Ыий находим из выражения [8]

Nui = 0,025Re0'8Pr0

1+|7

0,7

где Ыи — число Нуссельта для высоты бака I; Ие — число Рейнольдса; Рг — число Прандтля.

ОДУ для изотермических узлов ТЗП имеет вид

С(ТТЗП )р(ТТЗП )§тзп дТ'м

Хтзп (ТТЗП)

Дт

ЧТЗП

(TТЗп _T AI ) +

+ Ei а

(TT3n )4 _ T4p _ акв (Ттзп , ТСр )(тзп _ Tcp).

Процесс теплообмена между силовым кольцом и внешней средой характеризуется свободным движением воздуха и таким образом критерий Нуссельта определяется произведением чисел Грасгофа Gr и Прандтля при обтекании вертикальной стенки свободномолекулярным потоком воздуха [8]:

Ыи/ = 0,75^г (Т )Рг (Т )]0,25;

Gri =

= g ßp2 (т тзп _Tcp))3

r|T2

ß =

TH

где XB — коэффициент теплопроводности теплоносителя.

ОДУ для первого и второго изотермических узлов гелия имеют вид

rlT He

С, (T1He )p(T1He ))Д- V^ =

= Ср (T1He )Gbx (H _ T1He )FHe +

I „вын /rrAl _T"He \ с

"Г1ЛК у Ii 11 j -Тконв >

r)T He

cP (T2He )p(T2He ))Д- VHe =

He He He

_ \ 12 / «выЦ 1 вых 12 ) -'He "т"

i авын /TAl _ THe \ F • "Г1Хк l-^i i2 ! гконв >

FHe = nd_• ^конв = %dl; ViHe = Kd h • i = 1,2,

4 4

где С, — изобарная теплоемкость теплоносителя.

Расход паров гелия на входе в бак

Gbx = л/2рвх (рвх _ Р ) ,

где рвх — плотность теплоносителя на входе; рвх — давление внешней среды на входе в бак; р — давление внутри бака

Расход паров гелия на выходе из бака

GBbIX

= Рдр Р

п

Р0

где Рдр — площадь дренажного отверстия; Я — газовая постоянная; Тптах — максимальная температура пара, Тптах = 293 К; у — показатель адиабаты паров гелия; р0 — давление внешней среды на выходе из бака.

Показатель адиабаты паров гелия определяется следующим образом:

где g — ускорение свободного падения; Р — коэффициент расширения пара; ^ — кинематическая вязкость; Т — температура теплоносителя.

Коэффициент расширения пара определяется выражением

1

f

( р ^

р0 Y ,—

Тогда коэффициент свободной конвективной теплоотдачи от стенки алюминиевого бака к ядру гелия

2y р0 2 Y р0 Y+1 Y

Y_1 l Р V l р V

Y

Y_1

Y

Y_1

Рис. 2. Расчетная тепловая схема алюминиевого бака, захолаживаемого газообразным гелием: 1 — алюминиевый бак; 2 — ТЗП; 3 — силовое кольцо; Тм — температура силового кольца; Твых — температура паров гелия на выходе из алюминиевого бака

Рис. 3. Зависимость температуры характерных

изотермических узлов стенки бака Ti ( ), Т2 (♦), Т3 ( ) и силового кольца Т300 ( ) от времени захолаживания т

Температура газообразного гелия на входе в бак в начальный момент времени

ТНе (0 ) = 20 К.

ОДУ для изотермических узлов стального силового кольца имеет вид

дТ ст

С, (Тст )р(Тст )-^8ст = = А(Т ^(Т?! — Тст) + А(Т )(Тст — ТА1) +

h Al

h Al

+ sfo

(Тст )4 - Тс4р ] - акв (Тст, Тср ) (Тст - Тср ).

Данная система уравнений решена методом Рунге — Кутты четвертого порядка точности [14]. Характеристики гелия взяты из работы [15].

Оценим тепловое состояние бака на примере его расчета в характерных изотермических узлах, показанных на рис. 2.

Результаты оценки теплового состояния конструкции и времени захолаживания в характерных изотермических узлах, показанных на рис. 2, приведены на рис. 3.

Как видно из рис. 3, температуры Т\ и Т3 практически совпадают. В точке соприкосновения (стыка) алюминиевого бака с силовым кольцом для первого из них температура составила 78 К, для второго — примерно 79 К. На расстоянии 300 мм от этого стыка температура силового кольца Т300 составляет 140 К, на расстоянии 1000 мм от стыка — 233 К.

Вывод

Анализ полученных результатов показал, что крупногабаритные топливные баки можно охладить до температуры 20 К с помощью газообразного гелия. Причем температура стенки бака достигает этой температуры не более чем за 50 000 с (~14 ч), что является приемлемым с точки зрения временных затрат на прочностные испытания.

Литература

[1] Паничкин Н.Г., ред. ЦНИИмаш. Центр исследований прочности. Историяразвития.

Королев, ЦНИИмаш, 2001. 342 с.

[2] Васюкова Д.А., Колозезный А.Э., Юранев О.А. Использование криогенной гелиевой

системы для имитации эксплуатационных температур при испытаниях на прочность баков жидкого водорода перспективных средств выведения. Космонавтика и ракетостроение, 2012, № 2, с. 179-186.

[3] Васюкова Д.А., Колозезный А.Э., Юранев О.А. Эффективный подход к проведению за-

четных прочностных испытаний криогенных баков перспективных средств выведения. Авиакосмическая техника и технология, 2013, № 1, с. 23-25.

[4] Васюкова Д.А., Колозезный А.Э., Юранев О.А. Квалификация способов расчета захола-

живания крупногабаритной испытательной сборки «криогенного» топливного бака РКН при свободной конвекции газообразного хладагента. Полет, 2015, № 7, с. 18-24.

[5] Юранев О.А. Исследования различных способов захолаживания криогенных топлив-

ных баков изделий ракетно-космической техники. Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение, 2018, № 3, с. 50-57, doi: http://dx.doi.org/10.18698/0236-3941-2018-3-50-57

[6] Горбачев С.П., Попов В.П., Славин М.В. Определение времени захолаживания крио-

генного бака. Известия высших учебных заведений. Машиностроение, 2006, № 5, с. 43-53.

[7] Александров А.А., Бармин И.В., Денисова К.И. и др. Инновационная модель примене-

ния жидкого азота для охлаждения ракетного топлива в емкостях заправочных систем наземных комплексов. Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение, 2017, № 3, с. 4-17, doi: http://dx.doi.org/10.18698/0236-3941-2017-3-4-17

[8] Денисов О.Е., Золин А.В., Чугунков В.В. Методика моделирования охлаждения

компонентов ракетного топлива с применением жидкого азота и промежуточного теплоносителя. Наука и образование: научное издание, 2014, № 3, URL: http://engineering-science.ru/doc/699941.html

[9] Горбачев С.П., Кириенко К.И. Моделирование режимов заправки криогенного ем-

костного оборудования. Вести газовой науки, 2015, № 1, с. 124-132.

[10] Александров А.А., Бармин И.В., Кунис И.Д. и др. Особенности создания и развития криогенных систем ракетно-космических стартовых комплексов «Союз». Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение, 2016, № 2, с. 7-27, doi: http://dx.doi.org/10.18698/0236-3941-2016-2-7-27

[11] Сафонов В.С. Обоснование режимных параметров технологических трубопроводов комплексов СПГ с учетом требований промышленной безопасности. Вести газовой науки, 2017, № 1, с. 83-99.

[12] Бермант А.Ф. Краткий курс математического анализа. Москва, Наука, 1965. 664 с.

[13] Мартыненко О.Г., Соковишин Ю.А. Свободно-конвективный теплообмен. Минск, Наука и техника, 1982. 399 с.

[14] Деккер К., Вервер Я. Устойчивость методов Рунге-Кутты для жестких нелинейных дифференциальных уравнений. Москва, Мир, 1988. 332 с.

[15] Сычев В.В., Вассерман А.А., Козлов А.Д. и др. Термодинамические свойства гелия. Москва, Изд-во стандартов, 1984. 320 с.

References

[1] Panichkin N.G., ed. TsNIImash. Tsentr issledovaniy prochnosti. Istoriya razvitiya [TsNII-

mash. Center of strength studies. History of development]. Korolev, TsNIImash Publ., 2001. 342 p. (In Russ.).

[2] Vasyukova D.A., Kolozeznyy A.E., Yuranev O.A. Using cryogenic helium system for exploi-

tation temperatures imitation in strength tests of liquid hydrogen tanks of advanced launch vehicles. Kosmonavtika i raketostroenie [Cosmonautics and Rocket Engineering], 2012, no. 2, pp. 179-186. (In Russ.).

[3] Vasyukova D.A., Kolozeznyy A.E., Yuranev O.A. Effective approach to control strength test

works for cryogenic tanks of advanced launch vehicles. Aviakosmicheskaya tekhnika i tekhnologiya [Aerospace technology], 2013, no. 1, pp. 23-25. (In Russ.).

[4] Vasyukova D.A., Kolozeznyy A.E., Yuranev O.A. Qualification of method for prediction of

gas free convection cooling of the full-scaledtest article of a launcher cryogenic tank. Polet [Flight], 2015, no. 7, pp. 18-24. (In Russ.).

[5] Yuranev O.A. Investigating various ways of cooling cryogenic fuel tanks of aerospace equip-

ment. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Mashinostr. [Herald of the Bau-man Moscow State Tech. Univ., Mechan. Eng.], 2018, no. 3, pp. 50-57, doi: http://dx.doi.org/10.18698/0236-3941-2018-3-50-57 (in Russ.).

[6] Gorbachev S.P., Popov V.P., Slavin M.V. Determination of cooling time for cryogenic tank.

Izvestiya vysshikh uchebnykh zavedeniy. Mashinostroenie [BMSTU Journal of Mechanical Engineering], 2006, no. 5, pp. 43-53. (In Russ.).

[7] Aleksandrov A.A., Barmin I.V., Denisova K.I., et al. Innovative model of liquid nitrogen ap-

plication for cooling rocket fuel tanks filling the ground-based systems. Izvestiya vysshikh uchebnykh zavedeniy. Mashinostroenie [BMSTU Journal of Mechanical Engineering], 2017, no. 3, pp. 4-17, doi: http://dx.doi.org/10.18698/0236-3941-2017-3-4-17 (in Russ.).

[8] Denisov O.E., Zolin A.V., Chugunkov V.V. Simulation methods of rocket fuel refrigerating

with liquid nitrogen and intermediate heat carrier. Nauka i obrazovanie: nauchnoe izdanie [Science and Education: Scientific Publication], 2014, no. 3, URL: http://engineering-science.ru/doc/699941.html (in Russ.).

[9] Gorbachev S.P., Kirienko K.I. Simulation of the filling modes of cryogenic capacitive equip-

ment. Vestigazovoy nauki, 2015, no. 1, pp. 124-132. (In Russ.).

[10] Aleksandrov A.A., Barmin I.V., Kunis I.D., et al. Characteristic features of creating and developing cryogenic systems of space-rocket launch complex "Soyuz". Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Mashinostr. [Herald of the Bauman Moscow State Tech. Univ., Mechan. Eng.], 2016, no. 2, pp. 7-27, doi: http://dx.doi.org/10.18698/0236-3941-2016-2-7-27 (in Russ.).

[11] Safonov V.S. Substantiation of operating conditions for industrial pipe-lines of liquefied natural gas complexes with provision for safety requirements. Vesti gazovoy nauki, 2017, no. 1, pp. 83-99. (In Russ.).

[12] Bermant A.F. Kratkiy kurs matematicheskogo analiza [Short course of mathematical analysis]. Moscow, Nauka Publ., 1965. 664 p. (In Russ.).

[13] Martynenko O.G., Sokovishin Yu.A. Svobodno-konvektivnyy teploobmen [Free convective heat transfer]. Minsk, Nauka i tekhnika Publ., 1982. 399 p. (In Russ.).

[14] Dekker K., Verwer J.G. Stability of Runge-Kutta methods for stiff nonlinear differential equations. Elsevier Science, 1984. 307 p. (Russ. ed.: Ustoychivost' metodov Runge-Kutty dlya zhestkikh nelineynykh differentsial'nykh uravneniy. Moscow, Mir Publ., 1988. 332 p.).

[15] Sychev V.V., Vasserman A.A., Kozlov A.D., et al. Termodinamicheskie svoystva geliya

[Thermodynamics properties of helium]. Moscow, Izd-vo standartov Publ., 1984. 320 p. (In Russ.).

Информация об авторах

БОРЩЕВ Никита Олегович — инженер второй категории отдела тепловых режимов КА и воздействия ФКП. АО «ЦНИИмаш» (141070, Королёв, Российская Федерация, ул. Пионерская, д. 4, корп. 22, e-mail: [email protected]).

ЮРАНЕВ Олег Александрович — начальник лаборатории отдела температурно-статической прочности. АО «ЦНИИмаш» (141070, Королёв, Российская Федерация, ул. Пионерская, д. 4, корп. 22, e-mail: [email protected]).

Статья поступила в редакцию 17.05.2021 Information about the authors

BORSCHEV Nikita Ologovich — Engineer of the second category, Department of Spacecraft Thermal Modes and the Impact of Space Flight Factor. JSC "TsNIIMash" (141070, Korolev, Russian Federation, Pionerskaya St., Bldg. 4, Block 22, e-mail: [email protected]).

YURANEV Oleg Aleksandrovich — Head of the Laboratory of the Department of Temperature and Static Strength. JSC "TsNIIMash" (141070, Korolev, Russian Federation, Pionerskaya St., Bldg. 4, Block 22, e-mail: [email protected]).

Просьба ссылаться на эту статью следующим образом:

Борщев Н.О., Юранев О.А. Теоретическая оценка времени захолаживания бака жидкого водорода при испытании на прочность. Известия высших учебных заведений. Машиностроение, 2021, № 12, с. 83-89, doi: 10.18698/0536-1044-2021-12-83-89

Please cite this article in English as: Borschev N.O., Yuranev O.A. Theoretical Estimate of Cooling Time of a Liquid Hydrogen Tank during Structural Tests. BMSTU Journal of Mechanical Engineering, 2021, no. 12, pp. 83-89, doi: 10.18698/0536-1044-2021-12-83-89