Отыскал ту самую модель AlexAV о которой говорил сегодня днем.
Вот наша переписка тогда:
Aex_Semenov:Они обещают 85% энергии взрыва вложить в кинетическую энергию плазмы из вольфрама. Но насколько я понимаю, их не интересуют БОЛЬШИЕ скорости истечения.
~100 Км/с. В "Орионе" это всех устриавает. Но не нас.
То есть они могут иметь маленькой мощности бомбу в очень тяжелом футляре.
Улавливает?
AlexAV: Интересный вопрос. Набросал здесь простенькую модель разлёта чёрнотельного сгустка (при этом предполагал профиль плотности гауссовым, газ идеально теплопроводным (в каждый момент времени температура всех точек в среде одинакова), эффективный радиус излучающей сферы 2s, где s-дисперсия профиля плотности). В такой модели конечная скорость разлёта и доля энергии, уходящая в излёчение, будет зависеть только от начальной плотности и удельной плотности энергии мишени.
Вот что получилось для облака урана с нормальной начальной плотностью.
По оси абсцисс начальная плотность энергии в килотоннах т.э. на кг. По оси ординат на первом графике средняя конечная скорость разлёта в км/c. На втором – доля энергии, которая уйдёт в излучение.
Какие можно сделать выводы.
1) Получить скорости разлёта большие 160 км/c невозможно для такой системы ни при каком энерговкладе.
2) Оптимальный энерговклад находится где-то в области2.5 тонн т.э. на кг, при этом скорость разлёта составит около 130 км/c, а потери на излучении – 20%.
3) Поскольку любое термоядерное устройство по схеме Теллера-Улама имеет непрозрачную оболочку, то, по всей видимости, для системы типа “Орион” величина удельного импульса более 100-200 км/c принципиально недостижима.
4) Практически вне зависимости от материала греть сгусток выше температуры 1-4 кэВ бессмысленно, т.к. практически вся энергия сверх этой уйдёт в излёчение.
5) Поскольку разлёт чёрнотельного облака в первом приближении не слишком сильно зависит от материала, то материалы для которых критерий Лоуренца в системах прозрачных для излучения не может быть достигнут (например B-11 + p или Li-7 + p) для создания высокоимпульсного ТЯРД с инерционным удержанием непригодны. Т.е. выбор топлива для таких систем ограничен T+D, He-3 + D, D+D, и, быть может трёхкомпонентной системой Li-6+D+T. Если отбросить варианты со скоропортящимся тритием, то вообще остается только два.
Aex_Semenov: Интересно...
А как это согласуется с реальными экспериментами? То есть ядерными взрывами в космосе
AlexAV: Кажется, разобрался. Температура излучающей поверхности будет ниже, чем средняя температура среды из-за ограниченной теплопроводности. Очевидно, что для достаточно мощного устройства отношение температур внутри объёма и поверхности T/T0~1/R~1/E^(1/3), E – мощность устройства.
Но интереснее другое. Даже если мы устремим мощность устройства к бесконечности, а значит сделаем излучение через поверхность относительно малым, это всё равно нам поможет весьма ограниченно. Дело в том, что среда будет находится в равновесии с фотонным газом, причём энергия последнего существенно превышает энергию массивных частиц. При расширении такой системы охлаждении происходит не столько за счёт ускорения вещества, сколько в результате расширения фотонного газа. В итоге перекачка энергии среды в кинетическую энергию расширяющегося газа происходит медленно и на начальных стадиях расширения перекачивается менее 20% (при плотности энергии более 1 килотонны т.э. /кг). На поздних стадиях ускорение массы вещества становится совсем неэффективным, завися от приращения объёма системы по закону медленнее логарифмического, и даже незначительного излучения через внешнею поверхность оказывает достаточно, чтобы полностью остановить его.
В итоге получаем следующую зависимость конечной скорости вещества и доли энергии перешедшей в излучение от энерговклада (обозначения такие же, как и на предыдущем графике):
Какие можно сделать выводы.
1) Как видно, даже если предположить, что температура излучающей поверхности много ниже, чем средняя, то не смотря на это для зарядов с высокой плотностью энергии (более 1 кт.т.э) потери энергии через излучения оказываются очень значительны и превышают величину 80% всей энергии взрыва.
2) Величина конечной кинетической энергии вещества оказываются намного больше, чем в предыдущем случае, и оно может достигать скорости до 1600 км/c (что согласуется с приведенными вами данными).
3) По всей видимости, для схем типа Орион поднимать энергетику устройства выше 0.1 кт т.э./кг не имеет особого смысла. При этом достигается удельный импульс до около 850 км/c, что вполне соизмеримо с величиной 1600 км/c для устройства с энергетикой 6 кт т.э./кг, зато мощность рентгеновского излучения окажется почти в 500 раз ниже, чем в последнем случае.
4) По мере снижения мощности устройства (при заданной удельной энергетике), по всей видимости, относительные потери на излучение будут увеличиваться, что делает оправданным использование устройств мегатонного класса.
Картина уже не такая печальная, как в предыдущем случае, однако всё равно максимальный УИ 800 - 1500 км/c (реально меньше, ведь импульс разлетающейся плазмы ещё надо преобразовать в импульс направленной струи) для МП кажется несколько маловатым, позволяющим в лучшим случае развивать скорости 0.005с – 0.01с.
Да и даже 14% (для удельной энергии 0.1 кт т.э./кг), не говоря уже о 80%-95% (для 1-6 кт т.э./кг) мощности взрыва мегатонного класса в виде рентгеновского излучения – серьёзная проблема.
Aex_Semenov: Алекс, пожалуйста. Расжуйте мне вашу модель (последние кривули) подробней. Это очень важно!
AlexAV: Хорошо. Берём гауссов профиль начальной плотности. При этом предположим, что излучение сосредоточено в центральной области радиусом 2s, где s – дисперсия распределения. Также предположим, что градиенты температур в центральной области малы и сосредоточены в разреженной “короне”. Кроме того, положим, что температура области, с которой фотоны излучаются во внешнее пространство низка, и этим излучением можно пренебречь. Плазма рассматривалась как идеальный газ. В этом случае при расширении сгустка форма профиля его плотности будет сохраняться, в профиль скорости будет иметь вид V(r,t) = r F(t), где F(t) – некоторая функция, зависящая только от времени. Это позволяет свести задачу Коши для уравнения в частных производных к задаче Коши для системы из 3-х дифференциальных уравнений:
Уравнение для дисперсии:
ds/dt = s F
Уравнение для функции скоростей F:
dF/dt = 2RT/(mu s
2) – F
2 (R – универсальная газовая постоянная, mu – средняя молекулярная масса частиц плазмы с учётом электронов)
Уравнение для температуры:
dT/dt = - (Wad(F,T) + Wr(s,F,T))/C(s,T)
где
С(s,T) = 3/2 (M/ mu)R + (512/3)pi (st/c) s
3 T
3 (M – масса сгустка, st – постоянная Стефана – Больцмана, с – скорость света) – теплоёмкость системы с учётом вклада фотонного газа.
Wad(F,T) = 3 (M/mu) R F T – функция охлаждения в результата адиабатического расширения плазма.
Wr(s,F,T) = 128 pi (st/c) s
3 F T
4 – функция охлаждения в следствии расширения фотонного газа.
Начальные условия
s0 = (M/(pi^(3/2) ro))
1/3 (ro – начальная плотность)
F0 = 0
U(s0, T0) = M H (H – начальная плотность энергии)
U(s0, T0) = 3/2 (M/mu)RT + (128/3)pi(st/c) s
3 T
4 – энергия сгустка.
Модель конечно приближённая, по-хорошему надо правильно учитывать теплопроводность и профиль температур, да и уравнение состояние лучше брать как минимум в форме Томаса-Ферми, а не идеальный газ. Но не думаю, что учёт этих факторов изменит картину качественно.
P.S. Попробую сделать одномерную гидродинамику разлёта чёрнотельного облака с учётом вышеперечисленного, но в ближайшие дни не обещаю.