РФЯЦВНИИТФ Эксперимент и математическое моделирование в механике сплошных сред Всероссийская конференция «Взрыв в физическом эксперименте» 16 – 20 сентября. - презентация

1

2 РФЯЦВНИИТФ Эксперимент и математическое моделирование в механике сплошных сред Всероссийская конференция «Взрыв в физическом эксперименте» 16 – 20 сентября 2013 г. г. Новосибирск В.Ф. Куропатенко

3 РФЯЦВНИИТФ Механика взрыва – наука о поведении вещества под действием динамических нагрузок - получила мощный импульс в 50-х годах прошлого века в связи с разработкой ядерного оружия. При ударном сжатии или импульсном нагреве вещества потоком излучения возникают большие концентрации энергии. Они приводят к перестройке кристаллической структуры вещества, плавлению, испарению, диссоциации и ионизации. Необходимость моделировать свойства вещества в широком диапазоне изменения давления, плотности и температуры привела к бурному развитию компонентов триады эксперимент – модель – прогноз или от частного к общему, от общего к частному. 2

4 РФЯЦВНИИТФ К началу XX века были созданы основы МСС: Понятие сплошной среды. Закон сохранения массы Закон сохранения количества движения Закон сохранения энергии где Уравнения состояния или Определяющие уравнения (упругость, пластичность, вязкость). Инвариантность к преобразованию Галилея. 3

5 РФЯЦВНИИТФ До середины XX столетия, т.е. до появления ЭВМ, развитие МСС шло по пути аналитических решений. Резкое различие свойств агрегатных состояний веществ определило разделение МСС на ряд направлений (гидродинамика, газодинамика, аэродинамика, теория упругости, теория пластичности, теория детонации и др.). В каждом из направлений строились свои аналитические решения, с помощью которых получались требуемые промышленностью прогнозы характеристик физических процессов и оценивались параметры разрабатываемых машин. Несмотря на то, что аналитические решения получались для упрощённых постановок задач и упрощенных уравнений состояния, их ценность несомненна. В этот период были получены прекрасные аналитические решения, которые были востребованы при создании новой техники в авиации, кораблестроении, космонавтике, при решении оборонных задач. Это решение задачи о точечном взрыве, автомодельные решения о распространении ударных и детонационных волн и движениях оболочек, о разлете продуктов взрыва, о течениях несжимаемой вязкой жидкости и ряд других. 4

6 РФЯЦВНИИТФ Успехи механики сплошных сред, достигнутые к середине XX века в области создания моделей и аналитических решений, ярко изложены в следующих работах: 1.Я.Б.Зельдович. Теория ударных волн и введение в газодинамику. Изд.АН СССР, Н.Е.Кочин, И.А.Кибель, Н.В.Розе. Теоретическая гидромеханика. Гостехиздат, Г.Курант, К.Фридрихс. Сверхзвуковое течение и ударные волны. ИИЛ А.Ферри. Аэродинамика сверхзвуковых течений. Техтеорлит Л.Д.Ландау, Е.М.Лифшиц. Механика сплошных сред. Гостехиздат К.П.Станюкович. Неустановившиеся движения сплошной среды. Гостехиздат Ф.А.Баум, К.П.Станюкович, Б.И.Шехтер. Физика взрыва. Физматгиз, Л.И.Седов. Введение в механику сплошной среды. Физматгиз, В.Прагер. Введение в механику сплошных сред. ИИЛ

7 РФЯЦВНИИТФ Появление ЭВМ в середине XX столетия – это начало новой эры в механике. Для разработки ядерного оружия и освоения космоса срочно потребовались модели нового поколения, в которых учитывались бы реальные свойства вещества. Стала интенсивно формироваться Вычислительная механика сплошных сред. 6 УРС. Теор. модели УРС. Полуэмпирика Уравнения МСС. Законы сохранения Вычислительная механика сплошных сред ЭВМ Модели МСС Кинетические процессы Методы вычислений УРС. Эксперименты

8 РФЯЦВНИИТФ Появление ЭВМ создало новые условия для развития МСС. Появилась потребность в новых более сложных моделях среды, в новых УРС и определяющих уравнениях, в новых экспериментальных данных. Стала бурно развиваться физическая составляющая МСС: Теоретические модели свойств веществ (Томаса-Ферми, Саха, Хартри-Фока и др.) Полуэмпирические УРС, пригодные для расчётов конструкций. Кинетические модели фазовых переходов, разрушения. Новые задачи инициировали новые эксперименты: Ударноволновые и детонационные эксперименты для получения информации о свойствах веществ и процессов. Эксперименты по изучению перемешивания и турбулентности, разрушения и фазовых переходов. 7

9 РФЯЦВНИИТФ Достигнутый за полвека новой эры уровень моделей МСС и экспериментальных методов изучения свойств веществ отражен в литературе.Обзоры 1.Л.В. Альтшулер. Применение ударных волн в физике высоких давлений // УФН. 1965, т. 85, вып В.Б. Минцев, В.Е. Фортов. Взрывные ударные трубы // ТВТ. 1982, т. 80, 4. 3.А.В. Бушман, В.Е. Фортов. Модели уравнения состояния вещества // УФН. 1983, т. 140, вып Л.В. Альтшулер, Р.Ф. Трунин, К.К. Крупников, Н.В. Панов. Взрывные лабораторные устройства для исследования сжатия веществ в ударных волнах // УФН. 1996, т. 166, 5. Монографии 1. Цянь Cюэ-сень. Физическая механика. М.: Мир г. 2.Я.Б. Зельдович, Ю.П. Райзер. Физика ударных волн и высокотемпературных гидродинамических явлений. М.: Наука г. 3.Л.И.Седов. Механика сплошной среды, т.1, т.2, Наука, 1970 г. 4.С.Г. Андреев, А.В. Бабкин, Ф.А. Баум и др. Физика взрыва. т.1 и т. 2. Под редакцией Л.П. Орленко. М.: Физматлит г. 8

10 РФЯЦВНИИТФ К началу XXI века были экспериментально измерены зависимости между скоростью ударной волны и скоростью вещества за фронтом для большинства элементов и огромного количества соединений, сплавов и смесей. По ним с помощью законов сохранения строятся ударные адиабаты. Данные, полученные в ударно волновых экспериментах в период гг., изложены в справочнике: Р.Ф. Трунин, Л.Ф. Гударенко, М.В. Жерноклетов, Г.В. Симаков. Экспериментальные данные по ударно-волновому сжатию и адиабатическому расширению конденсированных веществ. – Саров.: РФЯЦ-ВНИИЭФ г. 9

11 РФЯЦВНИИТФ Новая наука Вычислительная МСС создалась и развивалась очень быстро из-за потребностей обороны. При разработке ядерного оружия специфика работ определялась тремя факторами: Невозможностью провести натурные испытания ядерного заряда в лабораторных условиях. Хроническим отставанием отечественной вычислительной техники по основным параметрам от ЭВМ США. Необходимостью создания конструкций ядерного оружия, не уступающих по своим характеристикам создаваемым в США. точности моделирования динамических процессов в ЯЗ Целью всех проводимых работ было создание системы моделей и методов для обеспечения высокой точности моделирования динамических процессов в ЯЗ. 10

12 РФЯЦВНИИТФ В Вычислительной МСС ищутся и находятся решения не дифференциальных, а разностных законов сохранения. Переход к разностным уравнениям порождает большой перечень обязательных исследований их свойств: 11 Характер аппроксимации Диссипативность Критерии устойчивости Монотонность Дистракцию Консервативность Вычислительная МСС состоит ФАКТИЧЕСКИ из двух блоков: разностных методов и моделей МСС. Основа вычислительной МСС – математическая модель – возможна только при условии тесного взаимодействия разностных методов и моделей МСС.

13 РФЯЦВНИИТФ 12 Адиабатическое ядро уравнения Эйлера-Гельмгольца и УРС.

14 РФЯЦВНИИТФ 13 В основе современных моделей сплошной среды лежит принцип аддитивности внутренней энергии где Каждый вид энергии определяет компоненту давления Энергия девиаторов определяется уравнением

15 14 Для каждой i-ой фазы вещества задаются сколь угодно сложные зависимости P i = P i (V i,E i ), T i = T i (V i,E i ), удовлетворяющие всем термодинамическим условиям. Границы фаз – это слабые разрывы изотерм и изэнтроп плазма плотная плазма жидкость газ тв. тело PXPX V K Ф4Ф4 Ф3Ф3 Ф1Ф1 P 0 Ф2Ф2 Схема агрегатных состояний в плоскости P, V Уравнение состояния железа в плоскости P, V (Дремов В.В.)

16 15 Каждый сильный разрыв заменяется непрерывным течением в слое конечной ширины (дистракция). В этом слое растет энтропия (диссипация энергии). Известно 4 принципиально разных механизма диссипации энергии. Год первой публикации Авторы метода расчета Страна Механизм диссипации энергии 1950Д. Нейман, Р. Рихтмайер [1]СШАПсевдовязкость 1954П. Лакс [2]СШААппроксим. вязкость 1957С.К. Годунов [3]СССРРаспад произв. разрыва 1960В.Ф. Куропатенко [4]СССРУравнения Гюгонио 1. Neumann J., Richtmayer R. A method for the numerical calculation of hydrodynamical shocks //J. Appl. Phys. – – V.21, 3 – pp Lax P.D. Weak solution of nonlinear hyperbolic equations and their numerical computations //Comn. Pure and Appl. Math. – – V.7 – pp Годунов С.К. Разностный метод расчета ударных волн//УМН т 12, вып.1.-с Куропатенко В.Ф. Метод расчета ударных волн //ДАН СССР.– – В.3, 4. – С

17 16 В 1965 г. в РФЯЦ-ВНИИТФ был создан неоднородный метод расчета нестационарных движений сплошных сред, в котором выделяются все существенные разрывы. Его применение устранило образование энтропийных следов возле контактных границ и в зонах взаимодействия ударных волн и резко повысило точность расчетов в целом. На рисунке приведена картина сильных, слабых и контактных разрывов в переменных r,t в сферической многооболочечной конструкции после взрыва наружного слоя ВВ.

18 РФЯЦВНИИТФ Сегодня быстродействие ЭВМ приближается к экзафлопному. Мат. моделирование применяется во всех естественных науках, в большинстве отраслей промышленности, в экономике, биологии, медицине и других областях. На ЭВМ рассчитываются: Новые конструкции самолетов, надводных и подводных кораблей, автомобилей, ракет и других образцов техники, Химические реакции в различных технологических процессах, Взаимодействие тел при высокоскоростных и низкоскоростных соударениях с целью отыскания средств защиты и оптимизации средств поражения, Горение топлива в камерах двигателей, Процессы в тепловых, атомных и других электростанциях с целью оптимизации существующих и поиска новых источников энергии, Движение крови в кровеносных сосудах. Проще назвать те области человеческой деятельности, где ЭВМ не применяются. 17

19 Любая модель является приближенной, тем более математическая модель. Отсутствие указания на приближенность модели в Вычислительной МСС порождает формально-математический подход к интерпретации результатов расчетов. Выдавая поток цифр, ЭВМ отрывает их от модели, как бы экранирует модель от пользователя и избавляет его от понимания её достоинств и недостатков. Доступность ЭВМ сделала мат. моделирование массовым занятием. Освоив ЭВМ, легко ввести в неё любую «новую» модель. Подобрав значения эмпирических параметров, можно даже описать несколько экспериментов, но на этом пути нельзя получить средство для прогнозирования других экспериментов. Такие результаты взаимодействия ЭВМ и пользователя вызывали и продолжают вызывать тревогу. Седов Л.И. (1976): «Теория моделирования деформируемых сред у нас и за границей заполнена мутными потоками дилетантской путаницы, которые сбивают с толку многих специалистов и учащихся» (ПМТФ, т. 40). Новожилов В.В. (1983): «Продолжают появляться без связи с потребностями практики всевозможные «мыслимые» модели сплошной среды, для упрощённых уравнений решаются краевые задачи, доказываются теоремы существования и единственности. Торжествует «механикоподобная математика». (Препринт 125 ЛГУ). Ситуация сегодня не улучшилась! 18

20 Все математические модели содержат эмпирические параметры, которые определяются не из теории, а из условия совпадения расчетов с базовым экспериментом. Конечная погрешность расчетов имеет вид: где - погрешность разностных уравнений - погрешность УРС, - погрешность определяющих уравнений, - погрешность состояния образца, - погрешность постановки расчета. Экспериментальные данные также имеют погрешность. Т.о. определение эмпирических параметров модели возможно только при взаимной компенсации указанных погрешностей. Для систем, далеких от базовой, результаты расчетов и экспериментов могут отличаться. 19

21 РФЯЦВНИИТФ В результате применения однородных методов, «размазывающих» ударные волны и слабые разрывы, в решении возникают энтропийные следы: При взаимодействии сильных разрывов друг с другом, При взаимодействии сильных разрывов с контактными границами, При распространении разрывов по сетке с переменным шагом. Энтропийные следы необратимы. Взаимодействие последующих ударных волн с энтропийными следами полностью искажает распад произвольного разрыва. Конечная погрешность возрастает. Она не уменьшается при увеличении числа точек сетки. Многие однородные методы сильно диссипативны, т.е. на непрерывных решениях энтропия меняется сильнее, чем на слабых ударных волнах. 20

22 РФЯЦВНИИТФ Современные УРС описывают фазовые переходы. Границы фаз неподвижны в пространстве термодинамических переменных. При моделировании процессов с фазовыми переходами «механическая» скорость движения вдоль изэнтропы НЕ коррелирует со скоростью образования новой фазы. Эта скорость должна рассчитываться в прикладной программе с помощью кинетического уравнения, когда вещество пересекло границу фазы. Примерно 90% применяемых в МСС программных комплексов этого делать не умеют. Это очень распространенная причина конечных ошибок. 21

23 22 На рисунке изображена изэнтропа разгрузки ABCD, пересекающая границу жидкости. В уравнении неразрывности отсутствует информация о фазовом состоянии вещества. В момент t вещество находится в точке A с параметрами P A, V A. Через время Δ t из разностного аналога уравнения (1) находится значение. Если алгоритм «умеет» определять точку B пересечения изэнтропы с бинодалью, то разгрузка пойдёт по изэнтропе в точку F (P F,V F ). Если нет, то будет найдено значение P E, и дальнейшая разгрузка пойдёт по другой изэнтропе EG, где S E < S F. Возникла конечная ошибка.

24 23 1. Одиночный откол. 2. Множественный откол. 3. Откол с дроблением вещества. + х + + х х Рис. 1 Рис. 2 Рис. 3 Однородные методы с «размазыванием» разрывов в принципе не могут описать ни одно из указанных решений. Погрешность математического моделирования откольного разрушения состоит из погрешности модели физ и погрешности численного метода мат = физ + мат. Для оценки мат используются аналитические решения для трёх возможных типов отколов :

25 24 Однородные методы: сильные и слабые разрывы «размазываются» (Пунктирные линии на Рис. 1.2). Неоднородные методы: сильные и слабые разрывы НЕ «размазываются» (Сплошные линии на Рис. 1.2). Рис. 1. Момент откола Рис. 2. Через t после откола При P_ < 0 на фронте УВ внутренняя энергия уменьшается Поскольку P = 0, V – V_< 0, P_< 0, то E – E_ < 0. В однородных же методах внутренняя энергия растёт Поскольку V – V_ 0, 0,5 P_ + q > 0, то E – E_ > 0. Любой однородный метод рассчитывает откол неправильно.

26 25 Фрагмент X-t-диаграммы, однородный метод

27 26 X-t-диаграмма. Формула линии уровня

28 27 Фрагмент X-t-диаграммы зоны дробления. Аналитическое решение. Расчет однородным методом с N= 1000.

29 28 Фрагмент X-t-диаграммы зоны дробления. Неоднородный метод. Выделяются слабые разрывы и ударные волны в точке разрушения. Сравнение результатов расчета с аналитическим решением

30 где Если сетка задана по геометрической прогрессии со знаменателем 2, то и при интегральная погрешность равна примерно 6%. 29 Все свойства разностной схемы математики проверяют на равномерных сетках. Все расчёты реальных систем проводятся на неравномерных сетках. На них при предельное решение отличается от точного решения. Пример. При, при на границе при x=0 Точное решение Аппроксимируем исходное уравнение разностным уравнением

31 30 Бытует великое заблуждение: Если хочешь повысить точность расчетов, увеличь число N точек сетки. Рис. Сходимость t p В действительности предельное решение, получающееся при всегда отличается от точного решения. При увеличении N результат приближается к предельному решению. На Рис. изображена зависимость времени образования откола t p в задаче об уединенном отколе от числа точек сетки N в расчете однородным методом.

32 РФЯЦВНИИТФ Пора прекратить подгонку под эксперимент с помощью моделей структурно однородных сред с простейшими УРС. Вычислительная механика должна прогнозировать количественные характеристики любого физического процесса, т.е. предсказывать результаты эксперимента без эмпирических параметров. Это возможно только с помощью моделей нового поколения. Двигаться к этой цели можно двумя путями. 1.Модели, основанные на гипотезе сплошной среды, далеко не исчерпали своих возможностей. Для уменьшения погрешности математического моделирования нужно совершенствовать каждую компоненту математической модели в отдельности. 2. Использовать модели, сочетающие континуальное описание вещества с дискретным: Модели частиц в ячейках Модели молекулярной динамики Модели мезоуровневой пластичности Модели многокомпонентных сред. 31

33 1. Внедрение неоднородных методов, выделяющих сильные, слабые и контактные разрывы в одномерной, двумерной и трехмерной постановках и моделирующих взаимодействие разрывов под произвольными углами. 2. Создание уравнений состояния веществ, не содержащих эмпирических констант. 3. Настройка вычислительных алгоритмов на особенности УРС: пересечение ударных адиабат и изэнтроп с границами фаз, расщепление ударных волн и волн разрежения, образование и распространение ударных волн разрежения. 4. Создание кинетических уравнений для моделирования фазовых переходов, химических реакций в ВВ, развития поврежденности кристаллической структуры материалов. 5. Создание уравнений, моделирующих взаимодействия компонентов смеси при различных комбинациях агрегатных состояний. 32 Смесь не имеет уравнения состояния! Каждый компонент смеси имеет свой УРС!

34 РФЯЦВНИИТФ Гипотеза о непрерывности сплошной среды дала возможность записать законы сохранения в виде дифференциальных уравнений в частных производных и создать механику сплошной среды, а затем и вычислительную механику сплошной среды. Модели МСС с эмпирическими параметрами были эффективно использованы при создании кораблей, самолетов, ракет, при разработке технологических процессов во многих отраслях промышленности. Модель сплошной среды верно служит людям больше 200 лет. 33 Сегодняшние возможности ЭВМ позволяют начать создание моделей нового поколения без эмпирических параметров. Жизнь заставляет нас не описывать эксперимент, а прогнозировать результат.

35 34 Благодарю за внимание

36 Поздравляю Владимира Михайловича с юбилеем!!!