Перейти к публикации

Не повышается температура после рабочего колеса центробежного компрессора


seyka

Рекомендованные сообщения

Добрый день!

Столкнулся со следующей проблемой: рассчитываю ступень центробежного компрессора. Моделирую её работу на расчетном режиме работы. Рабочая среда - природный газ. Его свойства задаю таблицами. 

Начальное давление 6,8 МПа, начальная температура - 15 С - параметры заданные на входе в проточную часть. На выходе задавалась массовая скорость. Вращение рабочего колеса задавалось по средствам подвижного тела. 

Суть проблемы. Произведен расчет, при этом контролируемые параметры (расход на входе, температура и давление на выходе из проточной части) вышли на горизонтали. Проблема заключается в том, что температура повышается незначительно (на десятые градуса). При этом повышение давления происходит ожидаемо и степень повышения давления равна 1,21.

С чем может быть связано поведение температуры после элемента в котором подводиться энергия?

Ссылка на сообщение
Поделиться на других сайтах


какую вы ожидаете температуру?

энергия ушла в давление, объем постоянным считать нельзя.

 

хотя от этих задач я далек...

 

посмотрите еще такие параметры как опорное давление и опорная температура.

М.б. у вас есть неувязка с этим, и с вашими граничными условиями?

Ссылка на сообщение
Поделиться на других сайтах

энергия подводимая в рабочем колесе срабатывает на повышение давления и температуры. Помимо этого в проточной части происходит перераспределение потенциальной энергии в термическую + температура газа повышается за счет вязкостного трения газа. Таким образом температура может повышаться до 15-20 гр С. Эта величина зависит выходного угла наклона лопатки рабочего колеса, условного числа Маха для рабочего колеса и сжимаемой среды. В моей расчетной и экспериментальной практике минимальное значение которое я получал было 5 гр С.

 

М.б. у вас есть неувязка с этим, и с вашими граничными условиями?

чет я на понял о чем Вы говорите( 

Ссылка на сообщение
Поделиться на других сайтах

Добрый вечер! Учтен ли теплоперенос в расчете? По моему во вкладке где выбирается конвекция и теплоперенос.

Ссылка на сообщение
Поделиться на других сайтах

Добрый вечер! Учтен ли теплоперенос в расчете? По моему во вкладке где выбирается конвекция и теплоперенос.

без включенного теплопереноса не задать температуру в этом пакете.

 

В 3,09,03 есть такая штука, как опорные значения расчета - самый верх дерева пре-пост-процессора.

Далее все величины - начальные и на ГУ - задаются относительно них. 

Поэтому если вы задали начальное давление 6,8 МПа, начальная температура - 15 град, то в реальности у вас 7,8Мпа и 315 град.

Смотрите раздел "общие установки"

 

Вы можете выложить ссылку на архив с клиентской частью проекта - м.б. я смогу его открыть и посмотреть.

Ссылка на сообщение
Поделиться на других сайтах

проект смогу выложить в понедельник вечером, поскольку он у меня на работе находиться. 

По поводу опорных величин. Да действительно давление порядка 6,8 МПа (точную величину не помню на память), а вот температура у меня традиционно ставится 273,15. Следовательно выдаваемые величины в результате расчета являются градусами цельсия:)

По поводу включения теплопереноса - да включен.

В добавок по версиям, есть и 3,09,4 - но я ее мало тестил с подвижным телом, поэтому для данной задачи решил использовать проверенную 3,09,03

 

Спасибо за не равнодушие к моей задаче:)

Ссылка на сообщение
Поделиться на других сайтах

обновляйтесь до 3,09,05 - там задачи типа ротор-статор можно решать в секторной постановке с использованием периодического ГУ.

Ссылка на сообщение
Поделиться на других сайтах

именно для этой задачи необходимо решать через подвижное тело - необходимо учитывать теретечки через пазухи у рабочего колеса. Хм 05 версия мне уже будет недоступна техническая поддержка закончилась и дальше не проплатили её(((

Ссылка на сообщение
Поделиться на других сайтах
Рабочая среда - природный газ.
Начальное давление 6,8 МПа, начальная температура - 15 С -
При этом повышение давления происходит ожидаемо и степень повышения давления равна 1,21.
какую вы ожидаете температуру?
Это Газпромовская компрессорная станция. Должно скакнуть не менее чем на 30-40 градусов.
Ссылка на сообщение
Поделиться на других сайтах

Задайте условия исключительно адиабатического процесса сжатия газа в турбине.

Для этих  условий программа обязана выполнить вычисления в соответствии 

с объединенным газовым законом (p*V)/T=const.

Должны измениться все три параметра и температура в том числе.

Для этих условий изменение температуры будет максимальным.

Затем добавьте условия теплообмена, изменение температуры д.б. 

меньше, чем для адиабатического процесса.

Хотя если считать, что в турбине и корпусе установилось тепловое равновесие,

то это будет означать квазиадиабатические условия сжатия.

Ссылка на сообщение
Поделиться на других сайтах

А вообще, если известно давление, температура и расход "до", и известно давление "после", массовый расход не меняется - то можно найти плотность "после", а по плотности и давлению - найти температуру "после" ;) 

Ссылка на сообщение
Поделиться на других сайтах
Это Газпромовская компрессорная станция. Должно скакнуть не менее чем на 30-40 градусов.

на сколько я помню такие ступени в газпроме есть - гл в каждом компрессора по 2 - 3 ступени для линейных станций. 30-40 это для всего компрессора, а для одной ступени где то 12-15. 


А вообще, если известно давление, температура и расход "до", и известно давление "после", массовый расход не меняется - то можно найти плотность "после", а по плотности и давлению - найти температуру "после" ;) 
 

да я могу взять результаты эксперимента и считать не стоит. Основная цель расчета - посмотреть величину перетечки через пазуху между рабочим колесом и статорными элементами проточной части. Эспериментально довольно противно это делается...

Ссылка на сообщение
Поделиться на других сайтах
проект смогу выложить в понедельник вечером, поскольку он у меня на работе находиться. 

как дела? 

Ссылка на сообщение
Поделиться на других сайтах

https://cloud.mail.ru/public/MgSH/viiG8ByrR

клиент расчета

1) зачем заменили плотность с закона идеального газа на таблицу?

2) отрицательную массовую скорость задавать нельзя - выполните урок из туториала по задаче ротор-статор

2.1) в целом надо ставить другое гу на выходе..., не вход/выход....

3) задавать нулевые параметры турбулентности очень не рекомендуется. нужно хотя бы 0,001 везде поставить

 

это первое, что попалось - завтра еще гляну.

Ссылка на сообщение
Поделиться на других сайтах
1) зачем заменили плотность с закона идеального газа на таблицу?

при больших значениях давлений уравнение состояния идеального газа все больше не подходит для расчетов (растет погрешность вычисления его компонентов) по этому были насчитаны массивы плотности по уравнению состояния, разработанного для смесей углеводородов (уравнение ВБР).

2) отрицательную массовую скорость задавать нельзя - выполните урок из туториала по задаче ротор-статор

хм - даж не знаю что ответить. Данный подход решения задачи, а именно задания ГУ разрабатывался совместно с технической поддержкой FV. Более кого, данное задание ГУ много кратно было подтверждено верификацией с экспериментальными данными. Все было шикарно просто!

3) задавать нулевые параметры турбулентности очень не рекомендуется. нужно хотя бы 0,001 везде поставить

Тут согласен!

 

Спасибо за уделяемое время данной проблемы! 

Ссылка на сообщение
Поделиться на других сайтах

при больших значениях давлений уравнение состояния идеального газа все больше не подходит для расчетов (растет погрешность вычисления его компонентов) по этому были насчитаны массивы плотности по уравнению состояния, разработанного для смесей углеводородов (уравнение ВБР).

тогда точность таблицы должна быть выше - это влияет на расчет тепловых задач 

 

хм - даж не знаю что ответить. Данный подход решения задачи, а именно задания ГУ разрабатывался совместно с технической поддержкой FV. Более кого, данное задание ГУ много кратно было подтверждено верификацией с экспериментальными данными. Все было шикарно просто! 

да, я ошибся, нельзя задавать отрицательную скорость, а массовую скорость отрицательную можно.  Нюанс в том, что и скорость, и массовая скорость работают внутри примерно одинаково - оба ГУ работают, по факту, как массовая скорость, поэтому очень интересно, почему в одном месте (-) задать можно, а в другом - нет.

 

еще посмотрю задачку.

еще посмотрите - сейчас есть два теплопереноса - через h и через H.

один для М<0.5 а второй для >0.5.

 

попробуйте переключить.

Ссылка на сообщение
Поделиться на других сайтах

еще нюансы:

 

опорным давлением лучше задать статику на входе (или выходе), чтоб относительное давление на входе (или выходе) было = 0. Это повышает точность/скорость немножко.

Сейчас у вас опора 100000, и полное на входе +7241000, в результате статика на входе +7100000.

А надо : опора 7200000.

после такого изменения задачу надо полностью пересчитывать!

 

сетка грубая. но можно еще в 2 раза загрубить (будет 50тыс вместо 300тыс ячеек), получить сошедшееся решение, и потом проадаптировать 1м уровнем по объему и 2м уровнем по лопаткам.

Для этого лопатки ротора и статора надо было еще в CADе покрасить в другой цвет...В уже сделанном проекте это не так удобно делать. Хотя можно сделать замену геометрии.

 

Далее - шаг по времени некорректно задавать фиксированным. У вас подвижное тело, такие задачи считаются с поверхностным CFL=1 или меньше. (конвективное ставим 100).

 

у меня сейчас на выходе +23 градуса, но это не сошедшееся решение. Но жду, что сойдется примерно к этой цифре, а потом можно проадаптировать сетку и уточнить решение.

 

Вот слева динамика температуры в моем варианте, справа - в вашем

post-26004-0-63329700-1475709023_thumb.jpg

Ссылка на сообщение
Поделиться на других сайтах

Присоединяйтесь к обсуждению

Вы можете опубликовать сообщение сейчас, а зарегистрироваться позже. Если у вас есть аккаунт, войдите в него для написания от своего имени.
Примечание: вашему сообщению потребуется утверждение модератора, прежде чем оно станет доступным.

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставить изображения напрямую. Загрузите или вставьте изображения по ссылке.

  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу.




  • Сообщения

    • Борман
      Если еще не поздно... я и сам погрузился бы. Физика такая, для вала с одним колесом... F - это понятно, центробежное смещение центра масс колеса M - момент от центробежных сил, который не появится, если колесо не повернется при изгибе вала. Обрати внимание, что BW-момент выворачивает ось в сторону действия силы F, т.е. при обратной прецессии они действуют сообща, и в этом опасность.   Отсюда будут разные собственные частоты при прямой и обратной прецессиях.  Теперь осталось понять, что такое прямая и обратная прецессии :), и почему момент при них направлен в разные стороны. По первой части вопроса см. мою картинку в этой теме от 27 августа 2013. По второй части вопроса - нужно включать воображение и усердно рисовать на бумажке. Мне пока сложно.    Физика именно такая. Нужно как то оценить величину этого момента и приложить его в софте, который не умеет ДК. А с тестовыми значениями можно и так порешать, и на выходе должна быть вилка. Т.е. F - задает типа "биссектрису" вилки (престресс-модальник), а момент - угол вилки.  
    • boomeeeer
      @anykeyto Я не оператор. Настраивал работу этой функции на токарно-фрезерном, где нужно было позиционировать шпиндель, для нарезания пазов. Что такое G07 не знаю. Вот такой пример для фрезерного станка нашёл   G108 B0; переключение на ось B для управления положением; G68.2 X10 Y10 Z20 I0 J90 K0 Положение заготовки (10,10,20) Угол Эйлера для определения характеристической системы координат; G53.2 Управление направлением инструмента   G109 B0 - отключит режим
    • anykeyto
      Ага, спасибо... А в связке с G07.1 это как работает? Есть детали ,кольца, на которых надо делать по 4 паза. Сейчас программист написал прогу под шаровую фрезу. Там рабочие движеня происходят через одновременное смещение по оси Х и Z... фреза убивается буквально на третьем кольце, да и обработка получаестя фиговой.. На сколько я понимаю надо стремиться к тому, чтобы на 5ти осном станке фреза всегда была перпендикулярна касательной ли самой поверхности. Мне посоветовал как раз цилндрческую интерполяцию, когда ось С преобразуется в плоскость. В мануале к HNC описано эт криво, но в примерах программ как раз сначала включается CTOS и потом G07.1. но дальше как запрограммировать понимания пока нет...  P.S. хочу добавить, что на 5-снике раньше работал только как оператор по готовым программам  на Хенденхайне. А сейчас на новом предприятии на новой для меня стойке нужно все делать с нуля и самостоятельно. С одной стороны интересно, с другой стороны спросить не у кого, в инете инфы по стойке нет от слова совсем, уж тем более видео....  чем это отличается от просто G0 C180? Или это уже будет шпиндель позиционироваться, а не стол?
    • AlexKaz
      Различают жёсткие и гибкие ротора. В данном случае вал жёсткий. Для него соответствуюшая формула. 4000 об/мин = 66.67 Гц, это ни о чём, эффекты для такой частоты почти не скажутся на жёстком роторе. Но проблема в том, что в лопаточных машинах (а насос на скрине наверное лопаточная?) собственные кратны числу лопаток. Я не разбирался в вопросе чтобы что-то здесь советовать. Лучше не спешить и внимательно почитать про модальный в лопаточных. Заодно проверить, а какое нагружение идёт в деталях насоса чтобы корректно выполнить престресс.
    • ДОБРЯК
      Хотите еще и матрицу масс изменить. :=) И сколько килограмм воды будете присоединять?  Опыт правильного моделирования, опыт расчета таких конструкций накапливается десятилетиями. Маловероятно, что кто-то на форуме будет вам объяснять как правильно моделировать, какие КЭ использовать, какие граничные условия задавать и т. д. Вы же всегда уклоняетесь от разговора на эти темы, даже для простых моделей.  Даже если перейдете на Ансис и не разберётесь в том, что он считает это не поможет вам написать грамотный отчет и объяснить грамотному заказчику, что вы там насчитали. Считает расчетчик, а программа это инструмент. Это только в рекламе пишут, что купите программу трехмерного моделирования, сделайте трехмерную модель, потом нажмите одну кнопку и на обычном игровом компьютере все посчитаете. Только в жизни все немножко иначе...:=)
    • Jesse
      Вот что у Биргера: Центробежные силы и индуцированный гироскопический момент зависят от квадрата частоты вращения вала   Хотя уже сама критическая частота зависит только от момента инерции колеса (альфа, бэта, гамма - это прогибы, углы поворота и прочая шелуха), т.е. частота вращения в формулу не входит, а всё определяется только моментом инерции Jm Ну у него рассматривается простая одномассовая система (фактически грузик на пружинке). Думаю поискать публикации.. мб где то раскопали аналитическую зависимость первой и высшей СЧ от частоты вращения для сложных систем. Хотя если выяснится что зависит,  всё равно обоснование придумывать придётся.  Эхх...
    • Fedor
      https://www.litres.ru/book/v-n-faddeeva/vychislitelnye-metody-lineynoy-algebry-65999878/  Ну и классику жанра :) 
    • Jesse
      4000 об/мин при наружном диаметре колеса 380 мм. Думаете, на больших оборотах вилка всё равно может разойтись и могут появиться гироскопические эффекты? Даже несмотря на вышеприведенные массово-инерционные оценки и расположение опор?   да я влияние престресс (вращение) особо не увидел... Думал обычным модальником прогнать, Ланцошом (он пошустрее).   Ещё мои раздумья направлены в сторону присоединённой массы воды, что по идее должно занижать частоты. Хотя для такого массивного вала как у меня даже на это дело можно забить, наверно...
    • Handrusik
      Maxx2000 - благодарю Павлуха - это не фриланс Утомлённый солнцем - приятно, что меня помнят
    • Борман
      Я так понимаю, что обороты не особо большие. ДК - это чисто высокооборотистая тема.   Просто мысли вслух... Если вилка будет сильно расходится, то всегда будет точка пересечения любого луча с кривой BW. Причем точка будет на меньших оборотах по сравнению с обычным престресс-модальником. Какие дальше действия ? С чем тут нужно дополнительно бороться по сравнению с тем, с чем нужно было бы бороться по результатам обычного престресс-модальника Если вилка не сильно расходится (от кривой обычного престресс-мобальника), то смысла стоить строгую ДК большого нет.    
×
×
  • Создать...