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

[a_schelyaev 366]

Занятная задачка.

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

Сейчас новый решатель все это проглатывает не кашляя.

В чем проблема. Проблема в том, что волна давления доходит до сужения сопла, где площадь сечения критики меньше по площади чем оставшаяся часть сечения. ВОЛна отражается от стенки конфузора и идет обратно. Доходит в виде ударной волны до входа, где задано давление и с этого места начинается "Страх и ненависть в Лас-Вегасе".

Все существующие методики, которые были тут перечислены, годятся для дозвукового истечения с ГУ.

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

Если есть оценка этой скорости, то нужно вместо ГУ Давление задаь Давление с Нормальной скоростью. И не забыть указать скорость.

Продожаем считать.

[Bonusfrag 104]

всё-таки я бы не отказывался от ГУ "давление"...

просто поднять его надо плавно - P=P(t)

поднять настолько плавно, чтоб не было ударной волны

[a_schelyaev 366]

Поднимал уже.

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

[a_schelyaev 366]

Все, я зашел в тупик.

А что вообще за устройство? Можно по подробнее описать, как осуществляется подача газа в сопло?

[Bonusfrag 104]

А что вообще за устройство? Можно по подробнее описать, как осуществляется подача газа в сопло?

...посмотри файл-отчет, который я посылал в конце года (Сопло.pdf)...

Не думаю, что это сопло от моего принципиально отличается.

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

вообще, заинтриговали) люблю проблемки) сейчас ковырну задачку)

[Bonusfrag 104]

Проект прилагаю.

в этом проекте мне бросилось в глаза:

1) температура вакуума = 300К .... многовато. поставил 3К

2) вернул опорные значения в 300К@100000Па (соотв-но в пространстве заданы -99999,8Па@-297К)

3) свободный выход : заменил сверхзвук на статическое давление (равное НУ). статическое давление автоматом переключится на сверхзвук в нужных ячейках.

4) вход/выход - Статическое давление на входе заменил на Полное давление. Статику можно задать только на входе (т.е. если вдруг пойдет за обратной волной масса, то ничего хорошего не будет). И еще момент - статика там меняется (на мой взгляд) хитрее, чем полное давление, поэтому, по идее, полное - корректнее.

5) перекроил геометрию, т.к. выглядела странно...

попробую посчитать)

UPD по п.4.:

"Статическое давление на входе" и "Полное давление" - это дозвуковые ГУ.

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

либо шаманить, меняя тип ГУ.

Вариантом тут является Нормальная скорость с давлением. Это ГУ допускает и вход и выход, на сверхзвуке близко к прибиванию скорости и давления, на дозвуке близко к прибиванию давления. В целом, логика зашита странная (не очевидная) - через массовую скорость, но в д.с. это может сыграть на руку.

в общем, пока оставил на нем считать. поставил маленькую скорость и рост давления.

[a_schelyaev 366]

4) вход/выход - Статическое давление на входе заменил на Полное давление.

Это можно делать, только если вход дозвуковой.

Все остальное я уже пробовал.

в общем, пока оставил на нем считать. поставил маленькую скорость и рост давления.

Я уже пробовал.

:)

[ilang 0]

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

[Bonusfrag 104]

Я уже пробовал.

А плотность таблично задавать?

Начальная плотность очень высокая получается.

М.б. стоит её понизить, задав таблично?

[a_schelyaev 366]

Я пока сделал финт ушами.

Поднял опорное давление до 10 Па, а шаг сверху ограничили 1e-9.

Картина течения другая. Теперь поток плавно доходит до конфузора и от него не отражается, а вдоль него образуется особая пристеночная зона течения.

Сейчас снял ограничение по шагу интегрирования и он вырос до 5e-8.

как получу струю, то опорное давлению опущу до исходного и посмотрим что получилось.

[Bonusfrag 104]

А разве можно менять опорное давление между шагами??

Мне говорили, что так нельзя делать.

а то, что шаг надо ограничивать - это понятный момент. я бы еще рекомендовал CFL 0.45 ставить.

и еще - там сужение резкое. в итоге у стенки вектор направлен вниз. Но ячейки вытянуты вдоль Х.

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

[a_schelyaev 366]

Сетку я уже давно переделал. Она у меня внутри сопла везде квадратная - выше изображение приводил.

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

Если перестало считать, то значит звоночек.

P.S. Да нет, один хрен отражается и бьет во вход.

[Bonusfrag 104]

P.S. Да нет, один хрен отражается и бьет во вход.

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

[a_schelyaev 366]

Это нормально. Как раз раньше мы тут сетку адаптировали.

Сейчас решалка все это проглатывает.

[Bonusfrag 104]

Это нормально. Как раз раньше мы тут сетку адаптировали.

Сейчас решалка все это проглатывает.

фишка в том, что на эти несколько ячеек не подействовала, например, установка переменной.

в общем, они стремные...

пробую одновременно поднимать давление на входе и сразу во всем докритическом объеме (с помощью установки переменной)

т.е. и там, и там одна и та же формула стоит.

тогда, по идее, не будет волны.

[a_schelyaev 366]

Мы не знаем какой режим должен быть на входе.

Мальчики ждут девочек - показаний автора.

[Bonusfrag 104]

Мы не знаем какой режим должен быть на входе.

Мальчики ждут девочек - показаний автора.

формально, да.

реально нужен стационар выхода )

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

по идее, это баллон с клапаном.

[Bonusfrag 104]

Это нормально. Как раз раньше мы тут сетку адаптировали.

Сейчас решалка все это проглатывает.

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

скорость большая, ячейка маленькая => наименьший шаг по времени там.

И что - можно множитель шага поднимать?

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

Занулить бы там модификатором скорость и считать дальше спокойно.

[a_schelyaev 366]

Если не срабатывает модификатор, то это повод для разбора.

В этом месте решалка должна прокашляться и поскакать дальше.

[Bonusfrag 104]

Если не срабатывает модификатор, то это повод для разбора.

В этом месте решалка должна прокашляться и поскакать дальше.

модификатор на таком месте вроде бы сработал 1 раз из трех.

И такие точки - они мигрируют и/или появляются в новых местах....(см.картинку) (скорость там упирается в ограничитель).

Правда, это дело происходило на 1м порядке точности. я этот момент не отследил в проекте, который скачал.

Сейчас поставил 2й порядок и 128 бит - лаги эти сдуло. Правда, они очень скоро и в более жестком развальном виде проявились в других местах ))

жесть )