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

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

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

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

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

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

 

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

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

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

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

 

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


UnPinned posts
18.02.2023 в 22:04, Jesse сказал:

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

 

Есть такая вязкость "схемная".  Она зависит от порядка точности схемы. 

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

Есть такая вязкость "схемная".  Она зависит от порядка точности схемы. 

Вы про вязкость как численные ошибки?

https://en.m.wikipedia.org/wiki/Numerical_diffusion

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

Да. Погрешность апроксимации 

да, скорей всего из-за этого и появилась завихренность, которой не должно быть в жидкости с заданной вязкостью 10^-40

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • Danil89nur
      Доброе время суток. Столкнулся с проблемой на новой работе. Много деталей у который в свойствах отсутствует раздел "Материал" и приходится каждый раз вручную заходить туда, создавать раздел и выбирать там ""SW-Material@Name.SLDPRT". Возможно кто-то знает, либо сможет помочь с написанием простенького макроса, при нажатии на который в свойствах детали будет создаваться необходимый раздел и в нем сразу написан выбранный материал (желательно, чтоб там было указано сразу наименование материала, а не формула для его получения ""SW-Material@Name.SLDPRT"). Пытался самостоятельно через ChatGPT, но так как в програмировании я 0, то ничего не получилось)
    • Bot
      Коллеги, Благодарю вас за интерес к проекту электронной энциклопедии в области инженерного программного обеспечения PLMpedia, который более 15 лет развивался группой компаний ЛЕДАС. Представленный ниже скриншот передаёт содержание и структуру базы данных энциклопедии, накопленной за эти годы: именно такую начальную страницу, можно было до сегодняшнего дня увидеть, придя по адресу https://plmpedia.ru/ В последние несколько лет ЛЕДАС особенно интенсивно и результативно вкладывает все ресурсы в развитие своего основного бизнеса: см. ниже «О компании ЛЕДАС». Сегодня я с некоторым сожалением констатирую, что у компании ЛЕДАС более нет возможности отвлекать ресурсы своих сотрудников на некоммерческий (фактически волонтерский) проект PLMpedia. Вся база данных ресурса может быть безвозмездно передана любому юридическому или физическому лицу с единственным условием: наш преемник должен будет публично объявить о состоявшейся передаче и взятии на себя [...] View the full article
    • gudstartup
      @NordCom  у вас наверное это так происходит читайте этот раздел itnc_530.pdf
    • fenics555
      интересный индус. А нового то чего?
    • gudstartup
      откудова качаете мне интересно?
    • ДОБРЯК
      Совсем забыл про собственные вектора       собственный тон частота           0.013                        1   Расчет собственных форм и частот   1         1   7.0711E-01   0.0000E+00   0.0000E+00   0.0000E+00   0.0000E+00   0.0000E+00   1         2   7.0711E-01   0.0000E+00   0.0000E+00   0.0000E+00   0.0000E+00   0.0000E+00       собственный тон частота       15910.75                         2   Расчет собственных форм и частот   1         1   7.0711E-01   0.0000E+00   0.0000E+00   0.0000E+00   0.0000E+00   0.0000E+00   1         2  -7.0711E-01   0.0000E+00   0.0000E+00   0.0000E+00   0.0000E+00   0.0000E+00        
    • gudstartup
      так я это и написал
    • Viktor2004
      Обычно если сразу по всем осям, тогда это вентилятор БП. А на старых БП вентилятор может быть и на плате управления
    • gudstartup
      @NordCom какого года чпу? посмотрим более древние руководства а в принципе надо setup.elf запустить или у вас с windows?
    • Deunet
      Спасибо Вам! Но на пикчер ссылки не рабочие, пришлите пожалуйста, если у вас есть.
×
×
  • Создать...