開発事例

1.　表面張力モデルの精度向上

標準のOpenFOAMでは曲率の計算精度が低い等の理由により、インクジェットの液滴など、微小スケールで表面張力が支配的となる計算では満足のいく計算精度が出せないことが多々あります。
これらの問題を解決するために、interFoamにおける曲率の計算精度向上を図り、 疑似流速の発生を抑える手法を組み込みました。

2.　多相流(interFOAM等)を含む着色機能の追加

一般的にパッシブスカラーと呼ばれる機能で、流体の一部を着色することで、 流れの様子を可視化する事がきます。
当社では、interFoamにスカラー移流拡散方程式を組み込み，液相（あるいは気 相）側だけに存在するスカラー移流拡散の解析を可能にしました。
また、スカラー依存で物性を変化させることも可能なので、kamalモデルや Macoskoモデルを組み合わせた熱硬化性樹脂の充填解析も可能です。

3.　二相流における熱計算機能の追加

interFOAMにエネルギー式の導入をおこない、温度場の解析を可能にしました。

4.　凝固・融解解析機能の開発

エンタルピー法を用いた凝固融解解析オプションを開発しました。凝固・融解解析では液体から固体への凝固、あるいは固体から液体への融解の相変化を解析することができます。
エンタルピー法ではエネルギー方程式に潜熱を考慮し、セルの液相割合を示す液相率（あるいは固相割合を示す固相率）が導入されます。
合金などの鋳造解析で生じるいわゆるMushy領域に対しては、液相率を用いて表わされるKozeny-Carman式から透過率を算出し、Mushy領域を多孔質体として扱います。
右の図は、初期温度が20℃の水が-10℃の壁面によって冷やされ、氷になっていく様子のアニメーションです。水は約4℃で密度が最大となるため、自然対流の解析には体積膨張率の温度依存性を考慮しないと精度のよい解が得られません。本オプションでは水の自然対流の解析を簡単に設定できるようにしてあります。
また、凝固・融解解析は自由表面解析（VOF法）と同時に実行できます。したがって、鋳造や溶接のシミュレーションも可能になっています。

5.　レーザー溶接・溶断解析ソルバーの開発

レーザー溶接では高パワーのレーザーを材料に照射します。高パワーレーザーを材料表面に局所的に照射することによって、材料の表面は急速に融解し蒸発します。溶けて蒸発した材料表面には反発力と呼ばれる高い圧力が作用し、材料にキーホールと呼ばれる狭く深い穴を生じさせます。レーザーは材料表面で反射を繰り返し、キーホールの底部へ到達します。これはまたレーザーから材料へのエネルギー伝達が増加することを意味します。このようにレーザー溶接の解析では、レーザーの材料表面での多重反射の考慮が必要になります。
（図：レーザーの多重反射）

これらレーザー溶接・溶断に固有の現象を捉えるために本ソルバーでは、以下の機能を実装しています。
　・蒸発反発力を考慮したCSF（体積力）モデル
　・レーザー多重反射の計算のための光線追跡アルゴリズム
また、本ソルバーはMPIによるコードの並列化をおこなっており、並列計算機での計算にも対応しています。

右の図（動画）は金属材料へのレーザー照射によるキーホールの生成過程の解析結果を示しています。
この解析では、エネルギー方程式、自由表面（VOF）、表面張力、凝固・融解が流れ場（質量保存式、運動量保存式）とともに解析されます。レーザー溶接・溶断解析ソルバーは、設定ファイル中で各解析機能のオン・オフに加えてさまざまな設定ができるソルバーとなっています。

6.　発泡材を考慮したソルバーの開発

シート成型等、発泡剤を考慮したソルバーを開発しました。反応率依存で、密度や体積、樹脂の粘度も計算する事ができるため、発泡現象を再現する事ができます。
また、空気の流動も解いているため、型内にトラップされるボイド位置も予測することができます。

7.　DEMと流体計算のカップリングソルバー開発

粒子計算(DEM)ソルバーと流体計算ソルバーの計算を一つのソルバーでおこない、流体と粒子の相互作用を考慮した、カップリングソルバーを作成しました。
DEMは粒子同士の衝突を解くため、流線を可視化する機能であるマーカー粒子・質量粒子と異なり、実際の砂のような粒子に近いふるまいを再現できます。
例えば、水中で堆積していく砂を模擬したり、多相流ソルバーと組み合わせて、水・空気・砂の解析も可能です。

8.　ろう付けを再現するオプション機能の開発

2層流解析用のソルバーであるinterFoamをカスタマイズして、ろう付けの計算が可能なソルバーを作成しました。
標準のinterFoamに湧き出し機能と多孔質体モデルを組み込み、さらに、毛管圧を考慮することで、溶融/未溶融状態のロウが混在する、接合面でのロウの流動を制御できるようにカスタマイズしました。
毛管圧などのパラメータは、計算条件の設定ファイルであるtransportPropertiesで変更可能な仕様にしました。
ロウ溶融プロセスのシンプルなモデル化と、非構造格子メッシュの利点を生かし、薄くて広い領域を伴う解析対象でも、実用的なＣＰＵ時間で計算可能です。

9.　新しい粘度式の追加

ユーザー様よりご提供いただいた任意の実験式を、新しい粘度モデルとして組み込みました。式中のパラメータは、transportPropertiesファイルで定義可能な仕様にしました。
複数の粘度モデルを追加可能で、どの粘度モデルを使用するかもtransportPropertiesファイルで選択できるようにしました。

10.　新しい境界条件追加

標準にはない新しい境界条件を追加しました。
境界条件が適切でないと、保存するべき量が無限に増え続けるなど、非現実的な結果となってしまうことがあります。
例えば、ある物理量が保存するように、特定の境界から流出する量とリンクして、他の境界から流入する量を決定する境界条件を作成しました。
分布を与える、実験データを反映させる、前ステップの結果から相対的に値を変える等、任意の境界条件も追加可能です。