Методы математического моделирования при численно-аналитических решениях задач для уравнений волнового теплопереноса тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат физико-математических наук Селин, Илья Александрович

• Селин, Илья Александрович
• кандидат физико-математических наук
• 2010

Оглавление диссертации кандидат физико-математических наук Селин, Илья Александрович

Введение

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

1.1. Существующие методы построения математических моделей физических процессов переноса потенциала при высокоинтенсивном воздействии.

1.2. Новый метод построения математических моделей на основе представления законов переноса потенциала в виде ряда по малому параметру.

1.3. Математическое моделирование процессов переноса при высокоинтенсивных воздействиях на примере волнового теплопереноса.

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

2.1. Теорема об эквивалентности описания волнового переноса на основе уравнений гиперболического типа и квазилинейных уравнений параболического типа.

2.2. Метод аналитического исследования явлений переноса в нелинейных анизотропных средах на основе квазилинейных уравнений' параболического типа.

2.3. Теорема об условии возникновения и распространения возмущений в нелинейных анизотропных средах в виде бегущих волн.

2.4. Численное моделирование распространения и возникновения ударных волн в нелинейных средах на примере волн теплопереноса.

3. Численное моделирование волнового теплопереноса в анизотропной среде.

3.1. Интенсивный температурный нагрев анизотропного тела.

3.1.1. Постановка задачи.

3.1.2. Конечно-разностная аппроксимация и метод численного решения.

3.1.3. Исследование аппроксимации и устойчивости конечноразностной схемы метода переменных направлений с экстраполяцией.

3.1.4. Результаты численного исследования.

3.2. Интенсивный нагрев анизотропного тела тепловым потоком.

3.2.1. Постановка задачи.

3.2.2. Конечно-разностная аппроксимация и метод численного решения.

3.2.3. Результаты численного исследования.

4. Комплекс программ.

4.1. Среда разработки и назначение программного комплекса.

4.2. Структура программного комплекса.

4.3. Описание интерфейса.

Введение диссертации (часть автореферата) на тему «Методы математического моделирования при численно-аналитических решениях задач для уравнений волнового теплопереноса»

При исследовании процессов переноса потенциала при высокоинтенсивном воздействии (при взаимодействии лазерного излучения с твёрдой поверхностью, входе гиперзвуковых летательных аппаратов в плотные слои атмосферы, сверхзвуковых контактных взаимодействиях и т.п.) гипотезы о пропорциональности плотности потока вектору градиента потенциала, построенные на феноменологических представлениях, приводят к бесконечной скорости распространения возмущений, что противоречит фундаментальным законам естествознания. К законам, построенным на данной теории относятся следующие: х,у,2,1) = -Х§гайТ{х,у,2,{) - закон Фурье, Я - коэффициент теплопроводности, плотность теплового потока, Г-температура;

0:,y,z,t') = -DgradC(<x,y,z,t) — закон Фика, И-коэффициент диффузии, ц - плотность потока диффузии, С - концентрация;

7О,.у,г,0 = -К%гас1Н{х,у,2,Г)- закон Дарси, К -коэффициент фильтрации, д- поток объёмного расхода или скорость фильтрации, Я- давление.

Все они имеют одинаковую форму, и в общем виде их можно записать так ¿7 О, .у, г,/) = -AgradU{x,y,z,t) — обобщённый закон переноса, где А - коэффициент переноса, д- плотность потока, и - потенциал. Дифференциальное уравнение переноса, получаемое из данного обобщённого закона переноса в одномерном случае имеет вид ди д2и

-= а——. д( дх2

Достаточно продифференцировать фундаментальное решение данного нестационарного уравнения переноса х,0 = -Д=ехр л/4 ла/

Г ^ ч и устремить время к нулю, откуда будет видно, что скорость переноса потенциала в начальный момент времени рано бесконечности. Данный парадокс возникает из-за специфического рассмотрения твёрдой среды. Законы переноса получены из рассмотрения тела, как сплошной среды и, как следствие, предположения о диффузионном распространении потенциала. Однако, при рассмотрении среды на молекулярном уровне, её нельзя считать сплошной, так как между молекулами находится свободное пространство. Приближение сплошной среды, подразумевающее у неё отсутствие внутренней структуры, означает, что в интегральных законах сохранения для этой среды можно совершать предельный переход при стремлении объёма к нулю. Такой предельный переход позволяет получить уравнение сохранения энергии в дифференциальной форме. С физической точки зрения эта процедура некорректна, так как среда всегда состоит из отдельных элементов и имеет свою внутреннюю дискретную структуру [82].

Для устранения парадокса о бесконечной скорости распространения возмущений в рамках теории теплопроводности К. Каттанео в 1958 году из молекулярно-кинетических представлений, используя гипотезу о конечности продолжительности удара молекул и представления о длине свободного пробега молекул, получил новый закон теплопроводности, в котором да появилось дополнительное слагаемое г—, учитывающее дискретную дt молекулярную структуру среды и отвечающее за инерционность тепла [88].

В этом слагаемом т — время релаксации, время установления термодинамического равновесия между тепловым потоком и градиентом температуры. Очевидно, что обобщённый закон переноса записать в следующем виде т^~- + q = -AgradU. (1)

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

При решении дифференциального уравнения, полученного из этого закона, наблюдается разрыв первого рода потенциала, распространяющийся от источника. Таким образом, закон (1) описывает возникновение волн при высокоинтенсивном воздействии, который приводит к локальной неравновесности системы. Эффекты локальной неравновесности чаще всего наблюдаются при воздействии на тело короткими импульсами энергии, в ударных волнах, при низких температурах. В этих случаях время релаксации системы к локальному равновесию сравнимо со временем самого процесса [67].

Локальное равновесие справедливо для моментов времени, значительно превышающих время релаксации.

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

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

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

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

Работ, посвященных волновому переносу, мало. Все их можно разделить на 2 группы. Первая представляет собой исследование волнового переноса на основе гиперболического уравнения, а ко второй группе относятся те работы, в которых волновой перенос описывается квазилинейным уравнением. Большинство работ, относящихся к первой группе, хоть и затрагивают волновой перенос на основе закона с учётом времени релаксации, но сами тепловые волны не исследуют, а решают прикладные задачи [1, 4, 5, 36, 41, 68, 81].

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

В [81] показано, что при воздействии лазерного излучения возникающие тепловые волны влияют на прочностные характеристики тел. Расчёт температурных полей по закону Фурье может давать критические ошибки для теплового напряжения, что может привести к разрушению нагреваемого тела.

В работе [1] рассматривается нагрев пластины импульсным лазерным излучением. Используется гиперболическое уравнение теплопроводности. Показано, что в области малых частот (m>«1) излучения получаются зависимости, согласующиеся с теорией Фурье. В пределе больших частот (tg>»1) гиперболическое уравнение теплопроводности описывает процесс, при котором температурные возмущения распространяются с конечной скоростью. Что показывает существенную неравновесность процесса, в котором необходимо учитывать появление тепловых волн.

Среди авторов, занимающихся исследованиями тепловых волн на основе закона Максвелла-Каттанео можно выделить Бубнова В.А., Корнеева A.C., Рубину Л.И., Соболева С.Л., Шашкова А.Г.

В работах, посвященных данной тематике рассматривается в основном одномерная область без учёта граничных условий [14, 16; 18, 23].

В [49] показано, что при высокоинтенсивном теплообмене в результате воздействия на тело теплового потока температура в начальный момент времени, рассчитанная по закону Максвелла-Каттанео оказывается больше, чем температура, рассчитанная по закону Фурье, на несколько сотен градусов.

Наиболее полной представляется работа [82]. Но в данной работе, не смотря на всесторонние исследования не показаны сами тепловые волны и нет аналитических решений краевых задач.

Нет исследований в нелинейных средах на основе гиперболического уравнения переноса; Работы по исследованию волнового переноса в анизотропных пространствах отсутствуют.

Исследованию процессов волнового теплопереноса на основе квазилинейного уравнения теплопроводности посвящены работы Зельдовича Я.Б., Баренблатта Г.И., Калашникова A.C., Ладыженской O.A., Самарского A.A. В работах Зельдовича было установлено, что скорость распространения тепла для процессов нелинейной теплопроводности является конечной, в отличие от задач, описываемых линейным уравнением теплопроводности и допускается существование разрывов производных на фронте тепловой волны.

В [8] предложен новый метод решения подобных задач.

В [28, 40, 45, 46, 80] разобраны общие вопросы существования и дифференцируемости решений, асимптотика, оценки.

Очень подробно, с большим количеством примеров волновой теплоперенос описан в [61, 63]. В [62] предлагаются численные методы для расчёта тепловых волн.

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

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

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

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

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

4. Модифицирован существующий численный метод для моделирования волнового теплопереноса в анизотропном пространстве.

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

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

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

В параграфе 1.2 описывается новый метод математического моделирования волнового переноса на основе разложения в ряд по малому параметру. Для полученного ^уравнения, ставятся' начально-краевые задачи.

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

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

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

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

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

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

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

Глава 3 посвящена численному моделированию теплопереноса на основе уравнения волнового теплопереноса в анизотропной среде.

В параграфе 3.1 численно решается первая краевая задача для уравнения волнового теплопереноса. Описана конечно-разностная схема, исследована её устойчивость и аппроксимация.

В параграфе 3.2 численно решается вторая краевая задача для уравнения волнового теплопереноса. Описывается модификация конечноразностной схемы с помощью интегро-интреполяционного метода A.A. Самарского.

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

В главе 4 описан программный комплекс, созданный для решения задач теплопроводности на основе методов математического моделирования, описанных в диссертации, в котором реализованы рассмотренные численные методы.

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

Классическая теория переноса — на основе градиентных законов — справедлива, если характерная скорость данного процесса много меньше скорости распространения возмущений в среде. При высокоинтенсивных процессах, классическая теория процессов переноса становится несправедлива и следует пользоваться локально-неравновесными методами описания таких систем [67]. В теории теплопроводности закон, представляющий взаимосвязь между тепловым потоком и градиентом температуры, предложенный сначала П. Верноттом, потом Каттанео и A.B. Лыковым, относится к одному из таких методов. Данный закон в общей форме является общим законом волнового переноса потенциала, которым необходимо пользоваться при< высокоинтенсивных воздействиях. Высокоинтенсивное воздействие подразумевает сообщение телу большого возмущения потенциала за промежуток времени t 0. При таком воздействии, поле потенциала, рассчитанное по параболическому уравнению переноса, в первые моменты времени не является достоверным. Кроме того, модель Вернотта-Лыкова учитывает конечную скорость распространения тепла [7], соответственно скорость получается конечной и при переноса любого потенциала, перенос которого построен на аналогичном законе.

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

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Селин, Илья Александрович

ЗАКЛЮЧЕНИЕ

На основе проведённых исследований в диссертационной работе получены следующие научные результаты:

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

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

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

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

5. Разработан программный комплекс, позволяющий решать задачи теплопереноса при высокоинтенсивном нагреве с учётом волнового теплопереноса.

Список литературы диссертационного исследования кандидат физико-математических наук Селин, Илья Александрович, 2010 год

1. Александров А.Н., Голубев Е.В. Нагрев бесконечной металлической пластины импульсным лазерным излучением// Вестник ЮУрГУ. Серия Математика, физика, химия. 2005. Вып. 5. № 2 (42).

2. Алексашенко A.A. Аналитическое исследование тепло и массопереноса с учётом конечной скорости переноса. Канд. дис. ИТМО. Минск, 1968.

3. Алексашенко A.A., Алексашенко В.А., Селезнёв Н.В. Решение уравнений тепло- и массопереноса для тел для тел классической формы с учётом конечной скорости капиллярного движения. В кн.: Строительная теплофизика. M.-JL: Энергия, 1966.

4. Анисимов С.И., Имас Я.А., Романов Г.С., Ходыко Ю.В. Действие излучения большой мощности на металлы. М.: Наука, 1970.

5. Аттетков A.B., Волков И.К., Тверская Е.С. Оптимальная толщина охлаждаемой стенки с покрытием, подверженной локальному импульсно-периодическому нагреву. ИФЖ, 2001. Т. 74, №6.

6. Баренблатт Г.И., Виши М.И. О конечной скорости распространения в задачах нестационарной фильтрации жидкости и газа. Прикладная математика и механика, 1956. Т. 20. Вып. 3.

7. Баумейстер К., Хамилл Т. Гиперболическое уравнение теплопроводности. Решение задачи о полу бесконечном теле.— Теплопередача, 1969, №4.

8. Баутин С.П. Аналитическая тепловая волна. М.: Физматлит, 2003.

9. Бахвалов Н.С. Численные методы. М.: Наука, 1975.

10. Бейтмен Г., Эрдейи А. Высшие трансцендентные функции. М.: Наука, 1973.

11. Беккер Р. Теория теплоты. М.: Энергия, 1974.

12. Березин И.С., Жидков Н.П. Методы вычислений. М.: Гос. издат. физико-математической лит-ры, 1960.

13. Берман Р. Теплопроводность твёрдых тел. М.: Мир, 1979.

14. Бубнов В.А. О тепловых волнах. Теплофизика высоких температур, 1982, Т.20, №5.

15. Бубнов В.А. Замечания к волновым решения нелинейного уравнения теплопроводности. Инженерно-физический журнал, 1981, Т. XL, №5.

16. Бубнов В. А. О характере теплообмена в акустической волне. ИФЖ. 1976 Т. 31, №3.

17. Бураков В.А., Берцун В.Н., Крицкий O.JI. Сравнительный анализ численных методов решения нестационарной задачи анизотропной теплопроводности. Томск: Пеленг, 2001.

18. Волосевич П.П., Леванов Е.И. Автомодельные решения задач газовой динамики и теплопереноса. М.: МФТИ, 1997.

19. Годунов С.К., Рябенький B.C. Разностные схемы (введение в теорию). М.: Наука, 1977.

20. Гуревич А. В., Минц Р. Г. Тепловые автоволны в нормальных металлах и сверх проводниках.— М.: ИВТАН СССР, 1987.

21. Двайт Г.Б. Таблицы интегралов. М.: Наука, 2005.

22. Дёч Г. Руководство к практическому применению преобразования Лапласа и Z-преобразования. М.: Наука, 1971.

23. Динариев О. Ю., Николаев О. В. О релаксационных процессах в низкопроницаемых пористых материалах. ИФЖ. 1990. Т. 58.

24. Жданов И.С. Применение функций Грина к решению задач математической физики. М.: Наука, 1995.

25. Зарубин B.C. Инженерные методы решения задач теплопроводности. М.: Энергоатомиздат, 1983.

26. Зарубин B.C., Кувыркин Г.Н. Математические модели термомеханики. М.: Физматлит, 2002.

27. Зарубин B.C., Кувыркин Г.Н. Математические модели теплопроводности. Тр. 2-ой российской национальной конференции по теплообмену. М.,1998.

28. Зельдович Я. Б., Баренблатт Г. И., Либрович В. Б., Махвиладзе Г. М. Математическая теория горения и взрыва.— М.: Наука, 1980.

29. Зельдович Я.Б., Райзер Ю.П. Физика ударных волн и высокотемпературных гидродинамических явлений. М.: Наука, 1966.

30. Змитренко Н.В., Михайлов А.П. Инерция тепла. М.: Знание, 1982.

31. Калашников A.C. Об уравнениях типа нестационарной фильтрации с конечной скоростью распространения возмущений. Вестник МГУ. Сер 1. Математика, механика. 1971, №6.

32. Калашников A.C. О характере распространения возмущений в задачах нелинейной теплопроводности с поглощением. Журнал вычислительной математики и математической физики. 1974. Т 14, №4.

33. Камья Ф.М. Импульсная теория теплопроводности. М.: Энергия, 1972.

34. Карслоу Г, Егер Д. Теплопроводность твёрдых тел. М.: Наука, 1964.

35. Карташов Э.М. Аналитические методы в теории теплопроводности твёрдых тел. М.: Высшая школа, 2001.

36. Кудинов В.М., Вовк A.A., Даниленко В.А., Макаренко A.C. Влияние процессов релаксации на горение и взрыв. Киев: Институт геофизики, 1983.

37. Колесников П.М. Простые и ударные волны при нелинейном высокоинтенсивном нестационарном процессе тепло массо переноса, ИФЖ 1968, Т.15, №3.

38. Корнеев А. Гиперболическое уравнение теплопроводности // Известия РАН. Энергетика. 2001. - № 4.

39. Краевые задачи теории теплопроводности // Сб. статей института математики УССР. Киев. 1975.

40. Крылов Н.В. Нелинейные эллиптические и параболические уравнения второго порядка. М.: Наука, 1985.

41. Крылович В.И., Дербан В.И. Термическая генерация упругих колебаний с учетом конечной скорости распространения тепла ИФЖ 1975. Т. 29, №3.

42. Кутателадзе С.С. Основы теории теплообмена. М.: Атомиздат, 1979.

43. Кутателадзе С.С, Накоряков В.Е. Тепломассообмен и волны вгазожидкостных системах. Новосибирск: Наука, 1984.

44. Кузнецова Е.Л., Селин И.А., Формалев В.Ф. Задача типа Стефана с произвольным числом подвижных границ фазовых превращений. Вестник МГТУ им. Н.Э.Баумана. 2010, №2 (37), с.49-58.

45. Ладыженская O.A., Солонников В.А., Уральцева H.H. Линейные и квазилинейные уравнения параболического типа. М.: Наука, 1967.

46. Леванов Е.И., Сотский E.H. // Математическое моделирование. Нелинейные дифференциальные уравнения математической физики.— М.: Наука, 1987.

47. Леванов Е.И., Сотский E.H. Некоторые свойства процесса теплопереноса в неподвижной среде с учётом релаксации теплового потока. Инженерно-физическй журнал. 1986. Т1, №6.

48. Леванов Е.И., Сотский E.H. Теплоперенос с учётом релаксации теплового потока. Математическое моделирование (нелин диф ур мат физ) 1987.

49. Лихт М.К. О распространении возмущений в задачах, связанных с вырождающимися квазилинейными уравнениями параболического типа. Дифференциальные уравнения, 1966. Т 2, №7.

50. Лыков A.B. Теория теплопроводности. М.: Высшая школа, 1966.

51. Лыков A.B. Тепломассообмен: (Справочник). М.: Энергия, 1978.

52. Маурер М.Ж., Томпсон Х.А. Эффекты отклонения от модели Фурье при высоких тепловых потоках. Теплопередача, 1973, №2.

53. Морс Ф, Теплофизика М : Наука, 1968.

54. Пирумов У.Г. Численные методы. М.: Вагриус, 2004.

55. Постригач Я.С., Коляно Ю.М., Обобщенная термомеханика, 1976.

56. Полянин А.Д. Справочник по линейным уравнениям математической физики. М.: Физико-математическая литература, 2001.

57. Рубина Л. И. О распространении слабых разрывов для системы уравнений магнитной газодинамики // Прикладная математика и механика. Т. 33. Вып. 5.

58. Рубина Л.И. О затухании и разрушении слабых разрывов, распространяющихся по области центрированных волн // Численные методы механики сплошной среды. 1972.1. Т. 3, № 4.

59. Рубина Л. И. О распространении слабых разрывов для квазилинейных систем // Прикладная математика и механика. 1972. Т. 36. Вып. 3.

60. Самарский A.A. Теория разностных схем. М.: Наука, 1983.

61. Самарский A.A., Галактионов В.А., Курдюмов С.П., Михайлов А.П. Режимы с обострением в задачах для квазилинейных параболических уравнений. М.: Наука, 1987.

62. Самарский A.A., Соболь И.М. Примеры численного расчёта температурных волн. Журнал вычислительной математики и математической физики, 1963, Т.З, №4.

63. Самарский А. А., Соболь И. М. Ц ЖВМ и МФ. 1963. Вып. 3, № 4.

64. Самарский A.A., Гулин A.B. Устойчивость разностных схем. М.: Наука, 1973.

65. Селин И.А. Моделирование тепловых волн в анизотропных твёрдых средах. Тезисы докладов восьмой международной конференции по неравновесным процессам в соплах и струях. 2010.

66. Скородинский В. А. Тепловой удар на поверхности полупространства с учетом конечной скорости распространения тепла.— Некоторые вопросы прикладной математики, 5. Киев, Изд-во Ин-та математики АН УССР, 1971.

67. Соболев C.J1. Процессы переноса и бегущие волны в локально-неравновесных системах. Успехи физических наук, 1991, Т. 161, №3.

68. Ступоченко Е. В., Лосев С. А., Осипов А. И. Релаксационные процессы в ударных волнах.— М.: Наука, 1965.

69. Титов С. С. О движении фронта нелинейной диффузии // Прикладная механика и техническая физика. 1996. Т. 37, №4.

70. Формалев В.Ф. Численное исследование двумерных нелинейных задач теплопроводности в анизотропных телах. ТВТ. 1988. Т.26. №6.

71. Формалев В.Ф. Метод переменных направлений с экстраполяцией по времени для параболических задач со мешанными производными. Вычислительные технологии. 1996. Т.1. №2.

72. Формалев В.Ф., Ревизников Д. Л. Численные методы. М.: ФИЗМАТ ЛИТ, 2004.

73. Формалёв В.Ф., Селин И.А. О тепловых волнах в нелинейных твёрдых средах// В тр. 16-ой межд. конф. по вычислительной механике и современным прикладным программным системам, 2009 г. с.732.

74. Формалёв В.Ф., Селин И.А. О тепловых волнах в твёрдых телах// В тр. XV Межд. симпозиума «Динамические и технологические проблемы механики конструкций и сплошных сред» им. А.Г.Горшкова. 2009 г, Т.1. сЛ 57-158.

75. Формалёв В.Ф., Селин И.А., Кузнецова Е.Л. Возникновение и распространение тепловых волн в нелинейном анизотропном пространстве. Известия РАН «Энергетика», 2010, №3, с. 136-141.

76. Формалёв В.Ф., Селин И.А. Тепловые волны в нелинейном анизотропном пространстве. Труды МАИ, 2010, №37, с. 12-13.

77. Формалёв В.Ф., Селин И.А. Моделирование тепловых волн в нелинейном анизотропном пространстве. Вестник Самарского Государственного Технического Университета, 2010, №1 (20), с.239-244.

78. Формалёв В.Ф., Селин И.А., Колесник С.А. Анализ тепловых волн в анизотропных пространствах. // В тр. 5-ой Российской национальной конференции по теплообмену. 2010.

79. Фридман А., Уравнения с частными производными параболического типа. М.: Мир, 1968.

80. Харитонов В.В., Субботин В.И., Гришутин П.А., Тимонин A.C.

81. Динамические эффекты при импульсном нагреве лазерных зеркал. Теплофизика высоких температур, 1983, Т.21, №6.

82. Шашков А.Г., Бубнов A.B., Яновский С.Ю. Волновые явления теплопроводности. М.: УРСС, 2004.

83. Шашков А. Г. Системно-структурный анализ процессов теплообмена и его применением.: Энергоатомиздат, 1983.

84. Штер И. М.Плоские гармонические волны в неограниченной термоупругой среде с конечной скоростью распространения тепла ИФЖ, 1973, Т.24, №4.

85. Шабловский О. Н. К исследованию нелинейных задач высокоинтенсивного нестационарного теплопереноса. ИФЖ.1987 Т. 52, № 2.

86. Шабловский О. Н. О нелинейных задачах плавления и испарения материалов под действием интенсивных потоков энергии с учетом тепловой релаксации. ИФЖ. 1988. Т. 55, № 3.

87. Яненко H.H. Метод дробных шагов решения многомерных задач математической физики. Новосибирск: Наука, 1967.

88. Cattaneo С. Comptes Rendus. 1958. Vol. 247, N 4.

89. Maxwell J. С., Philos. Trans. Roy. Soc. London, 1867.

90. Peletier L.A. A necessary and sufficient conditions for the existence of an interface in flows through porous media // Arch. Ration. Mech. Anal.1974 . V.56.