Регуляризация схем второго порядка точности

Таким образом, применение соотношений типа (3.111) основано на том, что элемент, представляемый явной схемой Эйлера в методе Рунге — Кутта, заменяется на неявный элемент, разрешаемый Ньютоновскими итерациями. Конкретный выбор значений параметров в (3.111) определяется процедурой регуляризации, состояш ей в установлении соответствия между численным решением и формальным разложением в ряд Тейлора с заданным порядком точности по к (порядок не может быть больше второго). Применяя формулы вычислительного процесса У п+1 = ФУп к исходному уравнению у = —Ку, всегда можно удовлетворить требованию ф < 1 выбором значений параметров в (3.111). Другие параметры выбираются либо пз сообра-жеиий простоты процедуры, либо регуляризацией иного типа, наделяющей численную схему дополнительными желательными свойствами. Таким образом, вычислительный процесс (3.102) легко управляем и является балансным, однако не имеет свойства положительности, т. е. в решении возможно появление отрицательных концентраций, продемонстрированное на примере (3.83). [c.188]

Справочник химика 21

Химия и химическая технология