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

Тюлькина Елена Юрьевна

УСОВЕРШЕНСТВОВАННЫЙ МОМЕНТНЫЙ МЕТОД РЕШЕНИЯ

КИНЕТИЧЕСКОГО УРАВНЕНИЯ И ЕГО ПРИЛОЖЕНИЕ

К ЗАДАЧАМ ТЕПЛОПЕРЕНОСА В МОЛЕКУЛЯРНЫХ ГАЗАХ

Специальность 01.04.02 – теоретическая физика

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

Москва 2010 2

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

Научный руководитель: доктор физико-математических наук, доцент Савков Сергей Анатольевич

Официальные оппоненты: доктор физико-математических наук, профессор Кузнецов Михаил Михайлович кандидат физико-математических наук, доцент Алехин Евгений Иванович

Ведущая организация: Ярославский государственный университет им. П. Г. Демидова

Защита диссертации состоится " 14 " октября 2010 г. в 15 час. 00 мин.

на заседании диссертационного совета Д 212.155.07 в Московском государственном областном университете по адресу: 105005, Москва, ул. Радио, д. 10а.

С диссертацией можно ознакомиться в научной библиотеке Московского государственного областного университета.

Автореферат разослан " " сентября 2010 г.

Ученый секретарь диссертационного совета, кандидат физ.- мат. наук, доцент Барабанова Н.Н.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

Определяющую роль при теоретическом описании указанного явления играет число Кнудсена Kn = L, здесь – средняя длина свободного пробега молекул газа, L – характерный размер задачи. При L >> состояние газа определяется уравнениями динамики сплошной среды, для решения которых разработан широкий арсенал аналитических и численных методов. В случае, когда средняя длина свободного пробега молекул газа сравнима или больше характерной длины, фигурирующей в задаче, необходим учет дискретности строения газа, что требует рассмотрения кинетического уравнения.

Впервые математически корректный способ решения этого уравнения в приложении к задачам теплопереноса был предложен Лизом [1]. Основу данного метода составляет идея о сведении кинетического уравнения к системе уравнений переноса, для замыкания которой Лиз использовал двухстороннюю (четырехмоментную) функцию распределения. Такой алгоритм позволяет удовлетворить всем необходимым законам сохранения при использовании в функции распределения минимального числа моментов и делает возможным решение задачи в аналитической форме. Однако стандартный подход к реализации этого метода дает заниженное значение коэффициента скачка температуры. Другим принципиальным недостатком метода Лиза является произвол в выборе моментов функции распределения с одной стороны и составлении системы моментных уравнений – с другой. Можно показать, что использование различных наборов функций скорости может приводить к разным, а иногда и просто бессмысленным результатам. Более последовательным в данном отношении следует признать использование одного и того же набора разрывных функций скорости для составления функции распределения и моментных уравнений, как в методе полупространственных моментов [2].

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

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

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

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

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

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

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

Научная новизна работы.

1. Развит метод решения кинетического уравнения в приложении к задачам вычисления потока тепла в молекулярном газе для случая плоской, сферической и аксиальной геометрии.

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

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

4. Впервые из решения уравнения Ван Чанга – Уленбека с интегралом столкновений в форме Хансона – Морзе получены значения коэффициентов скачка температуры для конкретных молекулярных газов.

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

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

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

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

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

Положения, выносимые на защиту:

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

результаты решения кинетического уравнения в задачах теплопереноса в плоском, сферическом и цилиндрическом слоях молекулярного газа во всем диапазоне числа Кнудсена;

соотношения, определяющие зависимость потока тепла от характера аккомодации энергии;

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

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

Результаты диссертации представлялись на IX Международной научнопрактической молодежной конференции "Человек и космос. 50-летие космической эры". (Украина, Днепропетровск, 2007 г.); X, XI Международных молодежных научно-практических конференциях "Человек и космос" (Украина, Днепропетровск, 2008, 2009 гг.); 14 и 15 Всероссийских научных конференциях студентов-физиков и молодых ученых (ВНКСФ – 14, 15), (Уфа, 2008, 2009);

Всероссийской конференции по математике и механике (Томск, 2008); Международной конференции "Современные проблемы математики, механики и их приложений, посвященной 70-летию В.А. Садовничего" (Москва, 2009). А также докладывались и обсуждались на заседаниях и на научных семинарах кафедры “Теоретической физики и математического моделирования” Орловского государственного университета, кафедре “Высшей математики” Орловского государственного технического университета и кафедре “Теоретической физики” Московского государственного областного университета.

Публикации. По теме диссертации опубликовано 15 работ, включающих 9 статей в научных журналах и сборниках научных работ, 6 тезисов докладов на научных конференциях как всероссийского, так и международного уровня; в том числе работ, опубликованных в ведущих рецензируемых журналах, определенных ВАК РФ, – 3.

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

Структура и объем диссертации. Диссертация состоит из введения, трех глав, заключения, списка литературы из 127 наименований и приложений. Содержит 26 рисунков и 4 таблицы. Объем работы составляет 104 страницы.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

Первая глава диссертации посвящена рассмотрению процесса теплопереноса между параллельными пластинами.

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

Рассматривается слой газа толщиной d, заключенный между двумя неподвижными пластинами, на поверхности которых поддерживается постоянная температура Ts1 > Ts2. Полагается, что перепад T s = T s1 T s2 достаточно мал для того, чтобы ограничиться линейным приближением.

Состояние газа описывается уравнением Здесь V – вектор собственной (тепловой) скорости молекул газа, r – радиусвектор рассматриваемой точки пространства, J st – интегральный оператор столкновений, f – функция распределения.

В силу линейности поставленной задачи, решение уравнения (1) может быть представлено в виде где f 0 = n 2kT exp( C ) – равновесная функция распределения, – поправка, определяемая из решения соответствующего (1) линеаризованного уравнения:

Здесь C = V m 2kT0 – безразмерное значение скорости поступательного движения молекул газа, I st – линеаризованный оператор столкновений.

Конкретные расчеты проведены для БГК-модели:

За единицу длины выбрана величина и – коэффициенты температуро- и теплопроводности, T0 и n0 – некоторые, принятые за равновесные, значения температуры и концентрации молекул газа, k – постоянная Больцмана, m и s – масса и число степеней свободы молекул газа.

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

Значения nr и Trk определяются требованием отсутствия массового движения газа и характером аккомодации энергии Eik, Erk – обезразмеренные значения энергии, соответствующие приносимой падающими и уносимой отразившимися от поверхности k-той пластины молеk кулами; E s – энергия, которую уносили бы молекулы, если бы отражались с температурой Tsk.

Из (3) – (4) следует, что структура решения должна определяться соотношением = f1 ( z, C z ) + f 2 ( z, C z )C 2. Функции f1, f 2 могут быть разложены в ряд по любой полной системе ортогональных полиномов, при этом необходимо учитывать разрывный характер функции распределения на поверхности каждой из пластин. Сказанное позволяет представить искомую функцию в виде:

H ( x) = ( x + x ) 2 x – стандартная функция Хевисайда.

Коэффициенты ai±k ( z ) определяются из системы однородных дифференциальных уравнений, для составления которой кинетическое уравнение (3) следует последовательно умножить на входящие в (6) – (7) моменты и проинтегрировать по всему пространству скоростей.

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

На рисунках 1а) и 1б) представлены результаты расчетов в случае 1 = 2 = 1 и 1 = 2 = 0.826 соответственно. Там же приведены экспериментальные данные [4].

Рис. 1. Зависимость приведенного потока тепла от расстояния между пластинами. Линия 1 соответствует N = 1; 2 – N = 10 (при N = 2 – 9 кривые практически сливаются с линией 2); 3 – результаты стандартного метода Лиза; • – экспериментальные данные для аргона [4].

В режиме, близком к газодинамическому и условии полной аккомодации энергии, поток тепла можно представить в виде Здесь Ct – коэффициент скачка температуры, Kn = d, средняя длина свободного пробега определена соотношением:

Изложенный подход к решению кинетического уравнения дает Ct = 2.2193, 2.2076, 2.2057 и 2.2049 при N = 1, 2, 3 и 10 соответственно. Таким образом, при N = 1 отличие коэффициента скачка температуры от точного значения 2. [5] составляет порядка 0.7%, в случае N = 2 – 0.2%, N = 3 – менее 0.04%, а при N = 10 совпадает с результатами аналитического решения с точностью до приведенных в указанной работе числа цифр. Стандартный метод Лиза для данной модели кинетического уравнения дает заниженное (более чем на 14%) значение Ct = 1.875.

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

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

Конкретные расчеты проведены для модельного интеграла столкновений релаксационного типа [3]:

Рассмотрение этой модели представляет интерес в плане возможности сравнения значения коэффициента скачка температуры с известным аналитическим решением [3].

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

и nr – температура и концентрация молекул газа, отразившихся от поверхности k-той пластины. Значения nr, Trk, v и Trk, определяются требованием отk сутствия массового движения газа и характером аккомодации энергии Ev, i, E, i и Ev, r, E, r – обезразмеренные значения энергии поступательного и вращательного движения соответственно, приносимой падающими и уносимой отразившимися от поверхности k-той пластины молекулами; Ev, s, E, s – энергия, которую уносили бы молекулы, если бы отражались с температурой Tsk.

Поток тепла между пластинами задается выражением где В случае, когда обе пластины выполнены из одного материала ( v, 2 = v и, 2 = ), получено аналитическое выражение, задающее зависимость потока тепла от характера аккомодации энергии Здесь Qv и Q – значения потоков энергии поступательного и вращательного движения молекул газа у поверхности первой пластины, вычисленные при Соотношения (14) являются общими и не зависят от формы кинетического уравнения и способа его решения.

Для определения входящих в (14) параметров необходимо решать уравнение (3). Его решение также можно представить в виде (6), где Коэффициенты определяются из системы моментных уравнений, процесс составления которой описан выше.

В третьем параграфе рассматривается уравнение Ван Чанга – Уленбека [6] Здесь fl – функция распределения для молекул, находящихся в состоянии с энергией El.

С учетом линейности поставленной задачи, решение (15) может быть представлено в виде:

где – равновесная функция распределения, l = El kT0, = exp( l ).

Поправка l определяется из соответствующего (15) линеаризованного уравнения В качестве интеграла столкновений принята модель Хансона – Морзе [7]:

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

За единицу длины выбрана величина В качестве граничных условий принят закон чисто диффузного отражения молекул газа от поверхности каждой из пластин, что эквивалентно Значения nr, tr и int определяются требованием отсутствия массового движения газа и характером аккомодации энергии Искомый поток тепла может быть представлен в виде В случае, когда обе пластины выполнены из одного материала ( tr, 2 = tr и int2 = int ), получено аналитическое выражение, задающее зависимость потока тепла от характера аккомодации энергии где Под Qtr и Qint понимаются значения Q tr и Qint, вычисленные при tr = 1, Для определения входящих в (23) параметров необходимо решать уравнение (18). По аналогии с предыдущими параграфами решение представим в виде Коэффициенты ai±k определяются из системы дифференциальных уравнений, для составления которой кинетическое уравнение (18) необходимо последовательно умножить на все входящие в (24) моменты, просуммировать по всем l и проинтегрировать по пространству скоростей.

Четвертый параграф посвящен обсуждению результатов и сравнению с экспериментом.

Результаты расчетов в случае полной аккомодации энергии представлены на рисунках 2 и 3. На рисунке 3 также приведены результаты квадратурного метода, используемого в работе [9].

Рис. 2. Зависимость приведенного потока тепла от расстояния между пластинами. Линии 1 и 2 соответствуют результатам решения релаксационного кинетического уравнения усовершенствованным моментным методом при N = 1 и соответственно, 3 – стандартному методу Лиза. Для N = 2, 3 и 4 кривые почти полностью сливаются с линией 2;

Рис. 3. Значения приведенного потока тепла в слое азота при 1 = 2 = 1. Линия 1 соответствует решению кинетического уравнения в форме ХансонаМорзе моментным методом при N = 5; 2 – результаты квадратурного метода Как видно из рисунка 3, а также из приведенных в диссертации данных, в свободномолекулярном режиме результаты квадратурного и рассматриваемого методов совпадают. В случае d 1 отличие составляет порядка 1%. С ростом расстояния между пластинами погрешность квадратурного метода увеличивается, достигая 12% при d = 40. Очевидно, что в континуальном режиме можно говорить лишь о качественном характере результатов [9].

На рисунке 4 приведены результаты расчетов при различных коэффициентах аккомодации, соответствующим условиям эксперимента [4].

Рис. 4. Зависимость приведенного потока тепла от параметра A = Кривая 1 соответствует результатам решения уравнения в форме ХансонаМорзе в слое азота при tr = tr = int = int = 0.76, 2 – результатам решения релаксационного уравнения при тех же значениях коэффициентов аккомодации, 3 – результатам решения уравнения в форме Хансона-Морзе в слое воздуха при tr = tr = int = int = 0.95 ;, + – экспериментальные данные Теагана и Спрингера [4].

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

Таким образом, отличие коэффициента скачка температуры, полученного при решении кинетического уравнения релаксационного типа от точного значения 2.0579 [3] составляет порядка 0.6%, 0.2% и 0.009% при N = 1, 2 и 5 соответственно. Анализ представленных в диссертации результатов показывает, что как и в случае атомарного газа, указанная точность сохраняется во всем диапазоне значений числа Кнудсена.

Из таблицы и графиков следует, что использование интеграла столкновений релаксационного типа (10) и в форме Хансона – Морзе (18) для описания двухатомных газов дает фактически совпадающие результаты. В частности, в случае азота различие между значениями коэффициента скачка температуры, полученными при различных моделях, не превышает 0.4%, что оправдывает возможность использования аналога БГК-модели для описания данного класса явлений.

Вторая глава диссертации посвящена изучению процесса теплопереноса в сферическом слое газа.

В первом параграфе рассматривается двухатомный газ, заключенный между концентрическими сферами с радиусами R1 < R2, на поверхности которых поддерживается постоянная температура Ts1 > Ts2.

Введем сферическую систему координат, начало которой свяжем с центром сфер. Состояние газа определяется уравнением:

Cr проекция вектора C на направление r.

Искомый поток тепла может быть представлен в виде тепла с единицы поверхности внутренней сферы.

Учитывая разрывный характер функции распределения на поверхности каждой из сфер и тот факт, что с любой точкой в объеме газа связаны три инвариантных конуса в пространстве скоростей, границы которых молекулы пересекают только за счет столкновений между собой, решение уравнения (25) представим в виде Функции i,1 (первое приближение) содержат моменты, входящие в распределение Чепмена – Энскога, что обеспечивает переход решения в газодинамический режим. В i,2, i,3 и т.д. учтены моменты, которые получаются при подстановке предыдущего приближения в левую часть кинетического уравнения.

Во втором параграфе рассматривается уравнение Ван Чанга – Уленбека, которое в силу сферической симметрии задачи примет вид:

Решение (28) как и в случае двухатомного газа определяется выражением Искомый поток тепла может быть представлен в виде (26), где Для случая уединенной сферы, что эквивалентно R2, получено аналитическое выражение, задающее зависимость потока тепла от характера аккомодации энергии, которое совпадает по форме с (23), где Следует заметить, что при G = 1 соотношения (30) совпадают с приведенными в диссертации [8] выражениями для случая двухатомного газа.

Третий параграф посвящен обсуждению результатов.

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

На рисунке 6 приведены значения отношения Q Qgd, рассчитанные при R2 R1, N = 3 и условии полной аккомодации энергии для модели Хансона газодинамическом пределе. Там же приведены экспериментальные данные.

Рис. 5. Значения Q, вычисленные при условии полной аккомодации энергии для различного числа моментов, удерживаемых в функции распределения для кинетического уравнения релаксационного типа. Линия 1 соответствует N = 1;

2 – N = 2. В случае N > 2 графики сливаются с кривыми 2. Линия 3 – результаты метода Лиза.

Рис. 6. Зависимость отношения Q Qgd от параметра A = 6Kn 19 в слое воздуха; – экспериментальные данные [4].

Анализ рисунка свидетельствует о качественном согласии полученных результатов с данными экспериментов.

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

В первом параграфе рассматривается слой двухатомного газа, заключенный между коаксиальными цилиндрами с радиусами R1 < R2, на поверхности которых поддерживается постоянная температура Ts1 > Ts2.

Состояние газа определяется уравнением:

Искомый поток тепла задается выражением:

Решение уравнения (31) представим в виде (27), где вместо H1 и H При этом Следует заметить, что рассмотрение следующего приближения требует удержания большего, чем в задаче сферической геометрии, числа моментов, что приводит к существенному увеличению времени счета. Учитывая, что в случае сферического слоя отличие результатов, полученных при N = 2 и N = 3, составляет не более 0.1%, ограничимся учетом двух указанных приближений.

Во втором параграфе рассматривается уравнение Ван Чанга – Уленбека, которое в силу аксиальной симметрии задачи примет вид Искомый поток тепла определяется соотношением (32), где Решение (33) зададим выражением (29), где Третий параграф посвящен обсуждению результатов и сравнению с экспериментом.

На рисунке 7 показана зависимость приведенного потока тепла, рассчитанного для модельного кинетического уравнения релаксационного типа, от радиуса внутреннего цилиндра в случае полной аккомодации энергии.

На рисунке 8 представлены значения отношения Q Qgd, полученные при N = 2 для фиксированных R2 R1 и условии полной аккомодации энергии в слое азота и воздуха. Здесь Qgd = (3 + 2G )F1 R1 ln 2 – величина потоR1 T ка тепла, вычисленная в газодинамическом пределе. Там же приведены экспериментальные данные.

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

Рис. 7. Значения приведенного потока тепла при условии полной аккомодации энергии и фиксированных R1 R2, рассчитанные для релаксационного кинетического уравнения. Кривая 1 соответствует N = 1, 2 – N = 2; 3 – результаты стандартного метода Лиза.

Рис. 8. Значения отношения Q Q gd а) в слое азота; A = 5 ( 4 Kn ) ; б) в слое воздуха. Точками обозначены опытные данные [4, 10]. Линия 1 соответствует результатам решения уравнения Ван Чанга - Уленбека (33), 2 – кинетическому уравнению релаксационного типа (31).

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

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ

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

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

3. Впервые для модельного оператора Хансона – Морзе получены значения коэффициента скачка температуры для N 2, CO2, SO2 и воздуха.

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

5. Показано, что для азота в приложении к задачам теплопереноса в случае плоской, сферической и цилиндрической геометрии использование модели кинетического уравнения релаксационного типа [3] и модели ХансонаМорзе [7] дает практически совпадающие результаты.

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

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

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

ЦИТИРУЕМАЯ ЛИЕТРАТУРА

1. Lees L., Liu Chung-Yen. Kinetic theory description of conductive heat transfer from a fine wire // Phys. Fluids, 1962. – Vol. 5. – No 10. – P. 1137–1148.

2. Gross E.P., Ziering S. Heat flow between parallel plates // Phys. Fluids, 1959.

– Vol. 2, – No 6. – P. 701–712.

3. Латышев А.В., Юшканов А.А. Аналитическое решение задачи о скачке температуры с вращательными степенями свободы // Теор. и мат. физика, 1993. – Т. 95. – № 3. – С. 530–540.

4. Кошмаров Ю.А., Рыжов Ю.А. Прикладная динамика разреженных газов.

М., Машиностроение, 1977. – 182 с.

5. Латышев А.В. Применение метода Кейза к решению линеаризованного кинетического БГК уравнения в задаче о температурном скачке // Прикл.

матем. и механика. 1990. – Т. 54. – В. 4. – С. 581–586.

6. Wan Chang C.S., Uhlenbeck G.E., Boer J. The heat conductivity and viscosity of polyatomic gases. // Studies in Statistical Mechanics. – Amsterdam: North Holland Pablishing Company. 1964.

7. Hanson F.B., Morse T.F. Kinetic models for a gas with internal structure // Phys. Fluids, 1967. – Vol. 10. – No 2. – P. 345–353.

8. Савков С.А. Обобщенная теория теплопереноса в газовой среде при всех числах Кнудсена: дис. … доктора физ. мат. наук. – Москва, 2004. – 271 с.

9. Pazooki N., Loyalka S.K. Heat transfer in a rarefied polyatomic gas I. Plane parallel plates // J. Heat Mass Transfer, 1985. – Vol. 28. – No 11. – P. 2019– 10. Semyonov Yu.G., Borisov S. F., Suetin P.E. Investigation of heat transfer gases over a wide range of Knudsen numbers // J. Heat Mass transfer, 1984. – Vol. 27. – No. 10. – P. 1789–1799.

ОСНОВНОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ ОТРАЖЕНО

В СЛЕДУЮЩИХ РАБОТАХ

1. Савков С.А., Тюлькина Е.Ю. Об учете аккомодации энергии и вычислении потока тепла в плоском слое двухатомного газа. // ЖТФ, 2006. – Т.

76. – Вып. 2. – С. 25–29.

2. Савков С.А., Тюлькина Е.Ю. О решении кинетического уравнения при вычислении потока тепла в многоатомных газах. // ЖТФ, 2008. – Т. 78. – Вып. 7. – С. 16–20.

3. Савков С.А., Тюлькина Е.Ю. О решении кинетического уравнения в задаче вычисления потока тепла от сферической частицы в разреженном молекулярном газе // Письма в ЖТФ, 2009. – Т. 35. – Вып. 1. – С. 63–68.

4. Савков С.А., Тюлькина Е.Ю. О вычислении потока тепла между коаксиальными цилиндрами в разреженном молекулярном газе // Известия ОрелГТУ. Серия Фундаментальные и прикладные проблемы техники и технологии, 2009. – № 1 / 273 (559). – С. 35–40.

5. Савков С.А., Тюлькина Е.Ю. Об учете аккомодации энергии и вычислении потока тепла между параллельными пластинами // Вестник науки.

Сборник научных работ преподавателей, аспирантов и студентов физикоматематического факультета ГОУ ВПО «ОГУ». Выпуск 4. – Орел. 2005. – С. 151– 6. Савков С.А., Тюлькина Е.Ю. О решении кинетического уравнения в задаче вычисления потока тепла между параллельными пластинами // Вестник науки. Сборник научных работ преподавателей, аспирантов и студентов физико-математического факультета. Выпуск 6. – Орел: Издательство ОГУ, Полиграфическая фирма «Картуш», 2007. – С. 151–155.

7. Тюлькина Е.Ю. О вычислении потока тепла между концентрическими сферами в двухатомном газе // Вестник науки. Сборник научных работ преподавателей, аспирантов и студентов физико-математического факультета ГОУ ВПО «ОГУ». Выпуск 7. – Орел: Изд-во ОГУ, 2008. – С.

138–140.

8. Тюлькина Е.Ю. О вычислении потока тепла между концентрическими сферами в молекулярном газе // Альманах современной науки и образования. Математика, физика, строительство, архитектура, технические науки и методика их преподавания. – Тамбов: Грамота, 2008. – № 7(14). – С. 213–215.

9. Тюлькина Е.Ю. О вычислении потока тепла между коаксиальными цилиндрами в двухатомном газе // Альманах современной науки и образования. Математика, физика, строительство, архитектура, технические науки и методика их преподавания. – Тамбов: Грамота, 2008. – № 12 (19).

– С. 198–200.

10. Тюлькина Е.Ю., Савков С.А. К вычислению коэффициента скачка температуры в молекулярных газах. // Человек и космос. 50 – летие космической эры: сборник тезисов IX междунар. научн. – практич. конф. – Днепропетровск, 2007. – С. 31.

11. Тюлькина Е.Ю., Савков С.А. К вопросу о вычислении потока тепла между концентрическими сферами в двухатомном газе. // Человек и космос:

сборник тезисов X междунар. научн. – практич. конф. – Днепропетровск, 12. Тюлькина Е.Ю. К вопросу о вычислении потока тепла между концентрическими сферами в молекулярном газе. // Сборник тезисов, материалы Всероссийской научной конференции студентов – физиков и молодых ученых (ВНКСФ – 14): материалы конференции, тезисы докладов: В 1. – Т. 1 – Екатеринбург – Уфа: изд-во АСФ России, 2008. – С. 251–252.

13. Тюлькина Е.Ю. К вопросу о вычислении потока тепла от равномерно нагретой сферы в разреженном молекулярном газе // Всероссийская конференция по математике и механике, посвященная 130 – летию Томского государственного университета и 60-летию механико-математического факультета: (Томск, 22 – 25 сентября 2008 г.) – Томск: Томский государственный университет, 2008 г. – С. 174–175.

14. Тюлькина Е.Ю. К вопросу о теплопереносе между коаксиальными цилиндрами в разреженном молекулярном газе // Современные проблемы математики, механики и их приложений: материалы международной конференции, посвященной 70 – летию ректора МГУ академика В.А. Садовничего. – М.: Изд-во «Университетская книга», 2009. – С. 334–335.

15. Тюлькина Е.Ю. К вопросу о вычислении потока тепла через цилиндрический слой двухатомного газа при всех числах Кнудсена // Сборник тезисов, материалы 15 Всероссийской научной конференции студентовфизиков и молодых ученых (ВНКСФ – 15, Кемерово – Томск): материалы конференции, тезисы докладов: В 1. – Т. 1 – Екатеринбург – Кемерово:

изд-во АСФ России, 2009. – С. 263–264.