Scientific journal
Advances in current natural sciences
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,775

Математическая модель приземного слоя атмосферы представляет собой замкнутую систему уравнений, записанных в безразмерном виде:

f(1)

с заданными краевыми условиями

f                  (2)

где U,V - проекции скорости единичной массы воздуха на горизонтальной оси координат Ox,Oy,

θ - потенциальная температура,

b2-величина, пропорциональная средней кинетической энергии турбулентности,

l - масштаб турбулентности,

K - коэффициент турбулентного обмена,

ε - средняя скорость диссипации,

R0 - число Россби,

αD = f - отношение коэффициента турбулентного обмена для примеси и коэффициента турбулентного обмена для импульса,

αН = f - отношение коэффициента турбулентного обмена для теплоты и коэффициента турбулентного обмена для импульса,

Χ - постоянная Кармана,

с - безразмерная константа,

r - безразмерный коэффициент пропорциональности, зависящий от вида диффундирующей примеси.

Как видно из системы (1), для определения первых двух неизвестных необходимо знать значение коэффициента турбулентного обмена K, который можно найти из уравнения 6.

Значения средней кинетической энергии и масштаба турбулентности найдем, решая совместно уравнения 4 и 5 системы (1), преобразовав пятое уравнение этой системы:

f             (3)

Учитывая, что l= l(z,t),

f

после дифференцирования первого уравнения системы (3) по z, получим:

f                          (4)

подставим значения b и f во второе уравнение системы (3), и преобразуем его к следующему виду:

f             (5)

Представим производные в виде конечных разностей:

f  (6)

откуда

f   (7)

 f            (8)

Получаем явную разностную схему для вычисления значений l в узлах сетки.

Для проведения вычислений разработана программа, с помощью которой было установлено, что масштаб турбулентности убывает с ростом высоты.

Вычислив значения l в узлах выбранной сетки, находим значения b в этих узлах, воспользовавшись первым равенством из (3):

f                 (9)

Из системы (1) видно, что зная l и b, легко рассчитать коэффициент турбулентного обмена K и среднюю скорость диссипации - ε. Подставляя значение K в уравнения 1-3 системы можно найти скорости движения воздушных масс и потенциальную температуру. Таким образом, все неизвестные параметры замкнутой системы (1-2) будут найдены.

Приведем примеры вычисления значений l и b, с помощью указанной программы, если время наблюдения - T=5 минут, нижняя граница пограничного слоя - Z1=0 м, высота пограничного слоя - Z2=1000 м, шаг сетки по высоте - dz=100 м, шаг сетки по времени - dt=30 с, начальное значение масштаба турбулентности - l0=500 м, граничное значение - lz1= f, где i-номер шага по времени. Рассчитанные значения представим в соответствующих таблицах.

Таблица значений масштаба турбулентности (l)

T,с

z, м

0

30

60

90

120

150

180

210

240

270

300

0

0,052

0,074

0,09

0,104

0,117

0,128

0,138

0,148

0,157

0,165

0,173

100

500

692,04

839,41

963,64

1073,09

1172,04

1263,03

1347,73

1427,27

1502,51

1574,07

200

500

678,08

814,72

929,92

1031,41

1123,16

1207,54

1286,07

1359,84

1429,60

1495,96

300

500

665,13

791,83

898,65

992,76

1077,84

1156,08

1228,90

1297,30

1361,99

1423,52

400

500

653,12

770,61

869,66

956,92

1035,82

1108,37

1175,89

1239,32

1299,30

1356,36

500

500

641,98

750,93

842,77

923,69

996,85

1064,12

1126,74

1185,55

1241,17

1294,08

600

500

631,66

732,68

817,85

892,88

960,71

1023,09

1081,16

1135,69

1187,27

1236,33

700

500

622,08

715,76

794,73

864,31

927,21

985,05

1038,89

1089,46

1137,29

1182,78

800

500

613,20

700,07

773,29

837,81

896,14

949,77

999,70

1046,59

1090,94

1133,12

900

500

604,97

685,52

753,42

813,24

867,33

917,06

963,36

1006,84

1047,96

1087,07

1000

500

597,34

672,02

734,99

790,46

840,61

886,73

929,66

969,98

1008,11

1044,38

Таблица значений b - величины, пропорциональной кинетической энергии турбулентности

t,с

z, м

0

30

60

90

120

150

180

210

240

270

0

0,00018

0,00014

0,00012

0,00010

0,00009

0,00009

0,00008

0,00008

0,00007

0,00007

100

1,631

1,251

1,055

0,929

0,840

0,773

0,719

0,675

0,639

0,608

200

1,512

1,160

0,978

0,862

0,779

0,716

0,667

0,626

0,592

0,563

300

1,402

1,076

0,907

0,799

0,722

0,664

0,618

0,581

0,549

0,522

400

1,300

0,998

0,841

0,741

0,670

0,616

0,573

0,538

0,509

0,484

500

1,205

0,925

0,780

0,687

0,621

0,571

0,532

0,499

0,472

0,449

600

1,118

0,858

0,723

0,637

0,576

0,530

0,493

0,463

0,438

0,417

700

1,037

0,795

0,671

0,591

0,534

0,491

0,457

0,429

0,406

0,386

800

0,961

0,738

0,622

0,548

0,495

0,455

0,424

0,398

0,377

0,358

900

0,891

0,684

0,577

0,508

0,459

0,422

0,393

0,369

0,349

0,332

1000

0,826

0,634

0,535

0,471

0,426

0,392

0,364

0,342

0,324

0,308

СПИСОК ЛИТЕРАТУРЫ

  1. Наац И.Э., Семенчин Е.А. «Математическое моделирование динамики пограничного слоя атмосферы в задачах мониторинга окружающей среды», Ставрополь, 1995 г.
  2. Турчак Л.И., Плотников П.В. «Основы численных методов», М, 2002 г.