Конструктивный подход к численному решению квазилинейных уравнений переноса

1. Конструктивный подход к численному решению квазилинейных уравнений переноса А.П. Фаворский1, А.М. Галанина2, В.А. Исаков3 Московский Государственный Университет им. М.В. Ломоносова _________________________________ 1 МГУ им. М.В. Ломоносова, профессор, 117899, Москва, Воробьёвы горы, д.1, e-mail: galaninaanna@gmail.com 2МГУ им. М.В. Ломоносова, аспирант, 117899, Москва, Воробьёвы горы, д.1, e-mail: galaninaanna@gmail.com 3 МГУ им. М.В. Ломоносова, аспирант, 117899, Москва, Воробьёвы горы, д.1, e-mail: victorisakov88@gmail.com

2. Введение • Уравнения переноса: • составляют неотъемлемую часть математического описания многих явлений; • оказываются ключевым моментом при рассмотрении и построении решения многих задач (например, перенос инвариантов Римана в газовой динамике); • уравнение переноса суть простой и содержательный представитель уравнений гиперболического типа, удобный для обсуждения методов численного решения, включая квазилинейный случай;

3. Требования, предъявляемые к построению численного решения: • явность; • сохранение монотонности профиля решения; • приемлемая точность на «гладких» участках решения; • соблюдение законов сохранения в потоковой форме; • отсутствие специальных регуляризаторов (искусственной вязкости, например);

4. Постановка задачи Квазилинейное уравнение переноса (1) Дивергентная потоковая форма уравнения переноса: где (2) Проводится построение финитного решения на постоянном фоне в предположении отсутствия влияния граничных условий

5. Построение численного решения Проводятся линии параллельные фону. Решение в начальный момент времени заменяется кусочно - линейной функцией где Значения остаются постоянными, перемещаясь вдоль линий с постоянной характеристической скоростью

6. Свойства кусочно-линейной функции • Эволюция кусочно-линейной функции: • воспроизводится точно; • Замена произвольной достаточно гладкой исходной функции кусочно-линейной происходит с порядком при достаточно малых шагах по времени ; • Однозначное решение выстраивается до момента возникновения «градиентной катастрофы». В последующие моменты времени следует вводить разрыв функции исходя из условия Гюгонио на фронте волны. • Рассмотренный подход распространяется на квазилинейные уравнения более общего вида: монотонная функция

7. Построение локального, линейного сплайна на отрезке на момент времени

8. Функция u(x,t) на отрезке при заменяется локально-линейным сплайном где Такая кусочно-линейная функция не нарушает монотонности профиля на отрезке [1]. характеризует собой наклон сплайна. Величина соответствует построению схемы первого порядка

9. Расслоение сплайновой функции Для вычисления необходимых при построении разностной схемы интегральных потоков на границах ячейки линейный сплайн заменяется структурой ступенчатых функций («малых возмущений»), расположенных последовательно одна на другой, отсчитывая от фона.

10. Скорость перемещения «кирпича» вдоль оси x определяется по формуле где и - значения на соответствующих слоях. «Малые возмущения» распространяются вдоль характеристики. Если скорость , то соответствующий «кирпич» имеет возможность пересечь границу ячейки . В этом случае из ячейки в ячейку перейдёт интегральное количество функции равное

11. Вычисление интегральных потоков Интегральный поток за время через границу расчетной ячейки , соответствующий перемещению «кирпичей» из ячейки в ячейку , равен Фоновый поток определяется по формуле Результирующий интегральный поток через границу равен

12. Разностная схема • имеет второй порядок точности на гладких решениях; • схема консервативна; • не нарушает монотонность профиля волны, включая разрывные решения; • является устойчивой при соблюдении условия Куранта: • учёт задания граничных значений функции на границах не представляет затруднений; • схема не содержит искусственных регуляризаторов;

13. Результаты численных расчетов Квазилинейное уравнение переноса. Точное решение

14. Расчет по сплайн-схеме: синусоидальный импульс

15. Расчет по сплайн-схеме: ступенчатый импульс

16. Построение схемы для уравнений газовой динамики Одномерные уравнения газовой динамики в дивергентной форме где - любая из функций - потоки: Система уравнений замыкается уравнением состояния где - показатель адиабаты.

17. Построение разностной схемы Численное решение системы уравнений строится в классе сеточных функций , отнесенных к центрам ячеекравномерной по направлениям координат и прямоугольной сетки с шагами и. Границы пространственных ячеек проходят через полуцелые точки . Построение схемы проводится интегро-интерполяционным методом. Все уравнения системы последовательно интегрируются по прямоугольной пространственно-временной ячейке:

18. Интегральные балансные соотношения Выражают законы сохранения массы, импульса и энергии газа. Величины: - средняя по ячейке на момент времени объемная плотность; - объемная плотностьимпульса; -объемная плотностьполной энергии;

19. Величины - средняя по массе ячейки скорость; - средняя по массе ячейки внутренняя энергия; равны интегральным по времени потокам в сечениях - дискретный аналог уравнения состояния;

20. Аппроксимация интегрального потока Локально-линейная сплайн-реконструкция функций заменяем локально-линейным сплайном, моделирующим её поведение в пределах ячейки на момент времени : Функцию характеризует угол наклона сплайна и вычисляется по формуле: где или Величина где

21. Расслоение (разбиение) линейных сплайн-функций Горизонтальное сечение, ближайшее к среднему между значениями сплайн-функций и на границе раздела ячеек и , назовём общим постоянным фоном .

22. Акустическое приближение для малых возмущений В течение времени от до каждое малое возмущение распадается на волны Римана, бегущие каждая по своему фону в соответствии с известным решением уравнений акустики: где

23. Малые возмущения на каждом слое в соответствии с решением линеаризованной системы имеют следующий вид: где - соответствующие характеристики, а через обозначены отклонения плотности, скорости и давления от своих фоновых значений на момент времени .

24. Вычисление интегрального потока Интегральный по времени поток импульса через границу из ячейки в ячейку за время равен Здесь - поток, обусловленный общим постоянным фоном с параметрами Вклад от послойных возмущений связан с вычислением интегралов от ступенчатых функций возмущений

25. Представление алгоритма в переменных Лагранжа Алгоритм численного решения уравнений газовой динамики в переменных Лагранжа не отличается от изложенного выше и даже несколько проще, поскольку число волн Римана сокращается с трёх до двух.

26. Результаты численных расчётов Расчет задачи с гладким начальным профилем. Приводятся графики точного и приближенного решения (функция скорости) на последовательные моменты времени.

27. Пример расчета задачи о распаде произвольного разрыва (a) – функция скорости в эйлеровых переменных на фиксированный момент времени при различных шагах и ; (б) – функция давления в лагранжевых переменных в различные моменты времени при фиксированных шагах и .

28. Расчёт цилиндрического взрыва

29. Распространение газовых струй от двигателей самолёта в отсутствие экрана

30. Взаимодействие струи с отбойником экспериментатор И.К. Ермолаев

31. Расчет взаимодействия газовых струй, вылетающих из двигателя самолёта, с отражающим экраном

33. Свойства разностной схемы • схема является явной; • не содержит искусственных регуляризаторов; • воспроизводит гладкие решения со вторым порядком точности на реальных сетках; • численное решение слабо зависит от числа Куранта при соблюдении условия устойчивости: • где ; • схема монотонна на гладких решениях и квазимонотонна в окрестности разрывов, структура фронта которых занимает 2 – 3 интервала расчётной сетки; • схема допускает обобщение на случай двух и трёх пространственных измерений.

34. Список литературы Фаворский А.П., Тыглиян М.А., Тюрина Н.Н., Галанина А.М., Исаков В.А. Численное моделирование распространения гемодинамических импульсов // Мат. Моделирование. 2009. Т. 21 № 12. с. 21-34. Фаворский А.П., Тыглиян М.А., Тюрина Н.Н., Галанина А.М., Исаков В.А. Численное моделирование распространения акустических импульсов в гемодинамике // Дифференциальные уравнения. 2009. Т. 45 № 8. с. 1179-1187. Абакумов М.В., Галанина А.М., Исаков В.А., Тюрина Н.Н., Фаворский А.П., Хруленко А.Б. Квазиакустическая схема для уравнений Эйлера газовой динамики // Дифференциальные уравнения. 2011. Т. 47 № 8. с. 1092-1098.

35. Спасибо за внимание!