- •Оглавление
- •Предисловие
- •Введение
- •Р а з д е л 1.Общая характеристика учебных версий ппп Flow Vision и Gas Dynamics Tool
- •§ 1.Физико-математические модели
- •§ 2.Граничные и начальные условия
- •§ 3.Особенности численных расчетов
- •§ 1.Физико-математические модели
- •§ 2.Граничные и начальные условия
- •§ 3.Особенности численных расчетов
- •Р а з д е л 2.Решение учебных задач с использованием пакета Flow Vision г л а в а 1.Cостав и назначение основных моделей пакета
- •§ 1.Препроцессор
- •§ 2.Солвер
- •§ 3.Постпроцессор
- •Г л а в а 2. Алгоритм моделирования в пакете Flow Vision
- •§ 1.Геометрический препроцессор (Solid Works)
- •§ 2.Физико-математическая постановка задачи
- •§ 3.Подготовка к численному моделированию
- •§ 4.Моделирование с помощью солвера
- •§ 5.Подготовка к визуализации результатов
- •§ 6.Визуализация скалярных полей
- •§ 7.Визуализация отдельных числовых значений
- •§ 8.Визуализация векторного поля скорости
- •§ 9.Представление результатов и подготовка отчета
- •Г л а в а 3.Течение вязкой жидкости в прямом плоском канале § 1.Основные соотношения
- •§ 2.Постановка задачи
- •§ 3.Задание
- •Результаты расчета скорости и длины установления течения
- •§ 4.Представление и анализ результатов
- •Г л а в а 4.Обтекание круглого цилиндра вязкой несжимаемой жидкостью § 1.Основные соотношения
- •§ 2.Постановка задачи
- •§ 3.Задание
- •Результаты расчета и сравнение с теорией силы сопротивления, испытываемой цилиндром
- •§ 4.Представление результатов
- •Г л а в а 5.Течение жидкости в канале Переменного сечения § 1.Основные соотношения
- •§ 2.Постановка задачи
- •§ 3.Задание
- •Потери давления при сужении (расширении) канала
- •§ 4.Представление результатов
- •Г л а в а 6.Обтекание эллиптического цилиндра и плоской пластины идеальной несжимаемой жидкостью § 1.Основные соотношения
- •§ 2.Постановка задачи
- •§ 3.Задание
- •Расчет и сравнение с теорией силы и момента сил
- •§ 4.Представление результатов
- •Г л а в а 7.Удар воздуха о торец пластины § 1.Основные соотношения
- •§ 2.Постановка задачи
- •§ 3.Задание
- •Расчет и сравнение с теорией давления при ударе
- •Р а з д е л 3.Решение учебных задач с использованием пакета Gas Dynamics Tool г л а в а 1.Алгоритм моделирования в Gas Dynamics Tool
- •§ 1.Выбор параметров пакета
- •Набор параметров, задаваемых для расчета в пакете gdt
- •§ 2.Визуализация с помощью постпроцессора
- •§ 3.Проведение расчетов и представление результатов
- •Г л а в а 2.Ударная волна
- •§ 1.Основные соотношения
- •§ 2.Постановка задачи
- •§ 3.Особенности выполнения задания
- •Г л а в а 3.Истечение из сопла
- •§ 1.Основные соотношения
- •§ 2.Постановка задачи
- •§ 3.Задание и особенности его выполнения
- •Расчетные области и значения параметров, задаваемых в них
- •§ 4.Представление результатов
- •Параметры на оси сопла
- •Г л а в а 4.Сверхзвуковой диффузор
- •§ 1.Потери полного давления в ударной волне
- •§ 2.Постановка задачи
- •§ 3.Моделирование диффузора
- •Координаты для построения поверхностей
- •§ 4.Представление результатов
- •Г л а в а 5.Удар воздуха о торец пластины
- •§ 1.Моделирование в Gas Dynamics Tool
- •§ 2.Представление результатов
- •Расчет и сравнение с теорией давления при ударе
- •Г л а в а 6.Течения с подводом тепла и детонация (gdt) § 1.Основные соотношения
- •§ 2.Оценки параметров и диаграммное представление
- •§ 3.Постановка задачи
- •§ 4. Представление результатов
- •Сравнение с теорией параметров расчета при течении с подводом тепла
- •Заключение
- •Приложение применение теории функций комплексной переменной к решению задачи обтекания идеальной несжимаемой жидкостью эллиптического цилиндра и пластины
- •Список литературы
§ 1.Основные соотношения
Уравнения, описывающие изменения параметров в ударной волне, имеются в большинстве учебников по МЖГ, например [3, 4, 6, 8]. В пособии приводится краткий вывод уравнений импульсов и энергии.
Уравнение импульсов для стационарного течения
В случае
движения ударной волны в канале уравнения
удобно записывать в системе координат,
в которой ударная волна покоится. В этом
случае течение можно считать стационарным.
Уравнение импульсов получается из
второго закона Ньютона F = m
Выразим силу и массу через давление, плотность и единичную площадь:
;
;
.
(2.1)
Тогда уравнение движения запишется в виде
.
(2.2)
В канале постоянного сечения произведение плотности на скорость постоянно, поэтому получаем следующий интеграл уравнения импульсов для прямой ударной волны:
(2.3)
Уравнение энергии
В отсутствие внешних и химических источников тепла с использованием уравнения импульсов получаем из первого начала термодинамики:
,
(2.4)
(2.5)
где E ― внутренняя энергия, h ― энтальпия. Интегририруя (2.5), получим
const.
(2.6)
Уравнение (2.6) описывает адиабатический процесс (в том числе и необратимый). Уравнение энергии можно также записать с использованием второго начала термодинамики:
(2.7)
где S и Т соответственно энтропия и абсолютная температура.
Уравнения сохранения для ударной волны
Из уравнений (2.3) ― (2.7) можно получить соотношения для прямой ударной волны. Обозначим индексом «1» параметры до ударной волны, а индексом «2» параметры за ударной волной. Тогда в системе координат, связанной с ударной волной, получаются следующие выражения, представляющие собой уравнения сохранения:
,
(2.8)
,
(2.9)
.
(2.10)
Эти соотношения необходимо дополнить уравнением состояния, которое в данном случае будем рассматривать в виде уравнения для идеального газа:
(2.11)
Решение системы уравнений (2.8) – (2.10) приводит к следующему известному под названием адиабата Гюгонио соотношению для прямой ударной волны:
.
(2.12)
Здесь V ― удельный объем. Соотношение (2.12) можно переписать через число M1 набегающего потока и соответственно для отношения плотностей «после» и «до» ударной волны:
,
(2.13)
.
(2.14)
Отметим, что представление соотношения параметров на ударной волне через удельный объем (например, формула (2.12)), полезно для анализа так называемых PV-диаграмм; что касается расчетов, то удобнее пользоваться соотношениями, выраженными через плотность (формула (2.14)).
Полезно привести также уравнения сохранения для газа, покоящегося перед ударной волной:
1Us = 2(Us – U2), (2.15)
P1 + 1Us2 = P2 + 2 (Us – U2)2 , (2.16)
h1 + Us2/2 = h2 + (Us – U2)2/2. (2.17)
Для скорости газа за ударной волной имеем выражение
,
(2.18)
индекс “1” относится к параметрам перед ударной волной, индекс“2” относится к параметрам за ударной волной, индекс “s” относится к фронту ударной волны. , P, U, h, a, M – плотность, давление, скорость, энтальпия, скорость звука, число М соответственно.
