Методы исследования устойчивости

Существуют три стандартных метода, обычно используемых для исследования устойчивости или свойств неустойчивостей, связанных с данным равновесным состоянием плазмы.

Простейший из них —интуитивный подход, согласно которому равновесное состояние подвергается возмущению, изменяющему силы, действующие на плазму. Если эти измененные силы действуют так, что увеличивают первоначальное возмущение, плазма неустойчива. Примером применения интуитивного подхода служит анализ змейковой неустойчивости. Этот подход не позволяет получить количественную информацию об инкрементах, однако с его помощью удается вскрыть механизм неустойчивости, показать, какие равновесные ситуации требуют дальнейшего анализа, а также получить информацию о том, какие моды, по-видимому, должны быть неустойчивыми.

Второй метод основан на энергетическом принципе. Равновесное состояние плазмы с потенциальной энергией W устойчиво, если W представляет собой минимальное значение, и неустойчиво в противном случае. Для исследования устойчивости плазмы этим методом необходимо вычислить изменение потенциальной энергии плазмы в результате данного возмущения. Плазменное равновесие устойчиво, если изменение потенциальной энергии положительно при любых возмущениях, не противоречащих уравнениям плазмы. Если же существует хотя бы одно допустимое возмущение, для которого приращение потенциальной энергии отрицательно, плазма неустойчива. Этот метод удобен при определении равновесных параметров в случае неустойчивого равновесия. Он не пригоден для нахождения инкремента конкретной неустойчивости.

Третий метод исследования устойчивости равновесия плазмы — анализ собственных частот. В этом методе предполагается, что, если исследуемое равновесное состояние подвергнуто возмущению, линеаризованные уравнения плазмы, описывающие развитие возмущения во времени, могут быть решены при соответствующих граничных условиях в предположении временной зависимости вида . При этом мы получаем уравнение дляω, в которое входят равновесные параметры. Решения этого уравнениямогут быть действительными, мнимыми или комплексными. Если все решения дляω действительны, то любые переменные, описывающие возмущение, осциллируют гармонически и плазма устойчива. Если хотя бы одно из решений дляω имеет положительную мнимую часть, система неустойчива, поскольку соответствующее собственное колебание будет раскачиваться со временем.

Анализ собственных частот позволяет получить полную информацию о неустойчивостях, присущих конкретному равновесному состоянию плазмы. Развитие любого начального возмущения может быть прослежено до того предела, до которого справедливы линеаризованные уравнения. К сожалению, анализ собственных частот может быть применен лишь в тех случаях, когда равновесное состояние плазмы достаточно простое, чтобы дифференциальные уравнения плазмы могли быть решены для каждой из гибридных частот. Уменьшение показателя преломления с возрастанием частоты называется в оптике аномальной дисперсией, причем в наиболее характерных случаях аномальной дисперсии это уменьшение происходит резко, в узком интервале частот. Именно поэтому гибридные частоты и следует называть частотами аномальной дисперсии.

(источник:Н.Кролл, А.Трайвелпис, Основы физики плазмы, гл.5, §§1-3)

Пучковая неустойчивость

Проведем анализ одного важного следствия формулы (31.23) для декремента затухания Ландау плазменных колебаний. Из этой формулы следует, что затухание колебаний, обусловленное обменом энергий с резонансными частицами, реализуется только в том случае, когда функция распределения резонансных частиц есть монотонно убывающая функция скорости. Это имеет место, например, при максвелловском распределении частиц плазмы по скоростям, когда единственный максимум функции распределения соответствует скорости

При наличии в плазме достаточно интенсивного пучка электронов функция распределения по скоростям может иметь дополнительный максимум для значения , отличного от нуля. Этот случай изображен на рис. 1.17. Волна, фазовая скорость которой лежит вблизи второго максимума функции распределения между точками 1 и 2, раскачивается электронным лучком. Механизм раскачки — все тот же резонансный обмен энергией между волной и частицами, который в данном случае за выполнения условия , т. е. преобладания в исходном состоянии частиц, движущихся быстрее волны, приводит к росту амплитуды волны. Таким образом, из всего спектра волн, которые вследствие флуктуации всегда присутствуют в плазме, электронный пучок выделяет узкую спектральную область находящихся с ним в фазовом резонансе волн и интенсивно накачивает энергию в эти волны.

В определенном смысле электронный пучок в плазме можно рассматривать как классический мазер на продольных колебаниях, поскольку фактически механизм неустойчивости связан с инверсностью заселенности уровней для частиц со скоростями, близкими к фазовой скорости волны (условие ).

Скорость нарастания амплитуды плазменных волн, возбуждаемых электронным пучком, определяется формулой (1.148). Из этой формулы, в частности, следует, что для возникновения неустойчивости необходимо, чтобы интенсивность пучка в плазме превышала некоторую минимальную величину или чтобы пучок был сильно надтепловым, т. е. приходился на далекую область «хвоста» максвелловского распределения тепловых частиц плазмы. Очень слабый пучок, очевидно, не сможет повлиять на функцию распределения электронов по скоростям таким образом, чтобы ее производная могла изменить знак. Если же в области скоростей (резонансных с неустойчивыми волнами) можно пренебречь вкладом плазменных частиц, то очевидно, что функцию распределения в резонансной области скоростей можно представить в виде , где —плотность пучка; — тепловой разброс в пучке (ширина его функции распределения по скоростям). Производная от функции распределения , и в этом случае для инкремента нарастания находим из (31.23) формулу

которая часто используется для оценок.

Формулы (31.23) и (32.1) описывают неустойчивости пучков сбольшем тепловым разбросом. Дело в том, что за нестационарности волнового процесса, в данном случае роста амплитуды волны при неустойчивости, резонанс волна — частица имеет конечную ширину по скорости .Именно на такую минимальную ширину по фазовым скоростям размыт пакет волн из-за мнимой части частоты. Формула (32.1), при выводе которой пренебрегалось конечной шириной резонанса волна — частица, относится к случаю сильно размытых пучков, когда тепловой разброс в пучке существенно больше ширины резонанса и для каждой неустойчивой волны на функции распределения пучка выделяется своя небольшая группа резонансных частиц.

В обратном предельном случае весь пучок как целое находится в резонансе с неустойчивой волной. Именно в этом случае следует ожидать развития наиболее сильной пучковой неустойчивости. Дисперсионное уравнение для такой неустойчивости можно получить из общего уравнения. При вычислении интеграла по скоростям в этом уравнении следует выделить область тепловых скоростей, в которой интеграл по v вычисляется обычным путем (разложением по параметру ), и область скоростей пучка, при рассмотрении которой можно считать пучок имеющим δ-образное распределение по скоростям: Однако нагляднее вывести искомое уравнение по аналогии с дисперсионным уравнением неустойчивости Бунемана, считая электроны пучка своего рода отдельным сортом частиц со своей плазменной частотой

(32.2)

где — квадрат ленгмюровской частоты в пучке. Применим графический метод анализа корней этого уравнения. График правой части этого уравнения при фиксированном k, т. е. для возмущений с заданной длиной волны, приведен на рис. 1.18. Видно, что когда минимальное значение F оказывается меньше единицы, у уравнения (32.2) все четыре корня вещественны, что соответствует периодическому изменению возмущений со временем. Если же минимум F больше единицы, то уравнение (32.2) имеет только два вещественных корня, а два корня являются комплексно сопряженными . Эти корни, очевидно, соответствуют изменению амплитуды колебаний по закону , и, таким образом, один из этих корней приводит к росту амплитуды со временем, т. е. к неустойчивости. Из (32.2) имеем, что и, следовательно, неустойчивыми окажутся достаточно длинноволновые возмущения, для которых

Коротковолновые колебания с устойчивы. Примерный график функции можно представить следующим образом. Неустойчивы только волны, сохраняющие фазовый резонанс с пучком . При малых инкремент нарастает с волновым числом: