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

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

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

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




  • Сообщения

    • Ninja
    • Ninja
      В зимнее время водители устанавливают на колёса специальные цепи против скольжения. В мире разработано множество видов цепей и способов их установки. Японцы применяют "цепи Абэ". Кто этот великий инженер? Это знаменитый японский писатель Кобо Абэ. Для уединенного писательского труда Абэ купил домик в горах. А в горах зима ранняя, снежная, затяжная. Без цепей на колёсах можно и улететь с обрыва. Абэ периодически заводил машину и спускался в долину. Почту получить/отправить, в магазине продукты купить, ящик рисовой водки. (Известное дело: насухую великую книгу не напишешь). В долине приходилось цепи снимать. По правилам цепи запрещается применять при отсутствии на дороге снежного покрова. В те времена для установки/снятия цепи требовалось домкратом поднимать каждое колесо. В машине 4 колеса. Товарищу Абэ приходилось корячиться с домкратом. При этом в самых нелитературных выражениях комментируя политику правящей партии и все домкраты в мире. Купил в магазине бухлишко-закусь. Поехал домой в горы. А там снег, скользко. Цепи требуются. Опять писатель и драматург Абэ домкрат крутит. Всё это насмерть надоело Абэ. Обратился к компаниям-производителям: придумайте способ установки цепей без домкрата! Но производители ничего не придумали. Или не захотели. Тогда Абэ сам придумал способ установки цепей без необходимости поднимать колесо. Послал изобретение на конкурс. И неожиданно получил призовое место и патент. Новый способ установки цепей быстро завоевал популярность. Так и пошло у японцев - "цепи Абэ". Кобо Абэ не получил Нобелевскую премию по литературе. Уже выдвигался, уже все были уверены в том, что он её получит. Но не успели, писатель умер. А это всё дурацкие цепи виноваты. Сколько времени он потерял в установке/снятии цепей? Сколько времени, вместо того чтобы писать, он крутил домкрат? Вот и не успел. https://ru.ruwiki.ru/wiki/Абэ,_Кобо
    • gudstartup
      приемлимо. температуру мотора посмотрите. если вентилятор нормально вращается и его скорость соответствует оригинальному а также контакт в разъеме нормальный а F все равно появляется то даже не знаю что у вас может быть так как привод у вас новый
    • Koels
      Хмм, я думал раз она вылазит даже в простое, то Z координата висит на тормозе и дело не в этом. Спасибо за мысль, щас узнаю. Вообще у нас рядом стоит точно такая же макина, точь в точь и там таким проблем нет. Много заказов и отсутствие специалистов сделало свое дело. :) @gudstartup, 58% нагрузка на координату
    • mrVladimir
      Что-то я немного засомневался. Если решим брать новое ЧПУ, то искать точно такое же необходимо по номеру сзади корпуса (пластикового пластмассового) - так? А если брать только плату (материнскую или как ее лучше назвать...), на которой установлена (запаена) микросхема 1 на моих фото выше, то искать плату необходимо по номеру, который указан на самой плате - так? В моем случае : ЧПУ : A02B-0321-B520. Код на плате : A20B-8201-0081/01A. Если, допустим, мы покупаем ЧПУ A02B-0321-B520 - будет ли это гарантией того, что в нем установлена плата A20B-8201-0081/01A. И можно ли будет его считать таким же. И с отдельно приобретаемой платой тоже самое - если номер на ней полностью совпадает с нашим, значит ли это, что она полностью идентична. P.S. : извиняюсь за , возможно, навязчивые вопросы. Просто не хотелось бы попасть впросак из-за своей некомпетентности. такой файл есть, но мне говорили, что после снятия архивов, его лучше вообще никогда нигде ни использовать. А на каком этапе он нам может понадобиться?. Если, как говорил Виктор, мы зальем архив SRAM на старый модуль (ROM-SRAM), установленный в новое ЧПУ, то и файл OPRMINF не нужен. Правильно же?
    • Slavdos
      Доброе. делюсь 1 внедренным китайцем. купили у ЛЛС MARVEL PRO 6000-3015.HGT , 2 шт , с автоматизацией. станки неплохие, интерфейс русские, достаточно дружелюбен. за автоматизацию зря переплатили, китай похоже в этом сильно уступает европе. из неожиданностей- резка воздухом дает неудаляемый грат, по сути необходимы зачистные станки.
    • Killerchik
      Я 4 комфорки на скоростной шпиндель поставил, в этом же качестве)))))
    • gudstartup
      я предлагал операторам чай кипятить чем больше и интенсивнее работают тем чаще чайку попить можно нельзя же чтобы тормозной кипятильник перегрелся
    • Viktor2004
      и этот человек смеялся когда я использовал кипятильник вместо тормозного резистора
    • gudstartup
      @mrVladimir если у вас нет копии oprminf обязательно сделайте
×
×
  • Создать...