Расчет отраженных ударных волн в задаче о схлопывании пустой полости Текст научной статьи по специальности «Математика»

index for surgical treatment of intracranial aneurysms // J. of Neurosurgery. 1989. V. 45. P. 119-129.

32. Sforza D., Putman C., Cebral J. Hemodynamics of cerebral aneurysms // Annu. Rev. Fluid Mech. 2009. V. 41. P. 91-107.

33. Tateshima S. et al. Three-dimensional blood flow analysis in a wide-necked internal carotid artery-

ophthalmic artery aneurysm // J. Of Neurosurgery. 2003. V. 99. P. 526-533.

34. Torii R. et al. Fluid-structure interaction modeling of aneurismal conditions with high and normal blood pressures // Computational Mechanics. 2006. V. 38. P. 482-490.

УДК 532.5:533.6.011.5

РАСЧЕТ ОТРАЖЕННЫХ УДАРНЫХ ВОЛН В ЗАДАЧЕ О СХЛОПЫВАНИИ ПУСТОЙ ПОЛОСТИ

В.С. Кожанов

Саратовский государственный университет, кафедра вычислительного эксперимента в механике E-mail: [email protected]

Рассматривается автомодельная задача о захлопывании пустой цилиндрической или сферической полости в сжимаемой жидкости с отношением удельных теплоемкостей y. Обсуждаются два возможных варианта течения после схлопывания, связанные с поведением энтропии при переходе через отраженную ударную волну. Проведенные расчеты показывают, что основные отличия в поведении параметров течения на стадии отражения носят количественный характер. Степень сжатия отраженной ударной волны, характеризуемая отношением р2/pi, убывает при увеличении y для обоих вариантов отражения.

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

Calculation of Outgoing Shock Waves in The Empty Cavity Collapse Problem

V.S. Kozhanov

Saratov State University,

Chair of Computational Experiment in Mechanics E-mail: [email protected]

The self-similar problem about a collapse of an empty cylindrical or spherical cavity in compressible fluid with adiabatic exponent Y is considered. Two possible variants of the flow after collapse are discussed. The variants are connected with the entropy behavior through the outgoing shock. The calculations show that the main difference in the flow quantities behavior at reflection stage have a quantitative character. Outgoing shock compression ratio, characterized by relation p2 /pi, decreases for both variants of the reflection when Y is increase.

Key words: one-dimensional flows, self-similar flows, homentropic flows, cavitation, bubble, outgoing shock, reflected shock.

ВВЕДЕНИЕ

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

Интерес к классическим работам [1, 2] связан с публикацией [3]. Ее авторы заявляют об осуществлении эксперимента, в котором, возможно, наблюдались признаки термоядерных реакций внутри схлопывающихся пузырьков газа в жидкости. Выполненный эксперимент напоминает опыты по изучению сонолюминисценции. Ультразвук с частотой 19.3 кГц воздействовал на жидкий ацетон, в котором атомы водорода были замещены на атомы дейтерия. Синхронно с ультразвуком ацетон С3В6О облучался импульсами высокоэнергетических нейтронов, что является новым элементом по сравнению с предшествующими экспериментами. Нейтроны инициировали зарождение в жидкости очень маленьких пузырьков, которые росли до относительно больших размеров и затем интенсивно схлопы-вались под действием ультразвука. О термоядерных реакциях свидетельствовало появление трития и нейтронов с энергиями 2.5 МэВ, которые должны сопутствовать синтезу трития из атомов дейтерия. Результат эксперимента стал большой неожиданностью, так как, согласно расчетам, при схлопывании пузырьков не может достигаться температура более 11000 К, а для термоядерного синтеза необходима температура на 2-3 порядка выше. По этой причине многие специалисты сомневаются в корректности

В.С. Кожанов. Расчет отраженных ударных волн в задаче о схлопывании пустой полости

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

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

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

1. Энтропия при переходе через отраженную УВ возрастает скачком до некоторого конечного значения и далее остается постоянной во всей области течения за УВ;

2. Энтропия при переходе через отраженную УВ не меняется и во всей области течения за УВ остается равной своему первоначальному значению.

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

Предполагается, что объекты, имеющие отношение к первой/второй модели отражения, наделены соответственно индексами I и II.

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

В работах [5, 6] при расчете стадии отражения используется другой подход. Полагается, что отраженная УВ обеспечивает скачкообразный рост энтропии. В отчете [5] расчет выполнен при значении 7 = 3 для случая сферически симметричного потока. Стоит отметить, что авторы [5] наряду со значением положения ударного перехода, дают некоторые гидродинамические характеристики УВ: коэффициент сжатия и отношение энтропий.

Статью Лазаруса [6] саму по себе можно считать новаторской, так как она представляет собой успешную попытку объединить в рамках одной постановки две задачи со сходящимися течениями: задачу о схождении УВ и задачу о схлопывании полости. В работе [6] расширяются результаты, данные в работе [5], на случай двух симметрий (цилиндрической и сферической) и большого числа значений показателя адиабаты 7. Однако среди параметров, характеризующих течение после схлопывания, фигурирует только положение отраженной УВ.

Целью данной работы является уточнение и обобщение результатов [1, 5, 6], касающихся течения после схлопывания. Основное внимание направлено на расчет таких гидродинамических параметров, характеризующих отраженную УВ и поток за ней, как коэффициент сжатия, перепад давлений, отношение энтропий, число Маха за УВ. Проводится сравнение зависимостей (от 7) перечисленных параметров для двух моделей отражения. Результаты представлены для случая сферически симметричного течения.

1. ПОСТАНОВКА ЗАДАЧИ

Чтобы описать процесс схлопывания пустой цилиндрической (V = 2) или сферической (V = 3) полости в среде с отношением удельных теплоемкостей 7, необходимо решить уравнения неразрывности

и сохранения количества движения:

dc2 dc2 , _ 2 _+„ + (7 - 1)с2

du , „,u —+ (v - 1)-dr r

— 0 du + u^u + 1 dc2 (i)

' dt dr y — 1 dr '

где £ — время, г — координата, и = и (г, £) — скорость частицы жидкости, с2 = с2(г, £) — квадрат скорости звука.

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

s —

р_ Р7

(2)

сохраняет свое значение не только вдоль траектории частицы, но и во всей области течения (для всех траекторий).

Плотность среды р = р(г, £) связана с квадратом скорости звука уравнением

2 Р

c — Y- . Р

Граничные условия на границе полости имеют вид

при r — Rc : u — Rc

c2 — 1,

(3)

(4)

где г = Яс(£) — уравнение движения границы.

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

Задача допускает автомодельное решение, которое будем искать в форме

— —Anr1-1/nF (n), c2 — A2 n2 r2-2/nG(n), П — At/r1/n, A — const.

(5)

Рис. 1. Плоскость rt

Здесь п — независимая автомодельная переменная, а ^ = ^(п) и С = С(п) — автомодельные представители (по г) скорости частицы жидкости и и квадрата скорости звука с2 соответственно, п — показатель автомодельности. Они описывают поведение и и с2 при фиксированном г в зависимости от £.

Подставляя (5) в (1) и разрешая полученные уравнения относительно производных, получим систему двух обыкновенных дифференциальных уравнений (ОДУ) относительно функций ^ и С:

dF dn

1 A

(Y — 1)n ' Ao '

dG

dn

G A2 1 + nF ' Ao '

(6)

Ах = (7 - 1)^п - 1)^Ао - Д5, Д2 = -2(1 - п)^До + Да, До = (1 + п^)2 - п2^, Д5 = (7 - 1)(1 + п^) [(V - 1)пп^ + V™ - 1] ^ + 2(1 - п)пС. Граничное условие (4) в автомодельных переменных принимает вид

при п = -1 : F(-1) = 1, С(-1) = 0.

Требование ограниченности и и с2 при г ^ налагает на F и С условие [1]

при n — 0 : (—n)1-nF(n) — 0, (—n)2-2nG(n) — 0.

2 —2n,

(7)

(8)

Таким образом, чтобы получить решение автомодельной задачи о схлопывании пустой полости для стадии схлопывания, необходимо решить краевую задачу для системы ОДУ (6) с граничными условиями (7) на границе полости и (8) на бесконечности.

u

B.C. Кожанов. Расчет отраженных ударных волн в задаче о схлопыванпп пустой полости

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

Исследование задачи на фазовой плоскости переменных V = —nF(n), Z = n2G(n) (рис. 2) сводится к изучению поведения интегральных кривых одного нелинейного ОДУ

- = (Y — 1)Z ¿2, dV W ; Д i'

Д1 = [(y — 1)vnV — 2(1 — п)]Д 0 + (1 — V )Д 5, Д 2 = 2пД 0 + Д 5, Д о = (1 — V )2 — Z, Д5 = (Y — 1)(v — 1)nV2 — [(y — 1)(vn — 1) + 2(1 — n)] V + 2(1 — n).

(9)

,z=(i-v):

-0,2

0,2 0,4 0,6

Рис. 2. Фазовая плоскость VZ

1,0 V

Уравнение (9) в общем случае имеет шесть особых точек в конечной части плоскости. Перечислим их:

1) точка Рь V = 1, И = 0;

2) точка Р4: V = 0, И = 0;

2

3) точка Р5: V =

(a + 2)n' (а + 2)2n2' 4-5) точки P2 и P3:

, Z = ,(Y, 0!2)а2, где а = (y — 1)v;

V = (—Bi ± vB2 — 4B2Bo) / (2B2), Z = (1 — V)2,

Bi = b + 2d, B2 = b +(y — 1)d, Bo = 2d, b = (y — 1)(vn — 1),

d =1 — n.

Эти точки являются образами предельной характеристики. Присвоим обозначение Р2 точке с меньшим значением переменной V;

6) точка с координатами V = 1/п, И = 0. Эта точка не участвует в решении.

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

удаленная особая точка

6 •

V=

2(1 — n) (Y — 1)vn'

Z =

(10)

Подробный качественный анализ ОДУ (9) и перечисленных особых точек выполнен в работах [2, 6].

Граничным условиям (7), (8) отвечают соответственно точки фазовой плоскости Р1 и Р4. В настоящей работе рассматривается решение задачи, когда полость обладает цилиндрической или сферической симметрией для двух [2, 6] интервалов значений показателя адиабаты: 71 < 7 < 72

и 72 <7 .В этих диапазонах существует так называемое физическое решение задачи. Это решение характеризуется тем, что при требуемом значении показателя автомодельности п интегральная кривая уравнения (9), соединяющая точки и Р4, проходит через узел Р3, если 71 < 7 < 72, или узел Р2, если 72 < 7, вдоль уса отдельного направления.

Величины 71 и 72 являются трансцендентными числами и определяются численно. Их значения

Таблица 1 приведены в табл. 1.

Для каждого конкретного набора параметров

(V, 7) показатель п выбирается из некоторого промежутка [птщ,птах]. В этом случае точки Р2 и Р3 будут узлами. При птщ = птах точки Р2 и Р3 совпадают.

Метод определения показателя автомодельности заключается в отыскании интегральной кривой уравнения (9), которая при искомом п выйдет нужным образом из особой точки Р3 (если 71 < 7 < 72) или Р2 (если 72 < 7) и попадет в особую точку Рх, т.е. будет второй сепаратрисой седла Рх и, следовательно, решением. Чтобы определить нужный показатель автомодельности будем численно интегрировать уравнение (9), выбирая п из интервала приемлемых значений [птт,птах], соответствующих заданному набору параметров (V, 7). Координаты особой точки, из которой выходит интегральная кривая, брать в качестве начальных условий интегрирования нельзя. В малой окрестности особой точки решение можно представить в виде ряда

Z = Z * + г г 5\ V = V * + 5<< 1, (11)

где (V— координаты одной из особых точек Р2, Р3 (той, которая учавствует в решении).

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

Обозначим пс некоторое текущее значение показателя автомодельности. Установлено, что если п < пс < птах, то получаемая интегральная кривая проходит слева от искомой интегральной кривой. Если же птт < пс < п, то получаемая интегральная кривая отклоняется вправо. Эти закономерности служат ориентирами при поиске требуемого п.

Вычисления необязательно вести на фазовой плоскости VZ. Описанный метод можно перенести на систему ОДУ (6) и исследовать поведение ее решения в зависимости от п. Для системы ОДУ (6) определение искомого значения п означает возможность аналитического прохода функций F, С через предельную характеристику (образами которой являются особые точки Р2 и Р3 на плоскости VZ). Условие на границе полости п = -1 при этом также выполняется (это соответствует приходу интегральной кривой в седло Рх на плоскости VZ). Аналогично, чтобы определить нужный показатель автомодельности, будем численно интегрировать систему ОДУ (6), выбирая п из интервала приемлемых значений [птт,птах], соответствующих заданному набору параметров (V,7). Начальные условия находятся из тех же соображений, что и (11).

2. ТЕЧЕНИЕ ПОСЛЕ СХЛОПЫВАНИЯ

В настоящей работе рассматриваются две модели отражения. Для обеих моделей справедливо описание, которое представлено ниже.

Стадия схлопывания завершается фокусировкой полости в центре в момент времени £ = 0. Распределения скорости, квадрата скорости звука, плотности и давления по координате в этот момент имеют вид

1—п 2(1-п) _ 2(1_п) _ 2(1-п)т

и (г, 0) ~ г_ п , с2 (г, 0) ~ г_ п , р(г, 0) ~ г (7_1)п , р(г, 0) ~ г _1)п . (12)

Решение, описывающее схлопывание и действительное для £ < 0, можно продолжить и на положительные значения времени £ > 0, то есть из области отрицательных в область положительных п. На стадии схлопывания была проинтегрирована система (6) для F и С в интервале -1 < п < 0. Поскольку F и С конечны при п = 0 и уравнения системы (6) не имеют особенности в этой точке, интегрирование можно продолжить и в область положительных значений п. При продолжении решения в область п > 0 возникает особенность при некотором п = пн, так как при этом выполняется

Значения 71 и 72 для двух симметрий течения

V 71 72

2 2.9776 16.8693

3 2.4056 8.4635

В.С. Кожанов. Расчет отраженных ударных волн в задаче о схлопыванпп пустой полости

равенство До = 0. Эту особенность, которая соответствует предельной линии, выходящей из центра О на плоскости гЬ, следует устранить. Для этого необходимо ввести разрыв в величинах Р и С при некотором значении п = nJ, меньшем значения пн. Физически это означает, что в точке г = 0, Ь = 0 образуется отраженная УВ, которая распространяется от центра навстречу набегающему на нее потоку. Поток впереди УВ является гладким продолжением течения, найденного на стадии схлопывания. Следовательно, значение показателя автомодельности п в автомодельном решении, описывающем течение за отраженной УВ, должно быть тем же самым, что и в решении до момента фокусировки, чтобы обеспечить то же самое поведение (12) при Ь = 0. Таким образом, закон распространения отраженной УВ имеет вид nJГ1/n = АЬ. Формально интегральная кривая ОДУ (9) на плоскости VZ проходит через начало координат и продолжается в область V < 0, И > 0 до тех пор, пока она не пересечет параболу И = (1 — V)2 в точке Н, отвечающей значению пн. Этот интеграл теряет физический смысл за точкой Н, так как функции Р и С должны быть однозначными функциями п. Поскольку отраженная УВ выходит из точки О на плоскости гЬ, не существует ни одной характеристики, которая выходила бы из точки О, и поэтому выражение Д0 меняет знак при п = пJ, испытывая при этом разрыв. Соответственно этому интегральная кривая на плоскости VZ перескакивает через параболу И = (1 — V)2, также претерпевая разрыв.

Значение nJ будем называть положением ударного перехода. Для каждой из двух рассматриваемых моделей отражения при одном и том же наборе параметров (V, 7) это положение различно. Различие связано с тем, что для моделей отражения используются разные соотношения на УВ, которые обусловлены, в свою очередь, выбором поведения энтропии на УВ.

На плоскости VZ точка 1, отвечающая значению пJ, характеризует автомодельные представители, которые принадлежат течению непосредственно перед УВ.

Найдем вид течения в области nJ < п < за отраженной УВ. Эта область имеет нулевые размеры при Ь = 0. Течение за отраженной УВ, как отмечено ранее, является автомодельным с прежним показателем автомодельности. Решение удобно искать в той же форме (5), что и на стадии схлопывания. При этом остается справедливой система ОДУ (4).

В статье [1] показано, что для течения за отраженной УВ существует единственно верная возможность: полость остается захлопнутой при Ь > 0. Это обстоятельство определяет граничные условия при п ^ то есть в центре схлопывания. Скорость частицы и должна быть равна нулю, то есть и(0, Ь > 0) = 0, а квадрат скорости звука с2(0, Ь > 0) должен принимать положительное значение. Следовательно, граничные условия для функций Р и С имеют вид

при п ^ : Р(+го) = 0, = сопб^ (13)

Таким образом, для того, чтобы получить автомодельное решение в области за отраженной УВ, необходимо проинтегрировать систему (6) для функции Р и С с граничными условиями при п = пJ на отраженной УВ и (13) в центре схлопывания.

На плоскости VZ течение за отраженной УВ описывается интегральной кривой ОДУ (9), выходящей из бесконечно удаленной особой точки типа седло Р6 и заканчивающейся в некоторой точке К, лежащей внутри области И > (1 — V)2. Точка К, как и точка 1, отвечает значению nJ, однако характеризует автомодельные представители, которые принадлежат течению непосредственно за УВ, т.е. точки 1 и К являются образами берегов отраженной УВ.

Окончательно, полное решение для течения после схлопывания на плоскости VZ состоит из (О) — участка Р41 кривой, прошедшей через начало координат и лежащей целиком в области И < (1 — V)2,

(b) — кривой, соединяющей точку Р6 с точкой К и лежащей целиком в области И > (1 — V)2,

(c) — ударного перехода, соединяющего части (а) и (Ь) в соответствии с соотношениями на УВ.

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

ния на плоскость VZ используются участки одних и тех же интегральных кривых. Отличие состоит в том, что точки 1 и К имеют на соответствующих кривых разное местоположение. Это следствие того, что разные модели отражения имеют для общего набора параметров (V, 7) разные положения ударного перехода nJ.

Определим вид соотношений при переходе через УВ. Для каждой из двух рассматриваемых моделей отражения эти соотношения различны и учитывают заложенное в модель поведение энтропии.

Обозначим параметры течения перед фронтом УВ индексом 1, а за фронтом — индексом 2. Тогда уравнения сохранения массы и количества движения для УВ можно записать в виде

Р2(и2 - Я) = Р1 (и1 - Я), Р2 [с2 + 7(и2 - Я)2] = Р1 [с2 + 7(и1 - Я)2] . (14)

Если рассматривается первая модель отражения, когда энтропия возрастает скачком на УВ, уравнения (14) необходимо дополнить уравнением сохранения энергии на УВ

2с2 + (7 - 1)(и2 - Я)2 = 2с1 + (7 - 1)(и1 - Я)2. (15)

Если же предполагается, что энтропия при переходе через УВ не меняется, то есть во всем течении справедливо единственное соотношение между давлением р и плотностью р (вторая модель отражения), то к уравнениям (14) нужно добавить соотношение (2).

В уравнениях (14) и (15) Я — скорость распространения УВ. Учитывая (5), выразим скорость Я в автомодельных переменных

Я = ^ = пг = кпг1_1/п п_1.

Автомодельные представления уравнений (14), (15) и соотношения (2) имеют соответственно вид

Д2(1 + nF2) = Д1 (1 + п^), Я2 (п2С + 7(1 + пF2)2) = (п2С + 7(1 + пFl)2) ; (16) 2п2С + (7 - 1)(1 + nF2)2 = 2п2С + (7 - 1)(1 + п^)2; (17)

Я = (п2С)1/(7_1) . (18)

При получении (18) было использовано выражение (3).

Таким образом, для определения значений параметров течения за УВ, когда известны значения параметров перед УВ, можно составить два набора условий. В случае первой модели отражения необходимо использовать уравнения (16) и (17), а в случае второй модели отражения — уравнения (16) с учетом (18). Течение для всего двухстадийного процесса схлопывания-отражения определяется решением системы двух ОДУ (6). Поэтому для осуществления ударного перехода через отраженную УВ достаточно только двух соотношений, которые должны связывать значения F и С перед и за УВ. С помощью преобразований получаем

1+ nF2 = 1-1(1 + п*)2 + ^ С2 = С1 + I-1 .(^^^(Цп^; (19)

7 + 1Ч ' у 7 + 11+ п^' 2 п2

^2/(7_1) (1 + nF2) = ^1/(7_1)(1 + nFl),

ч1/(7_1) ^2п I I „р ^2) _ гч1/(7_1) („2п \ (Л I ^Т? \2)

(20)

(п2^2 + 7(1 + nF2)2) = С^7"^ (п2^1 + 7(1 + п^)2) .

Системы уравнений (19) и (20) представляют собой окончательные варианты условий на отраженной УВ для первой и второй моделей отражения соответственно. Соотношения (19) и (20) должны выполняться при п = nJ. При этом важно еще раз отметить, что положение ударного перехода nJ различно для каждой из рассматриваемых моделей течения после схлопывания.

Максимальную вычислительную сложность задачи представляет определение значения показателя автомодельности п, которое проводится на первом этапе решения. Основные вычисления при построении решения после схлопывания связаны с нахождением положения ударного перехода nJ, но эти вычисления имеют меньшую трудоемкость.

3. РЕЗУЛЬТАТЫ РАСЧЕТОВ

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

Расчеты проведены для случая сферически симметричного течения (V = 3). Исследованы все диапазоны значений показателя 7, для которых возможно построение аналитического решения. Результаты расчетов представлены в виде графиков.

В.С. Кожанов. Расчет отраженных ударных волн в задаче о схлопывании пустой полости

На рис. 3 представлены распределения автомодельных представителей скорости, плотности, квадрата скорости звука и энтропийной функции на обеих стадиях течения для двух значений показателя адиабаты: 7 = 3 и 7 = 7. Автомодельные представители обладают разным поведением для первой и второй модели отражения (рис. 3, кривые I и II соответственно). Однако эти различия в большей степени носят количественный характер.

1,0-1 ^

1,51

5 у=7

Н

у=3

-L-J-L

0,5-

II П

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

в

г

Следует отметить, что приводимые автомодельные функции при изменении независимой автомодельной переменной п от —1 до описывают поведение соответствующего размерного параметра от времени Ь при фиксированном значении координаты г.

Значения автомодельных представителей Р и С на берегах отраженной УВ для четырех значений 7 (V = 3) сведены в табл. 2.

Таблица 2

Параметры Модель Точка П ) С

7 = 3.0 п = 0.70856126 I 31 0.613874 0.425515 1.35655

I К1 0.613874 -0.271603 3.73506

II 3п 0.670408 0.375532 1.47137

II Кп 0.670408 -0.245560 3.30371

7 = 7.0 п = 0.55521127 I 31 0.222869 0.348165 6.15235

I К1 0.222869 -0.542505 29.6114

II 3п 0.263572 0.275848 7.04475

II Кп 0.263572 -0.445235 22.6999

7 = 20.0 п = 0.47708490 I 31 0.0755734 0.342964 28.9678

I КI 0.0755734 -0.746679 298.737

II 311 0.106180 0.256646 35.7991

II Кп 0.106180 -0.498699 167.747

7 = 200 п = 0.42443812 I 31 0.00714694 0.314150 890.776

I К1 0.00714694 -1.01801 37890.2

II 3п 0.0159785 0.242132 1207.77

II Кп 0.0159785 -0.391378 9075.54

На рис. 4-5 показаны кривые, характеризующие зависимости (от 7) положения ударного перехода nJ, коэффициента сжатия р2/р1, перепада давления р2/р1, отношения энтропий 52/^1 и числа

б

Рис. 4. Изменения (в зависимости от 7) а — положения ударного перехода ЦJ, б — коэффициента сжатия

В.С. Кожанов. Расчет отраженных ударных волн в задаче о схлопывании пустой полости 6 *

Р\

ЗАКЛЮЧЕНИЕ

4-

у=2.4056

14

20 У

2,0-

1,6

1.2-

0,8-1

у =2.4056

у= 8.4635

14

20 Y

0,15-1

0,10-

0,05-

14

20 Y

Рис. 5. Изменения (в зависимости от 7) а — перепада давлений, б — отношения энтропий, в — числа Маха после перехода через УВ

Каждое из решений двухпараметрического семейства (V, 7), описывающее схлопывание полости и справедливое для £ < 0, может быть закончено для £ > 0. При этом решение для положительных моментов времени зависит от того, каков закон изменения энтропий на отраженной УВ. Однако основные отличия в поведении газодинамических параметров носят больше количественный, чем качественный характер.

Автор благодарит И.А. Чернова за внимание к работе. Библиографический список

1. Хантер К. О захлопывании пустой полости в воде // Механика: Периодический сб. переводов иностр. ст. 1961. № 3 (67). С. 77-100.

2. Брушлинский К.В., Каждан Я.М. Об автомодельных решениях некоторых задач газовой динамики // УМН. 1963. Т. 18, вып. 2 (110). С. 3-23.

3. Taleyarkhan R.P., West C.D., Cho J.S. et al. Evidence for Nuclear Emissions During Acoustic Cavitation // Science. 2002. V. 295, № 5561. P. 1868-1873.

4. Lord Rayleigh. On the Pressure Developed in a Liquid During the Collapse of a Spherical Cavity // Phil. Mag. 1917. V. 34, № 200. P. 94-98.

а

б

в

5. Richtmyer R.D., Lazarus R.B. Singularity Fitting in Hydrodynamical Calculations II. Los Alamos Sc. Lab. Report LA-6108-MS, 1975. 19 p.

УДК 539.3

ИЗГИБ СОСТАВНОЙ АНИЗОТРОПНОЙ ПЛИТЫ НОРМАЛЬНОЙ НАГРУЗКОЙ

В.И. Копнина, Е.Ю. Крылова*

Саратовский государственный университет, кафедра математической теории упругости и биомеханики; "Саратовский государственный технический университет, кафедра математика и моделирование E-mail: [email protected]

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

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

6. Lazarus R.B. Self-Similar Solutions for Converging Shocks and Collapsing Cavities // SIAM J. Numer. Anal. 1981. V. 18, iss. 2. P. 316-371.

Bend of Composite Anisotropic Slab Under Normal Loading V.I. Kopnina, E.Yu. Krylova*

Saratov State University,

Chair of Mathematical Theory of Elasticity and Biomechanics; * Saratov State Technical University, Chair of Mathematics and Modelling E-mail: [email protected]

In this contribution, deflected mode of a thin slab under bending is investigated by the Lehnicky method of complex potential. The slab is composed of two elliptical rings; they are embedded in each other without tension. Material of the rings is anisotropic and different.

Key words: elliptical thin slab,, anisotropic, bend, complex potential by Lehnicky.

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

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

Обозначим внешний контур плиты [кольца] ¿00) [¿01) ], внутренний контур плиты [кольца] — ¿10) [I]; полуоси контурных эллипсов — аП^, бПТ) (п = V = 0.1). Здесь и далее величины, характеризующие упругое состояние плиты, имеют индекс (0) сверху, а кольца — (1).

Задачу об изгибе такой плиты можно привести [2] к определению комплексных потенциалов С.Г. Лехницкого Ш3(V) ^) (здесь и далее ] = 1, 2) из граничных и контактных условий [3]:

2Re ± W< v)' (j >) = — ^

j = 1

2 dW ( v)

2Re£ j W( V) (tjv)) = — ^W^ на L (0)(v = 0), l(v = 1); (1)

2

2Re£ [dg'wf (j) — jwf' (*<»)] = ft (k =14 на if [4«] , (2)

3=1

где «(V) = -(V) = , Ы .(V) = Р^ .(V) = (V). . = дШ(0) + дШ(1) . = дШ(0) + дШ(1)

где «31 = I, 3 = ^ , Й33 = -ру, Й34 = Яз ; X1 = — + , /2 = о^ + ,

/3 = — Сх + С1, /4 = 6Я (ж, у) ^20) — v21)) х + Су + С2. Здесь ьЗ^ — аффиксы точек на контурах соответствующих эллипсов, расположенных в областях изменения ^3^, которые получаются из области, занятой срединной плоскостью составной плиты, при помощи известных афинных преобразований [2];