Слошинг в CFX

[soklakov 1 475]

19 часов назад, karachun сказал:

Возможно это было вызвано численными ошибками - на некоторых шагах была плохая сходимость и в результате давления получались ниже чем должны были быть. Прямо как с подбрасыванием монетки на линейке в МКЭ: слишком большой шаг+контакт=грубое нарушение закона сохранения энергии.

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

[piden 2 874]

И как ты планируешь задачи сопрягать?

[soklakov 1 475]

1 минуту назад, piden сказал:

Твой пик низкого давления тоже весьма кратковременный - успеет ли он оказать существенное влияние на НДС бака?

 

В 03.03.2020 в 22:12, karachun сказал:

Конечно можно было бы кинуть всю массу на одну стенку и посчитать консервативный вариант. Или поискать упрощенные методики. Но клиент хочет именно слошинг.

 

1 минуту назад, piden сказал:

И как ты планируешь задачи сопрягать?

one-way-fsi вполне должно хватить.

[piden 2 874]

1 minute ago, soklakov said:

one-way-fsi вполне должно хватить.

Можно соединить блоки в ВБ, можно сделать каплинг, можно просчитать cfd и экспортировать давление на стенках на каждом шаге, потом зачитать их в транзиенте. Все это 1-way-fsi. Но все работают с разной скоростью.

Для сложной модели, по крайней мере, последний был многократно быстрее.

[soklakov 1 475]

1 минуту назад, piden сказал:

Но все работают с разной скоростью.

каков личный опыт эксплуатации?

[piden 2 874]

Just now, soklakov said:

каков личный опыт эксплуатации?

выше добавилось.

[karachun 2 038]

16 минут назад, piden сказал:

Да, я видел, и попробовал..

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

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

Попробую во флуенте сделать.

Класс, я пока курю мануалы, может к вечеру тоже сделаю простую модель. Можно будет потом "раскрыть карты" - сравнить кто что насчитал.

18 минут назад, piden сказал:

И как ты планируешь задачи сопрягать?

Только Фемап только костыли) 

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

В CFDшной задаче есть их аналог?

Может сохранения массы и количества движения. Запускал модель с очень малым шагом по времени - 1е-8 сек. И вместо хорошей сходимости у меня в баке были колоссальные осциляции давления и задача вообще разваливалась. Возможно это тоже из этой оперы.

[karachun 2 038]

@piden Задача нормально заработала если выделить три ядра и разбить домен на три горизонтальных партишна, получается что вся свободная поверхность находится внутри одного партишна.

 

[piden 2 874]

11 minutes ago, karachun said:

заработала

т.е. только начала считаться? А что будет, скажем, после 0,3 s?)

[karachun 2 038]

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

А что будет, скажем, после 0,3 s?)

Я просчитывал до 2 секунд (на одном ядре) и больше там никаких пиков по нагрузке не было. Так что все интересующие меня процессы происходят до 0.3 секунд, даже до 0.1.

Дальше течение просто медленно устаканивается.

Изменено 5 марта 2020 пользователем karachun

[piden 2 874]

20 minutes ago, karachun said:

Так что все интересующие меня процессы происходят до 0.3 секунд, даже до 0.1.

Ок) Но я к тому, что уже к 0,3 с граница раздела фаз должна попасть в другой партишинг - как себя тогда солвер поведет?

[karachun 2 038]

3 минуты назад, piden сказал:

к 0,3 с граница раздела фаз должна попасть в другой партишинг

За это время граница успевает размазаться по всему домену. В любом случае запущу сегодня расчет секунд на 5, просто потому что это красиво.

https://i.imgur.com/PYpql2W.mp4

[karachun 2 038]

@piden, ты оказался прав. Расчет развалился на 0.46 с и к тому же результаты не сохранились. Получается что для того чтобы быть уверенным в расчете надо запускать на одном ядре. 

[piden 2 874]

5 minutes ago, karachun said:

и к тому же результаты не сохранились.

Разве ты .trn не сохранял? Можно ведь рестартовать с последнего нормального?

On 3/5/2020 at 10:59 AM, piden said:

Попробую во флуенте сделать.

Эммм... я пока ковыряю, но для своей задачи)

[karachun 2 038]

@piden Эврика, я задал координаты (выбрал точку внутри объема воздуха) для Pressure Level и вот эта нижняя ветка графика просто исчезла. При этом результаты по давлению не поменялись! Я пока не понимаю почему оно работает и надо еще прогнать на более подробных сетках, но выглядит обнадеживающе.

[piden 2 874]

3 minutes ago, karachun said:

я задал координаты (выбрал точку внутри объема воздуха) для Pressure Level

в смысле, для buoyancy model?

 

UPD: уже нашел, в advanced option...

Не знал о ней))

[karachun 2 038]

@piden Я включил вот эту штучку в Solver Control.

[piden 2 874]

Just now, karachun said:

Я включил вот эту штучку в Solver Control.

А как теперь стало с multicore?)

Попробовал - вроде получше!)

[karachun 2 038]

8 минут назад, piden сказал:

А как теперь стало с multicore?)

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

[piden 2 874]

5 minutes ago, karachun said:

Надо будет проверить.

С партишинингом по-умолчанию:

На 4х ядрах пока полет нормальный (>400 tsteps), на 7 - вылет на 20м шаге