Возврат на первую страницу веб-сайта ВЦ

in English

ОТДЕЛ МЕХАНИКИ СПЛОШНЫХ СРЕД

Заведующий отделом — д.ф.м.н. Чарахчьян Александр Агасиевич

Сектор аэродинамики, зав. сектором — д.ф.м.н. Чарахчьян Александр Агасиевич
      (E-mail: chara@ccas.ru)
Сектор моделирования климата, зав. сектором — к.ф.м.н. Пархоменко Валерий Павлович
      (E-mail: parhom@ccas.ru)
Сектор асимптотических методов, зав. сектором — д.ф.м.н. Жук Владимир Иосифович
      (E-mail: zhuk@ccas.ru)


Содержание страницы



Историческая справка

      Отдел механики сплошных сред возник одним из первых в Вычислительном центре Российской академии наук. С момента образования и до 1991 года его возглавлял крупный ученый в области аэрогидродинамики профессор Ю.Д.Шмыглевский. С 1992 по 2002 год отдел возглавлял известный специалист в области механики сплошных сред Ученый секретарь института д.ф.м.н. Л.В.Шуршалов.

      В отделе разрабатываются математические модели, вычислительные методы и алгоритмы для исследования физических процессов в сплошных средах.



Ведущие научные сотрудники отдела

д.ф.м.н. Ю.Д.Шмыглевский
д.ф.м.н. Л.В.Шуршалов
д.ф.м.н. В.М.Борисов
д.ф.м.н. В.И.Жук
д.ф.м.н. В.И.Зубов
д.ф.м.н. И.Е.Михайлов
д.ф.м.н. А.М.Тарко
      к.ф.м.н. А.В.Шипилин (1938-2003)
к.ф.м.н. Б.Н.Азаренок
к.ф.м.н. В.Н.Котеров
к.ф.м.н. В.М.Кривцов
к.ф.м.н. В.А.Люлька (1937-2006)
к.ф.м.н. В.П.Пархоменко
к.ф.м.н. И.В.Савенков
к.ф.м.н. А.А.Фролова

      Совместители: д.ф.м.н. В.Н.Диесперов (МФТИ), к.ф.м.н. А.А.Иванков (НПО им. С.А.Лавочкина).



Тематика отдела



Основные научные результаты за 1998-2003 гг.

  1. Уравнения Навье-Стокса для установившихся плоскопараллельных течений несжимаемой жидкости преобразованы к однородным формам. Изучены бесконечные вихревые системы в вязкой жидкости в окрестности угловых точек границы области течения. Показана возможность любых подходов разделяющих линий тока к обтекаемому контуру в плоскопараллельном потоке вязкой жидкости. Изучены стоксовские плоскопараллельные течения в каналах с отрывом потока от их стенок. Построено аналитическое представление некоторых плоскопараллельных отрывных течений вязкой несжимаемой жидкости у границ областей течения. Найден класс точных решений уравнений Навье-Стокса в осесимметричном случае. Класс содержит произвольную функцию и определяет, в частности, системы вихревых образований типа шарового вихря Хилла.

  2. Разработан программный комплекс «Теплозащита». С его помощью проведены расчеты теплозащиты космических аппаратов при их движении в атмосфере Венеры.

  3. Проведено численное моделирование газодинамических явлений, которые могут развиваться в протяженных подводных газопроводах как при штатных технологических мероприятиях (например, пуск газопровода), так и при нештатных ситуациях (разрыв или значительное повреждение трубопровода). Исследование проведено на основе специально разработанной математической модели, учитывающей неидеальность перекачиваемого газа (метан+высшие углеводороды+азот+углекислый газ) и большие перепады высот залегания трубопровода. Изучен вопрос возможного обнаружения места разрыва трубопровода по изменениям газодинамических характеристик потока на контрольных участках газопровода. Работа выполнена по двум хоз. договорам ВЦ РАН № 17 от 2 февраля 1998 г. и № 18 от 5 января 2000 г. с ВНИИГАЗ (головная научная организация ГАЗПРОМ) в рамках проекта сооружения магистрального газопровода Россия-Турция через акваторию Черного моря (проект «Голубой поток»).

  4. Разработана математическая модель и создан комплекс программ, предназначенный для расчета внутренних пространственных течений газа в многоступенчатых охлаждаемых осевых турбинах. Моделирование осуществляется путем численного интегрирования осредненных по Рейнольдсу трехмерных уравнений Навье-Стокса, дополненных полуэмпирической двухпараметрической qw-моделью турбулентности.

  5. Разработана математическая модель, предложен вычислительный алгоритм и создан комплекс программ для проведения численных исследований влияния параметров трехмерного эжекторного сопла на его тяговые и иные характеристики при наличии внешнего обтекания данного устройства. Работа выполнена в рамках хоз. договора ВЦ РАН с ОАО «ОКБ Сухого» № 18/2002 от 11 июля 2002 г.

  6. В рамках методологии быстрого автоматического дифференцирования получены новые формулы вычисления матрицы вторых производных сложных функций и продемонстрирована их эффективность на модельных задачах оптимального управления. Исследована задача оптимального управления процессом плавления и кристаллизации вещества по нахождению распределения мощности источника по времени, обеспечивающего минимум суммарной энергии при заданных ограничениях на массу расплавленного вещества и скорость кристаллизации. Ограничения, накладываемые на решение задачи, учитывались с помощью метода штрафных функций. Для вычисления градиента целевого функционала в дискретном случае использовалась методология быстрого автоматического дифференцирования. Выяснено влияние на оптимальное управление таких параметров задачи, как предельная скорость остывания, ширина источника, и размер области у оси симметрии, где допускается произвольная скорость остывания. Предложен новый численный алгоритм решения прямой задачи (определения температурного поля), который позволяет существенно ускорить получение решения.

  7. Разработан алгоритм и создана программа для многопроцессорного компьютера для расчета двумерных и трехмерных течений разреженного газа на базе нелинейного уравнения Больцмана.

  8. Численно исследовано сжатие углерода в конических свинцовых мишенях с помощью алюминиевых ударников, имеющих скорость 4 км/с. Постановка задачи включает в себя уравнения гидродинамики сжимаемых сред, полуэмпирические широкодиапазонные уравнения состояния материалов и кинетическую модель неравновесного превращения графита в алмаз, откалиброванную с учетом известных экспериментальных данных. Исследованы все стадии процесса от кумулятивного эффекта на оси симметрии мишени до разгрузки образца, вызванной приходом волны разрежения от тыльной поверхности ударника, учитывая явление откола.

  9. Численно исследована неустойчивость границы раздела сред после выхода на нее двух ударных волн, следующих одна за другой через некоторый промежуток времени. Рассмотрены свободная граница алюминия, граница алюминий — дейтериевый лед, воздух — гелий и фреон — воздух. Показано, что если к моменту выхода на границу второй ударной волны развитие возмущений находится на нелинейной стадии, то, в отличие от линейной стадии, изменение скорости роста амплитуды слабо зависит от длины волны начального возмущения. Предложена инженерная формула, позволяющая определять воздействие второй ударной волны на скорость роста амплитуды, и показана ее хорошая точность как для синусоидального начального возмущения, так и для негладкого возмущения с точками излома. Продемонстрирована возможность определения параметров второй ударной волны, обеспечивающих эффект замораживания неустойчивости, по скорости роста амплитуды независимо от длины волны начального возмущения.

  10. Численно исследовано бесструйное сжатие пластины на клине сильной ударной волной. Построено приближенное решение с присоединенной ударной волной, учитывающее фазовый состав материала пластины в волне разгрузки, и показана хорошая точность этого решения. Обнаружено значительное (в несколько раз) увеличение температуры за фронтом присоединенной ударной волны по сравнению с температурой за фронтом сжимающей волны. Показана принципиальная возможность инициирования термоядерной реакции при бесструйном сжатии пластины из дейтериевого льда сильной ударной волной.

  11. Исследована задача о нагреве сжимаемой жидкости постоянным потоком тепла. Обнаружено, что начиная с некоторого момента времени давление всюду за ударной волной мало отличается от некоторого постоянного значения. Получены приближенные аналитические формулы, из которых следует независимость давления от величины коэффициента теплопроводности и некоторые закономерности зависимости давления от величины потока тепла. На ряде примеров продемонстрировано хорошее соответствие между формулами и решением задачи.

  12. Численно исследованы течения пористых сред со сходящимися ударными волнами при наличии некоторого механизма ограничения кумуляции: ударное сжатие в конических твердотельных мишенях, одномерные сферические и цилиндрические волны с ограничением кумуляции в виде жесткого шарика или цилиндра небольшого радиуса. Обнаружено, что переход от сплошного вещества к пористому в этих течениях приводит к существенному увеличению не только температуры, но и давления. Такой характер изменения давления при изменении пористости качественно отличается от случая плоских волн постоянной интенсивности, для которых увеличение пористости всегда приводит к уменьшению давления.

  13. Развит метод построения регулярных адаптивных конечно-разностных сеток в применении к расчету неодномерных газодинамических течений с разрывами. Развит новый метод построения регулярных сеток для течений с сильно искривленными подвижными границами. Построен алгоритм решения задачи о распаде разрыва с произвольным уравнением состояния в табличной форме; предложен способ коррекции разностной схемы типа С.К.Годунова, направленный на подавление энтропийного следа. Разработан общий способ вычисления первых вариаций функционалов разной кратности. Разработан метод псевдохарактеристик для расчета многомерных уравнений гиперболического типа. Разработан алгоритм трехмерной консервативной интерполяции на гексаэдральных сетках.

  14. Разработаны варианты моделей циркуляции атмосферы, гидро-термодинамики океана и морского льда в арктическом регионе, эффективно аккумулирующие, адаптирующие и использующие существующие и новые данные о процессах взаимодействия атмосфера — океан — лед, конвективных и динамических процессах в океане, а также данные о структуре и эволюции морских льдов. Разработан вариант компьютерной системы для анализа, обработки и визуализации результатов численных экспериментов и данных наблюдений. Осуществлена численная реализация и создан программный комплекс для статистического анализа данных с помощью обобщенных линейных разделяющих функций.

  15. Получены многолетние (на сроки до 70 лет) распределения характеристик морского льда в виде карт, графиков как для средних по пространству и времени величин, так и в отдельных точках. Представлены также карты расчетных распределений толщины льда в Арктике для различных сезонов в установившемся режиме. Установлены новые результаты, связанные с сильной естественной межгодичной изменчивостью ледового покрова. Проведен анализ чувствительности модельных результатов к параметрам моделей, определяющих взаимодействие атмосферы, льда и океана. Установлен существенно неоднородный по пространству характер влияния на ледовый покров этих параметров. Так, уменьшение альбедо снега и льда приводит к уменьшению в среднем толщины ледового покрова, но существуют области, где толщина увеличивается.

  16. Проведено изучение проблемы распространения нефтяных разливов в океане при наличии морского льда для построения математических моделей. Предложены соответствующие модели.

  17. Разработана система математических моделей для прогнозирования воздействия хозяйственной деятельности на природные системы и климат. Система моделей биосферных содержит точечные и пространственные модели биогеохимических циклов углерода и азота в системе «атмосфера — растения — почва» и «атмосфера — океан». Пространственное разрешение 0.5х0.5 и 4х5 град. георафической сетки. Получены долгосрочные прогнозы изменения количества двуокиси углерода в атмосфере и климата, а также экологических процессов в результате действия индустриальных выбросов двуокиси углерода в атмосферу, вырубки лесов, эрозии почв. Получены расчеты баланса поглощения и выделения двуокиси углерода стран и регионов при различных сценариях экономической деятельности. Сделаны расчеты последствий ограничения выбросов СО2 в атмосферу — выполнения Киотского протокола к Рамочной конвенции ООН об изменении климата и его невыполнения, в частности, США. Разработана математическая трактовка выполнения принципа Ле-Шателье в биосфере. Определена способность экосистем суши и океана компенсировать антропогенные воздействия на биосферу. В пространственной модели системы атмосфера-океан рассчитана способность различных океанов и их частей поглощать выбросы СО2, а также оценена роль биоты океана поглощении выбросов и формировании пространственной и временной картины динамики полей СО2 и неорганического азота.

  18. Разработана система моделей действия атмосферных загрязнений на лесные экосистемы в регионе.

  19. Изучена тонкая структура областей взаимодействия пограничных слоев с внешним потенциальным потоком (локальный и глобальный отрыв, сверхзвуковые зоны с образованием ударных волн, генерация волн Толлмина-Шлихтинга, распространение возмущений солитонного типа). Проанализированы различные аспекты неустойчивости пограничных слоев и пристеночных струй. Названные исследования нацелены на решение исключительно актуальной общей проблемы ламинарно-турбулентного перехода (вызывающей интерес на протяжении более столетия), а также математических свойств уравнений Навье-Стокса при больших числах Рейнольдса.

  20. Подытожены исследования различных аспектов асимптотической теории неклассического пограничного слоя (монография В. И. Жука «Волны Толлмина-Шлихтинга и солитоны», М.: Наука, 2001 г.).



Основные публикации 1998-2003 гг.



Участие в проектах и хоз. договорных работах

  1. Проект РФФИ № 96-01-00869 «Методы и алгоритмы решения пространственных оптимизационных задач информатики и математической физики»
  2. Проект Миннауки РФ по подпрограмме «Перспективные информационные технологии» № 0201.06.233 «Информационные технологии, опирающиеся на использование высокопроизводительных вычислительных систем для задач гидрогазодинамики в фундаментальных исследованиях и технических приложениях»
  3. Проект Миннауки РФ № 201-23(00)-П, ФЦТНП «Исследования и разработки по приоритетным направлениям развития науки и техники гражданского назначения», подпрограмма «Перспективные информационные технологии»
  4. ФЦП «Мировой океан». Проект № 5.14/6-01 «Моделирование циркуляции атмосферы, океана и эволюции морского льда в арктическом регионе с целью изучения процессов взаимодействия атмосфера — океан — лед, конвективных и динамических процессов в океане», проект 5.14 «Исследование физических процессов, определяющих основные черты гидрологического режима арктических морей России»
  5. ФЦНТП «Исследования и разработки по приоритетным направлениям развития науки и техники». № 12/02, Тема № 9: «Математическое моделирование промышленных аварий и экологических чрезвычайных ситуаций» (шифр темы 01.37.03.01.09), Блок 1 «Ориентированные фундаментальные исследования», Раздел «Информационные технологии и электроника»
  6. Проект РФФИ № 96-01-00193 Математические модели и методы планетарной экологии
  7. Проект РФФИ № 99-01-00945 Математические модели и методы для исследования проблемы глобального потепления в применении к проблемам Арктики и Субарктики
  8. Проект РФФИ № 02-01-01040 Разработка математических моделей и методов для исследования и прогнозирования глобальных и региональных климатических и экологических последствий антропогенной деятельности в применении к проблемам Арктики и Субарктики
  9. Проект РФФИ № 97-01-00095 Математическое моделирование импульсного сжатия веществ в конических твердотельных мишенях
  10. Проект РФФИ № 00-01-00551 Численное исследование гидродинамических течений, вызванных ударным сжатием конденсированных сред
  11. Проект ВНШ № 96-15-96298 Математическое моделирование глобальных и региональных процессов (методология, реализация)
  12. Проект ВНШ № 00-15-96119 Математическое моделирование глобальных и региональных процессов (методология, реализация)
  13. Программа № 17 Президиума РАН «Параллельные вычисления на многопроцессорных вычислительных системах», тема Параллельные методы и алгоритмы решения больших задач оптимизации и численного моделирования. 2003 г.
  14. Программа № 3 фундаментальных исследований ОМН РАН «Вычислительные и информационные проблемы решения больших задач» (ГК № 10002-251/ОМН-03/026-024/240603-805 от 24.06.2003 г.)
  15. РФФИ 96-01-01928 «Нелинейное взаимодействие сильных разрывов и обтекание препятствий в проводящей среде магнитным полем применительно к космической плазме»
  16. РФФИ 99-01-00264 «Конструирование алгоритмов построения адаптивных сеток на основе теории гармонических отображений»
  17. РФФИ 02-01-00236 «Разработка современных технологий построения многомерных сеток и их приложения в задачах математической физики»
  18. Проект Института математики и механики УрО РАН, Екатеринбург, «Разработка алгоритмов и программ построения регулярных трехмерных сеток и интерполяция газодинамических полей»
  19. Проект РФФИ № 98-01-00517 «Теория и приложения метода быстрого автоматического дифференцирования»
  20. Проект РФФИ «Ведущие научные школы» № НШ 00-15-96080
  21. Проект РФФИ «Ведущие научные школы» № НШ-1737.2003.1
  22. Хоз. договор № 17 от 2 февраля 1998 г. с ВНИИГАЗ «Прогноз нестационарных технологических режимов работы газопровода и процессов в газопроводе при его разрыве»
  23. Хоз. договор № 18 от 5 января 2000 г. с ВНИИГАЗ «Разработка основных положений методики определения мест разрыва подводного газопровода по изменениям параметров технологического режима»
  24. Хоз. договор № 18/2002 от 11 июля 2002 г. с ОАО «ОКБ Сухого» «Разработка методики численных расчетов параметров течения внутри схематизированного плоского эжекторного сопла при наличии внешнего обтекания. Проведение предварительных параметрических расчетов и сравнение с экспериментальными данными»
  25. РФФИ № 98-01-00341: «Асимптотические модели вязких нестационарных трансзвуковых течений» (1998-2000 гг.).
  26. РФФИ № 01-01-00999: «Исследование различных механизмов передачи возмущений вверх и вниз по потоку в задачах теории свободного взаимодействия пограничного слоя при до-, сверх- и трансзвуковых скоростях движения» (2001-2003 гг.)
  27. Проект «Учебно-научный центр по фундаментальной и прикладной математике» (на базе МФТИ, МИ РАН, ВЦ РАН) целевой программы «Интеграция»


Участие в конференциях




Страница размещена на сайте 19 октября 2005 года.
Модифицирована 10 января 2007 года.
Предшествующий вариант страницы