КАРТОЧКА ПРОЕКТА ФУНДАМЕНТАЛЬНЫХ И ПОИСКОВЫХ НАУЧНЫХ ИССЛЕДОВАНИЙ,
ПОДДЕРЖАННОГО РОССИЙСКИМ НАУЧНЫМ ФОНДОМ

Информация подготовлена на основании данных из Информационно-аналитической системы РНФ, содержательная часть представлена в авторской редакции. Все права принадлежат авторам, использование или перепечатка материалов допустима только с предварительного согласия авторов.

 

ОБЩИЕ СВЕДЕНИЯ


Номер 21-71-00006

НазваниеПрямое численное моделирование волновой турбулентности в нелинейных средах со слабой дисперсией

РуководительКочурин Евгений Александрович, Кандидат физико-математических наук

Организация финансирования, регион Федеральное государственное бюджетное учреждение науки Институт электрофизики Уральского отделения Российской академии наук, Свердловская обл

Период выполнения при поддержке РНФ 07.2021 - 06.2023 

Конкурс№60 - Конкурс 2021 года «Проведение инициативных исследований молодыми учеными» Президентской программы исследовательских проектов, реализуемых ведущими учеными, в том числе молодыми учеными.

Область знания, основной код классификатора 01 - Математика, информатика и науки о системах, 01-309 - Турбулентность и гидродинамическая устойчивость

Ключевые словаПрямое численное моделирование, волновая турбулентность, нелинейные волны, волновые взаимодействия, прямое численное моделирование, спектры Колмогорова-Захарова, волновые резонансы

Код ГРНТИ29.27.21


СтатусУспешно завершен


 

ИНФОРМАЦИЯ ИЗ ЗАЯВКИ


Аннотация
Создание новых и эффективных инструментов прямого численного моделирования турбулентного движения является серьезным вызовом перед современной вычислительной физикой. Несмотря на взрывной рост производительности вычислительной техники, наблюдавшийся в последние несколько десятилетий, проблема численного описания турбулентности до сих пор актуальна. Это связано с тем, что для корректного статистического описания турбулентного движения необходимо проводить расчеты с высоким пространственным разрешением и в течение длительного интервала времени. Задача еще сильнее усложняется при описании сложного трехмерного движения сплошной среды, когда увеличивается число независимых переменных, требуемых для полного описания эволюции системы. В этом случае высокую эффективность демонстрируют расчетные технологии, основанные на параллельных вычислениях. В рамках настоящего проекта мы планируем использовать самые современные технологии распределенных вычислений на основе графических процессоров для описания волновой турбулентности, возникающей в нелинейных средах со слабой дисперсией. Волновая (или слабая) турбулентность наблюдается при взаимодействии нелинейных дисперсионных волн во многих физических процессах [1]. Примерами такой турбулентности являются оптическая турбулентность [2], турбулентность дисперсионных капиллярных и гравитационных волн на поверхности жидкости [3,4], акустическая турбулентность [5], МГД-турбулентность альфвеновских волн [6,7]. Важно, что слабая турбулентность возникает в условиях, когда в исследуемой системе есть малый параметр (амплитуда волн). Это позволяет применять подход, основанный на слабонелинейном разложении исходных уравнений движения. Фундаментальная физическая теория, описывающая статистический характер взаимодействия нелинейных волн, носит название теории слабой турбулентности. Она впервые предложена в конце 1960-х годов в работах ак. Захарова и соавт. [3-5]. Важнейшим достижением теории слабой турбулентности является то, что она аналитически предсказывает показатели спектров турбулентности в произвольных нелинейных волновых системах (т.н. спектры Колмогорова-Захарова). Этот факт отличает спектры слабой турбулентности от классического спектра Колмогорова-Обухова, полученного на основе качественного размерностного анализа для описания сильной вихревой турбулентности. К настоящему моменту времени спектры Колмогорова-Захарова с высокой точностью подтверждены численно и экспериментально для дисперсионных нелинейных волн на поверхности жидкости [8-10]. Отдельной задачей в теории слабой турбулентности является исследование волн, распространяющихся без дисперсии. К таким волновым процессам можно отнести акустические волны [11], в том числе волны второго звука в жидком гелии [12], магнитогидродинамические волны [13], и волны на поверхности жидкости во внешнем магнитном (электрическом) поле [14]. Главной проблемой в описании турбулентности бездисперсионных волн является тенденция к их опрокидыванию, приводящая к формированию сильно нелинейных ударных фронтов. В этом случае альтернативой теории слабой турбулентности является подход, предложенный в [15, 16], согласно которому акустическая турбулентность рассматривается как ансамбль газодинамических ударных волн. Таким образом, в настоящий момент времени существуют два различных подхода к описанию турбулентности нелинейных акустических бездисперсионных волн. Настоящий научный проект направлен на решение фундаментальной задачи о реализации спектров Колмогорова-Захарова в средах со слабой дисперсией. На основе методов прямого численного моделирования будут исследованы две физически различные системы. Первая задача посвящена моделированию трехмерного движения жидкости со свободной границей во внешнем магнитном (электрическом) поле. В работах [17,18] было экспериментально обнаружено, что на поверхности магнитной жидкости может развиваться особый тип турбулентности слабодисперсионных магнитогидродинамических волн. Такой тип турбулентности в плоско-симметричном случае был исследован теоретически в наших недавних работах [19,20]. В настоящем проекте мы впервые исследуем новый тип волновой турбулентности в рамках полной 3D геометрии. Вторая задача будет посвящена прямому численному моделированию акустической турбулентности в рамках трехмерного модельного уравнения [21]. Отметим, что в настоящий момент времени не существует научных работ, которые достоверно воспроизводят спектр Колмогорова-Захарова для бездисперсионных сред в рамках методов прямого численного моделирования. Таким образом, в рамках настоящего проекта будет впервые систематически исследована волновая турбулентность, возникающая в нелинейных средах со слабой дисперсией. Разработанные численные модели будут основаны на самых современных методах параллельных вычислений. Главным результатом работы будет являться непосредственное обнаружение спектров Колмогорова-Захарова для новых, ранее не исследованных слабодисперсионных волновых систем. 1. Zakharov V.E., Falkovitch, G. and V. S. L’vov, Kolmogorov Spectra of Turbulence I: Wave Turbulence (Springer, Berlin, 1992). 2. Picozzi A., et al, Physics Reports, v. 542, p. 1, (2014). 3. Zakharov V.E., Filonenko N.N. J. Appl. Mech. Tech. Phys. V. 4, p. 506, (1967). 4. Zakharov V. E. J. Appl. Mech. Tech. Phys. 9, 190 (1968). 5. Zakharov V. E., Sagdeev R. Z. Dokl. Akad. Nauk SSSR, V. 192, N. 2, p. 297–300 (1970). 6. Ирошников П.С., Астрономический журнал, т. 40, с. 742, (1963). 7. Kraichnan R.H., Physics of Fluids, v. 8, p. 1385, (1965). 8. Pushkarev A.N., Zakharov V.E. Phys. Rev. Lett. 76, 3320 (1996). 9. Kolmakov G.V., Brazhnikov M.Y., Levchenko A.A., Abdurakhimov L.V., McClintock P.V.E., and Mezhov-Deglin L.P. Prog. Low Temp. Phys. 16, 305 (2009). 10. Falcon E., Discrete Contin. Dyn. Syst. B 13, 819 (2010). 11. L’vov, V. S., L’vov, Y., Newell, A. C., & Zakharov, V. (1997). Phys. Rev. E, 56(1), 390. 12. Ganshin, A. N., Efimov, V. B., Kolmakov, G. V., Mezhov-Deglin, L. P., & McClintock, P. V. Phys. Rev. Letters, 101(6), 065303 (2008). 13. Andrés, N., Clark Di Leoni, P., Mininni, P. D., Dmitruk, P., Sahraoui, F., & Matthaeus, W. H. Physics of Plasmas, 24(10), 102314 (2017). 14. Zubarev, N.M. Phys. Lett. A. 333 (3-4), 284–288 (2004). 15. Кадомццев Б.Б., Петвиашвили В.И., ДАН СССР, в. 208, в. 4, с. 794, (1973). 16. Kuznetsov, E. A. JETP Lett., 80(2), 83-89 (2004). 17. Boyer, F. and E. Falcon. Phys. Rev. Lett. 101, 244502 (2008). 18. Dorbolo, S. and E. Falcon. Phys. Rev. E 83, 046303 (2011). 19. Kochurin, E.A. JETP Lett. 109 (5), 303-308 (2019). 20. Kochurin, E.A. J. Magn. Magn. Mater. 503, 166607 (2020). 21. Zakharov, V.E. Appl. Mech. Tech. Phys., 6(4), 22-24 (1965).

Ожидаемые результаты
Общий план работ по проекту будет разделен на два этапа. Первый год будет посвящен прямому численному моделированию волновой магнитогидродинамической турбулентности поверхности жидкости в полной трехмерной геометрии. По итогам первого этапа выполнения проекта будут получены следующие научные результаты: 1. Будет сформулирована математическая модель, описывающая нелинейную динамику свободной поверхности жидкости, помещенной во внешнее тангенциальное электрическое (магнитное) поле. Полученная система будет включать действие сил поверхностного натяжения, вязкости и энергетической накачки. В пределе сильного внешнего поля распространение волн на поверхности жидкости будет происходить без дисперсии. 2. Будет разработан программный комплекс для прямого численного моделирования волновой магнитогидродинамической турбулентности свободной поверхности жидкости в полной трехмерной геометрии. В основе численного алгоритма будут лежать псевдо-спектральные методы расчета пространственных производных и интегральных операторов. Для ускорения расчетов будет использована технология CUDA распределенных вычислений на графических процессах. 3. На основе разработанного программного комплекса будет впервые проведено прямое численное моделирование волновой турбулентности поверхности жидкости во внешнем магнитном поле. Будет показано, что под действием внешней вынуждающей силы (энергетической накачки) исследуемая система переходит в режим квазистационарной диссипации энергии. Спектр поверхностных возмущений в этом режиме приобретает степенной характер с показателем близким к аналитическому спектру Колмогорова-Захарова, полученному для бездисперсионных волн. 4. На основе данных прямого численного моделирования будут детально исследованы статистические свойства магнитогидродинамической турбулентности поверхности жидкости. В частности, будут построены функции плотности вероятности для амплитуды поверхности жидкости. Будет анализировано нелинейное дисперсионное соотношение, а также проведен корреляционный анализ, демонстрирующий ключевую роль трех-волновых резонансов в процессе перекачки энергии в малые масштабы. На втором этапе выполнения проекта будет проведено прямое численное моделирование слабой турбулентности нелинейных волн акустического типа. В основе модели будет лежать трехмерное уравнение Бюргерса с диссипацией и слабой дисперсией [Захаров, ПМТФ, (1965)]. Основные научные результаты, которые планируется получить в ходе выполнения второго этапа проекта: 1. Будет разработана вычислительная модель трехмерной волновой турбулентности слабодисперсионных волн, в основе которой будет лежать численное решение волнового уравнения с квадратичной нелинейностью. Численное решение модельного уравнения будет строиться на основе псевдо-спектральных методов и с использованием новейших средств параллельных вычислений CUDA. Такие методы позволят проводить вычисления с высокой точностью и скоростью. 2. Будет впервые проведено прямое численное моделирование волновой турбулентности слабодисперсионных волн в полной трехмерной геометрии. Ранее такие исследования проводились в рамках упрощенных уравнений, описывающих лишь эволюцию спектра волн, т.е. не позволяющих восстановить профиль волнового поля. Прямое численное моделирование, запланированное в рамках настоящего проекта, позволит непосредственно исследовать волновую турбулентность и ее статистические свойства. 3. Будет показано, что в режиме развитой волновой турбулентности волновой спектр приобретает степенной характер в инерционной области волновых чисел. Впервые будет проведено сравнение показателей численного спектра и аналитического спектра Колмогорова-Захарова (Захарова-Сагдеева), полученного для бездисперсионных звуковых волн. Несмотря на то, что акустическая турбулентность исследовалась многими авторами, в рамках прямого численного моделирования спектры слабой турбулентности для бездисперсионных систем до сих пор не были воспроизведены убедительно. Таким образом, настоящий научный проект предполагает систематические исследование турбулентности слабодисперсионных волн в различных физических системах. Главным результатом проекта будет непосредственное обнаружение спектра Колмогорова-Захарова, полученного для нелинейных волн, распространяющихся без дисперсии. Обнаружение аналитических спектров слабой турбулентности для новых физических систем – сложная и амбициозная задача, важная как с прикладной, так и с фундаментальной точки зрения. Решение поставленной задачи существенно расширит современное представление о волновой турбулентности, возникающей в средах со слабой дисперсией. Мы убеждены, что предполагаемые научные результаты находятся на передовом крае современной мировой науки.


 

ОТЧЁТНЫЕ МАТЕРИАЛЫ


Аннотация результатов, полученных в 2021 году
За отчетный период получены следующие научные результаты: 1. Сформулирована математическая модель, описывающая развитие турбулентности слабодисперсионных электро- или магнитогидродинамических волн на свободной поверхности жидкости во внешнем электрическом (магнитном) поле. Модель учитывает квадратично нелинейные эффекты и включает в рассмотрение силы поверхностного натяжения, вязкости, механической накачки и пондеромоторные силы, действующие со стороны внешнего поля. 2. Разработан программный комплекс для прямого численного моделирования волновой магнитогидродинамической турбулентности свободной поверхности жидкости в трехмерной геометрии. В основе программы лежат алгоритмы быстрого преобразования Фурье для расчета пространственных производных и интегральных операторов. Интегрирование по времени проводится явным методом Рунге-Кутта четвертого порядка точности. Для ускорения точности и времени расчета используется технология параллельных вычислений CUDA на графических процессорах NVidia. На разработанный программный комплекс получено свидетельство о государственной регистрации. 3. Проведено детальное численное моделирование волновой электрогидродинамической турбулентности свободной поверхности жидкости во внешнем электрическом поле, направленном тангенциально невозмущенной поверхности. Показано, что под воздействием внешней вынуждающей силы система переходит в режим квазистационарной диссипации энергии, когда действие внешней силы полностью компенсируется влиянием вязкости. 4. Исследованы статистические свойства поверхностной магнитогидродинамической волновой турбулентности. В квазистационарном состоянии построены функции плотности вероятности. Найденные распределения плотности вероятности близки к гауссовому распределению. Исследовано нелинейное дисперсионное соотношение в режиме развитой волновой турбулентности. Как и ожидалось волны, распространяющиеся вдоль внешнего поля, являются слабодисперсионными. Для возмущений, распространяющихся перпендикулярно полю, обнаруживается вторая ветвь дисперсионного соотношения, представляющая собой нелинейно связанный волновой пакет (т.н. связанные волны). Впервые проведен корреляционный анализ, демонстрирующий ключевую роль трех-волновых резонансных взаимодействии в развитии поверхностной МГД турбулентности. 5. Впервые исследована развитая волновая турбулентность поверхности жидкости в режиме сильного внешнего электрического (магнитного) поля. Показано, что в этом режиме движение жидкости сильно анизотропно: доминируют возмущения, распространяющиеся перпендикулярно внешнему полю. Показатель спектра турбулентности становится близок к значению −3, что согласуется с оценками на основе качественного размерностного анализа спектров слабой турбулентности. Отметим, что подобный спектр справедлив для классической МГД турбулентности альфвеновских волн в проводящей жидкости или плазме.

 

Публикации

1. Е. Кочурин, Г. Рикард, Н. Зубарев, Е. Фалькон Three-Dimensional Direct Numerical Simulation of Free-Surface Magnetohydrodynamic Wave Turbulence Physical Review E, - (год публикации - 2022)

2. Е.А. Кочурин,М.А. Гашков,Н.М. Зубарев Formation of a Small-Scale Chaotic Relief on the Free Surface of a Dielectric Liquid in a Tangential Electric Field Proceedings of IEEE Conference on Dielectric Liquids, - (год публикации - 2022)

3. Кочурин Е.А. Программный комплекс для трехмерного численного моделирования волновой магнитогидродинамической турбулентности свободной поверхности жидкости -, №2022616151, дата государственной регистрации 05 апреля 2022 года (год публикации - )


Аннотация результатов, полученных в 2022 году
За отчетный период получены следующие научные результаты: 1. Проведено детальное численное моделирование электрогидродинамической (ЭГД) турбулентности поверхности жидкости для широкого диапазона управляющих параметров: напряженности поля и уровня механической накачки системы. Нам удалось впервые провести моделирование перехода дисперсионной капиллярной турбулентности к бездисперсионной ЭГД турбулентности при увеличении внешнего поля. В режиме развитой ЭГД турбулентности полная электрическая энергия намного превышает энергию поверхностных капиллярных волн, что свидетельствует о доминирующем влиянии электрогидродинамических эффектов в эволюции системы. 2. Детально исследованы свойства анизотропии поверхностной электро- и магнитогидродинамической волновой турбулентности. Построены точные двумерные спектры турбулентности относительно фазового угла в Фурье пространстве. Оказалось, под влиянием случайного (некоррелированного) внешнего воздействия на границе жидкости может развиваться квази-изотропная МГД турбулентность. Измеренные функции плотности вероятности для амплитуды границы жидкости оказались близки к нормальному гауссовому распределению. Этот факт свидетельствует о развитом характере исследуемой турбулентности. 3. Проведены серии численных экспериментов для фиксированного поля. Было впервые показано, что коэффициенты спектров имеют корневую зависимость от скорости диссипации энергии. Этот факт однозначно свидетельствует о доминирующем влиянии трех-волновых резонансных взаимодействий в эволюции рассматриваемой турбулентности. Отметим, что гипотеза о доминирующем влиянии резонансных трех-волновых взаимодействий лежала в основе подхода к поиску аналитического спектра поверхностной МГД турбулентности. Численные данные с высокой точностью согласуются с полученным аналитическим спектром турбулентности. Получена численная оценка для постоянной Колмогорова-Захарова для поверхностной электрогидродинамической и магнитогидродинамической волновой турбулентности. 4. Исследованы процессы передачи магнитной (электрической) энергии в энергию капиллярных волн. Показано, что в первую очередь происходить накачка гармоник, ответственных за электрогидродинамические эффекты. После накачки электрической энергии происходит ее перенос в область капиллярных волн, распространяющихся перпендикулярно внешнему полю. Оценены характерные времена нелинейной перекачки энергии. Они оказались намного больше линейного времени (периода линейной волны). Этот факт подтверждает слабонелинейный характер эволюции системы. 5. Разработан программный пакет для прямого численного моделирования акустической турбулентности в трехмерной геометрии. В основе численной модели лежат высокоточные спектральные методы расчета пространственных производных. Для ускорения расчетов программа позволяет использовать технологию распределенных вычислений CUDA на графических процессорах NVidia.

 

Публикации

1. Дмитриев И.А., Кочурин Е.А.,Зубарев Н.М. Numerical study of free-surface electro-hydrodynamic wave turbulence IEEE Transactions on Dielectrics and Electrical Insulation, - (год публикации - 2023) https://doi.org/10.1109/TDEI.2023.3256350

2. Кочурин Е.А. Программный комплекс для моделирования трехмерной акустической турбулентности Свидетельство о государственной регистрации программы для ЭВМ, №2023616015 от 21.03.2023 (год публикации - 2023)

3. - Ученые смоделировали движение жидкостей, подобных космической плазме Indicator, ФИЗИКА 13:18, 29 ИЮНЯ 2022 (год публикации - )


Возможность практического использования результатов
Дальнейшее направление использования результатов проекта видится в возможности применения разработанных методов численного моделирования движения жидкостей. Создание новых и эффективных методов прямого численного моделирования течения жидкостей – важная задача для перехода к передовым цифровым технологиям в различных областях высокотехнологического производства. Сложное гидродинамическое течение возникает при плавлении материалов во многих физических процессах (к примеру, при горении вакуумной дуги, лазерной абляции вещества, плазменной и электронной обработке поверхностей и т.д.). Управление микроструктурами, возникающими в результате плавления материала, принципиально важно для повышения качества 3D печати. Разработанный в рамках проекта программный комплекс нацелен на создание новых методов управления хаотическим мелкомасштабным рельефом поверхности материалов, находящихся в жидкой фазе. Его главной задачей является обнаружение и исследование нового механизма турбулентного движения жидкости, который может быть использован, к примеру, для создания новых материалов с заданным микрорельефом поверхности и уровнем шероховатости. Таким образом, настоящий проект направлен на решение фундаментальных и прикладных задач, соответствующих направлению Н1 из Стратегии НТР РФ «Переход к передовым цифровым, интеллектуальным производственным технологиям, роботизированным системам, новым материалам и способам конструирования, создание систем обработки больших объемов данных, машинного обучения и искусственного интеллекта».