днеперсивности с увеличением пути переноса; метод требует ми-
нимальных затрат для моделирования переноса многокомпонент-
ных и взаимосвязанных мигрантов; при повторных расчетах оди-
наковых примеров отчетливо виден стохастический характер, так
что может быть полезнес, например, решать задачу 10 раз © 200
блуждающими точками, чём один раз с 2000
К недостаткам метода случайных блужданий можно отнести
необходимость сравнительно большого числа блуждающих точек,
тарантирующего 
достаточно 
гладкий 
результат, 
а 
также 
то,
что при выборе сравннтельно грубой дискретизации внутри обла-
сти_концентрация может оказаться больше граничных значений.
Примеры 
вычислительных 
программ 
по 
методу 
случайных
блужданий приведены_ в работе [13]_с использованием_для_ реше-
ния фильтрационной задачи неявной конечно-разностной аппрок-
симацни на прямоугольной сетке и неявного метода переменных
направлений:
Метод конечных элементов для моделирования миграции почти
исключительно использует алгоритм, базнрующийся на. формули-
ровке Галеркина. Хорошо обоснованной, отработанной и опробо-
ванной_считается программа РЕЕ(О\, представленная в работах
[50, 51]. Эта программа охватывает одновременно определение
поля скоростей и концентраций конечно-элементной аппроксима-
цией. Можно выбирать элементы различной формы (рис. 27), од-
нако` прежде всего оказываются пригодными изопараметрические
криволинейные четырехугольные элементы с биквадратным мате-
матическим выражением (восьмиточечный элемент, м, рис. 27, г).
В качестве первичных_неизвестных вводятся как Н, с, так н
естественные переменные Н, о, Фу, с. При формулировке фильтра-
ционной задачи в напорах необходимые в дальнейшем скорости
Ффильтрации получаются из поля напоров путем дифференцирова-
ния. Выводимые при этом математические функции для ©; 
И 
бу
таким образом сще только билинсарные, Поскольку такая форму-
лировка не обеспечивает устойчивого состояния ДЛЯ 0 н 0,, 0Со-
бенно при возможности значительных различий скоростей от эле-
мента к элементу, то рекомендустся формулнровка, включающая
в первичные переменные скорости ох н оу. В четырехугольном эле-
менте определяются дискретные величины Н, о, Фу, © В угловых
точках и значения ,, Фу» © В точках середины сторон. Такая ап-
@
6 
в
Рис. 27. Схемы изопараметрических конечных элементов плоского потока,
‚Элементы: 
—- треуголыимй е лицейными соотнопечнями; 6 — четырехугольный с билнней-
чыми 
соотношениями: 
@ — треугольный 
с 
квадратичными 
соотмОШеннямн 
е =— четырек-
Угольный © бикоадратимми соотношеннями
г
131