Файл: Чайлдс Э. Физические основы гидрологии почв.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

Просмотров: 253

Скачиваний: 5

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

12.6.Горизонтальная диффузия воды

вполубесконечную колонку с постоянной исходной влажностью

Рассматриваемый здесь случай сам по себе не представляет боль­ шого интереса, но его можно считать введением к более важному, но и более трудному случаю развития вертикального профиля влаж­ ности. Кроме того, как указывалось в параграфе 11.2, уравнения, описывающие горизонтальное движение влаги, приблизительно верны и для вертикального передвижения, если мала влагопровод­ ность, свойственная почве, как у тяжелых глин, или, что является более общим условием, если мала влажность почвы на всем рас­ сматриваемом участке профиля.

Конкретно мы рассмотрим следующую задачу. Имеется горизон­ тальная колонка пористого материала, ось которой совпадает с на­ правлением оси X, берущей начало на торце колонки. В направлении положительных значений х колонка считается бесконечной. Сначала влажность во всех точках колонки одинакова и равна с0. В некото­ рый момент t, который можно принять за нуль при отсчете времени, на торце колонки создается другая влажность С, которая и под­ держивается в дальнейшем. Этого можно добиться, например, за­ топлением свободной поверхности колонки сухой почвы и поддер­ жанием состояния насыщенности на этой поверхности. Можно под­ держать и другую постоянную влажность, хотя это связано с техническими трудностями, о которых мы здесь не будем говорить, но которые вполне преодолимы. Обсуждаемая задача состоит в том, чтобы рассчитать распределение влажности в колонке для любого момента времени.

Соответствующим уравнением неразрывности, решение которого содержит необходимую нам информацию, является уравнение (11.29)

dc/dt — d[D(dc/dx)\/dx.

(12.25)

Граничные условия, характеризующие данную задачу,

состоят

в том, что концентрация (влажность) с для всех значений х, отличных от а:0, в нулевой момент времени равна с0; при х, равном нулю, влаж­ ность равна С для любых t. Интуитивно можно понять, что после внезапного затопления поверхности требуется бесконечно большое время для того, чобы затопление ощутилось на бесконечно большом расстоянии от начала колонки; понятно также, что через бесконечно большой срок все точки, находящиеся на конечном расстоянии от начала колонки, будут иметь одну и ту же влажность С. Следова­ тельно, мы знаем, что при х, равном бесконечности, с равно с0 для всех конечных отрезков времени, а также что при бесконечно боль­ шом t с равно С для всех конечных х. Эти граничные условия можно сформулировать следующим образом:

с = с0

при £ = 0 (1 /г = оо) и £ > 0 ,

 

с = с0 при

1/х = 0 {х —оо) и t <

оо (1/t >

0 ),

с —С при х = 0 (1/х = оо)

и t~> 0 ,

 

с = С при

1/£ = 0(£ = оо) и X <

оо (1/х >

0).


Математическая эквивалентность нулевого расстояния и беско­ нечного времени, нулевого времени и бесконечно большого рас­ стояния настолько бросается в глаза, что это побудило Больцмана объединить обе независимые переменные х и t в единую перемен­ ную % соотношением

1

X = */**•

(12.26)

Подстановка этой переменной позволяет преобразовать дифферен­ циальное уравнение в частных производных (12.25) в обыкновенное дифференциальное уравнение. Из уравнения (12.26) следует, что

1

dxldx --=l/t2

(12.27)

д% /dt =

Теперь уравнение (12.25) можно переписать следующим образом:

{d%/dt) (dc/dx) = (дхідх)* d [D (dc/d%)]/dx.

Подставив соответствующие выражения для частных производ­ ных из уравнений (12.27), получим

( Y x/t 2 ] dc/dx = d [D(dc/dx)]/dx.

Инаконец, используя уравнение (12.26), найдем

X dc/dx + 2d [D(dc/dx)]/dx = 0.

(12.28)

Данное уравнение имеет только одну независимую переменную % и зависимую переменную с. Это означает, что его решением, если такое существует, является кривая или математическое выражение, представляющее с в виде функции переменной %. Рассматривая уравнение (12.26), можно заметить, что для данного фиксированного момента времени кривую зависимости с от х можно считать кривой зависимости с от х, где масштаб х зависит от выбранного момента времени; чем больше отрезок времени, тем больше расстояние х при данных значениях р с, т. е. тем глубже проникновение в ко­ лонку данной влажности с. Точно так же для фиксированного х

кривую можно рассматривать как график зависимости с от t i/2, где масштаб зависит от выбранного значения х.

В опытах по исследованию развития профиля влажности в раз­ личные моменты времени одна и та же влажность обнаруживается на различных расстояниях от начала колонки, однако если движение влаги соответствует уравнению (12.28), то отношение пройденного расстояния к корню квадратному из времени будет постоянным. Отсюда следует, что если такая инвариантность действительно наблюдается, она служит доказательством, что движение влаги находится в соответствии с уравнением диффузии, безотносительно

к тому, как зависит от

влажности D, поскольку вид этой зависи­

мости не играет роли

при выводе уравнения (12.28) из уравне­

ния (12.25).

 


Решение же уравнения (12.28), очевидно, зависит от характера связи между D и влажностью, и пока эта связь не уточнена, решение не может быть получено. Однако природу этого решения до некото­ рой степени можно пояснить хотя бы в самом общем виде. Прежде всего развернем уравнение (12.28) и запишем

 

 

 

de

 

 

 

 

 

О,

 

 

 

(12.29)

 

 

 

1 dÿ.

 

 

 

 

 

 

 

 

где

dDJdc — параметр, характеризующий

изменение

коэффициента

диффузии

с

влажностью.

В

частном

случае,

когда

коэффициент

(с-с0)/(С-с0)

 

 

 

 

диффузии постоянен и потому

 

 

 

 

не

зависит

от

влажности,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

второй

член

в

уравнении

 

 

 

 

 

 

 

(12.29)

исчезает,

а

из

остав­

 

 

 

 

 

 

 

шейся

части уравнения сле­

 

 

 

 

 

 

 

дует, что

если

первый член

 

 

 

 

 

 

 

отрицателен, то третий поло­

 

 

 

 

 

 

 

жителен, и обратно. Посколь­

 

 

 

 

 

 

 

ку % и D существенно поло­

 

 

 

 

 

 

 

жительны,

то

когда

dc/d%

 

 

 

 

 

 

 

положительно,

d2c/d%2

дол­

 

 

 

 

 

 

 

жно быть отрицательно. Сле­

 

 

 

 

 

 

 

довательно, график зависимо­

 

 

 

 

 

 

 

сти с от %должен

иметь вна­

 

 

 

 

 

 

 

чале

наибольшую крутизну,

 

 

 

 

 

 

 

а затем становиться все ме­

 

 

 

 

 

 

 

нее и менее крутым по мере

 

 

 

 

 

 

 

увеличения %. Для этого слу­

Рис.

12.8.

Распределение влажности с

 

чая

существует

аналитиче­

 

 

в функции от

 

 

 

ское

 

решение,

 

приводимое

Вода диффундирует от поверхности, имеющей

 

в соответствующих

руковод­

влажность

С,

в полубесконечный

профиль,

 

влажность которого вначале всюду равна с0. D

 

ствах (см.

[2 1 ]):

 

 

 

 

постоянный коэффициент диффузии.

 

 

 

 

 

 

 

 

 

 

гд е

 

 

С — С0=

(С — с0) erf C(X /2Z>‘/ 2),

 

 

 

 

(1 2 .3 0 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

00

 

 

 

 

 

 

 

 

 

 

 

 

 

erfc(z) = (2 /л1/г) J ехр( —г|2) dr]

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

— табулированная функция.

Решение

уравнения

(12.30)

предста­

влено на рис. 1 2 .8 в форме графика

зависимости

(с — с0)/(С — с0)

от (х/2D'/*),

демонстрирующего отмеченную выше особенность. Если

вспомнить, что при нахождении этого решения мы полагали D по­

стоянным, то рис. 1 2 .8 можно рассматривать как график зависи­ мости между влажностью и расстоянием от поверхности для данной длительности процесса и данного значения D, величиной которых определяется масштаб оси абсцисс. Выбранный на рисунке масштаб

непосредственно характеризует

х для всех комбинаций D u t , для

которых произведение 2

равно единице.


Можно задаться другими формами профиля влажности и пы­ таться выяснить, каким зависимостям коэффициента диффузии от влажности эти формы соответствуют. Простым случаем является такой, при котором влажность уменьшается линейно от С на поверх­ ности до с0 на расстоянии %0, а затем остается равной с0при всех больших значениях %, как показано на рис. 12.9. Это решение можно записать в форме

 

 

 

с —С— (С— с0) (х/Хо) при 0 <% <%о>

(12.31)

 

 

 

 

 

с = с„ При %>Хо.

(12.32)

Из

 

уравнения

(12.31)

 

 

 

имеем

 

 

 

 

 

 

 

deld%= (С— с0)/%о

 

 

 

d2c/dx2 = 0

 

(12.33)

 

 

 

 

 

 

 

 

 

 

 

а из уравнения (12.32) по­

 

 

 

лучаем в интервале % боль­

 

 

 

ших,

чем

%0,

 

 

 

 

 

dc/dx = d2c/dx2 — 0.

 

(12.34)

 

 

 

Подставляя

уравнение

Рис.

12.9. Гипотетическое распределение

(12.34)

в

уравнение

(12.28),

сразу

видим, что

уравнение

Такое

влажности с в функции от х-

диффузии

удовлетворяется

распределение характеризует

зависимость

коэффициента диффузии от влажности.

при

X

больших,

 

чем

 

 

 

Подставив уравнение (12.33) в уравнение (12.28), получим

 

 

 

 

 

 

dD/dx =

—(х/2),

(12.35)

после интегрирования которого имеем

 

 

 

 

 

 

D = A - x V 4,

(12.36)

где А — постоянная интегрирования, которую надо найти. В точке перехода от линейной части графика к постоянному значению с, а именно при х = Хо> величина d2c/dx2 стремится к бесконечности, поскольку длина зоны перехода становится бесконечно малой, и раз­ вернутой форме уравнения (12.28), т. е. уравнению (12.29), можно удовлетворить только в том случае, если D в этой точке равно нулю. Следовательно, для этой точки уравнение (12.36) принимает вид

0 = Л - х 2о/4,

или

-4 = Хо/4,

так что общее выражение для D при этом значении А , согласно урав­ нению (12.36), имеет вид

Д = (Х2о -Х 2)/4.

(12.37)

219