Jump to content

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


Guterfreund

Recommended Posts

Guterfreund

Всем привет!

 

Решаю относительно простую задачу: определяю перепад давления на арматуре.

 

На входе задан расход, на выходе постоянное статическое давление, модель k-e, шаг по времени - по умолчанию.

 

Обычно все сходится нормально, но бывает, что не получить качественного решения: на графике перепад давления-итерации синус, невязки по давлению и скорости также скачут между e-4 и e-5, уходя в горизонт.

 

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

 

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

 

Спасибо! 

Link to post
Share on other sites


UnPinned posts
karachun

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

 

Может образовался вихрь который ведет себя нестационарно?

 

Если перепад давления сошелся то можно считать что решение тоже сошлось, хоть и плоховато.

 

Можно для очищения совести просчитать в нестационаре и осреднить результат по времени.

  • Нравится 1
Link to post
Share on other sites
Guterfreund
48 минут назад, karachun сказал:

в постпроцессоре вывести их контуры

А если проблема внутри объема? Переменные, насколько я понимаю, выводятся на границе среды.

 

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

Если перепад давления сошелся то можно считать что решение тоже сошлось, хоть и плоховато.

Ясно. Но перепад давления в моем решении скачет в пределах 5-10%. Можно, конечно, усреднить, но хотелось бы понять как сделать правильно. 

Link to post
Share on other sites
karachun
3 минуты назад, Guterfreund сказал:

Переменные, насколько я понимаю, выводятся на границе среды.

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

 

3 минуты назад, Guterfreund сказал:

но хотелось бы понять как сделать правильно. 

Усреднить после нестационара. Это будет правильно.

Это как нахождение к-та сопротивления при обтекании цилиндра.

Link to post
Share on other sites
Guterfreund
2 минуты назад, karachun сказал:

Усреднить после нестационара. Это будет правильно.

Вероятно. Но решаемая задача вполне стационарная. Можно все усложнить, но с практической точки зрения это не удобно и не эффективно.

 

Мне кажется, есть способы решения такой задачи в стационарной постановке.

 

Проверю, где высокие невязки и поработаю с сеткой в той области.

 

Link to post
Share on other sites
piden
8 minutes ago, Guterfreund said:

Но решаемая задача вполне стационарная.

Что-то советы Карачуна мимо пролетели(

 

@Guterfreund , можешь и шаг по времени задрать (все равно в псевдо-транзиенте решаешь) - и тоже получишь сходимость. Только толку от нее? Чтобы графики невязок глаз радовали?

 

1 hour ago, Guterfreund said:

Решаю относительно простую задачу: определяю перепад давления на арматуре.

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

Spoiler

output.gif

 

  • Нравится 1
Link to post
Share on other sites
Guterfreund
4 минуты назад, piden сказал:

Чтобы графики невязок глаз радовали?

Нет такой цели. Хочу перепад давления - горизонтальную линию.

 

4 минуты назад, piden сказал:

Что-то советы Карачуна мимо пролетели(

Ничего не пролетает мимо ;). Просто задача на потоке и времени в обрез.

Link to post
Share on other sites
karachun

@Guterfreund Так вы определитесь - вы хотите решить правильно или вы хотите решить побыстрее?

Я бы вообще забил на осциляции если у них амплитуда 5-10%. Это не так много. Взял бы и тупо осреднил.

21 минуту назад, Guterfreund сказал:

Просто задача на потоке и времени в обрез.

Ну так тем более.

Просто не обращать внимание на эти осциляции и делов.

 

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

Но это не будет правильно. Вы просто сделаете красивую картинку и более красивые графики сходимости.

Link to post
Share on other sites
Guterfreund
2 минуты назад, karachun сказал:

Так вы определитесь - вы хотите решить правильно или вы хотите решить побыстрее?

Я бы вообще забил на осциляции если у них амплитуда 5-10%. Это не так много. Взял бы и тупо осреднил.

Скорее всего я и буду осреднять. 

 

Я не могу назвать себя опытным CFD расчетчиком. Меня просто интересует, почему 90% моделей сходится на ура, а 10% имеет проблемы. Будем считать, что это спортивный интерес. Если фантазировать, то ответ мог бы быть таким: товарищ, ты не знал, что такая плохая сходимость для таких задач решается включением такой-то галки в настройках сольвера. Стыдно, товарищ.

Link to post
Share on other sites
piden
9 minutes ago, Guterfreund said:

Если фантазировать

Зачем фантазировать, если уже был дан ответ:

 

2 hours ago, karachun said:

образовался вихрь который ведет себя нестационарно

1 hour ago, karachun said:

как .... при обтекании цилиндра

 

 

Link to post
Share on other sites
karachun
7 минут назад, Guterfreund сказал:

Меня просто интересует, почему 90% моделей сходится на ура, а 10% имеет проблемы.

От острых кромок, и от разных торчащих в потоке пиптиков получаются вихревые дорожки а они имеют нестационарную природу.

 

Вон, в теме про сужающее устройство

https://cccp3d.ru/topic/110410-расчёт-сужающего-устройства-су/

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

2020-07-29_20-37-13.png

sketch.jpg

  • Нравится 1
Link to post
Share on other sites
Guterfreund

@karachun пересчитал с выводом невязок. С ходу не могу понять, как в постпроцессоре вывести невязки по объему. На Default Domain Default я могу это прицепить, но выводятся лишь границы среды, а объема нет :(

Link to post
Share on other sites
karachun
6 минут назад, Guterfreund сказал:

С ходу не могу понять, как в постпроцессоре вывести невязки по объему.

Isoclip (или может isosurface) или Volume Rendering

19 минут назад, Guterfreund сказал:

товарищ, ты не знал, что такая плохая сходимость для таких задач решается включением такой-то галки в настройках сольвера

Это прямо влажные фантазии, что придет человек из Кемерово и молча поправит все.

Увы это физика а от нее кнопки еще не придумали.

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

Edited by karachun
Link to post
Share on other sites
soklakov
12 часов назад, Guterfreund сказал:

Я не могу назвать себя опытным CFD расчетчиком.

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

  • Нравится 1
Link to post
Share on other sites
Guterfreund

@karachun , не подскажете, есть ли возможность в CFX настроить монитор на отображение перепада давления между двумя произвольными сечениями, а не между Input и Output? В постпроцессоре как сечения вставить понятно, а в препроцессоре не найти такой кнопки :(

 

Спасибо!

Link to post
Share on other sites
karachun

@Guterfreund Вроде так нельзя. В препроцессоре нужно опираться на существующие поверхности. Можно взять и собрать трубу из нескольких тел, поверхности на границе тел можно добавит в именованную выборку и тогда возможно ее удастся выбрать для монитора.

 

А просто произвольную поверхность - роде нельзя.

 

Может еще @piden подскажет.

Link to post
Share on other sites
karachun

@Guterfreund В CFX Pre есть некая User surface. Вроде она должна помочь. Я помню что для одной задачи мне тоже нужно было еще в процессе расчета следить за одним параметром в точке внутри объема жидкости и я сделал малюсенькую треугольную User surface и мониторил давление на ней.

Edited by karachun
  • Нравится 2
Link to post
Share on other sites
piden
17 minutes ago, karachun said:

В CFX Pre есть некая User surface. Вроде она должна помочь.

+1

Создать в CFD-Post, потом экспортировать в Pre.

 

17 minutes ago, karachun said:

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

Для такого ведь вроде Monitor Point должно хватить?

  • Нравится 2
Link to post
Share on other sites
karachun
5 минут назад, piden сказал:

Для такого ведь вроде Monitor Point должно хватить?

Наверное я тоже так думал но не смог Monitor Point настроить и вместо этого использовал маленькую User Surface.

Edited by karachun
Link to post
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.
Note: Your post will require moderator approval before it will be visible.

Guest
Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

  • Recently Browsing   0 members

    No registered users viewing this page.



  • Сообщения

    • Aer
      поуродовался и построил через "сру курвес мех"....
    • kex1570
      Напишу в эту тему, не хочу новую плодить. Кто-нибудь сталкивался с тем, что лицензия на Ansys Student, которая выдаётся вроде как на год(Renewable, twelve-month lease.), заканчивалась раньше? Мб можно её продлить? спасёт ли переустановка? или нужно подчищать следы за собой? Ansys 21-й, последний.
    • Aer
    • Alecs_st
    • Nastyrnuy
      1  Итак, я не понимаю как можно быть толковым электриком(электронщиком) если не знать основы блоков питания которые обеспечивают стабилизированное питания различными напряжениями(+/-5 В/ и +-12В), ну плюс еще на всякий случай включить туда контроль тока для защиты( чтобы каждый раз не портить плату с выпайкой и перепайкой..). Эти схемы в учебной доке рассписаны от и до..! Ну есть конечно некоторые извраты связанные в свое время с ленью чинит БП Персоналок, так как проще было типа купить новый и помощнее например на 350-450 Вт вместо 200-250 ватников(да были времена соотношения 1 $ -33рубь...!). И у меня есть такие знакомые( и даже классные) на ардуине соберут ЧПУ, но БП чинить не хотят(тем более если нет схем,,). Но есть и другие которым не в лом и БПки чинить(так как ничего зазорного в этом не видят..подумаещь схем нет, а должно так работать и все. Главное есть откуда плясать Есть вход есть выход все и этоого достаточно!) так что скажем либо ленятся, либо не толковые, либо стимуляции нет.!! 2- помнится я говорил что у Ажи идет последовательная эволюция плат, и что  на ранних стадиях на станках вместо 20 LPS  были другие... 15  гугления(  я просто знал что и как искать..) вытыщил аналог функциональный и конструктивный LPS-03! Можете на него не только полюбоваться, но и сравнить элементную базу, и схемы поискать на него  и саму плту поискать.. дерзайте..
    • Viktor2004
      Мне нужен физический контакт. И время
    • gudstartup
      На картинке похож на ПО станков компании neway  
    • Viktor2004
      Когда мне попадется, обязательно посмотрю. Пока вживую не видел
    • Viktor2004
      Давайте перейдем сюда. А то там чужая тема   Ладно. Мне все равно долго вникать надо что к чему
    • gudstartup
      Вот китайцы например выпускают пульты оператора с аппаратной блокировкой на форуме была такая темка!
×
×
  • Create New...