Список использованных источников
1 Безопасность жизнедеятельности и промышленная безопасность: Учебное пособие / Под ред. проф. В.Д. Шантарина - Тюмень: ТюмГНГУ, 2001. - 308с.
2 Вестник" Самарского государственного аэрокосмического университета имени академика С.П. Королева, №2 - Материалы IV Всероссийской научно-практической конференции "Процессы горения, теплообмена и экология тепловых двигателей". Самара: СГАУ, 2002. - 111с.
3 Зенкевич О. Метод конечных элементов в технике: Пер. с англ. - М.: Мир, 1975. - 541с., ил.
4 Кашапов Р.С., Максимов Д.А., Скиба Д.В., Куликов С.В., Баштанников М.Н. Исследование автоколебаний давления в камере сгорания с предварительным смешением топлива // Газотурбинные технологии №4 (13). Рыбинск, 2001. С.34-37.
5 Латыпов Р.Ш. Вопросы рациональной эксплуатации газотурбинных установок: Учебное пособие. Уфа: УГНТУ, 2000. - 100с.
6 Методические указания по дипломному проектированию для студентов очного и заочного обучения специальности 220201 - Управление и информатика в технических системах. - Тюмень: ТюмГНГУ, 2006. - 28с.
7 Методические указания к оценке экономической эффективности автоматизированных систем в курсовом и дипломном проектировании для студентов специальности АСОиУ, АТП, УИТС дневного и заочного обучения - Тюмень: ТюмГНГУ, 2003. - 33с.
8 "Моделирование процессов горения в ANSYS CFX. Хитрых Д., "ANSYS Solutions. Русская редакция" № 3 (май) '2006.
9 Норри Д., де Фриз Ж. Введение в метод конечных элементов: Пер. с англ. - М.: Мир, 1981. - 304с., ил.
|
|
10 Справочная система CAE-пакета ANSYS.
11 Старикова Г.В., Столбова Н.В., Дорофеева Э.С. Охрана окружающей природной среды: Учебное пособие. - Тюмень: ТюмГНГУ, 2000. - 94с.
12 Чигарев А.В., Кравчук А.С., Смалюк А.Ф. ANSYS для инженеров: Справ. пособие. М.: Машиностроение-1, 2004, 512с.
13 www.ansys. spb.ru Санкт-Петербургский сайт лаборатории "Вычислительная механика", СПбГПУ, 2004г.
14 www.aviaport.ru ГТУ - основа будущей энергетики России: электронная статья, Олег Фаворский, 2005г.
15 www.cae. tgngu. tyumen.ru Компьютерный инженерный анализ. Современные технологии проектирования.
16 www.college.ru Электронный справочник по физике газов: справочные таблицы.
17 www.engin.ru Официальный сайт Тюменской ТЭЦ-1
18 www.esco-ecosys. narod.ru. Электронный журнал энергосервисной компании "Экологические системы", №1, январь 2006г.
19 www.gtt.ru Электронный журнал "Газотурбинные технологии", 2004г.
20 www.siemens.ru Официальный русскоязычный сайт компании "Siemens".
Приложения
Приложение А
Листинг программы
!***************************************************************************
!Программа газодинамического расчета камеры сгорания газовой турбины v64.3A *
!Автор: студент 5-го курса УИТС, Панышев Евгений Александрович *
|
|
!***************************************************************************
!Начало
!Импортирование модели формата IGES в ANSYS
/AUX15
IOPTN,IGES,NODEFEAT
IOPTN,MERGE,YES
IOPTN,SOLID,YES
IOPTN,SMALL,YES
IOPTN,GTOLER, DEFA
IGESIN,'chambre','IGS','.\GrafModel\'
LPLOT
/REPLOT,RESIZE
/AUTO, 1
/REP
/VIEW, 1 ,,,1
/ANG, 1
/REP,FAST
/REPLOT,RESIZE
FINISH
!Выбор единиц измерения
/units,si0020
!Задание констант моделирования
/PREP7
vvoz=1.54 !Скорость воздуха
tvoz=450 !Температура воздуха
vgas=1.57 !Скорость газа
tgas=350 !Температура газа
tflame=(vvoz*tvoz+vgas*tgas)/vgas-vvoz*vvoz*15
tsten=tflame*0.6 !Температура стенок(т.к. не учитываем излучения пламени)
!Масштабирование графической модели
FLST,2,28,4,ORDE,2
FITEM,2,1
FITEM,2,-28
LSSCALE,P51X, , ,0.01,0.01,0.01, ,1,1
!Сжатие номеров
NUMCMP,KP
NUMCMP,LINE
!Постройка площадей
FLST,2,26,4
FITEM,2,21 !Выбор линий первой (рабочей) плоскости
FITEM,2,20
FITEM,2,19
FITEM,2,18
FITEM,2,17
FITEM,2,10
FITEM,2,9
FITEM,2,8
FITEM,2,28
FITEM,2,4
FITEM,2,1
FITEM,2,3
FITEM,2,27
FITEM,2,5
FITEM,2,6
FITEM,2,7
FITEM,2,11
FITEM,2,12
FITEM,2,13
FITEM,2,14
FITEM,2,15
FITEM,2,24
FITEM,2,23
FITEM,2,25
FITEM,2,26
FITEM,2,22
AL,P51X !Постройка рабочей плоскости
FLST,2,6,4
FITEM,2,24 !Выбор линий второй плоскости
FITEM,2,23
FITEM,2,25
FITEM,2,16
FITEM,2,26
FITEM,2,22
|
|
AL,P51X !Постройка второй плоскости
!Подключение фильтров графического меню
/NOPR
/PMETH,OFF,0
KEYW,PR_SET,1
KEYW,PR_STRUC,0
KEYW,PR_THERM,1
KEYW,PR_FLUID,1
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0
KEYW,MAGEDG,0
KEYW,MAGHFE,0
KEYW,MAGELC,0
KEYW,PR_MULTI,2
KEYW,PR_CFD,1
/GO
!Выбор типа элементов
ET,1,FLUID141
!!!!!!!!Разбитие на элементы рабочей площади
!Разметка линий
FLST,2,26,4,ORDE,5
FITEM,2,1
FITEM,2,3
FITEM,2,-15
FITEM,2,17
FITEM,2,-28
LESIZE,P51X, , ,-1, ,1
FLST,5,15,4,ORDE,4
FITEM,5,11
FITEM,5,-15
FITEM,5,17
FITEM,5,-26
LSEL,S, , ,P51X
FLST,5,15,4,ORDE,4
FITEM,5,11
FITEM,5,-15
FITEM,5,17
FITEM,5,-26
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,0.03, , , , , , ,1 !Разбивка линий по расстоянию между узлами (0.03)
FLST,5,6,4,ORDE,2
FITEM,5,5
FITEM,5,-10
LSEL,S, , ,P51X
FLST,5,6,4,ORDE,2
FITEM,5,5
FITEM,5,-10
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,0.04, , , , , , ,1 !Разбивка линий по расстоянию между узлами (0.04)
FLST,5,5,4,ORDE,5
FITEM,5,1
FITEM,5,3
FITEM,5,-4
FITEM,5,27
FITEM,5,-28
LSEL,S, , ,P51X
FLST,5,5,4,ORDE,5
FITEM,5,1
FITEM,5,3
FITEM,5,-4
FITEM,5,27
FITEM,5,-28
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,0.03, , , , , , ,1 !Разбивка линий по расстоянию между узлами (0.03)
!Выбор вида элементов (треугольные)
MSHAPE,1,2D
MSHKEY,0
CM,_Y,AREA
ASEL, , , , 1 !Выбор площади для разбиения
CM,_Y1,AREA
CMSEL,S,_Y
AMESH,_Y1 !Наложение конечно-элементной сетки на выбранную площадь
|
|
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!!!!!!!!!!!!!!!!Задание начальных условий и нагружений
! Задание параметров стенок КС
/UI,MESH,OFF
ALLSEL,ALL
FLST,5,22,4,ORDE,6 !Выбор линий
FITEM,5,3
FITEM,5,-14
FITEM,5,18
FITEM,5,-23
FITEM,5,25
FITEM,5,-28
LSEL,S, , ,P51X
NSEL,ALL
NSLL,S,1
FLST,5,239,1,ORDE,10
FITEM,5,1
FITEM,5,-100
FITEM,5,112
FITEM,5,-218
FITEM,5,224
FITEM,5,-233
FITEM,5,240
FITEM,5,-249
FITEM,5,255
FITEM,5,-266
CM,_Y,NODE
NSEL,R, , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
CMDELE,_Y
D,_Y1,VX,0 !Задние скорости на границы КС
D,_Y1,VY,0
d,_Y1,temp,tsten !Задание температуры стенок КС
CMDELE,_Y1
! Задание скорости течения воздуха на левый патрубок
LSEL,S, , , 17
NSEL,ALL
NSLL,S,1
FLST,5,7,1,ORDE,3
FITEM,5,245
FITEM,5,249
FITEM,5,-254
CM,_Y,NODE
NSEL,R, , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
CMDELE,_Y
D,_Y1,VX,vvoz !Задание нагрузки
D,_Y1,VY,0
CMDELE,_Y1
!Задание температуры потока на левый патрубок
FLST,2,7,1,ORDE,3
FITEM,2,245
FITEM,2,249
FITEM,2,-254
/GO
D,P51X,TEMP,tvoz !Задание нагрузки
!Задание скорости потока на правый патрубок
LSEL,S, , , 15 !Выбор линии
NSEL,ALL !Выбор всех узлов на линии
NSLL,S,1
FLST,5,7,1,ORDE,3
FITEM,5,212
FITEM,5,218
FITEM,5,-223
CM,_Y,NODE
NSEL,R, , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
CMDELE,_Y
D,_Y1,VX,vvoz*sin(0.2618) !Задание нагрузки
D,_Y1,VY,-vvoz*sin(1.309) !в соответствии с глобальной системой координат
CMDELE,_Y1
!Задание темепературы потока на правый патрубок
FLST,2,7,1,ORDE,3
FITEM,2,212
FITEM,2,218
FITEM,2,-223
/GO
D,P51X,TEMP,tvoz !Задание значения нагрузки
!Задание скорости потока газа
LSEL,S, , , 24 !Выбор линий
NSEL,ALL !Выбор всех узлов линии
NSLL,S,1
FLST,5,8,1,ORDE,3
FITEM,5,228
FITEM,5,233
FITEM,5,-239
CM,_Y,NODE
NSEL,R, , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
CMDELE,_Y
D,_Y1,VX,vgas*sin(0.8727) !Задание значения нагрузки
D,_Y1,VY,-vgas*sin(0.6981) !в соответствии с глобальной системой координат
CMDELE,_Y1
!Задание температуры потока газа
FLST,2,8,1,ORDE,3
FITEM,2,228
FITEM,2,233
FITEM,2,-239
/GO
D,P51X,TEMP,tgas !Задание значения нагрузки
!Задание граничного условия на выход КС
LSEL,S, , , 1 !Выбор линии
NSEL,ALL !Выбол всех узлов линии
NSLL,S,1
FLST,5,13,1,ORDE,3
FITEM,5,91
FITEM,5,100
FITEM,5,-111
CM,_Y,NODE
NSEL,R, , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
CMDELE,_Y
/GO
D,_Y1,PRES,0 !Задание значения нагрузки
CMDELE,_Y1
SAVE !Сохранение конечно-элементной модели с нагружениями
!Создание поля температур для эмуляции пламени
FLST,5,1032,1,ORDE,84 !Выбор элементов поля
FITEM,5,443
FITEM,5,453
FITEM,5,459
FITEM,5,468
FITEM,5,-470
FITEM,5,475
FITEM,5,480
FITEM,5,697
FITEM,5,700
FITEM,5,-703
FITEM,5,705
FITEM,5,-706
FITEM,5,708
FITEM,5,721
FITEM,5,725
FITEM,5,-732
FITEM,5,925
FITEM,5,928
FITEM,5,-931
FITEM,5,933
FITEM,5,946
FITEM,5,949
FITEM,5,-958
FITEM,5,1142
FITEM,5,-1150
FITEM,5,1161
FITEM,5,-1162
FITEM,5,1165
FITEM,5,-1172
FITEM,5,1341
FITEM,5,1344
FITEM,5,-1345
FITEM,5,1348
FITEM,5,-1356
FITEM,5,1365
FITEM,5,-1367
FITEM,5,1369
FITEM,5,-1374
FITEM,5,1539
FITEM,5,-1540
FITEM,5,1543
FITEM,5,-1550
FITEM,5,1566
FITEM,5,-1573
FITEM,5,1645
FITEM,5,-1654
FITEM,5,1725
FITEM,5,1728
FITEM,5,-1735
FITEM,5,1749
FITEM,5,-1758
FITEM,5,1813
FITEM,5,-1836
FITEM,5,1896
FITEM,5,-1903
FITEM,5,1911
FITEM,5,-1921
FITEM,5,1943
FITEM,5,-1944
FITEM,5,1949
FITEM,5,-1951
FITEM,5,1957
FITEM,5,-1985
FITEM,5,2039
FITEM,5,-2045
FITEM,5,2050
FITEM,5,2053
FITEM,5,-2062
FITEM,5,2075
FITEM,5,-2117
FITEM,5,2123
FITEM,5,-2144
FITEM,5,2163
FITEM,5,-2172
FITEM,5,2175
FITEM,5,2177
FITEM,5,-2184
FITEM,5,2187
FITEM,5,-2266
FITEM,5,2268
FITEM,5,-2269
FITEM,5,2271
FITEM,5,2274
FITEM,5,-2914
NSEL,S, , ,P51X
FLST,2,1032,1,ORDE,84
FITEM,2,443
FITEM,2,453
FITEM,2,459
FITEM,2,468
FITEM,2,-470
FITEM,2,475
FITEM,2,480
FITEM,2,697
FITEM,2,700
FITEM,2,-703
FITEM,2,705
FITEM,2,-706
FITEM,2,708
FITEM,2,721
FITEM,2,725
FITEM,2,-732
FITEM,2,925
FITEM,2,928
FITEM,2,-931
FITEM,2,933
FITEM,2,946
FITEM,2,949
FITEM,2,-958
FITEM,2,1142
FITEM,2,-1150
FITEM,2,1161
FITEM,2,-1162
FITEM,2,1165
FITEM,2,-1172
FITEM,2,1341
FITEM,2,1344
FITEM,2,-1345
FITEM,2,1348
FITEM,2,-1356
FITEM,2,1365
FITEM,2,-1367
FITEM,2,1369
FITEM,2,-1374
FITEM,2,1539
FITEM,2,-1540
FITEM,2,1543
FITEM,2,-1550
FITEM,2,1566
FITEM,2,-1573
FITEM,2,1645
FITEM,2,-1654
FITEM,2,1725
FITEM,2,1728
FITEM,2,-1735
FITEM,2,1749
FITEM,2,-1758
FITEM,2,1813
FITEM,2,-1836
FITEM,2,1896
FITEM,2,-1903
FITEM,2,1911
FITEM,2,-1921
FITEM,2,1943
FITEM,2,-1944
FITEM,2,1949
FITEM,2,-1951
FITEM,2,1957
FITEM,2,-1985
FITEM,2,2039
FITEM,2,-2045
FITEM,2,2050
FITEM,2,2053
FITEM,2,-2062
FITEM,2,2075
FITEM,2,-2117
FITEM,2,2123
FITEM,2,-2144
FITEM,2,2163
FITEM,2,-2172
FITEM,2,2175
FITEM,2,2177
FITEM,2,-2184
FITEM,2,2187
FITEM,2,-2266
FITEM,2,2268
FITEM,2,-2269
FITEM,2,2271
FITEM,2,2274
FITEM,2,-2914
/GO
D,P51X,TEMP,tflame !Задание нагрузки
!Задние поля обратного потока в пламени
FLST,5,255,1,ORDE,106 !Выбор элементов
FITEM,5,728
FITEM,5,-730
FITEM,5,925
FITEM,5,931
FITEM,5,946
FITEM,5,951
FITEM,5,-956
FITEM,5,1144
FITEM,5,-1149
FITEM,5,1165
FITEM,5,1167
FITEM,5,-1170
FITEM,5,1349
FITEM,5,-1354
FITEM,5,1370
FITEM,5,-1374
FITEM,5,1545
FITEM,5,-1550
FITEM,5,1571
FITEM,5,-1573
FITEM,5,1730
FITEM,5,-1735
FITEM,5,1755
FITEM,5,-1758
FITEM,5,1896
FITEM,5,-1903
FITEM,5,1912
FITEM,5,1919
FITEM,5,-1921
FITEM,5,2040
FITEM,5,-2045
FITEM,5,2054
FITEM,5,2060
FITEM,5,-2062
FITEM,5,2163
FITEM,5,2165
FITEM,5,-2170
FITEM,5,2178
FITEM,5,2183
FITEM,5,-2184
FITEM,5,2279
FITEM,5,2281
FITEM,5,-2285
FITEM,5,2290
FITEM,5,2295
FITEM,5,-2298
FITEM,5,2387
FITEM,5,-2390
FITEM,5,2394
FITEM,5,-2396
FITEM,5,2398
FITEM,5,-2402
FITEM,5,2482
FITEM,5,2484
FITEM,5,-2487
FITEM,5,2492
FITEM,5,-2493
FITEM,5,2496
FITEM,5,-2499
FITEM,5,2507
FITEM,5,-2510
FITEM,5,2576
FITEM,5,2578
FITEM,5,-2582
FITEM,5,2584
FITEM,5,-2589
FITEM,5,2594
FITEM,5,-2600
FITEM,5,2624
FITEM,5,-2627
FITEM,5,2658
FITEM,5,-2667
FITEM,5,2671
FITEM,5,-2678
FITEM,5,2698
FITEM,5,-2702
FITEM,5,2728
FITEM,5,-2730
FITEM,5,2733
FITEM,5,-2736
FITEM,5,2739
FITEM,5,-2747
FITEM,5,2764
FITEM,5,-2769
FITEM,5,2790
FITEM,5,-2795
FITEM,5,2798
FITEM,5,-2806
FITEM,5,2820
FITEM,5,-2825
FITEM,5,2841
FITEM,5,-2843
FITEM,5,2846
FITEM,5,-2854
FITEM,5,2862
FITEM,5,-2868
FITEM,5,2878
FITEM,5,2880
FITEM,5,-2881
FITEM,5,2884
FITEM,5,-2891
FITEM,5,2894
FITEM,5,-2899
FITEM,5,2906
FITEM,5,2909
FITEM,5,-2914
NSEL,S, , ,P51X
FLST,5,255,1,ORDE,106
FITEM,5,728
FITEM,5,-730
FITEM,5,925
FITEM,5,931
FITEM,5,946
FITEM,5,951
FITEM,5,-956
FITEM,5,1144
FITEM,5,-1149
FITEM,5,1165
FITEM,5,1167
FITEM,5,-1170
FITEM,5,1349
FITEM,5,-1354
FITEM,5,1370
FITEM,5,-1374
FITEM,5,1545
FITEM,5,-1550
FITEM,5,1571
FITEM,5,-1573
FITEM,5,1730
FITEM,5,-1735
FITEM,5,1755
FITEM,5,-1758
FITEM,5,1896
FITEM,5,-1903
FITEM,5,1912
FITEM,5,1919
FITEM,5,-1921
FITEM,5,2040
FITEM,5,-2045
FITEM,5,2054
FITEM,5,2060
FITEM,5,-2062
FITEM,5,2163
FITEM,5,2165
FITEM,5,-2170
FITEM,5,2178
FITEM,5,2183
FITEM,5,-2184
FITEM,5,2279
FITEM,5,2281
FITEM,5,-2285
FITEM,5,2290
FITEM,5,2295
FITEM,5,-2298
FITEM,5,2387
FITEM,5,-2390
FITEM,5,2394
FITEM,5,-2396
FITEM,5,2398
FITEM,5,-2402
FITEM,5,2482
FITEM,5,2484
FITEM,5,-2487
FITEM,5,2492
FITEM,5,-2493
FITEM,5,2496
FITEM,5,-2499
FITEM,5,2507
FITEM,5,-2510
FITEM,5,2576
FITEM,5,2578
FITEM,5,-2582
FITEM,5,2584
FITEM,5,-2589
FITEM,5,2594
FITEM,5,-2600
FITEM,5,2624
FITEM,5,-2627
FITEM,5,2658
FITEM,5,-2667
FITEM,5,2671
FITEM,5,-2678
FITEM,5,2698
FITEM,5,-2702
FITEM,5,2728
FITEM,5,-2730
FITEM,5,2733
FITEM,5,-2736
FITEM,5,2739
FITEM,5,-2747
FITEM,5,2764
FITEM,5,-2769
FITEM,5,2790
FITEM,5,-2795
FITEM,5,2798
FITEM,5,-2806
FITEM,5,2820
FITEM,5,-2825
FITEM,5,2841
FITEM,5,-2843
FITEM,5,2846
FITEM,5,-2854
FITEM,5,2862
FITEM,5,-2868
FITEM,5,2878
FITEM,5,2880
FITEM,5,-2881
FITEM,5,2884
FITEM,5,-2891
FITEM,5,2894
FITEM,5,-2899
FITEM,5,2906
FITEM,5,2909
FITEM,5,-2914
CM,_Y,NODE
NSEL,R, , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
CMDELE,_Y
D,_Y1,VX,-(0.01*exp(tflame/1000)-0.04)*sin(0.8727) !Скорость противотока газов
D,_Y1,VY,(0.01*exp(tflame/1000)-0.04)*sin(0.6981)
D,_Y1,ENKE,-1
CMDELE,_Y1
EPLOT
SAVE !Сохранение конечно-элементной модели со всеми нагружениями
!Настройка параметров решателя
FLDATA1,SOLU,TRAN,0 !Используем стационарный анализ
FLDATA1,SOLU,FLOW,1 !Используем расчет потоков
FLDATA1,SOLU,TEMP,0 !Не используем температурный расчет
FLDATA1,SOLU,TURB,1 !Используем модель турбулентности (по умолчанию)
FLDATA1,SOLU,COMP,0
FLDATA1,SOLU,VOF,0
FLDATA1,SOLU,SFTS,0
FLDATA1,SOLU,IVSH,0
FLDATA1,SOLU,SWRL,0
FLDATA1,SOLU,SPEC,0
FLDATA1,SOLU,ALE,0
alls !В расчете участвуют все элементы
/prep7
flda,conv,iter,5 !Задание параметров сходимости
FLDATA21,CONV,PRES,1e-008,
!Задание параметров среды
/solu
flda,iter,exec,50 !Необходимо несколько итераций для начального расчета
flda,temp,nomi,350 !Начальная температура 350K
flda,prot,dens,air-si !Начальное решение для воздуха
flda,vary,dens,true
flda,prot,visc,air-si
flda,vary,visc,true
flda,prot,cond,air-si
flda,vary,cond,true
flda,prot,spht,air-si
save
solve !Провести 50 итераций для воздуха при 350K
finish !Подготовка к многокомпонентному переносу
!Задание параметров многокомпонентного переноса
/prep7
flda,prot,dens,cmix !Растворенная смесь для плотности
flda,prot,visc,cmix !Растворенная смесь для вязкости
keyopt,1,1,2 !Установка двух компонент смеси
flda,solu,spec,t !Включение опции многокомпонентного переноса
msdata,1 !Задание компонента 1 (GAS) в качестве алгебраического
!!!!!!! Задание свойств каждого из компонентов
!Для топливного газа
msspec,1,gas,16.04 !Молекулярная масса
msprop,1,DENS,GAS,0.7168,300,1.01325E+5 !Плотность
msvary,1,dens,t !Варьирование плотности
msnomf,1,.3 !Массовая доля вещества
msprop,1,VISC,GAS,3.786E-4,300,1.01325E+5 !Вязкость
msprop,1,mdif,GAS,1.601E-6 !Массовая диффузия
msvary,1,visc,t !Варьирование вязкости
msprop,1,cond,GAS,.01598 !Теплопроводность
mscap,1,1
msspec,2,voz,28.96 !Молекулярная масса
msprop,2,DENS,GAS,1.2928,300,1.01325E+5 !Плотность
msvary,2,dens,t !Варьирование плотности
msnomf,2,.7 !Массовая доля вещества
msprop,2,VISC,GAS,2.2067E-5 !Вязкость
msprop,2,mdif,GAS,2.149E-5 !Массовая диффузия
msvary,2,visc,t !Варьирование вязкости
msprop,2,cond,GAS,.02674 !Теплопроводность
mscap,2,1
!Степень стабилизации параметров
MSMIR,1,1,
MSMIR,2,1,
!Настройка решателя многокомпонентного переноса вещества
msrelx,2,1.0
msmeth,2,3
mssolu,2,,,2,1.e-9
msmeth,1,3
mssolu,1,,,2,1.e-9
alls
!!!!!!! Установка граничных условий смесей
lsel,s,,,15,17
nsll,,1
d,all,gas,0.0 !Для подвода воздуха
d,all,voz,1.0
lsel,s,,,24
nsll,,1
d,all,gas,1.0 !Для подвода топливного газа
d,all,voz,0.0
alls
!!!!!!! конец установки граничных условий
/prep7
flda,iter,exec,100 !Необходимо выполнить 100 глобальных итераций
flda,conv,outp,land !Настроить монитор сходимости
save
finish
/solu
solve !Провести 100 итераций с активированными компонентами
!Подготовка к решению уравнений энергии
flda,solu,temp,t !Включить решатель уравнения энергии вещества
flda,solu,flow,f !заморозить поле потока
flda,meth,temp,3 !Активировать решатель PCCR
flda,conv,temp,1.e-8 !Критерий сходимости для PCCR
flda,iter,exec,50 !Для решения необходимо несколько итераций
flda,relx,temp,1.0 !Отсутствует релаксация на температуру
save
solve !Решение уравнения энергии
!! Подготовка совместного решения
flda,iter,exec,50 !Число итераций на проведение совместного решения
flda,solu,flow,t !Разморозка поля потоков
solve !Совместное решение
flda,iter,exec,100 !Необходимо выполнить 100 дополнительных итераций
!для улучшения результатов расчета
save
solve !Решение
FINISH
!!Улучшение модели турбулентности и задание её параметров
/solu
FLDATA24,TURB,MODL,4 !Выбор модели турбулентности
FLDATA24,NKET,C1MX,0.43, !Настройка параметров
FLDATA24,NKET,C2,1.9,
FLDATA24,NKET,SCTK,1,
FLDATA24,NKET,SCTD,1.2,
save
solve !Решение с новой моделью турбулентности
finish
!!!!!!!Включение опции сжимаемого потока
!Расчет уравнения энергии
/SOLU
FLDATA1,SOLU,TRAN,0
FLDATA1,SOLU,FLOW,0 !Отключение опции переноса вещества
FLDATA1,SOLU,TEMP,1
FLDATA1,SOLU,TURB,1
FLDATA1,SOLU,COMP,1 !Сжимаемый поток
FLDATA1,SOLU,VOF,0
FLDATA1,SOLU,SFTS,0
FLDATA1,SOLU,IVSH,0
FLDATA1,SOLU,SWRL,0
FLDATA1,SOLU,SPEC,1
FLDATA1,SOLU,ALE,0
!Совместное решение
SOLVE
FLDATA1,SOLU,TRAN,0
FLDATA1,SOLU,FLOW,1 !Включение опции переноса вещества
FLDATA1,SOLU,TEMP,1
FLDATA1,SOLU,TURB,1
FLDATA1,SOLU,COMP,1 !Сжимаемый поток
FLDATA1,SOLU,VOF,0
FLDATA1,SOLU,SFTS,0
FLDATA1,SOLU,IVSH,0
FLDATA1,SOLU,SWRL,0
FLDATA1,SOLU,SPEC,1
FLDATA1,SOLU,ALE,0
SOLVE
FINISH
!************************************************************************
! Расчет модели камеры сгорания газовой турбины закончен *
!************************************************************************
!Подготовка к просмотру результатов решения
/POST1
SET,LAST
!Просмотр результатов ведется в интерактивном режиме
!Конец
Приложение Б
Диаграммы результатов экспериментов
Рисунок Б.4 - Диаграмма поля распределения температур в камере сгорания газовой турбины V64.3A (a = 8)
Рисунок Б.6 - Диаграмма векторного поля распределения скоростей (a = 8)
Рисунок Б.7 - Диаграмма поля распределения температур в камере сгорания газовой турбины V64.3A (a = 10)
Рисунок Б.8 - Диаграмма кинетической энергии турбулентности (a = 10)
Рисунок Б.9 - Диаграмма векторного поля распределения скоростей (a = 10)
Рисунок Б.10 - Диаграмма поля распределения температур в камере сгорания газовой турбины V64.3A (a = 4)
Рисунок Б.11 - Диаграмма кинетической энергии турбулентности (a = 4)
Рисунок Б.12 - Диаграмма векторного поля распределения скоростей (a = 4)
Рисунок Б.13 - Диаграмма поля распределения температур в камере сгорания газовой турбины V64.3A (a = 2)
Рисунок Б.15 - Диаграмма векторного поля распределения скоростей (a = 2)
Размещено на Allbest.ru
Дата добавления: 2018-09-22; просмотров: 213; Мы поможем в написании вашей работы! |
Мы поможем в написании ваших работ!