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

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

Продолжаю тему сингулярности. На этот раз цель получить сингулярность, т.е. отсутствие сеточной сходимости около тупого угла в потоке идеальной жидкости. 
Ожидал получить что-то типа такого
синг.png

Считал в floefd в 2D. Задал для воды вязкость 10e-40. На всякий включил гравитацию.
Сверху симметрия. Самая плотная сетка (решил не борщить, для ускорения процесса - много расчетиков прогонял)

mesh1.jpgmesh2.jpg
Пробовал обычный нестационарный k-e решатель в первом случае и ламинарный решатель (с вырубленной турбулентной вязкостью) во втором случае.

На эпюре скорости скачков вообще не видать и на разных сетках результаты похожи 
скорость.jpgскорость увелич.jpg

 

а вот эпюра завихренности оказалась многообещающей: тут для разных сеток с кол-вом ячеек в ~2 раза больше в каждом расчёте результаты получились ~0, 57 и 220!

ЗАВИХРЕННОСТЬ САМАЯ ГРУБАЯ СЕТКА.jpgЗАВИХРЕННОСТЬ СРЕДНЯЯ СЕТКА.jpgЗАВИХРЕННОСТЬ ПЛОТНАЯ СЕТКА.jpg

Вообще результаты довольно странные... можно начать с того, что завихренность бывает только в вязкой жидкости и отсутствует в идеальной (если её там нет в "начальных условиях")
Хотя возможно физичные результаты получились как раз на грубой сетке (там где завихренность == 0), а на плотной сетке начала проявлять себя численные особенности всякие.
Вообще, ожидал получить сингулярность по полю скорости:biggrin:
Мб такие задачи надо только в DNS решать?

Вообще кто-то что-то подобное делал?)

 

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


UnPinned posts
15 часов назад, Jesse сказал:

Вообще, ожидал получить сингулярность по полю скорости

Это все равно что получить сингулярность по перемещениям. Получается такое?

15 часов назад, Jesse сказал:

вот эпюра

А размеры вихрей не меняются с сеткой?

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

Это все равно что получить сингулярность по перемещениям.

а где это писано?) я импровизировал:biggrin:

 

4 минуты назад, soklakov сказал:

Получается такое?

неа

 

4 минуты назад, soklakov сказал:

А размеры вихрей не меняются с сеткой?

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

 

вот тут лабу делал человек (первая попавшаяся ссылка)
https://vunivere.ru/work11691
ориентировался на это. Плюс минус те же габариты и всё остальное. Правда у них стационарный ламинарный вихрь. А у меня они двигаются дальше уступа.

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

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

размеры вихрей не меняется.

Это нестационар?

9 часов назад, Jesse сказал:

а где это писано?

 Не знаю) но в мкэ мы ищем перемещения. А в мко скорости

А вот завихренность может и градиент какой-то))

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

Это нестационар?

на картинках да. А так пробовал и так  , и сяк..

 

12 часов назад, soklakov сказал:

в мкэ мы ищем перемещения. А в мко скорости

 

всё верно. Но кто сказал что в МКО мы не можем получить сингулярность по скорости? почему?)

 

12 часов назад, soklakov сказал:

А вот завихренность может и градиент какой-то

да, аналогия работает..

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

всё верно. Но кто сказал что в МКО мы не можем получить сингулярность по скорости? почему?)

Чтобы продать что-нибудь ненужное...

Чтобы получить сингулярность, надо что-нибудь на ноль поделить. Ну не ноль, на бесконечно малое. А поле скоростей, как и перемещений - мы находим в виде кусочно линейной (в случае линейных элементов конечно) функции. Там негде поделить на ноль и нет никаких бесконечных значений.

 

Вот градиент - то, что нужно, чтобы получить сингулярность при бесконечно малом элементе.

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

Ссылка на сообщение
Поделиться на других сайтах
  • 2 недели спустя...
18.02.2023 в 23:04, Jesse сказал:

Вообще кто-то что-то подобное делал?)

 

Пробовал считать течение вязкого газа в канале с квадратными призмами. Атмосферные условия, скорость на входе 50 м/с. Размеры области 0.2х0.11 м. Число Рейнольдса Re=25000. Использовал модель Смагоринского. Получились такие картинки для разных моментов времени.

 

1_.png

 

2_.png

 

Сетка 530000 элементов.

 

setka_.png

 

setka2_.png

 

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

Изменено пользователем Foksmen
Ссылка на сообщение
Поделиться на других сайтах
17 минут назад, Foksmen сказал:

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

что за программа? и картинки чего? Поля Завихренности? И почему всё чёрно-белое?..)

 

17 минут назад, Foksmen сказал:

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

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

Изменено пользователем Jesse
Ссылка на сообщение
Поделиться на других сайтах
4 часа назад, Jesse сказал:

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

 

Параметрические исследования можно провести конечно. Пока было бы интересно узнать можно ли получить подобные картинки с использованием k-e или других параметрических моделей в принципе?

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

Нет это не завихрённость, это линии тока. Программа своя ) 

понял. За свою программу респект!)

 

29 минут назад, Foksmen сказал:

можно ли получить подобные картинки с использованием k-e или других параметрических моделей в принципе?

так у меня и k-e. По идее это не должно влиять. А у вас прямой числ. счёт?

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

понял. За свою программу респект!

 

Спасибо )

17 минут назад, Jesse сказал:

так у меня и k-e. По идее это не должно влиять. А у вас прямой числ. счёт?

 

Прямой. Уравнения Навье-Стокса с моделью Смагоринского. Схема 2-го порядка точности по пространству.

Забыл написать. Адаптации сетки в процессе расчёта не производилось. Никаких возмущений на входе не задавалось.

Изменено пользователем Foksmen
Ссылка на сообщение
Поделиться на других сайтах
38 минут назад, Foksmen сказал:

Прямой. Уравнения Навье-Стокса с моделью Смагоринского.

Кажется это все таки не прямой, а LES. 

 

В прямом Смагоринского подставлять некуда.

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

Пока было бы интересно узнать можно ли получить подобные картинки с использованием k-e или других параметрических моделей в принципе?

Да. Только лучше всё-таки без осреднения по Рейнольдсу)

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

Кажется это все таки не прямой, а LES. 

Ну да LES. я так понял не заподозрил ли меня Jesse в мультипликации ) Есть такие программки. Нет, это численный расчёт. В реале. Без мультиков.

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

я так понял не заподозрил ли меня Jesse в мультипликации )

Вряд ли. 

А вы могли бы?? Знаете кого-то кто мог бы? "Подсовывать мультики вместо расчета".

Тут такого отродясь не видали.

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

Ну и хорошо, если нет. Здесь может и нет, но позднее дам ссылку. Кстати он не скрывает, что это мультик. Просто снимает мультики.

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

 

@Foksmen  мы все тут мультифизикой занимаемся, а не спец эффекты для голливудских блокбастеров в океане пилим)))

Ссылка на сообщение
Поделиться на других сайтах
  • 3 недели спустя...
Солнечный энергетик

В этих вопросах не голливуд рулит, а  Болливуд. Какой фильм не возьми, задумываться начинаешь, может на другой планете снято? 

Или в крутилках программ всеми физическими постоянными играются.

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

Вопрос философский. От задачи зависит. Может и играются. Получили правильные коэффициенты,  может и нормально на определённом этапе.

Изменено пользователем Foksmen
Ссылка на сообщение
Поделиться на других сайтах

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • Bot
      АСКОН, российский разработчик инженерного программного обеспечения и ИТ-интегратор, подвел итоги 2023 года. Выручка компании выросла на 47% и достигла 5,6 млрд рублей, штат сотрудников увеличился до 1250 человек. Клиентская база пополнилась тысячью предприятий и сегодня с решениями АСКОН работают почти 16 000 заказчиков из всех отраслей промышленности и строительства. Смотреть полностью
    • Sergey78
    • scrimline
      Уговаривать ему пришлось да и оскорбили его. Я вообще то сразу написал что постараюсь ладдер скинуть, а ты сарказмом со всех щелей хлестать начал. 
    • andrey2147
      Я подключал 532 , позже поищу ссылки, ВинПсином и кабелем там хитрость в настройках винпсина 7 бит надо ставить вам зачем ?
    • Guhl
      Надеюсь, в моих словах вы ничего оскорбительного не усмотрели Или это я?
    • Guhl
      Представляю. Я же не предлагал на полной скорости считать Получается, что до м19 считает правильно, а после - неправильно?
    • mannul
      Вы сначала проверьте точно ли режущий край пластины в оси детали. А то китайские державки и пластины бывают с такими косяками.
    • kkk
      Что-то вдруг резко понадобилась возможность посчитать площадь разной отделки в изделии. Есть куча многотельных деталей, какие-то поверхности красятся одним цветом, какие-то другим, какие-то третьим. Раскрасить в модели проблем нет, а нет ли возможности высчитать эту информацию автоматом? Есть ли какой-то способ наложения "отделки", чтобы Солид мог ей оперировать в дальнейшем.    Раньше, если изредка нужно было что-то подобное посчитать, просто выбирал нужные грани мышкой и смотрел сумму площадей в инструменте измерения. Но хочется больше интерактива, чтобы Солид сам отслеживал, что у него чем "покрашено".
    • andrey2147
      Пришлось уговаривать оператора что-бы выполнил просьбу Виктора, ещё и  оскорбили.  
    • gudstartup
      да не надо ничего паять система прекрасно считает импульсы на оборот в 445 параметре диагностики вот только беда в том что один раз войдя в режим контроля позиционирования после m19 она из него уже не выходит хотя сразу после включения находится в режиме управления скоростью и 445 параметр ничего не отображает и контроль числа импульсов отсутствует. Вопрос как выключить контроль позиционирования. он вводится по G70.6 но и после того как этот сигнал остается система продолжает контролировать кодер. в обще то понятно что это датчик так как на скорости 100 или 200 оборотов а также при ориентации все хорошо а потом видимо амплитуда сигнала падает ниже нормы там все ок вы частоту представляете при 5000 об 5000*4096
×
×
  • Создать...