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

Перепад давления на участке прямой гидравлически гладкой трубы


Guterfreund

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

Еще вопрос: ставлю для невязок значение 10^-5. Когда решение уже почти сошлось, невязки начинают расти (все). Т.е. до 100 итерации невязки дружно падали и почти достигли 10^-5, а потом начинается их рост и следующие 50 итерации они дружно растут. На 200-ой итерации у меня не хватило терпения и я тормознул решение.

 

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

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


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

https://cae-club.ru/publications/kogda-ostanavlivat-reshenie-v-cfx

Если задал маленький шаг по времени то возможно это все только начало и поток не установился.

Получается что если невязки большие то надо смотреть на качество сетки.

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

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

https://cae-club.ru/publications/kogda-ostanavlivat-reshenie-v-cfx

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

 

В 11.12.2018 в 12:40, karachun сказал:

Получается что если невязки большие то надо смотреть на качество сетки.

На качество в прямом смысле или на ее размер?

 

Вот еще думаю, я задаю Reference Pressure всегда 0 Па для своих расчетов. Пишут, для воды не должно быть критичным. Не может тут быть заковырки? Вот мой график сходимости на большом расходе и мелкой сетке (см. ниже). Данная задача на меньшем расходе и более крупной сетке (в соответствии с требуемым Y+) сходилось без проблем.

Monitor.jpg

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

На качество в прямом смысле или на ее размер?

На качество, покрученная сетка дает плохое решение.

4 часа назад, Guterfreund сказал:

Вот еще думаю, я задаю Reference Pressure всегда 0 Па для своих расчетов. Пишут, для воды не должно быть критичным. Не может тут быть заковырки?

Не должно, у воды по умолчанию нет свойств, зависимых от давления.

Выложите проект, может ошибка закралась в другом месте.

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

покрученная сетка дает плохое решение.

Можно поподробнее, что есть такое покрученная сетка?

 

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

Ну в смысле сетка с элементами плохого качества - острые углы, большое соотношение сторон (но не для подслоя, там поток направлен параллельно элементу и все хорошо).

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

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

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

 

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

 

Спасибо!

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

Подскажите, как в постпроцессоре отобразить плохие по качеству элементы

качество элементов можно посмотреть еще в пре-

на этапе генерации сетки.

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

качество элементов можно посмотреть еще в пре-

Ок, тогда копаем глубже. На данный момент я генерирую сетку стандартным модулем Mesh. Если не говорить о различного рода лопатках достаточно ли функционала данного модуля для построения качественных сеток для CFD задач или есть куда стремиться и оптимальнее использовать другие генераторы?

 

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

 

По качеству сетки: практически всегда имею некоторое количество плохих элементов по Mech Expansion Factor, вижу это в Output файле решателя. В окне Mesh, такого показателя качества сетки напрямую не нашел, поэтому и спросил про постпроцессор.

 

И чтобы окончательно разобраться. В Output файле решателя я вижу, что есть некоторое количество плохих элементов, допустим менее 1%. Насколько я понимаю, глобально на качество решения такое количество оказать влияние не может. Допустим, я смотрю перепад давления, эта величина усредненная по сечению, качество отдельных элементов навряд ли окажет влияние на данную величину. Какой же процент плохих элементов тогда допусти, 5, 10? Или же хорошим тоном является вообще отсутствие плохих элементов во всем объеме модели?

 

Спасибо!

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

Насколько я понимаю, глобально на качество решения такое количество оказать влияние не может.

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

7 часов назад, Guterfreund сказал:

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

сеточную сходимость никто не отменял.

7 часов назад, Guterfreund сказал:

На данный момент актуален выбор между гексаэдрической и тетраэдрической сетками?

 

 

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

сеточную сходимость никто не отменял.

По сеточной сходимости для CFD расчетов: имеется, допустим, труба. Я хочу определить коэффициент гидравлического сопротивления. Насколько я понимаю, для конкретного расхода я не смогу оценить сеточную сходимость? Размер пристеночного элемента определяется величиной Y+, т.е. скоростью и, в данном случае, диаметром. От размера пристеночного элемента пляшет вся сетка в поперечном сечении. Т.е. для конкретного Y+ я не могу замельчить сетку и проверить на ней сходимость. Таким образом, либо увеличивать Y+ в пределах, характерных для модели турбулентности, либо увеличивать, например, скорость и, автоматически, тем самым, менять размер пристеночного элемента.

 

Интересуют практические рекомендации по анализу сеточной сходимости для такого вида расчетов...

Изменено пользователем Guterfreund
Ссылка на сообщение
Поделиться на других сайтах

Если вы на стенке получаете сетку с Y+ не более 1, то дальше вам остается только измельчать сетку в объеме.

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

По сеточной сходимости для CFD расчетов

а что если посмотреть на это так: побить трубу кубиками или элементами с Aspect Ratio близкими к 1. уменьшать размер кубиков, пока не будет достигнута сходимость. никаких погранслоев, сплошные кубики.

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

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

Тут проблема вот в чем: я работаю с моделью турбулентности K-E, для которой Y+ должен быть в пределах 30-300. Моя стартовая величина Y+ = 40. Получается, если я на следующем этапе мельчу сетку, то выхожу на Y+<30, что влечет за собой либо переход на другую модель турбулентности, например SST c Y+ = 1-2 (насколько я понял эту тему), либо в пределах старой модели турбулентности я, вероятно, получу неточные результаты уже из-за достаточно мелкого для данной модели размера пристеночного элемента.

Ссылка на сообщение
Поделиться на других сайтах
  • 1 месяц спустя...
В 24.12.2018 в 10:32, Guterfreund сказал:

Тут проблема вот в чем

мдя... шляпа.

может взять стартовый у+ в районе 300? и проверить, что на пути к 30 сходимость будет наблюдаться уже в районе 100, а может раньше.

или отключить модель турбулентности) тем более труба прямая, можно скорость снизить до ламинарного течения.

но это так, в порядке бреда.

Ссылка на сообщение
Поделиться на других сайтах
В 19.02.2019 в 18:07, soklakov сказал:

может взять стартовый у+ в районе 300

можно попробовать, но при десяти пристеночных элементах как бы не схлопнулся диаметр пристеночной области, а еще и ядро хотелось бы заиметь :). Все-таки для K-E размер первого элемента уже исчисляется миллиметрами, а не их тысячными долями...

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • AlexArt
      Ну допустим, ты и на другом ресурсе это опубликовал. А не коммуниздил. Но вот продвигать воровство от государства, ворующее из Вики, это верх мерзости.
    • maxx2000
      Ах, да. Фильтры выбора добавили. Теперь можно выбрать только то что видно на первом плане, а не вместе с тем что с обратной стороны детали. В общем надо обновляться. Как раз работёнка на прессформу нарисовалась 
    • maxx2000
      Причина того - Кроилово. Кроилово всегда приводит к попадалову. Месяц простоял сколько мильонов деревянных потеряли? Вопрос риторический. И ещё будет стоять. Как памятник человеческой глупости и жадности.
    • AlexKaz
      "9 июля 1968 года на мышах был проведен самый знаменитый эксперимент американского ученого-этолога Джона Кэлхуна «Вселенная-25». Суть опыта заключалась в создании идеальных условий, где мыши могли бы жить и размножаться, не ведая никаких забот, вдали от хищников и в отсутствие эпидемий и заболеваний. Для этих целей ученый построил специальный загон, куда были помещены четыре пары белых мышей (самцов и самок). В распоряжении мышей всегда была чистая вода и еда в изобилии, специальные гнезда, где можно обустроить себе жилище ― гнезд в загоне хватало для проживания нескольких тысяч мышей. Температура в загоне в среднем составляла около 20 ℃ и была комфортной для мышей. Животные не подвергались никаким влияниям извне и жили в идеальных условиях в свое удовольствие. А дальше началось самое интересное. На первом этапе эксперимента мыши хорошо размножались, вели активный образ жизни, охотно играли. На следующей фазе эксперимента мыши стали есть меньше, перестали наедаться до отвала. На третьей фазе эксперимента, когда в загоне были уже сотни мышей, произошло распределение социальных ролей, стала ярко выраженной иерархия, клановость. Появились так называемые отверженные ― молодые особи, которых другие, взрослые мыши сгоняли в центр загона, не давали им вести нормальный образ жизни, причиняли физический вред. В природе такое, наверное, было бы невозможно, ведь эти мыши-агрессоры просто не дожили бы до старости: их бы съели хищники. Но в загоне Кэлхуна хищников не было, и взрослые мыши начали попросту издеваться над молодняком. Образовались две большие группировки: самцы-одиночки и самки-одиночки. При этом самки-одиночки отказывались спариваться <с менее статусными многочисленными молодыми самцами и с оставшимися старыми статусными> и отвергали ухаживания самцов. У мышей стал проявляться тотальный индивидуализм, мыши не стремились создать семью. На последней, четвертой стадии мышиная популяция стала сокращаться. Появились самцы, которых сам Кэлхун назвал «красивыми» (англ. beautiful ones), из-за отсутствия ран и рубцов. <В оригинале: They never engaged in sexual approaches toward females, and they never engaged in fighting, and so they had no wound or scar tissue. Thus their pelage remained in excellent condition. - Дословный перевод: Они никогда не прибегали к сексуальным подходам к самкам, и они никогда не участвовали в боях, и поэтому у них не было ран или рубцовой ткани. Таким образом, их шерсть сохранилась в отличном состоянии.> Эти мыши не вступали в борьбу за самок и территорию, не проявляли активности к размножению и только питались, спали и чистили шёрстку. У мышей стали проявляться различные формы девиантного поведения, вспышки агрессии. Самки стали проявлять агрессию, защищать себя сами, стали умерщвлять своих детенышей, а затем окончательно отказались размножаться. На пике эксперимента в загоне одновременно проживало чуть более двух тыс. мышей. Еды и гнезд было достаточно для дальнейшего роста популяции, но через четыре года после начала эксперимента Кэлхун остановил свой опыт, потому что в загоне осталось чуть более сотни мышей, и все они уже вышли из репродуктивного возраста. По итогам эксперимента Кэлхун пришел к выводу, что достижение определенной плотности населения и заполнение социальных ролей в популяции приводит к распаду общества" https://physicsoflife.pl/dict/pic/calhoun/calhoun.. https://scientificrussia.ru/articles/utopiya-dlya-mys.. https://ru.wikipedia.org/wiki/Кэлхун,_Джон_(этолог)
    • gudstartup
      @Koels вот в чем дело пока ds609 это предупреждение поэтому F может и не появится если sv601 это значит ошибка. возможно при нагреве радиатора серво определяет это как предупреждение или ваш вентилятор крутиться медленнее чем оригинальный и серва думает что он встал хотяпри этом обычно на экране в строке состояния FAN.мигает больше у меня вариантов нет....  
    • ДОБРЯК
      Решите любым алгоритмом. Тогда будет конструктивный разговор. :=)
    • Fedor
      https://en.wikipedia.org/wiki/List_of_numerical_analysis_topics#Eigenvalue_algorithms     :) 
    • Юрий К.Ф.
      Добрый день. Не нашёл тут тему по стойке Китайско Китайской)) Lynuc N3ME. Видать мне так повезло с её наличием)) Приобрели 5-ти осевой Китаец. В б/у состоянии после удара по оси Z. Отремонтировали по механике, заменили батареи на драйверах, выставили лимиты. Всё Ок. Пригласили со стороны людей которые бы разобрались по операторской части. Те два выходных ковырялись, после сказали покажут расскажут, но за огромные деньги. Не сошлись. После месяц станок простоял, когда включили перестал реагировать на регулировку скорости шпинделя. То есть в режиме Jog, включаем обороты, которые стандартно 2140-2149 об/мин. При регулировке процетности не меняются (сама процентность показывает на мониторе). Так же при включении оборотов через команду M03S300 или другое значение, скорость так же показывает 2140-2149 об/мин. Грешить на тех товарищей с которыми не сошлись по деньгам для обучения, как то не хочется. Поковырялся в настройках шпинделя, вроде всё в норме. Проводку на шпинделе прозвонил, целая. В чём причина, не понятна. Кто нибудь сталкивался с подобным, или с подобной стойкой? Может подсказать варианты причины подобного?
    • ДОБРЯК
    • Fedor
      https://en.wikipedia.org/wiki/Eigenvalue_algorithm   Вот , кстати, важное свойство которое раньше встречал у Стренга. :)    http://elibrary.bsu.edu.az/files/kitablar/1022.pdf   
×
×
  • Создать...