БЛОК-СХЕМА. α, Н0, N, М, λ(θ), с(θ),ρ(θ), θ1 Ввод начальных данных и

Начало

α, Н0, N, М, λ(θ), с(θ),ρ(θ), θ1

Ввод начальных данных и

описание термодинамических

характеристик грунта.

- - - - -

Вычисление параметров

разностной схемы и

начальной функций

Δt, Δh, У[l], l =0, …, N У1[0] = θ1

- - - - -

конец

J = 0, M - 1, 1

E, αN-1, β N-1

Начальные значения

коэффициента прогонки

- - - - - -

I = N-1, 1, -1

αi-1, β i-1

I = 0, N - 2, 1

рис. 14

Если M и N достаточно большие величины, то в оперативной памяти компьютера может не хватит места для массива . Чтобы избежать этого, вводятся одномерные массивы . Вместо массивов A_i, B _i, C _i используются идентификаторы A_, B_, C. Для отводятся одномерные массивы . В формуле (7.16) при определений A_i,C_i,B _i используется отношение . Если это выражение очень большое, то вычислительный процесс будет не устойчивой. В этом случае не выполняется теорема 2.

Теорема 2. Пусть решение дифференциальной задачи (7.4)-(7.7) θ обладает непрерывными производными до четвертого порядка и если , то решение разностной схемы (7.8) – (7.10) сходится к решению дифференциальной задачи (7.4) – (7.7).

Алгоритм реализации схемы расчета приведен на рис. 3.