Интерес к задачам о цилиндрически и сферически симметричном нестационарном сжатии стимулирует ряд приложений, включая проекты реализации инерциального управляемого термоядерного синтеза [1—6].
В связи с ограниченными возможностями традиционной энергетики и постепенным исчерпанием разведанных запасов топлива, на повестку дня выходит управляемый термоядерный синтез (УТС). Согласно современным представлениям звезды светят благодаря стационарной термоядерной реакции, и неконтролируемый термоядерный процесс удалось реализовать с помощью водородной бомбы.
В лабораторных условиях УТС реализуется, но, к сожалению, пока цена ядерной энергии слишком высока. Иными словами, затраты энергии превышают энергетический выигрыш.
Проблема термоядерного синтеза состоит в решении двух задач: нагрева вещества до необходимых температур и его удержания на время, достаточное для «сжигания» заметной части термоядерного топлива. Это время определяется критерием Лоусона [7]: лт > const где г — время’удержания высокотемпературной плазмы в системе, п — концентрация её частиц, константа зависит от реакции. В таком виде критерий удобен для перспективных реакторов непрерывного действия.
УТС с инерционным удержанием, или инерциальный УТС, основан на возможности получения положительного энергетического эффекта в реакторах импульсного типа. В инерциальном УТС небольшая масса, миллиграммы, термоядерного топлива (дейтериево-тритиевой смеси) сжимаются оболочкой, ускоряемой за счет реактивных сил, возникающих при испарении оболочки с помощью мощного облучения. Облучение термоядерной мишени осуществляется лучами лазера, пучками ионов или электронов [1−4, 8, 9]. Энергия выделяется в виде микровзрыва, когда в процессе сжатия в плазме мишени достигаются необходимые условия для термоядерного горения. Время жизни такой плазмы определяется инерционным разлетом смеси и поэтому критерий Лоусона для> инерционного удержания принято записывать в терминах произведения рг, где р — плотность реагирующей смеси и г — радиус сжатой мишени. Для того, чтобы за время разJ лета смесь успела выгореть, нужно, чтобы рг >/*,/* ~ 0.1 — 3 г/см [7]. Для о субмиллиметровых мишеней это отвечает плотностям в сотни г/см, т. е. необходимо достичь, плотности. термоядерного вещества, в тысячи и более раз превышающую плотность твёрдого дейтерия.
Из критерия Лоусона следует, что чем меньше масса, сжигаемой мишени, тем сильнее надо ее сжать. Поэтому реализация неуправляемого термоядерного синтеза требует во много раз меньшей степени сжатия вещества по сравнению с управляемым.
Попытки’реализовать, термоядерное горение с помощью взрывчатых веществ, или газодинамический термоядерный синтез, на данный момент не увенчались успехом-[10].
К настоящему времени энергоэффективную реакцию удалось получить в водородной бомбе. Она состоит из взрывателя — атомной* бомбы и термоядерного топлива, дейтерида лития-6 LiD, которое советские оружейники назвали «Лидочкой» [11].
Взрыв термоядерной бомбы доказал возможность синтеза с инерционным удержанием.
Недавно в Ливерморе (США) заработала мощная лазерная установка National Ignition Facility с энергией импульса до 1.8 МДж [12]. Лазерные лучи со всех сторон направляют на мишень, содержащую термоядерное топливо. Под действием светового давления и реактивной силы от испаряющегося с поверхности вещества происходит одновременно сильное сжатие и нагрев мишени. Использование лазеры с энергией порядка 1 МДж в установках Laser Megajoule, Бордо, Франция и National Ignition Facility, Ли-вермор, США, должно обеспечить условие термоядерного зажигания и энергоэффективные режимы УТС [7].
Взрыв полноценной водородной бомбы — неуправляемый термоядерный синтез, при котором выделяется слишком много энергии, что делает его непригодным для энергетических целей.
Предложен метод получения ядерной энергии путём взрывов атомных зарядов, инициирующих DD-реакцию. Согласно ему предлагается производить в камере котла вспышечного сгорания термоядерные взрывы большой мощности (а не микровзрывы, как винерциальном термояде) с целью получения энергии. [13].
Проблема реализации реактора типа котла вспышечного сгорания* состоит в большой мощности взрыва, а значит и большой опасности такого реактора, а также том, что все атомные взрывы, даже подземные, сейчас запрещены.
Сильное сжатие предполагает неограниченный^ рост давления, а при очень больших давлениях для вещества становится применимо уравнение состояния’идеального газа. Например, горные породы* при< прохождении через них сильной ударной, волны ведут себя какидеальный газ с показателем адиабаты у = 3, так называемый газ Ландау-Станюковича. Для замагниченной плазмы используется уравнение состояния' идеального газа с показателем адиабаты 2.
Вследствие пространственного усиления^ из-за уменьшения сечений трубок тока некоторые параметры по мере приближениям центру симметрии неограниченно возрастают, это приводит к тому, что начальные значения этих параметров перестают влиять на течение, т. е. «забываются», а само-те-чение выходит на автомодельный режим. Так происходит в задачах о схождении ударной волны к центру симметрии и схлопывании пустой полости, в которых из всех параметров вначале неавтомодельного течения не «забываются» только начальная плотность перед волной и начальная энтропия соответственно.
Негазообразные вещества при меньших давлениях можно приближенно заменить идеальным газом с большим эффективным показателем адиабаты, вплоть до у = оо.
В диссертации рассмотрены 4 задачи о сжатии цилиндрических и сферических объемов газа:
1) Автомодельная задача о сходящейся к центру симметрии сильной ударной волне для показателей адиабаты от 1.001 до 3 (задача Гу-дерлея).
2) Неавтомодельная задача о быстром сильном сжатии,.предполагающем одновременное получение сколь угодно больших температур и плотностей для показателей адиабаты. от 1.001 до 3. К особой характеристике задачи Гудерлея примыкает пучок волн сжатия, сфокусированный наконечном расстоянии от центра симметрии. Рост температуры обеспечивается интенсивностью сходящейся ударной волны, а рост плотности происходит за ней в центрированной волне сжатия.
3) Автомодельная задача о сходящейся к центру симметрии сильной ударной волне с изменением показателя адиабаты на фронте отраженной ударной волны.
4) Автомодельная задача о схлопывающейся к центру симметрии пустой сферической полости с изменением показателя адиабаты, на фронте ударной волны.
Автомодельная задача о схлопывающейся полости, впервые рассмотренная в [14] в сферически симметричном случае, позже рассматривалась как отдельно [15, 16, 17], так и совместно [18, 19] с одной из основных газодинамических задач атомного проекта — задачей о схождении ударной волны к центру симметрии? и отражении от него, впервые рассмотренной Г. Гудерлеем [20].
В литературе отсутствовало полное решение задачи Гудерлея для показателей адиабаты от 1 до 3.
С помощью сжатия сходящейся ударной волной^ или ударного сжатиятеоретически можно достичь сколь угодно высоких температур, определяемых интенсивностью ударной: волны. При этом критерий Лоусона для микромишеней термоядерного/ синтеза не будет удовлетворен, требуется больший рост плотности: Изэнтропическое или. безударное сжатие может дать сколь угодно высокую плотность [5, 21−24]. При. выполнении критерияЛоусона уровень температуры будет недостаточным для реализации У ТС [25].
Основные классические задачи о сжатии рассмотрены в обзоре [24].
По вышеуказанным причинам для достижения условий, необходимых для реализации энергоэффективного УТС, нужно объединить ударное сжатие с безударным. Для?минимизации, диссипативных эффектов-.сжатие должно быть как можно более быстрым [25].
Задача о быстром сильном сжатии была сформулирована А.Н. Крайко• «и решена для нескольких показателей адиабаты, 6/5, 7/5 и 5/3 [26]. Быстрое-сильное сжатие получается присоединением центрированной волны. сжатия к особой характеристике задачи Гудерлея, приходящей в центр симметрии одновременно со сходящейсяударнош волной: Позже оказалось, что «построенное решение го дитсяшоказателей’адиабаты, не превышающих Г. 9.
Третьяш четвертая задачшранее в литературе не рассматривались.
Обзор литературы по теме диссертации.
Автомодельная задача о схлопывающейся сферически симметричной полости в газе исторически не представляла большого самостоятельного интереса, поэтому ей посвящено относительно небольшое количество работ.
Будучи очень близкой по физическим и математическим особенностям, она рассматривалась совместно с задачей о схождении ударной волны к центру симметрии, рассмотренной Г. Гудерлеем [20] в 1942 г. Немецкий физик в то время занимался кумулятивными снарядами и побочно пришел к идее взрыва во внутрь (имплозии). В дальнейшем идея имплозии была подхвачена физиками-ядерщиками, участвовавшими в атомном проекте [27].
В работе [14] рассмотрена задача о схлопывающейся сферической полости в жидкости с эффективным показателем адиабаты у = 7. Рассмотрена движущаяся от центра симметрии ударная волна, считающаяся" гомоэнтро-пичной, — энтропия газа в ударной волне остается* неизменной. Возможно, такой подход упрощает решение уравнений газовой динамики, но не отражает сути явления. Поэтому от такого подхода следует отказаться.
В работе [28] рассмотрены возможные типы автомодельных течений, описывающих симметричное* течение жидкости в пустую сферическую полость, при показателях адиабаты Г < у < 7. Установлено, что автомодельные течения, при которых жидкость ускоряется в полость, существуют при показателях адиабаты у > 3/2. При* меньших значениях у существует решение с постоянной скоростью границы.
Подробное исследование особых точек и вопроса о количестве возможных автомодельных решений приведено в работе [18]. При этом отмечается, что, в соответствии с гипотезой Гельфанда [29], в случае, когда через особую точку на звуковой линии — узел — можно теоретически пройти неединственным способом по аналитической кривой без разрыва производных любого порядка во всех точках, практически реализуется единственное, особое решение с минимально* возможным показателем автомодельности. Интегральная кривая при этом является отдельным усом вышеупомянутой узловой особой точки. В работе [30] показано, что при выборе особых начальных условий может реализоваться альтернативное, неособое решение, при котором показатель автомодельности совпадает с найденным в [18] и интегральная кривая является усом общего направления особой точки на звуковой линии.
В работе Лазаруса [19] для множества показателей адиабаты от нижней границы существования особого решения до 100 приведены показатели ав-томодельности с точностью до 10 знака после запятой, проведено исследование особых точек, построены графики основных параметров, а также построены примеры решений с разрывом производных на особой характеристике. Позже выяснилось, что такие решения нестабильны относительно симметричных (одномерных) возмущений [31]. Основные результаты работ [19, 31] подтверждаются численным моделированием [32] перехода схлопывания полостей в автомодельный режим. При этом для значений У < У er, Jcr ~ 3/2 действительно реализуется решение с практически постоянной скоростью границы полости, с точностью до 0.3%. А при больших значениях у, вплоть до 2, область существования автомодельного решения с непостоянной скоростью границы полости составляет порядка 1% от начального радиуса полости.
Исследование [33] показало, что учет неавтомодельных эффектов на последней стадии схлопывания пустой полости в воде, хотя и меняет распределение параметров, практически не сказывается на траектории границы полости. Показано, что схлопывание полости сопровождается очень высокими давлениями при относительно высокой скорости коллапса.
В работе [34] построено глобальное аналитическое решение в области между поверхностью полости и звуковым фронтом в виде бесконечных рядов. Решение существует вплоть до момента схлопывания полости при у < 5/3 в сферическом, и у < 2 в цилиндрическом случаях, при больших значениях у ряды расходятся, начиная с определенного момента времени. Результаты работы находятся в соответствии с выводами [33]: при значении у «2 граница полости вначале движется с постоянной скоростью, а затем выходит на автомодельный режим [19].
Обзор некоторых последствий, связанных со схлопыванием пустых полостей, а также возможности применения этого явления для зажигания взрывчатых веществ, можно найти в работе [35].
Как уже отмечалось, ключевая газодинамическаяпроблема в термоядерном синтезе — задача о сходящейся к центру или оси симметрии ударной волне большой интенсивностиЕудерлей опубликовал решение в сферическом случае дляпоказателя адиабаты, у — 7/5: [20]. Плотность перед сходящейся волной? в классической постановке1 считается^ постоянной. Согласно этому решению интенсивность.волны.по мере движения к. центру симметрии увеличивается, а после отражения ударной волны в центре симметрии температура бесконечна! В" СССР" этой, задачей занимались. Лі Д. Ландау, К.ГГ. Станюкович [36] и ряд других ученых.
Исследователи^ занимающиеся задачей Гудерлея, сталкиваются?- с: множествомсложностей: Для полного решения необходимосоединитькачественный1. анализ: обыкновенных дифференциальных уравнений^ с численным исследованием их поведения: Ключевойпараметр автомодельной задачи. Гу-дерлея— показатель, автомодельности — является трансцендентным числом и. определяетсяшщроцессе решенияї Значения показателя-автомодельности для-различных показателей адиабаты. приведены в работе [19] с точностью до 10 знака после запятой. Віработе [37], рассмотрена^ задача со степеннымшрофи-лем плотности. В этой постановке: задачи дляосновных показателей. адиабаты в диапазоне от: 1.1 до оо приведены-показателиавтомодельности с точностьюдо- 15-го знака после запятой, а в задаче с постоянной плотностьюпри у = 1.4 для сферического и цилиндрического случая аж до 30-го знака1 после запятойПомимо этогопостроены графики распределеншгг параметров> за сходящейся ударной волной до момента ее прихода в центр симметрии.- Дополнительные сведения’о вкладе начального степенного профиля в решение приведены в работах [38, 39].
Приближенный метод вычисления показателя автомодельности (с точностью до нескольких знаков после запятой), основанный на мнимом устранении сингулярности на звуковой линии с помощью алгебраических преобразований, приведен в [40].
Большая часть научных публикаций по задаче Гудерлея затрагивает решение задачи до момента фокусировки ударной волны в центре симметрии, оставляя вопрос о проходящих после этого процессах открытым. В некоторых работах приводятся качественные сведения о том, что происходит после фокусировки ударной волны в центре симметрии [36]. Работа Лазаруса [19] содержит довольно подробное исследование задачи Гудерлея, включая отраженную ударную волну. В этой работе, помимо множества других графиков, для значений показателя адиабаты у = 7/5, 5/3, 2, 3, 6 приведены гра.
Рисунок 0.1. Сравнение показателей автомодельности, вычисленных Лазару-сом [19], сферический случай отвечает сплошной линии, цилиндрическийпунктиру, и Фуджимото и Мишкином [46, 47], показаны крестиками в сферическом и плюсиками в цилиндрическом случае. Рисунок из работы [45] фики основных параметров газа за отраженной ударной волной. Графики построены в абсолютных масштабах, что несколько затрудняет их сравнительное исследование. Также приведены неаналитические решения (с разрывом производных на особой характеристике), вопрос о стабильности оставлен для будущих исследований. Выполненное позднее исследование стабильности автомодельных сходящихся течений [31]. показало нестабильность неаналитических решений и в. случае сходящейся ударной волны.
Математически строгий способ решения задачи Гудерлея предполагает исследование особых точек в плоскости безразмерных скорости газа — скорости звука (или квадрата скорости звука) [18, 19] и должен обеспечивать прохождение интегральной кривой’через определенную особую точку. Также имеются приближенные методы решения задачи Гудерлея-, самым известным из которых, является метод Чиснелла-Честера-Уизема, известный также как. метод СС" ^ по первым английским буквам, фамилий авторов, разработавших его [41−43]. К сожалению, известные приближенные методы не дают инфорг мации о течении после отражения ударной волны от центра симметрии. Помимо этого, имеются работы с т.н. аналитическим решением' задачи Гудерлея, где за1 «довольноточное» или даже точное значение показателя' автомодельности берется-его первое приближение, используемое при итерационном вычислении показателя автомодельности [44]. При некоторых показателях адиабаты оно близко к истинному, при других далеко (см. рисунок 0.1 из работы [45]). И приводятся «доказательства» применимости такого подхода, основанные не на глубоком анализе проблемы перехода, через особую точку на звуковой линии, а на основе гипотезы, о максимуме давления, [46−49]. Работа [46] рассматривает цилиндрически, а. 47] - сферически симметричный случай. Лазарус отверг такой подход [50], но Мишкин не согласился с приведенными доводами [51]. Более подробная' критика взглядов Мишкина и Фуджимото содержится в работах [44, 45], поддержка — в,[52].
Внекотором диапазоне показателей адиабаты задача Гудерлея, как и задача о схлопывающейся сферической полости, допускает неединственное решение без разрыва производных [18], но согласно гипотезе Гельфанда [29] реализуется решение с наименьшим показателем автомодельности. Это подтверждается результатами исследования [53], в котором решение построено в виде многочлена 40-го порядка от времени.
Вопросы устойчивости схождения ударной волны рассматривались разными авторами [31, 54]. Сферическая сходящаяся ударная волна оказалась неустойчивой относительно радиальных возмущений [54]. А решения с разрывом производных на звуковой характеристике, примеры которых приведены в [19], оказались неустойчивыми относительно колебаний, направленных вдоль радиуса, т. е. безусловно неустойчивыми, а потому нереализуемыми в общем, случае [31]. Области устойчивости относительно различных колебаний найдены, в численном исследовании [55]'. В работе [56] исследована стабильность сходящейся ударной волны как линеаризованных уравнений метода Чиснелла, так и в нелинейном случае. Некоторые моды колебаний, стабильные для линеаризованных уравнений, оказались нестабильными в нелинейном. случае. Рост возмущений, оказался не экспоненциальным, а степенным. Рост амплитуды" не зависит от моды. При достаточно большой амплитуде возмущения форма волны испытывает излом, и это ведет к потенциальным потерям энергии сходящейся ударной волны.
В обзорной работе [24], в частности, приведены полные качественные решения как задачи Гудерлея со степенным профилем плотности перед ударной волной, так и задачи о схлопывающейся полости. Скорость звука на границе полости ошибочно считается бесконечной, что приводит к отличному от других авторов решению. Приведены некоторые количественные данные.
Анализ неавтомодельного схождения ударной волны методом интегральных соотношений [57], методами характеристик, случайной выборки, и конечноразностным второго порядка с искусственной вязкостью [58] показал, что течение приближается к автомодельному, но никогда не выходит на истинно автомодельный режим. Скорость приближения оказалась ниже, чем предполагалось. Приближение к автомодельному решению, особенно профиля плотности, быстрее при большей начальной энергии, вызвавшей сходящуюся ударную волну (рисунок 0.2 из [57]). «Уо — энергия инициирования Ш.
Рисунок 0.2. Распределение безразмерной скорости Г частиц около оси симметрии [57] сходящейся ударной волны,? — безразмерное расстояние, пройденное фронтом ударной волны. Начальной точке соответствует значение Е, = 0, моменту прихода волны на ось симметрии отвечает Е, = 1, г — расстояние до оси симметрии, Кз — радиус ударной волны.
Решение неавтомодельной задачи методом характеристик [59] показывает постепенный выход течения на автомодельный режим.
В одной из первых отечественных публикаций по задаче Гудерлея, в книге [60], задача рассматривается качественно, приводятся приближенные количественные оценки параметров задачи, в частности значения показателя автомодельности для нескольких показателей адиабаты.
В работе [61] течение со сходящейся к центру симметрии ударной волной рассчитано методом конечных разностей Лакса, для конечности давления на оси симметрии введена искусственная диффузия. Найдено, что рассчитанные параметры газа находятся в хорошем соответствии с расчетами по методу Чиснелла.
В работе [62] рассмотрена фокусировка ударной волны с учетом лучистой теплопроводности, определена ширина зоны прогрева. Показано, что при лучистом теплообмене неограниченная кумуляция не исчезает, а видоизменяется. Вместо конечной плотности и бесконечной температуры получается конечная температурами бесконечная плотность. Численное исследование процесса превращения автомодельного решения Гудерлея в изотермическую волну под действием электронной теплопроводности [63] показало, что теплопроводность очень сильно влияет на максимальные значения параметров в центральной области, но качественная картина течения остается неизменной. Сходящаяся ударная волна, если и вырождается в движущуюся с постоянной скоростью изотермическую волну вследствие влияния теплопроводности, то существует в настолько незначительной области, что попытки ее обнаружить не увенчались успехом. Этот результат дает подтверждение выводу работы [62] о неограниченной кумуляции возмущений, идущих к центру в теплопроводном газе.
При особой геометрии экспериментальной установки можно получить сверсферическую кумуляцию, при которой интенсивность сходящейся ударной волны нарастает быстрее, чем в сферическом случае [64].
Работа [65], в частности, рассматривает вклад в автомодельное решение задачи Гудерлея возмущений, связанных с противодавлением и однородным тепловыделением.
В работе [66] рассмотрено схождение ударных волн конечной интенсивности с помощью приближенного квазиавтомодельного метода Ошимы. Показано, что течение асимптотически приближается к автомодельному решению задачи Гудерлея, и в пределе совпадает с ним. Обнаружено хорошее соответствие результатов метода Ошимы с результатами других приближенных методов.
В работе [67] рассмотрена начальная стадия схлопывания цилиндрической или сферической волны. Получено аналитическое решение с точностью до членов четвертого порядка по пространству. Классическое автомодельное решение является нулевым приближением полученного решения. Исследование дает отличный способ выбора подходящих начальных условий для численного расчета. Работа [68] дополняет [67] в области, где результаты последней работы неприменимы, совпадая с ее результатами возле цилиндрической стенки. Исследование [69] показало, что учет двух последующих членов к решению Гудерлея-существенно расширяет область справедливости разложения-по сильным ударным волнам.
Схождение сильных ударных волн по среде, состоящей из газа и твердых несжимаемых частиц, рассмотрено * без учета, теплопроводности [70−72], при конечной теплопроводности [73] и при очень большой, так что течение1 заударной волной можно считать изотермичным [74, 75].
В работе [76] дано автомодельное представление для уравнений радиационной газовой динамики со сферической симметрией, подходящее для широкого разнообразия законов излучения. Построены численные решения для сходящейся ударной волны с радиационными потерями в оптически толстом пределе. Более подробно схождение ударной волны в излучающем газе рассмотрено в [77].
В работе [78] рассмотрено схождение ударной волны с гомотермичным течением за' ее фронтом. Показатель автомодельности, как и в задаче Гудерлея, определяется в процессе решения. Проведено сравнение пространственных распределений параметров в этой задаче и в задаче Гудерлея. Полная энергия в автомодельной области убывает более быстро, чем в изэнтропиче-ской постановке. При этом задача имеет единственное решение при у < 3 и неединственное при больших показателях адиабаты у [79].
В работе [80] проведен анализ оптимизированных имплозий сферических дейтерий-тритиевых мишеней, инициируемых импульсами лазера и ведущих к эффективному термоядерному горению. Получены параметры оптимального импульса для мишеней массой от 40 нг до 250 мкг.
В работе [81] исследовано влияние плотной оболочки термоядерных мишеней на имплозию. Рассмотрены мишени как со"сферической полостью в центре, так и без нее.
В работе [82] модифицированным методом Уизема теоретически исследован процесс генерации и. схождения, ударных волн произвольной формы.
В работе [83] экспериментально получена сходящаяся детонация в ацетилен-кислородной смеси. В пределе сильно пересжатая детонация приближается к ударной волне. Измеренные экспериментально значения давления детонации находятся в хорошем* соответствии с теоретически предсказанными^ по методу Чиснелла. Опыты со сходящимися ударными волнами, приведенные в [84] так же показали соответствие с теоретическими предсказаниями, в частности, радиусы сходящейся и отраженной ударной волны, измеренные по шлирен-фотографиям, являлись функциями времени, близкими1 к степенным.
В работе [85] обсуждается возможность получения условий’термоядерного горенияс помощью лазеров. Оценен требуемый энерговклад. Для начальной температуры за сходящейся ударной волной^ получена оценка в 7×106 К.
В работе [86] экспериментально получена цилиндрическая сходящаяся детонация в ацетилен-кислородной смеси, а так же описан метод получения цилиндрической имплозии в бомбе, а в работе [87] исследования продолжены. Из-за отсутствия вихревой крупномасштабных вихрей за отраженной ударной волне имплозия признана устойчивой. После отражения ударной волны от центра симметрии спектроскопические наблюдения зафиксировали температуру в 1.89×105 К.
В работе [88] описано получение технических алмазов в имплозийной камере иТ1А8 из графита с помощью взрывчатых веществ.
В работе [89] приведены результаты спектроскопических измерений имплозии в водород-кислородной смеси. Зафиксирована температура в 4500—6000К и давление примерно 12 500 атм при начальном давлении в 7 атм. В опыте не наблюдалось никаких больших неустойчивостей или нарушений.
Рисунок 0.3. Процесс развития нестабильности 4-й моды. После отражения ударной волны в центре области остаются вихри [92] симметрии.
В работе [90] представлен экспериментальный способ получения сходящихся цилиндрических волн. Усиление волн в эксперименте оказалось в хорошем соответствии с предсказаниями метода Чиснелла. Эксперименты [91] показали, что сходящаяся цилиндрическая ударная волна с возмущением сохраняет симметрию вплоть до малого радиуса, при котором в кривизне ударной волны внезапно появляется излом. Искусственно вызванные возмущения ведут к такому типу нестабильности. Во всех случаях нестабильность.
Рисунок 0.5 Меридиональные сечения обжатых металлических шаров [96].
Рисунок 0.6. Образование пустой полости в точке фокусировки сходящейся ударной волны в стальном шаре [97] проявляется в наличии вихревых пар за отраженной волной.
Экспериментальное исследование стабильности сходящихся цилиндрических ударных волн [92] показало, что развивается нестабильность 4-й моды (см. рисунок 0.3).
В экспериментах [93] в детонационной камере с размером волн поряд.
7 8 ка 500 мм наблюдалась температура от 10 до 10 К, плотность ионов при этом была от 1019 до 1021 см'3, что предполагает возможность реализации термоядерного синтеза в камере в 3 раза большего размера.
В ряде экспериментов по ударно-волновому или взрывному нагруже-нию металлических шаров наблюдается интересное явление: увеличение размеров образца после нагружения. При этом в центре образца наблюдается сферическая полость [94−97] (см рисунки 0.5−0.7). Также ссылки на другие подобные эксперименты можно найти в обзоре [10].
При этом в случае, когда точка фокусировки не совпадает с центром шара, образуются две полости, одна в точке фокусировки, другая в центре шара (см. рисунок 0.7).
Рисунок 0.7. Меридиональные сечения обжатых шаров [95].
Актуальность.
Вопросы, рассмотренные в диссертации — возможность достижения условий реализации инерциального УТС с помощью одной ударной волны, а также с помощью следующей за сходящейся’ударной волной непрерывной волны изэнтропического сжатияпостроение траектории поршня, осуществляющего близкое к оптимальному сжатие с одновременным получением сколь угодно больших плотностей и температурвлияние изменения показателя адиабаты на отражённой ударной волне на картину течения — являются важными, для проблемы инерциального УТС. Все вышеуказанные вопросы рассмотрены для широкого диапазона показателей адиабаты. Эти вопросы являлись предметом внимания многих исследователей в последние годы. Можно считать, что. и в будущем они также окажутся важными, в проблеме реализации инерциального УТС.
Результаты диссертации применимы в численном, моделировании сильного сферически или цилиндрически симметричного сжатия газа, а также схлопывания пустой сферической полости в газе, как для сравнения полученных решений с автомодельной асимптотикой, так и для. разрешения особенности около центра, или оси симметрии (далее — центра симметрии), возникающей в момент прихода ударной волны или границы*полости в него.
Цель и предмет исследований.
В диссертации исследуются процессы автомодельного схождения ударной волны к центру симметрии в сферически и цилиндрически симметричном случаях как в постановке с постоянным показателем адиабаты во всём поле течения, так и в постановке, допускающей изменение показателя адиабаты на отражённой ударной волнесхлопывания пустой сферической полости в постановке, допускающей изменение показателя адиабаты на ударной волне, идущей от центра симметриипостроение дополнительной к сходящейся ударной волне неавтомодельной изэнтропической волны сжатия, сфокусированной в центре симметрии.
Тема диссертации была предложена А. Н. Крайко. Под его руководством были выполнены все этапы работы. Первоначально планировалось рассмотреть сжатие газа сходящимися ударными волнами и дополнительного сжатия в центрированной волне сжатия, сфокусированной в центре симметрии. Позже было-принято решение рассмотреть также задачу о схлопыва-нии пустой сферической полости в газе в постановке, допускающей изменение показателя адиабаты на ударной волне, идущей от центра1 симметрии, как близкую к задаче о схождении ударной волны к центру симметрии по физическим и математическим особенностям.
1. Вфаботах [98, 101] построено полное решение автомодельной задачи о схождении ударнот волны к центру симметрии, для показателей адиабаты от 1.2 до 3. Цель построения1 полного решения состоит в том, чтобы проверить возможность достижения условийнеобходимых для реализации инер-циального УТС, с помощью сжатия сходящейся ударной-волной большой интенсивностиа также в использовании' полученного решения для построения траектории поршня, обеспечивающего одновременное достижение сколь угодно больших температур и плотностей.
2. В-работах [99, 102, 103, 105, 107−111, 115−117] рассмотрено цилиндрически и сферически симметричное быстрое сильное сжатие идеального совершенного газа с показателями адиабаты от 1.001 до 3. При быстром сильном сжатии одновременно достигаются сколь угодно большие плотности и температуры за время, много меньшее времени пробега звуковой волны через несжатый объем газа. Цель такого рассмотрения состоит в том, чтобы показать возможность достижения чисто газодинамическими средствами условий, необходимых для реализации инерциального УТС.
3. В работах [100, 104, 106, 112—114] рассмотрены автомодельные-зада-чи о схождении ударной волны к центру симметрии в сферически и цилиндрически-симметричных случаях и о схлопывании пустой сферической полости с изменением показателя адиабаты на ударной волне, идущей от центра симметрии. Цель такого рассмотрения состоит в том, чтобы выяснить поведение вышеуказанных нестационарных течений совершенного газа при-изменении показателя адиабаты на ударной волне, идущей от центра симметрии. Уменьшение показателя адиабаты приближенно моделирует различные физико-химические процессы (например, ионизацию и диссоциацию), а в задаче о схлопывании полости — фазовый переход (превращение жидкости в пар).
Структура и объём работы.
Диссертация состоит из введения, трёх, глав, заключения и списка литературы. Общий объём диссертации 137 страниц, 41 рисунок.
Список литературы
включает 128 наименований.
Заключение
.
Рассмотрены четырезадачи о сжатии цилиндрических и сферических объемов газаДля двух, из них рассматривавшихся: ранее, решение в более широком диапазоне параметров, в частности, для показателей адиабаты 1.001 < у < 3. Другие две задачи рассмотрены впервые.
1) Получено полное решение задачи о схождениигударной волны-к центру или оси симметрии для показателей адиабаты. 1.001 < у < 3 с шагом, непревышающим 0:1. Расчеты показали, что1 основное изменение: плотности, внутренней энергии, а, следовательно, и температуры за отраженной ударной волной происходит в малой окрестности центра или оси симметрии. Большой рост плотности? (более: чем-, втысячи раз)? при? сжатии, сходящейся ударной волной5 реализуетсятолько для газовс показателем адиабаты, близкимк 1. Всюду за. отраженной? ударношволной! кинетическая" энергия-газа составляет малую? долю от его? внутренней? энергииИзменениепо5радиусу давления за> отраженной ударной волной не превышает 1.5 раза.
2) Рассмотрена задача о быстром сильном сжатии в цилиндрически и сферически симметричном! случаях. Для показателей' адиабаты 1.001 < у < 3 построенырешения со сферически, и цилиндрическисимметричными иеав-томодельными-пучками воли сжатиях фокусировкой вблизи центра симметрии на особой характеристике и на отраженной ударной волне задачи Гудер-лея. Установлено, что зависимости радиуса сжимающего поршня и плотности на^нем.от времени близки к степенным.
3) Показано, что для всех рассмотренных. показателей (адиабаты, у чисто газодинамическими средствами возможно? получение плотностей" и температурнеобходимых для реализации управляемого инерциальнош термоядерного синтеза. .
4). Рассмотрены задачи о схождении ударной волны к центру или оси симметрии и о схлопывании пустой сферической полости в постановке, допускающей изменение показателя-адиабаты на фронте ударной волны, идущей от центра симметрии. Эти задачи обнаружили ряд неожиданных, интересных и важных моментов, возникающих в модели с изменением показателя адиабаты на ударной волне. В построенных решениях самые слабые («звуковые») ударные волны, которые при уменьшении показателя адиабаты на фронте ударной волны движутся по газу перед ними со скоростью звука, имеют конечную интенсивность. Более того, при у+ = 1 на такой звуковой ударной волне при небольшом (всего двукратном) повышении давления отношение плотностей, как на сильной ударной волне, бесконечно. Скорость ударных волн меньшей интенсивности дозвуковая, из-за чего они, будучи неэволюционными, запрещены.
5) В задачах о схождении ударной волны к центру или оси симметрии и схлопывании пустой сферической полости автомодельные решения с повышением показателя адиабаты на фронте ударной волны, идущей от центра или оси симметрии, запрещены из-за уменьшения энтропии. Типичны ситуации, в которых сохранение автомодельности возможно только при введении поршня, расширяющегося от центра симметрии. При его отсутствии течение за отраженной ударной волной становится неавтомодельным с момента отражения.