ПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ ТЕПЛОПРОВОДНОСТИ ВЫСОКОПОРИСТЫХ МАТЕРИАЛОВ В УСЛОВИЯХ ВОЗДЕЙСТВИЯ ФАКТОРОВ КОСМИЧЕСКОГО ПРОСТРАНСТВА Текст научной статьи по специальности «Физика»

УДК 629.78

doi:10.21685/2307-4205-2022-4-11

ПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ ТЕПЛОПРОВОДНОСТИ ВЫСОКОПОРИСТЫХ МАТЕРИАЛОВ В УСЛОВИЯХ ВОЗДЕЙСТВИЯ ФАКТОРОВ КОСМИЧЕСКОГО ПРОСТРАНСТВА

Н. О. Борщев

Астрокосмический центр Учреждения Российской академии наук Физического института имени П. Н. Лебедева, Москва, Россия moriarty93@mail.ru

Аннотация. Актульность и цели. Рассмотрен метод определения теплопроводности высокопористого материала для условий функционирования при орбитальном полете. Данный метод основан на классическом подходе минимизации среднеквадратичного отклонения между теоретическим и экспериментальным температурным полем в зонах установки датчиков температур. Материалы и методы. Для этого сперва решается «прямая» задача теплообмена для пористого образца при орбитальных начально-граничных условиях при предварительном задании сплайнов, учитывающим зависимость теплопроводности от температуры. Процесс теплообмена внутри пор будет сопровождаться лучисто-кондуктивным теплопереносом. Далее решается задача минимизации целевой среднеквадратичной ошибки. Для этого сначала необходимо найти компоненты градиента целевого функционала шага спуска, так как данная задача оптимизации решается методом спряженных направлений, позволяющим достигнуть требуемой сходимости за минимальное число итерационных приближений. Результаты и выводы. Результаты показали, что при данном уровне температур коэффициент теплопроводности будет варь-

Вт

ироваться в пределах 0,26-0,28 -. Разработанный метод и результаты расчета могут быть использованы в во-

мК

просах проектирования теплозащитных покрытий.

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

Для цитирования: Борщев Н. О. Параметрическая идентификация теплопроводности высокопористых материалов в условиях воздействия факторов космического пространства // Надежность и качество сложных систем. 2022. № 4. С. 98-109. doi:10.21685/2307-4205-2022-4-11

PARAMETRIC IDENTIFICATION OF THERMAL CONDUCTIVITY OF HIGHLY POROUS MATERIALS UNDER THE INFLUENCE OF SPACE FACTORS

N.O. Borshchev

Astro Space Center Institutions of the Russian Academy of Sciences P.N. Lebedev Institute of Physics, Moscow, Russia moriarty93@mail.ru

Abstract. Background. In this paper, a method for determining the thermal conductivity of a highly porous material for operating conditions during orbital flight is considered. This method is based on the classical approach of minimizing the standard deviation between the theoretical and experimental temperature field in the zones of installation of temperature sensors. Materials and methods. To do this, the "direct" heat transfer problem for a porous sample is first solved under orbital initial boundary conditions with preliminary spline construction taking into account the dependence of thermal conductivity on temperature. The heat exchange process inside the pores will be accompanied by radiant-conductive heat transfer. Next, the problem of minimizing the target root-mean-square error is solved. To do this, it is first necessary to find the components of the gradient of the target functional of the descent step, since this optimization problem is solved by the conjugate directions method, which allows achieving the required convergence in the minimum number of iterative approximations. Results and conclusions. The results showed that at this temperature level, the thermal conductivity coefficient will vary between 0,26-0,28 W/mK. The developed method and calculation results can be used in the design of heat-protective coatings.

Keywords: highly porous insulation, inverse thermal conductivity problem, conjugate directions method, thermal regime assurance system, radiation-conductive heat exchange

© Борщев Н. О., 2022. Контент доступен по лицензии Creative Commons Attribution 4.0 License / This work is licensed under a Creative Commons Attribution 4.0 License.

For citation: Borshchev N.O. Parametric identification of thermal conductivity of highly porous materials under the influence of space factors. Nadezhnost' i kachestvo slozhnykh sistem = Reliability and quality of complex systems. 2022;(4):98-109. (In Russ.). doi:10.21685/2307-4205-2022-4-11

Введение

При проектировании теплового режима космических аппаратов зачастую встает вопрос выбора теплоизоляционных материалов, удовлетворяющих заданному температурному режиму, как правило они должны быть низкотеплопроводными. Известно, что наилучшими теплоизоляционными свойствами обладают пористые материалы. Однако необходим адекватный метод определения его теплофизических свойств, поскольку внутри материала присутствуют все виды теплообмена и речь идет об эффективных характеристиках. Таким образом, определение теплопроводности теплоизоляционных высокопористых материалов является актуальной задачей на сегодняшний день.

Постановка задачи теплообмена

Первым шагом в восстановлении коэффициента теплопроводности является составление постановки «прямой» задачи теплообмена конструкции, моделирующей условия эксперимента.

Применительно к данной постановке задачи можно применить упрощенную одномерную постановку задачи теплопроводности при допущении о малом перераспределении теплового потока по поверхности прибора при его одномерном нагреве [1]:

( - П )СЭф (Т)

дТ(x,т) _ д

дт

дх

Ьэф (Т )

дТ (х, т)

дх

В работе [2] получена формула для расчета эффективной теплопроводности пористой системы методом обобщенной проводимости:

ь эф (Т )_ь ск (Т)

I + A

-+V

1 --

+

1+-+-

V2h

здесь

A _■

1

Ь к (Т) +Ь

Ьск (Т )

—^ ( h )2103 4кк km L

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

1) теплопроводностью через контакт;

2) теплопроводностью через среду, заполняющую поры и пространство выступающими шероховатостями контактирующих поверхностей;

3) лучистым теплообменом между поверхностями;

4) конвекцией газа в порах.

При выводе данной формулы были сделаны следующие допущения:

1) линии теплового тока параллельны вертикальным образующим элементарной ячейки;

2) структура пористой системы монодисперсная (средний размер частиц одинаков);

3) форма зерен шарообразная.

Выражение для контактной теплопроводности можно представить в виде

Ь к (Т)_ Л

xdx

Y

Ь „(Т )+Ь ск (Т)

Влияние излучения в порах может быть оценено из выражения [2]

к (Г ) = 2е(Г )2 aT (х, т)3 /п.

Таким образом, выражение для коэффициента теплопроводности через среду, заполняющую поры и пространство выступающими шероховатостями контактирующих поверхностей [2]:

mT ) = 4п (х )е(Г )сГ (x, т)3 1п .

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

Г (0) = Г0; ^ (Г = As (Г )9г + е(Г ^ + qm (т) — е(Г )сГ (4, т)4.

Представим коэффициент теплопроводности исследуемого образца в параметризированном виде, помноженным на соответствующие базисные функции, учитывающие их зависимость от температуры:

м

кэф (ГЬ Тр(Г),

т=1

где Ыт (Г) - базисные функции, описывающие зависимость искомых теплофизических характеристик от температуры. В данной работе используются линейно-непрерывные базисные функции, имеющие следующий вид [3-8]:

^т (Г ) =

0, Г < Гт-1,

Г - Г

Г - Г

т т—1

Г . < Г < Г ,

' т—1 т '

Г — Г

т+1 Г < Г < Г

Г Г ' т—1 т'

т+1 т

0, Г > Гт, т = 1,М.

Рассмотрим восстановление исследуемой характеристики на основе среднеквадратичного функционала невязки между теоретическим и экспериментальным полем температур [7, 8]:

1'

5(Ч)=1 ¡Т[Г(Ч)—Г(хт)] ^

0 '=1

В работе используется метод безусловной минимизации функционала 5 (кр) с помощью

ме-

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

Последовательный алгоритм метода сопряженных градиентов можно представить в следующем виде на примере параметризированного коэффициента теплопроводности [9-13]:

кп+1 = кп +ЛГ+1,

где

Лк "+1 =—вкР

Направление спуска определяется из

рп = вгаа5(кп) + впРР1; в = 0, Р(0)= grad5(к(0)); Ри = ^( ^2-

; grad5 (к(п—1))

RELIABILITY AND QUALITY OF COMPLEX SYSTEMS. 2022;(4) Критерием останова итерационного процесса является выражение

|gradS (Ь(°Г

I

p_1

BS ( дЬ p

<5s

где 5sum - погрешность входных данных, вычисленная в той же метрике, что и целевой функционал:

5 _5 + 5f + 5 ,

sum a f окр'

где 5 f — погрешность входных температур, определяемая следующим выражением:

''max M

5f _ j I5l (T)dт

0 1=1

Для определения компоненты градиента целевого функционала воспользуемся методом неопределенных множителей Лагранжа. Для этого запишем обобщенный функционал Лагранжа с

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

( \

S

Ь *,C *, А

V As Р

_2 J ![Т(ЬР)—Т(x,т)] dT+ J JV(x,т)[Сэф(Т)р

дТ (x, т) д

0 i_1

дт д

Ь эф (Т)

дТ (x, т)

Bx

dxd т +

+

Т , т)

ь эф (т —+£(Т )q+q^+е(Т К (lx, т)4

(т + jJ(0,т)[Т(0,т)-Т0]dт.

Для нахождения компонент сообщим элементарные приращения искомых рассматриваемым теплофизическим характеристикам. Получим следующие зависимости:

Ь ск (Т + ДТ )(1 + Д)Ь * Nm (Т + АТ )_(1 + Д)Ь *

/ ч Bn (т) N (Т) + —dT

m V / дТ

Таким образом, при подстановке данных компонент получим возмущенное уравнение теплопроводности:

Сэф (Т ) + dT

дТ

дТ (x, т) + ДТ (x, т) , \ , \ . .

1 ' -Т (0, т)_ Т0 (0, т) +ДТ (0, т);

дт

Ь (ТдКИ(Т

эф1 ' дТ Bx

дТ (lx, т) + AT (, т) Bx

д^

- А ^ - (('т)4 + (,т)3 (,т)) = > — + —— +--------.

к £ дТ т (1 + А)

Выражение для производной эффективной теплопроводности от температуры примет следующий вид:

( \

дЬэф (Т) _дЬск (Т)

дТ дТ

1

2 • + V2

'1—h v

^ I + А

+

1 + - + -

lx V2h

Ьск (Т )дА '' lY ^

^ I + А

0

здесь

к (Г)дМП+к (Г)Экск (Г) экг.з.(Г) Iь

дА = ск1 ' дГ Л ' дГ дГ 4ккк

дГ ~

103

(| кг.з.(Г ) Г Н }2103

кск (Г ) 4кккт У Ь Выражение для производной по контактной теплопроводности имеет вид

Эк к (Г) = дГ

'Экг.з.(Г) + Экск (Г)

дГ

дГ

^ г11

xdx

(к „(Г ) + к ск (Г ))2

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

Г - до Г ^ Л

к г.з.(Г ) = п (х )Ы„ 12п (х )е(Г )Г (х, т)2 + 3 ^ Г (х, т)

Для получения выражения относительно вариации температуры вычтем из возмущенного уравнения теплопроводности невозмущенное [14-17]. Получим

Сэф (Г )АГМПгП dг

Эт

дГ Эт

_э_

дх

" ( )ЭАГСх^ + дкэф(Г)дГ(x,т) дГ(x,т) + АГ(хт) [кэф (Г ) Эх дГ дх ^] Эх

г (0)=аг ; кЭф (г )дАГМ+дкдфГГ) дГМ г=Т ^ Г+*—4аГ ЗАГ (/х,т)

Эх

дГ

дх2

дГ

(1+ Л)

Теперь можно записать вариацию обобщенного функционала Лагранжа:

55

ттах 1х

к р, с р

= ||[Г (к р)— Г (х, т)] ЛГ + { у(х, т)[Сэф (Г)

дЛГ (х, т)

! р /

Эт

+

+ ЭСэф (Г) дГ (x,т) Г _Э_ дГ Эт дх

кф (Г )ЭЛГ (х,т)+дМП дГ(хт) Г ^ ' дх дГ Эт

dxd т +

д А (Г)

+

{ , т^-д^Г — d т + ^(0, т) | ЛГ (0, ту т. о дГ (1 + Л) 0

Избавимся от знака повторного дифференцирования для выражения

{у^т)

дх

к э ф (Г)

ЭАГ (х, т)

Эх

^ = т)к эф (Г)

ЭАГ(х,т) |1х VЭу(х,т)

Эх

Эх

-х

Хкэф (Г )дАГМ dx = у(х, т)кэф (Г )дАГхт) 10 — к эф (Г )ЛГ (х, т)|0

э ^т) к эф (Г )АГ (х, т) + М^ ^МЦ ^^ АГ (х, т

Эх2 Эф1 ' 1 ' ' дх дГ дх 1 '

+

+

т=1

0

0

Запишем теперь новый вид функционала с учетом сделанного преобразования:

^ (,) = 7 Ъ (и) - Т Т)]АТ + {„( X, т)р[СЭф (Т )ЭАТ(^ +

0 1=1

дт

эт^dT _

дт

дт

д ^т)х эф (T )ДТ (x, т^М^ дт (, т)

дх2

дх

дТ дх

dxd т —

„(*,т)А.ф (Т)~~Тд~х"") 10* А* (Т)АТ(х,т)|

дх

d т +

M

+ Iр

m=1

, ЧЛТ (( , т) дN (Т)Т (( ,т) . " N. (Т )-+ —mV У W ' ЛТ ((х, т)

Эх

дТ дх

+

'max M

J % - т)©.,

N (Т )ЛКк!) + N (Т )Т Сх ■т) ЛТ (( , т)

mV / -\„ -ITT -\„ Vх'/

+9,

Эх

д а (т )

—e-dT■

дТ

дТ дх

+ I^рN. (Т))^ +

дх

(1+^Лу]dт + ^(0,т) J ЛТ(0,т^т.

Для нахождения трех неопределенных множителей Лагранжа перегруппируем члены, стоящие перед вариацией температуры и производной от вариации по температуре, и приравняем их к нулю. Получим

ДТ (X, т): „(X, т)рд£±Ш дПхИ, (Т )+Мх1) дя(Т >3Т (*т)=о;

дТ дт Эх2 дх дТ дх

д А ( т )

ДТ (, т):-М^О я, Эф (т ) + ^, т)^^-= 0;

Эх

дТ

ДТ (0, т): Ç(0, т)+дЧСLT) À эф (Т ) = 0;

дДТ ((х ■ т):

Эх

: Ч((х ■ т) эф (Т ) = 0;

дДТ(0,т) :

Эх

: у(х, т)А эф (Т ) = 0.

Решая эту систему, находим неопределенные множители Лагранжа. Теперь перепишем наш функционал в другом виде:

( \ ттах 1х

5Б

и- и- Л р J J.....

0 0 т=1

М 7ЬТ{лт М

V

= ПдТЛА' + J J х (Т^^ +

0 0 ОА р

дТ (х, т)

:р

р дт

+ £Ар ^ЛТ(х,т) + In. (Т)Ар Щ + ÎÎNm (Т)ЛАр ^]}т-

m=. т.

+

1 дТ " дх .=1

дЛТ(х,т) , дч(х,т)

Ч(х, т)А эф (Т )

Эх

Эх

т

А эф (Т )ЛТ (х, т)|0

m=1

т

M

d т + р ^^ (Т ^^^ +

Эх

ттах A M

+ J ^(/х, т)[?ж [Л A IN. (Т ) +

0 e .=1

(Т ) e дТ

ЛТ(х,т)]-Л9(]dт + ^(0,т) {ЛТ(0,)dт.

Группируя члены относительно вариаций параметризированного коэффициента теплопроводности материала, получим аналитическое выражение компоненты градиента целевого функционала:

д^ Y !rM , дТ (х, т)

дА

= J £>(х,т)Nm (Тdxdт+ J ЩААрN.(Т)

р 0 0 .=1

0 .=1

Т ( ^ т)

дх

d т.

т

m=1

т

0

В качестве регуляризирующего параметра выбирается номер итерации, т.е. на п + 1 итерации целевой функционал должен быть минимален:

5 ( + Дк р) = 5 ( р) - а Л^]2; 5 ( + Дк р) = 1

дк р 2 т=1 0

т (к р )-а сп ёга^5

р

-т (х, т)]2

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

М Ттах 1

асп = V [ ——т. сп ¿1 дт(к,)

Эк,

Таким образом, можно сформулировать принцип работы данного алгоритма идентификации теплофизических параметров. Суть его заключается в следующем:

1. Решается «прямая» задача прогрева конструкции при реализации граничных условий, характерных для условий тепловакуумных испытаний изделия.

2. Далее после результатов тепловакуумных испытаний составляется среднерегуляризируемая интегральная ошибка между теоретическим и экспериментальным полем температур в зонах установки датчиков температур.

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

4. Ищется шаг спуска в методе «сопряженных» направлений, обеспечивающий минимум целевого функционала на следующей итерации.

5. Получаем значения целевых функций на следующей итерации. При выполнении условия останова полученные значения каждого из блоков домножаются на свои характерные базисные функции, считается, что зависимости получены, иначе процесс 1-5 выполняется повторно.

Реализация решения данного алгоритма проиллюстрирована на блок-схеме на рис. 1.

Рис. 1. Блок-схема алгоритма идентификации искомых теплофизических параметров

В численном эксперименте рассматривается образец в виде параллелепипеда, по толщине которого установлено шесть термопар. Все его поверхности теплоизолированные, кроме верхнего основания, на которое падает лучистый интегральный тепловой поток. Таким образом, реализуется одномерный прогрев по толщине материала, моделирующий заданную постановку задачи. В качестве источника теплового потока используется медный линейчатый нагреватель. Результаты определения падающего теплового потока и экспериментального температурного поля представлены на рис. 2, 3.

60

со

5С

О

5 40

с JS О

§ 30

а <и

20

Xi X Л

5 10 £

Г4*- •

-e[qzem+qpereiz] •Asqsoln

20

40 60

Время, мин

80

100

Рис. 2. Удельные падающие тепловые потоки на нагреваемую поверхность образца

Время,

Рис. 3. Температурное поле в местах установки датчиков температур

При итерационном уточнении параметризированной величины коэффициента теплопроводности материала будет также по итерациям восстанавливаться температурное поле, стремясь к своему экспериментальному аналогу. На рис. 4 приведена зависимость температур в точке замера температур для двух крайних моментов времени в зависимости от номера итерации.

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

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

Номер итерации

Рис. 4. Итерационное изменение температурного поля для двух моментов времени

ю

\

\ \

\ \ \

\ \ \

- -1- 1

1,5

2,5

Номер итерации

-М=1

-м=2

3,5

Рис. 5. Итерационное изменение среднеквадратичной ошибки между теоретическим и экспериментальным температурным полем в местах установки датчиков температур

Итерационные изменения коэффициента теплопроводности материала представлены для каждого из двух временных блоков на рис. 6, а изменение коэффициента теплопроводности от температуры - на рис. 7.

о------.---,--,-

1 1,5 2 2,5 3 3,5 4

Номер итерации

Рис. 6. Итерационное изменение коэффициента теплопроводности высокопористого материала от температуры

Рис. 7. Изменение коэффициента теплопроводности высокопористого материала от температуры

Заключение

1. Разработан метод параметрической идентификации коэффициента теплопроводности высокопористых материалов как функций от температуры методом итерационной регуляризации в приближении однонаправленного прогрева для орбитального участка полета.

2. Продемонстрированы результаты данного алгоритма на примере определения коэффициента теплопроводности образца высокопористого теплозащитного покрытия.

3. Результаты показали, что при данном уровне температур коэффициент теплопроводности

Вт

будет варьироваться в пределах 0,26-0,28 -.

мК

4. Данный алгоритм может быть использован и для более широкого температурного диапазона для определения коэффициента теплопроводности высокопористых материалов в изотропном приближении.

Список литературы

1. Залетаев В. М., Капинос Ю. В., Сургучев О. В. Расчет теплообмена космического аппарата. М. : Машиностроение, 1979.

2. Васильев Л. Л., Танаева С. А. Теплофизические свойства пористых материалов. Минск, 1971. С. 17-20.

3. Басистов Ю. А., Яновский Ю. Г. Некорректные задачи в механике (реологии) вязкоупругих сред и их регуляризация // Механика композиционных материалов и конструкций. 2010. Т. 16, № 1. С. 117-143.

4. Бакушинский А. Б., Кокурин М. Ю., Кокурин М. М. Прямые и обратные теоремы для итерационных методов решения нерегулярных операторных уравнений и разностных методов решения некорректных задач Коши // Журнал вычислительной математики и математической физики. 2020. Т. 60, № 6. С. 939-962.

5. Фанов В. В., Мартынов М. Б., Карчаев Х. Ж. Летательные аппараты НПО им. С. А. Лавочкина (к 80-летию предприятия) // Вестник НПО им. С. А. Лавочкина. 2017. № 2/36. С. 5-16.

6. Блох А. Г., Журавлев Ю. А., Рыжков Л. Н. Теплообмен излучением. М. : Энергоатомиздат, 1991.

7. Тулин Д. В., Финченко В. С. Теоретико-экспериментальные методы проектирования систем обеспечения теплового режима космических аппаратов. М. : МАИ-ПРИНТ, 2014. Т. 3. С. 1320-1437.

8. Цаплин С. В., Болычев С. А., Романов А. Е. Теплообмен в космосе. Самара : Самарский университет, 2013. 53 с.

9. Алифанов О. М., Артюхин Е. А., Румянцев С. В. Экстремальные методы решения некорректных задач. М. : Наука, 1988. 288 с.

10. Алифанов О. М. Обратные задачи теплообмена. М. : Машиностроение, 1988. 280 с.

11. Формалев В. Ф. Теплоперенос в анизотропных твердых телах. М. : Физматлит, 2015. 238 с.

12. Васин В. В. Модифицированный метод наискорейшего спуска для нелинейных регулярных операторных уравнений // Доклады Академии наук. 2015. Т. 462, № 3. С. 264.

13. Голичев И. И. Модифицированный градиентный метод наискорейшего спуска решения нелениаризован-ной задачи для нестационарных уравнений Навье-Стокса // Уфимский математический журнал. 2013. Т. 5, № 4. С. 60-76.

14. Формалев В. Ф., Ревизников Д. Л. Численные методы. М. : Физматлит, 2004. 400 с.

15. Формалев В. Ф. Анализ двумерных температурных полей в анизотропных телах с учетом подвижных границ и большой степени анизотропии // Теплофизика высоких температур. 1990. Т. 28, № 4. С. 715-721.

16. Формалев В. Ф. Идентификация двумерных тепловых потоков в анизотропных телах сложной формы // Инженерно-физический журнал. 1989. Т. 56, № 3. С. 382-386.

17. Формалев В. Ф., Колесник С. А. Аналитическое решение второй начально-краевой задачи анизотропной теплопроводности // Математическое моделирование. 2003. Т. 15, № 6. С. 107-110.

1. Zaletaev V.M., Kapinos Yu.V., Surguchev O.V. Raschet teploobmena kosmicheskogo apparata = Calculation of the heat exchange of the spacecraft. Moscow: Mashinostroenie, 1979. (In Russ.)

2. Vasil'ev L.L., Tanaeva S.A. Teplofizicheskie svoystva poristykh materialov = Thermophysical properties of porous materials. Minsk, 1971:17-20. (In Russ.)

3. Basistov Yu.A., Yanovskiy Yu.G. Incorrect problems in mechanics (rheology) of viscoelastic media and their regularization. Mekhanika kompozitsionnykh materialov i konstruktsiy = Mechanics of composite materials and structures. 2010;16(1):117-143. (In Russ.)

4. Bakushinskiy A.B., Kokurin M.Yu., Kokurin M.M. Direct and inverse theorems for iterative methods of solving irregular operator equations and difference methods for solving ill-posed Cauchy problems. Zhurnal vychislit-el'noy matematiki i matematicheskoy fiziki = Journal of Computational Mathematics and Mathematical Physics. 2020;60(6):939-962. (In Russ.)

5. Fanov V.V., Martynov M.B., Karchaev Kh.Zh. Aircraft of S.A. Lavochkin NGO (to the 80th anniversary of the enterprise. VestnikNPO im. S. A. Lavochkina = Bulletin of S.A. Lavochkin NGO. 2017;(2/36):5-16. (In Russ.)

6. Blokh A.G., Zhuravlev Yu.A., Ryzhkov L.N. Teploobmen izlucheniem = Heat exchange by radiation. Moscow: Energoatomizdat, 1991. (In Russ.)

7. Tulin D.V., Finchenko V.S. Teoretiko-eksperimental'nye metody proektirovaniya sistem obespecheniya teplovogo rezhima kosmicheskikh apparatov = Theoretical and experimental methods of designing systems for ensuring the thermal regime of spacecraft. Moscow: MAI-PRINT, 2014;3:1320-1437. (In Russ.)

8. Tsaplin S.V., Bolychev S.A., Romanov A.E. Teploobmen v kosmose = Heat transfer in space. Samara: Samar-skiy universitet, 2013:53. (In Russ.)

9. Alifanov O.M., Artyukhin E.A., Rumyantsev S.V. Ekstremal'nye metody resheniya nekorrektnykh zadach = Extreme methods for solving ill-posed problems. Moscow: Nauka, 1988:288. (In Russ.)

10. Alifanov O.M. Obratnye zadachi teploobmena = Inverse problems of heat transfer. Moscow: Mashinostroenie, 1988:280. (In Russ.)

11. Formalev V.F. Teploperenos v anizotropnykh tverdykh telakh = Heat transfer in anisotropic solids. Moscow: Fizmatlit, 2015:238. (In Russ).

12. Vasin V.V. Modified steepest descent method for nonlinear regular operator equations. Doklady Akademii nauk = Reports of the Academy of Sciences. 2015; 462(3):264. (In Russ.)

13. Golichev I.I. Modified gradient method of the steepest descent of the solution of the non-leniarized problem for unsteady Navier-Stokes equations. Ufimskiy matematicheskiy zhurnal = Ufa Mathematical Journal. 2013;5(4):60-76. (In Russ.)

14. Formalev V.F., Reviznikov D.L. Chislennye metody = Numerical methods. Moscow: Fizmatlit, 2004:400. (In Russ.)

15. Formalev V.F. Analysis of two-dimensional temperature fields in anisotropic bodies taking into account mobile boundaries and a large degree of anisotropy. Teplofizika vysokikh temperature = Thermophysics of high temperatures. 1990;28(4):715-721. (In Russ.)

16. Formalev V.F. Identification of two-dimensional heat flows in anisotropic bodies of complex shape. Inzhenerno-fizicheskiy zhurnal = Engineering and Physics journal. 1989;56(3):382-386. (In Russ.)

17. Formalev V.F., Kolesnik S.A. Analytical solution of the second initial boundary value problem of anisotropic thermal conductivity. Matematicheskoe modelirovanie = Mathematical modeling. 2003;15(6):107-110. (In Russ.)

References

Информация об авторах / Information about the authors

Никита Олегович Борщев

кандидат технических наук, ведущий инженер, Астрокосмический центр Учреждения Российской академии наук Физического института имени П. Н. Лебедева (Россия, г. Москва, Ленинский пр-т, 53) E-mail: moriarty93@mail.ru

Nikita O. Borshchev

Candidate of technical sciences, leading engineer, Astro Space Center Institutions of the Russian Academy of Sciences P.N. Lebedev Institute of Physics (53 Leninsky avenue, Moscow, Russia)

Авторы заявляют об отсутствии конфликта интересов / The authors declare no conflicts of interests.

Поступила в редакцию/Received 15.03.2022 Поступила после рецензирования/Revised 18.04.2022 Принята к публикации/Accepted 16.05.2022