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

Lamb–Oseen vortex


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



3 minutes ago, Bonusfrag said:

В чем делал?

В Математике

 

3 minutes ago, Bonusfrag said:

И можно же было сразу t=0 дать...

Да... Но у меня уже формула для V была забита. И попугать тебя хотелось)

Я ее все равно как массив чисел в CFX забирал)

Первая колонка - радиус в мм, вторая - давление в Па.

 

pressure_init.zip

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

спасибо, но у меня нет цилиндрических координат для НУ.

так что мне нужно будет P(x,y) делать - формулу или таблицу...

 

UPD вроде нашел как задать....

 

@piden

 

UPD2 - могу сказать, что мне начальное распреддаление не помогло. буквально за 0,1сек максимум с 0,02 уплывает в 0,035

 

 

Ссылка на сообщение
Поделиться на других сайтах
3 hours ago, Bonusfrag said:

буквально за 0,1сек максимум с 0,02 уплывает в 0,035

 


Time, [s]   Max Vel Radius, [mm]
0.0             22.52
0.010           24.41
0.05            25.07
0.11            27.08
0.20            30.67
0.30            31.85
0.40            31.85
0.50            33.84
0.60            33.84
...
1.00            37.50
2.00            44.06
3.00            49.53
4.00            53.77
5.00            55.47

 

https://i.imgur.com/FbERH3c.png

 

 

Сетка пока та же.

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

вот наложил... Только не понимаю, почему у тебя не из 0.02 график идет, если R0=0,02

ну, точнее, понимаю - 1,25^0,5  у тебя учтен...

compare.png

Ссылка на сообщение
Поделиться на других сайтах
18 minutes ago, Bonusfrag said:

Только не понимаю, почему у тебя не из 0.02 график идет, если R0=0,02

Оно то в условиях R0=0.02, но если в графике для скорости (можно в начальный момент времени) найти максимум, то он окажется в R=0.0224181.

 

https://i.imgur.com/EexuHw5.png

29 minutes ago, Bonusfrag said:

ну, точнее, понимаю - 1,25^0,5  у тебя учтен...

А я вот этого не понял) О чем ты?

Ссылка на сообщение
Поделиться на других сайтах
15 минут назад, piden сказал:

Оно то в условиях R0=0.02, но если в графике для скорости (можно в начальный момент времени) найти максимум, то он окажется в R=0.0224181.

 

 

А я вот этого не понял) О чем ты?

в вики есть эмпирический коэффициент 

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

у тебя сетка не имеет никакой симметрии относительно центра, из-за чего вихрь "колбасит" немного. 

У меня сетка - квадратики, с вершиной в 0.0.0 - т.е. какая-никакая симметрия есть. вихрь не колбасит.

Но лучше чем теор.кривая с вязкостью x50 я сделать не могу. 

FVvortex.png

Ссылка на сообщение
Поделиться на других сайтах
1 hour ago, Bonusfrag said:

Но лучше чем теор.кривая с вязкостью x50 я сделать не могу.

У меня рецепт счастья оказался таким:

1. Сетка (топология такая же) - замельчил в центре. Размер элементов: 0,5мм. Не хотелось с симметричной долго мудрить. С таким размером элементов и топологией оно строит ее минут 5... долго ждать, чтобы разные варианты пробовать))

 

2. Max CFL <= 5

3. Инициализация шага - по результатам предыдущего (без предсказания/экстраполяции)

 

Не знаю, насколько повлиял последний пункт - я его решил включить на всякий случай. Можно было еще double precision задействовать, но с ней скорость счета вообще грустная.

 

Этому варианту соответствует оранжевая линия на графике. Предыдущий (синяя линия) был c размером элементов 2мм, Max CFL <= 10, инициализация шага - значение по-умолчанию (нужно в справке смотреть что оно там химичит).

 

Резы:

 

https://i.imgur.com/Q5df0lD.png https://i.imgur.com/I80I2hk.png

 

Сетка:

 

https://i.imgur.com/o4AyO1d.jpg https://i.imgur.com/m0nAnhM.jpg

2 hours ago, Bonusfrag said:

у тебя сетка не имеет никакой симметрии относительно центра, из-за чего вихрь "колбасит" немного. 

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

2 hours ago, Bonusfrag said:

Но лучше чем теор.кривая с вязкостью x50 я сделать не могу.

Я б посоветовал попробовать новый проект с нуля создать. Может, пропустил что, или порядок в какой-нибудь константе потерял.

Могу тебе свой файл .cfx скинуть, если у тебя есть, где запустить.

Ссылка на сообщение
Поделиться на других сайтах
17 minutes ago, piden said:

2. Max CFL <= 5

3. Инициализация шага - по результатам предыдущего (без предсказания/экстраполяции)

+ еще критерий сходимости RMS Residuals < 1e-6 (в предыдущем случае - дефолтное 1e-4). По факту они оказались вообще ~5e-7 даже при одном цикле на шаг.

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

1. Сетка (топология такая же) - замельчил в центре. Размер элементов: 0,5мм. Не хотелось с симметричной долго мудрить. С таким размером элементов и топологией оно строит ее минут 5... долго ждать, чтобы разные варианты пробовать))

да, измельчение сетки в этой задаче приводит к улучшению результата. Это я проверял. На картинке у меня в центре 0,6 мм.

Цитата

2. Max CFL <= 5

С этим - наоборот: у меня чем больше шаг, тем меньше отклонение от теории. Но если совсем большой шаг сделать - то расколбашивает вихрь

Цитата

3. Инициализация шага - по результатам предыдущего (без предсказания/экстраполяции)

Это мне не понятно - о чем речь?

Цитата

Можно было еще double precision задействовать, но с ней скорость счета вообще грустная.

у меня есть имитация 128 бит, но она мне не давала никакой разницы

Цитата

+ еще критерий сходимости RMS Residuals < 1e-6 (в предыдущем случае - дефолтное 1e-4). По факту они оказались вообще ~5e-7 даже при одном цикле на шаг.

вот это попробую поменять. у меня по умолчанию 1е-05

 
Ссылка на сообщение
Поделиться на других сайтах
12 minutes ago, Bonusfrag said:

Это мне не понятно - о чем речь?

Ну, начиная считать шаг N, солвер смотрит на скорости/давления/всякое такое и их производные по времени из предыдущего шага N-1 и решает, с каких значений скорости/давления итп начинать итерации на этом шаге. Если "брать с предыдущего" - то инициализирует шаг значениями, полученными на предыдущем шаге) Если "экстраполировать" - то, зная производные, экстраполирует.



oYAIJve.png

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

@Bonusfrag , расскажи пока, ради чего вообще задача считается? Как ты будешь системную вязкость оценивать? И для чего это нужно?

Ведь, если на то пошло - то уже сейчас ты ее можешь оценить.. приняв результаты, полученные в ФВ как данность) 

 

И как ты значение циркуляции в начальных условиях выбрал? Наугад?)

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

Ещё: ты указывал ГУ как симметрию везде. У меня на внешней границе стоит стенка с проскальзыванием. Ведь симметрия подразумевает равенство нулю производной в направлении нормали.

Ссылка на сообщение
Поделиться на других сайтах
33 минуты назад, piden сказал:

Ещё: ты указывал ГУ как симметрию везде. У меня на внешней границе стоит стенка с проскальзыванием. Ведь симметрия подразумевает равенство нулю производной в направлении нормали.

да, это правильный шаг. хотя влиять не должен.

Цитата

И как ты значение циркуляции в начальных условиях выбрал? Наугад?)

если знать скорость и радиус - то она вычисляется по формуле. в pdf было же...

в принципе, это параметр, задающий макс.скорость

 

ради чего: была у меня в лабе задачка, где  DBD-разряд формирует вихрь. Но в CFD этот вихрь затухает быстрее, чем в эксперименте. Встал вопрос - почему.

Теперь понятно почему. Также видно,(спасибо!) что, если повозиться, то в ANSYS можно получить правильное решение базовой задачи.

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

Так что надо либо еще возиться, либо делать моделирование в другом ПО под эксперимент.

1 час назад, Bonusfrag сказал:

вот это попробую поменять. у меня по умолчанию 1е-05

не помогло....

 

2 часа назад, piden сказал:

Я б посоветовал попробовать новый проект с нуля создать. Может, пропустил что, или порядок в какой-нибудь константе потерял.

Могу тебе свой файл .cfx скинуть, если у тебя есть, где запустить.

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

 

 

 

UPD

хотя смотрю я на твою сетку, и понимаю, что с такой сеткой в 3D не перейти.

Ссылка на сообщение
Поделиться на других сайтах
1 hour ago, Bonusfrag said:

хотя смотрю я на твою сетку, и понимаю, что с такой сеткой в 3D не перейти.

Да.. С таким размером элемента и курантом это почти DNS получается)

Но если у вас там супер-комп имеется или нормальный кластер хотя бы :g: 

 

Вообще, мне странно что начальное распределение давления никакого влияния на твой результат не возымело. Не могу понять, как ФВ умудряется "самосогласованное давление" на первом шаге установить. Если не сложно, можешь показать картинку давления на нулевом шаге и на первом (у меня так в 16 посте).

 

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

 

Сейчас попробую сделать "красивую симметричную" сетку с размером элемента как в первом случае и запустить со всеми примочками. Узнаем, насколько тут сетка решает)

 

 

Ссылка на сообщение
Поделиться на других сайтах
15 hours ago, piden said:

Сейчас попробую сделать "красивую симметричную" сетку с размером элемента как в первом случае и запустить со всеми примочками. Узнаем, насколько тут сетка решает)

Попробовал. Сетка решает :biggrin:

 

 

https://i.imgur.com/XxfbbAB.png https://i.imgur.com/vr2cYVQ.png

 

 

"Почти структурированная" сетка. Размер элемента в центре - 2 мм.

https://i.imgur.com/JK1Lbks.png https://i.imgur.com/y9gSMM4.png

Ссылка на сообщение
Поделиться на других сайтах
20 часов назад, piden сказал:

Вообще, мне странно что начальное распределение давления никакого влияния на твой результат не возымело. Не могу понять, как ФВ умудряется "самосогласованное давление" на первом шаге установить. Если не сложно, можешь показать картинку давления на нулевом шаге и на первом (у меня так в 16 посте).

Нужно буквально 3-4 шага: давление после второго = -400, после 3го = -700, после 4го (время 0,2мс) :

 

press_vortex.png

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

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

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

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

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

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

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

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

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

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

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



  • Сообщения

    • Snake 60
      Скорее всего исправляли в конкретной сборке, для конкретных тел. Насколько я помню, надо в самом эскизе профиля заполнить свойство Description
    • Snake 60
      Делаете дырку в свае - вставляете в нее трубу и 2 человеко-силы крутят сваю ) Дёшево и сердито и не нужен никакой планетарный редуктор ))
    • gudstartup
      вы всегда ищите свой путь или будете делать то что вам советуют? с самого начала прежде чем приступать к экспериментам вам нужно было вынуть носитель и сделать образ. diskcopy делает копию гибкого диска на другой а вам надо сделать образ жесткого диска или карты памяти. делайте как советует @Viktor2004 иначе еще чего нибудь сломаете или вам станок не нужен и вы его приобрели для экспериментов ? вот а печатной машинке из прошлого нельзя а вот на синумерике можно - удивительные вещи в настоящем происходят!! на 100 % уверен что там не модуль а либо писимиси ай либо жд @mircomax81 хотите работать в дос почитайте помощь по командам 
    • PuKoLLleT
      Доброго дня всем. Поделитесь опытом,раньше такое не доводилось делать. На торце фланца нужно сделать  восемь отверстий ф17H7 на токарно-фрезерном станке с осью Y ,имеется ManualGuide.В моем понимании нужно просверлить,а затем торцевой фрезой расточить нужный диаметр.Получится такой фокус,или надо разверткой доводить нужный размер?  Стойка Fanuc 0i-Tf.
    • Snake 60
      Всё верно, я даже русификатор делал для него (ушло на это где-то около года). Хороший софт был, я всех знакомых на него подсадил в свое время :) У меня даже где-то образ на виртуалке лежит с SW2014+MechSoft, Жаль, что автодеск его поглотил, а не DS :( А по сути вопроса, проектировал двухступенчатый планетарный редуктор в свое время. Что было нужно: Справочник по планетаркам (автора не помню, если надо поищу), из софта MathCAD + GearTeq (софтина идет вместе с GearTrax, только рассчитывает и моделирует несколько зацеплений согласно выбранной схеме, в том числе планетарки) Ой и ошибаетесь. Звездочки - не равно зубчатые колеса в планетарке) Мы намучились с малой точностью изготовления на эл.эрозии, а Вы лазером ))) Был собран опытный образец редуктора с незакаленными деталями, чисто проверить геометрию, отладить сборку. Так потом начальство этот редуктор сказало испытать в бою, как я не противился)) Помер за пару дней интенсивных нагрузок ))) Материал, термообработка и точности изготовления решают )
    • Snake 60
      @Тихоход  Ещё один вариант - изучать программирование и писать свою программу/макрос. Вот пример конфигуратора двери: https://www.youtube.com/watch?v=wv4HryWQBSk
    • mircomax81
      Посредством Diskcopy можно сделать адекватную копию на флоппи диск?
    • Andrey_kzn
      Да, ведь советовали автору вытащить жёсткий диск и сделать образ. Неужели самому не интересно, как устроен этот старый промкомпьютер? Там может стоять как жёсткий диск, так и CF-карта, или же Disc-on-chip или Disc-on-module.  Как-то очень давно, на подобном промышленном оборудовании (большая печатная машина) я не смог перенести  ПО с родной CF - карты на другую, такого же размера но другой фирмы. Была также станочная оболочка под досом.  Образ я делал Нортон гостом, развернул его на другую CF-карту без проблем. Машина загрузилась нормально и работала ровно 2 дня, на неродной CF-карте,  потом зависла, и отказалась загружаться. Пришлось вызывать сервисника. От него мы узнали, что нельзя просто так взять и заменить карту  - носитель системного ПО, на другую. А полезли мы туда, так как руководство захотело печатать отчёты по работе машины. В станочном ПО такая функция была. Вот только принтер надо было купить у фирмы - производителя этой самой машины. В попытках заставить работать другой принтер мы редактировали файл  AUTOEXEС.BAT, на предмет загрузки драйвера принтера, ничего не получилось, зато перестали открываться некоторые страницы станочного ПО.  Образ у нас был давно снят, и мы развернули его на другую карту памяти. Родную карту я трогать не стал (как чувствовал, что нельзя ). Сервисник восстановил ПО с архивной дискеты на родную карту памяти, и всё заработало. Дискеты с архивами ПО шли вместе с документацией к машине.  Там были дискеты с Досом, станочной оболочкой, параметрами приводов, программой ПЛК. Он же и предупредил нас о предельно аккуратном обращении с родной картой памяти, и посоветовал немедленно сделать образы с архивных дискет, что мы сразу и проделали.. Вот такой сюрприз от машины эпохи Дос.
    • Борман
    • gudstartup
      а ведь все работало...
×
×
  • Создать...