Презентация на тему: Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,

Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,
1/32
Средняя оценка: 4.1/5 (всего оценок: 44)
Код скопирован в буфер обмена
Скачать (5716 Кб)
1

Первый слайд презентации

Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ, Долгопрудный к.ф.-м.н. Уткин Павел Сергеевич e-mail: utkin@icad.org.ru, pavel_utk@mail.ru (926) 2766560 Система уравнений газовой динамики и разностные схемы для ее решения Семинар № 5

Изображение слайда
2

Слайд 2

Краткое содержание предыдущих семинаров 2 Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5. Семинар № 1 (0 7.02.1 4 ). Некоторые разностные схемы для решения линейного уравнения переноса. Семинар № 2 ( 14.02.1 4 ). Построение схем для решения линейного уравнения переноса в пространстве неопределенных коэффициентов. Семинар № 3 (2 1.02.1 4 ). Понятие монотонности разностных схем. Теорема Годунова. Семинар № 4 (2 8.02.1 4 ). Обобщение на случай системы уравнений гиперболического типа.

Изображение слайда
3

Слайд 3

Система уравнений газовой динамики 3 Дивергентная форма записи (в форме законов сохранения) Вектор консервативных переменных Вектор потоков Полная энергия Внутренняя энергия (уравнение состояния) Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
4

Слайд 4

Анализ системы уравнений газовой динамики 4 Характеристическая форма Скорость звука Матрица Якоби Куликовский А.Г., Погорелов Н.В., Семенов А.Ю. Математические вопросы численного решения гиперболических систем уравнений. – М. : Физматлит, 2001. Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
5

Слайд 5

Гиперболическая система уравнений 5 Наличие конечной скорости распространения бесконечно слабых возмущений Возможность существования разрывных решений даже для гладких начальных данных Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
6

Слайд 6

Метод конечных объемов (1) 6 V m ( q 1 m, q 2 m, … ) 0 x y z Для произвольной компоненты вектора консервативных переменных : или Проинтегрируем по объему ячейки расчетной сетки и по времени : Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
7

Слайд 7

Метод конечных объемов (2) 7 V m ( q 1 m, q 2 m, … ) σ n σ 0 x y z Возьмем интеграл в первом выражении и применим теорему Остроградского-Гаусса к расчету интеграла во втором : S m ( q 1 k, q 2 k, … ) Аппроксимируем поверхностный интеграл через сумму интегралов по граням ячейки : Основной вопрос – как определять численный поток? Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
8

Слайд 8

8 Постановка задачи о распаде произвольного разрыва Задача Коши для системы уравнений газовой динамики с разрывом первого рода в начальных данных Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
9

Слайд 9

9 Возможные конфигурации решения Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
10

Слайд 10

Соотношения на разрыве 10 Соотношения Ренкина – Гюгонио : Разрывы Контактные – нет потока массы вещества через разрыв Ударные волны Обтекание тела сверхзвуковым потоком Ван-Дайк М. Альбом течений жидкости и газа. – М. : Мир, 1986. Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
11

Слайд 11

Элементарная теория ударных волн 11 Адиабата Гюгонио p p 0 адиабата Гюгонио адиабата Пуассона В ударной волне газ нельзя сжать больше, чем в ( γ + 1 ) / ( γ – 1 ) раз. Ударная волна бесконечно малой интенсивности распространяется относительно газа со скоростью звука. Фронт ударной волны распространяется относительно фона со сверхзвуковой скоростью. Теорема Цемплена : не существует ударных волн разрежения. Δ S = 0 Самарский А.А., Попов Ю.И. Разностные методы решения задач газовой динамики. – М. : Наука, 1992. Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
12

Слайд 12

Соотношения для ударных волн и волн разрежения 12 «левая» УВ «правая» УВ «левая» ВР «правая» ВР Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
13

Слайд 13

Метод Ньютона для нахождения давления на контактном разрыве 13 Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
14

Слайд 14

14 Система тестов Toro E.F. Riemann Solvers and Numerical Methods for Fluid Dynamics. – Springer, 19 99. Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
15

Слайд 15

15 Тест 1. Задача Сода. Анализ (p-v)- диаграммы. p L = 1.0 U L = 0.0 ρ L = 1.0 γ L = 1.4 p R = 0.1 U R = 0.0 ρ R = 0.125 γ R = 1.4 Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
16

Слайд 16

16 Тест 1. Задача Сода. Точное решение. Распределение плотности. Начальная плотность слева от разрыва Начальная плотность справа от разрыва Ударная волна (УВ) Контактный разрыв (КР) Волна разрежения (ВР) «Голова» ВР «Хвост» ВР p L = 1.0 U L = 0.0 ρ L = 1.0 γ L = 1.4 p R = 0.1 U R = 0.0 ρ R = 0.125 γ R = 1.4 x 0 = 0.5 t = 0.2 Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
17

Слайд 17

17 Начальная скорость слева от разрыва Начальная скорость справа от разрыва Ударная волна Скорость на контактном разрыве непрерывна Волна разрежения «Голова» ВР «Хвост» ВР p L = 1.0 U L = 0.0 ρ L = 1.0 γ L = 1.4 p R = 0.1 U R = 0.0 ρ R = 0.125 γ R = 1.4 x 0 = 0.5 t = 0.2 Тест 1. Задача Сода. Точное решение. Распределение скорости. Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
18

Слайд 18

18 Начальное давление слева от разрыва Начальное давление справа от разрыва Ударная волна Давление на контактном разрыве непрерывно Волна разрежения «Голова» ВР «Хвост» ВР p L = 1.0 U L = 0.0 ρ L = 1.0 γ L = 1.4 p R = 0.1 U R = 0.0 ρ R = 0.125 γ R = 1.4 x 0 = 0.5 t = 0.2 Тест 1. Задача Сода. Точное решение. Распределение давления. Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
19

Слайд 19

19 Начальная энергия слева от разрыва Начальная энергия справа от разрыва Ударная волна Волна разрежения «Голова» ВР «Хвост» ВР Контактный разрыв p L = 1.0 U L = 0.0 ρ L = 1.0 γ L = 1.4 p R = 0.1 U R = 0.0 ρ R = 0.125 γ R = 1.4 x 0 = 0.5 t = 0.2 Тест 1. Задача Сода. Точное решение. Распределение внутренней энергии. Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
20

Слайд 20

20 Алгоритм построения точного решения 1. Определяем по (p-v)- диаграмме конфигурацию, возникающую при распаде. 2. В результате решения нелинейного алгебраического уравнения методом Ньютона ищем давление на контактном разрыве. 3. Определяем оставшиеся параметры – скорости ударных волн и наклоны крайних характеристик, описывающих веер волны разрежения. Самарский А.А., Попов Ю.П. Разностные методы решения задач газовой динамики. – М. : Наука, 1992. Годунов С.К. и др. Численное решение многомерных задач газовой динамики. – М. : Наука, 1976. Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
21

Слайд 21

Тестирование схемы С.К. Годунова решения уравнений газовой динамики 21 Сергей Константинович Годунов род. 1929 г., академик РАН Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
22

Слайд 22

22 p L = 1.0 U L = 0.0 ρ L = 1.0 γ L = 1.4 p R = 0.1 U R = 0.0 ρ R = 0.125 γ R = 1.4 x 0 = 0.5 t = 0.2 Тест 1. Задача Сода. Метод Годунова. Распределение плотности. Δ x = 0. 01 Δ t = 0. 001 УВ размазывается на ~ 9 ячеек КР размазывается на ~ 15 ячеек Погрешности в описании ВР Схема С.К. Годунова : 1-ый порядок аппроксимации Монотонность Физичность результатов (есть исключения) Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
23

Слайд 23

23 p L = 1.0 U L = 0.0 ρ L = 1.0 γ L = 1.4 p R = 0.1 U R = 0.0 ρ R = 0.125 γ R = 1.4 x 0 = 0.5 t = 0.2 Тест 1. Задача Сода. Метод Годунова. Распределение скорости. Δ x = 0. 01 Δ t = 0. 001 Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
24

Слайд 24

24 p L = 1.0 U L = 0.0 ρ L = 1.0 γ L = 1.4 p R = 0.1 U R = 0.0 ρ R = 0.125 γ R = 1.4 x 0 = 0.5 t = 0.2 Тест 1. Задача Сода. Метод Годунова. Распределение давления. Δ x = 0. 01 Δ t = 0. 001 Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
25

Слайд 25

25 p L = 1.0 U L = 0.0 ρ L = 1.0 γ L = 1.4 p R = 0.1 U R = 0.0 ρ R = 0.125 γ R = 1.4 x 0 = 0.5 t = 0.2 Тест 1. Задача Сода. Метод Годунова. Распределение внутренней энергии. Δ x = 0. 01 Δ t = 0. 001 Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
26

Слайд 26

26 Тест 1. Задача Сода. Сравнение методов Годунова и КИР (1). Распределение плотности Распределение скорости Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
27

Слайд 27

27 Тест 1. Задача Сода. Сравнение методов Годунова и КИР (2). Распределение давления Распределение внутренней энергии Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
28

Слайд 28

28 Тест 2. Метод Роу ( Roe ). Toro E.F. Riemann Solvers and Numerical Methods for Fluid Dynamics. – Springer, 19 99. Проблемы с законом неубывания энтропии Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
29

Слайд 29

29 Тест 2. Метод Роу с энтропийной коррекцией. Toro E.F. Riemann Solvers and Numerical Methods for Fluid Dynamics. – Springer, 19 99. Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
30

Слайд 30

30 Тест 1. Метод Лакса-Вендроффа. Sod G.A. A Survey of Several Finite Difference Methods for Systems of Nonlinear Hyperbolic Conservation Laws // J. Comp. Phys. – 1978. – V. 27. – P. 1 – 31. Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
31

Слайд 31

31 Тест 1. Метод Русанова. Sod G.A. A Survey of Several Finite Difference Methods for Systems of Nonlinear Hyperbolic Conservation Laws // J. Comp. Phys. – 1978. – V. 27. – P. 1 – 31. Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда
32

Последний слайд презентации: Нелинейные вычислительные процессы 7 марта 201 4 г., МФТИ,

32 Монотонные схемы повышенного порядка аппроксимации метод Годунова 1-го порядка аппроксимации метод Годунова повышенного порядка аппроксимации точное решение Уткин П.С. Нелинейные вычислительные процессы. Семинар № 5.

Изображение слайда