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

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


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

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

post-26004-1336391868_thumb.png

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

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


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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Спасибо!

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

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

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

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

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

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

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

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

:)

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

Так сможете ?

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

:)

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

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

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

Релизная?

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

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

Релизная?

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

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

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

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

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

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

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

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

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

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

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

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

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

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

post-26004-1336416340_thumb.pngpost-26004-1336416349_thumb.png

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

forFsapr2000.zip

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

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

:)

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

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

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

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

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

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

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

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

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

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

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

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

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

UPD:

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

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

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

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

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

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

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

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

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

post-35213-1336546717_thumb.jpg

post-35213-1336546732_thumb.jpg

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

fsapr2000.pdf

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

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

post-30625-1336583665_thumb.jpg

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

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

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

post-30625-1336583877.jpg

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

post-30625-1336583916.jpg

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

post-30625-1336583936.jpg

в близи:

post-30625-1336583949.jpg

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

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

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

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

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

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

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

:)

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

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

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

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

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

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

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

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

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

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

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




×
×
  • Создать...