CC BY
29
25 5 5 5 117
---8---2 +--1 +--12 +--3---0 = — = 058(3).
288 36 288 288 36 288 12
Результат совпадает с известным аналитическим решением. Заключение
Автор считает, что в данной работе новыми являются следующие положения и результаты. Новизна метода проявилась в том, что разработан алгоритм выбора узлов решетки при фиксированных коэффициентах, и эта задача является более трудной задачей. Для классов функций, зависящих от гладкостей по выбранным координатным осям, построенные формулы сохраняют полученный ранее порядок сходимости. Качество полученных методов проверены на контрольных задачах с известными решениями.
Литература
1. Соболев С.Л. Введение в теорию кубатурных формул. М.: Наука, 1974. 808 с.
2. Шойнжуров Ц.Б. Некоторые вопросы теории кубатурных формул в неизотропных пространствах С.Л. Соболева. // Докл. АН СССР, 1973. Т. 209/ № 5. С. 1036-1038.
3. Рамазанов М.Д. Лекции по теории приближенного интегрирования. Уфа: Изд-во Башкир. ун-та, 1973. 174 с.
4. Владимиров В.С. Уравнения математической физики. - 5-е изд., доп. М.: Наука, 1988.
512 с.
5. Никольский С.М. Приближение функций многих переменных и теоремы вложения. -2-е изд., перераб. и доп. М.: Наука, 1977. 456 с.
6. Шойнжуров Ц.Б. Некоторые вопросы теории кубатурных формул в пространстве
Wpm (En ) // Сб.Теория кубатурных формул и приложения функционального анализа к задачам
математической физики. Новосибирск: Наука, 1980. С.302-306.
7. Юмова Ц.Ж. Вычисление параметров функционалов погрешностей кубатурных формул. - Saarbrücken: LAP LAMBERT Academic Publishing GmbH&Co. KG, 2011. 108 c.
On the cubature formula with the involving derivatives
Yumova Tsyrenkhanda Zhembeevna, Candidate of Physics and Mathematics, Associate professor of the Department of High Mathematics,
East Siberian State University of Technology and Management, Ulan-Ude, Kluchevskaya st, 40V.
The algorithm of replacing some other lattice nodes to minimize the error functional norm is developed under the concordance condition of the convergence order with lattice step and the smoothness of the function along the selected coordinate directions. The results obtained by the method tested in the control examples.
УДК 532.517.4
ДИНАМИКА ЗАГЛУБЛЕНИЯ ПЕРЕМЕШАННОГО СЛОЯ ЖИДКОСТИ В ТРЕХМЕРНЫХ РАСЧЕТАХ
Лидия Алексеевна Компаниец, канд. физ.-мат. наук, доцент, с.н.с. тел.: 8 391 249 8811, e-mail: kla@icm.krasn.ru ФГБУН ИВМ СО РАН http://icm.krasn.ru Татьяна Валерьевна Якубайлик, канд. физ.-мат. наук, н.с. тел.: 8 391 249 8811, e-mail:ytv@icm.krasn.ru ФГБУН ИВМ СО РАН http://icm.krasn.ru
Людмила Владимировна Гаврилова, канд. физ.-мат. наук, доцент тел.: 8 391 249 7147, e-mail: lv.gavrilova@gmail.com ФГАОУ ВО "Сибирский федеральный университет" http://sfu-krasn.ru
С использованием трехмерной численной модели выполнены расчеты заглубления перемешанного слоя под действием ветровой нагрузки в стратифицированном водоеме прямоугольной формы. Проведено сравнение результатов для различных моделей турбулентности.
Ключевые слова: трехмерное численное моделирование, турбулентное замыкание, заглубление перемешанного слоя.
Введение
Задача численного исследования заглубления турбулентного слоя жидкости рассматривалась в ряде работ [1;2;3;4] и опирается на натурные измерения, приведенные в
работе [3]. Натурные эксперименты проводились для водоема глубины 0,3 ?'"к м. Область наших интересов — иссле-
дование основных гидрофизических характеристик неглубоких стратифицированных озер. В частности, нами ' применялся пакет GETM [5] для расчета течений в озере Шира, глубина которого 25 м. Поэтому в дальнейшем мы будем опираться на статью [2], в
Л.В. Гаврилова кот°р°й результаты статьи [3] те°ре- Т В Я(<убайлик тически (численно) обобщаются на
водоемы глубины 40 м, и результаты работы новосибирских ученых О. Ф. Васильева, Т. Э. Овчинниковой, Г. Г. Черных. Нашей задачей было сравнение результатов указанных авторов и результатов трехмерных расчетов по программе GETM.
Модель GETM — это наиболее динамично развивающаяся модель, специализированная для расчета течений в различных водоемах (в отличие, например, от универсальной MITgsm, которая позволяет, наряду с водоемами, рассчитывать и движение в атмосфере и за счет этого имеет некоторые ограничения в вариативности). В GETM возможен выбор из четырех вариантов сетки для дискретизации по вертикали, достаточно большой выбор разностных аппроксимаций для расчета адвекции (по сравнению, например, с POM). Также она позволяет подключать много разных моделей вертикального турбулентного перемешивания (для расчета коэффициента вертикального турбулентного обмена).
1 Постановка задачи
Постановка задачи по [4] следующая.
Уравнение состояния записывается в виде: Р = Pö (1 -а(( - T0 ) + ), где
—3 —3
а = 0,253547-10 1/°C и ß = 0,717031-10 л/г получены путем линеаризации уравнения состояния, предложенного Ченом и Миллеро [6], To = 3,98 °С — температура максимальной плотности для пресной воды.
Распределение температуры по вертикали постоянное (гомотермия) T = 3,5 °С. Имеем линейную стратификацию плотности по глубине. На рис. 1 приведены распределения солености и плотности по глубине в начальный момент времени.
В начальный момент времени скорости нулевые. На поверхности задается граничное условие, позволяющее учесть влияния ветра: К2 ди/& = тw|р, где тк —
напряжение ветра. Для солености на поверхности задаются нулевые потоки. На дне задаются условия проскальзывания с трением для скорости и постоянное значение для солености. Кроме того, задаются начальные данные и краевые условия для уравнений, позволяющих учесть вертикальный турбулентный обмен. На рис. 2 изображены значения относительной скорости, энергии турбулентности е, коэффициента вертикального турбулентного обмена и плотности. Рисунки предоставлены Г. Г. Черных и Т. Э. Овчинниковой.
z/H 1 -
0.8 —
0.6
0.4
0.2
0
0 0.5 U/Um„
Рис. 1. Начальное распределение солености и плотности
1
0.8 — 0.6 0.4 0.2 0
1
0.8 — 0.6 0.4 — 0.2 0
1
0.8 0.6 0.4 0.2 0
0.0000 0.0002 0.0004 e [m2/s2]
0 0.02 0.04 v [m2/s]
1000.2 1000.4 1000.6 р [kg/m3]
Рис. 2. Вертикальные распределения параметров для моментов времени I = 7,5 часов, t = 15 часов, t = 22,5 часа, t = 30 часов
30 -1
20 -
10 -
200
400
600
Рис. 3. Динамика заглубления пере-
u в
мешанного слоя hi = Nq hju зависимости от безразмерного времени t\ = Not
На рис. 3, взятом из [7], представлены экспериментальные данные [3] и результаты расчетов из статьи [4]. Точками обозначены результаты натурного эксперимента, сплошной линией — результаты численных расчетов.
12
Здесь Nq =
)t=o F
частота
Вяйсяля - Брента, h — толщина перемешанного
* / / * *
слоя, u = -утw/ р , р — начальное значение плотности жидкости на поверхности. Толщина перемешанного слоя определялась следующим образом: h = H - zh, где zh — величина, удовлетворяющая
условию e(t, zh ) = 0,02emax (t), «max (t) = max e(t, z).
0
0
No t
z
40 30 20 10
0
40 30 ^ 20 10 0
Нашей целью было исследовать различные модели турбулентности в расчетах задачи перемешанного верхнего слоя в модельном прямоугольном бассейне. При расчетах использовались пакеты программ ОЕТМ, ООТМ [8].
Система уравнений, используемая в этой программе, является модификацией хорошо известной и часто используемой для расчетов российскими и зарубежными учеными, системой уравнений гидрофизики, основанной на приближениях Буссинеска и гидростатики. По сравнению с системой уравнений, выписанной, например, в [9], из уравнений не исключается вертикальная молекулярная вязкость.
Коэффициент вертикального турбулентного обмена рассчитывается с использованием моделей турбулентного замыкания разного уровня сложности. Существует возможность расчета с постоянным коэффициентом вертикального турбулентного обмена. Запрограммировано несколько вариантов к -в модели. Также присутствует возможность расчета по моделям второго порядка. Коэффициенты горизонтальной турбулентной диффузии берутся постоянными.
Далее приведены результаты расчетов по двум моделям турбулентности.
1. Вариант турбулентного замыкания, который в ООТМ именуется как модель с вычислением ТКЕ и масштаба длины (далее модель 1). При вычислении ТКЕ используется динамическое уравнение типа Меллора-Ямады [1]. Величина I вычисляется через динамическое уравнение скорости диссипации. Функция стабильности постоянна.
2. Модель 2-го порядка с динамическим уравнением типа к -в модели (далее модель 2).
2 Тестовые расчеты
Тестовые расчеты проводились на трехмерной задаче о заглублении верхнего перемешанного слоя в прямоугольном бассейне шириной 7,5 км и длиной 10 км с постоянной глубиной 40 м.
Моделирование проводилось для стратифицированной по солености жидкости с использованием ОЕТМ. В этом случае скорости зависят от времени, глубины и еще двух пространственных переменных. Наш опыт вычислений показывает, что если в модели ОЕТМ коэффициент вертикального турбулентного обмена не слишком большой, то влияние бокового обмена на величины скоростей в центре бассейна невелико и можно попытаться сравнить одномерный модельный расчет [4] со специальной усовершенствованной моделью турбулентности [1].
На рис. 4 представлена динамика заглубления перемешанного слоя И = ЩИ/и * в зависимости от безразмерного времени ^ = Not. Графики выведены в центральной точке
Центр прямоугольного бассейна на экваторе
Центр прямоугольного бассейна на 54° северной широты
200
400
600
800
1000
1200
Рис. 4. Динамика заглубления перемешанного слоя
бассейна.
Толщина перемешанного слоя определялась следующим образом:
a) модель 1, К2 (г, zh )< 10
b) И = 1513 и *Щ-2^1/3
И = Н - 1и , где 1и — величина, удовлетворяющая условию
ч-4 2/
м /с,
с) И = (б/5)14и*щЛ1^
ё) к = 214 и * N-1¥/2
е) модель 2, к = Н - 1к , где 1к —
К2 (t, )< 10 -4 м2/с.
Соотношения Ь) - ё) взяты из статьи [2].
На рис. 5 приведен типичный график коэффициента вертикального турбулентного обмена в центре бассейна, который качественно совпадает с графиком этой величины, приведенном в [7].
Результаты, приведенные на рис. 4-5 получены при коэффициенте горизонтального обмена Кх = 10 м2/с.
Кроме того, исследовалось влияние коэффициента горизонтального турбулентного обмена на величину заглубления турбулентного слоя. На рис. 6 представлены результаты расчетов с Кх = 200 м /с.
величина, удовлетворяющая условию
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016
К^, м/с
Рис. 5. Распределение коэффициента вертикального турбулентного обмена по глубине
Центр прямоугольного бассейна на экваторе
20
200 400 600 800 1000 1200
Центр прямоугольного бассейна на 54° северной широты
20
Рис. 6. Динамика заглубления перемешанного слоя, Кх = 200
Заключение
Физика задачи о заглублении перемешанного верхнего слоя в замкнутом прямоугольном бассейне с жесткими стенками при постоянных граничных условиях обуславливает выход на квазистационарное течение, что демонстрируют рис. 4 и 6. Аналитические формулы Ь) - ё) со степенной зависимостью величины заглубления перемешанного верхнего слоя от времени выхода на стационар не предусматривают.
Мы пришли к следующим выводам:
1) получаемый результат зависит от того, каким взят параметр Кориолиса (т.е. от места расположения бассейна);
2) в начальный момент времени в реальных расчетах перемешивание идет быстрее, чем в расчетах по аналитическим формулам;
3) с течением времени в реальных расчетах перемешивание практически прекращается, что может быть связано с выходом решения на стационар.
4) увеличение коэффициента горизонтального турбулентного обмена приводит к ускорению процесса перемешивания.
Авторы благодарят Г. Г. Черных и Т. Э. Овчинникову за постановку задачи и плодотворное обсуждение результатов.
25
20
15
10
5
40
0
0
40
0
0
200
800
200
N t
Литература
1. Mellor G.L., Yamada T. Development of a turbulence closure model for geophysical fluid problems // Reviews of geophysics and space physics. 1982. V. 20. No. 4. Pp. 851-875.
2. Deleersnijder E., Luyten P. On the practical advantages of the quasi-equilibrium version of the Mellor and Yamada level 2.5 turbulence closure applied to marine modelling // Appl. Math. Modelling. 1994. V.18. Pp.281-287.
3. Kato H., Phillips O.M. On the penetration of a turbulent layer into stratified fluid // J. Fluid. Mech. 1969. V. 37, part 4. Pp. 643-655.
4. Васильев О.Ф., Овчинникова Т.Э., Черных Г.Г. Математическое моделирование заглубления турбулентного слоя в стратифицированной жидкости // Доклады академии наук. 2012. Т. 443. № 5. С. 578-582.
5. GETM [Электронный ресурс]. URL: http://getm.eu (дата обращения: 22.02.2016).
6. Chen C.T., Millero F.J. Precise thermodynamic properties for natural waters covering only the limnologies range // Limnol. Oceanogr. 1986. Vol.31(3). Pp. 657-662.
7. Васильев О.Ф, Овчинникова Т.Э., Черных Г.Г. Численные модели заглубления турбулентного слоя в устойчиво стратифицированной жидкости // Математическое моделирование. 2015. Т. 27. № 5. C. 52-64.
8. GOTM [Электронный ресурс]. URL: http://www.gotm.net/ (дата обращения: 22.02.2016).
9. Марчук Г.И., Саркисян А.С. Математическое моделирование циркуляции океана. М.: Наука, 1988. 304 c.
Dynamics of fluid mixed layer penetration in 3d calculations
Lidiya Kompaniets, PhD, Associate Professor, Senior Researcher, ICM SB RAS
Tatiana V. Yakubaylik, PhD, Researcher, ICM SB RAS
Lyudmila Gavrilova, PhD, Associate Professor, Siberian Federal University
The wind-influenced deepening of the mixed layer in a stratified rectangular reservoir has been calculated using three-dimensional numerical model. A comparison of results has been made for various turbulence models.
Keywords: three-dimensional numerical simulation, turbulent closure, mixed layer penetration.
УДК 519.714
ВЕРИФИКАЦИЯ ЛОГИЧЕСКИХ ОПИСАНИЙ С ФУНКЦИОНАЛЬНОЙ НЕОПРЕДЕЛЕННОСТЬЮ НА ОСНОВЕ МОДЕЛИРОВАНИЯ
Людмила Дмитриевна Черемисинова, д.т.н., главный научный сотрудник Тел.: 375 17 284 2 82, e-mail: cld@newman.bas-net.by Объединенный институт проблем информатики НАН Беларуси
http://uiip.bas-net.by
Рассматривается подход к верификации логических описаний комбинационных устройств для случая, когда исходное функциональное описание проектируемой схемы не полностью определено. Предлагаемый подход основан на моделировании поведения схемы для общего случая, когда входные стимулы представляются троичными векторами.
Ключевые слова: автоматизация проектирования, верификация, моделирование.
Одним из наиболее эффективных методов повышения качества и надежности оборудования является создание средств автоматизации, гарантирующих корректность проектируемых устройств. Это может быть обеспечено путем верификации проектных решений на всех стадиях проектирования. По мере возрастания сложности проектируемых устройств функциональная верификация становится все более необходимым и дорогим этапом процесса проектирования. Официальные издания проектных компаний
