Jump to content

Recommended Posts

Jesse

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

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

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

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

 

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

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

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

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

 

Link to post
Share on other sites


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

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

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

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

вот эпюра

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

Link to post
Share on other sites
Jesse
4 минуты назад, soklakov сказал:

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

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

 

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

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

неа

 

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

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

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

 

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

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

Link to post
Share on other sites
soklakov
9 часов назад, Jesse сказал:

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

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

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

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

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

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

Link to post
Share on other sites
Jesse
12 часов назад, soklakov сказал:

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

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

 

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

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

 

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

 

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

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

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

Link to post
Share on other sites
soklakov
17 минут назад, Jesse сказал:

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

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

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

 

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

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

Link to post
Share on other sites
  • 2 weeks later...
Foksmen
18.02.2023 в 23:04, Jesse сказал:

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

 

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

 

1_.png

 

2_.png

 

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

 

setka_.png

 

setka2_.png

 

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

Edited by Foksmen
  • Нравится 1
Link to post
Share on other sites
Jesse
Posted (edited)
17 минут назад, Foksmen сказал:

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

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

 

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

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

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

Edited by Jesse
Link to post
Share on other sites
Foksmen

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

Link to post
Share on other sites
Foksmen
4 часа назад, Jesse сказал:

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

 

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

Link to post
Share on other sites
Jesse
4 часа назад, Foksmen сказал:

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

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

 

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

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

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

Link to post
Share on other sites
Foksmen
17 минут назад, Jesse сказал:

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

 

Спасибо )

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

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

 

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

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

Edited by Foksmen
Link to post
Share on other sites
soklakov
38 минут назад, Foksmen сказал:

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

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

 

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

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

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

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

Link to post
Share on other sites
Foksmen
26 минут назад, soklakov сказал:

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

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

  • Нравится 1
Link to post
Share on other sites
soklakov
2 минуты назад, Foksmen сказал:

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

Вряд ли. 

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

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

  • Нравится 1
Link to post
Share on other sites
Foksmen

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

Link to post
Share on other sites
Jesse

 

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

Link to post
Share on other sites
Foksmen

Ну... я так на всякий случай ) Вдруг сомнения какие )

 

Link to post
Share on other sites
  • 3 weeks later...
Солнечный энергетик

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

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

Link to post
Share on other sites
Foksmen

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

Edited by Foksmen
Link to post
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.
Note: Your post will require moderator approval before it will be visible.

Guest
Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

  • Recently Browsing   0 members

    No registered users viewing this page.




  • Сообщения

    • миталник
      Шпонка, передавая крутящий момент, работает на срез и смятие. Исходя из этого и подбирается (рассчитывается) размер шпонки.
    • AlexKaz
      Под результатами ессно надо понимать экономию времени и денег. Если покупка более дорогой и худшей по качеству услуги выливается в тотальный геморрой и не приводит ни к какой экономии ресурсов - нафиг не нужны такие услуги...
    • AlexKaz
      IMHO, не по Сеньке шапка. Софт англосаксов - это следствие диверсификации их экономик. Софт - это сфера услуг, призванная экономить деньги. А вот деньги зарабатывать можно только на товарном производстве (рынок акций, депозитов, и прочие финансы - это только обслуга реального сектора). Толкая софт впереди паровоза экономики как её драйвер есть риск не добиться ни-че-го, кроме довольных рож софтописателей. Во-первых, эти редиски ломят цену за софт внутри РФ словно уже живут за бугром и получают з.п. как бугром. Во-вторых, рынок РФ они не ставят приоритетным, РФ для них - это 1/50 мирового пирога - так себе полигон для отработки технологий. Соответственно, эти люди не понимают в массе что нужно миру - потому как в РФ экономика откровенно не выдающаяся и всех тонкостей и нужностей силами спецов РФ не создать и не описать. Как следствие, на выходе имеем большой пшик - деньги влиты, а результаты мир не устраивают.   Лунная же миссия и будущий Марс - это резервирование человеческой популяции, и чем раньше сделать бэкап - тем лучше. В глобальном смысле. В локальном же смысле эти проекты есть уйма заказов на железо. А железо можно и без софта создать или использовать примитивные софтины... Чёт Королёв без CAD/CAE обходился (но не без вычислитетелей). Вот ведь странный человек - ни одной CAD-шараге не забашлял...
    • Александр1979
      Такое бывает если фазы на синхронный двигатель перепутаны.
    • lublu_abs
      понял посмотрите для одного диаметра вала b и h   Посмотрите на b и h для одного диаметра вала, например 20 и 45. гост на конические концы валов Гост шпоночные соединения: Как не сложно заметить, размеры отличаются  
    • irbitsky
      Могут быть любыми. Имеется ввиду на условном виде сверху. Это касается узких торцов паза. А вот  Они идентичны    Или я просто не заметил разницы. Укажите конкретно, какие размеры отличаются. Обоими. Могут быть любыми. Имеется ввиду на условном виде сверху. Это касается узких торцов паза. А вот  Они идентичны    Или я просто не заметил разницы. Укажите конкретно, какие размеры отличаются. Обоими.
    • Viktor2004
      Ну, рассказываю все по порядку Купили манагеры глобусный стол но без двигателей. Так дешевле Потом где-то по объявлению купили два движка и привод. И сказали собирать. Тут уже в качестве юмора постили картинку где вместо шин постоянного тока кусок кабеля ВВГ примерно 25 кв мм. Так вот я вместо шин скрутил в несколько слоев кабель ВВГ 2,5 кв мм. Оптические пачкорды у меня были, а вот соединители CXA2A --- CXA2B пришлось поискать. Вроде в ЧИП и ДИПе они на сайте есть, но в реалии не оказалось. С трудом достали. Там же достали разъемы кабелей 4-5 осей. Осталось сварганить кабели. Купили контрольный электрический кабель. Спаяли, подключили, а он не работает. SV0368 Заземление эффекта не дало. Поехал искать витую пару, да такую, что бы каждая жила была многожильной. Нашел на рынке. Витую пару подключил только на линию данных. 5-6 контакт энколера. Остальные провода оставил от контрольного кабеля. У него сечение побольше. Ошибка SV0368 ушла. Затем оказалось что цоколевка разъемов энкодера 4 и 5 оси разная. Конечно на двигателях и на приводе она одинаковая, а на промежуточном разъеме оси она разная. Пришлось перепаивать. Подключил. Тут оказалось что энкодер одного двигателя не обнуляется и не отслеживает положение вала. Менял местами двигатели, отдельно энкодеры, не работает именно тот энкодер. Стал пытаться запускать все по очереди. Забил параметры как указал выше и стал проверять сначала на 4-й оси, потом на 5-й. Поочередно отключая то одну, то другую. При движении от маховичка вначале ось страгивается с места, потом начинает гудеть, а двигатель греться. Иногда дергается так что аж подпрыгивает. На обеих каналах ведет себя одинаково что один, что другой движок, когда исправный энкодер на них вешаю.    
    • lem_on
      че за бред, сам сначала влез, а потом я ему пишу оказывается.
    • статист
      Никогда этой версии не ставил. Попробуйте 21 версию поставить. С ней таких проблем не наблюдалось
    • Maik812
      Я переименования файла с разными свойствами решал так.  Все го то сделать нужно плагин который писать будет через апи свойство в фаил. На прямую он пишет но ссылка то на имя старое. Создать библиотеку с именами в файлах со свойствами различными в любом порядке и тд за пару кликов.. На данном фото переименовании по габаритам и по имени папки где лежат файлы. Все сохраняется по шаблону и просто выбирается.
×
×
  • Create New...