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

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

УДК 532.517.532.546 Логвинов Олег Анатольевич Особенности неустойчивого вытеснения вязкой жидкости из ячейки Хеле-Шоу при больших числах Пекле 01.02.05 – механика жидкости, газа и плазмы

АВТОРЕФЕРАТ

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

2 Москва – 2011

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

М. В. Ломоносова Научные руководители: доктор физико-математических наук, профессор Н. Н. Смирнов доктор физико-математических наук, О. Е. Ивашнев

Официальные оппоненты: доктор физико-математических наук, профессор А. Г. Петров кандидат физико-математических наук, А. А. Осипцов

Ведущая организация: РГУ нефти и газа им. И. М. Губкина

Защита состоится 21 октября 2011 года в 15 часов 00 мин. на заседании Диссертационного совета Д 501.001.89 при Московском государственном университете имени М. В. Ломоносова по адресу: 119992 ГСП-2, Москва, Ленинские горы, Главное здание МГУ, механико-математический факультет, аудитория 16-

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

Автореферат разослан « » сентября 2011 года.

Ученый секретарь диссертационного совета доктор физико-математических наук А. Н. Осипцов

Общая характеристика работы

Актуальность темы. Ячейка Хеле-Шоу представляет собой две параллельные пластины, разделенные малым зазором, существенно меньшим, чем ширина пластин. Подобные ячейки используются для аналогового моделирования фильтрационных течений в пористых средах. Предметом настоящего исследования является вытеснение вязкой жидкости, находящейся в ячейке, другой – менее вязкой, подающейся в ячейку. Вытесняющая жидкость стремится прорваться через слой вытесняемой, образуя в ней каналы, называемые вязкими пальцами (viscous fingers, Saffman P. G., Taylor G., 1958). В результате развивается так называемая неустойчивость Саффмана – Тейлора.

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

Экспериментально установлено, что критерием подобия, определяющим характерный размер вязких пальцев в случае, когда вытесняющая и вытесняемая жидкости не смешиваются, является модифицированное капиллярное число (Park C. W., Homsy G. M., 1985). В случае смешивающихся жидкостей таким критерием выступает число Пекле (Tan C. T., Homsy G. M., 1986). Однако существуют эксперименты по скоростному (числа Пекле больше 1000) вытеснению водоглицериновой смеси водой (смешивающиеся жидкости), в которых этот критерий не срабатывает (Smirnov N.N., Nikitin V.F., Ivashnyov O.E., etc., 2002–2004). Образующиеся пальцы имеют конечную ширину, не зависящую от скорости вытеснения.

Целями работы являются: установление критерия подобия и выделение стабилизирующего фактора вытеснения вязкой жидкости из ячейки Хеле-Шоу менее вязкой жидкостью, смешивающейся с вытесняемой, при больших числах Пекле, когда эффекты диффузии малы по сравнению с проникновением одной жидкости в другую за счет развивающейся неустойчивости;

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

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

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

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

постановке (неустойчивость типа Саффмана – Тейлора).

Об устойчивости границы раздела при сдвиговом течении в линейной 2.

постановке (неустойчивость типа Кельвина–Гельмгольца).

Численное моделирование вытеснения вязкой жидкости из ячейки Хеле 3.

Шоу.

Научная новизна заключается в том, что автором впервые:

Выявлен стабилизирующий фактор, определяющий структуру фронта 1.

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

Показано, что за счет инерции вытесняющей жидкости на боковой 2.

поверхности пальцев развивается неустойчивость. По схожести механизмов обнаруженная неустойчивость аналогична неустойчивости Кельвина – Гельмгольца (Ламб Г., 1947). Силы вязкости в плоскости ячейки приводят к стабилизации боковой поверхности.

Достоверность результатов подтверждается качественным соответствием с экспериментальными данными и численным моделированием.

Научная и практическая ценность заключается в том, что:

Осреднением уравнений Навье – Стокса по зазору ячейки Хеле-Шоу 1.

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

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

Хеле-Шоу при больших числах Пекле.

На защиту выносятся:

Осредненные уравнения движения вязкой жидкости в ячейке Хеле-Шоу, 1.

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

Критерий подобия вытеснения вязкой жидкости из ячейки Хеле-Шоу при 2.

больших числах Пекле.

Объяснение механизма развития неустойчивости на боковой поверхности 3.

пальцев.

Метод моделирования вытеснения вязкой жидкости из ячейки Хеле-Шоу 4.

при больших числах Пекле.

Публикации и апробация диссертации. Основные результаты достаточно полно отражены в 10 публикациях, 3 из которых представлены в журналах перечня ВАК. Все положения, вынесенные на защиту, получены автором самостоятельно.

Результаты работы докладывались на следующих научно исследовательских семинарах и конференциях: семинаре кафедры газовой и волновой динамики под руководством академика Р. И. Нигматулина (МГУ им.

М. В. Ломоносова, 2010 г.), семинаре по многофазным средам под руководством профессора Н. Н. Смирнова (МГУ им. М. В. Ломоносова, 2008– 2010 гг.), семинаре кафедры аэромеханики и газовой динамики под руководством академика Г. Г. Черного (МГУ им. М. В. Ломоносова, 2009 г.), семинаре по механике сплошной среды под руководством академика РАН А. Г.

Куликовского, профессора В. П. Карликова и член-корр. РАН О. Э. Мельника (НИИ механики МГУ, 2011 г.);

конференции «Ломоносовские чтения» (МГУ им. М. В. Ломоносова, 2007–2009 гг.), конференции–конкурсе молодых ученых НИИ механики МГУ (2008 г.), Всероссийской конференции «Современные проблемы механики сплошной среды», посвященной столетию со дня рождения академика Л. И. Седова (Москва, 2007 г.), международной конференции «Современные проблемы газовой и волновой динамики», посвященной столетию со дня рождения академика Х. А. Рахматулина (Москва, 2009 г.), Одиннадцатой международной школе-семинаре «Модели и методы аэродинамики» (Евпатория, 2011 г.) Структура и объем. Диссертация состоит из введения, четырех глав, выводов и списка литературы (92 наименования);

изложена на 117 страницах с 27-ю рисунками и 2-мя таблицами.

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

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

В рамках классической модели (Saffman P. G., Taylor G., 1958) анализируется устойчивость первоначально плоского фронта вытеснения, движущегося по своей нормали. Модель включает уравнение неразрывности и аналог закона Дарси для ячейки Хеле-Шоу. Решается задача на малые возмущения в бесконечной области: слева менее вязкая жидкость («1»), справа более вязкая («2», рис. 1). На границе раздела жидкостей ставятся условия равенства нормальных к фронту скоростей и равенства давлений. На бесконечности – условие постоянства скорости.

12 µ 12 µ p p u v = 2 u, = 2 v, + = 0;

x y x y d x = ( y, t ) : p1 = p2, u1 = u2 = ;

dt x + : u2 = U 0, v2 = 0, x : u1 = U 0, v1 = 0.

u, v, p – осредненные по зазору ячейки компоненты скорости и давление, µ – вязкость жидкости.

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

Линейный анализ приводит к дисперсионному соотношению между длиной волны возмущения и интенсивностью его роста h µ2 h M 1, = h= M=, h=,. (1) M +1 µ1 U h = h ( ) для M = 100 показана на рис. 2 кривой 1.

Зависимость Видно, что при M 1 процесс вытеснения более вязкой жидкости менее вязкой абсолютно неустойчив, декремент затухания строго положителен для всех длин волн. Причем чем меньше длина волны возмущения, тем больше интенсивность его роста, так что возмущения бесконечно малой длины нарастают с бесконечной скоростью. Таким образом, на начальной стадии формирования вязких пальцев мгновенно должно образоваться огромное число тонких «каналов». Это представляется логичным, так как классическая модель учитывает только дестабилизирующий фактор – силу трения о пластины ячейки. При вытеснении более вязкой жидкости менее вязкой она приводит к развитию неустойчивости Саффмана – Тейлора.

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

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

2 3 / 2 p1 = p2 1 + +.

y y – коэффициент поверхностного натяжения между жидкостями.

С учетом поверхностного натяжения получается дисперсионное соотношение U 0 µ M 1 1 M h= Ca =,. (2) M + 1 12 3 ( M + 1)Ca Зависимость h = h ( ) для Ca = 0.1 и M = 100 показана на рис. кривой 2. Учет поверхностного натяжения в плоскости ячейки привел к стабилизации коротковолновых возмущений. Существует критическая длина волны *, интенсивность роста которой наибольшая. Считается, что вязкие пальцы имеют ширину наиболее быстро растущего возмущения. Для нее получена аналитическая оценка.

Рис. 2. Зависимость темпа роста возмущений от длины волны: 1 – по решению (1) при отсутствии поверхностного натяжения ( Ca = );

2 – по решению (2) с учетом поверхностного натяжения ( Ca = 0.1 ). В обоих расчетах M = 100.

В предположении, что течения в двух ячейках подобны, если в них образуется одинаковое число пальцев, получается критерий подобия Uµ l l Ca = 0 * = = const,.

, ( µ 2 µ1 )U 0 */2 l – ширина ячейки.

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

Если жидкости смешиваются, структура фронта вытеснения определяется молекулярной диффузией. Коротковолновые возмущения стабилизируются, так как прорвавшиеся тонкие пальцы диффундируют в вытесняемую жидкость, U 0l растворяясь в ней. Критерием подобия выступает число Пекле Pe =, D– D коэффициент диффузии.

В работе анализируются эксперименты с иным механизмом стабилизации.

Водоглицериновая смесь вытесняется водой (смешивающиеся жидкости) при числах Пекле больших тысячи (Smirnov N.N., Nikitin V.F., Ivashnyov O.E., etc., 2002–2004). Пальцы имеют конечную ширину, не зависящую от скорости вытеснения, поэтому число Пекле как критерий подобия не срабатывает.

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

В главе II уравнения Навье – Стокса несжимаемой вязкой жидкости осредняются по зазору ячейки. Интегрирование по вертикальной координате, перпендикулярной плоскости ячейки, приводит к уравнениям «типа Рейнольдса» (Шлихтинг Г., 1974) с осредненными параметрами течения и отклонениями от средних значений. Они вычисляются в предположении, что профиль скорости между пластинами параболический. Получаются осредненные уравнения, которые в отличие от классической модели содержат члены, характеризующие инерционные свойства жидкости и силы вязкости в плоскости ячейки.

2u 2u 12 µ u u p u t + u x + v y + x = 2 u + µ x 2 + y 2, 2v 2v 12 µ v v p v = 2 v + µ 2 + 2, + u + +v t y y x y x u v + = 0, x y /2 /2 / 1 1 / 2 / vdz, / 2pdz.

u= udz, v = p= 2 (черточки, означающие осреднение, в дальнейшем опускаются).

Показано, что учет инерции жидкости не приводит к стабилизации фронта вытеснения. Соответствующая зависимость h = h ( ) имеет вид U M = 84, D= = 1, Re = 2 0 = 1.2 ).

гиперболы (рис. 3 кривая 1, 1 µ Выбранные значения безразмерных параметров соответствуют экспериментам Smirnov N.N., Nikitin V.F., Ivashnyov O.E., etc., 2002–2004.

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

Поэтому в нее добавляются силы вязкости в плоскости ячейки (приходим к уравнениям «типа Бринкмана», Brinkman H. C., 1947). На границе раздела ставятся условия равенства векторов скоростей и векторов напряжений.

2u 2u 2v 2v 12 µ 12µ p p = 2 u + µ 2 + 2, = 2 v + µ 2 + 2, x y x y x y u v + = 0;

x y x = ( y, t ) : w1 = w2, 1n = 1n ;

x + : u 2 = U 0, v2 = 0, x : u1 = U 0, v1 = 0.

Линейный анализ приводит к дисперсионному соотношению ) ( M + 18 3 12 2 + 1 M h=. (3) M + 1 2 4 + 6 + 1 (18 + 12 + 1) 12 2 + 1 18 2 4 M 1 Рис. 3. Зависимость темпа роста возмущений от длины волны: кривая 1 – с учетом инерции жидкости, для M = 84, D = 1, Re = 1.2. Кривые 2–4 по решению (3) с учетом вязкости в плоскости ячейки для M = 3, 9, соответственно.

Зависимость h = h ( ) для M = 3, 9, 84 показана на рис. 3 кривыми 2–4. Учет вязкости в плоскости ячейки привел к стабилизации коротковолновых возмущений. Существует критическая длина волны возмущения *, интенсивность роста которой наибольшая. Ее зависимость от отношения вязкостей показана на рис. 4. Увеличение отношения вязкостей M не приводит к увеличению ширины пальцев. Она оценивается как 1.25.

Эта оценка и тенденция, чем больше зазор ячейки, тем шире образующиеся пальцы, качественно согласуется с реальной шириной пальцев на начальной стадии их образования (Smirnov N.N., Nikitin V.F., Ivashnyov O.E., etc., 2002– 2004).

Рис. 4. Зависимость критической длины волны от отношения вязкостей.

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

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

Пусть установилась структура фронта вытеснения из чередующихся пальцев, заполненных вытесняющей жидкостью «1», и промежутков между ними с вытесняемой жидкостью «2». Жидкость в пальцах движется со скоростью U 1, жидкость в промежутках – с меньшей скоростью U 2. Таким образом, возле границы раздела реализуется сдвиговой поток, когда скорости обеих жидкостей ей параллельны (рис. 6, а).

Рис. 5. Постановка задачи о сдвиговом течении двух вязких жидкостей в ячейке Хеле-Шоу.

Рис. 6. Профили скоростей: а – без учета сил вязкости, б – с учетом вязкости в плоскости пластин ячейки.

На первом этапе возникающая неустойчивость анализируется в рамках классической модели. Решается задача на малые возмущения в бесконечной области: сверху менее вязкая жидкость («1»), снизу – более вязкая («2», рис. 6, а). На границе ставятся условия равенства давлений и равенства между нормальной к границе скоростью и соответствующей субстанциональной производной вдоль нее. На бесконечности – условия постоянства скоростей.

12 µ 12 µ p p u v = 2 u, = 2 v, + = 0;

x y x y d d y = ( x, t ) : p1 = p2, v1 = 1, v2 = 2 ;

dt dt y + : u1 = U 1, v1 = 0, y : u 2 = U 2, v2 = 0.

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

[ ] [ ] Re h ( ) = hx ( ) = 0, Im h ( ) = hy ( ) =.

( M + 1) Таким образом, классическая модель не позволяет описать экспериментально наблюдаемую неустойчивость боковой поверхности вязких пальцев.

Высказывается гипотеза, что она аналогична неустойчивости Кельвина – Гельмгольца (Ламб Г., 1947), которая развивается за счет инерции жидкости.

Добавление этих сил в классическую модель приводит к следующему дисперсионному соотношению 6 M +1 D hx = + 2 + 2 +, (4) Re D + 1 M 18 ( M + 1) D 6 M 1 ( D M )D 2 D M = 2 + 2 2, =, ( D + 1) M Re ( D + 1) M Re M ( D + 1) 2 M U D = 2, Re = 2 1.

1 µ hx = hx ( ) Зависимость показана на рис. 7 кривой 1 при M = 84, D = 1, Re = 1.2. Возмущения бесконечно малой длины неограниченно возрастают, следовательно, боковая поверхность пальцев начинает разрушаться сразу же при их формировании. Дестабилизирующим фактором, приводящим к развитию неустойчивости данного типа, являются малые силы инерции. По схожести механизмов обнаруженная неустойчивость аналогична неустойчивости Кельвина – Гельмгольца (Ламб Г., 1947).

Рис. 7. Зависимость темпа роста возмущений от длины волны: кривая 1 – по решению (4) без учета вязкости в плоскости ячейки для Re = 1.2. Кривые 2– по неявному решению с учетом вязкости в плоскости ячейки для Re = 0.1, 0.7, 1.2 соответственно. Во всех расчетах M = 84, D = 1.

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

Эти силы добавляются в модель. На границе раздела жидкостей ставятся условия равенства векторов скоростей и векторов напряжений. На бесконечности – условия постоянства скоростей.

2u 2u 12µ u u p u = 2 u + µ 2 + 2, + u + v + t y x x y x 2v 2v 12 µ v v p v = 2 v + µ 2 + 2, + u + v + t y y x y x u v + = 0;

x y y = ( x, t ) : w1 = w2, 1n = 2 n ;

y + : u1 = U 1, v1 = 0, y : u 2 = U 2, v2 = 0.

Добавление вязких сил приводит к появлению стационарного экспоненциального профиля скорости (рис. 6, б), поэтому невозмущенное решение зависит от координаты y. Линейный анализ приводит к дисперсионному соотношению в неявном виде: f (, h ) = 0, f – комплексная функция. Вместо нее вводится действительная функция F = [ Re ( f )] + [ Im( f )], принимающая нулевые значения при тех же 2 значениях параметров. Корни действительной функции ищутся с привлечением, Re, M, D средств MATLAB. Параметры фиксируются, строится F (hx, h y ) = 0, h = hx + ih y, минимумы (нули) трехмерная поверхность которой определяются методом наискорейшего спуска. При всех рассмотренных значениях параметров существует только один минимум с положительной действительной частью.

Его зависимость от длины волны показана на рис. 7 при Re = 0.1, 0.7, 1.2 кривыми 2–4 соответственно. Видно, что учет вязкости в плоскости ячейки привел к стабилизации коротковолновых возмущений.

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

считается, что палец разрушается наиболее быстро растущим возмущением.

Для условий экспериментов Smirnov N.N., Nikitin V.F., Ivashnyov O.E., etc., 2002–2004, это несколько миллиметров, что намного меньше длины ячейки.

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

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

В главе IV вытеснение вязкой жидкости из ячейки Хеле-Шоу при больших числах Пекле моделируется методом сквозного счета. Движение обеих жидкостей описывается единой системой уравнений Навье – Стокса несжимаемой жидкости с переменной вязкостью, зависящей от объемной концентрации : = 0 соответствует вытесняющей жидкости («1»), = 1 – вытесняемой («2»). Зависимость вязкости от концентрации предполагается степенной: µ ( ) = µ1 (1 ) + µ 2. В расчетах использовались = 1 и = 1.

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

u v +u +v = 0, + = 0, (5) t x y x y u u v 12µ u u p u x = 2 u + x 2µ x + y µ y + x, + u +v + (6) t x y u v v v p 12µ v v = 2 v + µ + + 2µ.

+ u + v + (7) t y y x y x y y x Уравнения (6)–(7) решаются методом прогонки по неявной схеме, а уравнения (5) – по явной. Для тестирования метода проводится сравнение с аналитическим решением задачи о стационарном течении вязкой жидкости в прямоугольном канале между двумя пластинами. Получено полное соответствие.

На рис. 8 и рис. 9 представлены результаты двух расчетов, отличающихся только величиной зазора между пластинами ячейки: на рис. 8 = 1.2 мм, на рис. 9 = 3.7 мм. Остальные параметры одинаковы и соответствуют экспериментам Smirnov N.N., Nikitin V.F., Ivashnyov O.E., etc., 2002–2004.

Габариты ячейки: длина – 20 см, ширина – 10 см, отношение вязкостей M = 84, скорость вытеснения U 0 = 2 см/c. Результаты экспериментов даны в два момента времени. В начальный момент на границе между жидкостями задается полигармоническое возмущение типа sin (1/y ), одинаковое в обоих случаях. На входе ячейки поддерживалось постоянное давление Pвход. Заметно, что на начальном этапе характерный размер пальцев пропорционален зазору между пластинами. Изменение отношения вязкостей и скоростей вытеснения мало влияет на картину образования пальцев. Характерная ширина диктуется нелинейными эффектами.

Рис. 8. Результаты численного моделирования для = 1.2 мм : а – t = 1.5 c, б – t = 4c.

Рис. 9. Результаты численного моделирования для = 3.7 мм : а – t = 3 c, б – t = 5c.

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

Выводы Осреднением общих уравнений Навье – Стокса по зазору ячейки Хеле-Шоу 1.

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

Выявлен новый стабилизирующий фактор, определяющий структуру 2.

фронта неустойчивого вытеснения, – силы вязкости в плоскости ячейки.

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

Показано, что инерция вытесняющей жидкости приводит к развитию 3.

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

Проведено моделирование вытеснения вязкой жидкости из ячейки Хеле 4.

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

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

образующихся при вытеснении жидкости из ячейки Хеле-Шоу // Вестн.

Моск. Ун-та. Сер. 1, Матем. Механ. 2011. № 2. 40– 2. Logvinov O. A., Ivashnyov O. E., Smirnov N. N. Evaluation of viscous fingers width in Hele-Shaw flaws // Acta Astronautica 2010. 67. 53– 3. Звягин А. В., Ивашнев О. Е., Логвинов О. А. О влиянии малых параметров на структуру фронта неустойчивого вытеснения вязкой жидкости из ячейки Хеле-Шоу // Изв. РАН. Механ. жидкости и газа.

2007. № 4. 27–37.

4. Логвинов О. А. Особенности неустойчивого вытеснения вязкой жидкости из ячейки Хеле-Шоу при больших числах Пекле // Мат. 11ой межд. шк.-сем.

«Модели и методы аэродинамики». Евпатория. 2011. 5. Ивашнев О. Е., Логвинов О. А. Вытеснение вязких жидкостей из ячейки Хеле Шоу // Тез. докл. науч. конф. «Ломоносовские чтения». Секц. Механ.

Москва. 2009. Логвинов О. А., Ивашнев О. Е. Исследование устойчивости вытеснения 6.

вязких жидкостей из ячейки Хеле-Шоу // Тез. докл. междунар. конф.

«Современные проблемы газовой и волновой динамики», посвященной столетию со дня рождения академика Х. А. Рахматулина. Москва. 2009. 7. Ивашнев О. Е., Логвинов О. А. Неустойчивость при скоростном вытеснении вязкой жидкости из ячейки Хеле-Шоу // Матер. конф. «Современные фундаментальные и прикладные задачи математики и механики» в рамках всероссийского конкурса инновационных проектов. Москва. 2009. 8. Ивашнев О. Е., Логвинов О. А. Неустойчивость при фильтрационных течениях вязких жидкостей // Тез. докл. науч. конф. «Ломоносовские чтения». Секц. Механ. Москва. 2008. 9. Логвинов О. А., Ивашнев О. Е. Неустойчивое вытеснение вязкой жидкости из ячейки Хеле-Шоу // Тез. докл. всеросс. конф. «Современные проблемы механики сплошной среды», посвященной столетию со дня рождения академика Л. И. Седова. Москва. 2007. 114– Ивашнев О. Е., Логвинов О. А. Неустойчивое вытеснение вязкой жидкости 10.

из ячейки Хеле-Шоу // Тез. докл. науч. конф. «Ломоносовские чтения». Секц.

Механ. Москва. 2007.