Труды КНЦ вып.38 (ГЕЛИОФИЗИКА вып. 4/2016(38))

дри dpu2 dpuw dt dx dz dpw dpuw dpw2 ~дГ dP dx д д -č (.-) - ^ dx VZdx,. 1 dx dz dP dPu dPw dP д „/ \ dw ------+— £ (z )— , dz dxz dx,. dz -P (V v )- d ( d dT ----- к - z ) ----- dx,. dx,. ( ( z ) dvk d^ d^,. dx,. - + y - 1Vd?t dx Q (z) z - ) ^ r To{z ) dz dz В уравнениях: p — плотность; u и w — горизонтальная и вертикальная компоненты скорости; Р — давление; g — гравитационная постоянная, у — коэффициент адиабаты; £(z); k(z) — коэффициенты вязкости и теплопроводности; T0(z) — фоновая температура. Предполагается отсутствие волн в начальный момент времени, поэтому плотность и температура берется при отсутствии волн, давление определяется состоянием идеального газа. Условия на горизонтальных границах области моделирования являются периодическими, тем самым имитируется большая область. Эмпирические данные модели MSIS-90 используются для задания начальной фоновой температуры и плотности. Гидродинамические уравнения решаются с помощью конечно-разностного метода с применением явных и неявных схем. Алгоритм численного интегрирования уравнений был описан в работах [8, 9]. Рис. 1. Экспериментальные данные вариаций атмосферного давления Ap(x0; t) 10-11 апреля 2006 г., наблюдаемые на инфразвуковой станции IS17 (6:70N; 4:90W) При расчетах были использованы оцифрованные экспериментальные данные давления на поверхности Земли (рис. 1) во время грозовых явлений 10-11 апреля 2006 г., которые были получены на инфразвуковой станции IS17 (6:70N, 4:90W). Эти данные используются при моделировании в качестве нижнего граничного условия. 164