close

Вход

Забыли?

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

?

Применение метода галерейна к расчету вертикального распространения примеси в нестационарном устойчивом пограничном слое атмосферы.

код для вставкиСкачать
Пусть функция Z i x , y , i
условиям:
Zi ! 0 в : , Zi
0, n , удовлетворяет
0 на w: i , ’Zi
1 на w: i ;
тогда функция
Z
Z0 љ D Z1 љ D Z2 љ D ... љ D Zn
Здесь f
ECf 0 ? продолжение в область :
функции f 0 ( EC ? оператор продолжения граничm 1
ных условий); ) 0
, ) i , i 1, n , ? неопределенные компоненты структуры;
wZ w
wZ w
wx wx wy wy .
D Z1
(18)
будет удовлетворять таким условиям:
Z ! 0 в : , Z 0 на w: , ’Z
1 на w: .
В соответствии с методом R-функций можно построить структуры решения краевых задач (14), (15)
и (16), (17).
Структура решения краевой задачи (14), (15) имеет
вид
u 0m1
Z
fZ01
n
¦ Zk1
Z2)0m1
~ Z2) m1
M
0
0
а структура решения краевой задачи (16), (17) имеет
вид
n
n
Zi1 / ¦ Z k1 ZD Z1 (Zi1 / ¦ Z k1 ) Z 2 ) i
k 0
k 0
~ Z2 ) .
M
i
i
(20)
УДК 519.673+519.3
ПРИМЕНЕНИЕ МЕТОДА
ГАЛЕРКИНА К РАСЧЕТУ
ВЕРТИКАЛЬНОГО
РАСПРОСТРАНЕНИЯ ПРИМЕСИ В
НЕСТАЦИОНАРНОМ
УСТОЙЧИВОМ ПОГРАНИЧНОМ
СЛОЕ АТМОСФЕРЫ
ВЫХОДЦЕВ Е.И.
Рассматривается задача расчета вертикального распространения примеси от приземного источника в
нестационарном устойчивом пограничном слое атмосферы. Для приближенного решения этой задачи
предлагается использовать метод Галеркина. Предлагаемая вычислительная схема тестируется на модельной задаче, для которой известно точное решение.
Введение
Развитая структура и значительные объемы производства, транспорта и хранения низкотемпературных веществ представляет повышенную опасность
для окружающей среды и человека. Существующие
методики расчета последствий аварийных ситуаций, связанных с разливами, выбросами и утечками
низкотемпературных веществ в окружающую среду, часто дают сильно завышенные результаты, что
приводит к дополнительным затратам при составлении плана ликвидации последствий аварии.
Поэтому разработка методик, которые давали бы
44
Литература: 1. Лойцянский Л.Г. Механика жидкости и
газа. М.: Наука, 1987. 840 с. 2. Pearson C.E. A computational
method for viscous flow problems // J. Fluid Mech. 1965.
21, №4. P. 611-622. 3. Вабищевич П.Н. Метод фиктивных
областей в математической физике. М.: Изд-во МГУ,
1991. 156 с. 4. Рвачев В.Л. Теория R-функций и некоторые ее приложения. К.: Наук. думка, 1982. 552 с.
,(19)
k 0
ui
Для аппроксимации неопределенных компонент в
(19), (20) можно воспользоваться любым проекционным методом.
Поступила в редколлегию 19.11.2002
Рецензент: д-р физ.-мат. наук, проф. Колосов А.И.
Сидоров Максим Викторович, аспирант каф. прикладной математики ХНУРЭ. Научные интересы: математическое моделирование, математическая физика,
теория R-функций и ее приложения. Увлечения и
хобби: история культуры. Адрес: Украина, 61166, Харьков, пр. Ленина, 14, тел. (0572) 702-14-36.
более точные расчеты, а также позволяли проводить
прогностические расчеты для более полного класса
внешних условий, является актуальной.
1. Постановка задачи
Вертикальное распространение оседающей примеси от мгновенного источника мощности Q , находящегося на высоте h , описывается полуэмпирическим уравнением турбулентной диффузии [1]
wq
wt
wq ·
wq
w §
Ёkz
ёw
,
wz ©
wz №
wz
(1)
где q ? средняя концентрация примеси; w ?
скорость гравитационного оседания частиц; k z ?
вертикальный коэффициент турбулентной диффузии. Уравнение (1) дополняется следующими
начальными и краевыми условиями:
qt 0
Q G zh ,
(2)
wq
wq
w E q
0 , kz
0 , (3)
wz
wz z H
z z0
здесь E ? скорость сухого осаждения (характериkz
стика взаимодействия частиц примеси с подстилающей поверхностью); z 0 ? параметр шероховатости; H H t ? высота верхней границы нестационарного устойчивого пограничного слоя атмосферы; h ? высота источника примеси.
Наиболее полный обзор известных точных решений задачи приведен в [1].
РИ, 2003, № 1
Для неоседающей и неосаждающейся примеси
( w E 0 ), а также для неоседающей, но осаждающейся примеси ( w 0 , E z 0 ), при постоянных
вертикальном коэффициенте турбулентной диффузии k z и высоте верхней границы устойчивого
пограничного слоя атмосферы H точное решение
задачи (1)-(3) может быть легко получено методом
разделения переменных. Для первого случая оно
имеет вид ( z 0 0 )
f Q 2Q
¦e
H H j1
q z, t
k z S 2 j2
H2
t
cos
Sjh
Sjz
cos
, (4)
H
H
Под обобщенным решением задачи (1)-(3) будем
понимать функцию q z, t , которая почти при всех
t Џ 0, T принадлежит H A , имеет производную
wq
Џ L 2 0, T u : и при произвольной v z Џ H A
wt
почти всюду на 0, T удовлетворяет соотношениям
§ wq
·
Ё , v ё t q, v t Kq, v t
© wt
№
q z, 0 , v
Q¦
q z, t
)j h )j z
)j z
j 1
где Tj t
e
)j z
k z O jt
,) j z
2
cos O j z ,
(5)
E
sin O j z ,
kz
H kz E
k E
z
sin 2H O j 2k z
4k z O j
E
E
cos 2H O j 4k z O j
4k z O j ,
2
k
1 §
Ё Sj arctg z
2 Ё
E
H ©
Oj
Подставив (8) в (6) и (7) для определения коэффициентов c i t , приходим к системе ОДУ вида
N
dc j
j 1
dt
¦
Ku { w
du
dz
с областями определения
где V
Ѕ°
dq
0ѕ,
dz z H °
ї
v : v Џ W21 : , Av Џ L 2 : ; D K ‹ D A .
Через H A обозначим соответствующее оператору
A энергетическое пространство. Скалярное произведение и норма в H A имеют соответственно вид
u, v
H
і kz
z0
u
H
du dv
dz E w uv z z ,
0
dz dz
2
§ du ·
2
і k z Ё ё dz E w u
z z0 .
© dz №
z0
РИ, 2003, № 1
0 (9)
QM k h , k 1, N . (10)
Можно доказать, что задача (9)-(10) однозначно
разрешима.
В соответствии с общей методикой построения
структур решения краевых задач [3, 4] базисные
функции возьмем в виде
·
ёW z ё j
№
H z z z 0 2z z 0 H dW j z
dz ,
H z0 2
Mj z
z 0 , H функция. Введем в рассмотрение опе-
0, kz
j 1
j 1
раторы
­°
dv
DA ®v : v Џ V, kz w E v
dz
°?
z z0
N
M j , M k ¦ c j M j , M k KM j , M k
с начальными условиями, получаемыми из СЛАУ:
Пусть k z ? положительная и ограниченная на
d §
du ·
ё,
Ёkz
dz ©
dz №
(8)
j 1
N
Аналитические решения (1)-(3) получены лишь для
узкого класса задач, что делает актуальной проблему разработки эффективных численных методов.
Воспользуемся для приближенного аналитического решения задачи (1)-(3) методом Галеркина [2].
Au { N
¦ c j 0 M j z , Mk
2. Метод решения
:
(7)
где M j z ? полная система линейно-независимых
функций из H A .
2
·
ёё .
№
Qv h .
¦cj t Mj z ,
q N z, t
Tj t
(6)
Приближенное решение задачи (1)-(3) будем искать в виде
а для второго случая ? вид
f
0,
2
§
Ё1 E w H z z z 0
Ё
kz H z0 2
©
где W j z ? полная система функций (полиномы
Лежандра, сплайны и др.).
3. Численные результаты
Численное моделирование было проведено в три
этапа. Во всех случаях источник располагался на
таких уровнях: h
OH , где O
0,
1 1 3
, , ,1;
4 2 4
H 150 ? . Скорость сухого осаждения варьирова3
лась в пределах от 5 � 10
до 5 � 10 2 м/с. Для
оседающей неосаждающейся примеси полагалось
3
, 2 . Решение искалось в виде (8),
2
в качестве системы W j z были выбраны кубические
E
Jw , где J 1,
сплайны. Компьютерное моделирование выполнялось с помощью пакета Mathematica 4.0©.
45
3.1. Неоседающая и неосаждающаяся примесь
( w E 0)
4
Профили концентрации примеси для этого случая
представлены на рис. 1. Как видно, при постоянном
(усредненном по высоте) коэффициенте диффузии
с течением времени примесь равномерно распределяется по всей толщине пограничного слоя, что
соответствует теоретическим результатам (см. выражение (4)). Значения, отличающиеся от предель-
3
Q
ного значения
менее чем на 10 % , устанавливаH
ются примерно начиная с момента времени
W ??? 50 ??? . Было также установлено, что равномерное распределение наиболее быстро устанавливается при расположении источника в середине
пограничного слоя.
3.2. Неоседающая осаждающаяся примесь
( w 0, E z 0)
Профили концентрации примеси для этого случая
представлены на рис. 2. При учете скорости сухого
осаждения появляется осадок на земле. Часть примеси, находящейся в пограничном слое, с ростом
времени диффузии уменьшается и с увеличением
E быстро осаждается на землю.
В случаях 3.1 и 3.2 полученное приближенное
решение сравнивалось с точным решением вида (4)
и (5) соответственно. На сетке сплайнов с шагом 10
отклонение от точного решения в фиксированные
моменты времени по норме L 2 z 0 , H составило
менее 15%.
3.3. Оседающая осаждающаяся примесь
(w z 0 , E z 0)
Профили концентрации примеси для этого случая
представлены на рис. 3. Учет гравитационного
оседания приводит к сосредоточению примеси у
поверхности земли. Анализ численного решения
показывает, что сухое осаждение приводит к сглаживанию профиля концентрации и уменьшает
время установления линейного профиля практически по всей толщине пограничного слоя.
4
1
20
40
60
80
100
120
140
Рис. 2. Эволюция вертикального профиля
концентрации ( w 0, E 0,01 , h 0,25H )
4
3
2
1
20
40
60
80
100
120
140
Рис. 3. Эволюция вертикального профиля
концентрации ( w E 0,01 , h 0,25H )
Выводы
Таким образом, предложен и обоснован приближенный метод расчета вертикального распространения примеси в нестационарном устойчивом пограничном слое атмосферы. Этот метод был протестирован на модельной задаче, для которой известно точное решение. Также получено приближенное решение задачи, для которой точное решение не
известно. Полученные результаты хорошо согласуются с теоретическими и экспериментальными
данными [1].
Литература: 1. Берлянд М.Е. Прогноз и регулирование
загрязнения атмосферы. Л.: Гидрометеоиздат, 1985.
277 с. 2. Михлин С.Г. Численная реализация вариационных методов. М.: Наука, 1966. 432 с. 3. Рвачев В.Л.
Теория R-функций и некоторые ее приложения. К.:
Наук. думка, 1982. 552 с. 4. Выходцев Е.И. Построение
структур решения краевых задач для уравнения турбулентной диффузии // Радиоэлектроника и информатика. 2002. № 3. С. 55-57.
3
2
Поступила в редколлегию 29.10.2002
1
Рецензент: д-р физ.-мат. наук, проф. Дикарев В.А.
20
40
60
80
100
120
Рис. 1. Эволюция вертикального профиля
концентрации ( w E 0 , h 0,25H )
46
2
140
Выходцев Евгений Иванович, научный сотрудник кафедры ПМ ХНУРЭ. Научные интересы: математическая физика, теория R-функций, компьютерное моделирование. Адрес: Украина, 61058, Харьков, ул. Данилевского, 8, кв.130, тел. 43-87-84.
РИ, 2003, № 1
1/--страниц
Пожаловаться на содержимое документа