цыган и открытопенка

[geladen]

К предыдущему: что я, собственно, загоношился. До сих пор считал я себе ДТК в Solidworks Flow Simulation. Считал-считал, насчитал даже кое-какую годноту, но в какой-то момент заебался. В солидворкце встроена без вариантов только одна модель турбулентности – k-ε – вполне стабильная и не очень требовательная к ресурсам, но известная залихватским враньём и паническими метаниями в случае сильных градиентов давления, срывов потока, и резких углов (всё, что есть ДТК). Какие-то достойные доверия результаты получаются только с очень тонкими сетками, от миллиона ячеек и более, и тут колом встаёт другая проблема солидвокцовой вычислительной гидродинамики – überговняная масштабируемость. Более 4 ядер вкладывать – бессмысленно, и на черновой просчёт первых 0.0002 секунд потока на вусмерть оптимизированной сетке нужно больше суток счётного времени.

Есть предел терпению; решил я хромой солидворкцовой считалке изменить с другой, молодой и ебливой, хорошо масштабируемой, и понимающей, например, k-ω-SST или то же -DES, и скормить ей мою геометрию, и посмотреть что будет.

Вариантов, на самом деле, не так много.

FlowVision, судя по документации -- офигенный пакет, там грамотные инженеры серьёзно поработали над автоматизацией, отчего получается в большой степени мышеклац. Но когда я осведомился за цены, мне выкатили об 27 килофранков за лицензию на 8 ядер; в габарит "самофинансирующегося хобби" не проходит вообще никак. К тому же, доступ просто так за бабло не дают, изволь взъебни недельный курс в мск, за отдельные деньги (ага, щас, в редкие понаезды в нерезиновую, мне там, типа, заняться нечем, а писать и читать толковую и полную документацию нынче не модно).

Т.е. понятно, что кустарь-маниак-одиночка – не их целевая аудитория. К чести фловижына, гражданин, который мне отписал за цены – в высшей степени конструктивный, вежливый и грамотный. Прекрасно и сразу поняв, что по итогам нихера я у них не куплю, вник тем не менее в вопрос, не спрыгнул с переписки и высказал пару ценных соображений. Жаль, что я не НИИ, захомячил бы не глядя.

Из других вариантов – Ансис в облаке, с почасовой или поденной оплатой (чай не авиаконструкторское кб, не так мне много надо). Но ансисы (и их перепродавцы) с места в карьер начали мощно ебать мозг, про то как мне необходимы консалтинговые услуги, и сертификация, и что дескать только после того как они взглянут на мою задачу (разумеется, небесплатно), то решат что мне надо (им виднее), и вообще – допускать ли меня до их überсчиталки. Тогда же, может быть, и цены скажут. Да пошли они сами на хуй, вот что.

OpenFOAM, по сути – единственный выход, к тому же единственный идеологически правильный (код должен быть свободен!). Там, в принципе, есть все ништяки и плюшки, и все решатели, и все модели турбулентности, и масштабируемость на сотни процессоров, и лазерный дальномер, и ерьбаг со стороны пассажира, и кофе с утра в постель умеет подносить, если правильно попросить, но задокументировано это всё говнястенько, приходится лезть в сорцкод, чтобы понять что к чему, и/или мудировать методом тыка без гарантии результата. Времени это занимает уйму, и нервы не успокаивает.

Вот и шарюсь теперь в поисках гуру из которого проистекали бы, аки пороховые газы из патронника, сакральные Знания. Из фрилансеров по опенфому только индусы попадаются "ёр прайс из май прайс май френд", по предыдущим работам судя – полные невменько, да чехи, судя по всему – грамотные, но 10 килоевро "за бонжур".

Похоже, раз вышел на путь левой руки – придётся идти до конца. Надо же когда-то прекращать симулировать понимание.

"Блядзь, Пеця", где ж на всё время найти?

( Comments |Comment on this)

Comments

[fon-rotbar.livejournal.com]

Средь моих френдов один БиГлебовский рубит в гидравлике и газодинамике...

[just-reaper.livejournal.com]

Я, конечно, может крамолу скажу, но купить (в ознакомительных целях и на короткое время, конечно же) FlowVision на торрентах - не? Оно вроде как есть там.

[geladen.livejournal.com]

не то, чтобы я так чтил, прости Господи, "авторское" "право" (моё отношение см. https://sashnik.livejournal.com/180773.html -- срать я на "авторское" "право" хотел говном из жопы), но крутить нечто вялокоммерческое на пиратском софте, вкладываться в изучение, подсеть изолировать, крутить непатченные системы -- не уверен, что оно мне так уж необходимо.
если бы по деньгам в габарит проходило, непременно скачал бы сперва пиратку, чтобы убедиться что оно мне надо, а так...

[glockmeister.livejournal.com]

поговори с одним толковым парнем профильным
он в свое время считал и разрабатывал дтк, которые пошли в серию
может подскажет чего

контакт скину

[geladen.livejournal.com]

мерси! буду очень благодарен.

[(Anonymous)]

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

Если хочется, чтобы расчет не был медленным и печальным, то надо смотреть не только в стороны масштабируемости софтины, но и в сторону того, умеет ли она в динамическую адаптацию расчетной сетки, при чем не только с измельчением, но и с загрублением обратно (и желательно сразу по нескольким критериям)

[geladen.livejournal.com]

спасибо большое!

re: считать вообще в невязкой постановке --
мерси! попробую, посмотрю на что похожи результаты, и насколько разнится счётное время.

а когда дело дойдёт до турбулентности, есть ли какие-нибудь критерии оценки валидности модели для конкретного случая? что думаете про k-ω-SST? k-ω-SST/DES? иное?

re: умеет ли она в динамическую адаптацию расчетной сетки --
о! как раз об этом читал в мануале опенфома: dynamicRefineFvMesh does not morph the mesh shape. Instead, it performs topological refinements to the mesh, based on the values of specified fields. This is useful for cases like sonic shockwaves. Situations where you have a large region with relatively little change in field variables, and a small section with rapid and extreme change in field variables.

то есть оно в принципе имеется -- и измельчение и загрубление. теперь попытаюсь понять как это задать на градиенте давления (мануал сходу заявляет -- только на одном параметре), с синтаксисом, похоже, есть тонкости.

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

[(Anonymous)]

Да, потом можно попробовать k-w SST, а потом что-то более сложное гибридное. Простой DES сейчас уже не носят )) в моде более сложные варианты - DDES (Delayed DES), IDDES (Improved Delayed DES) и т.д. Но они все по большому счету отличаются тем, как идет определение «а где у нас тут погранслой?» чтоб ненароком не влезть в него LES частью. Но для ДТК КМК хорошее описание погранслоя не критично (если тепловое состояние не трогать), так что по идее пойдет любой вариант, что есть в софтине и с чем нормально считается.

С адаптацией сетки будет масса интересных моментов. Хотя градиент давления кажется очевидным выбором, с ним всё не так однозначно. Во-первых, численное значение градиента зависит от того, какая сетка у нас уже есть. Во-вторых, ударные волны от разных щелей и отверстий могут быть сильно разной интенсивности. И отсюда возникает вопрос - какие пороговые значения ставить для измельчения и загрубления. Поэтому народ пишет всякие ученые статьи на тему «мы придумали новый shock detector, который детектит лучше, чем раньше». Можно поискать в Google Academy по shock detection - там мноооого чего будет. Кроме градиента давления еще используют градиент градиента (вторую производную) давления, градиент числа Маха, произведение градиента давления на текущий размер элемента.

[geladen.livejournal.com]

ого, спасибо огромное! вот эти два абзаца сэкономили мне несколько суток колупания по тырнетам и учебникам, и, наверное, не одну неделю счётного времени.

посмотрел в доках опенфома, в модных вариантах DDES унд IDDES идёт только Spalart Allmaras, который, если я правильно понимаю, заточен под аэродинамику и не очень дружит со срывом потока. к-омега-сст-дес, однако, имеется -- на него пока и буду ориентироваться, мерси.

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

re: выбор параметра для адаптации сетки --
и вновь -- спасибо за науку!
из того, что сейчас понимаю, в опенфоме "из коробки" без допиливания сорцкода учесть размер сетки не получится. но -- (а) вторая производная давления, (б) градиент числа маха -- очень ценно.

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

[geladen.livejournal.com]

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

[drw-grinder.livejournal.com]

Я сварщик не настоящий и в CFD не умею, но есть вот https://www.autodesk.com/products/cfd/overview с бесплатным trial на 30 дней который в свою очередь интегрируется с CAD в виде https://www.autodesk.com/campaigns/fusion-360-for-hobbyists который сам по себе бесплатный для хобби и мелкого производства. Кадовская часть опробована man-of-motley и он вроде ей доволен. Насколько оно умеет именно то что вам надо - вопрос открытый. Но за попробовать денег не берут.

[geladen.livejournal.com]

в принципе, судя по описанию, там с функционалом ок, но после 30 дней -- 12K$ в год, увы.

[drw-grinder.livejournal.com]

Для студентов вроде бесплатно. https://www.autodesk.com/education/free-software/featured#free-software Неужели нет скубента в окружении чтобы на него зарегистрировать?