Параметры сопла Лаваля
Расчет сопла Лаваля с помощью газодинамических функций: проектирование дозвукового и сверхзвукового участков. Параметры течения газа по соплу. Расчет крыльевого профиля в среде Gas2. Определение профиля методом скачков уплотнения и волн разряжения.
Рубрика | Физика и энергетика |
Вид | курсовая работа |
Язык | русский |
Дата добавления | 07.08.2013 |
Размер файла | 1,1 M |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Федеральное агентство по образованию
Государственное общеобразовательное учреждение
высшего профессионального образования
Тульский государственный университет
Кафедра «Ракетостроение»
Курсовая по предмету аэрогазодинамика
Выполнил ст. гр.121362 Киричек И.В.
Проверил Худяков М.В.
Тула 2008
Расчет сопла Лаваля с помощью газодинамических функций
Параметры газа в камере сгорания
№ варианта |
P (МПА) |
T(K) |
Dвых |
h(км) |
k |
|
5 |
8 |
2000 |
0.2 |
1 |
1.4 |
Расчет параметров во входном критическом и выходном сечениях сопла.
Подсчитаем площадь выходного сечения сопла:
Запишем газодинамические функции
P0, T0, 0, параметры торможения (в камере)
FВХ, FКР, FA- площади входного выходного и критического сечений
RВХ, RКР, RA- радиусы входного выходного и критического сечений
DВХ, DКР, DВЫХ- радиусы входного выходного и критического сечений
Рассчитаем параметры для критического сечения сопла
;;
;
;
;
Рассчитаем параметры для входного сечения сопла
Принимаем
;;
;
;
;
;
Проектирование дозвукового участка сопла Лаваля
Рассчитаем длину дозвукового участка:
Спрофилируем дозвуковую часть сопла по формуле Виташинского
X-координата сечения
X |
0 |
0.04 |
0.08 |
0.12 |
0.16 |
0.2 |
|
R |
0.094 |
0.074 |
0.053 |
0.042 |
0.038 |
0.037 |
Проектирование сверхзвукового участка сопла
Примем согласно экспериментальным данным полуугол раструба сопла ,тогда найдем длину сверхзвуковой части (конического сопла). Она равна:
.
Спроектируем профилированное сопло, которое обеспечивает получение на выходе параллельного сверхзвукового потока с заданным числом M.
проведем через точку на стенке сопла , соответствующей критическому сечению ,прямую под углом к оси и продолжим ее до пересечения с осью (точка O)
из точки O проведем близко расположенные радиальные прямые ,получим
отложим от критического сечения на оси расстояние равное длине конического сопла (точка А)
проведем из точки А характеристики. Они проведены под углом равным15'
найдем пересечение радиальной прямой с характеристикой (точка)
Для сечения проходящего через А1 используя газодинамическую ф-ю из таблицы по значению находят значение , и из точкипроводим характеристики под углом (откладывается угол от) аналогично для других точек
Расчет параметров течения газа по соплу
Расчет параметров на дозвуковом участке сопла ведем по следующей схеме:
ранее дозвуковая часть была разбита на 6 частей, согласно данному разбиению выбираем соответствующие сечения;
считаем для каждого из полученных сечений;
для полученного значения находим из таблицы газодинамических функций число
по найденным числам можно найти параметры ,,, через газодинамические функции по сечениям
Расчет параметров на сверхзвуковом участке сопла ведем следующим образом:
на характеристике в точкахзначение числа Маха уже было найдено ранее; согласно данным значениям по газодинамическим функциям можно найти значение параметров ,,,в точках
значение параметров в точках и точках будут соответственно идентичны, так как соответствующие точки (например) будут лежать на одной характеристике, по построению, поэтому параметрам на стенке сопла будут соответствовать X точек
Найдем величину скорости звука в зависимости от выбранного сечения по формуле:
Величину скорости находим из зависимости
В результате вычислений получили следующие данные:
0.04 |
0.8 |
0.12 |
0.16 |
0.2 |
0,23 |
0,28 |
0,34 |
0,44 |
0,63 |
||
7. 896 |
7. 496 |
6.624 |
5.552 |
4.224 |
1.328 |
0.712 |
0.432 |
0.1192 |
0.088 |
||
13.76 |
13.27 |
12.14 |
10.71 |
8.81 |
3.86 |
2.46 |
1.74 |
0.97 |
0.57 |
||
0.137 |
0.304 |
0.525 |
0.74 |
1 |
1.827 |
2.23 |
2.544 |
3.070 |
3.595 |
||
1992 |
1962 |
1894 |
1802 |
1672 |
1198 |
1008 |
870 |
692 |
556 |
||
123 |
270 |
458 |
630 |
820 |
1268 |
1418 |
1504 |
1618 |
1700 |
||
895 |
888 |
872 |
851 |
820 |
694 |
636 |
591 |
527 |
473 |
Построим по полученным данным графики распределения давления, плотности, температуры, безразмерной скорости и скорости звука по длине сопла.
Расчет сопла Лаваля в среде Gas2
Структурная схема:
Сначала строим контур сопла. Далее мы разбиваем область внутри контура на ячейки с помощью линейной аппроксимации. На входном сечении сопла задаем границу втекания потока, а на выходной границу вытекания.
Так как, задача осесимметричная берем цилиндрическую систему координат.
На границе втекания задаем давление и температуру во входном сечении сопла. Учитываем в расчетах вязкость.
Схема распределения давления газов по соплу.
Схема распределения температуры газов по соплу.
Схема распределения плотности газов по соплу.
Высчитываем погрешность для распределения давления:
Высчитываем погрешность для распределения температуры:
Высчитываем погрешность для распределения плотности:
Расчет крыльевого профиля методом скачков уплотнения и волн разряжения.
Профиль симметричный с хордой в = 150 мм и максимальной толщиной с= 14 мм.
Параметры потока обтекающего крыловой профиль
№ варианта |
№профиля |
M |
P(МПА) |
T(K) |
k |
Угол атаки |
Угол атаки |
Угол атаки |
|
13 |
4 |
3.2 |
0.1 |
373 |
1.4 |
0 |
4 |
-2 |
Размещено на http://www.allbest.ru/
Рис 1 Торможение потока на скачке уплотнения.
Размещено на http://www.allbest.ru/
Рис 2 Ускорение потока на волне разрежения.
Расчет обтекания крыльевого профиля при заданных параметрах набегающего потока выполняется точным методом скачков - волн разрежения.
Схемы потоков и используемые параметры газа показаны на ( Рис 1,2)
1 Определяем параметры потока на каждом прямолинейном участке
рассчитываем параметры потока на скачках уплотнения
рассчитываем параметры потока на волнах расширения
Расчет параметров потока на скачках уплотнения проводим следующим образом: при заданном угле атаки и профиле находим угол отклонения потока(Рис 1),далее находим угол косого скачка , который зависит от величины скорости Маха и от угла (Cм газодинамические ф-ии ),затем ищем отношение давления после отклонения потока к давлению до него(),
Зная значение , находим значение ,далее находим и по газодинамическим функциям находим . Найдём значение скоростного напора непосредственно перед обтеканием крыла
V=Ma где -скорость звука
0 |
3.2 |
4 |
2110' |
2.9 |
1.105 |
0.1105 |
0.0105 |
0.015 |
|
4 |
3.2 |
0 |
1813' |
3.2 |
1 |
0.1 |
0 |
0 |
|
4 |
3.2 |
8 |
2440' |
2.8 |
1.394 |
0.1394 |
0.0394 |
0.055 |
|
-2 |
3.2 |
6 |
2250' |
2.95 |
1.32 |
0.132 |
0.032 |
0.045 |
|
-2 |
3.2 |
2 |
1945' |
3.1 |
1.076 |
0.1076 |
0.0076 |
0.011 |
Расчет параметров потока на волнах разрежения проводим следующим образом: при заданном угле атаки , профиле и скорости Маха находим углы ,(смотри Рис 2),далее находим фиктивный угол (из таблицы газодинамических функции)
Затем прибавляем угол ,далее находим по таблице скорости Маха отношение давления после отклонения потока к давлению до него ()будет равно
Зная значение , находим значение ,далее находим. Вычислим значение скоростного напора непосредственно перед обтеканием крыла
(-1) |
|||||||||||
0 |
2.9 |
2006' |
48 |
5634' |
3,34 |
1724' |
0.529 |
0.058 |
0.0525 |
0.073 |
|
4 |
3.2 |
1800' |
54 |
6234' |
3,74 |
1530' |
0.492 |
0.049 |
0.051 |
0.071 |
|
4 |
2.8 |
2048' |
46 |
5434' |
3,23 |
1800' |
0.539 |
0.075 |
0.0644 |
0.09 |
|
-2 |
2.95 |
1945' |
49 |
5734' |
3,40 |
1705' |
0.524 |
0.069 |
0.063 |
0.088 |
|
-2 |
3.1 |
1842' |
52 |
6034' |
3,60 |
1606' |
0.502 |
0.054 |
0.0536 |
0.075 |
На конце профиля:
(-1) |
|||||||||||
0 |
3,34 |
1724' |
56 |
86 |
6,188 |
918' |
0.033 |
0.0019 |
0.0561 |
0.078 |
|
4 |
3,74 |
1530' |
62 |
92 |
7,184 |
800' |
0.021 |
0.001 |
0.048 |
0.067 |
|
4 |
3,23 |
1800' |
54 |
84 |
5,875 |
948' |
0.037 |
0.0028 |
0.0722 |
0.101 |
|
-2 |
3,40 |
1705' |
57 |
87 |
6.321 |
906' |
0.03 |
0.0021 |
0.0669 |
0.093 |
|
-2 |
3,60 |
1606' |
60 |
90 |
6,845 |
824' |
0.025 |
0.0013 |
0.0527 |
0.074 |
2 Вычисление нормальной и продольной сил для профиля при рассматриваемых углах атаки
;
при = 0
Так как профиль симметричен и находится под нулевым углом атаки то подьемная сила.Осевая же сила
при = 4
Находим подъемную силу
Осевая сила
при = -2
Находим подъемную силу
Осевая сила
3Найдем коэффициенты для нормальной и осевой сил для профиля при рассматриваемых углах атаки
;
при = 0 ;
при = 4 ;
при = -2 ;
4 Найдем моменты тангажа для профиля при рассматриваемых углах атаки
при = 0
при =4
при =-2
5Найдем координаты относительного центра давлений при рассматриваемых углах атаки.
при =4
при =-2
6 Вычисляем коэффициенты лобового сопротивления и подъемной силы при рассматриваемых углах атаки.
при = 0
при = 4
при = -2
7. Оформление чертежа
Вычертить картину обтекания профиля при заданных углах атаки . Для каждого участка течения указать значения М, Р, Р0.
Ниже вычертить эпюру распределения по профилю избыточных давлений с указанием их величины. Рядом указать значение нормальной и продольной сил, их аэродинамических коэффициентов, а также указать графически координату центра давления.
Расчет крыльевого профиля в среде Gas2
Структурная схема:
Вначале строим внешний контур в виде квадрата, внутри которого строится контур крыльевого профиля. Далее мы разбиваем область вокруг крыльевого профиля на ячейки с помощью ступенчатой аппроксимации. На внешнем контуре задаем границы втекания и вытекания потока.
Система координат берем плоскую.
На границе втекания задаем скорость набегающего потока и его температуру. А для разных углов атаки профиля скорость набегающего потока раскладываем на соответствующие оси Vx и Vy. Кроме того в расчетах учитываем вязкость.
1. Распределение давления при г=0?.
2.Распределение давления при г=4?.
сопло газодинамический крыльевой уплотнение
3. Распределение давления при г=-2?.
Вычисление подъемной и продольной сил для профиля при рассматриваемых углах атаки
;
при = 0
Так как профиль симметричен и находится под нулевым углом атаки то подъемная сила.Осевая же сила
при = 4
Находим подъемную силу
Осевая сила
при = -2
Находим подъемную силу
Осевая сила
3.Найдем коэффициенты для нормальной и осевой сил для профиля при рассматриваемых углах атаки
;
при = 0 ;
при = 4 ;
при = -2 ;
4 Найдем моменты тангажа для профиля при рассматриваемых углах атаки
при = 0
при =4
при =-2
5Найдем координаты относительного центра давлений при рассматриваемых углах атаки.
при =4
при =-2
6. Вычисляем коэффициенты лобового сопротивления и подъемной силы при рассматриваемых углах атаки.
при = 0
при = 4
при = -2
Подводя анализ полученного результата можно отметить, что избыточное давление полученное в среде Gas2 и с помощью метода скачков уплотнения и волн разрежения в целом схожи, но заметно различаются на конце профиля. Подъемная сила для углов атаки отличных от нуля рассчитанная в Gas2 получилась несколько больше чем при использовании метода скачков уплотнения и волн разряжения.
Относительная погрешность для подъемной и осевой силах:
При = 0
При = 4
При = -2
Список литературы
Абрамович Г.Н. Прикладная газовая динамика. М. ,Наука ,.1976г. 824 с.
Алемасов В.Е. Теория ракетных двигателей М. Машиностроение 1980г. 533с.
Карафоли Е. Аэродинамика больших скоростей M изд. Академия наук СССР 1960г. 740c.
Краснов Н.Ф. Аэродинамика в вопросах и задачах М. Высшая школа 1985г.759c.
Лойцянский Л.Г. Механика жидкости и газа M.,Наука,1973г.847c.
Размещено на Allbest.ru
Подобные документы
Сопло Лаваля как техническое приспособление, служащее для ускорения газового потока. Рассмотрение основных особенностей построения графика газодинамических функций давления, скорости. Этапы расчета параметров течения воздушного потока в сопле Лаваля.
контрольная работа [394,1 K], добавлен 10.01.2013Задачи ориентировочного расчета паровой турбины. Определение числа ступеней, их диаметров и распределения тепловых перепадов по ступеням. Вычисление газодинамических характеристик турбины, выбор профиля сопловой лопатки, определение расхода пара.
курсовая работа [840,0 K], добавлен 11.11.2013Роль одномерного анализа при решении технических задач. Уравнения Бернулли для идеальной и реальной жидкостей. Выражение скорости звука через термодинамические параметры. Изоэнтропийное течение, критический расход. Сопло Лаваля и принцип его действия.
реферат [962,8 K], добавлен 07.01.2014Расчет технологической системы, потребляющей тепловую энергию. Расчет параметров газа, определение объемного расхода. Основные технические параметры теплоутилизаторов, определение количества выработанного конденсата, подбор вспомогательного оборудования.
курсовая работа [112,1 K], добавлен 20.06.2010Изучение механизма работы человеческого уха. Определение понятия и физических параметров звука. Распространение звуковых волн в воздушной среде. Формула расчета скорости звука. Рассмотрение числа Маха как характеристики безразмерной скорости течения газа.
реферат [760,2 K], добавлен 18.04.2012Расчет профиля диффузии сурьмы в кремнии, определение основных параметров этого процесса. Использование феноменологической модели диффузии. Влияние параметров на глубину залегания примеси. Численное решение уравнения диффузии по неявной разностной схеме.
курсовая работа [4,7 M], добавлен 28.08.2010Определение состава газа, расчет горения топлива. Расчет нагрева металла. Основные параметры, тепловой баланс, основные размеры печи, выбор горелок. Расчет рекуператора, гидравлический расчет трассы. Применение бетонов и волокнистой изоляции.
курсовая работа [1,4 M], добавлен 22.10.2010Параметры упругих гармонических волн. Уравнения плоской и сферической волн. Уравнение стоячей волны. Распространение волн в однородной изотропной среде и принцип суперпозиции. Интервалы между соседними пучностями. Скорость распространения звука.
презентация [155,9 K], добавлен 18.04.2013Начальные параметры ударной волны, образующейся движением пластины. Параметры воздуха на фронте ударной волны в момент подхода волны к преграде. Расчет параметров продуктов детонации в начальный момент отражения от жесткой стенки и металлической пластины.
курсовая работа [434,5 K], добавлен 20.09.2011Расчет конструктивных и технологических параметров поперечно-струйной турбины, водоводов и водоприемника. Количество вырабатываемой электроэнергии за год и объем плотины для гидроэлектростанции, работающей при расходе воды Qн=0,8м/c2 и напоре сопла Нс=6м.
контрольная работа [1,4 M], добавлен 24.09.2013