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

На правах рукописи

Кротов Герман Сергеевич

Модельные представления процесса теплопроводности в области с движущейся границей

Специальность 05.13.18 —Математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

31 СКГ 2013

Москва 2013

005536774

005536774

Работа выполнена на кафедре «Высшей и прикладной математики» ФГБОУ ВПО Московского Государственного Университета тонких химических технологий имени М.В. Ломоносова.

Научный руководитель

Официальные оппоненты:

Ведущая организация:

Заслуженный деятель науки РФ, доктор физико-математических наук, профессор Карташов Эдуард Михайлович

Чернова Татьяна Александровна

доктор технических наук, профессор кафедры «Электроники и информатики» «МАТИ - Российского государственного технологического университета имени К.Э. Циолковского (МАТИ)»

Ломовской Виктор Андреевич доктор физико-математических наук, профессор, заведующий Лабораторией структурообразования и дисперсных систем Института Физической химии и электрохимии РАН

Московский авиационный институт (Национальный исследовательский

университет).

Защита состоится 21 ноября 2013 года в 14.00 на заседании диссертационного совета Д.212.110.08 при «МАТИ - Российском

государственном технологическом университете имени К.Э. Циолковского (МАТИ)» по адресу: 121552, г. Москва, ул. Оршанская, д.З.

С диссертацией можно ознакомиться в библиотеке «МАТИ - Российского государственного технологического университета имени К.Э. Циолковского (МАТИ)» по адресу: 121552, г. Москва, ул. Оршанская, д.З.

Автореферат разослан 15 октября 2013 г.

Учёный секретарь диссертационного совета Д.212.110.08 кандидат физико-

математических наук

Спыну Марина Валерьевна

Общая характеристика диссертации.

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

В диссертации сосредоточено основное внимание на модельных представлениях процесса теплопроводности в областях с движущимися во времени границами, функциональных конструкциях в качестве аналитических решений для конкретных законов движения границы, простых примерах иллюстративного и содержательного характера. В связи с этим в диссертации не рассматриваются вопросы качественной теории уравнений параболического типа, асимптотические подходы. Развитые в работе подходы позволили рассмотреть важную для физической теории прочности материалов проблему теплового удара в терминах динамической термоупругости при наличии движения во времени границы области термонапряженного состояния по линейному закону y(l) = l + vt и закону >'(/) = ß-v/F, где /, v, р — постоянные. Построены многочисленные графики и произведён сравнительный анализ термонапряженных состояний для этих типов областей.

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

Актуальность работы

До последнего времени в работах, посвященных данной тематике, предполагалось, что геометрия тела остается неизменной с течением времени. Однако во многих практически важных ситуациях это предположение о неизменности во времени геометрии тела не выполняется и возникает необходимость развивать модельные представления задач теплопроводности и термомеханики в областях с движущимися во времени границами, то есть находить температуры, напряжения и деформации в этих областях. Стоит отметить, что краевые задачи для этих случаев чрезвычайно сложны, и для их аналитического решения классические методы математической физики оказываются неэффективными, так как не удается согласовать решение параболического уравнения с движением границы. В результате остро актуальной становится проблема создания нового математического аппарата, позволяющего находить аналитические решения вышеуказанных модельных задач для различных законов движения границы. Указанная проблема определила направление исследований. И хотя работы академиков Лыкова A.B., Тихонова А.Н., Самарского A.A.,

Гринберга Г.А. и его сотрудников (Чекмарева О.М., Косс В.А.), Карташова Э.М. и других заложили в этой области прочный фундамент, мы все же ещё находимся в начале пути.

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

Стефана

дх

Л

= А —>0, выражающее баланс энергии при переходе среды из г-Н<) л

„ . дТ(х,А

одного агрегатного состояния в другое, или типа Стефана —-—-

= Ч,(/),/>0 в более

"НО

дх

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

Цель работы

Цель диссертации — развитие модельных представлений процессов теплопроводности и термоупругости в областях с движущимися во времени границами. Это достигается путём развития методов нахождения точных аналитических решений краевых задач теплопроводности и термоупругости для этих указанных областей, получением новых функциональных конструкций в качестве аналитических решений для конкретных законов движения границы, построением точных решений в областях с границей, движущейся по закону у{1) = Рл/7, проведением соответствующих численных экспериментов, развитием численных методов и построения комплекса программ для нахождения корней ряда сложных трансцендентных уравнений, возникающих при нахождении аналитических решений в областях с границей, движущейся по закону >•(/) = Р>/7.

Научная новизна

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

Получены точные аналитические решения указанных краевых задач теплопроводности в области П, = {х > > о} с помощью метода функции Грина.

Численно решена проблема нахождения нецелых корней />„ уравнения Ор{г) = 0 (г = Ор(д) — специальная функция параболического цилиндра) при фиксированном

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

В результате рассмотрения проблемы теплового удара при температурном нагреве упругих полупространств z>l + vt, t>0, v = const и z > [i4t, t >0, /? = const в квазистатическом случае выведена точная формула для этих напряжений, когда граница движется по закону z > Pit.

Практическая значимость

В диссертации разработан математический аппарат, позволяющий получать важную информацию об особенностях теплового и термонапряженного состояния, возникающего в упругих твердых телах, подвергающихся резким термическим воздействиям, в том числе при изменении геометрических размеров тел. Развитые подходы позволили рассмотреть важную для физической теории прочности проблему теплового удара в терминах динамической термоупругости при наличии движения границы. Изученные в диссертации закономерности могут быть использованы при разработке методов применения лазеров в технологических операциях и при производстве самих лазероактивных материалов; при исследовании синтеза и свойств высокопрочных термостабильных полимеров и эластомеров; при изучении термоупругих и динамических эффектов в проводниках и диэлектриках; при исследовании термических напряжений, возникающих в космических аппаратах при вхождении в плотные слои атмосферы и разработке термостойких покрытий для них, способных выдерживать экстремальные термические воздействия, в геологических породах, нефтеносных пластов; при изучении распространения термоупругих волн в мантии Земли; в реакторостроении и ядерной энергетике — как при проектировании реакторов, так и при математическом моделировании аварий на них; а также в ряде других фундаментальных и прикладных исследований. Большую практическую ценность работа имеет для исследования термических эффектов при распространении трещин, когда на берегах распространяющейся трещины (область с движущейся границей) задается постоянная температура, приводящая к разрушению материалов, механизмов или летательных аппаратов; при исследовании процессов плавления электрических контактов; при изучении процесса горения, воздействия электрической дуги на контакты; при исследованиях промерзании растворов, грунтов; при исследовании кинетического роста кристаллов. Разработанные методы решения краевых задач могут найти применение в разных областях техники, связанных с авиацией и космическими исследованиями, в химической технологии и других сферах деятельности

человека, например, в МФТИ, ИПМ РАН, РНЦ «Курчатовский институт» и других научных центрах и организациях, а также при обучении магистров и студентов, обучающихся по направлению (специальности). Результаты диссертации, касающиеся краевых задач в областях с движущимися границами, используется в учебном курсе «Прикладная термомеханика» в ФГБОУ ВПО Московского Государственного Университета тонких химических технологий имени М.В. Ломоносова при обучении специалистов по направлению «Прикладная математика» (специализация «Математическое моделирование»).

Апробация работы

Материалы диссертации докладывались на Международной конференции "Математические модели физических процессов" (Таганрог, 2007 г.), на Международной научно - технической конференции "Повышение эффективности теплообменных процессов и систем" (Вологда, 2005 гг.), на городском семинаре по проблемам тепломассопереноса в МИТХТ им. М.В. Ломоносова (2000-2011 гг.), на 3-й Российской национальной конференции по теплообмену (Москва, 2002 г.), на 5-й Российской национальной конференции по теплообмену (Москва, 2010 г.).

Личный вклад автора

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

Публикации

По теме диссертации опубликовано 9 работ.

Структура и объем диссертации

Диссертация состоит из введения, шести глав, заключения и списка цитируемой литературы. Работа содержит 181 страницу машинописного текста, 32 рисунка, список литературы из 410 наименований.

На защиту выносятся следующие основные положения:

1. Использование метод функции Грина для областей цилиндрического и нецилиндрического типа.

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

3. Решена проблема нахождения нецелых корней рп уравнения Dp(z)= 0 при фиксированном значении аргумента. В результате рассмотрения получена Таблица первых 10 корней этого уравнения с точностью до седьмого знака после запятой для каждого z с шагом 0,1 в интервале от 0 до 10.

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

5. Рассмотрены математические модели термоупругости и проблема теплового удара в области, границы которой неподвижны во времени, области G, = {z > l + vr,t > 0} и

области П, = {х: > pV7,/ > о}.

6. Установлено влияние движения границы на поведение основной компоненты тензора напряжения а_ для упругого полупространства, граница которых неподвижна или движется по линейному закону.

7. Рассмотрено поведение компонент <т„ = ггуу тензора напряжений в полуограниченных областях с неподвижной границей и с границей, движущейся по линейному закону. Приводятся их точные решения. Проведен анализ компонент ст„ и <та= <т}у в этих областях.

8. Дано описание термической реакции упругих полупространств z>l + vt, t>0, v = const и z> p~Jt, t> 0, p = const в квазистатическом случае и выведены аналитические формулы для различных компонент тензора напряжений при температурном нагреве в частично ограниченной области, когда граница движется по закону z > pit.

2. Метод отыскания функции Грина в области

Определение её

функции в области iJ, = {с > р4t,t > о}.

Краткое содержание диссертации

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

В первой главе даётся постановка краевой задачи нестационарной теплопроводности в нецилиндрической области, приводится обзор современных аналитических методов решения краевых задача нестационарной теплопроводности. Приведены две теоремы об основных свойствах тепловых потенциалов, причем первая доказывается для общего случая движения границы области, а вторая для случая y(t) = I+vt.

Во второй главе рассматривается метод функции Грина как для канонических областей, так и для областей нецилиндрического типа. Объясняется разница в постановках краевых задач для функций Грина для этих областей. Показывается, что в областях с движущимися границами не сохраняется эквивалентность в записи граничных условий в постановках задач для функции Грина С?(х,/,х',т) по (x,t) и (х',т) в отличие от цилиндрических областей.

Для температурной функции T(M,t) в цилиндрической (канонической) области, удовлетворяющей условиям указанной ниже задачи:

^МЛ = a&T{M,t)+ — f{M, l\ М(х,у, z) <е D, t > О, at cp

T(M,t\_0 = <t>0(M\M(x,y,z)eD, (1)

дТ(М><) + PiT(M ,)= f),(p(M,tlM(x,y,z)eS,t > 0, on

рассмотрены 1-я, 2-я и 3-я краевые задачи в зависимости от условий на границе: если параметрам Pj, Р2 и Р3 задать значения Pj= О, Р2= ~Р3= 1 — в случае первой краевой задачи, Р2= 0, pj= Р3= 1 — в случае второй краевой задачи, Pj= 1, р2= Р3=-й (Л — относительный коэффициент теплообмена) в случае третьей краевой задачи.

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

Г(М,0 = ¡¡¡<b0(P)G(M,P,t,0)dVp +ajjj G-^-Г^-) dxdar +

D 0 S V dip JrtS

1 ' (2) + — J jjj F(p, t)G(M, P, и T)dVPdT,

CP 0 D

если известна соответствующая функция Грина G{M,P,t, т) для данной области как решение более простой задачи для однородного уравнения, соответствующего задаче (1) с однородными граничными условиями того же вида, что и для исходной задачи

8G(MP,t,т) aA t(G(M;Pt/;T)i M(x,y,z) е D, t >х, (3)

от

G(M,P,/,T)|i=t =5(Л/,Р), M(x,y,z) e D, P(x,y,z) e Д (4)

P, aG(A/'P'f'X) +P2G(M,P,f,t) = 0, M(W)eS,Ot. (5)

Здесь

— дельта-функция Дирака. Если рассматривать t, т) как функцию

точки Р и времени т, то для нахождения функции следует решать

эквивалентную задачу для сопряженного уравнения

8G(M, P,t,x) + дД ри т) = 0> p(x,y,z) ей, к/, (6)

дг

G(A/,P,i,x)|(=T =5(М,Р), P(x,y,z) е Д M(x,y,z) е D, (7)

р, 5g(m'''т) + p2g(m, л г, т) = 0, P(x,y,z) £ S, х <t. (8)

дпр

Если область D ограничена, то функция Грина имеет вид

yt;(,_г) = G(MJ>,t-x) = 2-,-и—р-е , (9)

"=> IKII

где (М) и у„ — собственные функции и собственные значения соответствующей однородной задачи для (6)—(8) (задача Штурма-Лиувилля)

Пусть теперь область nt имеет вид n,= {{x,t): х >){t), i > 0}, где у{1) — непрерывно дифференцируемая функция. Тогда T(x,t) — температурное поле в области П( может быть найдено в результате решения задачи

(УР д^Р

— = a—+f(x,t), х > y(t), t> 0; (10)

at дх

ГОгД,0 =Ф0(х), x>y(0), y(0)>0. (11)

с граничными условиями первого рода T(x,t) = <[>(')>' > 0, второго рода

dT(x,t) |

8х

8P(x,t)

= <p(/),i > 0 или третьего рода —;—-

дх

= h(T(x,t)-^0) t > 0. Для всех

случаев |7"(х,/)| < оо. Здесь /(х,1) е С°(Г2,); Ф0(*) е С'(Г2,); ^>(/)еС0[0;+оо); искомое решение Т(х,I) е С2(П,)ПС°(Д); £гас1, Т(х,1) е С°(Ц).

Функция Грина С(х, 1,х',х) — функция температурного влияния мгновенного точечного теплового импульса мощности () = ср в точке х' в момент времени т. Для функции Грина С(х,1,х', т) по переменным (х',т) имеем следующую задачу

дй 8гС , , .

= , х' > >(т), т <Г, (12)

С1 дх

С(х,1,х\т^г1=3(х -х), х'>у(т). (13)

с краевыми условиями в случае первой краевой задачи б(дг,/,х',т)|х, } = 0,т<?; в случае

(дв 1 Оу, ----(

\дх а ск

второй краевой задачи —----С I =0,т</ или в случае третьей краевой задачи

^ + — — = 0, т < Г . Во всех случаях |С(х,/,х',т)|<оо.

& V а л; Л'=,м

Основная формула, дающая интегральное представление аналитического решения уравнения (10) в области П, = {(*,/): х > >{/),' > 0} для всех типов краевых условий имеет вид:

МО)

■4МННН.....

' «'^(г) о >•(»•)

Построение функции Грина С(х,Г,х',т) для второй и третьей краевых задач, если исходить из постановки задачи по переменным , г), представляет серьёзные трудности и для большинства законов движения границы технически не осуществимо. Если исходить из постановки задачи по переменным (*,/), то ситуация существенно упрощается. Действительно, рассмотрим, функцию С(х,1,х' ,г), определяемую условием

8в 82в

-= а —

3/ дх

= х>уЦ\ >>г, (15)

т)|м=8(л-*'), х>у(х). (16)

и

в случае первой краевой задачи С(х,Г,х',т)| =0 > т ; в случае второй краевой задачи

= 0,Г >т.Во всех

8G дх

.дх

X=yU) 4

= 0,/ > т; в случае третьей краевой задачи

x-yW

случаях |с?(*,/,л:',т)| < +оо. Во §2 второй главы доказывается, что С(х,г,х',г) = С(х,/,х',г).

Таким образом, функция С(х,/, х', г) может быть получена как решение какой-либо из двух эквивалентных задач для уравнений (12) и (15) с соответствующими им краевыми условиями, приведенными выше. Из сказанного следует, что в областях с движущимися границами не сохраняется эквивалентность в записи граничных условий в постановках задач по (х,1) и (х\ т) в отличие от классических цилиндрических областей.

Искомые интегральные представления для аналитических решений Т(х,{) в ограниченной области П, = {(;с,г):О<л:<><(0,г >0} и полуограниченной области О, = {(х,/):х> >■((),; >0} имеют вид

2 I

T(x,t) = \T(x'fi)G(x,t,xfi)dx'

л(0) 1=0 0

дТ(х',х) , _,, , , a,i 8х, +арТ(х,х)

х| 7,1 DG(X^,X',X) +7aG(x.t,x',x) [[■ dx,

(17)

где

ал = у,2 = 0;а,2 = (- = 1 — первая краевая задача;

а12 — у„ =0;ал =(-l) ';у/2 =1 —вторая краевая задача; а,, = (- 1У;а,2 = й,;у,, = 0;y¡2 = 1 — третья краевая задача. Из (17) нетрудно получить аналогичные представления и для области [у(<), оэ); t > 0; х'= _у2 (т) —> м с учетом того, что при этом функции Г и G и их производные по х' стремятся к нулю.

Формула (17) позволяет в следующей главе построить точные аналитические решения первой и второй краевых задач в области П, = > р4t,t > о}.

В третьей главе приведены постановки первой и второй краевых задач нестационарной теплопроводности и их решения в области П, = (г > Р-Jt,t > о}, записанные в интегральной форме через функции Грина. Указывается метод отыскания функций Грина, а затем приводится их аналитический вид для первой и второй краевых задач. Действительно, функция Грина для этих задач удовлетворяет условиям

зс аб „-#=-

— = 0-- Л-'' ; (18)

¿я йг -

с(* ' * = Л^Н (19)

а,С(л:,/,*',т)+ а2 'т) | = 0,/ >т>0; |С(х,/,х',т)| < ас,р77 <х<а>,Г >0. (20)

При а, = 1,а2 = 0 получаем первую краевую задачу, при а, = 0,а2 = 1 — вторую краевую задачу. Формулируется и доказывается теорема. Теорема. Задача

с! 1Г(х)-Ш{х) ,„, ч л, ч., т=

——1— + а -^ + д Мх)= ^(гуЦг-г) г>

ах -Ж- - --

где а, = сог^ (/ = ),

эквивалента следующей

<1 1Г(х)-гЛУ(х) , ,

ах -ас- -

<т{х)

л

¿1Г(х)

сЫ

= А{т),

(21) (22)

(23)

(24)

(25)

то есть решения задач (21)—(22) и (23)—(25) совпадают.

Используя эту теорему, в главе 3 выводится аналитический вид функций Грина задач (18)—(20) и показывается, что они могут быть представлены в двух формах: в интегральной форме и форме ряда.

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

хл

8а! т I

й

41а

е ' ар

или в виде ряда

4—1 - ■

Ц I I . -я.-И рг~

дР

хО,

где />„ — нецелые корни уравнения

а) —

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

.т

в(х,1,х\т)=е \ Г(р + 1

Л гг! • / /*1 1Т-Г •

4тыатп ,

Ф

или в виде ряда

С{х I .

где рп — нецелые корни уравнения

О.

4~а

В §5 главы 3 рассмотрена проблема нахождения нецелых корней рп уравнения

О,(г) = 0

при фиксированном значении аргумента 2. Как несложно заметить, зная корни уравнения (32) при фиксированном значении аргумента г, можно легко получить корни уравнений (28) и (31), значения которых необходимо знать, чтобы построить функции Грина первой и второй краевых задач в виде ряда.

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

О»

-п."

Рисунок 1. График функции параболического цилиндра Ор(г) в зависимости от действительного р в интервале от 0 до 25 и действительного г в интервале от 0 до 10.

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

Для отыскания корней уравнения (32) использовался метод половинного деления. Используя этот метод, автором была написана программа для отыскания корней в программной среде МаЛсас! 11.0а.

rootD(a,b,z,e) := n

while Ial — ЫI >-2

middle«

al + bl

al <- middle if sign(D(al ,z)) = sign(D(middI^z)) Ы middle if sign(D(bl,z)) = sign(D(middl^z)) ' middle Л

al bl D(al,z) D(bl,z))

Здесь £>(p,z)— функция параболического цилиндра, sign(2) — функция знака. Аргументами определяемой функции rootD{a,b,z,e) являются: а — левая граница отрезка, где находится отыскиваемый корень уравнения (32), Ь — правая граница отрезка, где находится отыскиваемый корень уравнения (32), z — фиксированное значение аргумента функции параболического цилиндра Dp{z), при котором решается уравнение (54), е — погрешность вычисления. Значением функции тооЮ[а,Ь,г,е) является вектор размера 1x6. Его первой координатой является середина последнего интервала, получающегося в результате делений отрезков пополам, в котором находится отыскиваемый корень. Второй координатой является количество делений исходного отрезка, третьей — левая граница последнего отрезка, четвертая — правая граница последнего отрезка, пятая — значение функции параболического цилиндра Dp(z) на левой границе последнего отрезка, шестая — значение этой функции на правой границе последнего отрезка.

На рисунках 2 и 3 приведены графики функции Dp(z) при z = 0.5 и при z = .

Рисунок 2. График функции параболического цилиндра Е>р(г) при г = 0.5 в зависимости от действительного р в интервале от 0 до 25.

ОМ

1

\У 7 2 > 10 11 12 13 14 1) 1« 17 19 а 31 1 зз г- г

Рисунок 3. График функции параболического цилиндра Ор(:) при г = 1 в зависимости от действительного р в интервале от 0 до 25.

Далее в этой главе построена таблица первых 10 корней уравнения (32) с точностью до седьмого знака после запятой для каждого г с шагом 0,1 в интервале от 0 до 10.

2 Р Р Р Р Р Р Р

0,1 1,0817496 3,1216944 5,1516245 7,1765615 13,2360156 15,2527293 17,2683983 19,2831968

0,2 1,1674515 3,2474369 5,3073019 7,3571803 13,4760908 15,5095182 17,5408562 19,5704531

0,3 1,2571485 3,3772190 5,4670482 7,5418695 13,7202351 15,7703756 17,8173819 19,8617765

0,4 1,3508809 3,5110918 5,6308788 7,7306421 13,9684580 16,0353102 18,0979836 20,1571750

0,5 1,4486867 3,6490658 5,7988093 7,9235111 14,2207688 16,3043307 18,3826696 20,4566563

9,7 29,9669726 35,4689880 40,1547642 52,2501949 55,9258545 59,4875355 62,9670319 66,3068749

9,8 30,4988793 36,0323776 40,7430228 52,8964401 56,5882642 60,1653036 63,6537034 67,0354127

9,9 31,0356542 36,6005608 41,3360246 53,5473278 57,2552160 60,8463446 64,3479075 67,7439982

10,0 31,5772987 37,1735398 41,9337720 54,2028892 57,9276483 61,5360962 65,0576091 68,5563119

В §6 приведен пример использования метода Грина для конкретной задачи, аналитические решения и графики, которой приведены в безразмерных переменных

а!

I I л/а где / — масштабная единица длины.

Постановка этой задачи имеет следующий вид

-= —т-,ул/Ро <г< оо,Ро> 0:

дРодг2

=> = К>Г< 1; < +00'

а её решением является функция

¡Г (г, Ро) =

Рое^ '---:-г-^-,0<Ро<1

О

" д_ дР

2 ~Р,

■Ро2 ,Ро>1

(33)

(34)

(35)

(36)

Графики решения (36) имеют вид

Цг(2,¥о)

Рисунок 4. Изменение безразмерной температурной функции 1У{г,¥о) при у = 1 в

зависимости от 2 при: - Бо = 0.3;

.....Ро = 0.4; — Ро = 0.5.

Рисунок 5. Изменение безразмерной температурной функции ^(г^о) при у = 1 в

зависимости от 2 при:- Бо = 3;

.....Ро = 5; — Ро=7.

Рисунок 4 описывает распределение температуры во время нагрева на границе: чем ближе точка к границе, тем выше её температура. Рисунок 5 указывает на распределение температуры, когда нагрев на границе прекращен. При удалении от границы виден рост температуры до некоторого максимального значения, а затем стабилизация температуры со стремлением к нулю на бесконечности. Вероятно, такой эффект обусловлен наличием движущейся границы.

Стоит отметить, что эта глава является одной из ключевых в диссертации.

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

*=+ ОО А — 11

^млф" — в случае области О, = р<х<р77,/ >о|,

п ~т

X Мк — в случае области О, = {»: > Рл/7,Г > 0}, где £ = / = 7-1.

*=-«. ~ л/2 аг

Функциональные коэффициенты Nк (¡;) и Л/, (¡;) в находятся из удовлетворения

граничным условиям.

В первом параграфе данной главы рассматривается метод обобщенных рядов при решении краевых задач нестационарной теплопроводности в области 0 < х < у у/2 а/,? > 0, а во втором — для области х > у42хй,1 > 0.

Показано, что в частных случаях эти ряды сводятся к уже известным решениям, но имеющим иные функциональные конструкции. Например, если в области х > у ^2а/,/ > 0 в случае первой краевой задачи граничная функция является одночленом = Л/", то, используя метод обобщенных рядов, получим решение в виде

" ■ Д < {—1 г(х>,) = е"^• 07)

С помощью (38}—(39)

е2 />_„_, (в) + е 2 £>_„_, (-¿г)

(38)

(39)

42л

решение (37) приводится к

А ■ ;'"ег1с — т( \ \2yjat I ;

П*>') =-/ х >2- (40)

'ЧЙ

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

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

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

dt дхJ v ;

= A^-,t>0; (42)

■=>(') dt

(43)

^= ' > О (у(0) = О); Г(х,^ =ф(4'>0, (44)

при т = 0 — в декартовом, от = 1 — в цилиндрическом, т = 2 — в сферическом случае.

Решение задачи Стефана (41)—(44) в декартовом случае с граничными условиями П>.')|,=0 = Гс=соп,51 и ч/(/)=7п может быть получено в различных эквивалентных формах.

Например, решая эту задачу в декартовых координатах (т = 0) с помощью обобщенного интегрального преобразования

T(p,t)= ¡T(x,t)st{xjp)dx,

(45)

кг- я

где р=<у+на - комплексное число,--< arg Jp < —, получим решение в виде

4 4

(-1/1 а"42

к

. ^ . тсс ^р2» П V Р т!

Г(дг,/) = Х«„5ш-Т, где -^--X т' .... (46)

„=1 ЯП *! т=* а (2/л + 1)!

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

получаем закон движения границы области >'(/) = р^/7, а также решение в виде

Ь^) 2 „4-оал(2»+1)!е/" I1 -1 -Я1

Тп-^.^еЦф(-Т=11 = Фо + — Т^Фf-i=l.

2

Р

V2-Ja) \2-Jrt)

2 \2-Jai)'

(47)

Г 4а J Р Р2"'1»11 Используя соотношение sitae "Ф —¡= - ¿^ ——-—, (46) приводится к виду (47).

\2\!а J „„о а (2т +1)!

В §4 построены графики безразмерной температурной функции сначала при фиксированных временах, а затем при фиксированных значениях пространственной переменной и проведен их сравнительный анализ. Также построена с помощью вычислительных методов таблица возможных значений р в законе движения y(t) = р4t.

В шестой главе изучается проблема теплового удара для упругого полупространства z > / + V/ (для общности / * 0), t> 0, v = const в динамическом и квазистатическом случаях, а также для упругого полупространства z >pjt, t> 0, р = const. Рассмотрение

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

законы близки (рисунком 6), т.е. р-\/7 к ^ 7 = v/.

Ж)

Рисунок 6. Поведение графиков функций = л/7 - и у(г) = VI ----- в

окрестности точки 0.

Сначала даётся общая математическая постановка задачи термоупругости для произвольной области С, ={МеС,;(> 0}, затем рассматривается интересующая нас математическая модель о тепловом ударе в терминах динамической термоупругости

дга2! 1 дга2! (1 + 1/) с>2Т(1,1) , , -г^-^г-= ~-~атР-г>шI г >0

а.г V* а2 (1-у) т аг

! \ „ 8ст„(г,{) „ =0, -'

= 0,2 > >{о), «тв(г>/]| = = 0

(48)

(49)

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

0Г дгТ <\ п

Т(2,1)[_0=Т0, 2>у{0) и рассмотреть 3 режима: Т(г,1)= ТС, I >0 —температурный нагрев;

(50)

(51)

Яг

6Т(2,1)

д2

-цт, / > 0 — тепловой нагрев;

6Т(2,1)

~у(у)

дг

Ы')

Л(71г->(,) - ТС ) ' > 0 — нагРев

средой. Во всех случаях |Г(г,/)| < +°о, 2 > у((\ (>0.

Задача (48) - (49) содержит большое количество параметров, что неудобно для численного анализа, поэтому эту задачу целесообразно рассматривать в обобщённых переменных.

В §3 приведено решение динамической задачи для компоненты а_ тензора напряжений в полуограниченной области, граница которой движется по закону >>(<)=/+у/ при температурном нагреве

a„(z,Fo) = a£,'(r,Fo) +

О,

при 0<Fo<—,

Яо

a„(z,Fo) при — < Fo < —, а„ Ре

(52)

где

2[а02-Ре2 I, 2VFb

-Ре-а„

«о-Pe

-■e_a°2erfc

-Ре + а„

а„-Ре

- • e°°'erfc

z-2 • а0 • Fo 2-n/Fo

/ z + 2 • а0 • Fo 2,/Fo

(53)

-г Ре-Ял

. . ^(«o-^XFo-ÍFoJgjg.

«oj

k/Fo-<5Fo

-Pe + a0

2__^«„(«o+PíXFo-ÍFo),

a„ +Pe

erfc

Pe 2 "

—+ a„ ]VFo-¿)FO

(54)

Проводится сравнительный анализ поведения компоненты тензора напряжения а для областей, граница которых неподвижна, и когда граница движется по линейному закону. В §4 рассматривается вопрос поведения компонент = ап тензора напряжений в полуограниченных областях с неподвижной границей и с границей, движущейся по линейному закону. Приводятся точные решения этих задач и рассматривается вопрос о поведении компонент <тв = <т тензора напряжений для упругого полупространства в

зависимости от движения границы. Также проводится сравнительный анализ компонент а.. и а„ = <т>у в этих областях, который даёт следующий результат: напряжения вдоль оси z опаснее напряжений вдоль осей х и у, так как превосходят их примерно в два раза.

В §5 рассматривается вопрос термической реакции упругих полупространств z>l+vt и z>/W7, при |>0 и v,(i = const в квазистатическом случае при температурном нагреве. Выводится точная формула для напряжений, когда граница движется по закону z > -JVo

erfcf-

v-1

aJz,Fo)=a„ (z,Fo)=^—(55)

erfc

7

0 10 20 р0 Рисунок 7. Изменение напряжения сг„ (55) при температурном нагреве, при движении границы по корневому закону в квазистатическом случае в зависимости от безразмерного времени Ро в интервале от 0 до 30, V = 0.3, Ре = 0.000893 в сечениях: г = 3 - ; г = 5 -----; г = 7 — .

Один из важных новых результатов, рассмотренных в этой главе, — поведение компонент аа = а>у тензора напряжения при движении границы по корневому закону в квазистатическом случае, что хорошо иллюстрирует рисунок 7. Напряжение всё время остаётся сжимающим, растёт не резко, а плавно, незначительно меняется при росте времени. Из этого можно сделать вывод, что квазистатический случай является менее опасным с точки зрения напряжений, возникающих при температурном нагреве.

Рассмотрено поведение компоненты а при малых временах в различных областях.

О >10"' 1>10"' 640'' 1»10'' 001

Рисунок 8. Изменение напряжения оа (58) при температурном нагреве, при движении границы по корневому закону в квазистатическом случае в зависимости от безразмерного времени Ко в интервале от 0 до 0.01, 1/ = 0.3, Ре = 0.00083 в сечениях:

г = 0.01 - ; г = 0.015 .....;

г = 0.02

р0

Рисунок 9. Изменение напряжения аш при температурном нагреве, при движении границы по линейному закону в квазистатическом случае в зависимости от безразмерного времени Бо в интервале от 0 до 0.01, V = 0.3, Ре = 0.00083 в сечениях:

г = 0.01- ; 2 = 0.015 .....;

г = 0.02

Легко заметить практически полное совпадение рисунков 8 и 9, построенных при разных видах движения границы: линейном и корневом.

Рисунки 8 и 9 иллюстрируют другой важный новый результат, приведённый в этой главе, который заключается в следующем. Предполагалось, что, так как при малых временах именно линейный закон наилучшим образом приближает корневой закон, то и точное решение задачи о нахождении напряжений в случае линейного закона будет наилучшим образом приближать решение этой задачи в случае корневого закона движения границы. Это предположение оказалось верным в квазистатическом случае.

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

Основные выводы и заключение

В диссертации последовательно рассмотрен и решен следующий комплекс вопросов:

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

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

зависит от разности (I -т), а в нецилиндрической области для функции Грина определяющим будет / и т, т.е. не только время действия (/-х), но и т — момент возникновения теплового импульса мощностью () = ср. Также показывается, что в областях с движущимися границами не сохраняется эквивалентность в записи граничных условий в постановках задач по (*,/) и (дг',т) в отличие от цилиндрических областей. Эти свойства функции Грина являются настолько важными, что упустив их из вида можно получить неверные решения. Приведены постановки ряда краевых задач нестационарной

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

3. Решена проблема нахождения нецелых корней рп уравнения £)Дг) = 0 при фиксированном значении аргумента и построена таблица первых 10 корней этого уравнения с точностью до седьмого знака после запятой для каждого 2 с шагом 0,1 в интервале от 0 до 10 с помощью программы, написанной в среде МаЛсас!.

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

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

теплопроводности и их решения в области

записанные в

линейному закону. Приводятся точные решения этих задач и рассматривается вопрос влияния движения границы на поведение компонент а„ = ега тензора напряжений. Также проводится сравнительный анализ компонент <т_ и = в этих областях, который даёт важный для практики результат. Рассматривается вопрос термической реакции упругих полупространств z>l + vt, t> 0, v = const и z> p4t, / > 0, /? = const в квазистатическом случае. Выводится точная формула для напряжений, когда граница движется по закону z Показывается, что напряжение всё время остаётся

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

Список опубликованных работ

По теме диссертации опубликовано 9 работ.

1. Карташов Э.М., Кротов Г.С. Аналитическое решение однофазной задачи Стефана // Математическое моделирование, 2008, том 20,3, стр.77-86.

2. Карташов Э.М., Кротов Г.С. Функции Грина в задачах нестационарной теплопроводности в области с границей, движущейся по корневой зависимости // Известия РАН, Энергетика, №4,2006, стр. 134-149.

3. Карташов Э.М., Кротов Г.С. Проблема Стефана в вырожденной области // Математическое модели физических процессов. Таганрог: 2007. Том 1. С. 103-109.

4. Кротов Г.С., Карташов Э.М. Функция Грина для области с движущейся границей // Математические модели физических процессов. Таганрог: Издательство Таганрогского государственного педагогического института. 2005. Том 1. С. 89-95.

5. Кротов Г.С. Корни функции параболического цилиндра при фиксированном значении аргумента // Ученые Записки МИТХТ. Москва. 2005. Выпуск 14. С. 41-48.

6. Кротов Г.С., Карташов Э.М. Задачи нестационарной теплопроводности в области с границей, движущейся по корневой зависимости // МНТК. Вологда. 2005.

7. Кротов Г.С. Метод решение краевых задач теплопроводности для области, с границей, движущейся по параболическому закону // Третья Российская национальная конференция по теплообмену. Москва. 2002. Т. 7, стр. 161-163.

8. Кротов Г.С., Карташов Э.М. Метод решение краевых задач теплопроводности для области, с границей, движущейся по параболическому закону // Вопросы теории и расчета рабочих процессов тепловых двигателей. Межвузовский научный сборник №18. Уфа. 2000. С. 113-119.

9. Кротов Г.С., Ремизова О.И. Численное и аналитическое решения краевой задачи для уравнения параболического типа в нецилиндрической области // Вестник МИТХТ. Том VIII. №5. 2013

Кротов Герман Сергеевич Модельные представления процесса теплопроводности

в области с движущейся границей

Формат 60x90/16 Тираж 100 экз. Подписано в печать 15.10.2013 Заказ № 98 Типография ООО «Генезис» 8 (495) 434-83-55 119571, г. Москва, пр-т Вернадского, 86

Министерство образования и науки РФ ФГБОУ ВПО Московский Государственный Университет тонких химических

технологий имени М.В. Ломоносова

на правах рукописи

0420136406?

Кротов Герман Сергеевич

Модельные представления процесса теплопроводности в области с движущейся границей

Специальность: 05.13.18 — Математическое моделирование, численные методы и

комплексы программ

Диссертация на соискание ученой степени кандидата технических наук

научный руководитель Засл. деятель науки России, д.ф.-м.н., проф. Э.М. Карташов

Москва —2013

ОГЛАВЛЕНИЕ

Введение к работе........................................................................................................................5

ГЛАВА 1. Современные аналитические методы решения краевых задач нестационарной теплопроводности в областях нецилиндрического типа

Введение.......................................................................................................................................8

§1. Постановка краевых задач нестационраной теплопроводности............................................11

§2. Метод тепловых потенциалов...................................................................................................14

§3. Метод обощенного интегрального преобразования...............................................................18

§4. Метод дифференциальных рядов..............................................................................................21

§5. Метод функциональных преобразований. Новые законы движения....................................22

Выводы к главе 1..............................................................................................................................32

ГЛАВА 2. Интегральные представления аналитических решений краевых задач нестационарной теплопроводности в канонических областях и в областях с движущейся границей

Введение............................................................................................................................................34

§1. Метод функции Грина для цилиндрических областей канонического типа........................35

§2. Развитие метода функции Грина для областей нецилиндрического типа. Постановка

краевых задач. Интегральные представления........................................................................37

2.1. Интегральные представления аналитических решений краевых задач

теплопроводности в области 0 < х < ^ > О........................................................................38

2.2. Интегральные представления аналитических решений краевых задач

теплопроводности в области х>у(0> * ^ 0................................................................................43

Выводы к главе 2..............................................................................................................................47

ГЛАВА 3. Построение функции Грина и интегральных представлений решений задачи теплопроводности для полуограниченной области с границей, движущейся по корневой зависимости

Введение............................................................................................................................................49

§ 1. Постановка первой и второй краевых задач............................................................................49

§2. Об одном приёме сведения основной задачи к эквивалентным............................................51

§3. Построение функции Грина первой краевой задачи...............................................................54

§4. Построение функции Грина второй краевой задачи...............................................................60

§5. Корни функции параболического цилиндра при фиксированном значении аргумента.....63

Таблица корней функции параболического цилиндра при различных фиксированных

значениях аргумента г.............................................................................................................69

§6. Интегральные представления первой и второй задачи для области О, = {х > /?л/м > о}

и численные эксперименты......................................................................................................76

Выводы к главе 3..............................................................................................................................82

ГЛАВА 4. Метод обобщенных рядов при решении задача теплопроводности с границей, движущейся по параболическому закону [о,ул/2а7] и |ул/2а/,ао)* > 0

Введение............................................................................................................................................83

§ 1. Метод обобщенных рядов при решении краевых задач нестационарной

теплопроводности в области 0 < х < > 0 .............................................................83

§2. Метод обобщенных рядов при решении краевых задач нестационарной

теплопроводности в области X > /V, t > 0..................................................................93

Выводы к главе 4............................................................................................................................102

ГЛАВА 5. Проблема Стефана

Введение..........................................................................................................................................103

§1. Метод дифференциальных рядов............................................................................................105

§2. Решение одномерной задачи Стефана методом дифференциальных рядов.......................110

§3. Краевые задачи стефановского типа со свободной границей..............................................114

§4. Численные эксперименты........................................................................................................118

Выводы к главе 5............................................................................................................................122

ГЛАВА 6. Проблема теплового удара и динамическая термоупругость

Введение..........................................................................................................................................124

§ 1. Проблема теплового удара и математическая формулировка задач динамической

темроупругости........................................................................................................................126

§2. Задача темоупругости для упругого полупространства [/,со).............................................131

§3. Задача динамической темоупругости для упругого полупространства с равномерно

движущейся границей z>l + vt, t> 0, v = const ................................................................135

§4. Компоненты ctxx(z, Fo) и <5yy(z,Fo) тензора напряжений при температурном нагреве

для упругих полупространств z>l и z>l + vt, t> 0, v = const ......................................143

§5. Задача термоупругости для упругих полупространств z>l + vt, t> 0, v = const и

z > ß~Jt, t > 0, ß = const в квазистатическом случае........................................................150

Выводы к главе 6............................................................................................................................157

ВЫВОДЫ К ДИССЕРТАЦИИ...............................................................................................159

Литература...............................................................................................................................161

Введение к работе

Краевые задачи для уравнений в областях с движущихся во времени границами представляют неисчерпаемое поле для исследований. С годами их поток не уменьшается, охватывая все новые математические объекты и все большее число самых разнообразных приложений. Помимо подробных обзоров прошлых лет [372] по данной тематике, в настоящей диссертации приводится обзор последних работ в этом направлении.

Решению задач в нецилиндрических областях, как плоских, так и многомерных, по-прежнему, уделяется большое внимание как со стороны российских, так и зарубежных исследователей [30], [66]—[69], [71], [94], [96]—[97], [107], [157], [181], [207], [230]—[231], [249], [251]—[252], [258], [268], [292], [299]—[300], [303]—[305], [307], [318]—[319], [324], [326], [329], [333], [335]—[336], [338], [342]—[346], [349], [351], [354]—[356], [361]. Здесь стоит отметить ряд работ [35], [347], посвященных построению комплексов программ и численных методов, для решения задач в нецилиндрических областях. Большое количество приложений в экологии [26], [93], [294], криомедицине [28], задачах кинетического переохлаждения [48], фильтрации [59], теплопереноса в процессе кристаллизации [145], о промерзании грунта [331], теории горения [311], биологии [302], медицине [307], [310], физиологии [321], делают решение задач по данной проблематике одной из важнейших целей современной прикладной математики. В работах [65], [70], [230], [235]—[237] рассматриваются вопросы динамического напряжения упругой нецилиндрической области.

Важное место занимают работы, посвященные качественной теории по данной проблематике [106], [136], [166], [252], [293]. Среди них можно выделить работы, рассматривающие вопросы существования и единственности решения одномерных и многомерных задач, [104], [309], [312], [330], [341]; устойчивости и асимптотического поведения решения задач [13], [195]—[196], [211], [310], [319]; устойчивости и

асимптотического поведения подвижной границы или поверхности [79], [199]; корректности постановки задачи [105]; регулярности решений и свободной границы [17]—[18], [323].

Нельзя не отметить работы, посвященные развитию новых методов решения задач в нецилиндрических областях [233], [252]—[253], таких как метод функции Грина [67], [268] и получение новых интегральных соотношений [126].

Как и прежде, большое внимание уделяется задачам о фазовых превращения (так называемой проблеме Стефана) и задачам типа Стефана [24], [27]—[28], [33], [48], [53], [75], [85], [108], [195], [200], [206], [215], [229], [288]—[290], [293], [295], [297], [306], [317], [322], [325], [332]—[334], [336], [348], [350], [359]—[360]; обратным задачам Стефана [54]—[56]. Уделяется внимание задаче Стефана с движущимся источником холода [84], а также исследуется асимптотическое поведение решения задачи о фазовом переходе и свободной границы [352].

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

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

закону = рл// • На примере задачи Стефана показывается, этот закон является очень важным для практики и наблюдается во многих природных явлениях.

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

ГЛАВА 1. Современные аналитические методы решения краевых задач нестационарной теплопроводности в областях нецилиндрического типа

Введение

Первые результаты при решении краевых задач с движущимися границами были получены для частных случаев задания краевых условий (нулевые или постоянные) за счет применения искусственных приемов или удачных находок при проведении замены переменных или тождественных преобразований [128], при этом рассматривались простейшие законы движения границы в сочетании с узкими граничными условиями.

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

В 70-е был разработан аналитический метод решения краевых задач нестационарной теплопроводности в областях с движущимися границами, основанный на специальном функциональном преобразовании, переводящем данную краевую задачу в подвижную систему координат таким образом, что преобразованное уравнение теплопроводности допускало точное аналитическое решение [376]. Это позволило расширить класс решаемых задач.

Параллельно с развитием практических методов решения краевых задач нестационарной

теплопроводности в областях с движущимися границами, шло исследование теоретических проблем, связанных с условиями существования, единственности и устойчивости решений этих задач. Эти задачи рассматривались в рамках качественной теории дифференциальных уравнений параболического типа в нецилиндрических областях фазового пространства [109], [110].

С математической точки зрения краевые задачи для уравнений параболического типа в области с движущейся границей принципиально отличны от классических (для областей цилиндрического вида). Вследствие зависимости границы области от времени к этому классу задач в общем случае не применимы методы разделения переменных и интегральных преобразований Фурье-Хенкеля-Лапласа [86], [117]—[120], [125], [221], так как, оставаясь в рамках классических методов математической физики [12], [14], [29], [44], [50], [90], [140], [150], [153], [173], [183], [198], [208], [219], [250], [255], [259], [267], не удается согласовать решение уравнения теплопроводности с движением границы области теплопереноса. Развитие этой проблемы шло, по-видимому, следующим образом. С одной стороны, точные решения задач подобного типа удалось получить при помощи удачных догадок, искусственных приемов, причем для весьма ограниченного числа случаев движения границы (вначале линейного в

области x>l + vt, t> 0, затем параболического х > ß-ft, t> 0, [128]) и для частного вида граничных условий: постоянных и первого рода. С другой стороны, на этих задачах при весьма их общей постановке отрабатывались классические методы решения краевых задач для дифференциальных уравнений математической физики (и их модификации) [115]: тепловых потенциалов; контурного интегрирования; продолжений; функций Грина; вариационные; разложения искомой функции в ряд по обобщенным степеням; обобщение конечных интегральных преобразований на нецилиндрические области; функциональных преобразований Гринберга, а также методы, основанные на использовании интегральных, интегро-дифферен-

циальных или обыкновенных дифференциальных уравнений, разностные, асимптотические и численные [21], [45], [64], [81], [98]—[99], [101], [169]—[171], [174]—[176], [184]—[185], [187], [190]—[191], [194], [197], [214], [222]—[223], [243]—[244], [247]—[248], [269]—[271]. Объясним был и тот факт, что решения одного и того же класса задач применялись различные подходы. Это объясняется тем, что решение одной и той же тепловой задачи можно искать в различных классах функций, определяемых аналитическим подходом при решении задачи. Эти функции должны быть таковыми, чтобы они, во-первых, достаточно легко находились и, во-вторых, обеспечивали сходимость процесса настолько хорошо, чтобы можно было сделать требуемые в задаче заключения о свойствах полученного решения. Представление аналитического решения задачи в эквивалентных функциональных формах (тождественных в смысле числа) имеет большую практическую ценность, так как позволяет варьировать решением в зависимости от постановки задачи: например, решение в форме типа ряда Фурье, удобной для больших времен (находится интегральным преобразованием Фурье), или в виде формулы суммирования Пуассона, более подходящей для малых времен (находится операционным методом). Для областей с движущимися границами это обстоятельство имеет особенно важное значение, учитывая широкое применение этого класса задач. В то же время следует подчеркнуть, что несмотря на хорошо развитую аналитическую теорию нестационарного тепло- массопереноса и близких направлений [23], [58], [72], [86]—[87], [112], [117]—[120], [125], [133], [158]—[161], [186], [210], [221], [277], [74], [88]—[89], [102], [131], [134]—[135], [138]—[139], [143], [146]—[148], [154]—[155], [202], [205], [213], [216], [218], [239], [272], [274], [278], [280]—[281], [287], достигнутые за последние два десятилетия успехи в нахождении точных аналитических решений задач для различных законов движения границы весьма незначительны. По-видимому, одной из причин такого положения являются технические трудности вычислительного характера при нахождении точного аналитического

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

§1. Постановка краевых задач нестационарной теплопроводности

Пусть Q, — нецилиндрическая область в фазовом пространстве (п +l) измерений, сечение которой плоскостью-характеристикой t = const > /0 > 0, есть выпуклая область д(д е Л") изменения М(х1,х2,...,хп), S, — кусочно-гладкая поверхность, зависящая от времени t > 0 и о�