CC BY
30
УДК 51-73
РАЗРЕШИМОСТЬ КРАЕВЫХ ЗАДАЧ, ОПИСЫВАЮЩИХ ДИФФУЗИЮ АТОМОВ ПЛЕНКИ В ПОДСТИЛАЮЩЕЙ ПОВЕРХНОСТИ ПРИ ОБРАЗОВАНИИ ТОНКОПЛЕНОЧНЫХ СТРУКТУР
Тарасенко Елена Олеговна,
канд. физ.-мат. наук, доцент кафедры прикладной математики и математического моделирования Института естественных наук и математики Северо-Кавказского Федерального университета, Россия, 355000, г. Ставрополь, ул. Пушкина, 1, корп. 2. E-mail: galail@mail.ru
Гладков Андрей Владимирович,
ст. преподаватель кафедры прикладной математики и математического моделирования Института естественных наук и математики Северо-Кавказского Федерального университета, Россия, 355000, г. Ставрополь, ул. Пушкина, 1, корп. 2. E-mail: gavandrew@mail.ru
Маликова Наталья Владимировна,
учитель математики и информатики средней общеобразовательной школы № 1, Россия, 356250, Ставропольский край, с. Грачевка, ул. Советская, 47. E-mail: nmalik_@mail.ru
Актуальность работы. Добыча георесурсов требует создания новых технологичных решений их производства, например, покрытие буровых шнеков средствами для защиты от коррозии. Решение указанной задачи возможно при использовании диффузионного нанесения тонкой пленки (антикоррозийное вещество) на подстилающую поверхность (шнек). Математизация такого физического процесса, как диффузионный рост тонких пленок на подстилающей поверхности, в настоящее время является малоисследованной. При математическом моделировании часто возникает вопрос о существовании и единственности решения краевых задач, описывающих указанный физический процесс. Многие отечественные и зарубежные ученые проводили исследования по решению аналитическими и численными методами начально-граничных задач, в которых изначально явно или неявно допускается, что решение рассматриваемой задачи существует и единственно. Как правило, авторы публикаций, посвященных различным проблемам математического моделирования диффузии, либо вообще не затрагивают этот вопрос (о существовании и единственности решения), либо без должного на то основания ссылаются на классические работы. Поэтому исследования на разрешимость краевых задач, проводимые в данной работе, являются актуальными.
Цель исследования. Разработать критерии разрешимости (существование и единственность) краевых задач, возникающих при математическом моделировании роста тонкопленочных структур на подстилающей поверхности, в различных пространствах. Методы исследования. Достижение поставленной цели основывается на корректном использовании результатов и методов уравнений математической физики, интегральных уравнений, математического анализа, уравнений в частных производных, физики твёрдого тела, кристаллографии.
Результаты. Проведено исследование на разрешимость краевых задач, описывающих диффузионный рост тонких пленок на подложках; разработаны критерии существования и единственности решения указанных задач в различных пространствах. Выводы. В ходе проведения исследований при математическом моделировании диффузионного роста тонкой пленки на подстилающей поверхности были разработаны теоремы (критерии), обеспечивающие разрешимость (существование и единственность решения) начально-граничных задач. Рассмотрены краевые задачи для случаев полного отражения и поглощения атомов пленки подстилающей поверхностью. Настоящая статья представляют значительный интерес в прикладных исследованиях, позволяет ответить на вопрос: можно ли сразу приступать к численному (или, возможно, аналитическому) решению конкретно рассматриваемой краевой задачи, описывающей диффузионный рост тонкопленочных структур на подложках, или дополнительно проводить исследования по ее регуляризации.
Ключевые слова:
Разрешимость, краевая задача, тонкая пленка, подстилающая поверхность, подложка, диффузионный рост
Необходимость обработки больших массивов информации в реальном масштабе времени ставит задачу создания устройств функциональной электроники, объединяющих функции ввода, преобразования и вывода информации для последующей её обработки в цифровом коде с помощью традиционных принципов. Создание таких устройств упирается в изучение физического процесса - диффузии, возникающей в процессе роста тонких пленок на подложках. Исследованию такого процесса
посвящено значительное число работ [1-13], как отечественных, так и зарубежных исследователей. Задача диффузионного роста тонких пленок на подложках имеет широкое практическое применение в опто- и микроэлектронике, является прикладной задачей технологий георесурсов.
При математическом моделировании процесса диффузионного роста тонкой пленки на подстилающей поверхности возникает вопрос о разрешимости таких моделей.
В теории и практике современных исследований процесса диффузии, возникающей при образовании тонкопленочных структур, используются две начально-граничные задачи [3, 5, 14]. Одна из которых в самой общей постановке имеет вид:
дд 1 ^ дд 1 ^ ^ д ^ дд
^ + £ + aq — — К^ + f, ^ ^ ' дх, T1 j—'dx. " dx
Kj —
,=1 j—1 дх,. i, j = 1,2,3,
(1)
du.
У u —¡- = 0,
t! . dx,. '
(2)
t e[0,T], 0 < T <+<», (x, X, X) eG, q(t,x,X,X) - 0, t e[0,T), (Xj, x2, x3) eG, (3)
q(0, x, x, x) =ф( xi, x, x),
(4)
д(/, х1, х2, х3) |дс = х, х2, Х3). (5)
Вторая начально-граничная задача представляется в виде (1)-(4) и граничного условия
^ + ßq
dv
— F,
(6)
где q(t,x1,x2,x3) - функция, значения которой в каждый момент времени tе[3,Г) совпадает со средним значением концентрации атомов пленки в связной области О, дО=дО0^дО1идО2, дО0 - нижняя, дО! - боковая, дО2 - верхняя части границы дО, О-=ОидО; ц¡=ц¡(t,x1,x2,x3), ¿=1,2,3 - функции, значения которых совпадают со значениями средней скорости горизонтального переноса в момент t в точке (x1,x2,x3) соответственно вдоль осей Ox1, Ox1, Ox1 (рассматривается декартова прямоугольная система координат); a=a(t,x1,x2,x3) - функция, характеризующая взаимодействие атомов пленки с окружающей средой (или их радиоактивный распад) в момент t в точке (x1,x2,x3); K¡¡=K¡¡(t,xl,x2,x3), ¿,7=1,2,3 - элементы матрицы коэффициентов диффузии атомов пленки в подстилающей поверхности; f=f(t,xl,x2,x3) - функция, моделирующая источник атомов пленки (функция источника); ф=ф(t,x1,x2,x3) - функция, значения которой в точке (x1,x2,x3)eО в момент времени ^ совпадает со значениями концентрации атомов пленки (функция, описывающая фоновую концентрацию); ^ -
производная по внутренней нормали dG0: dq(t,x) = lim £ К (t, x)cos( N, x) dq(t,x)
dv
dv(t, x)
j=1
dv,
(7)
х = Х(Х1, ^ Х3), у = у(у^ у^ Уз), где N - внутренняя нормаль к дО0 в точке xeдО0•, К - конечный, замкнутый конус с вершиной xeдО0, который содержится в G+{x}; Дt,x1,x2,x3) -функция, характеризующая гравитационное осаждение атомов пленки на дО0; F(t,x1,x2,x3) - скорость осаждения атомов пленки на дО0, tе[3,Г),
Функция источника атомов пленки f задается в виде [14, 15]:
f (1, х, X, Х) = Q(t, Х\, X, Х Ж^, х1, х2, х3), (8)
где Q(t,x) - мощность источника атомов пленки (масса атомов пленки, выбрасываемых в области О в момент t в точке xeО), 5(t,x) - дельта-функция Дирака. При этом, если источник является - момент начала действия источника): _
1) точечным, сосредоточенным в точке x0=(x0,x20,x30)eО,
1.1) мгновенного действия, то ^^^^сопя^
f (/, х, X, Х) = в&(08(Х - х0)5(х - х°)8(х - х3°);
1.2) непрерывного действия, то
о а, х)=),
_/"(/, Х1, Х2, хз) =
= е(О^Жх - Х0Жх - х20)5(^ - х3°);
2) линейным, сосредоточенным на интервале [а,Ь] числовой прямой, параллельной оси Ox2 и пересекающей ось Ox3 в точке (0,0,x30),
2.1) мгновенного действия, то
_ , [0(x2), х2 ЬЬ
0(1, Х1, Х2, Х3) = 1 п „г ,,
[ 2, х £[а,Ь],
у(Х Х1,x2,хз) = е(хШШх1 - х10)5(хз - Х32);
2.2) непрерывного действия, то [0(1, х), ( е[2, Т], х е[а, Ь],
2, х £[а, Ь],
У(/,х,Х, Х) = 0(1,XХ - х°)5(х -х3°);
3) площадным, сосредоточенным на площадке 5, лежащей на плоскости x1Ox2 и пересекающей ось Ox3 в точке (0,0,x3),
3.1) мгновенного действия, то
[е(х,х), х1,х 2, х, Х £ ^,
у(/, Х1,x2, хз) = б(Х1,х2)5(/)5(хз - хз2);
3.2) непрерывного действия, то
^ , [°(1,Х1,x2), 1 е[2,Т] Х1,Х2 ^ ОМ, х,, х, Х) = 1
1 2, х, Х £
У(/, х, Х, Х) = 0(1, Х\, ХХ - х3°);
4) поверхностным, сосредоточенным на поверхно-
Q (t, xj, ^^^, x) —
Q(t, x1, x2, x3) —
Q(t, x1, x2, x3) —
сти Sn тела П,
4.1) мгновенного действия, то
|Q x2 , (X1, x2, x3) eSП,
0, (x1, x2, x3) ,
f (t, x, x2, x3) — Q( X, x2, x3)5(t);
4.2) непрерывного действия, то
Q (t, xj, ^^^, x) —
|Q(t,x,x,x), t e[0, T], (xj,x, x) eSn, 0, (x, x, x) ^ Sn,
f (t, x1, x2, x3) — Q(x1, x2, x3).
Уравнение (1) можно переписать в эквивалентном виде:
3 я„ 3 3 я2.
Ё1 + У И. * + ад = у У
д( £ ' дх. 4 и £ ' дх.дх.
И У дК . .
И.=и. -У -' = 1,2,3.
' ' ^ дх .=1 (Ж.
+/,
(9)
(10)
Следует заметить, что уравнение (9) (а значит и (1)) совпадает с уравнением из [16]:
Ьи = /ф, (11)
п п
Ьи = У У а
д и -А ди ди ' У"н--+ У Ь--+ си--
=1 .=1 дхрх '=1
дх,.
дг
(12)
при п=3 с точностью до знака у/ и /ф: / и /ф имеют противоположные знаки. Этот факт будет учитываться в приводимых ниже результатах исследования.
Уравнению (1) поставим в соответствие уравнение ^ + У и,^ +аЧ = у у — Кн^ + е, (13)
дг дх , дх, дх.
уравнению (9) - уравнение
д' +±*% + аЧ = У±К,] -дд- + е, (14)
дг
'=1 н= 1
] дх,дх..
отличающиеся, соответственно, от (1) и (9) лишь видом функции /: вместо /, задаваемой выражением (8), здесь рассматривается мощность источника атомов пленки е.
Рассмотрим следующую задачу: найти (указать) условия, при выполнении которых задачи (1)-(5) и (1)-(4), (6) имеют единственное решение.
Несмотря на очевидную необходимость проведения таких исследований (решению аналитическими и численными методами указанных начально-граничных задач посвящено значительное число работ, в которых изначально явно или неявно допускается, что решение рассматриваемой задачи существует и оно единственное), подобных исследований в этом направлении до настоящего времени не проводилось. Как правило, авторы публикаций, посвященных различным проблемам математического моделирования диффузии, либо вообще не затрагивают этот вопрос (о существовании и единственности решения), либо без должного на то основания ссылаются на классические работы [16-18]. Ниже можно будет убедиться: достаточно ясное и четкое освещение данного вопроса не является тривиальным и требует скрупулезного анализа результатов из [16, 17]. Исключение составляет монография [14], однако в этой работе найдены лишь достаточные условия единственности решения задач типа (1) ((9))-(5) и (1) ((9))-(4), (6). Вопрос о существовании их решения в [14] не затрагивался.
Теорема 1. Пусть коэффициенты щ, K¡j, а,
ра —т
i,j=1,2,3, принадлежат классу Н 2(Б4) и ограничены на Dт^, кроме того, щ, Kjj непрерывно диффе-
ренцируемы по xi, i,j=l,2,з в Dт, е>о, ф>о, е, ф ограничены, е удовлетворяет условию Гёльдера с показателем р, ф непрерывна в Dl, удовлетворяет условиям Ляпунова. Тогда задача (1)-(5) при у=0 имеет единственное решение, совпадающее с решением задачи (13), (2)-(5).
Доказательство. Так как щ, K¡j, -,;'=1,2,3, непрерывно дифференцируемы по xi в DT, то уравнение (1) эквивалентно уравнению (9), уравнение
(13) - уравнению (14).
Рассмотрим задачу (13), (2)-(5), которая эквивалентна задаче (14), (2)-(5). При выполнении условий данной теоремы, очевидно, выполняются условия:
• теоремы 16.2 из гл. 4 § 16 [17], а значит, решение задачи (14) ((13)), (2), (4), (5) существует;
• теоремы 5.2 из гл. 4 § 5 [17], а значит, решение задачи (14) ((13)), (2), (4), (5) единственное и оно представимо в виде
д(г, х1, х2, х3) =
г
= I¿тШ ^ ^ х2, х3; ^ у^ У2, у^у&Ж 3 +
0 о
+ Ш Чо , X1, X2, х3; Т, Уl, У2, Уъ)ф(Уl, У2, У3^У^3 (15)
о
(см. соотношение (16.17) из [17]), где q0(t,xl,x2,xз;т,yl,y2,yз) - функция -Грина для задачи
(14) ((13)), (2), (4), (5) в области DT^, т. е. ^ удовлетворяет уравнениям
± да* ^ ^ +
дг '=1 ' дх , , =1 ] =1 " дх, дх]
+5(г - т)5(х - У1)5(х2 - У2)5(х3 - У3),
( €%_ дг
У
да0 и ' + дх..
аа0 =у у ± +
]=1 дх1 дх
Л
-т)8(х1 - У1)5(х2 - У2)5(х3 - У3) )
и условиям:
ао(t, xl, х2, x3;T, Уl, У2, У3) = °
(16)
(17)
(18) (19)
а о (и xl, х2, х3; т у l, У2, У3) |( х ^ ^) Едо = 0.
Кроме того, в условиях данной теоремы выполняются условия теоремы 2.1 и следствия 2.1 из гл. 1 § 2 [17]. А тогда
Ч (г, х, х2, х) > 0, (г, х, х, хъ) еОт. (20)
Из (20) и условий е>0, ф>0 следует, что если выполнены условия данной теоремы, то решение задачи (14) ((13)), (2)-(5) существует и единственно.
Учитывая равенство (8) и используя свойства 5-функции Дирака [19], соотношение (15) можно переписать в эквивалентном виде:
а(г, х, х2, х) =
г
= Iа0(t, X1, х2, XЪ;T, Уl, У2, У3) /^У1^У2^У3 +
о о
+Шао(г,х1,х2,х3;т,У1,У2,У3)ф(У1,У2,У3)^1^2¿Уъ. (21)
о
Снова воспользовавшись свойствами 5-функ-ции, непосредственным подсчетом, можно убедиться, что функция (21) удовлетворяет уравнению (9), а следовательно, и уравнению (1). Учитывая, что q0(t,x1,x2,x3•т,y1,y2,y3) является функцией Грина для задачи (13), (2)-(5) (т. е. решением задачи (17)-(19)), заключаем, что функция (21) удовлетворяет условиям (3), (4), (5). Значит, решение (1)-(5) существует и единственно. Теорема доказана.
Теорема 2. Пусть выполнены все условия теоремы 1. Тогда задача (1)-(5) при f=0, ф=0 имеет единственное решение.
Доказательство. Если выполнены условия данной теоремы, то:
1) выполняются все условия теоремы 16.1 из гл. 4 § 16 [17], а значит, решение задачи (9), (2), (4), (5) приf=0, ф=0 существует;
2) выполняются все условия теоремы 5.2 из гл. 4 § 5 [17], а значит, решение задачи (9), (2), (4), (5) приf=0, ф=0 единственно;
3) выполняются все условия теоремы 2.1 из гл. 1 § 2 [17], а значит, согласно следствию 2.1 из этой теоремы, решение задачи (9), (2), (4), (5) приf=0, ф=0 неотрицательно, т. е. выполняется условие (3).
Уравнение (9) эквивалентно уравнению (1). А тогда, согласно 1)-3), решение задачи (1)-(5) при f=0, ф=0 существует и единственно. Теорема доказана.
Теорема 3. Пусть выполнены все условия теоремы 1. Тогда (1)-(5) имеет единственное решение, совпадающее с решением (13), (2)-(5).
Доказательство. Все условия теоремы 3 те же, что и условия теорем 1 и 2. Обозначим через q1(t,x1,x2,x3) решение задачи (1)-(5) при у=0, через q2(t,x1,x2,x3) - решение этой задачи при f=0, ф=0. Непосредственным подсчетом легко убедиться, что q(t,x 1^x3) q1(t,x1 ,x2,xз)+q2(t,x 1 ,X2,Xз)является решением (1)-(5). Согласно теореме 1, q1(t,x1,x2,x3) -единственное решение задачи (1)-(5) при у=0. Согласно теореме 2, q1(t,x1,x2,x3) - единственное решение задачи (1)-(5) при f=0, ф=0. Очевидно, что q(t,x1,x2,x3) будет единственным решением задачи (1)-(5). Теорема доказана.
Перейдем к анализу задачи (1)-(4), (6). Заметим, что уравнение (1) совпадает с уравнением из [16]:
Lq =
,=i j=i
д q , dq dq
a--Xh--aq--,
дх.дх. dx dt
I (t, x) - K (t0, x0)|< A(|x - x0|+|t -10|),
| h,.(t, x) - ht (t, x0)|< A | x - x0f,
| a(t, x) -a(t, x°)|< A | x - x0 f,
i, j = 1,2,3, A = const > 0.
Теорема 4. Пусть выполнены условия а), б), и, К непрерывно дифференцируемы по_Х;, ¿=1,2,3 в DT, граница dGeC1+a, ф непрерывна в G и равна нулю в некоторой G - окрестности границы dG, f, v непрерывны на dG0x[0,r]. Тогда решение задачи
(1)-(4), (6) существует и единственное. Доказательство. По условию, щ, Кф ¿=1,2,3
непрерывно дифференцируемы по xt в DJ. Поэтому уравнение (1) эквивалентно уравнению (9), уравнение (13) - уравнению (14).
Аналогично тому, как мы это делали при доказательстве теоремы 1, рассмотрим задачу (13), (2)-(4), (6), эквивалентную задаче (14),
(2)-(4), (6).
Если выполнены условия данной теоремы, то, очевидно, выполняются условия теоремы 2 из гл. 5 § 3 [16] (случай n=3). Откуда следует, что решение q задачи (14) ((13)), (2), (4), (6) существует, единственно и представимо в виде
q(t, x1, x2, x3) =
t
= i Г(t,x1,x2,x3; т,4, ц,в)г(т, 4, ц,в)dS +
0 dO
+Ш r(t, x, x2, x ;т,4, Ц,в)ф(т,4,Ц,в) d^d^dd +
о
t
+JdT JJJ F(t,x, x, x; т, 4, Я,4, Ц,в)d^d^dd, (22)
0о
где г(т,4,ц,в) непрерывная на dGx[0,T] функция, являющаяся решением интегрального уравнения (22), которое представимо в виде (3.10) из гл. 5 § 3 [15], dS - элемент поверхности dG, Г - фундаментальное решение уравнения Lq=0,
тг, д(.) д(.) + L(.) = + L h+ a(.) -
dt "j dx,.
-XX K
i=i a=i
д2(.)
дх,. дх,.
= 0,
fL(.) = f + 1 f- + a() -
дt дх..
определенное в Предположим, что Ь удовлетворяет следующим условиям: -
а) Ь - параболический коэффициент в Б44, т. е. при всех ^^еЛТ и любого вещественного вектора е=(е1,е2,е3)^(0,0,0) выполняется условие
п п
С,х)е,е, >0;
,=1 ] =1
б) коэффициент Ь непрерывен в Б4 и для всех (t,x)eD¡, (О^еДТ и некоторого ^=сопя1;, 0<^<1, удовлетворяет условиям:
1=11=1 дх i
Kj =о
a дха
(23)
(24)
Обозначим в (7) (а значит, и в (6))
3
Ь (/, х) = £ К (/, х)соб( N, х).
]=1
В условиях нашей теоремы выполняются (для случая п=3) все условия теоремы 2.2 из гл. 1 § 2 [17] (принцип максимума). А тогда, согласно этой теореме, q(t,x) удовлетворяет неравенству
q(t1, х) > sup min
X>a0
Щ-0)
n • ve 0, min
3-D1
eXtmin u( х,0);
b
X- a.
-min(geX(t1-to))
при
(х, х, х) ^ 3G, (4, г), в) е 3G,
(25)
|ь| ^ Ь2+ь22 + Ь1,
а0 = тах(-1| К(г, х) ||),
х = (х, х, х), Ч > Ч > 0. Рассмотрим выражение
а(г, х1, х2, х3) =
г
= I ^ ^ (г, х, х2, х ;Т, 4, ц,в) Г(Т, 4, ц,в) +
0 эо
Г(/, X, х2, х3; т, ц, в)ф(т, ц, в) в +
о
г
+1(т||| Г(/, X,х2, х3;т, ц,в)/(т, ц,в)С&цСв, (26)
0о
полученное из (22) путем замены в последнем слагаемом (22) е на/. Из (8) и свойств 5-функции Дирака следует, что (26) эквивалентно (22) в DT, т. е.
а(г, х,х, х) = а(г, х,X, х), (г, х^ х, х) £В4т. (27)
Из (27) вытекает, что q(t,x1,x2,x3) удовлетворяет условиям (3), (4), (6), так как этим условиям в Dт^ удовлетворяет q(t,x1,x2,x3).
Подставим формально (26) в (1), т. е. вычислим Lq(t,x1,x2,x3), (t,x1,x2,x3)еDT, (t,x1,x2,x3)<íдG, где L(.) имеет вид (24). Согласно (26), (24),
Ьа(г, х, х, х) =
г
= Ь\(т|||г (г, х, х, х; т, 4, ц, в)г( т, ц, +
0 до
+Ь ||| Г(г, х1, х2, х3;т,£,ц,в)ф(т,£,ц,в) ¡£с(ц((в +
о
г
+Ь1 ¡тЩ Г(г, х, х, х; т, 4, ц, в)/(т, £ ц, в)¿£сСцс(в. (28)
0о
Вычислим выражение в правой части (28), упростив для этого каждое слагаемое:
г
Ь | (тШ Г(/, х, х, х; т, 4, ц, в)г(т, £,, ц, в) ¡8 =
0 до
||| Г(/, х, х, х; т, 4, ц,в) г(т, ц, в) ¡8 +
до
г
I(т||| [ЬГ(/, х, х, х; т, 4, ц,в)]г(т, ц, в)¡8 = 0, (29)
0 до
так как по условию Г(t,x1,x2,x3;т,£,ц,в) - фундаментальное решение уравнения L)=0, а согласно определению Г
Г(г, х, х, х;т, 4, ц, в) = 0
Lr(t, х, х, х ;т,4,),в) = о. (30) LIII r(t, х, X, X;т, 4, Л,в)ф(т, 4, ),в)d%d)de =
G
Ш([Lr(t,х,,х2, х3; т,),в)] х\ | ^ 2' 3' ¡=0, (31)
LfdTfff(Г(t,Xl,X2,хз';Т,4,г),в) xl =
= III r(t, х, х, х; т, 4,), в) f (т, 4,), в) d4d)de +
G
+fdT i i if[LГ(t, х1, X2, 4,),в)] xl =
+1 TG vxf(T,4,),B)d4d)dO
= f (t, x, х2, х3). (32)
Равенство (32) вытекает из равенств
T(t, х1, х2, х3; т, 4,), в) = 8(х - 4)5(y - y - в), III f (т,4, п,в)8( х-4)8( y - П)8( y-в) d4d^в =
G
= f (t, х1, х2, х3)
и равенства (30).
Из (28), (29), (31), (32) следует, что
Lq(t, х, х, х) = f (t, х, х, х). Из (33) заключаем, что g(t,xj,x2,x3) является решением уравнения (1), из (27) - что решение единственно и оно удовлетворяет условиям (3), (4), (6) (в силу того, что этим условиям удовлетворяет q(t,xj,x2,x3)). Теорема доказана.
Полученные результаты принимаются для анализа математических моделей, используемых на практике.
Убедимся, что для основных начально-граничных задач, используемых в прикладных исследованиях диффузионного роста тонкопленочных структур на подложках, выполняются все условия теорем 1-4. В этих задачах обычно полагают, что коэффициенты и функции в задачах (1)-(5), (1)-(4), (6) имеют следующий вид [20] (для полного соответствия с обозначениями, используемыми на практике, будем считать, что x=x, x2=y, x3=z):
U = u (z) = С lnС = const > 0, u2 = щ = 0 (33)
(этот случай означает, что ось Ox сориентирована по направлению вектора скорости горизонтального переноса, а скорость ветра вдоль оси Oz изменяется по логарифмическому закону);
f Kf Ф 0, i = j, K„ =\ J i, j = 1,2,3
j [ K j = 0, i* j,
т. е. в матрице коэффициентов диффузии учитываются только диагональные элементы, а все элемен-
ты, не расположенные на главной диагонали, считают равными нулю; при этом К11=К22=с2и1, с2=сопя!;>0, где и1 задается выражением (34), К33=с32+с4, с3=сопя1>0, с4=сопя1>0.
В качестве О часто выбирают прямой круговой цилиндр высотой Н с достаточно большим радиусом Я основания, расположенного на подстилающей поверхности 2=0 [14]. Такой способ задания О удобен при аналитических (если это возможно в отдельных случаях [15]) и численных решениях рассматриваемых начально-граничных задач.
Функции а, Q, ф, у, ¡, V задают таким образом, что выполняются условия теорем 1-4. Чаще всего полагают, что эти функции являются постоянными величинами в О.
При данном выборе и, ¿=1,2,3, условие (2) выполняется тождественно.
Для указанных и, а, Кф ф, у, ¡, V, О условия теорем 1-4 выполняются. Поэтому используемые в прикладных исследованиях задачи вида (1)(5) и (1)-(4), (6) всегда имеют (и при том одно) решение.
Обратим внимание на следующее обстоятельство. При численном решении задач (1)-(5) и (1)-(4), (6) в уравнении (1) часто f заменяют на Q
СПИСОК ЛИТЕРАТУРЫ
1. Болтакс Б.И. Диффузия и точечные дефекты в полупроводниках. - Л.: Мир, 1972. - 362 с.
2. Верлань А.Ф., Москалюк С.С. Математическое моделирование непрерывных динамических систем. - Киев: Наукова думка, 1988. - 362 с.
3. Галай Е.О. Математическая модель образования плёнок на подложках // Обозрение прикладной и промышленной математики. - 2005. - Т. 12. - Вып. 4. - С. 932.
4. Самарский А.А., Гулин А.В. Численные методы математической физики. - М.: Научный мир, 2000. - 358 с.
5. Семенчин Е.А., Тарасенко Е.О. Моделирование процесса диффузии в полупроводниках // Обозрение прикладной и промышленной математики. - 2007. - Т. 14. - Вып. 2. - С. 352-353.
6. Тарасенко Е.О., Гладков А.В. Аналитические и численные решения некоторых обратных задач в рамках математической модели роста тонких пленок на подложках, уравнение которой допускает решение гауссового вида // Параллельная компьютерная алгебра и её приложения в новых инфокоммуника-ционных системах: сборник научных трудов I Международной конференции. - Ставрополь, 2014. - С. 93-98.
7. Тарасенко Е.О., Гладков А.В. Аналитическое решение краевой задачи о восстановлении мощности источника в математической модели роста тонких пленок на подложках // Параллельная компьютерная алгебра и её приложения в новых инфоком-муникационных системах: сборник научных трудов I Международной конференции. - Ставрополь, 2014. - С. 99-103.
8. Тарасевич Ю.Ю. Математическое и компьютерное моделирование. - М.: Едиториал УРСС, 2002. - 144 с.
9. Яценко Ю.П. Интегральные модели систем с управляемой памятью. - Киев: Наукова думка, 1991. - 284 с.
10. Aubry S. Toward a Rigorous Molecular Theory of Metastability // Ferroelectrics. - 1980. - V. 24. - P. 53-59.
без должного на то обоснования. Однако результаты численных расчетов в этом случае удовлетворительно согласуются с экспериментальными данными. Объяснить этот факт можно следующим образом. Из доказательств теорем 1-4 следует, что вид q(i,xj,x2,x3) не зависит от выбора в уравнении (1) в качестве свободного члена f или Q (см. (15) и (21), (22) и (26), (27)). Поэтому и результаты численных расчетов (при замене f на Q в уравнении (1)) всегда будут хорошо согласованы с экспериментальными данными (если только, конечно, сама модель (1)-(5) или (1)-(4), (6) адекватно экспериментальным данным описывает изменения значений q(t,x1,x2,x3) в G).
Заключение
В работе представлены результаты исследований на существование и единственность решения краевых задач, описывающих диффузионный рост тонкой пленки на подстилающей поверхности. Найдены критерии, обеспечивающие разрешимость (существование и единственность решения) в соответствующих пространствах. Данные критерии имеют важное значение в прикладных исследованиях.
11. Metiu H., Kitahara K., Ross J. Statistical mechanical theory of the kinetics of phase transitions // Fluctuation phenomena. Studies in statistical mechanics / Eds. Montroll E.W., Lebowitz J.L. North-Holland Publish. Comp., Amsterdam. - 1979. - V. 7. -P. 229-291.
12. Nishiama H. Martensitic Theory Transformations. - New-York: Academic, 1981. - 315 p.
13. Robledo A. The liquid-solid transition of the hard sphere system from uniformity of the chemical potential // J.Chem.Phys. -1980.- V. 72. - P. 1701-1712.
14. Марчук Г.И. Математическое моделирование в проблеме окружающей среды. - М.: Наука, 1982. - 320 с.
15. Семенчин Е.А. Аналитические решения краевых задач в математической модели атмосферной диффузии. - Ставрополь: СКИУУ, 1993. - 141 с.
16. Фридман А. Уравнения с частными производными параболического типа. - М.: Мир, 1968. - 427 с.
17. Ладыженская О.А., Уральцева Н.Н. Линейные и квазилинейные уравнения эллиптического типа. - М.: Наука, 1973. -576 с.
18. Ладыженская О.А. Краевые задачи математической физики. М.: Наука, 1973. - 406 с.
19. Афанасьев В.Н., Колмановский В.Б., Носов В.Р. Математическая теория конструирования систем управления. - М.: Высшая школа, 2003. - 614 с.
20. Берлянд М.Е. Современные проблемы атмосферной диффузии и загрязнения атмосферы. - Л.: Гидрометеоиздат, 1975. -448 с.
Поступила: 09.12.2015.
UDC 51-73
RESOLVABILITY OF BOUNDARY PROBLEMS DESCRIBING FILM ATOM DIFFUSION IN UNDERLYING SURFACE AT FORMATION OF THIN-FILM STRUCTURES
Elena O. Tarasenko,
North Caucasian Federal University, bld. 2, 1, Pushkin Street, 355000, Stavropol, Russia. E-mail: galail@mail.ru.
Andrey V. Gladkov,
North Caucasian Federal University, bld. 2, 1, Pushkin Street, 355000, Stavropol, Russia. E-mail: gavandrew@mail.ru.
Nataliya V. Malikova,
Secondary school no. 1, 47, Sovietskaya street, Grachevka village, Stavropol region, 356250, Russia. E-mail: nmalik_@mail.ru
Relevance of work. Extraction of geo assets requires the development of new technological solutions for their production, for example, the coating of drill screws with anticorrosion agents. A solution to this problem is possible using diffusion lacquer coating (corrosion inhibitor) on underlying surface (auger). The mathematization of such physical process as diffusion growth of thin films on the underlying surface is currently unexplored. In mathematical models the question on the existence and uniqueness of the solution of boundary-value problems describing the specified physical process often arises. Many domestic and foreign scientists have studied the analytical and numerical methods for solving the initial-boundary value problems, in which it is originally explicitly or implicitly assumed that the solution of the problem exists and it is unique. As a rule, the authors of publications devoted to various problems of mathematical modeling of diffusion, either do not address this question at all (about the existence and uniqueness of the solution) or refer to the classic works without good reason. Therefore, the studies on solvability of boundary value problems carried out in the paper are relevant. The aim of the research is to develop the criteria of resolvability (existence and uniqueness) of the boundary problems arising at mathematical modeling of the thin-film structure growth on underlying surface in various spaces.
Research methods. The achievement of a goal is based on correct use of results and methods of the equations of mathematical physics, the integrated equations, the mathematical analysis, the equations in private derivatives, physics of a solid body, a crystallography. Results. The authors have studied the resolvability of the boundary problems describing the diffusive growth of thin films on substrates; developed the criteria of existence and uniqueness of the solution of the specified tasks in various spaces. Conclusions. At mathematical modeling of diffusive growth of a thin film on underlying surface the authors developed the theorems (criteria) providing resolvability (existence and uniqueness of the decision) of initial-boundary tasks. The paper considers the boundary problems for cases of full reflection and absorption of atoms of a film by the underlying surface. The present article is of considerable interest in applied research, and allows answering a question:: is it possible to proceed immediately to a numerical (or possibly analytic) solution of the specific boundary value problem describing the diffusion growth of thin-film structures on substrates, or further carry out research on its regularization.
Key words:
Resolvability, boundary problem, thin film, underlying surface, substrate, diffusive growth.
REFERENCES
1. Boltaks B.I. Diffuziya i tochechnye defekty v poluprovodnikakh [Diffusion and dot defects in semiconductors]. Leningrad, Mir Publ., 1972. 362 p.
2. Verlan A.F., Moskalyuk S.S. Matematicheskoe modelirovanie ne-preryvnykh dinamicheskikh system [Mathematical modeling of continuous dynamic systems]. Kiev, Naukova dumka Publ., 1988.362 p.
3. Galay E.O. Matematicheskaya model obrazovaniya plenok na pod-lozhkakh [Mathematical model of formation of films on substrates]. Obozrenie prikladnoy i promyshlennoy matematiki, 2005, vol. 12, Iss. 4, p. 932.
4. Samarskiy A.A., Gulin A.V. Chislennye metody matematicheskoy fiziki [Numerical methods of mathematical physics]. Moscow, Nauchnyy mir Publ., 2000. 358 p.
5. Semenchin E.A., Tarasenko E.O. Modelirovanie protsessa diffuzii v poluprovodnikakh [Diffusion modeling in semiconductors]. Obozrenie prikladnoy i promyshlennoy matematiki, 2007, vol. 14, Iss. 2, pp. 352-353.
6. Tarasenko E.O., Gladkov A.V. Analiticheskie i chislennye reshe-niya nekotorykh obratnykh zadach v ramkakh matematicheskoy modeli rosta tonkikh plenok na podlozhkakh, uravnenie kotoroy
dopuskaet reshenie gaussovogo vida [Analytical and numerical solutions of some return tasks within mathematical model of growth of thin films on substrates which equation allows the solution of a Gaussian type]. Parallelnaya kompyuternaya algebra i ee prilozheniya v novykh infokommunikatsionnykh sistemakh. I Mezhdunarodnaya konferentsiya, sbornik nauchnykh trudov [Parallel computer algebra and its appendices in new infocommunica-tion systems. The I International conference, collection of scientific works]. Stavropol, 2014. pp. 93-98.
7. Tarasenko E.O., Gladkov A.V. Analiticheskoe reshenie kraevoy zadachi o vosstanovlenii moshchnosti istochnika v matematiches-koy modeli rosta tonkikh plenok na podlozhkakh [The analytical solution of a boundary problem on restoration of a source power in mathematical model of thin films growth on substrates]. Parallelnaya kompyuternaya algebra i ee prilozheniya v novykh info-kommunikatsionnykh sistemakh. I Mezhdunarodnaya konferent-siya, sbornik nauchnykh trudov [Parallel computer algebra and its appendices in new infocommunication systems. The I International conference, collection of scientific works]. Stavropol, 2014. pp. 99-103.
8. Tarasevich Yu.Yu. Matematicheskoe i kompyuternoe modelirovanie [Mathematical and computer modeling]. Moscow, Editorial URSS Publ., 2002. 144 p.
9. Yatsenko Yu.P. Integralnye modeli sistem s upravlyaemoypamy-atyu [Integrated models of systems with the operated memory]. Kiev, Naukova dumka Publ., 1991. 284 p.
10. Aubry S. Toward a Rigorous Molecular Theory of Metastability. Ferroelectrics, 1980, vol. 24, pp. 53-59.
11. Metiu H., Kitahara K., Ross J. Statistical mechanical theory of the kinetics of phase transitions. Fluctuation phenomena. Studies in statistical mechanics. Eds. Montroll E.W., Lebowitz J.L. Amsterdam, North-Holland Publish. Comp., 1979. Vol. 7, pp. 229-291.
12. Nishiama H. Martensitic Theory Transformations. New-York, Academic, 1981. 315 p.
13. Robledo A. The liquid-solid transition of the hard sphere system from uniformiry of the chemical potential. J. Chem. Phys., 1980, vol. 72, pp. 1701-1712.
14. Marchuk G.I. Matematicheskoe modelirovanie v probleme okruzhayushchey sredy [Mathematical modeling in an environment problem]. Moscow, Nauka Publ., 1982. 320 p.
15. Semenchin E.A. Analiticheskie resheniya kraevykh zadach v ma-temeticheskoy modeli atmosfernoy diffuzii [Analytical solutions of boundary problems in mathematical model of atmospheric diffusion]. Stavropol, SKIUU Publ., 1993. 141 p.
16. Fridman A. Uravneniya s chastnymi proizvodnymi paraboliche-skogo tipa [The equations with private derivatives of parabolic type]. Moscow, Mir Publ., 1968. 427 p.
17. Ladyzhenskaya O.A., Uraltseva N.N. Lineynye i kvazilineynye uravneniya ellipticheskogo tipa [Linear and quasilinear equations of elliptic type]. Moscow, Nauka Publ., 1973. 576 p.
18. Ladyzhenskaya O.A. Kraevye zadachi matematicheskoy fiziki [Boundary problems of mathematical physics]. Moscow, Nauka Publ., 1973. 406 p.
19. Afanasev V.N., Kolmanovskiy V.B., Nosov V.R. Matematiches-kaya teoriya konstruirovaniya system upravleniya [Mathematical theory of designing control systems]. Moscow, Vysshaya shkola Publ., 2003. 614 p.
20. Berlyand M.E. Sovremennye problemy atmosfernoy diffuzii i za-gryzneniya atmosfery [Current problems of atmospheric diffusion and pollution of atmosphere]. Leningrad, Gidrometeoizdat Publ., 1975. 448 p.
Received: 9 December 2015.
