расчет сопла Лаваля

[HFL 34]

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

Красивая картинка, но лично я вижу небольшую дисперсию за УВ висящей на срезе сопла хотя и говорят что это "дефект визуализации" :)

[Bonusfrag 104]

Красивая картинка, но лично я вижу небольшую дисперсию за УВ висящей на срезе сопла хотя и говорят что это "дефект визуализации" :)

Думаю, что не визуализации.

Думаю, она и есть. Ради этого и выкладывал.

Я (т.к. не спец в моделировании и численных методах) обозвал это "жабрами" и жаловался на подобное.

Но как-то всё не получается нормально разобрать этот момент.

[Bonusfrag 104]

А какие параметры турбулентности у вас были? Или все по умолчанию?

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

и ещё вопрос, что же выходит все таки сетку, построенную адаптацией в данном случае использовать не получится? Никак нельзя адаптировать ячейки только по двум осям?

нет. пользуйтесь сгущением и (максимум) адаптацией 1го порядка на стенке.

когда получите стационар - продатаптируете на скачках. (если потребуется).

и еще - задача валится, когда поток выходит за сопло и огибает кромку. я сейчас считаю этот момент на 0,01 CFL. Пока это помогает.

Т.е. суть - развалилось, откатились назад, на момент до развала, уменьшили CFL и считаете дальше.

И сделайте все НУ/ГУ нормально (например, не надо ставить 1 атм в объеме)

[ksitdikov 0]

когда получите стационар - продатаптируете на скачках. (если потребуется).

Скажите, а можно изменять сетку прямо в процессе расчета?

[Bonusfrag 104]

Скажите, а можно изменять сетку прямо в процессе расчета?

остановил расчет, создал локальную адаптацию на объекте (параллелепипед / конус/ сфера/любое импортированное тело),

и для всех адаптаций общее условие надо включить - "одноразово по шагам" и указать номер шага (следующего за текущим)

а базовая сетка одна на весь счет. её трогать нельзя.

[ksitdikov 0]

остановил расчет, создал локальную адаптацию на объекте (параллелепипед / конус/ сфера/любое импортированное тело),

Спасибо!

[a_schelyaev 366]

Это заблуждение когда речь идёт о сверхзвуковых задачах, связанных с УВ и вообще с разрываными задачами.

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

Спасибо за пищу для размышлений.

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

Я ориентируюсь на внешний вид картинок, которые даёт Tecplot, а в нём визуализация сделана достаточно точно.

Да нет. Я то про конкретную ситуацию в задаче. По картинкам судить не привык. Нужно убедиться что там с визуализацией и сеткой.

Не больше того.

:)

Я конечно не заказчик а скорее наоборот но вопрос чисто ради любопытства

Так сможете ?

Ха! А я как раз думал что за изображение на аватаре.

:)

Крест целовать не готов, но готов попробовать. Тем более, если это течение над клином.

и еще - задача валится, когда поток выходит за сопло и огибает кромку. я сейчас считаю этот момент на 0,01 CFL. Пока это помогает.

Ты на какой версии считаешь?

Релизная?

[Bonusfrag 104]

Ты на какой версии считаешь?

Релизная?

картинка посчитанного сопла, где помогало только CFL < 0.5 была релизная.

сейчас я на 05.04.2012 (т.к. в 25.04.2012 баг, замедляющий счет (+50% времени на итерацию))

Под "валится" я подразумеваю, например, что локальный мах стал >1000 в одной конкретной ячейке

(т.е. до полного развала не стал ждать, снизил CFL)

[a_schelyaev 366]

Чудеса. У меня все штатно было.

Ты задавал НУ так, чтобы стенки сопла их не разрезали?

[Bonusfrag 104]

Проблемная ячейка далеко от того места, где НУ задано в сопле.

Мож, конечно, у нас разные понятия о "штатно"?

невязка по давлению пульсирует от 2000 до 0,002

вот эта ячейка:

скорость (i-й шаг) и ЧислоМаха (i+1 шаг).

Т.е. они чередуются.

Солверная часть с одним шагом:

forFsapr2000.zip

[HFL 34]

Ха! А я как раз думал что за изображение на аватаре.

:)

Крест целовать не готов, но готов попробовать. Тем более, если это течение над клином

Если быть точным это ударная труба с 15 град. клином. Похоже на ютюбовский вариант ссылкой ниже но чуть более сложный. Сначала там получается двойное маховское отражение на клине а потом все это добро бьётся в торец трубы и получаем вот такую вот двойную лямбду на клине с замысловатым каскадом УВ. Конкретно эта картинка считалась в домене 2000х1000 но показана примерно половина

[a_schelyaev 366]

вот эта ячейка:

Это я понял. Раньше торец сопла адаптировали, чтобы эту встречу потоков внутри и снаружи сгладить.

На работе еще раз запущу этот свой вариант.

[ilang 0]

Спасибо, сейчас пробую так делать. Я раньше ставила 1 атм на НУ, потому что так решение не валилось сразу. Но теперь вроде подшаманила и с 0,13 все путем пока.

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

[Bonusfrag 104]

я думаю, вы до них не досчитываете.

покажите распределение плотности хоть для одного досчитанного до конца варианта.

и график давления на срезе сопла от времени (из статуса по характеристике)

UPD:

На работе еще раз запущу этот свой вариант.

Я, кажется, понял, в чем дело. У меня ограничитель по плотности установлен очень маленький.

Поднял его, и стало считать на большем CFL.

[ilang 0]

Вот что насчитала за ночь.

График давлние(время) из маленького параллелепипеда на срезе сопла.

Распределение плотности.

Графики для время=0.0001953, шаг=1953. Дальше решение начало разваливаться, поэтому я вернулась обратно, это видно на графике сходимости.

По поводу ограничителей, не понимаю,почему при увеличении значения можно ставить CFL больше?

[a_schelyaev 366]

По поводу ограничителей, не понимаю,почему при увеличении значения можно ставить CFL больше?

По набору имеющихся данных мы экстраполируем следующую расчетную точку/итерацию.

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

Т.к. в FV применяется схема расщепления, то каждое уравнение фактически живет само по себе.

Стоит одному из них начать "брыкаться" и оно тащит за собой все остальные.

Чтобы этого не происходило использую разные подходы:

1. Снижение шага интегрирования или числа Куранта.

2. Выставить разумные ограничители на физические переменные, чтобы во время расчета они не выбегали за эти ограничители.

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

4. Использовать совместные подходы, когда, например, давление и скорость записаны через одно выражение и решаются в одной матрице.

У нас пока что реализованы только первые три пункта.

Дальше идут хорошие свойства веществ, модели физических процессов, уравнения состояния.

[Bonusfrag 104]

По поводу ограничителей, не понимаю,почему при увеличении значения можно ставить CFL больше?

потому что ограничения не дают выйти значениям за машинные пределы.

Но это касается очень маленькой плотности, которую я поставил вместо вашей. И CFL всё равно выше 1 поднимать нельзя.

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

Как видите - никаких скачков уплотнения нету. просто гладкие распределения и разлет во все стороны...

Там я еще плотность не воткнул. В диапазоне 0-1 кг/м3 плотность выглядит также, как давление на картинке справа снизу на 1й странице.

fsapr2000.pdf

[ksitdikov 0]

Ого! Вот это вещь

Теперь сам ей пользуюсь =)

Я конечно понимаю, что задавать подобный вопрос в этой совсем не корректно с моей стороны, но расчет в FV меня пока не радует стационарностью хоть какой то. Поэтому вопрос по fluent'у: что с ним такое может быть?

Распределение скорости:

а вот и пятно вблизи:

Вот распределение давления:

в близи:

Вот что это такое? Опять проблемы решалки как в FV было с адаптацией сетки? Извините еще раз, и спасибо за ответы)

[Bonusfrag 104]

Поэтому вопрос по fluent'у: что с ним такое может быть?

Вот что это такое?

это высокий шаг по времени. далее с вопросами по fluent в соответствующий раздел переходите.

[a_schelyaev 366]

Я конечно понимаю, что задавать подобный вопрос в этой совсем не корректно с моей стороны, но расчет в FV меня пока не радует стационарностью хоть какой то. Поэтому вопрос по fluent'у: что с ним такое может быть?

Ого! Кто мокрицу в сопло засунул?

:)