• Search Research Projects
  • Search Researchers
  • How to Use
  1. Back to previous page

Study on Liquid Film Dynamics and Forming and Breaking-up of Foam by using Large-scale Numerical Simulation

Research Project

Project/Area Number 19H00746
Research Category

Grant-in-Aid for Scientific Research (A)

Allocation TypeSingle-year Grants
Section一般
Review Section Medium-sized Section 19:Fluid engineering, thermal engineering, and related fields
Research InstitutionTokyo Institute of Technology

Principal Investigator

青木 尊之  東京工業大学, 学術国際情報センター, 教授 (00184036)

Co-Investigator(Kenkyū-buntansha) 肖 鋒  東京工業大学, 工学院, 教授 (50280912)
白崎 実  横浜国立大学, 大学院環境情報研究院, 准教授 (50302584)
Project Period (FY) 2019-04-01 – 2020-03-31
Project Status Discontinued (Fiscal Year 2019)
Budget Amount *help
¥45,370,000 (Direct Cost: ¥34,900,000、Indirect Cost: ¥10,470,000)
Fiscal Year 2019: ¥17,680,000 (Direct Cost: ¥13,600,000、Indirect Cost: ¥4,080,000)
Keywords液膜 / 非ニュートン性流体 / 気液二相流 / 弱圧縮性流体計算 / AMR / GPU / FENE-CRモデル / 泡沫 / 非ニュートン性気液二相流
Outline of Research at the Start

気液二相流における液膜(液体自由膜)や泡沫は比較的安定であり、表面張力、粘性、非ニュートン性などの物性に強く影響され、液膜内の流動状態が大きく変わる。数値シミュレーションで解明するために、完全陽解法による弱圧縮性流体計算手法を導入し、粘弾性モデルを導入した計算を行う。さらに気液界面に細かい格子を動的に適合させるAMR法を用い、高解像度計算が必要な領域のみを細かい格子で計算する。GPUを大量に用いる液膜・泡沫の超大規模シミュレーションを実現し、マクロなスケールでの液膜や泡沫の生成・崩壊のダイナミクスの本質を解明し、液膜中の流動や熱伝達、物質拡散を明らかにする。

Outline of Annual Research Achievements

開発してきた2次元弱圧縮性気液二相流シミュレーション・コードを液膜・泡沫現象に適用するためにFENE-CRモデルに基づいた粘弾性流体モデルの導入し、気泡上昇に対してベンチマーク計算をよく再現する結果を得た。さらに3次元AMRコードへの拡張とGPU実装を完了させた。FENE-CRモデルの時間発展方程式を分散相側のみで解き、気相側に外挿することで計算の安定性を向上させ,複雑なトポロジー変化を含む流れへの適用性を向上させた。粘弾性モデルを導入した非ニュートン性流体中の気泡上昇計算では、気泡が気液界面に到達して生成された液膜の崩壊が、ニュートン性流体よりも遅れる結果が得られ、粘弾性効果が液膜を維持する効果があることを確かめた。
3次元計算で表面張力が液膜崩壊に与える影響を調べるために、シャボン玉形成シミュレーションを実行した。表面張力項の評価で重要となる曲率の精度を向上させるために、レベルセット関数により勾配の大きさを評価することで、識別変数として用いている低精度のPhase Field変数が液膜内部で使用されることを防ぎ、数値的液膜崩壊の一因を排除した。最も細かい格子幅が0.07825 mmとなる計算(均一格子の場合には4096×2048×2048格子相当)を実行したが、解像度不足が原因して長時間の液膜維持ができないことが明らかになった。そこで、2次元AMRコードにより高解像度化を行い必要な格子解像度を調べた結果、最小格子幅0.0195mmを用いた計算(均一8192×8192格子相当)では初期直径のt = 5.0 secまでの計算で約7.8倍まで液膜を膨張させることができ、非常に薄い液膜を含む流れの計算ができることが明らかになった。
本基盤研究(A)の内容を包含する基盤研究(S)が採択されたため、以降の研究遂行を基盤研究(S)の中で行う。

Research Progress Status

令和元年度が最終年度であるため、記入しない。

Strategy for Future Research Activity

令和元年度が最終年度であるため、記入しない。

Report

(2 results)
  • 2019 Comments on the Screening Results   Annual Research Report
  • Research Products

    (3 results)

All 2019

All Presentation (3 results)

  • [Presentation] AMR法を適用したphase-fieldデンドライト凝固計算の複数GPU並列化2019

    • Author(s)
      坂根 慎治, 高木 知弘, 大野 宗一, 澁田 靖, 青木 尊之
    • Organizer
      第24回計算工学講演会
    • Related Report
      2019 Annual Research Report
  • [Presentation] AMR法を導入した格子ボルツマン法を用いた噴流によるピンポン玉の浮遊シミュレーション2019

    • Author(s)
      渡辺勢也,青木尊之, Yang Kai, 小林宏充
    • Organizer
      第24回計算工学講演会
    • Related Report
      2019 Annual Research Report
  • [Presentation] 界面に適合するAMR法と等温過程に基づく弱圧縮性解法による高解像度気液二相流計算2019

    • Author(s)
      松下真太郎,青木尊之
    • Organizer
      第24回計算工学講演会
    • Related Report
      2019 Annual Research Report

URL: 

Published: 2019-04-18   Modified: 2021-01-27  

Information User Guide FAQ News Terms of Use Attribution of KAKENHI

Powered by NII kakenhi