ACU-T:6102 AcuSolve - EDEMの双方向連成を使用した風力選別機での粒子分離

このチュートリアルでは、AcuSolveおよびEDEMを使用して基本の双方向連成(二方向)のシミュレーションを設定して実行するワークフローを紹介します。このチュートリアルを開始する前に、HyperWorksの入門チュートリアルであるACU-T:1000 HyperWorksユーザーインターフェースおよびACU-T:6100 Altair EDEMを使用した風力選別機での粒子分離をすでに完了している必要があります。また、HyperWorks CFDAcuSolveEDEMの基本を理解しているものとします。このシミュレーションを実行するには、ライセンス供与済みバージョンのHyperWorks CFDAcuSolve、およびEDEMにアクセスできる必要があります。

このチュートリアルを実行する前に、HyperWorksCFD_tutorial_inputs.zip<Altair_installation_directory>\hwcfdsolvers\acusolve\win64\model_files\tutorials\AcuSolveから作業ディレクトリにコピーします。 ACU-T6102_windshifter.hm をHyperWorksCFD_tutorial_inputs.zipから取り出します。

問題の説明

ここで解析する問題を、acu_6102_intro_cfd_r.htm#acu_6101_intro_cfd_r__fig_cjm_1kg_kmbに概略的に示しています。これは風力選別機のモデルで、底部の入口から非常に高速(20m/sec)で領域内に空気が入り、出口から出ていきます。粒子は粒子入口から領域内に取り込まれ、軽い粒子は空気によって遠くへ運ばれ、重い粒子は底部の開口部から出ていきます。


図 1.

このモデルは45度に曲がった円筒管で構成されます。管の半径は0.25mで、粒子の入口は管全体の長さの中ほどに配置されます。

AcuSolve-EDEMの一方向連成シミュレーションのワークフローを以下に示します。


図 2.

これに応じて、このチュートリアルは次の2つの部分で構成されます。

  1. AcuSolveの設定と形状のエクスポート
  2. EDEMの設定とシミュレーション

AcuSolveモデルはHyperWorks CFDを使用して設定します。AcuSolveの設定が完了したら、EDEMのデック(形状を含む)がHyperWorks CFDからエクスポートされます。この入力デックはEDEMで開かれ、EDEMの設定を実行するために使用されます。EDEMデックの設定の終了後、連成シミュレーションを起動します。

EDEMシミュレーションで使用される2つの異なるBulk Materialsとそれぞれのプロパティを以下に示します。

名前 密度(kg/m3 粒子のサイズ(m) 個々の粒子の平均重量(kg) 生成率(1秒あたりの粒子数)
重い粒子 900 0.03 0.1 100
軽い粒子 100 0.03 0.01 100
粒子の流体抗力は、Ergun-Wen-Yuモデルを使用して計算されます。AcuSolveは粒子の位置と速度の情報をEDEMから受け取り、粒子の抗力を計算して流体フィールドを更新します。この流体力情報がEDEMに戻され、これを使用して粒子の位置と速度が更新されます。このループがシミュレーションの最後まで繰り返されます。


図 3.

パート1 - AcuSolveシミュレーション

HyperWorks CFDの起動とHyperMeshデータベースのオープン

  1. Windows のスタートメニューからスタート > Altair <バージョン> > HyperWorks CFDをクリックして HyperWorks CFDを起動します。
  2. HomeツールのFilesツールグループからOpen Modelツールをクリックします。


    図 4.
    Open Fileダイアログが開きます。
  3. モデルファイルの保存先ディレクトリを参照します。HyperMeshファイルACU-T6102_windshifter.hmを選択し、Openをクリックします。
  4. File > Save Asをクリックします。
  5. windshifter_bidirectionalとして、他の入力ファイルと同じディレクトリに保存してください。
    このディレクトリが解析用ディレクトリになり、シミュレーションに関連するすべてのファイルがこの場所に保存されます。

形状の検証

Validateツールは、モデル全体をスキャンし、サーフェスおよびソリッド上でチェックを実行して、形状に不具合(フリーエッジ、閉じたシェル、交差、重複、スライバーなど)があればフラグ付けします。

シミュレーションの物理パートに集中するために、このチュートリアルの入力ファイルにはすでに検証済みの形状が含まれています。GeometryリボンのValidateアイコンの左上隅に青色のチェックマークが表示されていることを確認します。これは、形状が有効で、フロー設定に進めることを示しています。


図 5.

流れのセットアップ

一般的なシミュレーションパラメータの設定

  1. FlowリボンからPhysicsツールをクリックします。


    図 6.
    Setupダイアログが開きます。
  2. Physics modelsの設定で、Multiphase flowラジオボタンを選択します。
  3. Multifluid typeをBidirectional EDEM Couplingに変更します。
  4. Eulerian materialドロップダウンメニューをクリックし、リストからMaterial Libraryを選択します。
    材料ライブラリに新しい材料モデルを作成できます。
  5. Material LibraryダイアログでEDEM 2 Way Multiphaseを選択し、My Materialタブに切り替えて、をクリックして新しい材料モデルを追加します。
  6. マイクロダイアログで、左上隅にあるEDEM Bidirectional Materialをクリックして、名前をAir-Particleに変更します。
  7. Carrier fieldをAirに設定します(まだ設定されていない場合)。
  8. 下図に示すように、抗力モデルと抗力係数を設定します。


    図 7.
  9. 材料モデルマイクロダイアログを閉じて、Material Libraryダイアログを閉じます。
  10. Setupダイアログで、Eulerian MaterialをAir-Particleに設定します。
  11. Time step sizeを0.001に、Final timeを1にそれぞれ設定します。TurbulenceモデルでSpalart-Allmarasを選択します。以下のようにGravityが0, 0, -9.81に設定されていることを確認します。


    図 8.
  12. Solver controls設定をクリックして、Minimum stagger iterationsを2に、Maximum stagger iterationsを4にそれぞれ設定します。


    図 9.
  13. ダイアログを閉じてモデルを保存します。

材料プロパティの割り当て

  1. FlowリボンからMaterialsツールをクリックします。


    図 10.
  2. Air-Particleが材料として割り当てられていることを確認してください。
  3. ガイドバーで、 をクリックしてツールを終了します。

流れ境界条件の定義

  1. FlowリボンのProfiledツールグループからProfiled Inletツールをクリックします。


    図 11.
  2. 下図でハイライトされているInlet面をクリックします。マイクロダイアログで、Average velocityに20m/sという値を入力し、Carrier fluid volume fractionを1.0に設定します。


    図 12.
  3. ガイドバーをクリックすると、コマンドが実行されてツールが終了します。
  4. Outletツールをクリックします。


    図 13.
  5. 下図でハイライトされている面を選択し、ガイドバーをクリックします。


    図 14.
  6. No Slipツールをクリックします。


    図 15.
  7. 下図でハイライトされている面を選択します。


    図 16.
  8. 境界の凡例で、WallをダブルクリックしてParticle_inletという名前に変更します。
    このサーフェスを参照形状として使用し、EDEMでparticle factoryを作成します。このため、このサーフェスを別の形状グループに配置します。
  9. ガイドバーをクリックしてします。
  10. モデルを保存します。

メッシュの生成

  1. MeshリボンからBatchツールをクリックします。


    図 17.
    Meshing Operations ダイアログが開きます。
  2. Average element sizeを0.3に設定します。


    図 18.
  3. Meshをクリックします。
    Run Statusダイアログが開きます。解析が終了すると、ステータスが更新され、ダイアログを閉じることができます。
    ヒント: メッシュジョブを右クリックし、View log fileを選択してメッシングプロセスの概要を表示します。
  4. モデルを保存します。

節点出力の設定

メッシングが完了すると、自動的にSolutionリボンに移動します。

  1. SolutionリボンからFieldツールをクリックします。


    図 19.
    Field Outputダイアログが開きます。
  2. Write Initial Conditionsを選択します。
  3. Time step intervalを10に設定します。


    図 20.

ソルバーデックのエクスポート

  1. メニューバーで、File > Export > Solver Deckに移動します。
  2. ファイルにwindshifter_bidirectionalという名前を付け、ファイルタイプにAcuSolve (*.inp)が選択されていることを確認します。
  3. Saveをクリックします。
AcuSolveの入力ファイルとEDEMの入力デック(形状セクションを含む)が作成されます。このEDEMデックを使用して、DEMシミュレーションデックを作成します。

パート2 - EDEMシミュレーション

Windowsのスタートメニューからスタート > Altair 2021 > EDEM 2021をクリックして、Altair EDEMを起動します。

EDEMの入力デックを開く

前述したとおり、AcuSolveシミュレーションの起動時に、HyperWorks CFDが一連のEDEMファイルを作業ディレクトリ内に作成しました。このEDEMの入力デックを開き、DEMシミュレーションを設定します。

  1. EDEMのCreatorタブで、File > Openに移動します。
  2. ダイアログで、AcuSolveの作業ディレクトリに移動し、EDEMフォルダにあるwindshifter_bidirectional.demファイルを開きます。
    形状が読み込まれます。
  3. Creator TreeでEnvironmentタブをクリックし、Auto Update from Geometryチェックボックスをオフにして再度オンにし、境界内に形状をフィットさせます。


    図 21.


    図 22.

Bulk Materials と Equipment Materialの定義

この手順では、重量および軽量のBulk Materials とEquipment Materialの材料モデルを定義します。

  1. Creator Treeで、Bulk Materialを右クリックしてAdd Bulk Materialを選択します。
  2. 材料の名前をHeavyに変更します。
  3. Creator Treeで、Solids Densityプロパティを900kg/m3に設定します。
    このチュートリアルの他のプロパティにはデフォルト値を使用します。


    図 23.
  4. Interactionの下のをクリックして、重量粒子間の衝突に関する相互作用のプロパティを定義します。ダイアログでOKをクリックします。
  5. Creator Treeで、Heavyを右クリックして、Add Shape from Library > Dual Sphere Shapeを選択します。
  6. 粒子の名前をHeavy particleに変更します。
  7. Heavy particleの下のPropertiesをクリックします。
  8. Heavy particleのSpheresパネルで、両方の球体のPhysical Radiusを0.03mに設定し、Enterを押します。


    図 24.
  9. Creator Treeで、Calculate Propertiesをクリックします。


    図 25.
  10. Creator Treeで、Bulk Materialを右クリックしてAdd Bulk Materialを選択します。
  11. 材料の名前をLightに変更します。
  12. Creator Treeで、Solids Densityプロパティを100kg/m3に設定します。
    このチュートリアルの他のプロパティにはデフォルト値を使用します。


    図 26.
  13. Interactionの下のをクリックして、重量粒子間の衝突に関する相互作用のプロパティを定義します。ダイアログで、Heavyを選択してOKをクリックします。
  14. もう一度をクリックして、軽量粒子間の衝突に関する相互作用のプロパティを定義します。ダイアログで、Lightを選択してOKをクリックします。
  15. Creator Treeで、Lightを右クリックして、Add Shape from Library > Dual Sphere Shapeを選択します。
  16. 粒子の名前をLight particleに変更します。
  17. Light particleの下のPropertiesをクリックします。
  18. Light particleのSpheresパネルで、両方の球体のPhysical Radiusを0.03mに設定し、Enterを押します。


    図 27.
  19. Creator Treeで、Calculate Propertiesをクリックします。


    図 28.
  20. Creator Treeで、Equipment Materialを右クリックしてAdd Equipment Materialを選択します。名前をSteelに変更します。
  21. Densityを7800kg/m3に設定します。
  22. Interactionの下のをクリックして、重量粒子間の衝突に関する相互作用のプロパティを定義します。ダイアログで、Heavyを選択してOKをクリックします。
  23. もう一度をクリックして、軽量粒子間の衝突に関する相互作用のプロパティを定義します。ダイアログで、Lightを選択してOKをクリックします。
  24. モデルを保存します。

Particle Factoryの作成

  1. Creator TreeタブでGeometriesを展開します。次に、Particle inletサーフェスグループを右クリックして、Copy Geometry > Single Copyを選択します。
  2. 新しい形状セクションの名前をParticle_factoryに変更します。
  3. Creator TreeでDefault Wall形状セクションを右クリックして、Merge Geometry(s)を選択します。
  4. Merge Geometryダイアログで、Particle_inletを選択してOKをクリックします。
  5. Default Wall形状セクションをクリックし、TypeをPhysicalに、MaterialをSteelに設定します(まだ設定していない場合)。


    図 29.
  6. Creator Treeで、Inletセクションをクリックし、TypeをVirtualに変更します。
  7. 同様に、OutletおよびParticle_factoryセクションのTypeをVirtualに変更します。
  8. Particle_factoryで、Transformをクリックします。X-Positionを0.035mに設定します。


    図 30.
    注: Opacityの値を0.2に設定し、パイプ形状内部の変換後のサーフェス位置を表示します。

    これは、粒子が流体領域内に生成されていることを確認するためです。

Particle Factoryの定義

これでBulk Materials 、形状セクション、Equipment Materialが定義されたので、次は粒子を生成するためParticle Factoryを作成する必要があります。Bulk Materials ごとに1つのファクトリを作成します。

  1. Creator Treeで、Particle_factoryを右クリックしてAdd Factoryを選択します。
  2. 新しいファクトリの名前をHeavy factoryに変更します。
  3. 下図に示すように粒子の生成パラメータを設定します。


    図 31.
  4. Velocityの横のをクリックし、X-velocityを1m/sに設定してOKをクリックします。
  5. 同じパラメータを使用して手順1~4を繰り返し、Light factoryという別のファクトリを作成します。ただし、MaterialにはLightを設定します。

Environmentの定義

この手順では、EDEMのシミュレーション領域の範囲と重力加速度の方向を定義します。

  1. Creator Treeで、Environmentをクリックします。
  2. Auto Update from Geometryのチェックボックスをアクティブにします(未選択の場合)。
    移動粒子が領域(環境)の境界面に接触すると、シミュレーションから取り除かれます。
  3. Gravityをアクティブにして、Z値を-9.81m/s2に設定します。
  4. EDEMデックを保存します。

シミュレーション設定の定義

  1. 左上隅のをクリックして、EDEMのSimulatorタブに移動します。
  2. Simulator Settingsタブで、Time Integration schemeをEulerに設定し、Auto Time Stepチェックボックスをアクティブにします(未設定の場合)。
  3. Total Timeを1sに設定し、Target Save Intervalを0.01sに設定します。
  4. Cell Sizeを4Rminに設定します。
    一般的に、最適なセルサイズとして3~6Rminの範囲の値が推奨されます。EDEMのセルサイズがシミュレーションの精度に影響を与えることはありません。実行時間にのみ影響します。
  5. Selected EngineをCPU Solverに設定し、可用性に基づいてNumber of CPU Coresを設定します。


    図 32.
  6. シミュレーション設定の定義が完了したら、モデルを保存します。

連成シミュレーションの実行

  1. EDEMCoupling Serverをクリックして、連成サーバーを起動します。


    図 33.
    連成サーバーがアクティブになると、アイコンが変化します。


    図 34.
  2. HyperWorks CFDに戻ります。
  3. SolutionリボンからRunツールをクリックします。


    図 35.
    Launch AcuSolveダイアログが開きます。
  4. additional arguments欄に以下のテキストを入力します。-tlog -lprobe.
    これにより、AcuSolveに、シミュレーションの進行状況に応じた結果のモニターに使用されるAcuTailAcuProbeウィンドウの起動を指示します。
  5. Parallel processingオプションをIntel MPIに設定します。
  6. オプション: プロセッサーの数を、利用環境に合わせ、4または8に設定します。
  7. Default initial conditionsを展開し、Pre-compute flowチェックボックスをオフにして、velocity valuesを0に設定します。Pre-compute Turbulenceチェックボックスをオフにします。
  8. RunクリックしてAcuSolveを起動します。


    図 36.
    AcuSolveEDEMの連成が成功すると、最初の時間ステップの前にAcuSolveログファイルにメッセージが出力されます。


    図 37.
    解析が進むに伴い、AcuTailウィンドウとAcuProbeウィンドウが自動的に開きます。AcuTailウィンドウには、残差指標と時間増分指標の情報がシミュレーションの進捗として出力されます。最後にシミュレーションのサマリーが出力され、シミュレーションの完了が示されます。


    図 38.
  9. 解析が完了した後、AcuTailウィンドウとAcuProbeウィンドウを閉じます。

結果の分析

AcuSolveのポスト処理

  1. 解析の完了後、Postリボンに移動します。
  2. メニューバーで、File > Open > Resultsの順にクリックします。
  3. 作業ディレクトリでAcuSolveログファイルを選択し、ポスト処理の結果を読み込みます。
    ソリッドとすべてのサーフェスがPostブラウザに読み込まれます。
  4. view cubeでLeft面を選択し、モデルをx-z平面に合わせます。


    図 39.
  5. Postブラウザで、流れ境界の横のアイコンをクリックすることで、境界サーフェスの表示をオフにします。


    図 40.
  6. Slice Planesツールをクリックします。


    図 41.
  7. モデリングウィンドウで、画面と平行の平面(x-z平面)をクリックします。
  8. スライス平面のマイクロダイアログで、をクリックしてスライス平面を作成します。
  9. 表示プロパティマイクロダイアログで、表示をvelocityに設定し、Legendのトグルスイッチをアクティブにします。
  10. 凡例のUpper Boundを36に変更します。
  11. をクリックして、Colormap nameをRainbow Uniformに設定します。


    図 42.
  12. ガイドバーをクリックしてします。
    シミュレーションの最初の速度コンターは、次に示すようになものになります。


    図 43.
  13. 下部のスライダーを26、51、76、101番目のフレームまでドラッグすることで、0.25s、0.5s、0.75s、1.0sのコンターをプロットします。


    図 44.


    図 45.


    図 46.


    図 47.
  14. PostブラウザSlice Plane 1を右クリックして、Editを選択します。
  15. マイクロダイアログで、Display variableをvelocityからvolume fraction edem particleに変更し、凡例のUpper Boundを0.186に設定します。
  16. ガイドバーをクリックしてします。


    図 48.
  17. 0.5s(ステップ=51/101)、0.75s(ステップ=76/101)、1.0s(ステップ=101/101)のコンターをプロットして、流体領域内のEDEM粒子分布の体積分率を表示します。


    図 49.


    図 50.


    図 51.

EDEMのポスト処理

  1. EDEMのシミュレーションの完了後、左上隅のをクリックして、EDEMのAnalystタブに移動します。
  2. Analyst Treeで、Display > Geometriesを展開してDefault Wallをクリックします。
  3. Display ModeがFilledに設定されていることを確認し、Opacityを0.2に設定します。


    図 52.
  4. Analyst Treeで、Particlesを展開してHeavy particleをクリックします。
  5. 表示色をMagentaに変更します。


    図 53.
  6. Light particleをクリックして、表示色をGreenに変更します。


    図 54.
  7. メニューバーで、以下をクリックすることで時間を0に設定します。


    図 55.
  8. View planeを+Yに設定します。


    図 56.
  9. Viewerウィンドウで、Playback Speedを0.1xに設定し、再生アイコンをクリックして粒子流のアニメーションを再生します。


    図 57.
  10. ドロップダウンメニューをクリックして、別の時間ステップの結果をプロットすることもできます。


    図 58.
    現在のケースでは、0.25s、0.5s、0.75s、0.98sの結果をプロットして、入口と出口の流体領域内の粒子分布の体積分率を表示します。


    図 59.
  11. メニューバーのCreate Graphアイコンをクリックします。
  12. Analyst Treeでをクリックして、プロットタイプをラインに変更します。
  13. X-axisタブをクリックし、下図のとおりに値が設定されていることを確認します。


    図 60.
  14. Y-axisタブをクリックします。下図のとおりに値を設定することで、シミュレーション時にわたる重い粒子の平均滞留時間のプロットを作成します。


    図 61.
  15. をクリックして別のY軸を追加し、TypeをLight particleに設定します。


    図 62.
  16. その他のオプションはすべてそのまま変更せずに、Create Graphをクリックして両方の粒子の平均滞留時間のプロットを作成します。


    図 63.

    緑色のプロットは重い粒子を表し、重い粒子の平均滞留時間の方が軽い粒子(青色)よりも長いことが確認できます。

要約

このチュートリアルでは、基本的なAcuSolve-EDEM双方向(二方向)連成の問題を設定し、実行する方法を知ることができました。最初のパートでは、HyperWorks CFDAcuSolveモデルを設定し、形状をエクスポートしました。次に、HyperWorks CFDによって作成されたEDEM入力ファイルをインポートし、EDEMモデルを設定して、連成シミュレーションを実行しました。連成シミュレーションが完了した後は、HyperWorks CFD PostとEDEM Analystの両方でアニメーションとブロットを作成する方法を学習しました。