CFD Module

数値流体力学(CFD:Computational Fluid Dynamics)のマルチフィジックスシミュレーション

CFD Module

2D近似流場とバッフル、乱流リアクタの3Dモデリングとの比較。

流体に関係する全てのアプリケーション

CFDモジュールは、複雑な流体モデリングを含むデバイスとシステムのシミュレーションを行うプラットフォームです。COMSOL Product Suiteの全てのモジュールと同様、 CFDモジュールは、グラフィカル・ユーザ・インタフェース(GUI)を介してモデルを入力できる定義済みフィジックスインタフェースを提供し、また、これらの入力をモデルの方程式を定式化するために使います。CFDモジュールが備えている特別なフィジックスインタフェースを使用して大半の流体の相をモデリング可能で、層流および乱流の領域、圧縮性、非等温、非ニュートン、二相流、多孔質流れなどの解析を含みます。CFDモジュールは、数値流体力学 (CFD)シミュレーションの一般的なツールとして使用可能で、また、流体が重要になるマルチフィジックスシミュレーションとしてCOMSOL Product Suiteにある他のモジュールとの連成が可能です。

CFD モジュールのGUIでは、モデリングプロセスの全ステップにアクセス可能です。下記のステップが含まれます。

*例えば、単相、二相、層流、あるいは乱流などのような流れの適切な状態の選択。 *モデルのジオメトリの作成やインポート。 *流体のプロパティの定義。 *流入口、流出口の項の追加、必要の応じて流体モデリングの重要な方程式の編集 *メッシュエレメントの選択、および、多様な位置におけるメッシュ密度の調整 *必要に応じて、ソルバの選択と調整


事例紹介

  • 混合の度合いがが粒子トレーシングによって観察可能な静的ミキサ内の流れ   混合の度合いがが粒子トレーシングによって観察可能な静的ミキサ内の流れ  
  • 回転機械インタフェースを使用した、撹拌されたバッフルミキサ内の流れ場 回転機械インタフェースを使用した、撹拌されたバッフルミキサ内の流れ場
  • 非等温流:換気装置の換気量のシミュレーションで、温度の等値面がプロットされます。モデリングでは、k-εモデルを使用して、非等温流や乱流をシミュレートします。 非等温流:換気装置の換気量のシミュレーションで、温度の等値面がプロットされます。モデリングでは、k-εモデルを使用して、非等温流や乱流をシミュレートします。
  • 二相流:回転流動層は、Euler-Eulerモデルを使用してモデリングされ、同時に固体微粒子は気体により流動化され、垂直ライザにより移送されます。 二相流:回転流動層は、Euler-Eulerモデルを使用してモデリングされ、同時に固体微粒子は気体により流動化され、垂直ライザにより移送されます。
  • 非ニュートン流:全て圧力の関数として、ポリスチレン溶液のせん断速度、動粘度、体積流量および等価なニュートン流の体積流量。 非ニュートン流:全て圧力の関数として、ポリスチレン溶液のせん断速度、動粘度、体積流量および等価なニュートン流の体積流量。
  • 乱流:翼の周辺の流れを示すSST乱流モデルのベンチマークの検討。図は、翼周辺の圧力係数の数値計算(ライン)と実験(マーク)によるそれぞれの結果を示しています。低圧側での実験データのみを表示しています。[参照:N. Gregory and C. L. O’Reilly, “Low-Speed Aerodynamic Characteristics of NACA 0012 Aerofoil Section, including the Effects of Upper-Surface Roughness Simulating Hoar Frost,” A.R.C., R. & M. No. 3726, 1970]

    乱流:翼の周辺の流れを示すSST乱流モデルのベンチマークの検討。図は、翼周辺の圧力係数の数値計算(ライン)と実験(マーク)によるそれぞれの結果を示しています。低圧側での実験データのみを表示しています。[参照:N. Gregory and C. L. O’Reilly, “Low-Speed Aerodynamic Characteristics of NACA 0012 Aerofoil Section, including the Effects of Upper-Surface Roughness Simulating Hoar Frost,” A.R.C., R. & M. No. 3726, 1970]

  • 非ニュートン流:全て圧力の関数として、ポリスチレン溶液のせん断速度、動粘度、体積流量および等価なニュートン流の体積流量。 非ニュートン流:全て圧力の関数として、ポリスチレン溶液のせん断速度、動粘度、体積流量および等価なニュートン流の体積流量。

全ての流体に最適なフィジックスインタフェースを

流体の多種多様な定義に適用するツールは搭載済みで、使いやすいフィジックスインタフェースとして提供されます。これらのインタフェースは、流体を定義する推進力、質量、そしてエネルギー方程式の不変性を定義し、他のフィジックスと連成するマルチフィジックスから構成されています。さらに、これらの方程式の安定フォームを定式化します。それらは時間依存問題もしくは、固定用の時間導関数の有限差分、空間の有限要素離散化を作成するため、COMSOLを使用します。安定した定式化は、流体種類のソルバコンフィグレーションとソルバ設定を示唆するフィジックスインタフェースにより流体プロパティとしての機能性および選択定義を適用します。オーダーメイドのフィジックスインタフェースは、以下の種類の流体において使用可能です:

  • Single-Phase Flow:  CFDモジュールは、流れ全ての速度領域をモデリングするための多種多様なNavier-Stokes方程式で解析します。低速度の流体、ゆっくりとした流れ(Stokes流れ)のモデリング、層流かつ弱圧縮性の流れ、そして乱流が含まれます。乱流は、k-ε, low-Reynolds k-ε, k-ω, SST (Shear Stress Transport), そしてSpalart-Allmaras 乱流モデルを含む、Reynolds-Averaged Navier-Stokes (RANS) 方程式を使用してモデリングされます。

    Navier-Stokes方程式の全ての変数と乱流モデリングの条件を操作するオプションを使用できます。例えば、連成した他のフィジックスインタフェースのモデルの変数を基にした方程式を含めることも可能です。乱流モデリングのソリューションプロセスの補助として、多数の追加ツールが内蔵されています。これらのツールには、壁関数 、境界のレイヤーメッシング、ハイブリッドメッシュの仕様、そしてメッシュ密度とメッシュ配置に適用する他のツールが含まれます。
  • Nonisothermal Flow:  伝熱との連成時において、温度‐浮力は、層流‐乱流内ではデフォルトで考察されます。CFDモジュールは、非等温流と共益熱伝達の定義済みマルチフィジックスインタフェースを含みます。このモジュールは、例えば、マッハ0.3以下の流れのような、弱圧縮性流を定義するために任意のマルチフィジックス連成が可能です。
  • Compressible Flow: CFDモジュールは、マッハ0.3以上の圧縮性流をモデリングすることも可能です。伝熱、圧縮、もしくは摩擦により生じた温度変化は、例えば、衝撃波のような圧縮性効果をもたらします。COMSOL Multiphysics内蔵のアダプティブメッシングは、衝撃波と液体流れのプロフィールの広範囲にわたる変化の解析に役立ちます。
  • Two-Phase Flow:  フィジックスインタフェースとそれがサポートする方程式により、二相流モデリングを自由に作成可能です。これが、2種以上の不混和流体を分離する移動境界の追跡を含む場合、CFDモジュールは、フェーズフィールド法とレベルセット法を使用します。

    CFDモジュールは、気泡流、混合モデル、およびEuler-Euler モデリング法を通じて、多数の粒子、液滴または気泡の懸濁液を含んでいる流れ記述可能な、分散二相流モデリングのフィジックスインタフェースも含んでいます。このEuler-Euler 法では、頻繁に衝突し、相間で相対速度が大きく変化する高濃度の気泡を解析可能です。伝熱モジュールには、空気の凝縮と湿度のモデリングインタフェースが含まれ、相の変化はCOMSOLMultiphysicsで定義済みのステップ関数を使用して記述されます。  
  • Porous Media Flow:  CFDモジュールを使用して、またDarcyのDarcyのBrinkmanの拡張を活用することにより、多孔質媒体内の単相流、二相流の輸送モデリングが可能です。Darcyの法則は、多孔質媒体の細孔が粘性の効果を打ち消すのに十分なほど小さい場合に適していて、一方でBrinkman方程式は粘性効果を考慮しており、流れは圧力差により発生します。内部条件では、フリーチャンネル流体流れと多孔質媒体の界面をモデリング可能です。  
  • Rotating Machinery:回転機械インタフェースは、流体ドメイン内で回転するプロペラの羽根や、ミキサーの翼などのような大規模なジオメトリ変更が可能で、回転部品を定義するモデリングツールを含んでいます。冷凍ローターインタフェースでは、シミュレーション中のジオメトリ変更を実行するのではなく、流体方程式に追加条件を課すことによって回転を解析します。このフィジックスインタフェースは、実際の回転現象を解く時よりも少ない計算リソースを使いながら、遠心力やコリオリ力を定常Navier-Stokes 方程式に追加します。そしてタービン、遠心分離機やミキサーのようなアプリケーションのモデリングによい近似値を提供します。旋回流インタフェースは、回転流れのモデリングでも使用可能です。このフィジックスインタフェースにおいて、面外旋回速度成分、完全な3Dモデルと比べて、計算量が低減された2Dジオメトリで軸対称モデリングが行えます。
  • Thin-Film Flow: CFDモジュールに含まれる専門的なフィジックスインタフェースは、例えば潤滑モデリングなどのように、2表面間もしくは、1表面上の薄層内に制限された液体流れ、もしくは気体流れをモデリングします。 
  • Non-Newtonian Flow:  CFD モジュールには、Carreauと、べき法則のモデリングが含まれるだけでなく、独自作成の方程式を定義可能で、ポリマーや他の非ニュートン性流体の粘性とせん断率を定義している外部データの読み込みも可能です。例えば、この方法で粘弾性モデリングを定義するとします。ステップ関数は、COMSOL Multiphysicsに内蔵され、例えばBingham流体を定義する時のように、流体プロパティの大規模変更もしくは突然の変更をモデリングする時に使用することができます。
  • Flow Through Thin Screens: 多孔板、グリル、そして金網を含むモデリングプロセスは、内蔵された薄層機能を使用して容易に作成可能です。これは薄層を通過する流れ、また薄層を介した層流もしくは乱流流れの影響を考慮した屈折と抵抗係数の相互関係を含んでいます。
  • Fluid Flow and Heat Transfer: CFDモジュールは、層流と乱流の流体流れを含む、完全連成された固体と流体の伝熱を定義する共役伝熱インタフェースを含みます。デフォルトでは、このモジュールは非等温流の解析を行い、たとえば伝熱モジュール内の表面輻射インタフェース、ジュール加熱、そして化学種輸送インタフェース内の反応熱などのフィジックスインタフェースとの連成が可能です。多孔質媒体の伝熱フィジックスインタフェースは、固体マトリクス伝導と流体相の伝導および対流を連成します。その一方で、液体により生じた蛇行パスと、これにより発生する散逸熱を組合わせます。
  • Reacting Flow:  希薄液と濃縮液内において化学種輸送により発生する層流-乱流、その両方と連成可能な専門的なインタフェースが含まれています。そしてChemical Reaction Engineering Moduleとの化学反応を定義しているインタフェースと連成することが可能です。

搭載済み関数の他にも、モデリング変数の任意関数式を定義すること、材料プロパティや境界条件を定義するソースや受信条件を移植するために編集フィールドに入力することも可能です。例えば、構造、温度、シェア率や他のモデリングとの連成の影響を定義している各自の関数を活用して、フィジックスプロパティを定義することも可能です。境界条件または、他のフィジックスモデリング変数の関数であり、他のフィジックスとも連成可能なソースや受信条件を定義可能です。これは、煩わしいユーザサブルーチンを使用することなく、GUI上で直接行うことができます。ドメインクロスセクション以上、もしくはドメインコントロール値内の平均的な値のようなモデリング変数の任意方程式を定義可能という特性により、レイノルズ、マッハ グラスホフ数の様々な代表値を使用する流体を特徴づけることができます。 

ソースや受信条件を追加したり、フィジックスプロパティをファンクションとして定義するステップとして、前回までにすでに定義されているフィジックスインタフェース内の基本的な方程式を操作することができます。それらは、流体定義の修正や他のフィジックスインタフェースとの規格外連成を実装するために行われます。

数多くのフィジックスインタフェースは、精巧な境界条件の定義をサポートします。定義されたSlip とNo Slipとの境界条件以外では、Sliding -Moving Wall条件をシミュレーションするためウォールを設定することが可能ですし、LeakingそしてOpenでも同様に可能です。そこでは、液体は自由な境界を作成すると想定されます。Wallファンクションとそれに関連する同調パラメータは、乱流モデルとして定義されます。出入口の速度や速度プロフィールは、Pressure, Stresses, Mass Flow条件とともに構成されます。これは、ある境界の出口フローから他の境界の入口フローにリンクされるPeriodic Flow /周期的流れ境界条件でも同様に構成されます。Periodic/周期的境界条件はは、反復的な単位セルアセンブリから構成されるジオメトリ内の単位セルをモデリングする時に大変便利です。  

マルチフィジックスと総合的シミュレーションの統合プラットフォーム

流体は多種多様なプロセスやアプリケーションの統合であり、他のプロセスの影響を考慮しながら、理解、最適化される必要があります。コンピュータハードドライブの効率的な冷却方法、加速度計の制動フィルムにおけるエネルギー分散方法、化学反応装置の異なる部分を通過した種の移動方法などが例として挙げられます。その中で流体は、他のフィジックスによって定義されたプロセスの主要因になっています。しかし実際には、電子装置から発生される熱は流体濃度に影響を及ぼします。また、加速度計の弾性は、流体を強制振動させ、その反応は化学成分を変化させます。潜在的に流体の流れに圧力を加えることになります。これが意味することは、全体的なプロセスの完全にそして正確な定義の影響を含めなければならないということです。

COMSOL MultiphysicsとCFDモジュールは、関連する全てのフィジックスとのシームレスな連成やGUI上のモデル方程式への無制限での直接的なアクセスも可能になり、プロセス定義をサポートします。両方間で自由に連成された流体構造連成 (FSI) 定式化を設定可能です。これにより、流体構造が変更される場合、また、この変更に対する構造反応が流体に影響を及ぼす場合、そのシナリオをモデリングすることが可能となります。 CFDモジュール内の全てのフィジックスインタフェースは、数値流体力学が考慮されるべきアプリケ―ションと一般的なプラットフォームを提供するCOMSOL Product Suite 内の他の全てのモジュールとを連成することも可能です。

COMSOLでは、CFDモジュールにおいて別の方法にて、フローモデリングを実装することが出来るモジュールをご提供します。その一方で、両方の方法からのメリットを活用して容易に連成することも可能です。一つの例は、Pipe Flow Moduleです。パイプ接合部分の定義として、エッジに沿った接線の平均速度と共に、エッジ要素を使用して発達の流れを2Dまたは3Dのパイプネットワークをモデリングします。これにより、タンクに接続中のパイプネットワーク内フローモデリングの実装が可能になります。また、これにより大規模な3Dメッシュが必要となるネットワーク内のパイプ接続部分のメッシングを回避することができます。COMSOLが備えている特徴は、完璧な2D、3Dジオメトリメッシュを作成し、パイプネットワークを接合するために、エッジから表面へそして表面からエッジへ、エッジから容積へそしてその逆へと、そのデータマッピングをシームレスに実装可能だということです。このようにして、パイプのネットワーク内の一定ユニットの数値流体力学プロパティを考慮しつつ、両方の接合部分の操作状態を調整することが可能になります。

全てのフィジックスは、標準的グラフィカルユーザインタフェースとワークフローを使用してモデリングされるため、CFDエンジニア達は、構造的、電気的もしくは化学的プロパティのようなのコンポーネントやプロセスにおいて、異なる特徴分析を行っているエンジニア達と、容易にコミュニケーションが可能です。ここですべきことは、ファイルを送付し、未解析分のフィジックスのスイッチを切り、他のフィジックスインタフェースを1つ、2つ追加し、そしてモデリングを続けるだけです。もちろん、これらの新しいフィジックスは、コンポーネントやプロセスのマルチフィジックスシミュレーション流体を定義したインタフェースとの連成を実現します。

最終的なCFDソリューションへのステップ

機器内もしくは、プロセス内で数値流体力学(CFD)をシミュレーションすることは、時として、最終的かつ正確なソリューションワークフローのステップになります。CFDモジュールは、そのワークフローの全ステップをサポートするための多様なツール、フィーチャー、設定そしてインタフェースが内蔵されています。

CAD Import Moduleもしくは、 LiveLink 製品のいずれかでは、サードパーティのCADソフトウェアでシミュレーションしようとするプロセス、コンポーネントやジオメトリのパートを取り込む作業のサポートをします。これらの製品は、CFDシミュレーションのメッシングを複雑にする、流れにとってあまり重要ではない小さなフィーチャーや構造を縮小するためのジオメトリ操作を可能にします。

CFDモジュール内で、3D CADで設計された形状を設定した後、即座に3Dシミュレーションを実行することはお勧めしません。COMSOL Multiphysicsでは、3Dジオメトリを使用した2Dモデリングワークスペース作成をサポートしてます。 代表的な横断面作成の2Dジオメトリでの作業を行うことにより、まずシミュレーション内の多数のパラメータに慣れることができます。3Dモデルの計算で必要となる大規模な計算リソースを使用することなく、以下のことが可能です:

*流体のプロパティがシミュレーション全体に及ぼす影響の調査 *使用すべき乱流モデルの決定 *適切なメッシングおよびレイヤーメッシングの境界設置の決定 *使用ソルバと設定の選択 *流体流れにおけるマルチフィジックス連成の影響の検討 *3Dモデル上で予測可能な正確性の見積もり    

システムの理解を進めながら、2Dモデリングから得た知識と最適化設定を使用することで、3Dシミュレーションを実装することが可能です。これは、対称もしくは線対称の3D CADで設計した形状の操作、全てを3Dでモデリングすることの回避、そして実質的な計算条件の削減を行う上でも大変役立ちます。

メッシングと強力なソルバにて柔軟さを提供するツール

メッシングは、デバイスやプロセス内での数値流体力学の大変重要なステップです。このメッシュでは、正確性が要求されますが、コンピュータリソースを流出するほどの微細さは要求されません。COMSOLマルチフィジックスは、液体シミュレーション上で最適なメッシュが実装できるよう、多様なツールをご提供しています。これには、構造化、非構造化およびスウェプトメッシュが含まれ、モデリング時におけるドメインのジオメトリサイズとそれらの割合のフレキシビリティを可能にします。また、フローの方向性への影響も含まれます。CFDモジュールでもまた、壁のような境界に沿ったメッシュの構造レイヤーを挿入するため、また全体としてハイブリッドメッシュを作成するために構造、非構造メッシュで覆った連成を実装するため、境界レイヤーメッシングを使用します 。

CFDモジュールでは、COMSOLマルチフィジックス内に含まれる、リニア、非リニア、時間依存、そしてパラメトリックソルバを利用します。2Dおよび小規模な3Dモデルを解析する直接的ソルバを含み、容易に収束可能な特性を備えています。また、大規模でより複雑なモデルの反復ソルバも含んでいます。プレコンディショニングおよびマルチグリッドソルバは、ソリューションを確定するために他のソルバと連成が可能です。クロスウィンドと流線拡散のような高性能ソルバ機能や平滑メソッドが使用可能で、これらの値は、他のソルバ設定と同様微調整することができます。CFDモジュールは、全く同一のシミュレーション内で別の要素を実装することも可能です。それは、低次要素を圧力のような変数解析で実装するために、またより高次要素を他の変数解析で実装するために適用することも可能ということです。

ソルバスキームによって、より近似値な解析プロセスの初期値を取得することができます。ある特定のモデリングドメインに内蔵されている層流場解析のように、フロー解析が容易に実装可能になるソルバスキーム設定が含まれています。また、乱流定義に対する最初の推測として、このソリューションを適用するソルバスキーム設定も含まれます。Frozen Rotor インタフェースを使用するソリューションは、回転モデリングドメインのシミュレーションの最初の予測数値として使用可能です。これにより多数の数値リソースを削減することができます。

CFDシミュレーションの正確性の拡張とデータ定義

CFDモジュールでは、下記のような流体のプロパティを算出します:フローパターン;圧力損失;流量-抵抗-揚力に依存する物体力;温度分布;系内の液体の多様な成分。さらに表面、流線、リボン、矢印、定性的粒子トレーシングのプロットだけでなく、アニメーションを含む定性的ポストプロセスを提供します。基本方程式の全てのパラメータを変数からのデータおよび追加条件は、他のパラメータもしくは変数に対する拡張やプロットにアクセスできます。これには、抵抗係数や揚力係数のような派生数値のポストプロセスも含まれます。CFDアプリケーションを解析する時に、粒子トレーシングモジュール内のフィジックスとカップリング、連成することで、粒子衝突やそれ自身の推進力を利用して、流体そのもの(Lagrange-Euler),およびお互いの粒子の影響の両方を考慮する事ができます。

Numerical Simulation-Based Topology Optimization Leads to Better Cooling of Electronic Components in Toyota Hybrid Vehicles

Modeling Optimizes a Piezoelectric Energy Harvester Used in Car Tires

The Burning Need for Modeling

Sea Floor Energy Harvesting

Silent Air Cooling: A New Approach to Thermal Management

Magnets Improve Quality of High-Power Laser Beam Welding

Cluster Simulation of Refrigeration Systems

A Smooth Optical Surface in Minutes

The Science of Water Screening

Restoration of Lake Water Environments

Heat Sink

Syngas Combustion in a Round-Jet Burner

Droplet Breakup in a T-junction

Circulated Fluidized Bed

Inkjet Nozzle - Level Set Method

Displacement Ventilation

Capillary Filling - Level Set and Phase Field Methods

Laminar Flow in a Baffled Stirred Mixer

Separation Through Electrocoalescence

Airflow over an Ahmed Body