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

• Зайнуллин, Рифат Гильметдинович
• кандидат физико-математических наук
• 2010

Оглавление диссертации кандидат физико-математических наук Зайнуллин, Рифат Гильметдинович

ВВЕДЕНИЕ.

Анализ работ по теме исследования.

1. Метод тепловых потенциалов.

2. Метод функций Грина.

3. Метод обобщенных рядов.

4. Метод функциональных преобразований.

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

6. Метод автомодельных преобразований.

7. Другие методы.

Глава 1. ОДНОМЕРНАЯ ДВУХФАЗНАЯ ЗАДАЧА ТЕПЛООБМЕНА

С ФАЗОВЫМ ПЕРЕХОДОМ С ОГРАНИЧЕННОЙ

ЗОНОЙ ТЕПЛОВОГО ВЛИЯНИЯ.

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

2. Нахождение собственных функций и закона движения свободной границы.

3. Вывод формулы для нормирующих множителей.

4. Определение собственных значений.

5. Построение температурного поля.

6. Анализ полученного решения.

Глава 2. ДВУМЕРНАЯ ОДНОФАЗНАЯ ЗАДАЧА ТЕПЛООБМЕНА

С ФАЗОВЫМ ПЕРЕХОДОМ.

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

2. Нахождение собственных функций по первой пространственной координате и закона движения свободной границы на горизонтальной оси.

3. Определение собственных значений по первой пространственной координате.

4. Нахождение собственных функций по второй пространственной координате и закона движения свободной границы.

5. Определение собственных значений по второй пространственной координате.

6. Построение температурного поля.

7. Анализ полученного решения.

Глава 3. РАЗРАБОТКА МЕТОДИКИ ЧИСЛЕННОГО РАСЧЕТА АНАЛИТИЧЕСКОГО РЕШЕНИЯ ОДНОМЕРНОЙ ЗАДАЧИ, ПОЛУЧЕННОГО МЕТОДОМ РАЗЛОЖЕНИЯ ПО СОБСТВЕННЫМ ФУНКЦИЯМ САМОСОПРЯЖЕННОГО ДИФФЕРЕНЦИАЛЬ

НОГО ОПЕРАТОРА.

1. Вывод вычислительных формул и оценка решения.

2. Построение графиков.

Выводы и заключения.

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

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

В работе изучаются одномерная двухфазная»задача Стефана с границами, перемещающимися по идентичному закону (с точностью до постоянного множителя), и двумерная однофазная задача Стефана. Задачи преобразуются к области с неподвижными границами с неоднородной правой частью и переменными коэффициентами в исходном уравнении, но с однородными краевыми условиями, и далее, для решения преобразованных задач строится конечное интегральное преобразование по каждой пространственной координате с неизвестным ядром, нахождение которого связано с постановкой и решением соответствующей спектральной задачи; причем собственные функции выражаются через вырожденные гипергеометрические функции (ВГГФ), а собственные числа в зависимости от величины введенного безразмерного параметра находятся тремя различными способами, два из которых асимптотические, а два — сочетаются с численными методами. При слишком больших значениях названного параметра асимптотика может быть использована повторно. Новый подход в развитии интегральных преобразований в конечном счете приводит к получению формулы обращения, что позволяет выписать приближенно-аналитическое решение задачи.

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

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

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

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

Поставленная цель достигается в результате решения следующих задач:

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

2. Доведение полученного приближенно-аналитического решения до компактных вычислительных формул с хорошо сходящимися рядами с помощью функциональных и рекуррентных соотношений и асимптотических представлений для* вырожденных гипергеометрических функций.

3. Оценка сходимости алгоритма численного расчета.

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

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

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

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

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

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

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

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

Апробация работы и* публикации. Результаты работы опубликованы в 15 научных статьях (2 из них в изданиях, рекомендованном ВАК, 4 в материалах конференций). Вклад соавторов публикации равнозначен. Список публикации по теме диссертации приведен в конце автореферата. Основные результаты диссертационной работы представлены и обсуждались:

- на конференции молодых учёных БФАН СССР "Исследования по математике, физике, механике и процессам управления" (Уфа, 1987г.);

- на IV Уральской региональной конференции "Функционально-дифференциальные уравнения и их приложения" (Уфа, 1989г.);

- на III Сибирском конгрессе по прикладной и индустриальной математике (Новосибирск, 1998г.); на-международной научной конференции "Моделирование, вычисление, проектирование в условиях неопределенности" (Уфа, 2000г.); на научно-практической конференции, посвящённой году здоровья^ и 70-летию БГМУ (Уфа; 2002г.);

- на международной математической конференции, посвященной памяти А.Ф. Леонтьева.(Уфа, 2007г.);

- на международной научной конференции "Дифференциальные уравнения и сложные проблемы" (Стерлитамак, 2008г.);

- на III международной научной, конференции' "Современные проблемы прикладной' математики и математического моделирования" (Воронеж, 2009г.);

- на VI Всероссийской научной4 конференции с международным участием "Математическое моделирование и кривые задачи" (Самара, 2009г.);

- на II международной научной конференции "Математическое моделирование и дифференциальные уравнения" (Минск, 2009г.);

- на научных семинарах кафедры дифференциальных уравнений Уральского,. Башкирского, Куйбышевского, Саратовского, Казанского, Московского госуниверситетов, института математики с ВЦ- УНЦ РАН, Стерлитамакского филиала Уфимского БГПУ (1990-2000г.г.). Основные положения, выносимые на защиту:

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

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

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

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

Классическим примером проблемы, связанной с решением уравнения теплопроводности в области, граница которой или ее часть вследствие имеющего место фазового перехода перемещается по неизвестному ранее закону, является задача о плавлении твердого тела, которое в начальный момент времени находится при температуре фазового перехода и заполняет полупространство 0 < к < X, а в области 0 < х < к это же вещество находится в жидкой фазе, где имеет место некоторое начальное распределение температуры ф(х) > 0; ф(/г) = Тп - температура плавления. На неподвижной границе х = 0 задается тепловой поток или температура. Требуется найти распределение температуры в жидкой фазе и движение поверхности раздела фаз. Задачи подобного типа возникают при изучении различных процессов, связанных с фазовыми превращениями, таких как, например, плавление или затвердевание, кристаллизация, испарение, диффузионные процессы и другие.

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

Используемый метод решения задач для- уравнения теплопроводности с подвижной границей, который для краткости принято называть методом ВГГП, основан на идеях профессора М. Н. Шафеева, изложенных в,работах [92]-[95]. Этот метод, подобно методам интегральных преобразований по пространственным- переменным, основывается на возможности разложения решения как элемента некоторого гильбертова пространства по собственным функциям соответствующего краевой задаче-оператора Штурма-Лиувилля.

Теоретической основой разложения по собственным функциям некоторой спектральной задачи является спектральная теория линейных операторов в гильбертовом пространстве. Весьма общие результаты по- спектральному разложению дифференциальных операторов, получены в работах [55], [60]-[62] и др. В работах [73] путем предельного перехода даны основные теоремы решения в* регулярном» и сингулярном случаях, исследован спектр. Одним из естественных конструктивных подходов можно назвать метод Э: Ч. Титчмарша [79], [80].

Автор выражает большую благодарность и глубокую признательность доктору технических наук, профессору М. Н. Шафееву за постоянную поддержку и за постоянное внимание к данной- работе, докторам физико-математических наук, профессорам А. В. Жиберу за руководство работой над диссертацией, Г.Т. Булгаковой, а также признательность участникам семинаров кандидату физико-математических наук- А. М. Абдрахманову и др. за ряд полезных советов и рекомендаций.

АНАЛИЗ РАБОТ ПО ТЕМЕ ИССЛЕДОВАНИЯ

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

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

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

Точные решения задач подобного типа удавалось получить с помощью удачных догадок, искусственных приемов, причем для весьма ограниченного числа случаев движения границы (линейного или параболического) и для частного вида граничных условий (постоянных и I рода) [47]. С другой стороны, на этих задачах при весьма их общей постановке отрабатывались классические методы решения дифференциальных уравнений математической физики и их модификаций: тепловых потенциалов, контурного интегрирования, функций Грина, обобщенных рядов, дифференциальных рядов, продолжений, «мгновенных» собственных функций Гринберга, а также методы, основанные на использовании интегральных, интегродифференциальных или обыкновенных дифференциальных уравнений, асимптотические и численные

46]. Эти методьь будут рассмотрены ниже более подробно. При этом точные решения аналитическимшутем удалось получить лишь для простейших законов движения? границы, а именно,^ равномерного и параболического:. Всякие-попытки; получить» аналитическим путем точное решение: краевой, задачи-обобщенного типа; в области с границей, движущейся; по произвольному, закону, приводили к системе интегральных уравнений Вольтерра II рода [86], разрешить^ которую не удавалось, вследствие сложности ядер уравнения системы. Устанавливались лишь качественные результаты для такого рода системы,,доказывались существованиерешенияи его единственность [46];, [86], [108].

Решение этой проблемы значительно? продвинулось. вперед в 70-х годах после выхода в свет серии фундаментальных работ чл.-корр. АН СССР F. А. Гринберга с сотрудниками [14], [15], [48], [69];, [89]'. [90]; Им, было* получено функциональное преобразование, переводящее краевую задачу теплопроводности ; обобщенного* типа, в подвижную систему координат, в которой преобразованное: уравнение теплопроводности допускало бы точное, решение классическим- методомфазделения переменных для весьма широкого класса новых законов движения границы, при, соответствующих граничных условиях .

В дальнейшем^; [50] было получено точное решение первой обобщенной краевой задачи в, конечной области с границей, движущейся по» произвольному закону в декартовой, цилиндрической и; сферической (радиальный случай) системах координат. Для полуограниченной области точное решение удалось получить при любом виде граничных условий, включая и случай произвольной зависимости коэффициента внешнего теплообмена от времени в законе Ньютона на границе, движущейся также по произвольному закону.

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

23], [81]. Постановка рассматриваемых в- главе задач теплопроводности в общем случае состоит в следующем. Пусть П - нецилиндрическая область в фазовом пространстве (п 4- 1) измерений (п - координат и I - времени, t > 0), сечение которой плоскостью-характеристикой t — const > 0 представляет собой пространственную область Gt (конечную или частично ограниченную) п измерений с границей St, зависящей от времени t. Пусть п - внешняя нормаль к границе (направлена изнутри Gy), a P(xj- точка п-мерного пространства Gt. Требуется найти в области Г2 дважды непрерывно дифференцируемое по пространственным координатам (лг19.,л:и) и непрерывно дифференцируемое по времени t решение Т = уравнения ы п xjy J

Г-непрерывна изнутри в области GyQ (t = 0), Р eGt,t> 0,

Т(Р,0) = \[/(Р), Р G GtQ ,а на боковой поверхности St непрерывна, в каждой ее точке (St) имеет предельное значение нормальной производной и удовлетворяет граничному условию: а3 -cp(P,i), t > 0.

Набор а{ 2?з определяет вид краевой задачи.

В декартовой системе координат область (0^2 гДе Уг (0 ~ законы одного класса, во многих случаях может быть сведена к области [0,у(?)] путем введения подвижной системы координат z = x — y-[{t) с соответствующим изменением исходного типа задачи. Заметим также, что с помощью несложных преобразований переменное по координате начальное распределение температуры можно свести к нулевому с соответствующим

-а1ЯШ+агПР, о on изменением граничных функций [54]. При этом характер краевой задачи остается тем же.

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

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

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Зайнуллин, Рифат Гильметдинович

Выводы и заключения

1. Методом разложения по собственным функциям самосопряженного дифференциального оператора получены приближенные аналитические решения нестационарных одномерной и двумерной задач теплообмена с фазовым переходом. Метод основан на конечных интегральных преобразованиях с ядрами, нахождение которых связано с постановкой и решением соответствующих спектральных задач через ВГГФ. Установлено, что свободная граница движется по параболическому закону. Для одномерной задачи при данных краевых условиях и теплофизических характеристиках среды положение свободной границы определяется соотношением ^(т) =0.512950"3 7т в Международной системе СИ единиц физических величин. Значение параметра А, когда т = 0, совпадает с результатом работы Г. Гребера, С. Эрка, У. Григулль, опубликованный для двухфазных задач Стефана с однофазным начальном состоянием среды при условии, если температура в промерзшем слое меняется по линейному закону.

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

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

00 на бесконечности абсолютно сходятся не хуже числового ряда времени погрешность определения скорости движения свободной границы 1 ^ для полученного решения имеет порядок малости О —¡= .

-JtJ

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

Как видно из (3.1.2) и (3.1.6), нами получены компактные вычислительные формулы, и они содержат ряды, которые хорошо сходятся, что является очень удобным при построении температурных полей. Аналогично одномерному случаю тем же самым способом можно было бы построить температурное поле и для двухфазной задачи, но при этом объем излагаемого материала заметно увеличится. Числовые результаты для двухмерной задачи получаются аналогично одномерной задаче, только решения будут содержать двойные ряды.

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

1. Абрамович, М. Справочник по специальным функциям. / М. Абрамович, И. Стиган. М.: Наука, 1979. - 832 с.

2. Авдонин, Н.А. Математическое описание процессов кристаллизации. — Рига: Зинатне, 1980. 180 с.

3. Арсенин, В.Я. Методы математической физики и специальные функции. М.: Наука, 1984. - 384 с.

4. Ахиезер, Н. И. Теория линейных операторов в гильбертовом пространстве. -М.: Наука, Т.1, 1977.-415 с.

5. Базалий, Б. В. Задача Стефана для уравнения Лапласа с учетом кривизны свободной границы // Укр. мат. ж., 1997. -49. №10. - С. 1299-1315.

6. Базалий, Б. В. О регулярности решения задачи со свободной границей дляуравнения utко хх ' > 1/ Б. В. Базалий, Н. В. Краснощек // Алгебраи анализ, 2000. 12. - №2. - С. 100-130.

7. Бейтмен, Г. Высшие трансцендентные функции. / Г. Бейтмен, А. Эрдейн. -М.: Наука, ч.1. 1973. - 294 е.; ч.2. - 1974. - 295 с.

8. Березанский, Ю. М. Разложение по собственным функциям самосопряженных операторов. — Киев: Наукова думка, 1965. — 798 с.

9. Бижанова, Г. И. О задачах со свободными границами для параболических уравнений второго порядка / Г.И. Бижанова, В.А. Солонников // Алгебра и анализ, 2000. т. 12. - №6. - С. 98-139.

10. Бородкин, М. А. Двухфазная контактная задача Стефана // Укр. мат. ж., 1995.-47.-№2.-С. 158-167.

11. Бородкин, М. А. Существование глобального классического решения в некоторой нелинейной параболической задаче со свободными границами // Доп. Нац. АН Украиш, 1999Г- №6. С. 7-12.

12. Гребер, Г. Основы учения о теплообмене. / Г. Гребер, С. Эрк, У. Гршулль ИЛ. - М.: 1958. - 566 с.

13. Гринберг, Г. А. О некоторых точных решениях уравнения Фурье, для-изменяющихся со временем областей. / Г. А. Гринберг, В; А. Косс // Прикл.матем. и мех., т.35. —№3. 1971.-С. 759-760.

14. Гринберг, Г. А. О решении задач диффузионного типа- для расширяющихся или сжимающихся областей,. форма которых меняется со временем без соблюдения подобия. / Г. А. Гринбёрг, В*: А. Косс // Прикл.матем;,и мех., т.ЗЗ. №4. - 19691- С. 755-756.

15. Гринберг, Г. А. О движении- поверхности раздела фаз в задачах Стефанского типа. / Гринберг Г.А., Чёкмарева ОМ:'// Журнал техн.физ., т.60: — №10. — 1970.-С. 2025-2031.

16. Гринберг, Г. А. О решении обобщенной;задачи Стефана о промерзании жидкости, а также родственных, задач- теплопроводности, диффузии- и других. // Журн.техн.физики, т.37. №9.- 1967. - С. 1598-1606.

17. Данилюк, И. И; О задаче Стефана. // Успехи матем. наук, т.40. вып.5. — 1985.-С, 133-185.

18. Зайнуллин, Р. Г. Решение одной задачи переноса тепла методом ВГГП. / Р. Г. Зайнуллин, М. Н; Шафеев. Деп. в ВИНИТИ 23.11.88. - №7030-В89. — Уфа: УАИ; 1989. -11с.

19. Зайнуллин, Р. Г. Решение одной нестационарной задачи Стефана методом ВГТП для несимметричной области. / Р. Г. Зайнуллин, М. Н.

20. Шафеев. Деп. в ВИНИТИ 19.11.90. - №5821-В90. - Уфа: УАИ, 1990. -18 с.

21. Зайнуллин, Р. Г. Решение одной- задачи переноса; тепла при наличии; движущихся границ. / Р. Г. Зайнуллин, И. А. Акимов, М. Н. Шафеев -Деп; в ВИНИТИ 26.03.91. №1207-В91. - Уфа: УАИ, 1991.-7 с.

22. Зайнуллин, Р. Б. Решение; одной сопряженной? задачи: теплообмена методом интегральных преобразований:. / Р. Г. Зайнуллин, И. А. Акимов- Дет в ВИНИТИ 26.03.91 №1308-В91 - Уфа: УАИ, 1991. - 5 с.

23. Зайнуллищ Р. Г. Решение одной двухслойной задачи теплообмена со свободными границами; / Р. Р. Зайнуллин, И; А. Акимов* М. Н. Шафеев- Деп. в ВИНИТИ-26Ю3.91. №1309^391 - Уфа: УАИ; 199Г. - 6 с.

24. Зайнуллин, Р. Г. Об одном приложении метода ВГГП. / Р. Г. Зайнуллин, М. Н. Шафеев. Деп. в ВИНИТИ 22.07.94. - №1944-В94. - Уфа: УГАТУ, 1994.-6 с.

25. Зайнуллин, Р. Р. Решение одной плоской нестационарной задачи' Стефана при обобщенных условиях. / Зайнуллин, М< Н. Шафеев. Деп. в ВИНИТИ 27.03.97. - №980-В97. - Уфа: УГАТУ, 1997. - 9с.

26. Зайнуллин, Р. Г. Решение одной двумерной задачи; переноса тепла со свободными; границами. // Актуальные проблемы математики; Математические методы в естествознании : Межвуз. науч. сб; — Уфа: УГАТУ, 1999. С. 120-125.

27. Зайнуллин, Р. Г. Об одном приложении метода ВГГП. / Р. Г. Зайнуллищ М. Н. Шафеев, Р. Г. Самигуллина // Межвузовский научный сборник

28. Вопросы теории и расчета рабочих процессов тепловых двигателей. Уфа: УГАТУ , 20001- С. 120-125.

29. Зайнуллин, Р. Г. Об одном аналитическом подходе к решению одномерной1 задачи теплообмена с движущимися границами. // Трудым/н научной конференции. 24-28.06.08, г.Стерлитамак Уфа: Гилем, 2008, Т1.-242 с.

30. Зайнуллин, Р. Г. Об одном аналитическом подходе к^ решению одномерной задачи переноса-тепла со свободными границами. // Изв. Вузов. Математика, №2, 2008. - С. 24-31.

31. Зайнуллин, Р. Г. Математическое моделирование процесса теплообмена с фазовым переходом // Вестник УГАТУ. — Сер. Управление, информатика, вычислительная техника. Т. 13, №2(35), 2009 С. 265-279.

32. Ильин, В. А. Спектральная теория дифференциальных операторов. Самосопряженные дифференциальные операторы. М.: Наука, 1991.-366 с.

33. Калиев, И. А. Некоторые задачи линейной термоупругости в теории фазовых переходов Гинзбурга-Ландау. // Прикладная механика и техническая физика. 2003. -Т.44. -№6(262). С. 140-147.

34. Калиев, И: А. Задача Стефана как коэффициентная обратная задача. / И. А. Калиев, Э. В. Вагапова // Современные проблемы физики и математики: Труды Всероссийской научной конференции, Стерлитамак, 16-18сент., 2004, т. 1. Уфа, 2004. - С. 43-49.

35. Каменемостская, С. Л. О задаче Стефана. // Матем.сб. — т.53. №4. -1961.-С. 489-514.

36. Карташов, Э. М. Аналитические методы в теории теплопроводности твердых тел. М.: Высшая школа, 1985 - 480 с.

37. Карташов, Э. М. Аналитические методы решения краевых задач уравнения теплопроводности в области с движущимися границами. / Э. М. Карташов, Б. Я. Любов // Изв.АН СССР, серия Энергетика и транспорт. №6. - 1974. - С. 83-111.

38. Карташов, Э. М. Метод обобщенного интегрального преобразования при решении уравнения теплопроводности в области с движущимися границами. // Инж.-физ. журнал т.52. - №3. - 1987. - С. 495-505.

39. Карташов, Э. М. Метод функций Грина при решении краевых задач уравнения теплопроводности обобщенного типа. // Изв. АН СССР, «Энергетика и транспорт». №2. — 1979. - С. 108-116.

40. Карташов, Э. М. Термокинетика процессов хрупкого разрушения полимеров в механических, температурных и диффузионных полях. // Автор.дисс. на соиск.уч.степ.док.физ.-мат. наук.Л., ИБС АН СССР, 1982.-54 с.

41. Карташов, Э. М. Аналитические методы решения краевых задач нестационарной теплопроводности в областях с движущимися границами // Изв. РАН, Сер. Энергетика, 1999. №5. - С. 3-34.

42. Карташов, Э. М. Проблема Стефана. / Э. М. Карташов, А. Г. Рубин, JI. М. Ожерелкова // Математические модели физических процессов: Сборник научных трудов 10 Международной конференции, Таганрог, 29-30 июня, 2004. Таганрог, 2004. - С. 88-92.

43. Квальвассер, В. И. Метод нахождения функции Грина краевых задач уравнения теплопроводности для отрезка прямой с равномерно движущимися границами. / В. И. Квальвассер, Я. Ф. Рутнер // Докл.АН СССР, т. 156. №6. - 1964. - С. 1273-1276.

44. Коддингтон, Е. А. Теория обыкновенных дифференциальных уравнений. / Е. А. Коддингтон, Н. Левинсон М.: ИЛ, 1958. - 474 с.

45. Костюченко, А. Г. Распределение собственных значений. / А. Г. Костюченко, И. С. Саргсян М.: «Наука», 1979. - 400 с.

46. Кошляков, Э. М. Уравнения в частных производных математической физики. / Э. М. Кошляков и др. М.: Высшая школа, 1970. - 712 с.

47. Курант, Р. Методы математической физики. / Р. Курант, Д. Гильберт -М.: Гостехиздат, т.1, 1957. 476 с.с

48. Курант, Р. Методы математической физики. / Р. Курант, Д. Гильберт -М.: Гостехиздат, т.2, 1957. 620с.

49. Ладыженская, О. А. Линейные и квазилинейные уравнения параболического типа. / О. А. Ладыженская, В. А. Солонников, Н. Н. Уральцева. М.: Наука, 1967. - 736 с.

50. Левитан, Б. М. Разложения по собственным функциям дифференциальных уравнений второго порядка. — М.: Гостехиздат, 1950. -159 с.

51. Левитан, Б. М. Операторы Штурма-Лиувилля и Дирака. / Б. М. Левитан, И. С. Саргсян-М.: Наука, 1988.- 431 с.

52. Левитан, Б. М. Введение в спектральную теорию. / Б. М. Левитан, И. С. Саргсян -М.: Наука, 1970. 671 с.

53. Лыков, А. В. Тепломассообмен. М.: Энергия, 1978. - 479 с.

54. Люстерник, Л. А. Об автомодельных решениях некоторых уравнений с частными производными. // Вестник МГУ, мат-мех., №9, 1974. С. 4054.

55. Люстерник, Л. А. Краткий курс функционального анализа. / Л. А. Люстерник, В. И. Соболев. -М.: Высшая школа,1982. 271 с.

56. Марченко, В. А. Спектральная теория операторов Штурма-Лиувилля. -Киев: Наукова Думка, 1972. 219 с.

57. Мейрманов, А. М. Задача Стефана. Новосибирск: Наука, 1986. - 240 с.

58. Наймарк, М. А. Линейные дифференциальные операторы. — М.:. Наука, 1969.

59. Нгуен Дин Уи. Об одной задаче со свободной границей для параболического уравнения. // Вестник МГУ, мат.-мех., №2, 1966. -С.40-54.

60. Никифоров, А. Ф. Специальные функции математической физики / А. Ф. Никифоров, В. Б. Уваров М.: Наука, 1984. - 344 с.

61. Олвер, Ф. Введение в асимптотические методы и специальные функции. -М.: Наука, 1978.-375с.

62. Олейник, О. А. Об одном методе решения общей задачи Стефана. // Докл.АНСССР, 1960. -т.135. -№5, С. 1054-1057.

63. Плеснер, А. И. Спектральная теория линейных операторов. М.: Наука, 1966.-624 с.

64. Прудников, А. П. Интегралы и ряды. Дополнительные главы. / А. П. Прудников, Ю. А. Брычков, О. И. Маричев. М.: Наука, 1986. -800с.

65. Рубинштейн, Л. И. Проблема Стефана. Рига: Звайгзне, 1967. - 457 с.

66. Рубцов, Н. А. Численное моделирование однофазной задачи Стефана в слое с прозрачными и полупрозрачными границами. / Н. А. Рубцов, С. Д. Слепцов, Н. А. Саввинова // Прикладная механика и техническая физика, 2006. т. 47. - №3. - С. 84-91.

67. Сильченко, Ю. Т. Одна краевая задача для области с подвижной границей. // Известия вузов. Математика, 1998. №3. - С. 44-46.

68. Слейтер Люси Дж. Вырожденные гипергеометрические функции. М.: ВЦ АН СССР, 1968. - 178 с.

69. Титчмарш, Э. Ч. Разложения по собственным функциям, связаны с дифференциальными уравнениями второго порядка. -М.:ИЛ,т.Г, 1960.-279 с.

70. Титчмарш, Э. Ч. Разложения по собственным функциям, связаны с дифференциальными уравнениями второго порядка. -М.: ИЛ, т.2, 1961.-556 с.

71. Тихонов, А. Н. Уравнения математической физики. / А. Н. Тихонов, А. А. Самарский -М.: Наука, 1966. 724 с.

72. Уиттекер, Э. Т. Курс современного анализа. / Э. Т. Уиттекер, Д. Н. Ватсон-М.: Физматгиз, т.1, 1963. -335 с.

73. Уиттекер, Э. Т. Курс современного анализа. / Э. Т. Уиттекер, Д. Н. Ватсон- М.: Физматгиз, т.2, 1964. — 467 с.

74. Ушакова, В. И. Поведение решений однофазной задачи Стефана при больших значениях времени. / В. И. Ушакова, А. В. Клочков // Известия вузов. Математика, 2006. № 11. - С. 55-60.

75. Федорюк, М. В. Асимптотические методы для линейных обыкновенных дифференциальных уравнений. М.: Наука, 1983. -352 с.

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

77. Хакимов, X. Р. Замораживание грунтов в строительных целях. М: Госстройиздат, 1962. - 257 с.

78. Хуснутдинова, Н. В. О поведении решения задачи Стефана при неограниченном возрастании времени. Динамика сплошной среды. -Институт гидродинамики СО АН СССР. — Новосибирск: Наука, вып.2, 1969.-С. 168-1177.

79. Цыбин, А. М. К решению задачи Стефана. // Журнал техн. физики, т.64, 1974.-С. 2441-2444.

80. Чекмарева, О. М. О движении поверхности фазового перехода при больших временах в осесимметричной задаче Стефана. // Журнал техн. физики, т.65, №2, 1975. С. 209-213.

81. Чекмарева, О. М. Решение задачи Стефана, когда движение поверхности фазового перехода происходит по закону VZ // Журнал техн. физики, т.64, №10, 1974. С. 2043-2050.

82. Шафеев, М. Н. Об одном аналитическом подходе к решению задачи теплообмена с фазовым переходом. Деп. в ВИНИТИ 25.11.86. -per. №8457-В86. - Уфа: УАИ, 1986. - 14 с.

83. Шафеев, М. Н. Плоская нестационарная задача Стефана для полуполосы. Деп. в ВИНИТИ 14.05.87. - рег.№3350-В87. -Уфа: УАИ, 1987.- 15 с.

84. Шафеев, М. Н. Решение одной нелинейной задачи методом ВГГП. // Известия вузов. Математика, 1980. -№12(233). С.73-75.

85. Шафеев, М. Н. Решение одной плоской задачи Стефана методом ВГГП. // Инж.-физ. журнал, 1978. т.34. - №4. - С. 713-722.

86. Янке, Е. Специальные функции. / Е. Янке, Ф. Эмде, Ф. Леш. -М.: Наука, 1977. 342 с.

87. Baconneau Oliver. Smooth solutions to a class of free boundary parabolic problems. / Oliver Baconneau, Alessandra Lunardi. // Trans. Amer. Math. Soc, 2004. v. 356. - №3. - P. 987-1007.

88. Bertsch Michiel, Dal Paso Roberta, Franchi Bruno // Math.Ann. 1992. - 294. -№3. —P. 551-578

89. Borisovich Andrei. Symmetry-breaking bifurcations for free boundary problems. / Andrei Borisovich, Avner Friedman. // Indiana Univ. Math.J., 2005. -v.54. -№3. -P. 927-947.

90. Buchholr, H. Die Konbeente hypergeometrische Funktion. Springer-Verlag, Berlin, 1953.-187 p.

91. Busher, S. Developing the solution of Stefan's problem. / S. Busher, V. Georgiev // Докл. Бълг.А.Н. 1994. - 47. -№3. -P. 9-12.

92. Stefan J. Uber einige Probleme der Theorie der Warmelcitung. Sitzber. Wien. Akad. Mat. naturwiss. 1989. b.98. - 11 a. - P. 473.

93. Calor Gabriel. Asymptotic expansion of the solution of an interfase problem in a polygonal domain with thin layer. / Gabriel Calor, Martin Costabel, Monique Dauge, Gregory Vial. // Asymptotic Anal. 2006. 50. -№1-2. - P. 121-173.

94. Cheng K.C. Историческое развитие теории теплоты и термодинамики: обзор и некоторые наблюдения. Historical development of the theory of heat and thermodinamics: Review and some observations. // Heat transfer Eng., 1992.-13.-№3.-P. 19-37. Англ.

95. Dancer, E. N. A uniqueness theorem for a free boundary problem. / E. N. Dancer, Yihong Du. // Proc. Amer. Math. Soc, 2006. -v.134. -№11. -P. 3223-3230.

96. Friedman, A. Variational principles and free-boundary problems. New-York: A Wiley-Interscience publ., 1982.-709 p.

97. Frolova, E. V. One-phase Stefan problem with vanishing specific heat. // Международная конференция «Дифференциальные уравнения исмежные вопросы», посвященная памяти И.Г. Петровского. Сборник тезисов. -М.: МГУ, 2007, С, 93-94. Англ.

98. Kortea, М. К. On the geometry of the free boundary of the one-phase Stefan problem // Int.Conf. Differ. And Funct. Differ. Equat. Moscow, Aug.16-21, 1999: Abstr. Б.М., 1999 - С.И.

99. Moraczewski Krzystof. Решение задачи типа Даниловской для полупространства с подвижной границей. Rozwiazanie zagadnienia tupu Danilowskiej dla polprzestrzeni sprezystej zrachomym. Pr.IPPTPAN., 1988. -№47. P. 1 — 47. - Пол.

100. Morikawa Kichizo. Анализ. теплопроводности в телах с движущейся границей. Analysis of heat condution fields with moving boundary./ Kichizo Morikawa, Kensuke Kawashimo. // Bull JSME, 1984. 27. - №231. -P. 1938-1943 (англ.).

101. Mucha Piotr Bogustaw. Stefan problem in 2D case. Colloq. Math., 2006. -105.-№1.-P. 149-169.

102. Tokuda Naoyuki. Решение-задачи Стефана с использованием разложения Лагранжа-Бюрманна. Stefan problem by Lagrange-Burmann expansions. -Adv.Phase Change Heat Transfer: Proc.Int. Synp.Chongqing,May 20-23,1988.-Oxford, 1989.-P. 115-119.-Англ.

103. Wilson, D. G. Moving boundary problems. / D. G. Wilson, A. D. Solomon, P. T. Boggs Proc. Symp. and Workshop, Galtinburg, Term., Sept., 1977. -New York: Academic Press, 1978.