ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 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