Max_7o

Расчет сетчатого фильтра

60 сообщений в этой теме
Pinned posts

Добрый день!

Коротко о том какая задача стоит. Есть фильтр и в него входит: сам корпус со входом и выходом и фильтрующий элемент (в моем случае это металлический перфорированный лист свернутый в цилиндр отверстиями 4 мм). Для этого фильтра надо рассчитать перепад давления. Максимально допустимое давление в системе известно 4 МПа, так же есть максимальный расход через всю систему в которой установлен фильтр это 41,6667 кг/с (перевел из другой единицы). В ходе гидравлических испытаний было установлено, что при чистом фильтрующем элементе перепад составляет 0,03 МПа. Среда - вода. 

По фильтрующему элементу: моделировать я каждое отверстие не стал, использовал пористое тело с такими настройками: volume porosity 0.4 (если я правильно понял, это площадь отверстий 40%). Далее по формуле CR2  = (K loss *p) /B2 нашел Quadratic Resistance Coefficient, у меня он получился 22,725 кг/м4. 

По ГУ: задал максимальный известный расход 41,6667 на входе; на выходе давление Averadge Statiс Pressure - 0. В результате расход на выходе  - минус 0,02 кг/с, а на входе -  3,2....кг/с.

Далее попробовал задать давление на входе 4 МПа, выход остался как и в первом варианте. Результат: давление на выходе упало до 60-70 кПа...

Далее я изменил пористое тело на Fluid, чтобы посмотреть как будет меняться расход. Гу: на входе расход 41,6667, на выходе Averadge Statiс Pressure - 0. Результат: на выходе расход  - минус 0,00798 кг/с; на входе - 0,109 кг/с.

Модель - SST.

Ansysом пользуюсь недавно. Может я в ГУ ошибку допустил. Не могу никак разобраться. Помогите пожалуйста)

 

Поделиться сообщением


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


UnPinned posts

Когда ты задал на входе 4 МПа а а выходе ноль, то и перепад получился 4 МПа, а по результатам испытаний 0,03 МПа, вот это давление и надо задавать как давление на входе. Для несжимаемой среды абсолютное давление не играет роли. По несовпадению расходов - очень странно, если ты задал на входе 40 кг/с то то 3,2 кг/с там быть уже не может.

С пористыми телами не работал, попробуй для начала вообще выкинуть пористое тело и посчитать сам корпус фильтра, добейся для начала баланса по массе. В сети можно найти туториалы по CFX, выполни пример 12. Flow in a catalytic Converter - там используется пористое тело. Еще можешь поискать конструкцию своего фильтра в справочнике Идельчика, там можно посчитать к-т потерь для всего фильтра, не Quadratic Resistance Coefficient а просто дзетту.

Ну и стандартные вопросы - делал ли ты сеточную сходимость (еще одна ссылка https://flowvision.ru/ru/support-menu-header-ru/blog-ru/ob-issledovanii-skhodimosti-po-setke) и достаточно ли далеко у тебя вход и выход от зоны обратных течений? Если выход слишком близко то в зависимости от типа ГУ это либо приведет к проблемам с решением либо воду будет подсасывать из выхода (если сделать выход типом Opening).

Поделиться сообщением


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

ID: 3   Опубликовано: (изменено)

27 минут назад, karachun сказал:

Когда ты задал на входе 4 МПа а а выходе ноль, то и перепад получился 4 МПа, а по результатам испытаний 0,03 МПа, вот это давление и надо задавать как давление на входе. Для несжимаемой среды абсолютное давление не играет роли. По несовпадению расходов - очень странно, если ты задал на входе 40 кг/с то то 3,2 кг/с там быть уже не может.

С пористыми телами не работал, попробуй для начала вообще выкинуть пористое тело и посчитать сам корпус фильтра, добейся для начала баланса по массе. В сети можно найти туториалы по CFX, выполни пример 12. Flow in a catalytic Converter - там используется пористое тело. Еще можешь поискать конструкцию своего фильтра в справочнике Идельчика, там можно посчитать к-т потерь для всего фильтра, не Quadratic Resistance Coefficient а просто дзетту.

Ну и стандартные вопросы - делал ли ты сеточную сходимость (еще одна ссылка https://flowvision.ru/ru/support-menu-header-ru/blog-ru/ob-issledovanii-skhodimosti-po-setke) и достаточно ли далеко у тебя вход и выход от зоны обратных течений? Если выход слишком близко то в зависимости от типа ГУ это либо приведет к проблемам с решением либо воду будет подсасывать из выхода (если сделать выход типом Opening).

Спасибо за ответ. 

По поводу сеточной сходимости есть ошибки P-mass находится ниже 10е-4, а как я вычитал из приведенной статьи такого быть не должно. Сейчас сделаю сетку более мелкой и уберу фильтрующий элемент.

На выходе outlet и я еще увеличил трубу на 5 диаметров от входа и выхода в фильтр.

Изменено пользователем Max_7o

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
26 минут назад, Max_7o сказал:

На выходе outlet

Можешь попробовать заменить его на Opening, через этот тип ГУ поток может и вытекать и втекать и для процесса решения это полезно.

28 минут назад, Max_7o сказал:

По поводу сеточной сходимости есть ошибки P-mass находится ниже 10е-4

Может выше? Но вообще это не так критично и смотреть надо не только на ресидуалс. https://cae-club.ru/publications/kogda-ostanavlivat-reshenie-v-cfx

Как я понял значение Quadratic Resistance Coefficient бралось из потолка. Можно сделать небольшое исследование - смоделировать малый участок этого фильтра, скажем 10х10 отверстий, можно даже сделать его плоским. добавит входной и выходной каналы а боковые стенки сделать симметрией. Вычисленный к-т сопротивления можно сравнить со значениями по Идельчику, если оба значения совпадут, то это добавит уверенности и будет намного точнее чем пальцем в небо.

35 минут назад, Max_7o сказал:

На выходе outlet и я еще увеличил трубу на 5 диаметров от входа и выхода в фильтр.

Нормально, если сетка позволяет, то можно увеличить до 20 диаметров, но это уже после того как сам расчет станет сходится, для повышения точности. Посмотри по результатам - поток должен успеть выпрямится, вихревое течение на выходе это плохо.

И по входным данным - перепад 0,03 МПа и макс. расход 41,67 кг/с могут быть на разных режимах работы и не совпадать, желательно найти расход именно во время испытаний.

Поделиться сообщением


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

ID: 5   Опубликовано: (изменено)

47 минут назад, karachun сказал:

Можешь попробовать заменить его на Opening, через этот тип ГУ поток может и вытекать и втекать и для процесса решения это полезно.

хорошо, попробую.

47 минут назад, karachun сказал:

Как я понял значение Quadratic Resistance Coefficient бралось из потолка

Не совсем, нашел на просторах порядок решение данного коэффициента.  Формулы во вложении. Не могу сказать правильные они или нет.

 

47 минут назад, karachun сказал:

И по входным данным - перепад 0,03 МПа и макс. расход 41,67 кг/с могут быть на разных режимах работы и не совпадать, желательно найти расход именно во время испытаний.

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

 

47 минут назад, karachun сказал:

Можно сделать небольшое исследование - смоделировать малый участок этого фильтра, скажем 10х10 отверстий, можно даже сделать его плоским. добавит входной и выходной каналы а боковые стенки сделать симметрией.

Как раз начал проводить данное исследование и наткнулся на формулы упомянутые выше, решил попробовать.

 

47 минут назад, karachun сказал:

Может выше?

По моему кривая стремилась к 10е-5. Могу сказать точно она отдалялась от других, тогда как остальные уже пересеклись и двигались в одном диапазоне.

 

Безымянный.png

Изменено пользователем Max_7o

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
2 часа назад, Max_7o сказал:

По моему кривая стремилась к 10е-5.

Это хороший показатель для невязок.

А вообще можно попробовать вариант по проще - не моделировать пористое тело а создать на месте фильтра интерфейс и задать к-т сопротивления на нем. Сам коэффициент можно найти по диаграмме 8-3 в справочнике https://dwg.ru/dnl/3036.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
В 01.02.2019 в 17:01, karachun сказал:

А вообще можно попробовать вариант по проще - не моделировать пористое тело а создать на месте фильтра интерфейс и задать к-т сопротивления на нем. Сам коэффициент можно найти по диаграмме 8-3 в справочнике https://dwg.ru/dnl/3036.

Я так и делаю. Назначил на фильтрующий элемент интерфейсы и указал что это пористое тело. Но вот к сожалению с диаграммой никак не могу разобраться.

В 01.02.2019 в 17:01, karachun сказал:

Сам коэффициент можно найти по диаграмме 8-3

данный коэффициент будет соответствовать  Quadratic Resistance Coefficient ? 

 

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

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
42 минуты назад, Max_7o сказал:

И возник вопрос, что задать на входе, чтобы увидеть какой перепад дает сам корпус?

Можно задать массовый расход, для которого был определен данный перепад.

43 минуты назад, Max_7o сказал:

данный коэффициент будет соответствовать  Quadratic Resistance Coefficient ?

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

54 минуты назад, Max_7o сказал:

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

Напиши параметры своей пластины - диаметр отверстий, толщина, живое сечение, расход через фильтр и габариты самой пластины (высота/ширина/диаметр).

58 минут назад, Max_7o сказал:

Я так и делаю. Назначил на фильтрующий элемент интерфейсы и указал что это пористое тело.

Мое предложение такое - заменить пористое тело одним интерфейсом нулевой толщины (пластина скорее всего тонкая), и а этом интерфейсе можно задать к-т потерь из Идельчика, не для пористого тела.

Вчера искал материалы по своей задаче и наткнулся на несколько тем на CFD-Online, у перфорированной пластины к-т сопротивления будет зависеть еще и от угла между потоком и пластиной, это можно учесть свойствами пористого тела но пока это задача со звездочкой. К ней можно будет приступать когда фильтр начнет считаться.

https://www.cfd-online.com/Forums/fluent/27490-perforated-plates.html

https://www.cfd-online.com/Forums/cfx/114670-perforated-sheet-wiremesh-resistance-coefficient.html

https://www.cfd-online.com/Forums/cfx/115973-perforated-sheet-hexagonal-holes.html

Поделиться сообщением


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

ID: 9   Опубликовано: (изменено)

22 минуты назад, karachun сказал:

Напиши параметры своей пластины - диаметр отверстий, толщина, живое сечение, расход через фильтр и габариты самой пластины (высота/ширина/диаметр).

Пластина: 150*150 мм (если я правильно понял, тебе нужна пластина с помощью которой я исследование провожу, а не фильтрующий элемент (пластина свернутая в цилиндр)), толщина 2 мм. Отверстия а - 4 мм, t - 6 мм. На картинке показано расположение. Расход через сам фильтр известно только максимальный (до сих пор не получил данные испытаний) - 41,6667 кг/с. 

22 минуты назад, karachun сказал:

Мое предложение такое - заменить пористое тело одним интерфейсом нулевой толщины (пластина скорее всего тонкая), и а этом интерфейсе можно задать к-т потерь из Идельчика, не для пористого тела.

Насколько я знаю для интерфейса нужен и вход и выход, т.е. один интерфейс это выход из корпуса и вход в фильтрующий элемент (цилиндр) и 2ой выход из цилиндра в корпус.

rv-form.jpg

Изменено пользователем Max_7o

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
В 01.02.2019 в 12:52, karachun сказал:

По несовпадению расходов - очень странно, если ты задал на входе 40 кг/с то то 3,2 кг/с там быть уже не может.

Сейчас закончился расчет перфорированного листа 150*150 (отверстия я моделировал). на входе задал расход 41,6667кг/с, на выходе Averadge Statiс Pressure - 0. Проверяю расход, на входе 0,0308 кг/с на выходе вообще минус 0,0012 кг/с. Может я не правильно результат смотрю? А смотрю так: выбираю Contour и в поле Locations - in/out; Variable выбираю Mass Flow; Range - Local.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
20 минут назад, Max_7o сказал:

Может я не правильно результат смотрю?

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

https://www.cfd-online.com/Forums/cfx/24529-how-monitor-mass-flow-rate-cfx.html

По расчету из диаграммы приложил таблицу. Дальше создаешь один интерфейс (выбираешь одну поверхность именованной выборкой), в свойствах выбираешь Interface Model -> Pressure Cgange. Чтобы перепад давления зависел от дин. напора надо писать выражение.

5c594eaeec637_1.JPG.233a145ebf4a30c8d7f0125501f29420.JPG

5c594eaf5fd29_2.JPG.86f6ed57f65aa553c93037ea61b49852.JPG

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

-7.4*1000*(massFlow()@REGION:interface/1000/area()@REGION:interface)^2/2

7.4 - к-т по Идельчику, 1000 - плотность воды, конструкция в скобках это расчет скорости как объемный расход деленный на площадь. Знак минус чтобы это было падение давления. Сама формула тоже взята из Идельчика, вот из такой формулы считается перепад давления.

5c595030cc327_.JPG.9bc9eeb4789dfcd6ee539451e3da06de.JPG

Также надо проверить число Рейнольдса для отверстий в плите - диаграмма 8-3 применима для Re>= 10e5. Возможно при таком расходе режим течения сильно отличается и все мои расчеты это филькина грамота.

Книга1.xlsx

Изменено пользователем karachun

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
3 часа назад, karachun сказал:

Дальше создаешь один интерфейс (выбираешь одну поверхность именованной выборкой)

Одну поверхность? Так разве можно? я думал у интерфейса есть свой вход и выход. Я так понимаю ты предлагаешь, убрать сам фильтр (он у меня в модели присутствует как отдельное тело) и назначить интерфейс на стенки которые имеют между собой расстояние в 2 мм или вообще избавиться от этих стенок и смоделировать одну? Скрин модели прикрепил.

За расчет отдельное спасибо.

 

3 часа назад, karachun сказал:

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

С расходом разобрался, спасибо)

Безымянный.png

Изменено пользователем Max_7o

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
1 час назад, Max_7o сказал:

или вообще избавиться от этих стенок и смоделировать одну

Вот это я имел в виду. Просто одна поверхность с перепадом давления, управляемым выражением.

И мой вариант и пористое тело это упрощения, при этом использовать один интерфейс проще и легче контролировать результаты.

Изменено пользователем karachun

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
17 часов назад, karachun сказал:

Просто одна поверхность 

Может конечно я чего то не знаю, но весь домен должен быть твердым телом...поверхность в DesignModeler не экспортировать. Вычесть только поверхность мне Solid не позволяет.

Я пытался в DM избавиться от лишних поверхностей, но так тоже не получается.

 

Изменено пользователем Max_7o

Поделиться сообщением


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

Какая у тебя версия ансиса?

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
1 минуту назад, karachun сказал:

Какая у тебя версия ансиса?

19

Поделиться сообщением


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

Вот пример, сетка вворкбенче, проект CFX в отдельном архиве.

Test.wbpz

Unnamed.zip

Посмотрел на сове выражение и оно мне не нравится, в нем берется средняя скорость для всей фильтрующей пластины но ее геометрия сложная и по этому такое предположение будет достаточно грубым. Надо смотреть скорости по отдельным элементам, но как это сделать выражением я пока не знаю (да и времени свободного у меня сейчас нет).

В общем не надо считать мой подход истинным а использование пористого тела неправильным. Просто я никогда не использовал porous domain и под рукой у меня только справочник Идельчика.

И еще - какой получился перепад давления для корпуса фильтра, без пластины.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
2 часа назад, karachun сказал:

В общем не надо считать мой подход истинным

Я попробую повторить то что ты сделал с одной поверхностью.

Теперь по давлениям, как я понимаю я снова не правильно смотрел результат. Я делал это через Contour и давление только на входе было правильным. 

На входе задал Static pressure - 4МПа, на выходе Opening c характеристиками opening pres. and Dirn - 0 Па ( даже не знаю почему этот параметр выбрал)

Результат, написал несколько выражений в expressions:

- areaAve(pressure)@in - 3.99904e+006 Па

- areaAve(pressure)@Out - 229.886 Па

- areaAve(pressure)@in-areaAve(pressure)@Out - 3.99881e+006 Па (это я вычитал из СFD, там это добавляли в монитор)

Изменено пользователем Max_7o

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
31 минуту назад, Max_7o сказал:

На входе задал Static pressure - 4МПа, на выходе Opening c характеристиками opening pres. and Dirn - 0 Па

0 Па это если выход в окружающую среду а у тебя там дальше сеть. При таком перепаде (40 атм) вода из этого фильтра будет вылетать со свистом, я бы напротив этого фильтра не стоял.

Еще раз. Тебе известно давление в сети 4 МПа, известен перепад на фильтре 0,03 МПа. Следовательно за фильтром давление 4-0.03=3.97 МПа, не ноль, разница огромна. Дальше начинаются особенности CFD если задать давление 4 МПа то на его фоне жалкие 0.03 МПа могут затеряться (все малые флуктуации давления съедаются округлением). Есть два варианта.

1) для сжимаемых сред задают "опорное давление" в него включают основную часть давления, если это течение в окружающей среде то берут атмосферу, у тебя же очевидно оно будет 4 МПа. В ГУ, на входе и выходе указывают давление без опорного 0.03 и 0 МПа соответственно. Свойства газа считают с опорным а сам расчет идет без него, течение определяет разница давлений а не абсолютные величины. Но все это только для расчета газов, тебе сейчас это не надо.

2) вода несжимамая, мы считаем что ее свойства не зависят от давления (в реальности зависят но мы можем эти пренебречь), получается что опорное давление и не нужно (по умолчанию в CDF-Post выводится давление без опорной части, все равно надо рыться в настройках чтобы его показать). Можно задать опорное 4 МПа, давления вход/выход 0.03/0 МПа или опорное 0 а на ГУ 0.03/0 или даже опорное 0 а на ГУ -3.97/-4 МПа (даже отрицательные, математика все стерпит). На расход это не повлияет, важна только разница. Но ошибки округления никто не отменял.

По этому я рекомендую на входе задать 0.03 МПа, на выходе 0, опорное 0 (для несжимаемых расчетов это правило хорошего тона). А на картинках результатов твердой рукой надписывать "для получения абсолютных величин давления прибавить 3.97 МПа".

Можно в принципе задать опорное и 3.97 МПа, в результатах расчета есть вроде бы "Absolute Total Pressure" и это полное давление+опорное. Но это не точно, надо разбираться.

Вообще понимание опорного давления это важная часть изучения CFD, без него начинающие часто путаются в результатах.

31 минуту назад, Max_7o сказал:

opening pres. and Dirn ( даже не знаю почему этот параметр выбрал)

Нормально, можешь им пользоваться и дальше.

 

33 минуты назад, Max_7o сказал:

это я вычитал из СFD, там это добавляли в монитор

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

Изменено пользователем karachun

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
19 минут назад, karachun сказал:

Вообще понимание опорного давления это важная часть изучения CFD, без него начинающие часто путаются в результатах.

C этим я действительно еще не разобрался.

В CFX это reference pressure? 

По задаче: на входе (Static или Total) pressure 0,03 МПа, на выходе opening pres. and Dirn - 0 Па? Сначала попробовать без фильтра?

 

22 минуты назад, karachun сказал:

Вообще в CFX множество задач делаются чрез выражения. Если присмотришься то почти все входные данные можно задать в виде выражений

А есть какой нибудь мануал по всем (ну хотя бы по основным) выражениям?

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
51 минуту назад, Max_7o сказал:

reference pressure

Да, на входе Total (stable) 0.03 МПа. По всему остальному тоже "да".

Выражения в CFX описаны в хелпе (там ссылки внизу). Обычно в том месте из которого скачивают ансис идет отдельный образ с документацией.

https://www.sharcnet.ca/Software/Ansys/16.2.3/en-us/help/cfx_ref/i1304009.html

https://www.sharcnet.ca/Software/Ansys/16.2.3/en-us/help/cfx_ref/i1304139.html

Много вопросов и ответов по выражениям есть на CFD-Online, если владеешь английским то можешь там зарегистрироваться и задавать вопросы.

Когда создаешь выражение жмешь правой кнопкой мыши на пустом поле и у тебя появляются выпадающие списки для самих функций (areaAve), переменных и объектов сетки, так меньше шансов сделать ошибку и сразу видно что тебе доступно.

Вообще по CFX есть еще такие ресурсы:

каналы с вебинарами:

https://www.youtube.com/channel/UCA-Xy7AVqmwEKX9uqMgELxg

https://www.youtube.com/user/CADFEMCIS

особенно вот это видео, это вообще маст хэв)

https://www.youtube.com/watch?v=Q6fMnTN5fvs

https://www.youtube.com/user/Ivanrips

https://www.youtube.com/channel/UCYbowzUPfsH3rA3srE618YA

есть еще https://cae-club.ru/ у них есть блог с полезными статьями с разными трюками в CFX и форум.

на cdf-online есть хорошая статья по разным вопросам

https://www.cfd-online.com/Wiki/Ansys_FAQ

И еще такая задача для самостоятельного изучения: для разных моделей турбулентности необходима разная высота первой ячейки сетки от стенки, это определяется параметром y+ (игрек плюс). Если ты еще не изучил этот вопрос то посмотри достаточно ли измельчена твоя сетка? Информация по необходимым значениям есть в документации к CFX (theory guide или как то так) и множество материалов есть в сети и на CFD online. Можно начать поиск с запроса "cfx sst yplus". В твоей задаче это значение может и не очень важное но в общем случае это тоже фундаментальная часть CFD.

https://www.cfd-online.com/Tools/yplus.php

Изменено пользователем karachun

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
8 минут назад, karachun сказал:

Да, на входе Total (stable) 0.03 МПа. По всему остальному тоже "да".

Выражения в CFX описаны в хелпе (там ссылки внизу). Обычно в том месте из которого скачивают ансис идет отдельный образ с документацией.

https://www.sharcnet.ca/Software/Ansys/16.2.3/en-us/help/cfx_ref/i1304009.html

https://www.sharcnet.ca/Software/Ansys/16.2.3/en-us/help/cfx_ref/i1304139.html

Много вопросов и ответов по выражениям есть на CFD-Online, если владеешь английским то можешь там зарегистрироваться и задавать вопросы.

Когда создаешь выражение жмешь правой кнопкой мыши на пустом поле и у тебя появляются выпадающие списки для самих функций (areaAve), переменных и объектов сетки, так меньше шансов сделать ошибку и сразу видно что тебе доступно.

Вообще по CFX есть еще такие ресурсы:

каналы с вебинарами:

https://www.youtube.com/channel/UCA-Xy7AVqmwEKX9uqMgELxg

https://www.youtube.com/user/CADFEMCIS

особенно вот это видео, это вообще маст хэв)

https://www.youtube.com/watch?v=Q6fMnTN5fvs

https://www.youtube.com/user/Ivanrips

https://www.youtube.com/channel/UCYbowzUPfsH3rA3srE618YA

есть еще https://cae-club.ru/ у них есть блог с полезными статьями с разными трюками в CFX и форум.

на cdf-online есть хорошая статья по разным вопросам

https://www.cfd-online.com/Wiki/Ansys_FAQ

И еще такая задача для самостоятельного изучения: для разных моделей турбулентности необходима разная высота первой ячейки сетки от стенки, это определяется параметром y+ (игрек плюс). Если ты еще не изучил этот вопрос то посмотри достаточно ли измельчена твоя сетка? Информация по необходимым значениям есть в документации к CFX (theory guide или как то так) и множество материалов есть в сети и на CFD online. Можно начать поиск с запроса "cfx sst yplus". В твоей задаче это значение может и не очень важное но в общем случае это тоже фундаментальная часть CFD.

https://www.cfd-online.com/Tools/yplus.php

Спасибо тебе большое за материалы)

 

9 минут назад, karachun сказал:

Обычно в том месте из которого скачивают ансис идет отдельный образ с документацией.

да, видел, в pdf файлы если не ошибаюсь.

По поводу видео, буду по возможности дома изучать. На работе все заблочено.

15 минут назад, karachun сказал:

для разных моделей турбулентности необходима разная высота первой ячейки сетки от стенки, это определяется параметром y+ (игрек плюс)

По поводу этого параметра, я где то читал. На стенке при построении сетки я задал пограничный слой в 5 уровней(слоев).

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
16 часов назад, karachun сказал:

Да, на входе Total (stable) 0.03 МП

На выходе opening - 0 Па; Опорное давление - 0. без фильтрующего элемента.

 

Вот что получилось:

- areaAve(Totalpressure)@in - 30005,5 Па

- areaAve(Totalpressure)@Out - 9625,21 Па

- areaAve(totalpressure)@in-areaAve(Totalpressure)@Out - 20380.2 Па 

Поделиться сообщением


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

Попробуй использовать Pressure вместо TotalPressure.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
13 минуты назад, karachun сказал:

Попробуй использовать Pressure вместо TotalPressure.

В выражениях? или снова расчет запустить с pressure вместо Total?

Изменено пользователем Max_7o

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
В 05.02.2019 в 16:09, karachun сказал:

Просто одна поверхность с перепадом давления

C моделью пока тоже не получается. Я закинул в DM фильтр без фильтрующего элемента. С помощью команды Extrude вытянул поверхность и закрыл ее другой поверхностью. В итоге у меня стало 3 тела. Закинул в mesh, сначала я попробовал строить сетку сразу на 3х телах, при переходе в CFX выдал ошибку: 

An error occured initializing CFX-Pre: Application error(s) occurred in CFX-Pre.  Please see the messages view for further details.
: Failed to export model.
There was a problem importing the mesh from the requested file.
The importing process reported the following problem:
  Model "Model" does not have a valid mesh associated with it.
 
  Body B126 is an unsuppressed surface body that does not enclose
  a volume. Only 3D volume meshes can be imported. If your model
  contains valid bodies as well as surface bodies then suppress
  the surface bodies and remesh before continuing.

 

Далее я попробовал создать для каждого тела сетку отдельно и вот что опять выдал:

 

An error occured initializing CFX-Pre: Application error(s) occurred in CFX-Pre.  Please see the messages view for further details.
It was not possible to extract any mesh from the file:
D:\Workbench\Filter_one_surface_files\dp0\SYS\MECH\SYS1.cmdb.
Please note that a Workbench installation is normally required to import this type of file.
 
Вывод: я вообще не умею с поверхностью работать в Ansys CFX...
 
А в том примере что ты скинул, я не могу понять у этих пластин есть толщина или нет? Если нет, то как ты их сделал solid?
Изменено пользователем Max_7o

Поделиться сообщением


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

В выражениях, то есть смотришь на статическое давление (Pressure) а не на полное (Total Pressure).

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

Как то так.

clip_image002.jpg

По ошибкам - по тексту мало что понятно, я с такими ошибками не сталкивался, лучше скинь архив проекта.

Изменено пользователем karachun

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
1 час назад, karachun сказал:

В выражениях, то есть смотришь на статическое давление (Pressure) а не на полное (Total Pressure)

Может тогда на входе задать статическое, а не полное?

Вот результат по давлению:

- areaAve(pressure)@in - 20762.1 Па

- areaAve(pressure)@Out - 1.06789 Па

- areaAve(totalpressure)@in-areaAve(Totalpressure)@Out - 20761.1 Па 

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

И как я уже говорил делал еще вариант с разными сетками для тела (solid) и для поверхностей, которые сам построил в DM.

Filter_one_surface.wbpz

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
31 минуту назад, Max_7o сказал:

на входе задать статическое

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

https://www.sharcnet.ca/Software/Ansys/16.2.3/en-us/help/cfx_mod/i1300679.html

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

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

Проект посмотрю, но позже.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
23 минуты назад, karachun сказал:

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

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

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

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
47 минут назад, karachun сказал:

И еще смотри результаты - нет ли вихрей у выхода

По поводу вихрей, подскажи как это сделать. С помощью Streamline? Если да, то все в порядке.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
2 часа назад, Max_7o сказал:

В дальнейшем будут еще фильтры, но уже с другими фильтрующими элементами.

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

2 часа назад, Max_7o сказал:

По поводу вихрей, подскажи как это сделать. С помощью Streamline?

Да линии тока, или вывести вектора или контур скорости.

Вот твой проект. Надо сперва добавить тело и поверхности в один парт а потом добавить Share topology (в опциях надо выбрать Allow Solid/Surface Parts - No). Тогда объем рассечется на две части но сетка будет общая (так как тела в одном парте).

И еще несколько замечаний, но только с точки зрения перфекционизма, как говорят автомеханики "на скорость это не влияет".

5c5c40101391b_4.thumb.JPG.2365bb685d344782b57546074646da30.JPG

1. Выкинуть из расчета этот маленький пиптик, там все равно выхода нет, только сетку лишнюю жрет.

2. Немного подрезать геометрию чтобы убрать острый угол, там все равно скорости низкие.

3. Скруглить края пластины, принять радиус скругления как 1/3-1/2 толщины пластины, она все равно в зоне интереса и сетка на ней подробная, без острых кромок решение будет сходится немного легче.

4. Хотел предложить сделать фаску 10-15 мм по всему контуру но это уже излишний перфекционизм.

Filter_one_surface_1.wbpz

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
В 07.02.2019 в 17:35, karachun сказал:

Надо сперва добавить тело и поверхности в один парт а потом добавить Share topology (в опциях надо выбрать Allow Solid/Surface Parts - No).

Я пытался назначить интерфейс на поверхность которая получилась Filter_surf (из проекта который ты скинул), но снова не получается. При назначении интерфейса нужны 2 поверхности или я все таки просто не умею это делать. Пытался назначить на этот цилиндр 2 поверхности через NS, CFX ругается, говорит что это одна и таже поверхность (что и следовало ожидать). 

И вопрос по 1му проекту который ты скидывал. В какой момент у тебя появились вот эти 2 интерфейса( test_interface,test_interface Shadow) в DM и Mesh их вообще нет и в CFX они появились?

В 05.02.2019 в 11:53, karachun сказал:

5c594eaeec637_1.JPG.233a145ebf4a30c8d7f0125501f29420.JPG

 

Поделиться сообщением


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

Скинь модельку с простой сеткой, мне проще показать на пальцах чем объяснять.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
30 минут назад, karachun сказал:

Скинь модельку с простой сеткой, мне проще показать на пальцах чем объяснять.

Любую модель? или именно фильтр со сгенерированой сеткой?

На фильтре у меня не сказал бы что простая сетка...

Поделиться сообщением


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

Построил сетку на той модели, которую ты скидывал раньше, действительно вторая поверхность с приставкой "_shadow" не появляется. По чему так происходит я не знаю. Но есть способ обойти это ограничение. Надо по очереди подавить сперва одно тело а потом другое и экспортировать две половинки сетки в разные файлы и потом закинуть их в CFX по отдельности. Тогда можно создать интерфейс между двумя сетками, на точность это не должно повлиять - сетка ведь узел к узлу. 

Filter_one_surface_2.zip

Возможно виноват Share Topology, но это не точно.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
1 час назад, karachun сказал:

Построил сетку на той модели, которую ты скидывал раньше, действительно вторая поверхность с приставкой "_shadow" не появляется. По чему так происходит я не знаю. Но есть способ обойти это ограничение. Надо по очереди подавить сперва одно тело а потом другое и экспортировать две половинки сетки в разные файлы и потом закинуть их в CFX по отдельности. Тогда можно создать интерфейс между двумя сетками, на точность это не должно повлиять - сетка ведь узел к узлу. 

Filter_one_surface_2.zip

Возможно виноват Share Topology, но это не точно.

Еще вопрос возник, а цилиндр весь получается перфорированный будет? просто дно у цилиндра в модели без отверстий. 

А если выбрать интерфейс, то весь цилиндр подсвечивается, или так задумано?

Поделиться сообщением


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

Выбери только боковую поверхность в Named selection. Тогда торец будет по умолчанию стенкой, но надежнее задать его стенкой вручную.

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

Поделиться сообщением


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

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

Также можно использовать симметрию и моделировать только половину задачи. Выгода очевидна.

5c6181a22c862_.thumb.JPG.da7cfc946083db12382bef950b7b6647.JPG

Filter_one_surface_3.wbpz

Еще касательно y+. В зависимости от выбранной пристеночной функции (wall function) требования будут разными. В таком случае можно выбрать automatic wall treatment и не делать сверх подробный подслой, сопротивление фильтра все равно будет зависеть в основном от перепада на интерфейсе/свойств пористого тела.

https://www.cfd-online.com/Forums/cfx/71606-yplus-sst.html#post242005

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах
16 часов назад, karachun сказал:

сопротивление фильтра все равно будет зависеть в основном от перепада на интерфейсе/свойств пористого тела.

Я вчера по проекту который ты скинул успел только сетку построить, но при этом запустил расчет с пористым доменом. Прошло всего 270 итераций, и я решил посмотреть что получилось. Я в монитор добавлял кривую перепада давления, заметил что она в пределах от 7 кПА до 8 кПа колеблется и остановил. Да и сразу скажу, что с опорным давлением пока не разобрался). Но когда я его задавал 0 атм для расчета с пористым телом, решатель падал в ошибку Overflow. Когда я вернул опорное давление 1 атм, расчет пошел. Но вот с результатом что то не так или я опять смотрю не правильно. Посмотри пожалуйста файл во вложении. Сравнивается результат с фильтрующим элементом и без. Но в варианте без фильтрующего элемента отсутствует опорное давление. Может попробовать его запустить с опорным? Или может еще надо итераций добавить?

Результат фильтр.xlsx

Изменено пользователем Max_7o

Поделиться сообщением


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

Создайте аккаунт или войдите для комментирования

Вы должны быть пользователем, чтобы оставить комментарий

Создать аккаунт

Зарегистрируйтесь для получения аккаунта. Это просто!


Зарегистрировать аккаунт

Войти

Уже зарегистрированы? Войдите здесь.


Войти сейчас

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

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