FuruRok

2.6K posts

FuruRok

FuruRok

@FuruRok

有限要素法ライブラリMFEMを用いた自作の三次元並列流体構造連成解析コードを趣味で開発しています。Space-Time法等、有限要素法(FEM;Finite Element Method)関連手法に興味あり。電機メーカー研究員。自分も勉強しながら、今どきのFEMの学び方をお伝えできればと思います。

เข้าร่วม Aralık 2022
982 กำลังติดตาม755 ผู้ติดตาม
FuruRok
FuruRok@FuruRok·
@dmbrkp_ がっ…がんばります
日本語
0
0
1
22
ダムP🌸🌟
ダムP🌸🌟@dmbrkp_·
Balanced-forceできたしば!うにょうにょしば!
日本語
2
0
13
252
FuruRok
FuruRok@FuruRok·
@dmbrkp_ こんなのがあったんですね、素晴らしいです!まさに必要としていたものでした。 ありがとうございます‼️
日本語
1
0
1
21
FuruRok
FuruRok@FuruRok·
@dmbrkp_ 私も勉強します!(非構造格子でできるかな)
日本語
1
0
0
23
FuruRok
FuruRok@FuruRok·
@dmbrkp_ 界面がすごくきれいです!
日本語
1
0
0
24
FuruRok
FuruRok@FuruRok·
@dmbrkp_ 2D-VOF-PLIC-構造格子-差分法 ですか?
日本語
1
0
0
35
FuruRok รีทวีตแล้ว
広江 克彦
広江 克彦@eman1972·
新刊発売予定です!|EMAN 『量子力学の〈新常識〉ー 意識・実在・情報が再定義する物理学ー』の紹介です。 note.com/eman/n/n33ede0…
日本語
0
45
111
30.5K
FuruRok
FuruRok@FuruRok·
@likehhkb 界面再構成が、非構造格子では難しいということでしょうか? 構造格子ならOK?
日本語
1
0
0
15
FuruRok รีทวีตแล้ว
Mathematica
Mathematica@mathemetica·
Sink (left, div < 0; everything collapses inward), Source (middle, div > 0; pure explosion outward), Incompressible (right, div = 0; swirls forever with zero net loss). Textbooks never made the flow rates feel this alive.
English
17
238
1.4K
34.9K
FuruRok
FuruRok@FuruRok·
たまには並列数変えて計算流そう、これ大事。
日本語
0
0
2
150
FuruRok
FuruRok@FuruRok·
ALE-FSI法、数ヶ月前には正しく動作していたものが、同じ条件とソルバーで途中で発散するようになった。色々調べると並列数を変えるとうまくいかなくなる場合があることが分かった。 ALEは完成したと思っていたがまだ修正箇所が残っているらしい。🫠
日本語
0
0
3
401
FuruRok รีทวีตแล้ว
ksa 🏴‍☠️
ksa 🏴‍☠️@kosa12m·
Best paper I've read so far this month: All elementary functions (sin, cos, tan, exp, log, powers, roots, hyperbolic functions, π, e, and even basic arithmetic) can be generated from just one binary operator: eml(x, y) = exp(x) − ln(y) …plus the constant 1.
ksa 🏴‍☠️ tweet media
English
88
1.7K
11.2K
1.4M
FuruRok
FuruRok@FuruRok·
入口流速境界を与えるところで、プロファイル計算に用いる寸法をハードコードしてしまっている部分があり、メッシュの変化に対応できていない箇所があり修整。理論解との比較が良好になった。🎉
日本語
0
0
1
180
FuruRok
FuruRok@FuruRok·
強形式の∇・uを小さくするには、grad-divの係数γを大きくする方法があり、実際にそうなることが確認できた。 しかし、γを大きくすると、方程式系が硬くなることや、有限要素解がずれてしまうという問題がある。これについても確認できた。 根本的には、有限要素空間を現状のH1から変更する等が必要。
日本語
0
0
1
120
FuruRok
FuruRok@FuruRok·
∇・u=0が満たされない件、これまで強形式残差の観点でモニターしていたが、FEMではこれは必ずしも満たされないので∇・uの値が大きいことが分かった。離散連続式残差の観点でモニターすると、きちんと小さくなっていることが確認できた。収束判定に使うなら後者と思われる。
日本語
1
0
8
377
FuruRok รีทวีตแล้ว
shenxiao
shenxiao@Shenxiao123971·
我要用这个视频告诉孩子们:「这个飞行姿势在空气动力学上太糟糕了」
中文
859
1.1K
20.2K
7.2M
FuruRok
FuruRok@FuruRok·
∇・uは入口流速がramp状に立ち上がる時間帯だけ大きく、rampが終わると1e-7くらいに小さくなるので、rampが終わってから収束判定を行うのが良さそう。
日本語
0
0
0
99
FuruRok
FuruRok@FuruRok·
∇・u(のノルム)、入口境界条件をポアズイユ流れのプロファイルにしたら劇的に下がった。🎉 強制的に与える入口プロファイルが、∇・u=0を満たしていなかった。 収束判定に使えるほど小さくはないけど。。。
FuruRok@FuruRok

∇・u(のノルム)、見直した結果1/10になった。🎉 ・出口圧力全体を0に拘束していたのを1点のみにした。 ・grad-divを入れていたつもりだったが、残差側には入っていたがヤコビ行列側には入っていなかったのを修整。 ただ、収束判定に使うにはまだ大きいので引き続き見直しする。

日本語
1
0
1
452
FuruRok
FuruRok@FuruRok·
∇・u(のノルム)、見直した結果1/10になった。🎉 ・出口圧力全体を0に拘束していたのを1点のみにした。 ・grad-divを入れていたつもりだったが、残差側には入っていたがヤコビ行列側には入っていなかったのを修整。 ただ、収束判定に使うにはまだ大きいので引き続き見直しする。
FuruRok@FuruRok

自作コード(FEM/VMS)では∇・u=0があまり満たされないんだよなぁ。初心に帰って見直すか。

日本語
0
0
2
741