Jump to content
Sign in to follow this  
karachun

Слошинг в Fluent

Recommended Posts

karachun

Здравствуйте!

Считаю тестовую задачу. Слошинг для двумерной геометрии - квадрат 256х256 мм. Домен наполовину заполнен водой. Вторая половина - воздух.

Initial_Volume_Fraction.PNG

На домен действует постоянное ускорение свободного падения по вертикали и переменное ускорение по горизонтали - бак трясеться. переменное ускорение изменяется по такому закону, амплитуда - 8g.

Load.PNG

Я уже посчитал эту задачу в CFX, теперь я хочу сделать то же самое во Флюенте. Я ознакомился с примером в туториалах (там где считаем каплю чернил). Еще я нащел туториал на похожую тему - слошинг в бензобаке.

Но у меня еще остались вопросы.

1. Я сделал Transient Table для переменного ускорения, таблица успешно считывается но в поле ускорения (там где Expression, UDF и все такое) не отображается.

2. На вкладке General есть поле для задания ускорений и в Operation Conditions тоже есть поле для ускорений. Мне надо дублировать все ускорения и туда и туда или постоянное в Operation Conditions а переменное в General? В туториале по баку ускорение вообще задавалось только в Operation Conditions.

3. В CFX я задавал гидростатическое давление в жидкости, нужно ли мне его задавать во Флюенте или оно само посчитается из НУ, в туториале гидростатику не задавали?

4. Как можно изменить допустимые уровни Residuals? В окне на картинке они неактивны, я жал все кнопки подряд но они так и не активировались.

Residual Monitors.png

Прикладываю файлы проекта (2019R3), сетку, файл с Transient Table и туториал по баку.

Acceleration.txt box_2D.cas.gz box_2D.msh Tutorial_Fuel_Tank_Sloshing_Introduction.pdf

  • Нравится 1

Share this post


Link to post
Share on other sites


UnPinned posts
karachun

С некоторыми вопросами разобрался.

Зависимость нагрузки от времени задал с помощью ужасного Expreddion. @piden должно понравиться)

IF(t<=0.02[s],8*g*(t/0.02[s]),0*g)+IF(AND(t>0.02[s],t<=0.05[s]),8*g,0*g) +IF(AND(t>0.05[s],t<=0.07[s]),8*g*(1-(t-0.05[s])/0.02[s]),0*g)+IF(t>0.07[s],0*g,0*g)

Как я понял - ускорения на вкладке General потом дублируются в Operating conditions.

Результаты в модели странные - тоже есть большое падение давления и сам график весь в зазубринах. График похож на то что получалось в CFX пока я не задал Pressure Level. Только во флюенте я задал давление в Operating conditions.

1344384771_.thumb.PNG.2bc6bc06295e362da2ec0e1e853b5a87.PNG

Сам процесс перемещения воды похож, но в конце есть различия.

https://imgur.com/a/PLiWqXC

Прикладываю обновленный файл проекта.

box_2D.cas

Share this post


Link to post
Share on other sites
karachun

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

https://imgur.com/a/LMzp7Lr

При этом результаты хорошо совпадают с CFX, пока давление не становится отрицательным.

559619609_3.thumb.PNG.e444edc82bcb3feeb7bc7bd8ab548a21.PNG

box.zip

Share this post


Link to post
Share on other sites
karachun

@piden, @soklakov Я разобрался с нефизичными давлениями. Оказалось что данных из Operating Conditions недостаточно для задания правильного уровня давления в домене. Я добавил к комену маленькую трубочку с Pressure Outlet на конце и все заработало!

 

Снимок.PNG

Снимок1.PNG

  • Нравится 3

Share this post


Link to post
Share on other sites
MotorManiac
В 23.03.2020 в 18:48, karachun сказал:

При этом результаты хорошо совпадают с CFX

Так наверное считалось по одним и тем же моделям просто

Share this post


Link to post
Share on other sites
karachun
Posted (edited)

@MotorManiac Задача вообще не сильно чувствительна к "тонким" настройкам. Модели турбулентности, поверхностное натяжение, гомогенная/негомогенная постановка многофазности почти не влияют на результат.

А вот сама постановка задачи - наличие данных о уровне давления - критично. Без него получаются нефизичные отрицательные давления. В CFX без задания Pressure Level был отрицательный горб на графике давления и при этом с уменьшением сетки горб все рос без намеков на сходимость.

Pressure_Level.png

Graph.PNG

Edited by karachun

Share this post


Link to post
Share on other sites
MotorManiac
3 минуты назад, karachun сказал:

А вот сама постановка задачи

видимо да, основа основ)  Вы не знаете как средствами post cfd осреднить значение переменной на выбранном интервале (на графике, скажем, с 0.15с до 0.2с)? или это нужно вытягивать данные, например, в эксель и там уже смотреть ?

Share this post


Link to post
Share on other sites
karachun

@MotorManiac В CFX можно задать осреднение еще до запуска, в настройках мониторинга. А если все посчитано то по моему мнению проще вытянуть в ексель. В CFD-Post есть Table Biewer, по идее в нем должен быть базовый функционал екселя, но я в нем не работал. Еще в посте есть возможность писать пользовательские скрипты, но это для меня темный лес. Так что ексель/маткад будет удобнее всего.

  • Нравится 1

Share this post


Link to post
Share on other sites
karachun

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

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

В архиве приложил .cas файлы. Баки наполовину наполнены водой, в моделях уже есть Cell Register для инициализации.

https://imgur.com/a/wbrgVmH

Chart.png

1.png

2.png

simplified_water_tank_NoBulkheads.zip

Share this post


Link to post
Share on other sites
an_rushin

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

Share this post


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.

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.

Sign in to follow this  
  • Recently Browsing   0 members

    No registered users viewing this page.



  • Сообщения

    • Борман
      Что не подтверждается то? Если подтверждается на месте - тогда редкие единичные плюсы. Если не подтверждается - везут в Москву и оттуда приходит "большое" подтверждение.
    • SHARit
      Нет, не совпадает с данными статистики по областям.
    • Горыныч
      День добрый! Пришлите чертежи, модели, количество на почту konevsky(собака)inbox.ru
    • Тигр
    • Zergus
      Добрый день! Для того, чтобы все получилось, необходимы два условия: 1) расположение неподвижных кронштейнов определено (не перемещается) и геометрически стержень может быть одновременно концентричен отверстиям поворотных вкладышей; 2) кронштейны - это сборочные единицы вставленные в Вашу сборку. Если да, то Вам необходимо добавить сопряжение выравнивания ось-в-ось стержня и отверстия поворотного вкладыша второго кронштейна. Siemens NX назначит это сопряжение (а также еще какие-то, ранее созданные) ошибочным. Далее необходимо в навигаторе сборки раскрыть первый экземпляр сборки кронштейна, нажать правую кнопку мыши на поворотном вкладыше и выбрать команду "Перезаписать позицию". Тоже самое повторить с поворотным вкладышем во втором экземпляре сборки кронштейна. Стержень должен будет выровняться вместе с вкладышами. Более подробно можно почитать в справке . Хотя там без примеров, не очень понятно, на мой взгляд. Писал по памяти, проверить сейчас негде по понятным причинам, поэтому прошу не пинать, если где-то не точно указал.
    • Тигр
      Это новое противовирусное средство? Оформляй патент, пока фармацевтические компании непрознали...
    • Борман
      Пойду покрышку запалю..
    • Тигр
      Судя по сайту стопкоронавирус.рф провели 639000 тестов по всей России, там же данные по заболевшим, выздоровевшим и умершим. Но данных по переболевшим ранее не зарегистрированным как больные нет, только один раз Голикова ляпнула про выборку которую они сделали. По идее мой регион должно было накрыть первым, мы граничим с Китаем, так у нас много отдыхающих в Китае, грузопоток с Китаем большой, ну и самих китайцев приезжает много по туризму. Но как-то тихо всё проходит. ТО что Москва очень большой город и рассадник всякой гадости эт давно понятно, но обычно кто ближе всех тот попадает первым. А собственно из всего этого не ясно самое главное есть ли на самом деле та самая эпидемия, то есть порог числа заболевших сообщающих о её наступлении, и насколько быстро шириться коллективный иммунитет как показатель приближения спада эпидемии и реальной тяжести её протекания. У меня такое ощущение что у нас народу пострадает больше от мер принятых противодействовать вирусу чем от самого вируса.
    • Тигр
      Это какие публичные от бабы Дуси или из оф источников?
    • Борман
      Какой еще расчет ? Данные публичные, расчет на глазок :)
×
×
  • Create New...