Отражение скачков уплотнеия

[dmiantosha 11]

Задача проста: в небольшом тунельчики обдувается сверхзвуком цилиндр (боковая поврехность), необходимо посчитать отражение скачков уплотнения от верхней и нижней стенки тунеля (квадрат). Если на снизу и сверху - Openning, все отлично, если снизу одна стенка - резульатат получаю, если две стены всерху и снизу решение разваливается. Как поступить в этой ситуации?

[BigBrother 1]

Картинку в студию. Геометрию и полученные решения. И параметры сетки тоже

[dmiantosha 11]

Вот сетка и картинка. Слева Inlet->Supersonic->NoremalSpeed=600 m/s

->Pressure=1 atm

Серху и снизу должны бать стены.

Справа Opening->StaticPressure=1 atm;

Результатов нет, так как нет решения.

[Dick 0]

А стенки с проскальзыванием ставили (free slip)?

[dmiantosha 11]

Нет, no slip, что надо free?

[BigBrother 1]

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

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

Ещё можно поставить не Auto Timescale, а Physical Timescale c величиной шага раз в 10 меньше, чем cfx давал при Auto Timescale

Попробуйте на выходе поставить сверхзвуковой outlet. А что на входе - это вопрос. Я не спец в этой области, но может быть происходит запирание канала? Тогда ставить на входе скорость нельзя. Попробуйте на входе задать некоторое избыточное давление. Несколько атмосфер.

Изменено 17 февраля 2007 пользователем BigBrother

[dmiantosha 11]

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

[Dick 0]

На стенках надо поставить проскальзывание, так как падающий скачок на стенке вызовет отрыв ПС. Когда стоит openning - отрыва не происходит, а имеет место отражение СУ. Вот вероятная причина неудач.

[dmiantosha 11]

Пробывал и такие и секие стены, решение не получаю.

[BigBrother 1]

Обождите. Что значит цилиндр летит? Он летит в этом канале?

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

Но действительно, тут проще поставить free slip. А если он летит в канале, то и нужно его ставить.

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

Все возможные комбинации ГУ на входе равнозначны. Только не надо наверное было задавать избыточное статическое давление 1 атм.

И ещё: The Opening boundary condition type is only available for Subsonic boundaries.

Можно попробовать расширить канал (или уменьшить диаметр цилинра :). Можно начать с небольшой скорости и увеличивать её.

[BadArtem 9]

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

Те с одной стенкой все-таки решает.

Задачка то симметричная когда стенка сверху и снизу - ну и надо обрезать ее по оси симметрии. Я думаю, все будет Ок.

[dmiantosha 11]

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

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

[BigBrother 1]

>Я думаю, все будет Ок.

Не будет.

Знакомые картинки

>Плохо, что нельзя сделать подвидным (реально) цилиндр, этот опыт хотелось бы повторить в CFX.

Что значит нельзя? Необходимо локализовать проблему, а для этого можно менять всё что потребуется. А там видно будет.

Изменено 18 февраля 2007 пользователем BigBrother

[BadArtem 9]

Кто здесь говорил, что не получится?

Этот рисунок демонстрирует, что расчет проходит - не более того.

Все остальное (проверка корректности решения и т.п.) к dmiantosha.

Изменено 18 февраля 2007 пользователем BadArtem

[BigBrother 1]

Дык! Раз решается половинка, тогда и целая будет. А если целая не решается, тогда и половина не будет

[BadArtem 9]

..и еще: (уже более серьезно)

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

Кто-нибудь сталкивался с подобной проблемой? И как решали ее? Этот эффект есть и в представленной выше картинке.

"А если целая не решается, тогда и половина не будет" - не согласен.

Кстати, на рисунке из атласа - ШАР. Его не корректно сравнивать с цилиндром.

[dmiantosha 11]

Кстати, на рисунке из атласа - ШАР. Его не корректно сравнивать с цилиндром.

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

to ВadArtem, так ты просто половинку взял и все, ГУ такие же как у меня?

[Krusty 199]

Знакомые картинки clap_1.gif

Могу выложить в рапиду, кому надо

[BadArtem 9]

dmiantosha, inlet: U=600 m/s, V=0, W=0, P=50000 Pa, T=100 K

Air Ideal Gas

насчет экватора - не думаю, что будет то же самое - вспомни расчет числа М на оси сопла (там зависимость от F) для плоского и осесимметричного получаются разные результаты

Насчет цилиндра в сверхзвуковом потоке - есть шлирная фотка в Чжене, том кажется II, страницы .... где то в начале.

Изменено 18 февраля 2007 пользователем BadArtem

[Tako 0]

dmiantosha, а что у Вас в направлении, перпендикулярном экрану?

Надеюсь, один слой Хекс и на обеих поверхностях симметричные или периодические ГУ.

У Вас на стенках СЛИШКОМ большие элементы, вот при таких скоростях и не сходится. Там пристеночный слой - самый ответстветный, а Вы пытаетесь все эффекты в один ряд элементов уместить? Попробуйте на стенках сделать рядов 5-10 очень мелких элементов. Лучше, если и на цилиндре тоже сделаете.

Вы же мне присылали пример по сверхзвуковому обтеканию крыла. Посмотрите там сетку вблизи крыла и сравните с Вашей.

Изменено 19 февраля 2007 пользователем Tako