диагоз срыва расчета

[an_rushin 128]

у меня вопрос скорее общего порядка

работаю, как правило, с очень большими моделями порядка 80-100 млн. ячеек. сетка строится в режиме "быстро & грязно", т.е. качество даже не плохое, а ужасное, скьюнес > 0.999 и сеточный генератор едва ли в состоянии улучшить качество. пользуюсь связкой Fluent - Harpoon (Sharc). часто случается, что расчет срывается из-за качества сетки. угадать эпицентр проблемы в такой модели помогут только интуиция и опыт, но что делать, если даже эти коварные спутники молчат? как бы вы поступили, диагностируя причину срыва? как бы искали? перезапустить расчет, выписать результат за несколько итераций до срыва, но место с очень высокой скоростью или очень низкой плотностью или... и .т.д. далеко не показатель настоящей причины. какую практическую рекомендацию вы бы дали?

[Bonusfrag 104]

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

во FlowVision сейчас расчет при развале сам останавливается, и проблема как правило находится в центре области развала (область неадекватных значений).

Там начинаю шаманить с сеткой или уменьшать шаг по времени.

 

Сетку-то как строите? полностью на автомате?

[HFL 34]

Задача/тип решателя хотя бы какие ? (сжимайка/несжимайка)

[an_rushin 128]

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

во FlowVision сейчас расчет при развале сам останавливается, и проблема как правило находится в центре области развала (область неадекватных значений).

Там начинаю шаманить с сеткой или уменьшать шаг по времени.

 

Сетку-то как строите? полностью на автомате?

 

 

с демонстрацией дело плохо, все, что делаю, абсолютно секретно.

сетка строится на полном автомате из cut-cell.

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

Задача/тип решателя хотя бы какие ? (сжимайка/несжимайка)

 

задача стационарная, решатель COUPLED (Fluent), среда сжимаемая, т.к. теплопередача и излучение.

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

[HFL 34]

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

 

 

У cut-cell известная болезнь это непредсказуемый размер резаных ячеек. Как следствие - Курант в "обрезках" может превышать всякие разумные пределы (хотя Fluent и умеет уменьшать Куранта автоматом но это скорее "предсмертная агония" решения)

[a_schelyaev 366]

 

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

 

 

У cut-cell известная болезнь это непредсказуемый размер резаных ячеек. Как следствие - Курант в "обрезках" может превышать всякие разумные пределы (хотя Fluent и умеет уменьшать Куранта автоматом но это скорее "предсмертная агония" решения)

 

 

В Струегляде на этот случай есть механизм слития "маленьких" ячеек с крупными соседями. Т.о. вопрос запредельно малого Куранта отпадает в начале, а в конце все это дело мы лакируем нашим всеядным линейным решателем, который успешно съедает такие разбросы по числу Куранта.

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

[HFL 34]

Недаром Саров скакал вокруг да около.

 

 

Зато у них с баблом всегда кучеряво  

 

А скакать вокруг того что им понравилось эт они умеют   

[a_schelyaev 366]

 

Зато у них с баблом всегда кучеряво

 

 

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

 

 

 

 

А скакать вокруг того что им понравилось эт они умеют    

 

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

[an_rushin 128]

 

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

 

 

У cut-cell известная болезнь это непредсказуемый размер резаных ячеек. Как следствие - Курант в "обрезках" может превышать всякие разумные пределы (хотя Fluent и умеет уменьшать Куранта автоматом но это скорее "предсмертная агония" решения)

 

 

 

 

ага... тогда зайдем с другого бока. что посоветуете, чтобы стабилизировать расчет? по крайней мере, повысить стабильность. в режиме coupled пробовал курант от 15 до 1, результат одинаков, ну или почти одинаков.

[Bonusfrag 104]

вы спрашиваете про черный ящик. поэтому отвечать всем тяжело.

 

расчет когда падает? в самом начале или в середине в произвольный момент? 

ну и не помешало бы всё-таки найти какой-то несекретный пример с проблемой.

Да хоть от руки какую-то близкую железку нарисуйте... и получите падение.  

 

Без объекта надо быть совсем ясновидящим, чтоб угадать решение для вас.

 

можно не сюда выложить, а на какой-нить гугл-диск на пару дней картинку, тут дать ссылку.

потом через день-два удалите ссылку и объект с диска.

Можно выделенным людям в личку ссылку кинуть.

 

Мне можно не кидать, я флюентом никогда не пользовался, и сетки там от FV отличаются.

[HFL 34]

 

 

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

 

 

У cut-cell известная болезнь это непредсказуемый размер резаных ячеек. Как следствие - Курант в "обрезках" может превышать всякие разумные пределы (хотя Fluent и умеет уменьшать Куранта автоматом но это скорее "предсмертная агония" решения)

 

 

 

 

ага... тогда зайдем с другого бока. что посоветуете, чтобы стабилизировать расчет? по крайней мере, повысить стабильность. в режиме coupled пробовал курант от 15 до 1, результат одинаков, ну или почти одинаков.

 

 

Общие рекомендации:

 

 -  1й-порядок дискретизации по пространству

 - ограничители по параметрам сверху-снизу (температура, давление, etc.)

 

PS: Тепло как вводится в систему ? Объемный источник или как то еще ?

 

C "резаной" сеткой  это только перемешить ну или использовать фичи мешеров как описано выше для FV

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

 

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

[a_schelyaev 366]

Я бы еще посоветовал на сегрегейт перейти.

[German_psv 20]

1. Для уменьшения скошенности следует уменьшить размер ячеек.

Сделайте сетку из небольших тетров.

Во Флуенте преобразуйте тетру в полихедру,

это в несколько раз сократит размер сетки и

несколько уменьшит остроту проблемы скошенности.  

 

2. Включите быстрый нестационар NITA.  

Полихедра+NITA будет решать не медленнее, чем в стационаре.

Кроме того, PISO автоматически корректирует скошенность.

 

3. Используйте переменный (во времени) шаг солвера.

Для этого: Time Stepping Method - Settings - Variable 

Опции переменного шага: число куранта 0.15, минимальный шаг 1e-08, остальное - по смыслу.

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

 

В этом случае срыва вычислений не будет

[a_schelyaev 366]

Курант 0,15???

[German_psv 20]

Transient, NITA - лучше с такого куранта начинать.

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

курант всегда можно поднять.

[a_schelyaev 366]

На дворе 21 век, а они все с Курантом меньше 1 считают.

[Bonusfrag 104]

На дворе 21 век, а они все с Курантом меньше 1 считают.

))) вы сами недавно исправились (для СЗ задач)

[HFL 34]

На дворе 21 век, а они все с Курантом меньше 1 считают.

Гиперзвук, детонацию и горение в нестационаре считают с Курантом < 1 чуть менее чем всегда   

[a_schelyaev 366]

 

На дворе 21 век, а они все с Курантом меньше 1 считают.

))) вы сами недавно исправились (для СЗ задач)

 

Извините, но всю жизнь с Курантом 1 считали, а не меньше 1.

 

На дворе 21 век, а они все с Курантом меньше 1 считают.

Гиперзвук, детонацию и горение в нестационаре считают с Курантом < 1 чуть менее чем всегда 

 

 

В общем, чудеса в решете.