WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

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

Степаненко Евгений Николаевич

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

01.02.05 — «Механика жидкости, газа и плазмы»

АВТОРЕФЕРАТ

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

Челябинск – 2012

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

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

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

Ведущая организация: ОАО «Государственный ракетный центр имени академика В.П. Макеева», г. Миасс

Защита состоится «19» апреля 2012г. в 10:00 на заседании совета Д212.296.02 в Челябинском государственном университете по адресу:

454021, г. Челябинск, ул. Братьев Кашириных, 129, конференц-зал.

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

Автореферат разослан "___" марта 2012 года.

Ученый секретарь диссертационного совета кандидат физико-математических наук, доцент С.Ф.Долбеева

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

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

Как показал анализ соответствующей литературы, указанной проблеме – проблеме создания новых моделей гетерогенной среды - и в отечественной литературе [31-36] и за рубежом [1-30] уделяется серьезное внимание, особенно актуальна разработка гиперболических моделей для описания динамических многофазных сред с произвольным числом компонент смеси. Так, например, в работе [13] проведено исследование трех систем уравнений разных моделей на гиперболичность. Авторами показаны существующие области комплексных решений корней характеристических уравнений для двух систем уравнений моделей и отсутствие таковых для строго гиперболической модели. Похожий анализ встречается и в других работах, например в [27].

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

Указанные выше модели, как и большинство опубликованных в литературе, могут работать только с двухкомпонентной средой. Причем дифференциальные уравнения моделей записаны в виде, не позволяющем провести обобщение на случай среды состоящей из произвольного числа компонент. Отдельно отметим модели, работающие с n-компонентной средой - Lallemand M.-H. и др., [17], Куропатенко В.Ф. [33] и Сурова В.С.

[39-40]. Так, в работах Сурова В.С. предлагается «обобщенно-равновесная модель» (ОР модель), лежащая в основе настоящего исследования.

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

1. модификацию ОР модели с учетом теплопроводности среды;

2. модификацию ОР модели с учетом вязкости среды;

3. многоскоростную модель среды на базе ОР модели;

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

Для достижения поставленной цели решались следующие задачи:

Модификация n-компонентной обобщенно-равновесной модели с учетом теплопроводности среды;

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

Разработка модели многоскоростной многокомпонентной среды;

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

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

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

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

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

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

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

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

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

3. Модель многоскоростной многокомпонентной гетерогенной среды.

Численные методы расчета для данной модели.

4. Модель течения грунтовых вод с модифицированным законом Дарси.

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

6. Модификация метода Годунова С.К., предназначенная для интегрирования уравнений односкоростной многокомпонентной адиабатической смеси.

7. Обобщение соотношений Ренкина – Гюгонио на участках терпящих разрыв в случае многоскоростной многокомпонентной смеси.

Научная новизна и практическая ценность.

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

2. Представлена модификация теплопроводной модели односкоростной многокомпонентной среды, в которой исключена искусственно введенная несущая фракция, описанная в работе [38].

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

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

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

6. Разработан модифицированный метод Годунова С.К., предназначенный для интегрирования уравнений односкоростной многокомпонентной адиабатической смеси, теплопроводной и вязкой ОР моделей.

7. Обобщены соотношения Ренкина-Гюгонио в рамках ОР модели среды.

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

Апробация работы. Результаты исследований, вошедшие в диссертацию, докладывались:

- на X Международной конференции «Забабахинские научные чтения» 1519 марта 2010 года. – г. Снежинск;

- на VII Международной конференции «Лаврентьевские чтения по математике, механике и физике» посвященной 110-летию со дня рождения академика М.А.Лаврентьева. 23-27 августа, г. Новосибирск, 2010;

- на VII Всероссийской научной конференции «Фундаментальные и прикладные проблемы современной механики», г. Томск, 12-14 апреля 2011;

- на XVIII Международной конференции по вычислительной механике и современным прикладным и программным системам 25-31 мая 2011 года, – г. Алушта;

- на X Всероссийском съезде по фундаментальным проблемам теоретической и прикладной механики 23-30 августа 2011. г. Нижний Новгород.

На семинаре под руководством д.ф.-м.н. А.П. Яловца в ФБГОУ ВПО «ЮУрГУ» НОУ (сентябрь 2011 г.);

На семинаре под руководством д.ф.-м.н. О.Н. Дементьева в ФГБОУ ВПО "ЧелГУ" (октябрь 2011г., ноябрь 2011г. ).

На семинаре под руководством д.ф.-м.н. чл.-корр. РАН, профессор Э.Е.

Сон в ФГАОУ ВПО «МФТИ (ГУ)» (февраль 2012г.).

Опубликованность результатов.

По теме диссертации опубликованы 2 статьи в рецензируемых журналах, рекомендованных ВАК. Из 5 публикаций 4 статьи в сборниках трудов конференций, в том числе 3 - в международных.

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

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

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

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

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

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

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

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

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

Дифференциальные уравнения односкоростной модели гетерогенной среды с гиперболичным адиабатическим ядром представлены системой u div u t u u u 1grad p t (1) p u p pc2 div u t 0 i u i i Gi div u t i u i i (1 Gi ) div u t j u div u j j t 1 i p i где, i = 1,…, m – 1, j = m +1,…, n, n – число Gi c2 0 0 i i i p фракций в рассматриваемой модели. Скорость распространения возмущений в среде описывается формулой 1 m1 n p p i p i i 1 0 0 0 0 (2) i i i i i 2 i j m1 j j i c m i i i 0 0 p i1 p i i i i Также, впервые обобщены соотношения Ренкина-Гюгонио на участках терпящих разрыв в случае односкоростной многокомпонентной смеси в виде векторного соотношения F D U , (3) обобщающим соотношения Ренкина-Гюгонио на односкоростные многокомпонентные смеси и которые выражают законы сохранения массы, импульса и энергии через поверхность разрыва.

В выражении (3) приняты следующие обозначения:

T 0 0 U= ,u,e,11,,m0,11u,,m0 u,11e1,,m0 em,m+1,,n, m m m 0 F= u, p u2,( p e)u,11u, m0 u,1( p 1e1),,m( p 0 em), m m T 1( p 1e1)u,,m( p 0 em)u,m+1u,,nu.

m Здесь T – оператор транспонирования, символ обозначает скачок соответствующей величины при переходе через поверхность разрыва, u sus 0uнапример, .

i mi V i Vi V – истинная (физическая) плотность, i,,, u приведенная физическая плотность, объемная доля (пористость) i-го компонента смеси и вектор скорости.

– удельная внутренняя энергия смеси, определяемая из n ii выражения i В третьей главе представлены новые гиперболические модели:

односкоростная многокомпонентная модель с учетом теплопроводности среды; односкоростная многокомпонентная модель среды с учетом вязкости среды; модель многоскоростной многокомпонентной среды;

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

Также приведены некоторые модельные расчеты.

В первом параграфе третьей главы рассматривается модификация ОР модели с учетом теплопроводности [А1].

Дифференциальные уравнения модели представлены системой:

div(u) t u (u )u grad p F t 1 2 1 2 p t 2 u u W F u 2 u div W (u )W grad T W t (4) n i0 div ii0u ij Jik i t k n i i (u )i i p ij Jik i (u )i t i k 1 t n (k i ) n 1 Rik Qik ) u ij Jik i k 1 k n div u) ij J j j t j k 1 jk n n 0 где, T, = i i, 1, i i( p,i ) Ti( p,i ) i ii i1 i, T T (, p,1,1,,m1,0,m1,,n).

(, p,1,1,,m1,0,m1,,n) m1 mПо сравнению с работой [38] в данной модификации удалось избавиться от искусственно вводимой несущей фракции. Дополнительное слагаемое, входящее в четвертое уравнение системы (4), представляющее собой модифицированный закон теплопередачи Фурье и учитывающую релаксационные процессы в смеси, позволяет получить гиперболическую систему уравнений. В решении отсутствуют волны, распространяющиеся с бесконечно большими скоростями. Так модель показывает, что возмущения в смеси распространяются со скоростями:

m1 n 1 c2 1 c1 kpH c4 H 2c2 kpH 4 k i Gi +k i (1 Gi ) k kp 0 k 1 j i ij 2 i1 jm , m1 n 1 c2 1 c2 kpH c4 H 2c2 kpH 4 k i Gi +k i (1 Gi ) k kp 0 k 1 j i ij 2 i1 jm где с1 – скорость звуковых волн в среде, с2 – скорость тепловых волн.

Систему дифференциальных уравнений (4) можно записать в векторно-матричном виде:

U U, (5) S t x T T где, S 0,,0,W , U , u, p, 10, 1,,m1, m1,m1,,n,W u 0 0 0 0 0 0 0 0 u 1 0 0 0 0 0 0 0 c2 u 0 0 0 0 0 0 H 0 1G1 0 u 0 0 0 0 0 0 0 0 0 0 1 1 G1 0 0 u A 0 0 Gm1 0 0 0 u 0 0 0 0 m 0 u 0 0 0 m1 1 Gm1 0 0 0 0 m1 0 0 0 0 0 u 0 0 n 0 0 0 0 0 0 u k 0 kp k k k k k k u 1 m1 m1 n 1 m Данная система будет гиперболичной, так как корни характеристического уравнения (4) - действительные числа ( ), и собственные векторы матрицы А, u c1, u c2, u,, u, u c2, u cсоответствующие корням характеристического уравнения, линейно независимы, поскольку матрицу А можно представить в виде.

Модельный расчет методом КИР для задачи распада разрыва в одномерном случае дает следующие результаты (рис. 1).

Начальные условия для данной задачи сведены в таблицу 1. Все величины приведены в системе СИ.

Таблица p0 u0 0 T0 10 (1) 20 1 2 X 0 0.15*106 0 0.073 3130 1.67 1.4 1000 2.58*102 60.2*102 1X 0 0.1*106 0 0.1 2930 1.19 1.4 1000 2.58*102 60.2*102 1Рис.1. Зависимость параметров течения при распаде произвольного разрыва к моменту времени t=0,4 мс в водно-воздушной смеси (штриховые кривые); при отсутствии теплопроводности (сплошные); с удвоенным значением (штрихпунктирные).

Во втором параграфе третьей главы рассматривается модификация ОР модели с учетом релаксационных сил вязкости в смеси [A4]. Дифференциальные уравнения модели представлены системой.

u t x p u u 1 u t x x (6) u u p u 2 x 2 t u u t x x ii0 ii0u t x i i p ii0 ii ii0 u u t x i0 t x j jj t x Данную систему (6) также можно представить в векторноматричном виде (5), где , S 0,0,0, , 0,, , U , u, p, , 10, 1,,m1, m1,m1,,n u 0 0 0 0 0 0 0 0 0 u 1 1 0 0 0 0 0 0 (7) 0 c2 u 0 0 0 0 0 0 0 0 0 u 0 0 0 0 0 0 0 0 1G1 0 0 u 0 0 0 0 0 0 1 1 G1 0 0 0 u 0 0 0 0 A 0 0 Gm1 0 0 0 0 u 0 0 0 m 0 m1 1 Gm1 0 0 0 0 0 u 0 0 0 m1 0 0 0 0 0 0 u 0 0 n 0 0 0 0 0 0 0 0 u 1 i p i Здесь G .

c2 i 0 0 i i i p Для общего случая матрицу А тоже удается записать в виде. Собственные значения матрицы (7) - действительные числа:

, где C c2 - скорость звука в среде.

u C,u,...,u,u C В рамках данной модели приведены модельные расчеты двух задач.

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

Начальные данные приведены в системе СИ в таблице 2.

Таблица p0 u0 0 10 (1) 20 1 2 3 X 0 1*106 10 0.1 1000 1.4 1,19 100 1000 5000 0,X 0 1*106 10 0.1 1000 1.4 1,19 100 1000 5000 0,Рис 2. Зависимости к моменту t = 0.4 мс для режима течения с двумя волнами разрежения при µ= 100 кг/(м·с) (кривые 1); 1000 (2); 5000 (3).

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

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

Начальные данные задачи приведены в системе СИ в таблице 3.

Таблица p0 u0 0 10 (1) 20 X 0 0,15*106 0 0.09 1000 1.4 1,19 1000 0,X 0 0,1*106 0 0.09 1000 1.4 1,19 1000 0,Рис. 3. Зависимости искомых параметров к моменту t = 0.6 мс для режима течения с ударной волной и волной разрежения при µ = 1000 кг/(м·с) (кривые 1); для невязкой среды (2).

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

многоскоростная модель [A2,A5]. Дифференциальные уравнения модели представлены в виде системы дифференциальных уравнений.

u t x u u 1 p u t x x p p u c2 u t x t x ui ui i p ii0 ui + i(u ui ) t x x i0 ii0ui (8) i t x i pui u iui (u ui ) i i i i i i0 ui2 i0 ui t x x j juj t x uj uj p 0 j j uj + (u uj ) j j t x x n ii iгде i – модельный коэффициент, структура которого может быть расписана по аналогии [13]:

3CD i ii u ui,CD 0,44.

4di Систему уравнений (8) можно представить в векторно-матричном виде (5), где , U , u, p, 10, u1, 1,,m1, um1, m1, um1 m1,, unn 1(u u1) m1(u um1) m1(u um1) n(u un) S 0,0,0,0, 11,0,0,,0, m1m1, 0, m1m1,0,, nn,0 u 0 0 0 0 0 0 0 0 0 0 0 1 0 u 0 0 0 0 0 0 0 0 0 0 c2 u 0 0 0 0 0 0 0 0 0 0 H1 K1 u1 L1 0 0 0 0 0 0 0 0 0 0 u1 G1 0 0 0 0 0 0 0 0 1H1 1K1 0 с1 u1 0 0 0 0 0 0 0 0 1 1 G1 0 .

Hm1 Km1 0 0 0 um1 Lm1 0 0 0 0 0 0 1 0 0 0 0 um1 Gm1 0 0 0 0 m 0 m1Hm1 m1Km1 сm0 0 0 0 um1 0 0 0 0 0 0 Gmm1 m 0 0 0 0 0 0 0 0 0 um1 Gm1 0 0 m 0 0 0 0 0 0 0 0 0 m1 um1 0 0 0 0 0 0 0 0 0 0 0 un Gn n 0 0 0 0 0 0 0 0 0 0 0 n un Собственными значениями матрицы А, являются:

, где u c, u, u1 c1, u1,, um1 cm1, um1,, um1 cm1,, un cn 1 m1 n p p i i i 1 p(i )2 0 0 0 , i i i i i jm1 j j i c m i i i 0 0 p i1 p i i i i pFi для i = 1, …, m – 1, c , i iGi iFi для j = m + 1, …, n, c p.

j j Матрицу А можно представить в виде A=-1 – следовательно, наша модель строго гиперболична.

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

Таблица p u u1 u2 g 1 2 g 1 X 0 0.12*105 0 0 0 1.19 1000 500 0.8 0.1 0.X 0 0.1 105 0 0 0 1.19 1000 500 0.8 0.1 0.p/pR u 1 1, 1,0,011,1,0,001,-40 -20 0 20 40 -40 -20 0 20 40 x x -40 -20 0 20 x uu1 2 0,01 1,2 1 0,6 0,010,0,0,000,0,0,000,0,x -40 -20 0 20 40 -40 -20 0 20 x -40 -20 0 20 x Рис. 4. Распад разрыва в трехкомпонентной смеси при j=0 (1), 3кг/(м3·с)(2) В четвертом параграфе третьей главы рассматривается ударная адиабата смеси без привлечения гипотезы об аддитивном ударном сжатии компонентов смеси. Получены уравнения, представляющие собой общую форму соотношений, связывающие параметры по обе стороны от поверхности разрыва. Эти соотношения обобщают известные соотношения Ренкина – Гюгонио на многоскоростные многокомпонентные смеси и выражают законы сохранения потоков массы, импульса и энергии через поверхность разрыва:

,, s(us D) 0(u0 D) uss(us D) ps u00(u0 D) p;

s(us D) s us2 2 psus 0(u0 D) 0 u02 2 p0u , i = 1, …, m – 1, is(uis D) i0(ui0 D), (9) isuis(uis D) is ps i0ui0(ui0 D) i0 p;

is(uis D) is uis2 2 isuis ps i0(ui0 D) i0 ui02 2 i0ui0 p , j = m + 1, …, n, (ujs D) (uj0 D) js j.

j ujs(ujs D) ps i0j uj0(uj0 D) pjs js jРассчитана задача о движении ударной волны в пузырьковой жидкости (рис 5.). Ударная волна перемещается по неподвижной (u0 = 0, p0 = 105 Па, OL=0.98) водно-воздушной смеси. Коэффициенты уравнения состояния (3.4.1) для воды равны: L=5.59, *L=1000 кг/м3, c*L =1515 м/с.

Для газовой составляющей – G=1.4, *G=1.24 кг/м3, c*G =0.

а) s, 1s б) в) us, D 1s 90,060,0640,0320,00 ps / p20 40 60 80 1ps / p0 ps / p0 20 40 60 80 100 0 20 40 60 80 1Рис. 5. Зависимости параметров течения,,,, за ударным фронтом от D s 1s us 1s интенсивности ударной волны ps pВ пятом параграфе третьей главы представлена модель течения грунтовых вод в пористой среде [A2], описываемая следующей системой уравнений:

h u h h u t x x u u p u u (10) t x x p gh const В диссертации показано, что система уравнений (10) строго гиперболична. Если же использовать оригинальный закон Дарси с =0, то скорости перемещения и возмущения становятся бесконечно большими.

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

В четвертой главе рассмотрены численные методы расчета течений для обобщенно-равновесной модели среды [A3].

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

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

На примере многоскоростной модели проведены расчеты и сравнение с экспериментальными данными [35]. Так было рассмотрено обтекание клина с различными углами (1-40, 2-100, 3-140) пузырьковой жидкостью со скоростью, превышающей скорость звука. Результаты представлены на рис.6.

D, м/с 3 1 I - - II 2III 10,0 0,3 0,6 0,p, МПа 0 5 10 15 g 0 gl 0 l MI 0.01 0.85 0.II 0.02 0.85 0.III 0.05 0.5 0.Рис. 6. Сравнение расчетных данных и Рис. 7. Зависимость скорости ударной волны экспериментальных при обтекании от давления за ее фронтом.

пузырьковой жидкости около клина. 1- ОР модель в адиабатическом приближении.

2- многоскоростная модель.

3- экспериментальные значения.

Применяя соотношения (9) полученные для многоскоростной среды, рассчитана скорость прохождения УВ по покоящейся однородной газо-жидкостной смеси. На рис. 7 приведены расчетные и экспериментальные [8] зависимости скорости распространения УВ в водно-глицериновом растворе с пузырьками воздуха от давления за ее фронтом при различной концентрации газа в смеси (ро=105 Па).

ЗАКЛЮЧЕНИЕ Представлена новая гиперболичная модель многоскоростной многокомпонентной среды, основанная на законах сохранения. Показано, что уравнения модели среды относятся к гиперболическому типу. С использованием численного метода Куранта – Изаксона – Риса рассчитан ряд задач о распаде произвольного разрыва в запыленном газе;

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

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

На базе ОР модели разработаны и уточнены соотношения Ренкина-Гюгонио для многоскоростной, многокомпонентной среды;

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

Необходимо отметить, что в рамках предложенных моделей хорошо зарекомендовал себя и метод Куранта – Изаксона – Риса;

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

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

СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ:

Список публикаций из списка ВАК:

А1. Степаненко Е.Н. Сеточный метод характеристик для расчета течений односкоростной многокомпонентной теплопроводной среды / Е.Н.

Степаненко, В.С. Суров // Вестник Челябинского государственного университета. 2010. №24 (205) Физика. Вып. 8. С. 15-22.

А2. Степаненко Е.Н. Новые гиперболические модели в механике многофазных сред / Е.Н. Степаненко, В.С. Суров // Вестник Нижегородского университета им.Н.И.Лобачевского. - 2011. - Вып.4(5).

- С.2518-2519.

Статьи в иных журналах и сборниках трудов конференций:

А3. Степаненко Е.Н. Модификация метода Годунова для расчета течений односкоростных многокомпонентных адиабатических смесей / Е.Н.

Степаненко, В.С. Суров // Тезисы докладов. Международная конференция. Лаврентьевские чтения по математике, механике и физике. Новосибирск. 2010. С.231-232.

А4. Степаненко Е.Н. Новые гиперболические модели в механике сплошных сред. Разработка и численная реализация / Е.Н. Степаненко, В.С. Суров // Забабахинские научные чтения: Сборник тезисов докладов X Международной конференции, 15-19 марта 2010 года. – РФЯЦ – ВНИИТФ, г. Снежинск. С. 259.

А5. Степаненко Е.Н. Динамика многофазных сред. Новый подход / Е.Н.Степаненко, В.С. Суров // Материалы XVIII Международной конференции по вычислительной механике и современным прикладным и программным системам, 25-31 мая 2011 года. – г.

Алушта. С.618-620.

СПИСОК ЦИТИРУЕМОЙ ЛИТЕРАТУРЫ:

1 Abgrall R., Karni S. Computations of compressible multifluids / J. Comp.

Phys. 2001. V.169. P. 594–623.

2 Allaire G., Clerc S., Kokh S. A five-equation model for the simulation of interfaces between compressible fluids / J. Comp. Phys. 2002. V.181. P.

577–616.

3 Andrianov N., Saurel R., Warnecke G. A simple method for compressible multiphase mixtures and interfaces, Int. J. Numer. Meth. Fluids 2003;

41:109–131.

Baer M., Nunziato J. Two-phase mixture theory for the deflagration to detonation transition (DDT) in reactive granular materials. International journal Multiphase Flows.12:861-889,1986.

5 Chung M.S., Lee S.J. A modified semi-implicit method for a hyperbolic twofluid model, Applied Numerical Mathematics 59 (2009) 2475–246 Cook T.L., Harlow F.H. Vortices in bubbly two-phase flow, Int. J.

Multiphase Flow Vol.12, № 1,pp.35-61,1986.

7 Drew D. A., Mathematical modeling of two-phase flow, Ann. Rev. Fluid Mech. 1983. 15:261-91.

8 Eddington R. B. Investigation of supersonic phenomena in a two-phase (liquid – gas) tunnel / AIAA J. 1970. V. 8, N. 1. P. 65–74.

9 Embid P., Baer M. Mathematical analysis of two-phase continuum mixture theory, Continuum Mech. Thermodyn. 4 (1992) 279-312.

10 Fedkiw R.P., Marquina A., Merriman B. An Isobaric Fix for the Overheating Problem in Multimaterial Compressible Flows, Journal of Computational Physics 148, 545–578 (1999).

11 Gavrilyuk S.L. Acoustic properties of a two-fluid compressible mixture with micro-inertia, European Journal of Mechanics B/Fluids 24 (2005) 397–412 Herard J.-M. A three-phase flow model//Mathematical and computer modeling. 2007. V. 45. P. 732 – 755.

13 Hudson J., Harris D. A high resolution scheme for Eulerian gas–solid twophase isentropic flow, Journal of Computational Physics 216 (2006) 494– 525.

14 Ivings M.J., Causon D.M., Toro E.F. On Riemann solvers for compressible liquids, Int. J. Numer. Meth. Fluids 28: 395–418 (1998).

15 Jayanti S., Valette M., Prediction of dryout and post-dryout heat transferat high pressures using a one-dimensional three-fluid model, International Journal of Heat and Mass Transfer 47 (2004) 4895–4910.

16 Kapila A. K., Son S.F., Bdzil J.B., Menikoff R., Stewart D. S., Two-phase modeling of DDT: Structure of the velocity-relaxation zone. Phys. Fluids (12), December (1997), 3885-3897.

17 Lallemand M.-H., Chinnayya A., Metayer O. Le, Pressure relaxation procedures for multiphase compressible flows, Int. J. Numer. Meth. Fluids 2005; 49:1–56.

18 Leea S.J., Chang K.S., Kim K. Pressure wave speeds from the characteristics of two Fuids, two-phase hyperbolic equation system, International Journal of Multiphase Flow 24 (1998) 855-866.

19 Lowe C.A., Two-phase shock-tube problems and numerical methods of solution, Journal of Computational Physics 204 (2005) 598–632.

20 Murrone A., Guillard H. A five-equation reduced model for compressible two phase flow problems / J. Comp. Phys. 2005. V.202. P. 664–698.

21 Ndjinga M., Kumbaro A., Vuyst F., Laurent-Gengoux P., Numerical simulation of hyperbolic two-phase flow models using a Roe-type solver, Nuclear Engineering and Design 238 (2008) 2075–2083.

22 Pokharna H.,Mori M., Ransom V.H., Regularization of Two-Phase Flow Models: A Comparison of Numerical and Differential Approaches.Journal of computational physics,№ 134,282-295 (1997).

23 Ransom V. H.,Hyperbolic Two-Pressure Models for Two-Phase Flow, Journal of computational physics,№ 53,124-151 (1984).

24 Romenski E., Toro E.F., Hyperbolicity and one-dimensional waves in compressible two-phase flow models, Shock Waves (2004) 13: 473–487.

25 Romenski E., Drikakis D., Toro E. Conservative Models and Numerical Methods for Compressible Two-Phase Flow, J Sci Comput (2010) 42: 68–26 Saurel R., Mtayer O., Massoni J., Gavrilyuk S. Shock jump relations for multiphase mixtures with stiff mechanical relaxation, Shock Waves (2007) 16:209–227 Thanh M.D. On a two-fluid model of two-phase compressible flows and its numerical approximation, Commun Nonlinear Sci Numer Simulat 17 (2012) 195–228 Trapp J. A., Welch S. W. J. An interface constitutive modeling problem for compressible flows with mass transfer//Int. J. Eng. Sci. 2010. V. 48. P. 18– 1906.

29 Wackers J., Koren B., A fully conservative model for compressible two-fluid flow, Int. J. Numer. Meth. Fluids 2005; 47:1337–1330 Zein A., Hantke M., Warnecke G. Modeling phase transition for compressible two-phase flows applied to metastable liquids//J. Comput.

Phys. 2010. V. 229. P. 2964 – 2998.

31 Волосевич П. П., Леванов Е. И., Северина Е. В. Решения типа бегущих волн с учетом гиперболического переноса / ИФЖ. 2008. Т. 81. № 2. С.

290 – 302.

32 Годунов С.К., Роменский Е.И., Элементы механики сплошных сред и законы сохранения. - Новосибирск: Научная книга,1998. - 280 с., ил.

(Университетская серия. Т. 4) 33 Куропатенко В.Ф., Новые модели механики сплошных сред. // Инженерно-физический журнал. 2011. Т. 84, №1. С. 74-92.

34 Ляхов Г. М. Ударные волны в многокомпонентных средах / Изв. АН СССР. Механика и машиностроение. 1959. №1. С. 46–49.

35 Наноряков В.Е., Покусаев Б.Г., Шрейбер И.Р. И др., Экспериментальное исследование ударных волн в жидкости с пузырьками газа // Волновые процессы в двухфазных системах.

Новосибирск. 1975. С. 54-97.

36 Нигматулин Р. И. Динамика многофазных сред. М.: Наука. 1987. Ч.1.

464 с.

37 Рахматулин Х. А. О распространении волн в многокомпонентных средах / ПММ. 1969. Т.33. Вып.4. С. 598–601.

39 Суров В. С. О локализации контактных поверхностей в многожидкостной гидродинамике. // ИФЖ. 2010. Т.83. №3. С. 518–527.

46 Суров В.С. Гиперболическая модель односкоростной многокомпонентной теплопроводной среды. // ТВТ. 2009. Т.47. №6. C.

905-913.

47 Суров В.С. Задача Римана для односкоростной модели многокомпонентной смеси // ТВТ. 2009. Т.47. №2. C. 283-291.

48 Суров В.С. Односкоростная модель гетерогенной среды с гиперболичным адиабатическим ядром // Журнал вычислительной математики и математической физики, 2008, том 48, №6, 1111-1125.

38 Суров В.С. Односкоростная модель многокомпонентной теплопроводной среды // ИФЖ. 2010. Т.83. №1. C. 132-141.

40 Суров В.С. Ударная адиабата односкоростной гетерогенной среды. // ИФЖ. 2006. Т. 79. №5. С. 46–52.

Подписано к печати 6.03.2012 г.

Формат 6084 1/16 Объем 1,0 уч.изд.л.

Тираж 100 экз. Бумага офсетная.

Отпечатано на ризографе в типографии МБОУ лицей №31 города Челябинска.

454080, г. Челябинск, ул. Володарского, 18.




© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.