по схеме Галеркина с линейными функциями. Вместе с тем неяв- ная схема, обратные разности членов конвекции и аккумуляции, а также формулировка МКЭ по Ритцу и по Галеркину с много- кратной коллокацией в значительной мере свободны от осцилля- ций. Вместе с тем чем «центральнее» численная схема, тем она точнее и чувствительнее к нарушениям. Поэтому применяемая на практике численная схема постоянно является компромиссом. меж- ду обенми тенденциями. `Наряду с погрешностями дискретности имеют значение также погрешности в стабильности, вытекающие из ограниченного коли- чества численных расчетов. Безусловно стабильной считается чис- ленная модель миграции, если численная погрешность (округле- ния) уменьшается от одного временнбго шага к другому, а условно стабильной — если это происходит только при определенных усло- внях. Эти условня для особых случаев аналитически представлены в работе [50]. Таким образом, сравнением © аналитическими ре- шениями фиксируется условие стабильности при заданной прост- ранственной дискретизации путем установления критической вели- чины временнбго шага через критические числа О1 или Сг. Без- условно стабильной_ является неявная схема решения с у== 1, при- чем с уменьшением у возрастает склонность к нестабильности. Численное решение составленной системы уравнений (матричное уравнение) также тант в себе возможности погрешностей. К очень большим погрешностям, сильно распространяющимся при услов- ном стабильном методе, может приводить решение системы' урав- нений с плохо выраженными условиями, у которых элементы ос- новных днагоналей матрицы коэффициентов_ в недостаточной сте- пени преобладают по сравнению с основными днагоналями кодиа- гоналей. Значительные погрешности в решении уравнений может вызы- вать решение всей системы уравнений с помощью частного метода шагов (например, неявного метода переменных направлений) и пе-- реноса элементов матрицы коэффициентов в правую Часть урав- мений путем.умножения на временные или итеративно зависимые переменные © обратной датировкой для создания ленточных мат- риц с незначительной шириной ленты (преимущественно триднаго- нальные матрицы коэффициентов). По этой причине следует тща- тельно проверять и контролнровать программы компьютера по- численному моделированию миграции, особенно путем сравнения © аналитическими решениями. `На_основе численного решения_производится первичное опреде- ление числа опорных точек в пространственно-временнбй сетке. Число опорных точек по времени или по размеру нтерационного- шага при нелинейном решении указывает количество определяе- мых локально-дискретных значений зависимых переменных (Н или иногда о;, Э,, ©) и таким образом влияет на число уравнений си- стемы. Затраты времени на одноразовое решенне Этой снстемы уравнений являются основной величиной для оценки затрат; они. зависят от типа ЭВМ, используемого метода для решения снстемы 124