close

Вход

Забыли?

вход по аккаунту

?

Двумерная математическая модель промерзания засоленного влажного грунта.

код для вставкиСкачать
В. В. Попов. ДВУМЕРНАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОМЕРЗАНИЯ ЗАСОЛЕННОГО ВЛАЖНОГО ГРУНТА
УДК 519.63
В. В. Попов
ДВУМЕРНАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОМЕРЗАНИЯ
ЗАСОЛЕННОГО ВЛАЖНОГО ГРУНТА
Проведено исследование двумерной математической модели промерзания грунта, содержащего соленую воду. Граница
раздела талой и двухфазной области определяется заданием температуры фазового перехода, определение температуры
в статье не рассматривается. Приводится зависимость решения задачи от этой температуры. Мерзлая область, строго
говоря, отсутствует. Но, как показывают расчеты, вблизи дневной поверхности (зимой) влажность практически можно
считать равной нулю.
Решение задачи сводится к решению системы уравнений, которое состоит из трех уравнений с разрывными
коэффициентами: уравнение диффузии, уравнение теплопроводности и уравнение термодинамического равновесия фаз.
Для численной реализации системы алгебраических уравнений применяется итерационный метод сопряженных градиентов.
Линеаризация проводится методом простых итераций. В области фазовых переходов распределение влажности вычисляется
из уравнения диффузии, а концентрация вычисляется из уравнения термодинамического равновесия.
Вычислительные расчеты показали эффективность примененного метода, что обосновано быстрой сходимостью
итерационного процесса и возможностью распараллеливания разработанного алгоритма.
Приведены сравнения результатов разностных аппроксимаций уравнения диффузии в дивергентном и недивергентном
представлении. Оказалось, что обе аппроксимации дают практически одинаковые результаты.
Ключевые слова: фазовый переход, диффузия, теплопроводность, теплоемкость, температура, концентрация, влажность,
водонасыщенность, теплота, пористость, аппроксимация, дивергентный, недивергентный, засоленный, грунт.
V. V. Popov
Two-Dimensional Mathematical Model of Frost of Saline Wet Soil
The research of two-dimensional mathematical model of soil freezing containing salt water is conducted. The interface and
the melt-phase area are defined by specifying the phase transition temperature. The determination of the temperature in this article
is not considered. The dependence of the solution of the problem of the temperature is given. The icy area, strictly speaking, is
absent. However, calculations show that near the surface (in winter) the humidity is almost zero.
Solution of the problem is narrowed down to solving a system of equations of the three equations with discontinuous
coefficients – diffusion equation, heat equation and the equation of thermodynamic equilibrium phases. For the numerical
implementation of the system of algebraic equations iterative conjugate gradients method is applied. Linearization is conducted
by the method of simple iterations. In the region of phase transitions moisture distribution is calculated from the diffusion equation,
and the concentration is calculated from the equation of thermodynamic equilibrium.
Numerical calculations have shown the effectiveness of the method used, which is justified by the rapid convergence of the
iterative process and the possibility of parallelization of the developed algorithm.
Comparisons of the results of difference approximations of the diffusion equation in divergent and non divergent view are
represented. It was found that both approximations give almost identical results.
Key words: phase transition, diffusion, thermal conductivity, specific heat, temperature, concentration, humidity, water saturation,
warmth, porosity, approximation, divergent, non divergent, saline, soil.
ПОПОВ Василий Васильевич – к. ф.-м. н., доцентисследователь кафедры вычислительных технологий ИМИ
СВФУ им. М.К. Аммосова.
E-mail: imi.pm.pvv@mail.ru
POPOV Vasily Vasilievich – Candidate of PhysicalMathematical Sciences, Associate Professor of the Department
of Computing Research, the Institute of Mathematics and
Informatics, the North-Eastern Federal University named
after M.K. Ammosov.
E-mail: imi.pm.pvv@mail.ru
Введение
Для описания процессов промерзания и протаивания
разработано множество математических моделей. В
работах [1-2] предложены математические модели с
фазовым переходом в протяженной области. Получены
автомодельные решения. Численные решения этих
моделей получены в работах [3-4]. В этих моделях
характерным является выделение границы фазового
перехода в виде уравнения, связывающего скорость
движения границы с разностью потоков. Математи19
ВЕСТНИК СВФУ, 2014, том 11, № 5
ческие модели без явного описания границ области
фазовых переходов разработаны в [5-6]. Многомерные
математические модели с фазовым переходом в
протяженной области еще недостаточно изучены и
численно не реализованы с применением современных
вычислительных методов и многопроцессорных систем.
Постановка задачи
Распределение температуры описывается уравнением теплопроводности [1-6]:
C
∂T ∂  ∂T  ∂  ∂T 
∂ν
 + mρw qw
= λ
,
 + λ


∂t ∂x  ∂x  ∂y  ∂y 
∂t
0< x < l1 , 0< y < l2 , t >0.
В области фазового перехода
концентрации
соли
описывается
уравнением диффузии [1-6]:
распределение
дивергентным
∂c  ∂c 
∂ν c
∂c  ∂c 
= D ν  + D ν  , T < T f , t >0,
∂y  ∂y 
∂t
∂x  ∂x 
и
выполняется
равновесия [1-6]
где полагаем
уравнение
термодинамического
на нижней границе:
T = T 0 , 0 < x < l1 , y = l2 ;
на боковых границах области:
∂T
=0, x=0, x = l1 , 0 ≤ y ≤ l2 .
∂x






C + mρw qw ν  ∂T = ∂ λ ∂T  + mρw qw D ∂ ν ∂T  +





T  ∂t ∂x  ∂x 
T ∂x  ∂x 
+
(4)
Теплоемкость и коэффициент теплопроводности
влажных пористых сред могут быть определены по
формулам [1-6]:
0
C = ( 1 - m)Cs ρs + mCl ρl (ν -ν )ρw / ρl + mCw ρwν ,
λ = ( 1 - m)λs + mλl (ν 0 -ν )ρw / ρl + mλwν ,
где нижние индексы s,l,w относятся соответственно к
скелету, льду и воде, а верхний индекс 0 – начальному,
m – пористости, ν – водонасыщенности.
Начальные условия:
T(x, y ,0 )= T 0 , c(x, y ,0 )= c0 , 0 £ x £ l1 , 0 £ y £ l2 .
(5)
Граничные условия для температуры на верхней
границе области:
(8)
где Тp , Т0 , Тa – соответственно температура под
зданием, на поверхности контакта с атмосферой,
температура воздуха, lp – длина здания.
Граничные условия для концентрации не требуются,
т. к. из уравнения диффузии будем определять
распределение влажности.
Численная реализация
Сначала приведем уравнения в виде, удобном
для проведения численных рассчетов. В уравнении
(1) исключим производную по времени от
водонасыщенности, используя для этого уравнение
диффузии (2):
c = -T / α , T < T f , t >0,
Тf =-c0a-e.
(7)
∂  ∂T 
D ∂  ∂T 
 + mρw qw
.
λ
ν
∂y  ∂y 
T ∂y  ∂y 
Дивергентный алгоритм
Уравнение диффузии (2) аппроксимируем
зависимости от направления потока:
 
ν ij cij -ν iJ cij
τ
Здесь
если
= ai +bi + aj +bj.
(9)
в
(10)
ci+1 j - cij <0,
ai = Dν ij
ci+1 j - cij
2
1
h
, bi = Dν i-1j
ci-1 j - cij
h12
,
(11)
,
(12)
в противном случае
ai = Dν i+1j
если
ci+1 j - cij
2
1
h
, bi = Dν ij
ci-1 j - cij
h12
cij+1 - cij <0,
cij+1 - cij
cij-1 - cij
aj = Dν ij
, bj = Dν ij-1
, (13)
2
h2
h22
в противном случае
l p p/ 2/ 2,< yx =
<0,
l1 / 2 +l p / 2, y = 0,
T = Tp , l1 / 2 − l p T
/ 2=<Txp ,<l1 l/1 2/ 2−+l
c - cij
cij-1 - cij
=−
T0l ,p0/ ≤
/ 2p −
l1 l/1 2, y+l
, y==D
0,ν ij+1 ij+1
, 0l≤/ x2 ≤
l1 T// 22,
2, lx1 ≤
/ 2l1+l
/ 2l p≤/ 2,
x≤
=p0,/ 2 ≤ x ≤ l1aj
, (14)
, bj = Dν ij
Tp , l1 / 2 − l p /T2 =<Tx0 <
+l
y
=
0,
2
1
p
h2
h22
или
∂T
T0 , 0 ≤ x ≤ l1 /λ2∂−Tl =
2,(Tl1 −
/ 2Ta+l
≤l1−
l1/,T2ya −
=
=≤α
(T
), l0,p0/ ≤
λ p /02≤
/ 2p −
l1при
),
xx≤
x≤
l/1 2, y+l
/ 2l1+l
/(6)
=p0;/ 2 ≤ x ≤ l1 , y = 0;
2, lx1 ≤
2l p≤/ 2,
p /α
∂y
∂y
Ti+1 j ,Tij ,Ti−1j ,Tij−1 ,Tij+1 < −α ⋅ c0 .
(15)
T
= α (T − Ta ), 0 ≤ x ≤ l1 / 2 − l p / 2, l1 / 2 +l p / 2 ≤ x ≤ l1 , y = 0;
y
20
В. В. Попов. ДВУМЕРНАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОМЕРЗАНИЯ ЗАСОЛЕННОГО ВЛАЖНОГО ГРУНТА
Недивергентный алгоритм
Уравнение диффузии (2) перепишем в недивергентном виде:
2
Полученное уравнение аппроксимируем с учетом
направления потока:
Ä
= ai + Dν ij
+aj + Dν ij
если
ci+1 j - 2 cij + ci-1j
h12
cij+1 - 2 cij + cij-1
2
2
h
(17)
,
h1
;
cij+1 - cij <0,
ν ij -ν ij-1 cij+1 - cij
;
aj = D
h2
h2
cij+1 - cij >0 ,
ν ij+1 -ν ij cij+1 - cij
h2
Cs ρs =1.86 e6; Cl ρl =1.80 e6; Cw ρw =4.19 e6;
(18)
qw =3.34 e4; m=0.2 -0.35; D=1.45 e - 9;
c0 =0.003,ν 0 =0.9; l p =2 ,T 0 = 4; Tp =4.
ν i+1 j -ν ij ci+1 j - cij
h1
Численная реализация алгоритма основана на
методе сопряженных градиентов [7]. Алгоритм
распараллелен с помощью системы передачи данных
MPI [8].
Результаты вычислений
Расчеты проведены при следующих значениях
параметров:
ρw =1000; ρl =910; α=66.7;
ci+1 j - cij >0,
aj = D
(22)
λs =2.09; λl =2.23; λw =0.58;
ci+1 j - cij <0,
ν ij -ν i-1 j ci+1 j - cij
;
ai = D
h1
h1
ai = D
Ti+1 j ,Tij ,Ti-1j ,Tij-1 ,Tij+1 < T f .
2
¶ν c
¶ν ¶c
¶ c
¶ν ¶c
¶ c
=D
+ Dν 2 + D
+ Dν 2 . (16)
¶t
¶x ¶x
¶x
¶y ¶y
¶y
 
ν ij cij -ν iJ cij
при
h2
Рис. 1. Распределение концентрации
при различных значениях e
(19)
(20)
;
(21)
2p
p
Ta = (T 0 − 25) sin( + t
) − 25
2
86400 ⋅ 210
.
На рис. 1 приведены результаты расчетов при e = 0,
0.1 и 0.2.
На рис. 2 приведены сравнения одномерной задачи с
двухмерной задачей при одинаковых условиях. На рис. 3
приведено сравнение распределений водонасыщенности
по глубине вблизи границы x = 0 на разных сетках
n = 200 ; n=400; n=800 . На рис. 4 приведены
распределения водонасыщенности по глубине около
границы x = 0 через 105 суток при m = 0.2; 0.4; 1 .
На рис. 5-6 приведены распределения водонасыщенности и температуры через 105.
Рис. 2. Распределение водонасыщенности
одномерной и двумерной задачи
21
ВЕСТНИК СВФУ, 2014, том 11, № 5
Рис. 3. Распределение водонасыщенности
через 105 суток на разных сетках
Рис. 5. Распределение водонасыщенности через 105 суток
Заключение
1. Вычислительные
расчеты
показали
эффекивность алгоритма.
2. Аппроксимация уравнения диффузии
в
дивергентном в недивергентном виде приводят к
одинаковым результатам.
22
Рис. 4. Распределение водонасыщенности через 105 суток
при пористости 0.2, 0.4, 1
Рис. 6. Распределение температуры
через 105 суток
T 0 = -22
Работа выполнена при финансовой поддержке
РФФИ р_восток_а - региональный конкурс «ДАЛЬНИЙ
ВОСТОК» № «12-01-98514».
В. В. Попов. ДВУМЕРНАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОМЕРЗАНИЯ ЗАСОЛЕННОГО ВЛАЖНОГО ГРУНТА
Литература
1. Ентов В. М., Максимов А. М., Цыпкин Г. Г.
Образование двухфазной зоны при промерзании пористой
среды. – М.: Препринт / ИПМех. АН СССР; № 269.
– 1986. – 55 с.
2. Максимов А. М., Цыпкин Г. Г. Образование двухфазной зоны при взаимодействии талых и мерзлых пород
с раствором соли. – М.: Препринт / ИПМех. АН СССР;
№ 305. – 1987. – 61 с.
3. Васильев В. И., Максимов А. М., Петров Е. Е.,
Цыпкин Г. Г. Математическая модель замерзания-таяния
засоленного мерзлого грунта // Прикладная механика и
техническая физика. – 1995. – Т. 36. – № 5. – С. 57-66.
4. Васильев В. И., Максимов А. М., Петров Е. Е., Цыпкин
Г. Г. Тепломассоперенос в промерзающих и протаивающих
грунтах. – М.: Наука, 1996.
5. Васильев В. И., Попов В. В. Численное решение
задачи промерзания грунта // Матем. моделирование. – 2008.
– Т. 20. – № 7. – С. 119-12.
6. Васильев В. И., Попов В. В. Численная реализация
математической модели взаимодействия талого грунта с
холодным раствором соли // Математические заметки ЯГУ.
– 2010. – Т. 17. – вып. 2. – С. 132-139.
7. Самарский А. А., Николаев Е. С. Методы решения
сеточных уравнений. – М.: Наука, 1978. – 562 с.
8. Корнеев В. Д. Параллельное программирование в
MPI. – Новосибирск.: ИВМ и МГ СО РАН, 2002. – 210 с.
References
1. Entov V. M., Maksimov A. M., Cypkin G. G.
Obrazovanie dvuhfaznoj zony pri promerzanii poristoj
sredy. – M.: Preprint / IPMeh. AN SSSR; № 269. – 1986.
– 55 s.
2.
Maksimov A. M., Cypkin G. G., Obrazovanie
dvuhfaznoj zony pri vzaimodejstvii talyh i merzlyh porod
s rastvorom soli. – M.: Preprint / IPMeh. AN SSSR;
№ 305. – 1987. – 61 s.
3. Vasil'ev V. I., Maksimov A. M., Petrov E. E.,
Cypkin G. G. Matematicheskaja model' zamerzanija-tajanija
zasolennogo merzlogo grunta // Prikladnaja mehanika i
tehnicheskaja fizika. – 1995. – Т. 36. – № 5. – S. 57-66.
4. Vasil'ev V. I., Maksimov A. M., Petrov E. E., Cypkin
G. G. Teplomassoperenos v promerzajushhih i protaivajushhih
gruntah. – M.: Nauka. 1996.
5. Vasil'ev V. I., Popov V. V., Chislennoe reshenie zadachi
promerzanija grunta // Matem. modelirovanie. – 2008. – T. 20. –
№ 7. – S. 119-12.
6. Vasil'ev V. I., Popov V. V. Chislennaja realizacija
matematicheskoj modeli vzaimodejstvija talogo grunta s
holodnym rastvorom soli // Matematicheskie zametki JaGU. –
2010. – T. 17. – vyp. 2. – S. 132-139.
7. Samarskij A. A., Nikolaev E. S. Metody reshenija
setochnyh uravnenij. – M.: Nauka, 1978. – 562 s.
8. Korneev V. D. Parallel'noe programmirovanie v MPI. –
Novosibirsk.: IVM i MG SO RAN, 2002. – 210 s.
23
Документ
Категория
Без категории
Просмотров
5
Размер файла
697 Кб
Теги
влажного, математические, промерзании, грунт, засоленного, двумерная, модель
1/--страниц
Пожаловаться на содержимое документа