Пиротехническая химия
Главная Начинающим пиротехникам Статьи Добавить статью Добавить материалы на сайт Поиск по сайту Карта книг Карта сайта
Книги в помощь
Военная история Изготовление и применение ВВ Пиротехника в военном деле Разное по пиротехнике Физика в пиротехнике Химия ВВ и составов
Новые книги
Суворов С. "Бронированная машина пехоты БМП -3 часть 1" (Военное дело)

Яковлев Г.П. "122 мм самоходная пушка образца 1944 г." (Военное дело)

Суарес Г. "Тактическое преимущество " (Военное дело)

Стодеревский И.Ю. "Автобиография записки офицера спецназа ГРУ " (Военное дело)

Семиколенков Н.П. "стрельба из танковых пулеметов " (Военное дело)
Взрывы и волны. Взрывные источники электромагнитного излучения радиочастотного диапазона - Прищепенко А.Б.
Прищепенко А.Б. Взрывы и волны. Взрывные источники электромагнитного излучения радиочастотного диапазона — М. : БИНОМ, 2008. — 208 c.
ISBN 978-5-94774-726-3
Скачать (прямая ссылка): vzriviivolni2008.djvu
Предыдущая << 1 .. 31 32 33 34 35 36 < 37 > 38 39 40 41 42 43 .. 61 >> Следующая

Область определения переменной s заранее известна — это полная масса жидкости, находящаяся в рассматриваемом объеме. Она не изменяется со временем. В эйлеровых координатах правая граница области определения переменна, причем закон ее зависимости от времени подлежит определению в процессе самого решения задачи.
Для решения системы дифференциальных уравнений (5.15) аппроксимируем ее полностью консервативной разностной схемой [5.9]:
pjt
= (,Д0,5)у(0,5) t rt =
(HH] f-1
{Sn) 5 S I P j
__„(0,5) (a)
+ fr'>
= _?(v).
p* 4n
^^(«)(,(0,5),(0,5)^ +q.
1
q = 2
\P*r J
e(v> +
/(+1)
4(0,5)
p* (+1)/-(+1)
e(v)(+l)
132 5. Излучающая ударная волна: схождение до последнего микрона
Здесь использованы следующие обозначения для сеточных функций:
Приняты следующие соответствия для сетевых величин: y(si9 tj) = у{; y(st ±К tj) = yj±l =y(±l);
y(si9tj+l) =yj+ =y; у
ґ h х іЛ M
1 2 J 2
h At
j;(g) =ay+(l-a)y.
Индекс і отвечает слою по времени, j — по пространству.
5.5.3.1. О методе раздельных прогонок для разностных схем магнитной гидродинамики
По сравнению с газовой динамикой система уравнений магнитной гидродинамики является более сложной, — здесь появляются дополнительные уравнения, описывающие электромагнитное поле; кроме того, в гидродинамических уравнениях появляются дополнительные члены. Поэтому для решения разностных схем в магнитной гидродинамике целесообразно использовать метод раздельных прогонок [5.9]. Естественно выделить уравнения электромагнитного поля в отдельную группу (магнитную часть) и рассматривать ее наряду с динамической и тепловой группами уравнений.
Схема метода раздельных прогонок для этого случая представлена на рис. 5.19. Кроме внутренних итераций в каждой из отдельных групп — динамической, магнитной и тепловой, а также внешних итераций, предусмотрены промежуточные итерации между группами уравнений (на рис. 5.19 такие промежуточные итерации указаны между магнитной и тепловой группами). Комбинируя в зависимости от характера задачи число внутренних, промежуточных и внешних итераций, удается достичь заданной точности за минимальное время. Электромагнитная сила, входящая в уравнение движения, и джоулев нагрев в уравнении энер-
5.5. Расчеты проведены.
133
Газодинамические уравнения
Уравнения МГД
Уравнения
теплопроводности
J+1
Итерации
Рис. 5.19. Схема метода раздельных прогонок для системы уравнений магнитной гидродинамики
гии вычисляются в магнитной части и во внутренних итерациях в группах I и III не участвуют. Точно так же скорость, плотность и электропроводность (зависящая от температуры и плотности), фигурирующие в уравнениях электромагнитного поля, во внутренних итерациях в части II считаются неизменными.
Заметим, что уравнения электромагнитного поля линейны, поэтому при их численном решении необходимость в итерационном процессе отсутствует.
Точность расчетов контролировалась путем вычисления на каждом временном шаге потока магнитного поля в монокристалле (как в разрушенной, так и в неразрушенной областях). При принятых граничных условиях на поршне величина потока должна сохраняться. Отличие рассчитанного потока от начального не превышало 5% и не имело тенденции к нарастанию.
5.5.4. Анализ результатов расчетов
5.5,4.1. Безразмерные параметры сжатия магнитного поля. Качественное описание их влияния на динамику процесса
Сформулированная задача определяется следующими размерными параметрами: с — скорость света;
R0 — начальный радиус монокристалла;
B0 — начальное магнитное поле;
P0 — плотность кристалла (р0 = 4,5 г/см3);
134 5. Излучающая ударная волна: схождение до последнего микрона
g0 — максимальная электропроводность;
и0 — константа, характеризующая внутреннюю энергию
кристалла (и0 = 7,6 • 106 м2/с2);
ир — массовая скорость за ударной волной, характеризующая интенсивность ударной волны (по постановке задачи ир — скорость толкающего поршня). Вместо ир можно ввести р — плотность за волной, и характеризовать ею интенсивность сжатия
Р/Ро-
В предположении, что ударная волна сильная, напряжением перед волной Gy пренебрегаем по сравнению с давлением за волной. Согласно тс-теореме [5.10] существуют три безразмерных параметра, определяющих исследуемый процесс. В качестве их возьмем:
Здесь параметр D характеризует отношение приложенной кинетической энергии к внутренней энергии кристалла. Величина Rem — магнитное число Рейнольдса, определяет соотношение между процессами диффузии магнитного поля по частицам и его конвективным переносом за счет движения (Rew >> 1 соответствует вмороженности магнитного поля в среду). Величина TV-параметр взаимодействия, характеризующий влияние магнитного поля на движение среды, и представляет отношение членов, связанных с действием магнитной силы, к инерционным (конвективным) членам. Вместо N можно рассмотреть число Альфвена А2 = N/Rem = B0 /(47Ip0Wo)5 которое представляет отношение энергии магнитного поля к внутренней энергии кристалла.
Указанные параметры характеризуют процесс в целом, а не отдельные его стадии, поскольку составлены по начальным значениям входящих в них величин. Процесс сжатия характеризуется конкуренцией трех процессов. • Кумуляции кинетической энергии ударно-сжатого вещества. Давление, а значит, и массовая скорость вещества за фронтом волны растут. Так, например, в газовой динамике на заключительной стадии схлопывания (задача Гудерлея [5.3])
Предыдущая << 1 .. 31 32 33 34 35 36 < 37 > 38 39 40 41 42 43 .. 61 >> Следующая
Реклама
 
 
Авторские права © 2010 PiroChem. Все права защищены.