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

Экспертное мнение по расчету модели горения non-premixed


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

Добрый день, господа! В недавней моей теме я спрашивал про программное обеспечение для расчета горения и по советам в итоге остановился на Fluent 17.1.

Поставленная задача довольно специфична, но не такая уж сложна на первый взгляд. Потратив практически две недели на просмотр уроков и чтение мануалов и туториалов приступил к выполнению задачи. Суть задачи состоит в расчете горения газа в вертикальном прямоугольном канале. Особенности модели таковы, что воздух для горения подается самотягой и порционно по высоте канала. Цель такого расположения - контроль распределения температуры (в реальной модели температуры стенки) по высоте канала. Некоторые результаты я получил, но так как сроки выдачи хоть какого-то мало-мальски приемлемого результата поджимают, я решил выложить на всеобщее обозрение и порицание саму модель и какие исходные данные я задавал. Возможно я где-то серьезно ошибаюсь и задаю в корне не верные параметры (т.к. по графику невязок решение сходиться не спешит, но разброс колебаний невязок от итерации к итерации не изменяется).

Вообщем сама модель предоставлена на рисунке ниже.

post-52196-0-98401700-1472206998_thumb.jpg

Вообщем так, на вкладке General выбрана модель Pressure Based, Вкладка Models включено Energy, модель турбулентности k-e с настройками по умолчанию, Radiation Model выбрана модель Diskrete Ordinates (в этом может быть причина несхождения), Spieces Model включена модель Non-Premixed combustion, включена опция Inlet Diffusion, значение  Fuel Stream Rich Flammability Limit установлено 0,8 (изначально здесь я совершил ошибку и поставил значение из туториала и при генерации PDF таблицы значение максимальной температуры равнялось начальной температуре газа и процесса горения как такового не происходило, после расчета стехиометрии и установки коэффициента 0,8 максимальная температура поднялась до 1800 К). На вкладке Boundaries задал мой реальный газ (CH4 - 0.01661, H2 - 0.0015, CO2 - 0.73729, CO - 0.2446 в массовых долях) и сгенерировал PDF таблицу. Температура газа 600 К, окислителя 300 К. На вкладке Materials в pfd-mixture изменил только значение Absorption coefficient на  wsggm-domain-based. Далее перешел к заданию граничных условий. Граничное условие gas-inlet, тип mass-flow-inlet, задал массовый расход на входе 0,033 кг/с, начальное давление 50 Па, параметры турбулентности интенсивность и гидравлический диаметр (применял на всех граничных условиях), параметр Main Mixture Fraction - 1. Граничные условия Air Inlet, тип pressure-inlet, давления по нулям. Граничные условия Outlet, тип pressure-outlet, давление минус 40 Па. Далее вкладка Solution Methods - Sheme - Coupled, Pressure - Presto. Далее Solution Controls, Flow Courant Number - 70, Under-Relaxation Factors, Density - 0.2, Bode Forces - 0.8 (честно признаюсь взял эти значения из туториала). Далее инициализация, число итераций я ставил 2000 и вперед.

Собственно результаты расчета в прикрепленных изображениях.

post-52196-0-45657900-1472207431_thumb.jpg

post-52196-0-58575000-1472207432_thumb.jpg

И мой вопрос на осуждение большинства - все ли я правильно сделал и если это не так, то как исправить либо возможно будут какие-либо рекомендации. Беспокоит еще тот факт, что по мере вычисления постоянно выскакивает предупреждение о возникновении обратного течения на выходе и то что сам флюент не пишет, что расчет сошелся а просто отсчитывает заданное количество итераций.

post-52196-0-93061100-1472207559.jpg

Буду очень благодарен за ответы и советы!

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


Diskrete Ordinates

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

включена опция Inlet Diffusion

не нужно.

значение  Fuel Stream Rich Flammability Limit установлено 0,8

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

начальное давление 50 Па,

не нужно.

тип pressure-inlet

должен быть outlet.

 

+ Не забыть включить гравитацию.

+ специфические настройки схемы

+ подгонка коэффициентов релаксации 

 

(с) DirHald

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

Про гравитацию забыл написать, что она включена. 

должен быть outlet.
Извиняюсь в каком именно граничном условии?

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

вызывает большие сомнения. перепроверьте свои расчеты. должно быть ближе к стехиометрии.
Я перепроверю конечно, но вроде у меня получалось 0,72 плюс рекомендуемые в туториале 10% 
Ссылка на сообщение
Поделиться на других сайтах

Могу конечно попробовать и в СFX но в своей предыдущей теме не дождался ответа, чем же они отличаются и что проще (лучше) из них двоих. Ну я и выбрал Fluent

Ссылка на сообщение
Поделиться на других сайтах
Legat сказал(а) 26 Ауг 2016 - 13:53: начальное давление 50 Па, не нужно.

Вот тут немного не понял. При заданной геометрии и массовом расходе давление определяет скорость истечения. Или я в чем-то ошибаюсь (применительно к расчету во Fluent)? 

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • Zergus
      Выкладывайте сюда свои "чудеса", может по модели получится понять, что в ней не так. А если есть возможность писать в техподдержку - кидайте им свой файл, думаю разберутся.
    • Богоманшин Игорь
      А есть подтверждение в каких либо нормативных документах? Или только из соображений "максимально минимальная себестоимость"?
    • ANT0N1DZE
      Это не поможет ответить на этот вопрос конкретными числами. Судя по всему, там металлические частички внутри материала. Подобные эластомеры существуют и применяются. 
    • gudstartup
      как перевести шпиндель снова первоначальный режим управления скоростью когда ему пофиг число импульсов на оборот? резьбу пока резать не будут а в подпрограмму смены можно что нить добавить или в реализацию m19  чтобы чпу из режима контроля позиционного кодера вышло?
    • gudstartup
    • mactepnew
      Вот я тоже так думал, пока не столкнулся с обратным) Как еще столкнусь, пришлю скриншот. Хаос еще в том, что порядок в дереве построения весьма условный. операция, которая стоит перед групой, имеет бОльший timestamp(пардон, не знаю как это в переводе), чем эта группа, и выполняется позже. При этом все операции в группе имеют timestamp на пару сотен меньше и выполняются раньше, чем две группы, стоящие перед ней. Логика?  Если в эскизе все элементы определены, то можно поставить еще один размер (хочешь переопределить эскиз - да пожалста!), изменить его - и NX самостоятельно выберет, какой размер или ограничение отпустить и весь эскиз летит к чертям. Это ли не хаос?   Назовете хоть один CAD, которым можно пользоваться так? Любой из встреченных мной имеет много неочевидных нюансов, которые конечно не проявляются при выполнении университетского или базового курса. В Creo я на протяжении 5 лет находил всё новые фишки, которые не прописаны в документации. Огромная благодарность Ruslan за пояснения тонкостей. Даже Компас и Solidworks, которые считаются самыми дружественными, оказалось непросто обуздать, когда проходил собеседования и делал тестовые. А ведь университетский курс знал на отлично!   Возвращаясь к NX - вот сегодняшний пример. Утром деталь была оранжевая, единая. Изменил пару уклонов и толщин, стала желтая с элементами оранжевого и частично твердотельной, частично листовой. Я даже не знаю, как сформулировать вопрос в поддержку. 
    • Mikhail_Podgornyj
      Добрый день, кто может поделиться SKF Bearing V1.6 (ANSYS App)? Заранее благодарен!
    • kkk
      Короче, как сделать автоматом не нашел. Раскрасил нужные грани нужными цветами. Затем на вкладке с цветом модели можно выбирать нужный цвет, тогда в модели выбираются грани, покрашенные этим цветом. Можно еще сохранить выбор этих граней, но смысла не вижу, т.к. все легко выбирается по цвету. Затем либо в инструменте Измерить смотрим площадь и как-то ее используем, либо можно макросом (что выкладывал Kir95 ) внести площадь в свойства (на вкладку Конфигурации). Можно еще внести выбранные грани в датчики, но как использовать их в свойствах не нашел.    Макрос прикрепляю.     Площадь выбранных граней в свойства.swp
    • IgnatvKudrin
      Для нашего производства требуются постоянно инструментальные пружины для штампов. Сейчас заказы обрабатываются довольно долго - до 2-х месяцев. А у нас сроки на штампы 1-2 месяца с испытанием, доработкой и т.п. Хотелось бы найти поставщика, у которого уже есть склад с данными пружинами или кто может доставить за 1 месяц.   Пружины стандартные по ISO, но на всякий случай приложил каталог. Вот список пружин, которые нужны очень срочно:   LR 32 x 44 - 1 шт. LR 25 x 51 - 2 шт. LB 63 x 89 - 7 шт. LB 32 x 102 - 10 шт. LB 32 x 89 - 20 шт. LB 32 x 76 - 12 шт. LB 32 x 51 - 6 шт. LB 25 x 127 - 9 шт.   Готов выслушать все предложения. Возможно даже есть те, кто сможет изготовить подобные пружины с такими же параметрами. 20231101.pdf
    • bubblyk
      Знатоки форума, подскажите, как смоделировать простую (на первый взгляд) задачу термомеханики в упругой постановке.   Есть брусок правильный формы ака кирпич, в котором генерируется НЕРАВНОМЕРНОЕ удельное тепловыделение (HGEN). Задняя и передняя поверхности конвективно охлаждаются с известными (но различными) коэффициентами теплоотдачи и температурами жидкости/газа. Свойства кирпича (теплофизика + механика) принимаются зависящими от температуры. Постановка - стационарная Собственно вопрос : как грамотно смоделировать температурное поле в кирпиче с учетом температурного расширения?   Исходный код приложил.     Надеюсь на коллективный разум и опыт! EPTH2D.mac p.s.   Если есть опыт решения этой задачи в Ansys WB, тоже можно сюда приложить командные вставки.
×
×
  • Создать...