Много глупых вопросов по FlowVision

[Avis 0]

Всем салют.

Ситуация такая: решил перейти на FlowVision 2.3 после того как разочаровался в адекватности решений во FloWorks.

Посему есть много глупых вопросов. Попрошу ногами не бить))) сразу)))

И так по порядку:

- как средствами FlowVision определить аэродинамические параметры крыла конечного размаха? Продуть то я его продул, а вот что дальше?

V=10 м/с

Среда - воздух из базы

Модель - несжимаемая жидкость.

На сайте в примерах я наблюдал поляры профилей полученные в FlowVision

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

Подскажите. Может я где то сильно ошибаюсь при задании параметров.

Хочу получить результат обтекания тела со скоростью 5 Мах

Параметры следующие:

V=331,46*5=1657,3 м/с

T=273,15

P=101 325 Па

Модель - полностью сжимаемая жидкость

Среда - воздух из базы

Заранее всем спасибо.

[a_schelyaev 366]

А дальше нужно взять слой Характеристика по Давлению, построенного на объекте Супергруппа, созданного на ГУ Стенка на твоем крыле.

Старая демка имела в себе баг, в ГУ Вход/Выход - что-то там с давлением на входе не так определялось.

Но можно на входе тоже Супергрупу сделать, а по ней Заливку или Изолинии и оценить велчину скорости на входе.

Или график вдоль кривой (по входу).

[Avis 0]

А дальше нужно взять слой Характеристика по Давлению, построенного на объекте Супергруппа, созданного на ГУ Стенка на твоем крыле.

Старая демка имела в себе баг, в ГУ Вход/Выход - что-то там с давлением на входе не так определялось.

Но можно на входе тоже Супергрупу сделать, а по ней Заливку или Изолинии и оценить велчину скорости на входе.

Или график вдоль кривой (по входу).

Не совсем понял про крыло. Т.е мы получаем график распределения давления по поверхности? А как дельше его интегрировать чтоб получить интегральную характеристику с тем чтоб получить Х и У ну или хотя бы R и угол наклона вектора R?

У меня интересные результаты получаются, если скорость на входе от 1657 м/с до примерно 2300 м/с число Маха на входе не превышеет значения 3.5 Мах. А если скорость на входе больше 2300 м/с число Маха сразу прыгает до 10 Мах. Где я мог напортачить? )))))

[a_schelyaev 366]

Практически везде.

График по крылу я строить не предлагал. Перечитай внимательней еще раз.

Также открой Учебник и пробегись по задачкам. Почитай FAQ в документации.

[sheyko.e 0]

a_schelyaev, помогите. пожалуйста разобраться с FlowVision...

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

Нужно узнать давление на выходе

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

Какое граничное условие на выходе надо задавать?

[a_schelyaev 366]

Открываем Учебник и делаем примеры.

Осмысливаем проделанное и приступаем к своей задаче.

[sheyko.e 0]

a_schelyaev, и учебник открывал (встроенное пособие и методичка МФТИ), и примеры делал...видимо проблема в осмысливании проделанного...

переосмыслил..все получилось. спасибо за внимание)

Изменено 20 октября 2009 пользователем sheyko.e

[nibulla 0]

а почету так сразу на FlowVision? Туда даже геометрию нормально не передашь, и модель турбулентности там тоже одна свего. Уж если перезодить так сразу на CFX или на Fluent. Вопрос где взять? так в пиринговой сети все есть - регестрируйся и качай! Успехов! P.S. Еще один совет занялся CFD - задудь про форточки - Linux, вот тема,

[a_schelyaev 366]

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

Под виндой все работает на ура, на линуксе тоже.

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

Вывод - знания периодически нужно освежать.

[ChancE 0]

В общем передо мной стоит задача сверхзвукового течения внутри допустим сопла (у меня более сложная геометрия).

Модель: полностью сжимаемая жидкость

На входе в сопло указываю (вход/выход):

Температура: значение на стенке (3056 К)

Скорость: нормальная скорость с давлением (1014; 5209000 Па)

ТурбЭнергия: пульсация (0,02)

ТурбДиссипация: масштаб турбулентности (0,0068м 1/10d входного)

На выходе из сопла указываю (свободный выход):

Температура:нулевой поток

Скорость: нулевое давление/выход (избыточное давления на выходе 0)

ТурбЭнергия: нулевой поток

ТурбДиссипация: нулевой поток

На стенки сопла (стенка):

Температура:нулевой поток

Скорость: стенка, логарифмический закон (шероховатость = 0)

ТурбЭнергия: значение в ячейке рядом со стенкой

ТурбДиссипация: значение в ячейке рядом со стенкой

Сетка: 100/100/100

Общие параметры:

Время конечное: 3 сек

Неявная схема; скошенная схема

Явное число куранта: 0,5

CFL 1 макс 1

Настройка: стационарная; убирать маленькие ячейки.

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

В чем моя ошибка, помогите разобраться? (может быть коряво все описал, если что-то непонятно, то задавайте вопросы).

[a_schelyaev 366]

Скорость: нормальная скорость с давлением (1014; 5209000 Па)

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

И температура тоже задается относительно опорной величины.

На выходе из сопла указываю (свободный выход):

Если на выходе св/зв истечение, то в свойствах данного ГУ для скорости укажи Свободную границу.

Неявная схема; скошенная схема

Явное число куранта: 0,5

Тут пока ничего не трогай. Оставь как было по умолчанию.

CFL 1 макс 1

На первых итерациях шаг интегрирования может скакнуть.

Поэтому Курант оставь 1, а сверху шаг ограничь 1e-5 для начала.

Потом ограничение ослабь.

Настройка: стационарная; убирать маленькие ячейки.

Если все вышележащее не поможет, то попробуй Совместную схему (2 итерации).

И в параметрах метода не забудь заптичить Использовать BsSolver для всех уравнений - расчет веселее пойдет.

Ну и желательно работать на последней версии.

P.S. Как крайняя мера - поставь в Параметрах метода 1й порядок аппроксимации для всех уравннений. Дойди до стационара и переключись на высокий порядок как и было.

[ChancE 0]

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

[a_schelyaev 366]

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

Хотя в новой версии и сверхзвук с курантом 100 считали, чтобы металл побыстрее прогрелся.

[ChancE 0]

Новая, это 3.06.00? Я сам работаю в 2.5.04 версии, хотя вроде видел где то примерчики и с большими курантом, к сожалению забыл где, недавно искал, но не нашел. Сам пробовал играть с числом Куранта, при увеличении числа невязки улучшались (правда значения искомых параметров мне не нравились все равно), но при достаточно больших числах решение вышибало с ошибкой (давно это было и диапозоны которые я брал не помню), может если время будет попробую снова. Проблема в моей задачи в том, что стационарность достигается очень не скоро, а при увеличении шага по времени решение вышибает, по ходу это связано с начальными параметрами, которые я задаю в области, думаю сделать одномерный расчет в другой программе этой же задачи и забить результат в начальные значения, тогда наверно стационарность наступит быстрее. Кстати, во FlowVision возможно решить задачу допустим движения шарика в трубе с жидкостью? Примеры с подвижными сетками (ради них собственно и перешел на FlowVision) только начал разбирать, вот возник такой вопрос. И ещё вопрос: почему с течением времени у меня меняется давление на входе, как это можно исправить,как сделать чтобы оно было постоянным?

Изменено 30 ноября 2009 пользователем ChancE

[fedigen 0]

Ребят, кто-нибудь считал во FlowVision c шагом по времени ~10^-9 с? Сможет ли он разрешить процесс с таким мелким шагом?

[a_schelyaev 366]

Новая, это 3.06.00? Я сам работаю в 2.5.04 версии, хотя вроде видел где то примерчики и с большими курантом, к сожалению забыл где, недавно искал, но не нашел. Сам пробовал играть с числом Куранта, при увеличении числа невязки улучшались (правда значения искомых параметров мне не нравились все равно), но при достаточно больших числах решение вышибало с ошибкой (давно это было и диапозоны которые я брал не помню), может если время будет попробую снова. Проблема в моей задачи в том, что стационарность достигается очень не скоро, а при увеличении шага по времени решение вышибает, по ходу это связано с начальными параметрами, которые я задаю в области, думаю сделать одномерный расчет в другой программе этой же задачи и забить результат в начальные значения, тогда наверно стационарность наступит быстрее. Кстати, во FlowVision возможно решить задачу допустим движения шарика в трубе с жидкостью? Примеры с подвижными сетками (ради них собственно и перешел на FlowVision) только начал разбирать, вот возник такой вопрос. И ещё вопрос: почему с течением времени у меня меняется давление на входе, как это можно исправить,как сделать чтобы оно было постоянным?

Новая это 3.07.

Шарик в трубе - искать в документации фильтр подвижного тела. Но мы не моделируем контакт со стенкой. Контакт рока только через ABAQUS. В 3.07.01 сделали ограничители по перемещениям и этот вопрос частично закрыли.

Давление меняется потому что течение устаканивается. Если нужно постоянное давление, то его и нужно задавать на ГУ. Если давление снимается с плоскости, то там значения берутся из ячеек, а не с ГУ.

Ребят, кто-нибудь считал во FlowVision c шагом по времени ~10^-9 с? Сможет ли он разрешить процесс с таким мелким шагом?

Ну в минус девятой что-то не припомню, но вот в минус восьмой точно було.

Профили_ударных_волн__рассчитанные_с_различными_пространственными.pdf

[fedigen 0]

Спасибо огромное! Буду пробовать считать с мелким шагом!

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

Постараюсь вкратце описать задачу...

На плоскости равномерно распределены "пятнышки" диаметром ~10мкм. В определенный момент времени они все начинают греться импульсно (~10 нс) до температуры ~1500 K. Необходимо построить распределение давления и температуры в плоскости расположения этих пятнышек.

На этом форуме уже поднимали вопрос о создании и разметке геометрии для FlowVision. Я сделал это следующим образом:

Я сделал цилиндрические выступы самой меньшей какую мог обеспечить солид высоты:

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

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

Может как-то разместить в объеме эти области? Я таким образом хочу посмотреть расширение маленьких объемчиков газа при его нагревании..., в 0-приближении на плоскости. Может есть другой способ локального нагрева вещества внутри счетного объема?

Вот такой глупый у меня вопрос :) Вообще, если откровенно, понравился мне FlowVision. Хорошая вещь!

[a_schelyaev 366]

В данной геометрической постановке альтернатив нет.

Можно еще инструментами CAD-системы создать на плоскости выделенные зоны, раскрасить их другим цветом и сохранить в VRML, т.к. он цвета понимает и они в разные ГУ в FV превратятся.

Или же просто внести через фильтры объемный источник тепла.

[ChancE 0]

Спасибо.

Но мы не моделируем контакт со стенкой

То есть во FlowVision нельзя смоделировать движение какого-либо тела в ограниченном объеме жикости?

Если давление снимается с плоскости, то там значения берутся из ячеек, а не с ГУ.

Да я понимаю. Но если я допустим на входе задаю давление, температуру, и скорость, а на выходе свободную границу, то почему он изменяет давление на входе (по сечению он меняет давление, это понятно, но вход у меня описывает только одна ячейка, откуда он там берет 20-30 значений непонятно, причем эти давления далеки от заданных мной), хотя скорость и температура во всей плоскости одни, и значения те которые я задал.

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

Вопрос в том, можно ли задать скорость и давление не изменными во входном сечении (желательно по всему сечению)? Для снятия данных я пользуюсь экспортом входной плоскости.

[a_schelyaev 366]

1. Движение тела смоделировать можно, но тело будет проходить через стенки. Если нужно моделировать контакт, то нужно использовать ABAQUS, с которым у FV реализована прямая связь для таких случаев. Или же дождаться очередного обновление для FV HPC версии 3.07.01, где будет реализован ограничитель на такие случаи.

2. В кишках FV алгоритм оперирует с понятием МАССОВЫЙ РАСХОД. Все что мы видим в ГУ используется лишь для того, чтобы определить расход на входе в кооперациями со опорными величинами и свойствами вещества. Поэтому, однозначен только массовый и объемный (для несжимаемой модели) расход.

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

:)