Уважаемые эксперты,
при решении задачи об истечении свободной затопленной струи (осесимметричная постановка, peudo transient, на входе и выходах задано давление) столкнулся с тем, что локальный дисбаланс значительно искажает картину решения. При приемлемом общем уровне невязок решение получается нефизичным. Поясню суть проблемы на иллюстрациях:
(автоматический шаг по времени)
(фиксированный шаг по времени, довольно маленький 1*10^-8с)
(сетка, вроде, нормальная)
(на измельченной сетке)
Вот так влияет на решение
Подскажите, пожалуйста, чем сие может быть вызвано, и как с этим бороться?
(картинки не вставились, повторяю)
Уважаемые эксперты,
при решении задачи об истечении свободной затопленной струи (осесимметричная постановка, peudo transient, на входе и выходах задано давление) столкнулся с тем, что локальный дисбаланс значительно искажает картину решения. При приемлемом общем уровне невязок решение получается нефизичным. Поясню суть проблемы на иллюстрациях:
(автоматический шаг по времени)
(фиксированный шаг по времени, довольно маленький 1*10^-8с)
(сетка, вроде, нормальная)
(на измельченной сетке)
Вот так влияет на решение
Подскажите, пожалуйста, чем сие может быть вызвано, и как с этим бороться?
Здравствуйте!
mesh/check ничего необычного не выдает? Пробовали густить сетку в продольном направлении? Сильно вытянутые ячейки могут давать локально большие невязки.
PS помимо mass imbalance можно через TUI включить для отображения невязки. Команду точно не помню, вроде где-то в solve/set
Здравстуйте!
Большое спасибо за ответ.
Измельчение сетки существенных результатов не дает, вот лучшее, чего удалось добиться
Похоже дело не столько в локальном массовом дисбалансе (абсолютные значения порядка +/-10^-9 кг/с), сколько в параметрах турбулентности.
Так решение по k (Turbulent Kinetic Energy) и Continuity (Mass Flow) имеет сравнительно меньшую устойчивость, что, вероятно, и проявляется на диаграммах, однако уровня 10^-5 добиться удается (см. рис.).
Другими словами, в какой-то момент интенсивность турбулентности внезапно падает практически до нуля, может срабатывают какие-то неочевидные ограничения?
А сетка там у Вас хорошая? Можете показать сами элементы?
Конечно
1653833 quadrilateral cells, zone 2, binary.
3303928 2D interior faces, zone 1, binary.
14 2D pressure-inlet faces, zone 5, binary.
2973 2D axis faces, zone 6, binary.
113 2D pressure-outlet faces, zone 7, binary.
2550 2D pressure-outlet faces, zone 8, binary.
622 2D pressure-outlet faces, zone 9, binary.
1204 2D wall faces, zone 10, binary.
1657572 nodes, binary.
1657572 node flags, binary.
Странно - сетка отличная.
Я не очень большой спец по флюенту но попробуйте другие варианты в Solution Methods -> Gradient.
Большое спасибо, помогло.
Пользуясь случаем хотелось бы спросить, как построить в постпроцессоре контур, например, созданной Custom Field Function?
Заранее благодарен.
Значит дело было в Gradient? Какая опция была изначально и на какой заработало?
Было Least Squares Cell Based,
стало Green Gauss Cell Based,
однако закономерность проследить пока не удается.
Например, при подключении Species Transport получается вот так:
и также Least Squares Cell Based
Добавить комментарий