Двухфазка во флюент

[cadcamcae 1]

Задачка. Дана конструкция труба в трубе диаметрами, например, 20 и 10мм. По внутренней трубе течет воздух, по внешней вода. Скорости обеих фаз 10 м/с. Давление во внутренней трубе 5 атм, во внешней 2 атм. Внутренняя труба имеет дефект, - сквозное отверстие диаметром 1мм. Необходимо рассчитать поле давления во внутренней трубе в окрестности дефекта.

Мое решение. 1) строю сетку одним доменом 2) задаю метод - vof, primary phase = вода, secondary= воздух 3) задаю граничные условия: следует отметить, что задаю на входе для смеси- скорость; для входа воздуха во внутреннюю трубу- volume friction=1. 4) далее устанавливаю инициализацию и метод решения 5) получаю решение - весь расченый домен заполнен primary phase, то есть водой, что на самом деле не так.

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

Спасибо.

[gebbels 1]

Задачка. Дана конструкция труба в трубе диаметрами, например, 20 и 10мм. По внутренней трубе течет воздух, по внешней вода. Скорости обеих фаз 10 м/с. Давление во внутренней трубе 5 атм, во внешней 2 атм. Внутренняя труба имеет дефект, - сквозное отверстие диаметром 1мм. Необходимо рассчитать поле давления во внутренней трубе в окрестности дефекта.

Мое решение. 1) строю сетку одним доменом 2) задаю метод - vof, primary phase = вода, secondary= воздух 3) задаю граничные условия: следует отметить, что задаю на входе для смеси- скорость; для входа воздуха во внутреннюю трубу- volume friction=1. 4) далее устанавливаю инициализацию и метод решения 5) получаю решение - весь расченый домен заполнен primary phase, то есть водой, что на самом деле не так.

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

Спасибо.

А на входе во внешнюю трубу задаёте правильно объёмную долю входящей воды? Если правильно, то поясните как именно вы проводите инициализацию? Вы используете разметку сеточной зоны и затем функцию патч чтобы установить начальное положение воды?

[cadcamcae 1]

А на входе во внешнюю трубу задаёте правильно объёмную долю входящей воды? Если правильно, то поясните как именно вы проводите инициализацию? Вы используете разметку сеточной зоны и затем функцию патч чтобы установить начальное положение воды?

Да, на входе во внешнюю трубу объемная доля воды=1. Граница воздуха зада на через патч. Инициализация: вначале инициализирую весь объем вода, затем патч воздух. Пуск на расчет .

Собственно расчет выполнен, вчера пошаманил с бубном...

Новый вопрос: как задать границу раздела фаз не патчем; с испозованием патча граница как бы получаемся предопределенной, но где она на самом деле?

Новая задачка вместе с тем: вместо внутренней трубы,- полый сосуд с некоторым объемом воздуха при высоком избытычном давлении. Сосуд также имеет дефект в виде отверстия малого диаметра. Воздух выходит во внешнюю среду. Необходимо смоделировать падение давления в сосуде за время истечения и собственно это время получить.

Если непонятно описана новая задача, пишите, поясню. Как ее делать, хз.

[gebbels 1]

Да, на входе во внешнюю трубу объемная доля воды=1. Граница воздуха зада на через патч. Инициализация: вначале инициализирую весь объем вода, затем патч воздух. Пуск на расчет .

Собственно расчет выполнен, вчера пошаманил с бубном...

Новый вопрос: как задать границу раздела фаз не патчем; с испозованием патча граница как бы получаемся предопределенной, но где она на самом деле?

Новая задачка вместе с тем: вместо внутренней трубы,- полый сосуд с некоторым объемом воздуха при высоком избытычном давлении. Сосуд также имеет дефект в виде отверстия малого диаметра. Воздух выходит во внешнюю среду. Необходимо смоделировать падение давления в сосуде за время истечения и собственно это время получить.

Если непонятно описана новая задача, пишите, поясню. Как ее делать, хз.

Задача у вас описана нормально, и как её делать тоже понятно. Вы помечаете для инициализации внутренний "баллон" как сеточную зону Adapt --> Mark, а затем опять же с помощью патча задаёте начальное давление (в первый момент времени) в этой области.

<noindex>http://www.youtube.com/watch?v=zxBZM55eujE</noindex>

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

"Новый вопрос: как задать границу раздела фаз не патчем; с испозованием патча граница как бы получаемся предопределенной, но где она на самом деле?"

Так где она на самом деле вы как раз определите в процессе расчёта. А так вы задаёте именно начальное приближение (нулевую иетрацию).

[cadcamcae 1]

Задача у вас описана нормально, и как её делать тоже понятно. Вы помечаете для инициализации внутренний "баллон" как сеточную зону Adapt --> Mark, а затем опять же с помощью патча задаёте начальное давление (в первый момент времени) в этой области.

<noindex>http://www.youtube.com/watch?v=zxBZM55eujE</noindex>

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

"Новый вопрос: как задать границу раздела фаз не патчем; с испозованием патча граница как бы получаемся предопределенной, но где она на самом деле?"

Так где она на самом деле вы как раз определите в процессе расчёта. А так вы задаёте именно начальное приближение (нулевую иетрацию).

Да, в первом приближении все понятно. Если будут вопросы реализации, я задам доп. вопросы. Спасибо за помощь.

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

Если дополонително задать, что у меня сосуд омывается водой с некой постоянной скоростью w,- какие граничные условия для этого использовать: inlet water скорость по компонентам+ volume fraction 1; а какая величина на outlet для воды будет? Ну pressure outlet water =0.

Короче вопроса два: как задать граничное учловие для воды( если в реале будет бинарная смесь воды и воздуха, поступающего из баллона)? Корректро ли задавать давление для outlet воды равным нулю?

Глобальное давление =1 атмосфере.

[gebbels 1]

Да, в первом приближении все понятно. Если будут вопросы реализации, я задам доп. вопросы. Спасибо за помощь.

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

Если дополонително задать, что у меня сосуд омывается водой с некой постоянной скоростью w,- какие граничные условия для этого использовать: inlet water скорость по компонентам+ volume fraction 1; а какая величина на outlet для воды будет? Ну pressure outlet water =0.

Короче вопроса два: как задать граничное учловие для воды( если в реале будет бинарная смесь воды и воздуха, поступающего из баллона)? Корректро ли задавать давление для outlet воды равным нулю?

Глобальное давление =1 атмосфере.

Ну так вода ведь выходит куда-то. Если вы зададите ноль относительно опорного давления,то будет выход в затопленное пространство, а если там отсос или что-то ещё, то тогда это и учитывайте. А вот объёмная доля воды - это хороший вопрос. Можно либо приконуть по массовым расходам примерно и задать, а можено просто поместить выход на удалении и задать там что угодно, расчёт всё равно не развалится даже если там будет 1 для воды.

[cadcamcae 1]

Ну так вода ведь выходит куда-то. Если вы зададите ноль относительно опорного давления,то будет выход в затопленное пространство, а если там отсос или что-то ещё, то тогда это и учитывайте. А вот объёмная доля воды - это хороший вопрос. Можно либо приконуть по массовым расходам примерно и задать, а можено просто поместить выход на удалении и задать там что угодно, расчёт всё равно не развалится даже если там будет 1 для воды.

Да, пока все так и делаю. Спасибо еще раз.

[cadcamcae 1]

+ банальный вопрос, под рукой нет флюента,- поэтому и спрашиваю: можно ли задать для воздуха, выходящего из объема сосуда, диаметр пузырьков? Таким образом, постараться соединить VOF и Эйлерову модели.

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

[cadcamcae 1]

+ банальный вопрос, под рукой нет флюента,- поэтому и спрашиваю: можно ли задать для воздуха, выходящего из объема сосуда, диаметр пузырьков? Таким образом, постараться соединить VOF и Эйлерову модели.

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

Да, видимо нельзя для vof диаметр пузырьков заать; модель только для большых пузырей...

[cadcamcae 1]

Гуру флюента и сфд в целом,- подскажите как диаметр пузырьков прикрутить к vof модели?

[cadcamcae 1]

Гуру флюента и сфд в целом,- подскажите как диаметр пузырьков прикрутить к vof модели?

Вопрос все еще актуален...