CC BY
18
УДК 51-74
Сибирский журнал науки и технологий. 2017. Т. 18, № 4. С. 719-726
ЧИСЛЕННО-СТАТИСТИЧЕСКИЙ МЕТОД ДЛЯ РЕШЕНИЯ ЗАДАЧ ТЕПЛООБМЕНА В ТЕПЛОЗАЩИТНЫХ КОНСТРУКЦИЯХ СОТОВОГО ТИПА
С. А. Гусев1, 2, В. Н. Николаев3*
1Институт вычислительной математики и математической геофизики СО РАН Российская Федерация, 630090, г. Новосибирск, просп. академика Лаврентьева, 6 2Новосибирский государственный технический университет Российская Федерация, 630073, г. Новосибирск, просп. Карла Маркса, 20 3Сибирский научно-исследовательский институт авиации имени С. А. Чаплыгина Российская Федерация, 630051, г. Новосибирск, ул. Ползунова, 21
*E-mail: [email protected]
Предложен метод для определения теплового состояния теплозащитной панели сотового типа. Применение теплозащитных материалов такого типа является перспективным направлением при конструировании высокоскоростных летательных аппаратов. Рассматриваемый процесс переноса тепла описывается краевой задачей для уравнения теплопроводности. При этом предполагается, что коэффициенты температуропроводности материалов, из которых состоит панель, являются постоянными величинами, а процесс передачи тепла осуществляется только посредством теплопроводности. Предложенный метод основан на вероятностном представлении решения краевой задачи в виде математического ожидания от функционала случайного процесса диффузионного типа и численном моделировании этого случайного процесса. В более ранних работах предлагалось проводить вычисления с использованием метода Эйлера. В связи с тем, что процесс теплообмена происходит в неоднородной среде, при численном моделировании траекторий случайного процесса используется сглаживание разрывных коэффициентов уравнения теплопроводности на основе интегрального усреднения. При этом в подобластях, в которых среда однородная, моделируемый случайный процесс является винеровским. Этот факт означает, что значительного ускорения счёта можно достичь путем использования в этих подобластях метода случайного блуждания по сферам. Предлагается для моделирования винеровского процесса внутри ячеек сотовой панели использовать метод блуждания по сферам, а для расчётов по каркасу панели и в его окрестности - метод Эйлера. Проведённые расчёты по данным физического эксперимента показали значительную эффективность предложенного комбинированного метода по сравнению с вычислениями только одним методом Эйлера.
Ключевые слова: сотовая панель, задача теплообмена, вероятностное представление, численно-статистическое моделирование, винеровский процесс, метод Эйлера, метод блуждания по сферам.
Siberian Journal of Science and Technology. 2017, Vol. 18, No. 4, P. 719-726
NUMERICAL AND STATISTICAL METHOD FOR SOLVING THE HEAT EXCHANGE PROBLEMS IN THE HEAT-PROTECTIVE CONSTRUCTIONS OF THE HONEYCOMB TYPE
S. A. Gusev1, 2, V. N. Nikolaev3*
institute of Computational Mathematics and Mathematical Geophysics SB RAS 6, Academika Lavrenjeva Av., Novosibirsk, 630090, Russian Federation 2Novosibirsk State Technical University
20, K. Marksa Av., Novosibirsk, 630073, Russian Federation
3SibNIA named after S. A. Chaplygin
21, Polzunovа Str., Novosibirsk, 630051, Russian Federation
E-mail: [email protected]
A method for determining the thermal state of a heat insulation panel of the honeycomb type is proposed in the paper. The use of heat-shielding materials of this type is a promising direction in the design of high-speed aircraft. The considered heat transfer is described by the boundary value problem for the heat equation. It is assumed that the thermal diffusivity coefficients of the materials of the panel are the given constants, and the heat exchange process occurs only due to thermal conductivity. The proposed method is based on the probability representation of the solution of the boundary value problem, which is an expectation of the functional of the random process of the diffusion type, and the numerical simulating this random process. In earlier papers the authors offered to perform the computations by the Euler method. As the heat exchange occurs in inhomogeneous medium, we use smoothing discontinuous coefficients
of the heat equation, based on the integral averaging, in numerical simulation of trajectories of the random process. The simulated random process coincides with Wiener process in subareas in which the medium is homogeneous. This fact means that a significant reduction in computational costs can be achieved by using the method of the random walks on spheres in these subareas. In the paper we propose to use the random walks on spheres in the cells of the honeycomb panel and the Euler method in moving on the frame and its vicinity. The calculations performed according to the data of the physical experiment have shown the high efficiency of the proposed combined method in comparison with calculations performed by only the Euler method.
Keywords: honeycomb panel, heat exchange problem, probability representation, numerical statistical modeling, Wiener process, Euler method, walks on spheres method
Введение. При конструировании авиационной и космической техники важнейшее значение приобретает устройство тепловой защиты летательных аппаратов. Одним из возможных вариантов решения этой проблемы является применение теплозащитных панелей сотового типа [1-5]. Сотовая теплозащитная панель представляет собой конструкцию, состоящую из параллельно расположенных пластин, между которыми заключен тонкий каркас в виде пчелиных сот, заполненных некоторым веществом с низкой теплопроводностью. Основными преимуществами таких средств тепловой защиты являются хорошие тепло-физические свойства и небольшой вес при относительно невысокой стоимости и достаточной прочности. Данная статья посвящена математическому моделированию процессов теплообмена в таких теплозащитных конструкциях и разработке метода расчёта для определения их теплового состояния в реальных условиях.
В работе [6] для расчёта теплопереноса в сотовой панели предлагался численно-статистический метод, основанный на вероятностном представлении решения краевой задачи для уравнения теплопроводности и численном решении стохастических дифференциальных уравнений (СДУ). При этом делалось сглаживание разрывных коэффициентов уравнения теплопроводности на основе интегрального усреднения, а численное решение СДУ осуществлялось методом Эйлера.
С помощью метода Эйлера возможно получение хороших статистических оценок решения уравнения теплопроводности, но большая трудоёмкость такого подхода даже при хорошей возможности распараллеливания и применении суперкомпьютерной техники не позволяет проводить реальные расчёты для достаточно больших интервалов полётного времени.
С другой стороны, если теплофизические свойства материалов, из которых состоит сотовая панель, незначительно изменяются в рассматриваемом диапазоне температур и внутри панели отсутствуют источники тепловой энергии, то возможна модификация алгоритма расчёта, позволяющая значительно уменьшить вычислительные затраты. Действительно, если коэффициенты температуропроводности материалов сотовой панели постоянны, а передача тепловой энергии внутри панели осуществляется исключительно теплопроводностью, то случайный процесс в однородной среде, который присутствует в вероятностном представлении, есть ни что иное, как винеровский процесс. Точка первого выхода винеровского процесса на границу шара имеет на ней равномерное распределение.
Следовательно, при моделировании траекторий винеровского процесса в подобластях с однородной средой можно двигаться не по траекториям, которые получаются методом Эйлера, а по границам шаров, находящихся в этих подобластях. В этом состоит основная идея метода случайного блуждания по сферам [7]. Но при решении нестационарных задач требуется помимо точки первого выхода винеровского процесса на границу шара ещё моделировать время первого достижения границы шара. Как известно [8], распределение времени первого выхода винеровского процесса на границу шара совпадает с распределением времени первого достижения уровня, равного радиусу шара, одномерным бесселевским процессом. Но плотность времени первого достижения бесселевским процессом заданного уровня выражается сложной формулой в виде ряда, содержащего функции Бесселя первого рода и их положительные нули. По этой формуле очень затруднительно построить эффективный алгоритм для моделирования времени первого выхода. С другой стороны, вместо заданного уровня, который требуется достичь, можно подобрать некоторые кривые [9], зависящие от t, такие, что получаются достаточно простые и удобные для моделирования формулы для плотности времени первого достижения бесселевским процессом этих кривых. Этот метод используется в данной работе, и он называется блужданием по движущимся сферам.
В качестве области изменения пространственных переменных х1, x2, х3 в работе принимается прямоугольный параллелепипед G = (-11,11)х (-12, 12)х (0, 13). При этом G есть объединение двух непересекающихся подмножеств: G = G1 и G2, где G1 - подобласть, соответствующая каркасу и пластинам, ограничивающим панель, а G2 есть объединение подобластей, соответствующих ячейкам, содержащим теплоизоляционный материал. Каждая ячейка представляет собой внутренность правильной шестиугольной призмы, боковая поверхность которой параллельна оси Xз. Предполагается, что рассматриваемый процесс теплопередачи происходит на отрезке времени [0,7] и описывается следующей краевой задачей для уравнения теплопроводности:
ди
(
— = У—
dt £ dXj
Л
где
b( x) =
иг \ ди b( x)^T
дхг
G
b2, x е G2
(1)
u(x,0) = ф(x) ,
du dx
du
dx2
= 0,
xi =-/i
= 0,
x2
du dx
du
dx2
= 0,
xi =i
= 0,
X2 = /2
du dx3
= a,(t) (u - ) ) =
x3 =0
-X
du dx3
= a2(t) (u - ^(t)).
(2)
(3)
(4)
(5)
(6)
x3 =/
T-h
свойством
J Jh i(u(x,t + h) -u(x,t))2dxdt ^ 0
h ^ 0, с нормой
uQ = max u (x, t) L +
1 qt 0<t <т" iil2 ;
J u2dxdt
\вт
при
(7)
ны можно, например, рассматривать его сглаживание в окрестности поверхностей разрыва с помощью интегрального усреднения с бесконечно дифференцируемым финитным ядром [11]:
Ь(т)(X) = | шрт (|х - у |)Ь(у)ду =
Iх-у| <Рт
:рт J «
Ут \x-y\<рт
i] x-y1 ^
b( y)dy.
(8)
В (1)-(6) использованы следующие обозначения: Ф - начальное распределение температуры в панели; X - коэффициент теплопроводности углепластика; Ь1, Ь2 - коэффициенты температуропроводности
углепластика и воздуха соответственно; а1, а2 -коэффициенты теплообмена между панелью и внешней средой; ^, у2 - температура внешней среды у нижней и верхней сторон панели соответственно.
Сглаживание коэффициентов. Поскольку теплозащитная панель состоит из двух материалов с различными теплофизическими свойствами, то в задаче (1)-(6) мы имеем уравнение теплопроводности с разрывным коэффициентом температуропроводности. Достаточно полное математическое исследование краевых задач для параболических уравнений с разрывными коэффициентами дано в книге [10]. Там же приведены условия существования и единственности обобщенного решения (теорема 5.1, гл. III) в пространстве функций К20,1(2г) ( = Ох(0,Т)), обладающих
При этом ю1(| ^ |) = 0 при | ^ |> 1; | ю^ \ |) = 1.
Предполагается, что рт ^ 0 при т ^да . Функции Ь(т) имеют все производные любого порядка [6]. Поскольку Ь е Ь (О) (q > 0), то для любой подобласти
О' с О , отстоящей от границы О на расстоянии, не меньшем рт, усреднение Ь(т) сходится к Ь в (О'),
т. е
. ||b(m)
- b < sup
"q,°' |v|<Pm
г q
I | b(x- v) -b(x) | dx
V G'
^ 0
при т ^да. Таким образом, аппроксимации Ь(т) сходятся к Ь в норме пространства Ь (О).
Известно также, что сходимость в Ь (О) влечет
сходимость по мере, а из сходящейся последовательности функций по мере можно выбрать подпоследовательность, сходящуюся почти всюду [12]. Причем в [12] дан конструктивный метод построения такой подпоследовательности. В связи со сказанным далее будем считать, что такая подпоследовательность
выбрана, и это есть последовательность Ь( т).
В результате замены в (1)-(6) коэффициента температуропроводности функцией Ь(т) приходим к следующей краевой задаче:
du(т) = b(m) (Л d2u(m) + Л db(m) du(m) dt dx2 dx, dx
i=i .,(m)
Как утверждается в [i0] (теорема 4.5, гл. III), обобщенное решение краевой задачи с разрывными коэффициентами в параболическом уравнении в норме пространства (QT) устойчиво относительно
вариаций коэффициентов. Из этой теоремы применительно к задаче (i)-(6) следует, что если при т ^го
равномерно ограниченная последовательность b(m) сходится почти всюду к b, тогда обобщенные решения задач с коэффициентами b(m) сходятся сильно в норме V20,1 (QT) к обобщенному решению задачи (i)-(6).
Таким образом, при достаточно больших значениях т приближенное решение u(m) задачи (i)-(6) с разрывным коэффициентом температуропроводности можно получить, если решать задачу, в которой в уравнении (i) этот коэффициент заменить на его
гладкую аппроксимацию b( m). В качестве такой заме-
ды(т >
dx
du(m)
(x,0) = ф( x)
= 0 du(m)
xi =-/i
dx2
= 0,
dx.
du(m)
= 0,
xi=/i
X
x2 = -/2 du(m)
dx2
= 0,
x2 = /2
dx3
= ai(t) (u - yi(t)),
x3 =0
-X
du(m)
dx3
= a2(t) (u - ^(t)).
(9)
(i0) (ii)
(i2)
(13)
(14)
Далее в работе будут использоваться следующие обозначения: дО - граница области О; хА - индикаторная функция множества А. Будем предполагать, что заданные в граничных условиях функции а1, а2,
0 G
Vj, у2 такие, что для некоторого 0 < а < 1 существует решение задачи (9)-(14) в пространстве Гельдера
3+а 3+а,-
H 2 (QT) [10] (теорема 5.3, гл. IV). Доказано [13], что при выполнении условий существования и единственности решений задач (1)-(6) и (9)-(14) для любого е > 0 существует натуральное число те такое, что при m > те выполняется неравенство
vrai sup | u(m) - u |< е . Следовательно, выбирая радиус
Qt
усреднения достаточно малым, можно с любой заданной точностью приблизить решение u(m) к u.
Вероятностное представление решения. Для решений краевых задач для параболических уравнений известны вероятностные представления [14; 15]. Чтобы получить вероятностное представление решения задачи (9)-(14), введём систему СДУ:
r
Xr = х + J ov x G ( Xv )dWv(3) +
T-t
r r
+ J av x g (Xv )dv + J n( Xv )x qg X )dKv, (15)
T-t T-t
r
Yr =1 + J a1(v)Yv Xx3 (v)=o dKv +
T-t
r
+ J a2(v)YvXx3(v)=z3dKv, (16)
T-t r
Zr = J XX3(v)=0dKv +
T-t
r
+ J а2(v)v2(v)YvXx3(v)=/3dKv, (17)
T-t
r
Kr =J Xqg (Xv )dKv, (18)
T-t
где х e G - начальная точка для Xv ; Wv(3) - трехмер-
ный винеровскии процесс; ctv =
(
вектор с координатами
Л
db(m) V dx, ,
где символом E,
x,T-t
обозначено математическое
ожидание относительно вероятностной меры Рх7 —,
соответствующей случайному процессу, стартующему в момент времени 7— из точки х, т. е. Х7= х.
Приближенные значения и (т)(х, г) можно получить методом Монте-Карло путем статистического моделирования траекторий системы СДУ (15)-(18). Для этой цели в работах [6; 13] предлагалось использовать метод Эйлера. Моделирование траекторий методом Эйлера для СДУ (15)-(18) осуществляется по следующей схеме:
1
Хгд+1 = х + И 2 о г £ г + Иаг, (20)
хг+: = Хгд+1 +(Аг+1 К)1кзу& (Хгд+1), (21)
У+1 = У (1 + (а1 )1 хдо(х )Д1+1К +
+ («2)1 Хдо(х )А1+!К), (22) X+1 = X + У (а1)1 О1 )1 Xдо(х1 )А1+1К +
+У (а2 )1 (^2 )1 XдО (х )А1+1К, (23)
к+! = к, +А1+1К, (24)
где 1, 1 + 1 - номера узлов сетки; - вектор трех независимых М(0,1) случайных величин;
Дг+1 K = d ( хД+, )
(2b(m)(Xv ))); av
i = 1, 2, 3; n(P) -
внутренняя нормаль в точке P е dG ; Kv - локальное время пребывания процесса Xv на границе, т. е. скалярный возрастающий процесс, который растет, только когда Kv edG.
Вероятностным представлением решения задачи (9)-(14) является математическое ожидание
и (x, t) = ExT_t (ф(XT )YT + ZT), (19)
расстояние от ХД до dG в случае выхода ХД+1 из области.
Применение методов Эйлера и случайного блуждания по движущимся сферам. Поскольку в работе предполагается, что теплофизические своИства материалов, из которых состоит сотовая панель, являются постоянными величинами, то случайный процесс Xv
в однородной среде, соответствующей каркасу или ячейке, с точностью до постоянного множителя совпадает с винеровским процессом.
Известно, что точка первого выхода из шара вине-ровского процесса, начинающего движение из центра этого шара, имеет на границе этого шара равномерное распределение. При этом эта точка не зависит от времени первого выхода на границу. В связи со сказанным предлагается для ускорения моделирования траекторий СДУ (15)-(18) при их прохождении внутри ячейки использовать метод случайного блуждания по сферам, а движение по каркасу и некоторой его окрестности моделировать с помощью метода Эйлера с малым шагом.
При решении нестационарных задач методом блуждания по сферам требуется каждый раз моделировать не только точку первого выхода из шара, но и время первого выхода. Для моделирования времени первого выхода из шара радиуса R в работе используется определение времени первого достижения уровня R бесселевским процессом, который является решением СДУ [8]:
St = So +
j S-1dr + Wt.
(25)
Одномерный случайный процесс (19) получается в результате применения формулы Ито к евклидовой норме трехмерного винеровского процесса с начальной точкой x e G. При этом выполняется равенство
Время первого выхода трехмерного винеровского процесса из шара радиуса Я и время первого достижения случайным процессом уровня Я имеют одинаковое распределение. Известны явные формулы для плотности времени первого достижения бесселевским процессом заданного уровня для начальных точек = 0 и 50 > 0 [9]. Но эти формулы представлены в виде рядов, содержащих функции Бесселя первого рода и их положительные нули, и поэтому численно-статистическое моделирование по ним крайне затруднительно. С другой стороны, вместо заданного уровня, который требуется достичь, можно подобрать некоторые кривые, зависящие от t, такие, что получаются достаточно простые и удобные для моделирования формулы для плотности времени первого достижения бесселевским процессом этих кривых. Каждая такая кривая определяется задаваемой положительной сигма-конечной мерой д. Определим функцию
1 М
г (V, х) = q(t,0, х)--Г q(t, у, х)д(ёу), (26)
а
0
где q(t, у, х) - плотность вероятности перехода бессе-левского процесса; а - числовой параметр.
Из свойства плотности q следует, что функции q и г удовлетворяют параболическому уравнению
ду = 1__д( 1 дt 2 дх2 дх I х
(27)
4
р* ^) = _ 4
(* (t)
| г^, х)йх
(28)
Если положить д(ф) = у 4у, то и р* (t) определяются по формулам [9]
* (t) =
2t 1п
V V
И)
1
^ 2
//
Р* С) =
2at
(29)
т* =
'((>2 И)
(30)
Пусть функция g(t) такая, что на кривой х = g(t) выполняется равенство г(t,x) = 0. Тогда эту кривую можно рассматривать в качестве граничного условия к уравнению (27). Причем с учетом того, что бессе-левский процесс стартует из заданной точки, решение полученной краевой задачи определяется единственным образом.
Обозначим т* время достижения бесселевским
процессом этой подвижной границы, т. е. т* = = М {/ > 0, = *(/)}. Можно показать [9], что плотность времени достижения границы х _ * ^) = 0 удовлетворяет уравнению
В этом случае получается простая формула для моделирования т*
где Н - гамма-распределенная случайная величина
5 2 В й
с параметрами —, 3. В ходе вычислений, после того,
как получено значение т* , радиус сферы для определения очередной точки положения случайного процесса определяется по первой формуле в (29). Выбор параметра а дает возможность определить радиусы моделируемых шаров так, чтобы они не выходили за границы текущей ячейки и чтобы они не были слишком малы. В случае, если расстояние до границы ячейки менее некоторой заданной величины е, вычисления производятся согласно методу Эйлера.
Результаты расчётов. Вычисления проводились для сотовой панели, каркас которой и ограничивающие пластины изготовлены из углепластика, а наполнителем служит воздух. Основные характеристики панели следующие: общая толщина панели - 0,035 м; толщина каждой из ограничивающих пластин панели -0,001 м; толщина стенок сотового каркаса - 6 -10_5 м; сотовыми ячейками являются правильные шестиугольники одинакового размера с длиной стороны, равной 0,0042 м; коэффициенты температуропроводности углепластика и воздуха равны соответственно 8,0 -10 _4 и 2,36 -10_5 м2/с.
Граничные условия для расчётов были сформированы на основании физических характеристик, полученных по данным лётного эксперимента в холодном климате и соответствуют первым 150 с взлёта самолета. На рис. 1 приведён график коэффициента теплообмена на внешней стороне обшивки.
На рис. 2 приведены полученные в результате расчётов при объёме выборки 2000 траекторий случайного процесса графики изменения значений температуры в середине панели и у края панели, обращённого внутрь салона.
В табл. 1, 2 приведены для сравнения затраты времени счёта предложенным комбинированным методом и только методом Эйлера с различными шагами по каркасу к^ и по объёму внутри ячеек ка.
Результаты расчётов позволяют сделать следующие выводы:
а) предложенный метод значительно эффективнее использовавшегося ранее метода Эйлера;
б) так как при увеличении шага в методе Эйлера эффективность комбинированного метода растет по сравнению с методом Эйлера, то основная часть вычислительных затрат приходится на движение по тонкому каркасу;
в) затраты в средней точке панели меньше, чем у края, поскольку в середине панели большая часть траектории приходится на движение по внутреннему пространству ячеек, где вычисления проводились с меньшим шагом.
2
Рис. 1. Коэффициент теплообмена на внешней стороне обшивки Fig. 1. Heat transfer coefficient on the outside of the skin
Рис. 2. Значения температуры внутри панели Fig. 2. Temperature values inside the panel
Затраты времени счёта в средней точке панели, с
Таблица 1
Метод ha = 2,5 • 10~7, hf = 2,5 • 10~8 ha = 5,0 • 10~7, hf = 5,0 -10~8 ha = 10~6, hf = 10~7
Комбинированный метод 1064 548 154
Метод Эйлера 5839 3070 1563
Таблица 2
Затраты времени счёта вблизи границы х3 = 0, c
Метод ha = 2,5 -10-7, hf = 2,5 •Ю-8 ha = 5,0 • 10-7, hf = 5,0 -10~8 ha = 10~6, hf = 10-7
Комбинированный метод 1394 479 188
Метод Эйлера 6193 4156 1854
Заключение. В работе предложен численно-статистический метод для оценки теплового состояния теплозащитных конструкций сотового типа. Данный метод основан на вероятностном представлении краевой задачи для параболического уравнения и комбинированном применении методов Эйлера и случайного блуждания по движущимся сферам. Получено значительное ускорение счёта по сравнению с моделированием траекторий только одним методом Эйлера.
Благодарности. Работа С. А. Гусева выполнена при финансовой поддержке РФФИ, грант № 17-01-00698 А.
Acknowledgments. S. A. Gusev's work was carried out with the financial support of the Russian Foundation for Basic Research, grant No. 17-01-00698 A.
Библиографические ссылки
1. Результаты экспериментальной отработки спускаемой капсулы космического аппарата «Фобос-Грунт» для доставки образцов грунта Фобоса на Землю / С. Н. Алексашкин [и др.] // Вестник ФГУП «НПО им. С. А. Лавочкина». 2011. № 5. C. 3-10.
2. Николаев В. Н., Гусев С. А. Экспериментально-теоретические исследования негерметизированного отсека летательного аппарата с сотовыми конструкциями обшивки // Авиакосмическое приборостроение. 2017. № 9. C. 10-19.
3. Николаев В. Н., Гусев С. А. Математическое моделирование теплового состояния отсеков и систем самолета с сотовыми конструкциями обшивки // Полёт: Общероссийский научно-технический журнал. 2017. № 5. С. 3-11.
4. Nikolayev V., Gusev S. Heat condition compartments of aircraft with a honeycomb structure. Saarbrucken : LAMBERT Publ., 2017. 133 p.
5. Николаев В. Н., Гусев С. А. Исследование теплового состояния отсеков пассажирского самолёта с сотовыми конструкциями фюзеляжа // Научный вестник Новосибирского государственного технического университета. 2016. № 1 (62). C. 146-167.
6. Gusev S. A., Nikolaev V. N. Calculation of heat transfer in heterogeneous structures such as honeycomb by using numerical solution of stochastic differential equations // Advanced Materials Research. 2014. Vol. 1016. Р. 758-763.
7. Muller M. E. Some continuous Monte Carlo methods for the Dirichlet problem // The Annals of Mathematical Statistics. 1956. Vol. 27, No. 3. Р. 569-589.
8. Revuz D., Yor M. Continuous Martingales and Brownian Motion. 3 ed. Springer Publ., 1999. 602 c.
9. Deaconu M., Herrmann S. Hitting time for the Bes-sel processes - walk on moving spheres algorithm (WOMS) // The Annals of Applied Probability. 2013. Vol. 23, No. 6. Р. 2259-2289.
10. Ладыженская О. А., Солонников В. А., Ураль-цева Н. Н. Линейные и квазилинейные уравнения параболического типа. М. : Наука, 1967. 736 с.
11. Соболев С. Л. Некоторые применения функционального анализа в математической физике. 3-е изд. М. : Наука, 1988. 336 p.
12. Колмогоров А. Н., Фомин С. В. Элементы теории функций и функционального анализа. М. : Наука, 1976. 336 с.
13. Гусев С. А. Применение СДУ к оценке решения уравнения теплопроводности с разрывными коэффициентами // Сиб. журн. вычисл. математики / РАН. Сиб. отд-ние. 2015. Т. 18, № 2. C. 147-161.
14. Гихман И. И., Скороход А. В. Стохастические дифференциальные уравнения и их приложения. Киев : Наукова думка, 1982. 612 с.
15. Кушнер Г. Дж. Вероятностные методы аппроксимации в стохастических задачах управления и теории эллиптических уравнений. М. : Наука, 1985. 222 с.
References
1. Alexashkin S. N., Martynov M. B., Pichkhadze K. M., Finchenko V. S. [Results of integrated development testing of "Phobos-Grunt" SC's descent capsule, intended for Phobos soil samples delivery to the Earth]. Vestnik FGUP NPO im. S. A. Lavochkina. 2011, No. 5, P. 3-10 (In Russ.).
2. Nikolaev V. N., Gusev S. A. [Experimental-theoretical studies of the unsealed aircraft compartment with honeycomb structures]. Aviakosmicheskoe priboros-troenie, 2017, No. 9, P. 10-19 (In Russ.).
3. Nikolaev V. N., Gusev S. A. [Mathematical modeling of the thermal state of compartments and aircraft systems with honeycomb structures]. Polet. Obshcherossiy-skiy nauchno-tekhnicheskiy zhurnal. 2017, No. 5, P. 3-11 (In Russ.).
4. Nikolayev V., Gusev S. Heat condition compartments of aircraft with a honeycomb structure. Saarbrucken, LAMBERT Publ., 2017. 133 p.
5. Nikolaev V. N., Gusev S. A. [Simulation of heat condition compartments of passenger aircraft with a honeycomb structure of the fuselage. Scientific]. Nauchnyy vestnik Novosibirskogo gosudarstvennogo tekhnicheskogo universiteta. 2016, No. 1 (62), P. 146-167 (In Russ.).
6. Gusev S. A., Nikolaev V. N. Calculation of heat transfer in heterogeneous structures such as honeycomb by using numerical solution of stochastic differential equations. Advanced Materials Research, 2014, Vol. 1016, P. 758-763.
7. Muller M. E. Some continuous Monte Carlo methods for the Dirichlet problem. The Annals of Mathematical Statistics, 1956, Vol. 27, No. 3, P. 569-589.
8. Revuz D., Yor M. Continuous Martingales and Brownian Motion. 3ed., Springer Publ., 1999. 602 p.
9. Deaconu M., Herrmann S. Hitting time for the Bessel processes - walk on moving spheres algorithm (WOMS).
The Annals of Applied Probability, 2013, Vol. 23, No. 6, P. 2259-2289.
10. Ladyzhenskaya O. A., Solonnikov V. A., Uralt-seva N. N. Lineynye i kvazilineynye uravneniya parabo-licheskogo tipa [Linear and quasilinear equations of parabolic type]. Moscow, Nauka Publ., 1967, 736 p.
11. Sobolev S. L. Nekotorye primeneniya funktsio-nal'nogo analiza v matematicheskoy fizike [Applications of Functional Analysis in Mathematical Physics]. Moscow, Nauka Publ., 1988.
12. Kolmogorov A. N., Fomin S. V. Elementy teorii funktsiy i funktsional'nogo analiza [Elements of the the-
ory of functions and functional analysis]. Moscow, Nauka Publ., 1988. 336 p.
13. Gusev S. A. Numerical Analysis and Applications. Springer Publ., 2015, No. 8(2), P. 22-34.
14. Gikhman I. I., Skorohod A. V. Stokhasticheskie differentsial'nye uravneniya i ikh prilozheniya [Sochastic differential equations and their applications]. Kiev, Naukova dumka Publ., 1982, 612 p.
15. Kushner H. J. Probability Methods for Approximations in Stochastic Control and for Elliptic Equations. Academic Press, 1977, 222 p.
© Гусев С. А., Николаев В. Н., 2017
