自己紹介(M. K)

自己紹介(M. K) -OpenFOAMへの期待現職
熱対策
熱対策・
による解析
熱対策・揮発性物質対策の
揮発性物質対策のCFDによる
による解析
(流体解析経験約
流体解析経験約1年
流体解析経験約 年)
普段使う
普段使う市販ソルバー
市販ソルバーでは
ソルバーでは、
では、限界を
限界を感じています。。
じています。。
市販ソルバー
など)
市販ソルバーであまり
ソルバーであまり対応
であまり対応していない
対応していない問題
していない問題(ElectroHydroDynamicsなど
問題
など
•ソルバー
ソルバーを
ソルバーを自作すると
自作するとデバッグ
するとデバッグ+管理
デバッグ 管理に
管理に時間がかかる
時間がかかる
並列ジョブ
並列ジョブを
ジョブを沢山流したい
沢山流したい
•ライセンス
ライセンス料
ライセンス料が膨大に
膨大に、、、。
OpenFOAMなら
なら、
なら、そこそこの
そこそこの手間で
手間で並列計算が
並列計算が可能な
可能なソルバー
を自作できる
の使い方や、ソルバー改造
自作できる、
できる、と期待しています
期待しています。
しています。OpenFOAMの
ソルバー改造
方法を
方法を勉強したいと
勉強したいと思
したいと思っています。
っています。このあたり情報共有
このあたり情報共有ができれば
情報共有ができれば、
ができれば、と期
待しています。
しています。
ここ3
でEHDソルバー
ソルバー作成中
ここ3ヶ月、業務の
業務の空き時間で
時間でOpenFOAMで
ソルバー作成中です
作成中です。
です。
OpenFOAMのソルバー改造
チュートリアルの
ケース
チュートリアルのchargedWireケース
のEHD解析
解析
基礎方程式
ρ
∇ φ =−
ε0
∂ρ
+ ∇ ⋅ (ρU − kρ∇φ ) = 0
2
∂t
ρ
∂U
1
2
+ U ⋅ ∇U = −
∇p + ν∇ U +
∇φ
ρ air
ρ air
∂t
∇⋅U = 0
非圧縮性ソルバー
非圧縮性ソルバーicoFoam
ソルバー
+静電場解析
静電場解析ソルバー
静電場解析ソルバーelectrostaticFoam
ソルバー
でコードの
コードの自作は
自作は簡単。
簡単。
OpenFOAMのソルバー改造
チュートリアルの
ケース
チュートリアルのchargedWireケース
のEHD解析例
解析例
圧力分布
電荷密度分布
誘電体用の
誘電体用の境界条件モジュール
境界条件モジュールを
モジュールを作成中
∂ σ
∂ t
=
j ⋅ n , E
n
=
σ
2 ε
En
0
j
+ + + + σ