Математическое моделирование комплексных процессов / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Вычисл. центр ; [редкол.: В. С. Мингалев (отв. ред.) и др.]. - Апатиты : Кольский филиал АН СССР, 1982.

где E q - внешнее постоянное электрическое поле, а Е обусловлено локальным нарушением квазинейтральности плазмы в волне. Для этой величины имеем d i v Е = 47 t e ( n . - n ) 4 1 е ' r o t Е = О. Учитывая, что при сильных электрических полях в верхней части Е-слоя может нарушаться неравенство со « і)е , в целях сохранения общности модели, слагаемое 3g^eY 9 t сохраним в уравнении дви жения. Оценки второго тока, сделанные в работе / 9 / , показали, что вклад электронов существенно выше вклада ионов в нелинейное взаи модействие волн Фали. По этой причине можно пренебречь нелиней ностями, входящими в уравнения ионной компоненты. Изотцрмичность электронов исключает из рассмотрения нелинейность g \ e / E g , т.е. в итоге в_качестве рабочей нелинейности модели будет исполь зоваться п е Е - нелинейность уравнения движения для электро нов. Такой шаг, как мы увидим ниже, существенно сокращает рас ход машинного времени и объема памяти. Однако при этом исключе ние температурных нелинейностей приводит и к нежелательному ис чезновению части квазилинейных механизмов стабилизации, к кото рым относятся коллективный нагрев электронов и ионов. Кроме это го, слагаемое gv е / E s обеспечивает эффект джоулева нагрева. Для компенсации этих недостатков необходимо ввести в модель учет температурных эффектов. Это можно сделать, например, путем ис пользования специальной сервисной подпрограммы. Для построения температурной подпрограммы можно воспользоваться результатами исследования квазилинейных процессов стабилизации / 1 0 / . Введение такой программы должно приблизить рассчитываемую модель к реальным условиям в полярной ионосфере и тем самым сократить необходимое число вариантов расчета. С физической точки зрения это выделение температурных вариаций можно связать с различием во временных масштабах процессов, т.е. с выполнением неравенства Введя обозначения 1 n « jn ( r , t ) e x p ( — i k r ) dr, g = j g ( r ,t ) * 4 e x p ( - i k r ) d r, z q = B j B ^ , перейдем к спектральному пред ставлению системы для электронов. Тогда имеем: Эп / х