328
Глава 17. Импульсные генераторы
с
электрическим взрывом проводника
§ 17.3 Магнитогидродинамический метод расчета схем с ЭВП
При магнитогидродинамических (МГД) расчетах в каждый момент времени оп-
ределяются, в зависимости от координат, плотность тока, напряженность электриче-
ского и магнитного полей, плотность, температура, давление и массовая скорость
вещества, а также поле излучения в веществе (если оно существенно в данной зада-
че). Одновременно расчет дает зависимость от времени токов и напряжений в каж-
дом элементе электрической схемы. Как правило, основу таких расчетов составляет
численное решение одномерных МГД-уравнений в однотемпературном приближе-
нии с учетом электронной теплопроводности. Перенос излучения в веществе описы-
вается в приближении спектральной диффузии. Система уравнений в лагранжевых
координатах для случая цилиндрической симметрии приведена
в
[12].
При подобного типа расчетах определяющую роль играет корректность описа-
ния свойств вещества, уравнений состояния, зависимостей коэффициентов элек-
тро- и теплопроводности от состояния вещества, а также спектральных коэффици-
ентов поглощения излучения. Описание свойств вещества в широком диапазоне
состояний, включающем область перехода от металла к плазме, представляет со-
бой самостоятельную и достаточно сложную физическую проблему. При расчетах
ЭВП используются различные методы такого описания. Приведем здесь один из
методов. В [10] для высокотемпературной области (Т> 1 эВ) использовались урав-
нения состояния, рассчитанные по модели Томаса-Ферми [13]. Для задания коэф-
фициентов электро- и теплопроводности применялись интерполяционные форму-
лы, построенные на основании расчетов по полуклассической теории переноса.
Спектральное поле излучения, по которому проводится осреднение коэффициен-
тов поглощения, определялось решением уравнений спектральной диффузии из-
лучения [12]. Коэффициенты поглощения при этом задавались либо таблицами,
полученными на основании детальных квантовомеханических расчетов, либо по
аналитическим формулам [14]. Процедура осреднения проводилась через опреде-
ленное количество шагов по времени.
Рассчитанные, исходя из термодинамических соотношений, значения гидроди-
намических величин на границе области смеси фаз в зависимости от относитель-
ной плотности вещества 8 (§ 5.3) представлены на рис. 4. Критическая точка ха-
рактеризуется значениями 8 = 0,322, е
т
= 5,28 кДж/г. Внутри области смеси фаз
гидродинамические величины получались интерполяцией по концентрации одной
из фаз между значениями этих величин на границе. Для определения зависимости
электропроводности от состояния вещества был применен расчетно-эксперимен-
тальный метод, при котором электропроводность выбирали так, чтобы расчеты с
достаточной точностью описывали ряд «опорных» экспериментов, проведенных
при значительно отличающихся друг от друга условиях взрыва проводников. На
рис. 5.5 главы 5 приведены полученные таким способом зависимости удельного
сопротивления меди от 8 для нескольких значений удельной тепловой энергии.
Как видно из сравнения рис. 4 и 5.5, для фиксированных значений е
т
, меньших
критического значения, удельное сопротивление при уменьшении 8 резко возрас-
тает вблизи значений 8, соответствующих границе области смеси фаз, и имеет
максимум при плотностях, меньших критической плотности. Это обстоятельство