Научный журнал КубГАУ, No133(09), 2017 года 
8
одномерным уравнением теплопроводности [1]:
дТ 
д’Т
Раба д_і = К 
д2261 `
(10)
ТСПЛОСМКОСТЬ, 
плотность 
и 
теплопроводность 
моренного 
материала
рассчитываются, исходя из представления о моренном субстрате‚ как о
смеси породы и воздуха (и, вообще говоря, воды) [1]:
са = с,(1- 
к) + с„ 
к
=р,(1-к)+р,к. 
(11)
Ка = К, (1-к)+ К„ 
к
Поскольку вторые слагаемые в правых частях (11) малы по сравнению с
тремя первыми, они могут быть опущены.
Уравнение 
(10) 
решается 
численно, 
безусловно, 
устойчивым
методом прогонки. Краевое условие (8) в конечно-разностной форме
запишется как
Тат:2 — Т$
С, = К- —————>|, 
(12)
Л:
где т — номер узла в конечно-разностном разбиении моренного слоя,
т=1,...„М ’ (рис. 3). Поскольку температура поверхности Т; входит в
уравнения (2) и (6) для расчета потоков тепла на поверхности моренного
чехла, то для определения граничного условия необходимо применить
итеративный подход. Для этого используется метод Ньютона-Рафсона:
та =т? -Л( 
)) Р() 
(13)
М-+1 
‚ 
М
где Т, — - значение температуры на итерации М+1, No=1,...‚„Миах; / (Тз ) - это
сумма функций в правой часть уравнения (1), а / ’(Т$No ) - соответственно ее
производная.
Бир://е).Кибаето.ги/2017/09/раг /78.раё