CC BY
11
В момент времени 0.001848 с в выходном сечении диффузора скорость движения продуктов детонации становится отрицательной, т. е. воздух начинает поступать через это сечение в двигатель. Расчёт проводился до этого момента времени. Из табл. 2 видно, что давление на дне камеры в этом промежутке времени монотонно уменьшается, давление в среднем сечении диффузора сначала уменьшается, затем его изменение принимает колебательный характер.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Сапунков Я.Г., Шиндяпин Г.П., Поршнев H.A., Федорец О.П. Математическая модель детонационного двигателя // Математика. Механика: Сб. науч. тр. Саратов: Изд-во Сарат. ун-та, 2001. Вып. 3. С. 178 - 181.
2. Сапунков Я.Г., Шиндяпин Г.П., Поршнев H.A., Федорец В.И. Вычислительный эксперимент по определению газодинамических параметров в камере детонационного двигателя // Математика. Механика: Сб. науч. тр. Саратов: Изд-во Сарат. ун-та, 2002. Вып. 4. С. 213- 216.
УДК 232.5; 232.135 Л. И. Сафрончик, М. И. Сафрончик
НЕУСТАНОВИВШЕЕСЯ «ЗАПАЗДЫВАЮЩЕЕ» ТЕЧЕНИЕ КУЭТТА ВЯЗКОПЛАСТИЧНОЙ СРЕДЫ МЕЖДУ ПАРАЛЛЕЛЬНЫМИ СТЕНКАМИ
Рассматривается неустановившееся течение вязкопластичной среды между параллельными плоскостями, одна из которых остаётся неподвижной, а другая начинает двигаться из состояния покоя с некоторой скоростью.
В отличие от вязкой и обычной Бингамовской жидкостей течение возникает лишь в области, непосредственно примыкающей к движущейся пластине, а остальная часть среды остаётся неподвижной. Область течения растёт, охватывая всё новые слои жидкости, пока не займёт всё пространство между пластинами.
Решение этой задачи связано с определёнными трудностями, так как она относится к классу не вполне корректно поставленных задач (область течения в начальный момент отсутствует и негде задавать начальное условие), а граница области течения изменяется во времени по неизвестному закону.
Для решения используется метод Колоднсра [1], по-видимому, единственный из известных на сегодня методов решения подобных задач.
Считая площадь пластин достаточно большой, течение можно рассматривать как плоскопараллельное. Для нахождения единственной отличной от нуля компоненты скорости необходимо решить следующую краевую задачу:
р^г^тт" (0<у<ко, о<t<n (1)
Ot dyi
Vx(0,t)=U(t), (2)
Hm = (3)
y-*h(l)-Q ay tl
lim Vx = 0, (4)
y->h(t)~0
A( 0) = 0, А(Г) = Я. (5)
Здесь p — плотность, г) - структурная вязкость, тс и тд - статическое и динамическое предельные напряжения сдвига, h(t) — нижняя граница зоны течения, Т- момент времени, когда вся область между пластинами придёт в движение. Начало декартовой системы координат расположено на движущейся пластине, а оси координат направлены вдоль движущейся пластины и внутрь слоя.
Следуя методу Колоднера, распространим определение решения на область h(t)<y <оо, полагая там Vx = 0. Сначала строится вспомогательное решение при дополнительных условиях:
к,(О,О=ф(0. (6)
Kr(y,0) = Vr(®,0 = 0, (7)
Ijm У*(У,О" lim Vx(y,t) = 0, (8)
Нга Hm = (9)
у-^А(|)+о су у-»А(1)-0 ay ч
Для единственности решения потребуем, чтобы оно было ограни-
3V,
ценным, а виде
ду
<™. Решение этой вспомогательной задачи записывается в
Ф
t
vx{y,t)=№Wc
о
(10)
ф(0(11) 2т1 Vn0 Jt-Z
Здесь U(t) - скорость пластины, v = —.
Р
Потребовав теперь, чтобы
lim Vx(y,t)=0, (12)
>->А(<)+0
ÖV
lim £Г»= 0, (13)
у->А(»)+о ду
получим два уравнения для определения неизвестной границы /?(г). В работе [2] показано, что если /г(/) является решением уравнения (12), то Ух = О в области /г(г)<у <ог! и, следовательно, уравнение (13) выполняется автоматически. Выпишем уравнения (11) и (12) более подробно:
Уравнения (14) и (15) относительно функции h{t) являются сложными интефо-дифференциальными уравнениями. Эффективных методов их решения нет. Для их решения применён метод обратных задач, заключающийся в том, что задаётся закон изменения функции h(t) и ищется функция U(t), соответствующая этому закону. Относительно функции U(t) уравнения (14) и (15) являются линейными интегральными уравнениями Воль-тсрра. Функция h(t) ищется в виде кусочно-ломаной линии, угловые коэффициенты наклона которой на каждом временном шаге определяются из условия совпадения значения истинной скорости в конце этого шага со значением скорости, найденной при решении обратной задачи.
Для оценки погрешности предлагаемого метода построено точное решение одной обратной задачи. Пусть А(/) = 2a4vt (такой закон изменения границы присущ всем известным автомодельным задачам стефанов-ского типа). Уравнения (14) и (15) в этом случае примут вид
(16)
(17)
Здесь А, В, С, D, Е обозначают найденные значения интегралов:
о V1-*
■а " crfc(a)erfc(- о) j,
1 Г~ 2 l—Jx Г—
Е = \}~\*пе" ^dx = я e°2erfc{a)erfc(- a)- —erf{v).
о 0-*) °
Единственным решением уравнения (16) является функция U(t) = М ,
■ • тс-тй 2Я-ЛС „
где М=-£—--, а уравнение (17) выполняется автоматически. Если
2г| С
ввести характерную скорость U0, то скорость пластины будет представляться формулой
(18)
Постоянная а в законе изменения границы будет решением уравнения
~2 ~.....т (19)
хд
= Гр
2 Ье-
Распределение скоростей в зоне течения задается формулой
(20)
Результаты расчётов по предлагаемому методу будут обсуждаться в наших дальнейших работах.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Kolodner J.J. l;ree boundary problem for the heat equation wich applications of change of phase // Comm. On Pure and Appl.Math. 1956. Vol. IX, № 1.
2. Сафрончик Л.И. Некоторые задачи неустановившегося течения вязкопластич-ных сред: Дис. ... канд. физ.-мат. наук. Ростов н/Д. 1962. 109 с.
УДК 232.5; 232.135
М. И. Сафрончик
НЕУСТАНОВИВШЕЕСЯ "ЗАПАЗДЫВАЮЩЕЕ" ТЕЧЕНИЕ ВЯЗКОПЛАСТИЧНОЙ СРЕДЫ ПО НАКЛОННОЙ ПЛОСКОСТИ
Данная статья посвящена реализации идеи сведения задачи о неустановившемся течении вязкопластичной среды к системе обыкновенных дифференциальных уравнений, изложенной в работе [1]. В качестве примера рассмотрен переход от одного стационарного режима течения по наклонной плоскости к другому при изменении угла наклона этой плоскости к горизонту.
