ACU-T：3311 代数オイラーモデルを使用した多相核沸騰

前提条件

このチュートリアルでは、管内の二相核沸騰の過渡シミュレーションを、代数オイラーモデルを使用して実行する手順を説明します。このチュートリアルを開始する前に、HyperWorks 入門チュートリアルである ACU-T：1000 HyperWorksユーザーインターフェースをすでに完了し、HyperMesh、AcuSolve、およびHyperViewの基本を理解しているものとします。このシミュレーションを実行するには、ライセンス供与済みバージョンのHyperMeshおよびAcuSolveにアクセスできる必要があります。

このチュートリアルを実行する前に、ここをクリックしてチュートリアルモデルをダウンロードしてください。 ACU-T3311_Steiner.hm をHyperMesh_tutorial_inputs.zipから抽出します。

HyperMeshデータベース（.hmファイル）には、メッシュ済みのジオメトリが含まれているため、このチュートリアルには、ジオメトリのインポートとメッシュ生成に関する手順は含まれません。

問題の説明

このチュートリアルで扱う問題は、図 1で図式的に示されています。ここでは例としてスタイナー問題を取り上げて、AcuSolveにおける多相核沸騰モデリングの機能を確認します。相変化を伴う代数オイラー（AE）モデルは、分散領域と搬送領域の間で見られる熱伝達と運動量交換のシミュレーションに使用します。

この問題では、搬送領域の材料は水であり、分散領域の材料は蒸気であると考えます。圧力が2バールで温度が95 ℃の水が、0.39m/secの平均速度で入口に流入し、130℃に保たれている加熱壁を通過します。


図 1. チャネルの概略図

HyperMeshモデルデータベースを開く

  1. HyperMeshを起動し、AcuSolveのユーザープロファイルを読み込みます。
    User ProfilesからAcuSolveを選択する方法については、HyperMeshの入門チュートリアルACU-T：1000 HyperWorksユーザーインターフェースをご参照ください。
  2. 標準ツールバーのOpen Modelアイコン をクリックします。
    Open Modelダイアログが開きます。
  3. モデルファイルの保存先ディレクトリを参照します。HyperMeshファイルACU-T3311_Steiner.hmを選択してOpenをクリックします。
  4. File > Save Asをクリックします。
    Save Model Asダイアログが開きます。
  5. 名前をMultiphase_NBとして新しいディレクトリを作成し、このディレクトリへ移動します。
    このディレクトリが作業ディレクトリになり、シミュレーションに関連するすべてのファイルがこの場所に保存されます。
  6. データベースのファイル名としてSteinerと入力するか、都合のいい名前を選択して入力します。
  7. 保存をクリックしてデータベースを作成します。

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

解析パラメータの設定

  1. Solverブラウザで01.Globalを展開してPROBLEM_DESCRIPTIONをクリックします。
  2. エンティティエディターで、Analysis typeをTransientに変更します。
  3. Abs. temperature offsetを200000に設定します。
  4. Turbulence modelをSpalart Allmarasに設定します。
  5. Temperature equationをAdvective Diffusiveに設定します。
  6. Multiphase equationをAlgebraic Eulerianに設定します。
  7. Global Gravity in the Z directionを-9.8に設定します。


    図 2.

ソルバー設定

  1. Solverブラウザで01.Globalの下の02.SOLVER_SETTINGSをクリックします。
  2. エンティティエディターで、Max time stepsを40に設定します。
  3. Initial time incrementを0.01に設定します。
  4. Min stagger iterationsを2に設定します。
  5. Max stagger iterationsを4に設定します。
  6. Relaxation factorを0に設定します。
  7. Flow、Turbulence、およびFieldの設定がOnになっていることを確認します。


    図 3.

節点出力の定義

  1. Solverブラウザで17.Ouputを展開してNODAL_OUTPUTをクリックします。
  2. エンティティエディターで、Time step frequencyを10に設定します。
  3. Output initial condition欄をオンにします。


    図 4.

材料モデルパラメータの設定

この手順では、まずVaporという名前の新しい材料を作成します。次に、多相材料モデルを設定します。

  1. Solverブラウザで、02.Materialsを展開します。
  2. FLUIDを右クリックし、Createを選択して新しい流体材料を作成します。これにVaporという名前を付けます。
  3. エンティティエディターで、Densityの値を1.128に設定します。
  4. Specific heatの値を2174.59に設定します。
  5. Latent heat typeをConstantに変更します。
  6. Specific heatの値を2256000に設定します。
  7. Latent heat temperatureの値を393.45に設定します。
  8. Viscosityの値を1.29e-05に設定します。
  9. Conductivityの値を0.02749に設定します。


    図 5.
  10. SolverブラウザのFLUIDの下のWater_HMをクリックし、材料プロパティを変更します。
    1. Densityを956.36kに変更します。
    2. Specific heatを4204.8に変更します。
    3. Viscosityを0.000314に変更します。
    4. Conductivityを0.675に変更します。
  11. Solverブラウザで、02.MaterialsのMULTIPHASEを右クリックし、Createを選択して新しい多相材料を作成します。VapourWater-Eulerian_HMという名前を与えます。
  12. エンティティエディターで、Field Interaction TypeがAlgebraic Eulerianに設定されていることを確認します。
  13. Carrier MaterialをWater_HMに設定します（まだ設定されていない場合）。
    この解析では水が搬送領域なので、このオプションは変更せずに、このままとします。
  14. Number of Dispersed Fieldが1に設定されていることを確認します。


    図 6.
  15. Edit Field Interaction Models横のチェックボックスを有効にします。
    Field Interaction Modelダイアログが開きます。
  16. このダイアログで、Dispersed Field 1 MaterialがVaporに設定します。
  17. Enable Phase change のチェックボックスを有効にします。
  18. Phase change typeをWall boilingに変更します（設定されていない場合）。
  19. Dispersed field diameter typeをBoiling Kurul exponentialに設定します。


    図 7.
  20. Save & Closeをクリックします。

境界条件と節点初期状態の設定

境界条件と節点初期状態の設定

この手順では、多相流体ボリュームに材料特性を割り当て、つづいてサーフェス境界条件を割り当てます。

  1. Solverブラウザで、11.Volumes > MULTIPHASEを展開します。
  2. FluidVolumeをクリックします。エンティティエディターで以下を設定します。
    1. TypeをMULITPHASEに変更します（設定されていない場合）。
    2. MaterialとしてVaperWater_Eulerian_HMを設定します。


    図 8.
  3. Solverブラウザで、12.Surfaces > INFLOWを展開します。
  4. Inflowをクリックします。エンティティエディターで以下を設定します。
    1. TypeをINFLOWに変更します（設定されていない場合）。
    2. Inflow TypeをAverage velocityに設定します。
    3. Average velocityを0.39m/secに設定します。
    4. Temperatureを368.15Kに設定します。


      図 9.
    5. No. of multiphase incoming fieldsを2に設定して、Enterを押します。
    6. ダイアログで、最初のIncoming FieldにWater_CF.VaporWater_Eulerian_HMを選択し、そのVolume Fractionを1に設定します。
    7. 2つ目のIncoming FieldでVapor_DF1. VaporWater_Eulerian_HMを選択し、そのVolume Fractionsを0に設定します。


      図 10.
    8. Save and Closeをクリックします。
  5. 同様にして、OUTFLOWを拡張し、Outflowコンポーネントをクリックします。エンティティエディターで、TypeをOUTFLOWに変更します（設定されていない場合）。


    図 11.
  6. 同様にして、WALLを拡張し、Side_MaxYコンポーネントをクリックします。エンティティエディターで、TypeがWALLに設定されていることを確認します。


    図 12.
  7. 同様にして、Side_MaxY、Top、BottomコンポーネントのTypeがWALLに設定されているのを確認します。
  8. HeatedWallをクリックします。エンティティエディターで以下を設定します。
    1. TypeがWALLに設定されていることを確認します。
    2. Temperature BC typeをValueに設定します。
    3. Temperatureの値を403.15に設定します。


    図 13.
  9. モデルを保存します。

節点初期状態の定義

  1. Solverブラウザで01.Globalを展開して03.NODAL_INTIAL_CONDITIONをクリックします。
  2. エンティティエディターで、x-Velocityの値を0.39に設定します。
  3. Temperatureの初期値を368.15に設定します。
  4. FieldタブでNo. of Fieldsを2に設定してEnterキーを押します。
  5. Fieldsダイアログで、値を以下のように設定します。


    図 14.
  6. Save & Closeをクリックします。
  7. モデルを保存します。

解析計算

この手順では、HyperMeshからAcuSolveを直接起動して解析を完了します。

  1. すべてのメッシュコンポーネントの表示をオンにします。
    解析を実行するには、アクティブなすべてのコンポーネントのメッシュを可視化した状態にする必要があります。
  2. ACUツールバーの をクリックします。
    Solver job Launcherダイアログが開きます。
  3. オプション: 解析時間を短縮するには、使用可能なプロセッサの数に応じて、使用するプロセッサの数に大きい値（4または8）を設定します。
  4. Output time stepsがAllに設定されているのを確認します。
  5. 他のオプションはデフォルト設定のままにし、Launchをクリックして解析プロセスを開始します。


    図 15.

HyperViewによる結果のポスト処理

AcuSolveの実行が終了したら、HyperWorks Solver View ダイアログを閉じます。HyperMesh Desktopウィンドウで、 AcuSolve ControlおよびSolver job Launcherダイアログを閉じます。次のステップでは、蒸気の体積分率のコンターをプロットします。

HyperViewインターフェースへの切り替えとAcuSolveモデルと結果の読み込み

  1. HyperMesh Desktopウィンドウで、グラフィックスウィンドウの左下隅のClientSelectorドロップダウンをクリックします。


    図 16.
  2. リストからHyperViewを選択します。
  3. 表示されたポップアップダイアログで、Yesをクリックします。
    インターフェースがHyperViewに変更されます。

    HyperViewを読み込むと、デフォルトでLoad model and resultsパネルが開きます。このパネルが表示されない場合は、File > Open > Modelの順にクリックします。

  4. Load model and resultsパネルで、Load modelの隣にある をクリックします。
  5. Load Model Fileダイアログで、作業ディレクトリに移動して、ポスト処理する解析実行のAcuSolve .Logファイルを選択します。この例で選択するファイルは、Steiner.1.Logです。
  6. Openをクリックします。
  7. パネル領域でApplyをクリックしてモデルと結果を読み込みます。
    読み込むと、モデルが形状で色分けされます。

蒸気の体積分率に関するコンターの作成

この手順では、蒸気の体積分率に関するコンターを表示します。
  1. これは過渡ケースであるため、最後の時間ステップにおける結果をプロットする必要があります。このためには、ResultsブラウザでTimeドロップダウンメニューをクリックし、リスト内の最後のオプションを選択します。


    図 17.
  2. ResultsブラウザでのComponentsリストを拡張表示します。
  3. Isolate Shownアイコンをクリックし、Ctrlキーを押しながらHeatedWallとBottomコンポーネントを選択して、必要なコンポーネントを除くすべてのコンポーネントの表示をオフにします。


    図 18.
  4. ResultsツールバーでをクリックしてContourパネルを開きます。
  5. パネル領域で、Result typeをVolume_fraction-1-Vapor_DF1.VaporWater_Eulerian_HM(s)に変更します。
  6. Componentsエンティティセレクターをクリックします。Extended Entity Selectionダイアログで、Displayedを選択します。
  7. Applyをクリックして、蒸気の体積分率のコンターをプロットします。
  8. パネル領域のDisplayタブで、Discrete colorオプションをオフにします。


    図 19.
  9. Legendタブをクリックし、Edit Legendをクリックします。表示されたダイアログで、Numeric formatをFixedに変更してOKをクリックします。
    次の図のようなコンタープロットが表示されます。


    図 20.

要約

このチュートリアルでは、代数オイラー多相モデルを使用し、過渡二相核沸騰流の問題を設定して解析するための基本的なワークフローを体験しました。まず、HyperMeshでモデルをインポートしました。次に、シミュレーションパラメータを設定して、HyperMeshから直接AcuSolveを起動しました。AcuSolveによる解析が完了した後、HyperViewを使用して結果をポスト処理し、体積分率のコンタープロットを作成しました。