RU2678155C1 - Способ прогнозирования помпажа в газовом компрессоре - Google Patents
Способ прогнозирования помпажа в газовом компрессоре Download PDFInfo
- Publication number
- RU2678155C1 RU2678155C1 RU2018105189A RU2018105189A RU2678155C1 RU 2678155 C1 RU2678155 C1 RU 2678155C1 RU 2018105189 A RU2018105189 A RU 2018105189A RU 2018105189 A RU2018105189 A RU 2018105189A RU 2678155 C1 RU2678155 C1 RU 2678155C1
- Authority
- RU
- Russia
- Prior art keywords
- compressor
- computer
- outlet
- inlet
- implemented method
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 238000012795 verification Methods 0.000 claims abstract description 4
- 239000012530 fluid Substances 0.000 claims description 9
- 238000005094 computer simulation Methods 0.000 claims description 8
- 239000000126 substance Substances 0.000 abstract 1
- 238000002485 combustion reaction Methods 0.000 description 13
- 238000011144 upstream manufacturing Methods 0.000 description 5
- 239000000446 fuel Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D27/00—Control, e.g. regulation, of pumps, pumping installations or pumping systems specially adapted for elastic fluids
- F04D27/001—Testing thereof; Determination or simulation of flow characteristics; Stall or surge detection, e.g. condition monitoring
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05D—INDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
- F05D2260/00—Function
- F05D2260/81—Modelling or simulation
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05D—INDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
- F05D2260/00—Function
- F05D2260/82—Forecasts
- F05D2260/821—Parameter estimation or prediction
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05D—INDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
- F05D2270/00—Control
- F05D2270/01—Purpose of the control system
- F05D2270/10—Purpose of the control system to cope with, or avoid, compressor flow instabilities
- F05D2270/101—Compressor surge or stall
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mechanical Engineering (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Control Of Positive-Displacement Air Blowers (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
Изобретение относится к способу прогнозирования точки помпажа компрессора. Технический результат заключается в автоматизации прогнозирования помпажа в рабочей характеристике газового компрессора посредством расчета CFD. Способ содержит этапы: генерирование множества сеток, сборка множества сеток, задание граничных условий вычислительной области, задание условий атмосферного давления на концевой площади выходного сечения сопла, вычисление потока вдоль компрессора и вычисление массового расхода на входе в компрессор и коэффициента давления компрессора, проверка достижения предопределенного предела численной стабильности, и если предел численной стабильности не достигнут, то: уменьшение размеров концевой площади выходного сечения сопла, повторное генерирование сетки выпуска компрессора, повторение этапов задания атмосферного давления на концевой площади выходного сечения сопла, этапа вычисления и упомянутого этапа проверки, если предел численной стабильности достигнут, то способ останавливается и точка помпажа компрессора прогнозируется как точка, соответствующая на карте компрессора заданной скорости и последнему вычисленному значению массового расхода на входе компрессора и коэффициенту давления компрессора. 7 з.п. ф-лы, 4 ил.
Description
ОБЛАСТЬ, К КОТОРОЙ ОТНОСИТСЯ ИЗОБРЕТЕНИЕ
Настоящее изобретение относится к реализуемому компьютером способу прогнозирования помпажа (пульсации) в рабочей характеристике газового компрессора посредством расчета CFD (ʺкомпьютерного моделирования динамики текучей средыʺ).
ПРЕДШЕСТВУЮЩИЙ УРОВЕНЬ ТЕХНИКИ
В процессах проектирования и разработки газового компрессора, наибольшую важность имеет возможность прогнозировать запас по помпажу с помощью эффективного по затратам расчета устойчивой CFD. Однако в современном уровне техники, прогнозирование запаса по помпажу с достоверностью может быть в значительной степени затруднительным. Это, в частности, обусловлено недостатками современных моделей CFD, т.е. областей, которые сопрягаются и решаются в расчете CFD. Как правило, известно применение методологий расчета CFD к моделям, содержащим:
- геометрию впуска, включающую в себя впускные направляющие лопатки;
- геометрию ступеней ротора;
- геометрию ступеней статора;
- геометрию выпуска, включающую в себя выходные направляющие лопатки.
Граничные условия среднего статического давления, применяемые непосредственно на выпуске выходных направляющих лопаток, изменяются для моделирования различных коэффициентов давления. Когда достигается численная нестабильность, предполагается, что достигнута оценка точки помпажа. Однако оценка реального компрессора за счет численной нестабильности может быть неточной, ввиду явно ограничительного характера граничного условия, применяемого настолько близко к выпуску выходных направляющих лопаток, тем самым приводя к искусственному разделению потока вблизи выхода, вызывающему численную неустойчивость. Фактически, за реальным компрессором обычно следует либо система сгорания, либо выходная нагнетательная камера, которые обычно не учитываются в известных моделях CFD, описанных выше.
Задача настоящего изобретения состоит в том, чтобы преодолеть эти недостатки, чтобы достичь лучшей оценки реальной точки помпажа компрессора посредством расчета CFD.
КРАТКОЕ ОПИСАНИЕ СУЩНОСТИ ИЗОБРЕТЕНИЯ
Чтобы решить основную задачу, определенную выше, предложен реализуемый компьютером способ в соответствии с независимым пунктом формулы изобретения. Зависимые пункты формулы изобретения описывают предпочтительные варианты осуществления и модификации изобретения.
В соответствии с настоящим изобретением, предложен реализуемый компьютером способ прогнозирования точки помпажа компрессора, содержащий следующую последовательность этапов:
- генерирование множества сеток, подходящих для расчета посредством ʺкомпьютерного моделирования динамики текучей средыʺ компрессора, причем сетки включают в себя:
сетку впуска компрессора,
по меньшей мере сетку ступени ротора компрессора,
сетку выпуска компрессора, представляющую множество выходных направляющих лопаток и выходное сопло и продолжающуюся до концевой площади выходного сечения сопла,
- сборку множества сеток для получения вычислительной области ʺкомпьютерного моделирования динамики текучей средыʺ,
- задание граничных условий вычислительной области,
- задание условий атмосферного давления на концевой площади выходного сечения сопла,
- вычисление посредством расчета ʺкомпьютерного моделирования динамики текучей средыʺ потока вдоль компрессора для расчета массового расхода на впуске компрессора и коэффициента давления компрессора,
- проверку, достигнут ли предопределенный предел численной стабильности, и если предел численной стабильности не достигнут, то:
- уменьшение размеров концевой площади выходного сечения сопла,
- генерирование заново сетки выпуска компрессора,
- повторение упомянутых этапов задания атмосферного давления на концевой площади выходного сечения сопла, упомянутого этапа вычисления и упомянутого этапа проверки.
Преимущественно, сопло предоставляется в вычислительной области в направлении вниз по потоку от выходных направляющих лопаток, чтобы постепенно ослабить давление до атмосферных уровней. Граничные условия атмосферного давления применяются на выпуске сопла независимо от моделируемого коэффициента давления. Тогда достигаются различные коэффициенты давления и массовые расходы путем изменения площади выходного сечения сопла до достижения численной нестабильности.
Предполагается, что коэффициент давления и массовый расход на впуске, которые рассчитаны при достижении точки численной нестабильности, соответствуют реальной точке помпажа компрессора. Экспериментальные испытания подтвердили, что точка помпажа, рассчитанная согласно настоящему изобретению, фактически близка к реальному помпажу компрессора с допустимым пределом погрешности.
В соответствии с возможным вариантом осуществления настоящего изобретения, длина сетки выпуска компрессора равна n-кратной длине выходных направляющих лопаток компрессора. Более конкретно, n может находиться в пределах от 10 до 15. Длина сетки выпуска выбирается так, чтобы поддерживать удобным образом и постепенно снижать давление вниз по потоку от выходных направляющих лопаток до атмосферных уровней.
В соответствии с другим возможным вариантом осуществления настоящего изобретения, длина сетки выпуска компрессора поддерживается постоянной на протяжении этапов способа, в частности, в течение этапа уменьшения размеров концевой площади выходного сечения сопла, если предел численной стабильности еще не достигнут. Предпочтительным образом, это позволяет работать только с одним граничным геометрическим параметром при модификации сетки выпуска компрессора между двумя последовательными выполнениями этапа вычисления.
В соответствии с другим возможным вариантом осуществления настоящего изобретения, сетка выпуска компрессора содержит секцию, в которой выходное давление компрессора определяется для расчета коэффициента давления компрессора.
Выходное давление затем используется для заполнения карты компрессора точками, определенными настоящим способом. При использовании точек на карте компрессора, вычисленных способом согласно настоящему изобретению, запас по помпажу компрессора, при эксплуатации компрессора, может прогнозироваться более точно.
В соответствии с другим возможным вариантом осуществления настоящего изобретения, этап задания граничных условий вычислительной области содержат подэтапы:
- задания атмосферных граничных условий на впуске,
- задания скорости вращающихся поверхностей,
- задания плоскости смешивания на последовательных границах раздела.
Предпочтительным образом, это обеспечивает реалистичную модель CFD компрессора.
КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ
Аспекты, определенные выше, и дальнейшие аспекты настоящего изобретения, будут понятны из примеров варианта осуществления, описанного ниже, и поясняются со ссылкой на примеры варианта осуществления. Изобретение будет описано более подробно ниже со ссылкой на примеры варианта осуществления, которые не ограничивают изобретение.
Фиг. 1 показывает продольный разрез газотурбинного двигателя, включающего в себя компрессор, характеристики которого могут прогнозироваться способом согласно настоящему изобретению,
Фиг. 2 показывает модель компрессора согласно фиг. 1,
Фиг. 3 показывает блок-схему последовательности операций способа согласно настоящему изобретению,
Фиг. 4 показывает, на карте компрессора, сравнение между экспериментальными точками и точками, рассчитанными в соответствии со способом согласно настоящему изобретению.
ПОДРОБНОЕ ОПИСАНИЕ
Фиг. 1 показывает пример газотурбинного двигателя 10 в разрезе. Газотурбинный двигатель 10 содержит, в последовательности потока, впускное отверстие (впуск) 12 для воздуха, секцию 14 компрессора, секцию 16 камеры сгорания и секцию 18 турбины, которые обычно расположены в последовательности потока и, как правило, вокруг и в направлении продольной оси или оси 20 вращения. Газотурбинный двигатель 10 дополнительно содержит вал 22, который может вращаться вокруг оси 20 вращения и который проходит в продольном направлении через газотурбинный двигатель 10. Вал 22 соединяет с возможностью привода секцию 18 турбины с секцией 14 компрессора.
При работе газотурбинного двигателя 10, воздух 24, который поступает через впускное отверстие 12 для воздуха, сжимается секцией 14 компрессора и подается в секцию 16 горелки секции сгорания. Секция 16 горелки содержит нагнетательную камеру 26 горелки, одну или более камер 28 сгорания и по меньшей мере одну горелку 30, закрепленную в каждой камере 28 сгорания. Камеры 28 сгорания и горелки 30 расположены внутри нагнетательной камеры 26 горелки. Сжатый воздух, проходящий через компрессор 14, входит в диффузор 32 и выпускается из диффузора 32 в нагнетательную камеру 26 горелки, откуда часть воздуха поступает в горелку 30 и смешивается с газообразным или жидким топливом. Затем воздушно-топливная смесь сжигается, и газообразные продукты 34 сгорания или рабочий газ из сгорания направляется через камеру 28 сгорания в секцию 18 турбины через переходный канал 17.
Секция 18 турбины содержит несколько несущих лопасти дисков 36, прикрепленных к валу 22. В настоящем примере, имеется два диска 36, каждый из которых несет кольцевой набор лопастей 38 турбины. Однако количество дисков, несущих лопасти, может быть различным, то есть, только один диск или более двух дисков. Кроме того, направляющие лопатки 40, которые прикреплены к статору 42 газотурбинного двигателя 10, расположены между ступенями кольцевых наборов лопастей 38 турбины. Между выходом камеры 28 сгорания и передними лопастями 38 турбины предусмотрены направляющие лопатки 44 впуска, которые направляют поток рабочего газа на лопасти 38 турбины.
Газообразные продукты сгорания из камеры 28 сгорания входят в секцию 18 турбины и приводят в движение лопасти 38 турбины, которые, в свою очередь, вращают вал 22. Направляющие лопатки 40, 44 служат для оптимизации угла подачи газообразных продуктов сгорания или рабочего газа на лопасти 38 турбины.
Секция 18 турбины приводит в действие секцию 14 компрессора. Секция 14 компрессора содержит осевую последовательность ступеней 46 лопаток статора и ступеней 48 лопастей ротора. Ступени 48 лопастей ротора содержат диск ротора, поддерживающий кольцевой набор лопастей. Секция 14 компрессора также содержит корпус 50, который окружает ступени ротора и поддерживает ступени 48 лопастей. Ступени направляющих лопаток включают в себя кольцевой набор радиально продолжающихся лопаток, которые смонтированы на корпусе 50. Лопатки предусмотрены для обеспечения потока газа под оптимальным углом для лопастей в данной рабочей точке двигателя. Некоторые из ступеней направляющих лопаток имеют переменные лопатки, где угол лопаток, вокруг их собственной продольной оси, может регулироваться для достижения угла в соответствии с рабочими характеристиками воздушного потока, которые могут возникать при различных условиях работы двигателя.
Первая ступень направляющих лопаток, непосредственно выше по потоку от первой в направлении вверх по потоку ступени лопастей ротора, идентифицируется как направляющие лопатки впуска (ʺIGVʺ) 57. Последняя ступень направляющих лопаток, непосредственно ниже по потоку от последней в направлении вниз по потоку ступени лопастей статора и выше по потоку от диффузора 32, обозначена как направляющие лопатки выпуска (ʺEGVʺ) 49.
Корпус 50 определяет радиально наружную поверхность 52 канала 56 впуска компрессора 14. Радиально внутренняя поверхность 54 канала 56 впуска по меньшей мере частично определяется барабаном 53 ротора, который частично определяется кольцевым набором лопастей 48.
Термины ʺвыше по потокуʺ и ʺниже по потокуʺ относятся к направлению воздушного потока и/или потока рабочего газа через двигатель, если не указано иное. Термины ʺвверх по потокуʺ и ʺвниз по потокуʺ относятся к общему потоку газа через двигатель. Термины ʺосевойʺ, ʺрадиальныйʺ и ʺокружнойʺ приводятся со ссылкой на ось 20 вращения двигателя.
На фиг. 2 показан пример модели 210 секции 14 компрессора в соответствии со способом 100 согласно настоящему изобретению.
В общем, согласно настоящему изобретению, любой компрессор может быть смоделирован, например, независимо от того факта, что компрессор может быть связан или не связан с турбиной. Секция 16 горелки и секция 18 турбины фактически не моделируются в соответствии со способом согласно настоящему изобретению, который выполняется для моделирования только рабочих характеристик компрессора и расчета рабочих точек компрессора, которые обычно определяет карту компрессора.
Могут также моделироваться одноступенчатые компрессоры.
Модели компрессора в соответствии с настоящим изобретением учитывают только окружную часть компрессора вокруг его оси вращения. Окружная протяженность модели считается вращательно-периодической.
Многоступенчатая модель 210 создается, в соответствии со способом согласно настоящему изобретению, посредством последовательности этапов, схематично представленных на блок-схеме на фиг. 3.
На первом этапе 110 способа 100 генерируется множество сеток, подходящих для расчета ʺкомпьютерного моделирования динамики текучей средыʺ (CFD) для компрессора 14. Множество сеток включает в себя:
- сетку 220 впуска компрессора, представляющую ступень направляющих лопаток 57 впуска;
- множество сеток 225 ступеней ротора компрессора представляющих, соответственно, множество ступеней 48 лопастей ротора;
- множество сеток 230 ступеней статора компрессора, представляющих, соответственно, множество ступеней 46 лопастей статора;
- сетку 235 выпуска компрессора, представляющую множество выходных направляющих лопаток 49 и, ниже по потоку от них, выходное сопло, продолжающееся до концевой площади 240 выходного сечения сопла.
Сетка 235 выпуска компрессора должна обрабатываться как одиночная сетка, представляющая как выходные направляющие лопатки 49, так и сопло, а не как множество сеток, разделенных искусственными границами раздела.
Сетка 235 выпуска компрессора имеет длину L, которая выбрана как n-кратное значение длины выходных направляющих лопаток 49 компрессора 14, причем n обычно является целым числом в пределах от 10 до 15. Выбрано первое значение площади A1 концевой площади 240 выходного сечения сопла.
Сопло вниз по потоку от выходных направляющих лопаток обеспечивает постепенное сопряжение между давлением выпуска вниз по потоку от выходных направляющих лопаток 49 и условиями атмосферного давления. На втором этапе 120 способа 100, множество сеток 220, 225, 230, 235 компонуются вместе, чтобы получить вычислительную область 250 ʺкомпьютерного моделирования динамики текучей средыʺ.
На третьем этапе 130 способа 100, задаются граничные условия вычислительной области 250. Третий этап 130 способа 100 содержит подэтапы:
- задание атмосферных граничных условий на впуске вычислительной области 250;
- задание скорости вращающихся поверхностей, то есть поверхностей, представляющих ступени 48 лопастей ротора;
- задание плоскостей смешивания, представляющих искусственные границы раздела между сетками 220, 225, 230, 235 вычислительной области 250, для учета переходов потока между неподвижными и вращающимися частями компрессора;
- определение, какие из поверхностей вращательно-периодической модели 210 должны рассматриваться как вращательно-периодические.
После третьего этапа 130, способ 100 содержит цикл 101 этапов, который включает в себя:
- четвертый этап 140 задания условий атмосферного давления на концевой площади 240 выходного сечения сопла,
- пятый этап 150 вычисления, посредством расчета CFD, потока вдоль компрессора 14 для расчета, в частности, скорректированного массового расхода на впуске компрессора и общего коэффициента давления компрессора между значениями давления в предопределенной секции 260 сетки 235 выпуска компрессора и на впуске вычислительной области 250. Предопределенная секция 260 представляет в сетке 235 выпуска компрессора секцию, расположенную непосредственно ниже по потоку от выходных направляющих лопаток 49 в компрессоре 14.
После пятого этапа 150, цикл 101 дополнительно включает в себя шестой этап 160 проверки того, достигнут ли предопределенный предел численной стабильности, и если предел численной стабильности не достигнут, то цикл 101 продолжается со следующими этапами:
- седьмой этап 170 заполнения карты М компрессора точкой, представляющей скорректированный массовый расход на впуске компрессора и общий коэффициент давления компрессора, вычисленный на предыдущем пятом этапе 150;
- восьмой этап 180 уменьшения размеров A1 концевой площади 240 выходного сечения сопла, при сохранении длины L постоянной. В частности, согласно варианту осуществления настоящего изобретения, A1 может уменьшаться с шагом 5%, то есть, например, если начальное значение A1 составляет 20 мм, второе значение составляет 19 мм, третье - 18 мм и так далее;
- девятый этап 190 повторной генерации сетки 235 выпуска компрессора.
После девятого этапа 190, повторяются четвертый этап 140, пятый этап 150 и шестой этап 160.
При каждом выполнении вычисления пятого этапа 150, сетка 235 выпуска компрессора немного отличается от предыдущего выполнения вычисления пятого этапа 150. Это приводит к различным значениям скорректированного массового расхода на впуске компрессора и общего коэффициента давления компрессора и, следовательно, к соответствующим различным точкам на карте M компрессора.
Цикл 101 заканчивается, когда на шестом этапе 160 достигнута численная стабильность. Точка на карте, рассчитанная во время последнего выполнения вычисления пятого этапа 150, считается репрезентативной для соответствующей точки помпажа компрессора.
Результаты, которые могут быть достигнуты с помощью настоящего изобретения, показаны на фиг. 4, где показана карта М компрессора, включающая в себя горизонтальную ось Х для скорректированного массового расхода на впуске компрессора и вертикальную ось Y для общего коэффициента давления компрессора. Для компрессора, вращающегося с заданной скоростью, треугольные метки 102 представляют экспериментальные тестовые данные. Среди множества меток 102, метка 102a, характеризующаяся более высоким общим коэффициентом давления и более низким скорректированным массовым расходом на впуске компрессора, представляет экспериментальную точку помпажа компрессора при заданной скорости вращения.
Для того же самого компрессора, множество точек, рассчитанных способом 100, представлено на карте M прямоугольными метками 103. Кривая 104, соединяющая прямоугольные метки 103, представляет на карте M компрессора кривую скорости при заданной скорости вращения. Другая аналогичная кривая может быть рассчитана с учетом другого значения скорости вращения во время выполнения способа 100, в частности, третьего этапа 130. Среди множества меток 103, метка 103a, характеризуемая более высоким общим коэффициентом давления и более низким скорректированным массовым расходом на впуске компрессора, представляет собой расчетную точку помпажа компрессора при заданной скорости вращения.
Расстояние между экспериментальной и вычисленной точками 102а, 103а помпажа как по значению коэффициента компрессора, так и по значению массового расхода на впуске, обеспечивает удовлетворительное подтверждение данного способа, будучи значительно меньшим, чем аналогичные расстояния, которые могут быть получены другими существующими способами CFD.
Claims (23)
1. Реализуемый компьютером способ (100) прогнозирования точки помпажа компрессора (14), содержащий следующую последовательность этапов:
- генерирование (110) множества сеток, подходящих для расчета ʺкомпьютерного моделирования динамики текучей средыʺ компрессора (14), причем сетки включают в себя:
сетку (220) впуска компрессора, представляющую ступень впускных направляющих лопаток (57) компрессора (14),
по меньшей мере сетку (225) ступени ротора компрессора, представляющую ступени (48) лопастей ротора компрессора (14),
сетку (235) выпуска компрессора, представляющую множество выходных направляющих лопаток (49) и выходное сопло и продолжающуюся до концевой площади (240) выходного сечения сопла сетки (235) выпуска,
- сборку (120) множества сеток (220, 225, 235) для получения вычислительной области (250) ʺкомпьютерного моделирования динамики текучей средыʺ,
- задание (130) граничных условий вычислительной области (250), включая задание по меньшей мере скорости поверхностей вычислительной области (250), представляющих ступени (48) лопастей ротора компрессора (14),
- задание (140) условий атмосферного давления на концевой площади (240) выходного сечения сопла,
- вычисление (150), посредством расчета ʺкомпьютерного моделирования динамики текучей средыʺ, потока вдоль компрессора (14) и вычисление массового расхода на входе в компрессор и коэффициента давления компрессора,
- проверку (160), достигнут ли предопределенный предел численной стабильности, и, если предел численной стабильности не достигнут, то:
- уменьшение (180) размеров (A1) концевой площади (240) выходного сечения сопла,
- повторное генерирование (190) сетки (235) выпуска компрессора,
- повторение упомянутых этапов задания атмосферного давления (140) на концевой площади (240) выходного сечения сопла, упомянутого этапа вычисления (150) и упомянутого этапа проверки (160),
- если предел численной стабильности достигнут, то способ останавливается, и точка помпажа компрессора (14) прогнозируется как точка, соответствующая на карте (M) компрессора (14) заданной скорости и последнему вычисленному значению массового расхода на входе компрессора и коэффициенту давления компрессора.
2. Реализуемый компьютером способ (100) по п. 1, в котором длина (L) сетки (235) выпуска компрессора равна n-кратной длине выходных направляющих лопаток (49) компрессора (14).
3. Реализуемый компьютером способ (100) по п. 2, в котором n представляет собой целое число от 10 до 15.
4. Реализуемый компьютером способ (100) по п. 2 или 3, в котором длина (L) сетки (235) выпуска компрессора поддерживается постоянной на упомянутом этапе уменьшения (180) размеров (A1) концевой площади (240) выходного сечения сопла.
5. Реализуемый компьютером способ (100) по любому из предыдущих пунктов, в котором сетка (235) выпуска компрессора содержит секцию (260), в которой определяется давление на выходе компрессора для вычисления коэффициента давления компрессора.
6. Реализуемый компьютером способ (100) по любому из предыдущих пунктов, в котором этап задания граничных условий вычислительной области (250) содержат подэтапы:
- задания атмосферных граничных условий на впуске,
- задания плоскости смешивания на границах раздела между сетками (220, 225, 230, 235) вычислительной области (250).
7. Реализуемый компьютером способ (100) по п. 1, причем способ (100) дополнительно включает в себя этап заполнения (170) карты (М) компрессора точкой, представляющий массовый расход на впуске компрессора и общий коэффициент давления компрессора, вычисленные на предыдущем этапе вычисления (150) потока вдоль компрессора (14).
8. Реализуемый компьютером способ (100) по п. 7, в котором, если предел численной стабильности достигнут, то точка на карте, вычисленная в течение последнего выполнения этапа вычисления (150) потока вдоль компрессора (14), считается репрезентативной для соответствующей точки помпажа компрессора.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP15181154.4 | 2015-08-14 | ||
| EP15181154.4A EP3131025A1 (en) | 2015-08-14 | 2015-08-14 | Method for the prediction of surge in a gas compressor |
| PCT/EP2016/068760 WO2017029131A1 (en) | 2015-08-14 | 2016-08-05 | Method for the prediction of surge in a gas compressor |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| RU2678155C1 true RU2678155C1 (ru) | 2019-01-23 |
Family
ID=54010864
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2018105189A RU2678155C1 (ru) | 2015-08-14 | 2016-08-05 | Способ прогнозирования помпажа в газовом компрессоре |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US11187235B2 (ru) |
| EP (2) | EP3131025A1 (ru) |
| JP (1) | JP6598983B2 (ru) |
| CN (1) | CN107924425B (ru) |
| RU (1) | RU2678155C1 (ru) |
| WO (1) | WO2017029131A1 (ru) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20180163736A1 (en) * | 2016-12-09 | 2018-06-14 | General Electric Company | Systems and methods for operating a compression system |
| US10344767B2 (en) * | 2017-05-01 | 2019-07-09 | GM Global Technology Operations LLC | Method for compressor surge detection to enable model base air estimation |
| CN109684750A (zh) * | 2018-12-27 | 2019-04-26 | 中国船舶重工集团公司第七0三研究所 | 一种用于氦气轮机的压气机模型建模方法 |
| CN109815624B (zh) * | 2019-02-21 | 2022-01-04 | 南京航空航天大学 | 一种考虑进气总压畸变影响的压气机稳定边界判断方法 |
| CN110005643B (zh) * | 2019-03-07 | 2020-10-02 | 北航（四川）西部国际创新港科技有限公司 | 基于面积律的跨声速轴流压气机机匣设计方法 |
| CN110728052B (zh) * | 2019-10-11 | 2023-09-22 | 中国航发沈阳发动机研究所 | 一种确定旋转盘腔相似试验边界条件的方法 |
| CN111079239B (zh) * | 2019-12-19 | 2023-07-21 | 中国航空发动机研究院 | 一种仿生压气机叶栅造型方法 |
| CN111737307B (zh) * | 2020-04-10 | 2024-04-09 | 中国电力工程顾问集团中南电力设计院有限公司 | 基于配电网典设边界条件全要素的自动校核方法 |
| CN114036652B (zh) * | 2021-08-06 | 2024-02-23 | 西安交通大学 | 一种燃气轮机叶片内外耦合的多尺度计算方法 |
| CN113836658B (zh) * | 2021-09-17 | 2024-02-27 | 中信建筑设计研究总院有限公司 | 一种y型铸钢节点抗压设计承载力的计算方法 |
| CN114810646B (zh) * | 2022-03-31 | 2023-03-24 | 清华大学 | 一种基于平行压气机改进模型的喘振边界判定方法 |
| CN116484771B (zh) * | 2023-06-21 | 2023-08-25 | 陕西空天信息技术有限公司 | 一种轴流压气机cfd网格生成的方法及装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5217348A (en) * | 1992-09-24 | 1993-06-08 | United Technologies Corporation | Turbine vane assembly with integrally cast cooling fluid nozzle |
| US20060034689A1 (en) * | 2004-08-11 | 2006-02-16 | Taylor Mark D | Turbine |
| RU2326271C2 (ru) * | 2002-10-04 | 2008-06-10 | Дженерал Электрик Компани | Способ (варианты) и система для обнаружения признаков замедления и помпажа компрессора |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6905303B2 (en) * | 2003-06-30 | 2005-06-14 | General Electric Company | Methods and apparatus for assembling gas turbine engines |
| GB0426439D0 (en) | 2004-12-02 | 2005-01-05 | Rolls Royce Plc | Rotating stall |
| JP4375369B2 (ja) | 2006-08-10 | 2009-12-02 | トヨタ自動車株式会社 | 過給機付き内燃機関の制御装置 |
| US8096756B2 (en) * | 2008-03-07 | 2012-01-17 | Pratt & Whitney Canada Corp. | Apparatus and method for controlling a compressor |
| CN101599104B (zh) * | 2009-07-16 | 2011-06-22 | 北京航空航天大学 | 一种航空涡轮发动机叶片颤振边界的模拟方法 |
| CN101718269A (zh) | 2009-11-20 | 2010-06-02 | 西安交通大学 | 基于频带能量的喘振先兆识别方法 |
| US8452515B2 (en) | 2011-09-15 | 2013-05-28 | General Electric Company | System and method for simulating a gas turbine compressor |
| JP6188069B2 (ja) * | 2013-10-17 | 2017-08-30 | 三菱重工業株式会社 | 圧縮機、及びガスタービン |
| GB201410180D0 (en) * | 2014-06-09 | 2014-07-23 | Rolls Royce Plc | Method and apparatus for controlling a compressor of a gas turbine engine |
| JP6866019B2 (ja) * | 2014-06-24 | 2021-04-28 | コンセプツ エヌアールイーシー，エルエルシー | ターボ機械の流動制御構造及びその設計方法 |
| US20160123175A1 (en) * | 2014-11-05 | 2016-05-05 | General Electric Company | Hybrid model based detection of compressor stall |
| US9932985B2 (en) * | 2015-02-03 | 2018-04-03 | Honeywell International Inc. | Gas turbine engine compressors having optimized stall enhancement feature configurations and methods for the production thereof |
-
2015
- 2015-08-14 EP EP15181154.4A patent/EP3131025A1/en not_active Withdrawn
-
2016
- 2016-08-05 EP EP16756957.3A patent/EP3335136B1/en active Active
- 2016-08-05 CN CN201680048441.6A patent/CN107924425B/zh active Active
- 2016-08-05 JP JP2018507557A patent/JP6598983B2/ja not_active Expired - Fee Related
- 2016-08-05 US US15/747,498 patent/US11187235B2/en active Active
- 2016-08-05 WO PCT/EP2016/068760 patent/WO2017029131A1/en active Application Filing
- 2016-08-05 RU RU2018105189A patent/RU2678155C1/ru active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5217348A (en) * | 1992-09-24 | 1993-06-08 | United Technologies Corporation | Turbine vane assembly with integrally cast cooling fluid nozzle |
| RU2326271C2 (ru) * | 2002-10-04 | 2008-06-10 | Дженерал Электрик Компани | Способ (варианты) и система для обнаружения признаков замедления и помпажа компрессора |
| US20060034689A1 (en) * | 2004-08-11 | 2006-02-16 | Taylor Mark D | Turbine |
Non-Patent Citations (3)
| Title |
|---|
| Stephanie Bergqvist. "Prediction of Turbo Compressor Maps using CFD Master's thesis in Applied Mechanics", 01.07.2014. * |
| Stephanie Bergqvist. "Prediction of Turbo Compressor Maps using CFD Master's thesis in Applied Mechanics", 01.07.2014. Yu Wang et al. "Simulation and performance analysis on centrifugal compressors of different dimensions and variable operation speed", 16.04.2011. * |
| Yu Wang et al. "Simulation and performance analysis on centrifugal compressors of different dimensions and variable operation speed", 16.04.2011. * |
Also Published As
| Publication number | Publication date |
|---|---|
| US20180216623A1 (en) | 2018-08-02 |
| EP3335136B1 (en) | 2019-10-16 |
| EP3335136A1 (en) | 2018-06-20 |
| CN107924425B (zh) | 2021-03-09 |
| US11187235B2 (en) | 2021-11-30 |
| EP3131025A1 (en) | 2017-02-15 |
| JP6598983B2 (ja) | 2019-10-30 |
| WO2017029131A1 (en) | 2017-02-23 |
| JP2018532178A (ja) | 2018-11-01 |
| CN107924425A (zh) | 2018-04-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2678155C1 (ru) | Способ прогнозирования помпажа в газовом компрессоре | |
| US20210209264A1 (en) | Modeling and calculation aerodynamic performances of multi-stage transonic axial compressors | |
| Pachidis et al. | Towards a full two dimensional gas turbine performance simulator | |
| Barbosa et al. | Influence of Variable Geometry Transients on the Gas Turbine Performance | |
| Léonard et al. | A quasi-one-dimensional CFD model for multistage turbomachines | |
| Chen et al. | Effects of nonuniform tip clearance on fan performance and flow field | |
| Dickmann et al. | Unsteady flow in a turbocharger centrifugal compressor: 3D-CFD-simulation and numerical and experimental analysis of impeller blade vibration | |
| Baturin et al. | Identifying the approach to significantly improve the performance of NK-36ST gas turbine power plant | |
| US11568099B2 (en) | System and process for designing internal components for a gas turbine engine | |
| Hegde et al. | Influence of Disc Modes and Sideband Excitations on the Mistuned Forced Response Behaviour of an Embedded Compressor Rotor | |
| Touyeras et al. | Aerodynamic design and test result analysis of a three stage research compressor | |
| US10552555B2 (en) | Method for the prediction of turbomachine performances | |
| Riegler et al. | Some aspects of modelling compressor behavior in gas turbine performance calculations | |
| Gao et al. | Flow interactions between shrouded power turbine and nonaxisymmetric exhaust volute for marine gas turbines | |
| Lotz | Aerodynamic Optimization Process for Turbocharger Compressor Impellers | |
| Hariharan et al. | Effect of inlet clearance gap on the performance of an industrial centrifugal blower with parallel wall volute | |
| Grigoriev et al. | Benchmarking CD-Adapco’s Star-CCM+ in a production design environment | |
| Patel et al. | CFD analysis of mixed flow pump | |
| Baghdadi et al. | Compressor stability analysis | |
| Gezork et al. | Influence of tip shroud cavity detailing on turbine blade forcing calculations | |
| Yamagami et al. | CFD Modeling effects on unsteady multistage simulation for a transonic axial compressor | |
| Li et al. | The impact of the multiple reference frame interface on modelling the interaction between IGVs and the impeller in turbocharger compressors | |
| Petrovic et al. | Part Load Behavior of the LP Part of an Industrial Gas Turbine | |
| Nigam et al. | Design of Parametric Modeller for Centrifugal Compressors | |
| Buske | A CFD-Based Prediction Method for Tip Clearance Losses and Deviations in Axial Turbines |