КАРТОЧКА
ПРОЕКТА ФУНДАМЕНТАЛЬНЫХ И ПОИСКОВЫХ НАУЧНЫХ ИССЛЕДОВАНИЙ,
ПОДДЕРЖАННОГО РОССИЙСКИМ НАУЧНЫМ ФОНДОМ
Информация подготовлена на основании данных из Информационно-аналитической системы РНФ, содержательная часть представлена в авторской редакции. Все права принадлежат авторам, использование или перепечатка материалов допустима только с предварительного согласия авторов.
ОБЩИЕ СВЕДЕНИЯ
Номер 22-21-00271
НазваниеОптимизационный анализ обратных задач для дифференциальных уравнений тепломассопереноса, магнитной гидродинамики и электромагнетизма
РуководительАлексеев Геннадий Валентинович, Доктор физико-математических наук
Организация финансирования, регион Федеральное государственное бюджетное учреждение науки Институт прикладной математики Дальневосточного отделения Российской академии наук, Приморский край
Период выполнения при поддержке РНФ | 2022 г. - 2023 г. |
Конкурс№64 - Конкурс 2021 года «Проведение фундаментальных научных исследований и поисковых научных исследований малыми отдельными научными группами».
Область знания, основной код классификатора 01 - Математика, информатика и науки о системах, 01-111 - Дифференциальные уравнения с частными производными
Ключевые словадифференциальные уравнения, гидродинамика, тепломассоперенос, магнитная гидродинамика, обратные задачи, задачи управления, единственность, оценки устойчивости, численные алгоритмы
Код ГРНТИ27.31.00, 27.35.00
СтатусУспешно завершен
ИНФОРМАЦИЯ ИЗ ЗАЯВКИ
Аннотация
Настоящий проект направлен на развитие математического аппарата решения обратных задач для дифференциальных уравнений с частными производными, основанного на оптимизационном методе, и применение разрабатываемого аппарата для решения новых обратных задач для дифференциальных уравнений тепломассопереноса, магнитной гидродинамики и электромагнетизма. На основе разрабатываемой теории исполнителями будет проведен теоретический анализ рассматриваемых обратных задач и выполнено строгое математическое обоснование применения оптимизационного метода для их решения; установлены дополнительные свойства оптимальных решений в зависимости от выбора функционала качества и множества управлений, на котором он минимизируется; выведены достаточные условия на исходные данные, обеспечивающие существование, единственность и устойчивость оптимальных решений; разработаны эффективные численные алгоритмы решения рассматриваемых обратных экстремальных задач. Разработанный математический аппарат будет применен для решения прикладных задач, возникающих при проектировании специальных функциональных устройств, обеспечивающих невидимость материальных тел относительно волновых или статических физических полей.
Научная новизна настоящего проекта будет заключаться в получении фундаментальных результатов в области теоретического анализа новых обратных задач для дифференциальных уравнений тепломассопереноса, магнитной гидродинамики и электромагнетизма с использованием оптимизационного метода.
Предполагаемые результаты проекта будут иметь высокую научную значимость как в области теоретических исследований, так и в плане возможных практических приложений. Разрабатываемый математический аппарат представляет самостоятельную ценность, причем в силу своей фундаментальности и универсальности он может быть применен в будущем для исследования других новых обратных задач, не вошедших в рамки данного проекта. Создание основ данного аппарата потребует решения сложнейших математических проблем теории экстремальных задач для дифференциальных уравнений в частных производных.
Прикладная значимость результатов проекта связана с возможными приложениями развиваемой теории при решении широкого множества проблем от маскировки материальных тел до высокотехнологичной медицинской диагностики. Решения рассматриваемых обратных задач, полученные с помощью разрабатываемых алгоритмов, могут послужить основой для создания реальных технических устройств управления физическими полями различной природы, отличающихся простотой реализации и использованием распространенных природных материалов
Ожидаемые результаты
В результате выполнения проекта предполагается получить следующие научные результаты:
1) Развить оптимизационный метод исследования обратных задач для стационарных уравнений магнитной гидродинамики вязкой несжимаемой жидкости, рассматриваемых при смешанных краевых условиях для магнитного поля. На основе разработанного метода доказать новые теоремы о глобальной разрешимости обратных экстремальных задач, исследовать дополнительные свойства оптимальных решений и установить на основе их анализа достаточные условия на исходные данные, обеспечивающие локальную единственность и устойчивость оптимальных решений. Полученные результаты будут обобщены на случай более сложной модели магнитной гидродинамики вязкой теплопроводной жидкости, рассматриваемой при смешанных краевых условиях для магнитного поля и температуры и неоднородном граничном условии Дирихле для скорости.
2) Разработать математический аппарат исследования краевых и экстремальных задач для нелинейных моделей тепломассопереноса в случае, когда коэффициенты, входящие как в уравнения, так и в граничные условия рассматриваемых моделей, нелинейно зависят от искомого решения. На основе разработанного аппарата доказать новые теоремы о глобальной разрешимости краевых и экстремальных задач для рассматриваемых моделей, обобщающих приближение Обербека-Буссинеска, и исследовать устойчивость решений экстремальных задач.
3) Сформулировать и исследовать новые обратные экстремальные задачи для двумерных и трехмерных моделей переноса тепла и масс и электромагнетизма, возникающие при разработке технологий дизайна специальных функциональных устройств, служащих для управления соответствующими физическими полями в сплошных средах. К таким устройствам, наряду с устройствами, обеспечивающими невидимость материальных тел, можно отнести концентраторы, служащие для накопления энергии, вращательные оболочки, служащие для вращения волновых либо статических полей, устройства иллюзиона или камуфляжа и ряд других специальных устройств. На основе теоретического анализа указанных задач установить дополнительные свойства оптимальных решений в зависимости от выбора функционала качества и множества управлений с целью использования их при разработке численных алгоритмов решения рассматриваемых задач.
4) Развить эффективные численные алгоритмы решения двумерных и трехмерных обратных экстремальных задач для рассматриваемых моделей переноса тепла или масс и электромагнетизма, возникающих при разработке технологий дизайна специальных функциональных устройств. Исследовать свойства разработанных алгоритмов, выполнить вычислительные эксперименты и на основе их анализа оценить влияние основных параметров, входящих в рассматриваемые задачи, на свойства полученных решений. Применить полученные результаты для решения прикладных задач, возникающих при проектировании специальных функциональных устройств, и, в частности, высокоэффективных концентраторов энергии.
Оригинальность и научная значимость результатов, полученных при выполнении проекта, будет заключаться в том, что для решения исследуемых в проекте задач, являющихся с математической точки зрения коэффициентными обратными задачами для дифференциальных уравнений рассматриваемых моделей, будет применен и получит дальнейшее развитие подход, основанный на оптимизационном методе решения обратных задач. Это позволит, в частности, учесть те дополнительные условия и ограничения, которые связаны с необходимыми условиями реализуемости отыскиваемых решений рассматриваемых задач, и тем самым получить решения, допускающие относительно простую техническую реализацию. Кроме того, это позволит при исследовании рассматриваемых обратных задач использовать некоторые факты, такие как принцип максимума или bang-bang принцип, присущие именно решениям экстремальных задач, и получить тем самым дополнительную полезную информацию о свойствах искомых решений.
Результаты, полученные при выполнении проекта, будут иметь важное прикладное значение. В частности, результаты решения задач дизайна тепловых функциональных устройств могут использоваться при создании высокоэффективных концентраторов тепла нового типа, состоящих из доступных для технической реализации природных материалов.
Результаты, полученные при решении задач дизайна устройств магнитной невидимости, могут найти потенциальные применения в медицинских технологиях, в которых используемые магнитные поля должны иметь заданную структуру и не допускать возмущений различного рода присутствующими магнитными объектами. Примерами таких медицинских технологий являются магнитно-резонансная томография и транскраниальная магнитная стимуляция.
Разработанные исполнителями алгоритмы предполагается зарегистрировать в базе авторских свидетельств о государственной регистрации программ для ЭВМ.
ОТЧЁТНЫЕ МАТЕРИАЛЫ
Аннотация результатов, полученных в 2022 году
Разработан математический аппарат исследования новых краевых и экстремальных задач для нелинейных моделей массопереноса, обобщающих приближение Буссинеска. В основе доказательства глобальной разрешимости краевых задач с неоднородным условием Дирихле для концентрации вещества лежит применение принципа Лере-Шаудера и лифтингов для концентрации вещества с малой L4-нормой. Разработан универсальный метод, позволяющий установить принцип максимума и минимума для концентрации вещества, а также температуры среды для нелинейных уравнений тепломассопереноса в случае наиболее распространенных граничных условий для концентрации вещества и температуры.
С использованием данного математического аппарата, доказана глобальная разрешимость и локальная единственность решений краевых задач для нелинейных моделей массопереноса в случае, когда коэффициент реакции и коэффициент массового расширения нелинейно зависят от концентрации вещества. Установлены достаточные условия на исходные данные рассматриваемых краевых задач, при которых имеет место локальная единственность их решений, и условия, при которых справедлив принцип максимума для концентрации вещества.
Получил дальнейшее развитие математический аппарат исследования экстремальных задач для нелинейных моделей массопереноса. Применив этот аппарат, мы доказали разрешимость задач граничного управления для нелинейной модели массопереноса для функционалов качества общего вида, а также в случае, когда коэффициент реакции достаточно произвольно зависит от концентрации вещества. В свою очередь, для нелинейного уравнения реакции-диффузии-конвекции, в случае степенной зависимости коэффициентов реакции и массобмена, для экстремальных задач выведены системы оптимальности. На основе анализа данных систем получены оценки локальной устойчивости оптимальных решений для конкретных функционалов качества.
Разработана математическая технология решения задач дизайна устройств невидимости, экранирования и других специальных устройств, служащих для управления физическими полями, для стационарных моделей тепломассопереноса и электромагнетизма.
Разработанная технология была применена для решения двух рассматриваемых в проекте задач. Первой задачей является обратная задача для трехмерной модели теплопереноса, связанная с проектированием маскировочного устройства сферической формы. В рамках разработанной технологии доказана разрешимость прямой и обратных экстремальных задач, построена система оптимальности, описывающая необходимые условия экстремума, установлены некоторые свойства оптимальных решений, являющиеся следствием структуры системы оптимальности.
Второй задачей является задача дизайна многослойных цилиндрических радиально-анизотропных и изотропных экранирующих оболочек. С использованием оптимизационного метода упомянутая задача дизайна сведена к конечномерной экстремальной задаче, для решения которой разработан эффективный численный алгоритм, основанный на методе роя частиц. Вычислительные эксперименты показали, что предложенный метод позволяет спроектировать многослойную экранирующую оболочку, обладающую высокой эффективностью и простотой технической реализации.
Публикации
1. Алексеев Г.В. Theoretical analysis of cloacking problem for 3D model of heat conduction Дальневосточный математический журнал, Т. 22, N 2, С. 143-149 (год публикации - 2022) https://doi.org/10.47910/FEMJ202214
2. Сарицкая Ж.Ю. Краевая задача для нелинейных уравнений массопереноса с условием Дирихле Сибирские электронные математические известия, Т. 19, N 1. С. 360-370. (год публикации - 2022) https://doi.org/10.33048/semi.2022.19.031
3. Сарицкая Ж.Ю., Бризицкий Р.В. Boundary control problems for nonlinear reaction-diffusion-convection model Дальневосточный математический журнал, - (год публикации - 2023)
4. Спивак Ю.Э. Computer design of cylindrical cloaking shell for the magnetostatics model Дальневосточный математический журнал, Т. 22, N 2, С. 238-244 (год публикации - 2022) https://doi.org/10.47910/FEMJ202232
Аннотация результатов, полученных в 2023 году
Исполнителями сформулированы задачи управления для стационарных уравнений магнитной гидродинамики при смешанных граничных условиях для магнитного поля и неоднородном граничном условии Дирихле для скорости. Роль управлений играют нормальная и тангенциальная составляющие магнитного поля на разных частях границы и плотность сторонних электрических токов внутри области течения. Целью исследования является построение математического аппарата изучения локальной единственности и устойчивости оптимальных решений. Достижение этой цели основывается на обобщении подхода, развитого руководителем проекта для исследования устойчивости оптимальных решений для стационарных систем Навье--Стокса и магнитной гидродинамики (МГД) при стандартных краевых условиях для магнитного поля. В этом подходе не требуется определять вторую производную минимизируемого функционала качества, поскольку он основан на анализе фундаментальных свойств системы оптимальности для общей задачи управления для рассматриваемой модели и использовании специальных оценок разности решений исходной и возмущенной задач управления. Метод достаточно прост, естественен и применим к моделям гидродинамики, тепловой конвекции, массопереноса и другим гидродинамическим моделям, основанным на системе Навье-Стокса. Использование данного подхода позволяет получить оценки устойчивости оптимальных решений по отношению к малым возмущениям как функционала качества, так и одной из заданных функций, входящих в уравнение состояния, в явном и достаточно простом для интерпретации виде.
Исполнителями был выполнен теоретический анализ описанных выше задач управления для стационарных уравнений магнитной гидродинамики, рассматриваемых в ограниченной области трехмерного пространства, при смешанных краевых условиях для электромагнитного поля. Роль управлений в рассматриваемых задачах управления играют нормальная и тангенциальная компоненты магнитного поля, заданные на разных частях границы области течения, и плотность сторонних электрических токов, сосредоточенных внутри области течения. В качестве функционалов качества были выбраны квадратичные функционалы для скорости, магнитного поля или давления. Разработана общая теория (или математическая технология) исследования сформулированных задач управления, которая была применена для доказательства локальной единственности и устойчивости оптимальных решений рассматриваемых задач управления для конкретных функционалов качества.
Применение этой теории позволило доказать разрешимость задач управления для широкого класса слабо полунепрерывных снизу функционалов качества, вывести системы оптимальности и на основе их анализа доказать теоремы о локальной единственности и устойчивости оптимальных решений для квадратичных (tracking–type) функционалов относительно скорости, магнитного поля или давления.
Одной из особенностей проведенного исследования является доказательство разрешимости задач оптимального управления при минимальных требованиях к нормальной и тангенциальной компонентам магнитного поля, заданных на разных частях границы и играющих роль граничных управлений. Этот результат позволяет использовать простые L^2-нормы для нормальной и тангенциальной компонент магнитного поля вместо норм обычно используемых пространств следов на соответствующих частях границы в качестве регуляризаторов Тихонова, служащих для обеспечения корректности рассматриваемых некорректных задач управления.
Уникальность и универсальность разработанной технологии состоит в том, что она одинаково применима для решения как задач управления, так и задач идентификации параметров для рассматриваемой модели МГД, которые сводятся к задачам управления с помощью оптимизационного метода. Отметим, что наличие смешанных граничных условий для магнитного поля существенно усложняет исследование не только задач управления, но и исходной неоднородной краевой задачи для системы МГД. Ввиду этого, разработанная технология существенно использует авторскую находку, доказывающую возможность ослабления гладкости магнитной компоненты решения при негладких граничных данных для магнитного поля с сохранением фундаментального свойства единственности решения. Именно использование этой находки позволило доказать единственность и устойчивость решений рассматриваемых задач управления при минимальных требованиях на используемые граничные управления.
Исследованы задачи управления для нелинейного уравнения реакции-диффузии-конвекции, рассматриваемого при смешанных краевых условиях для концентрации вещества. Роль управлений в рассматриваемых задачах управления играют две функции, зависящие от пространственных переменных: коэффициент диффузии и множитель при коэффициенте распада реакции, зависящем от концентрации. В качестве функционалов качества выступают L^2 либо H^1-норма отклонения искомой концентрации от заданной концентрации, измеряемой в некоторой подобласти области, где изучается процесс распространения вещества. В случае степенной зависимости коэффициентов модели от решения для экстремальных задач выведены системы оптимальности. С их помощью получены оценки локальной устойчивости решений конкретных задач управления относительно малых возмущений как функционалов качества, так и одной из заданных функций краевой задачи.
В рамках направления 3) п 1.2 плана работ разработанная исполнителями в прошлом году математическая технология решения задач дизайна устройств невидимости была применена для решения задачи сопряжения и экстремальных задач для 3D модели магнитостатики. Указанные задачи связаны с дизайном трехмерных магнитных маскировочных оболочек. В качестве маскировочной среды был выбран эллиптический материал, заполняющий область, топологически эквивалентную сферическому слою. Доказана разрешимость краевой и экстремальных задач, построена система оптимальности, описывающая необходимые условия экстремума, установлены некоторые важные свойства оптимальных решений
Разработанная исполнителями математическая технология была применена для решения задач проектирования многослойных цилиндрических и сферических анизотропных и изотропных маскировочных устройств для двумерных моделей электропроводности и магнитостатики и трехмерной модели магнитостатики. С помощью упомянутой математической технологии и оптимизационного метода рассматриваемые задачи проектирования были сведены к конечномерным экстремальным задачам, для решения которых разработаны эффективные численные алгоритмы, основанные на методе роя частиц. Вычислительные эксперименты позволили сделать вывод об эффективности разработанной математической технологии и предложенных алгоритмов. С их помощью были установлены важные геометрические и физические свойства оптимальных решений, в частности, принцип релейности. Их применение дает возможность проектировать высокоэффективные маскировочные оболочки, обладающие простотой технической реализации для обеих рассматриваемых моделей.
Исследована задача идентификации неизвестного коэффициента диффузии для модели массопереноса по дополнительной информации о решении в случае, когда коэффициент диффузии является функцией, зависящей от пространственных переменных. В соответствии с общей теорией экстремальных задач построена система оптимальности, на основе которой разработан алгоритм численного решения указанной обратной экстремальной задачи. Проведено исследование разработанного численного алгоритма с помощью пакета Scilab, в котором решение прямой и сопряженной задач выполняется с помощью метода конечных разностей.
Для реализации и апробации численного алгоритма написана программа для ЭВМ "Восстановление коэффициента диффузии для стационарной модели массопереноса". Указанная программа зарегистрирована в Реестре программ для ЭВМ, и на нее получено свидетельство N 2023681974. Правообладатель О.В. Соболева.
Публикации
1. Алексеев Г.В. Теоретический анализ задач магнитной маскировки с использованием эллиптических метаматериалов Дальневосточный математический журнал, Том 23, Номер 2., С. 152-160 (год публикации - 2023) https://doi.org/10.47910/FEMJ202313
2. Алексеев Г.В. Analysis of control problems for stationary magnetohydrodynamics equations under the mixed boundary conditions for a magnetic field Mathematics, V. 11, no. 12: 2610, P. 1-29 (год публикации - 2023) https://doi.org/10.3390/math11122610
3. Бахтин К.Е., Савинов П.А. Оптимизационный метод решения задач маскировки для 2D модели электропроводности Дальневосточный математический журнал, Том 24. Номер 1. (год публикации - 2024)
4. Бризицкий Р.В., Максимов П.А. Об устойчивости решений задач управления для нелинейной модели реакции-диффузии-конвекции Дифференциальные уравнения, Tом 59, Номер 3, стр. 409-421 (год публикации - 2023) https://doi.org/10.31857/S0374064123030123
5. Спивак Ю.Э. Численная оптимизация в задачах проектирования многослойных магнитных маскировочных оболочек Дальневосточный математический журнал, Том 23, Номер 2, С. 264-269. (год публикации - 2023) https://doi.org/10.47910/FEMJ202323
6. Савинов П.А., Бахтин К.Е., Киличов У.Ш. Оптимизационный метод решения задач маскировки для 2d модели электропроводности Вычислительные технологии и прикладная математика : Материалы II Международного семинара (12-16 июня 2023 г., Благовещенск) / отв. ред. А. Г. Масловская. - Благовещенск : Изд.-во Амурского гос. университета, 2023., с. 169-171 (год публикации - 2023) https://doi.org/10.22250/9785934933921_169
7. Соболева О.В. Исследование численного алгоритма восстановления неизвестного коэффициента диффузии для стационарной модели переноса вещества Вычислительные технологии и прикладная математика : Материалы II Международного семинара (12-16 июня 2023 г., Благовещенск) / отв. ред. А. Г. Масловская. - Благовещенск : Изд.-во Амурского гос. университета, 2023., с. 184-186. (год публикации - 2023) https://doi.org/10.22250/9785934933921_184
8. Соболева О.В., Осипова М.А. Восстановление коэффициента диффузии для стационарной модели массопереноса -, 2023681974 (год публикации - )
Возможность практического использования результатов
Результаты, полученные в ходе выполнения проекта, возможно применять при решении ряда важных прикладных задач. В качестве некоторых из них можно отметить задачи маскировки материальных тел, задачи дизайна медицинских магнитночувствительных устройств, магнитно-резонансная томография (МРТ), транскраниальная магнитная стимуляция (ТМС), задачи проектирования устройств беспроводной передачи энергии и данных.