Slide background

If you are interested in our current studies,
do not hesitate to contact us!

Slide background

A wide range of engineering application

Slide background

Physics based computer graphics

Slide background

Numerical modeling of industrial processes

Slide background

Advanced multi-physics simulations


In my group, discrete element modeling is extensively studied, specifically, modeling of granular and multi-phase flows in a complex shape, modeling of a solid-liquid flow involving a free surface, numerical simulation of a multi-phase flow involving heat transfer and development of the parallel computation technologies. If you are interested in these technologies and would like to comprehend them efficiently, please read my review articles [link1, link2].

1. DEM

離散要素法(Discrete Element Method (DEM))は、ラグランジュ的記述に基づく粉体シミュレーション手法であり、個々の固体粒子の挙動をニュートンの第2法則に基づいて模擬します。 DEMは、バネ、ダッシュポットおよびフリクションスライダーを用いたシンプルな手法ですが、粉体シミュレーションの世界標準の手法であり、極めて有用な情報が得られます。 酒井研究室では、DEMの産業応用を促進するために、DEMの任意形状壁モデルとして符号付距離関数(Signed Distance Function (SDF))を提案しています。 SDFの導入により、既存手法では極めて困難であった、スクリュー搬送[A1]、乾式ミル[A1]、二軸混練器[A2]、リボンミキサー[A3]、複雑形状の粉末金型充填[A4]、ポットブレンダー[A5]、三本ロールミル[A6]などの複雑形状容器内の粉体流動に関する数値シミュレーションがシンプルなアルゴリズムで実行できるようになりました。

The Discrete Element Method (DEM) is well employed in a numerical simulation of a granular flow. The DEM is a Lagrangian approach based on the Newton's second law of motion, where the contact force is modeled by springs, dash-pots and a friction slider. Although the DEM is a very simple model, it provides valuable information for better understanding of complexly granular dynamics. Originally the DEM has hardly been applied to a complex shape domain because arbitrary shape wall modeling is quite difficult. In order to solve this problem, a new arbitrary shaped wall boundary model which is referred to as the Sign Distance Function (SDF) has been developed in my group. Arbitrary shape wall boundary can be created easily by the SDF. The DEM/SDF makes it possible to perform numerical simulation of a granular flow in a complex domain, e.g., screw conveying [A1], a twin-screw kneader [A2], a ribbon mixer [A3], die-filling [A4], a pot blender [A5] and a three-roll mill [A6].

[A1] Y. Shigeto and M. Sakai, Chem. Eng. J., 231, 464-476 (2013)
[A2] M. Sakai et al., Chem. Eng. J., 279, 821-839 (2015)
[A3] G. Basinskas and M. Sakai, Powder Technol., 287, 380-394 (2016)
[A4] Y. Tsunazawa, Y. Shigeto, C. Tokoro, M. Sakai, Chem. Eng. Sci., 138, 791-809 (2015)
[A5] G. Basinskas and M. Sakai, Powder Technol., 301, 815-829 (2016)
[A6] M. Sakai et al., Phys. Fluids, 31, 033302 (2019)

2. DEM-CFD method

固体-流体連成シミュレーションとして、DEMと数値流体力学(Computational Fluid Dynamics (CFD))を連成したDEM-CFD法が広く使われています。 酒井研究室では、DEM-CFD法に関するモデリング&シミュレーションを研究しており、本分野において世界を牽引しています。 DEM-CFD法に壁境界モデルとしてSDFと埋込境界法(Immersed Boundary Method (IBM))を導入した手法(SDF/IBM)[B1-B3]を開発し、 これにより任意壁面形状の固体-流体連成シミュレーションをシンプルな操作で実行できるようになりました。 抗力項の陰解法アルゴリズム[B4]も開発し、これにより従来のDEM-CFD法の安定性を格段に向上させることができました。 Dual Grid法とIBM/SDFを融合した技術の開発により、薄板が含まれる固気混相流のDEM-CFDシミュレーション[B5]が実行できるようになりました。 これらの手法が妥当であることも実験結果との整合により示されています。 また、DEM-CFD法における抗力項の安定計算条件の導出にも成功し、これまで試行錯誤で決めていた時間刻みを理論的に与えることができるようになりました。

The discrete element method (DEM) coupled with computational fluid dynamics (CFD) has been widely employed for a solid-fluid interaction simulation. We are world leading group in the modeling and simulation for the DEM-CFD method, and has developed excellent models thus far. Combination of the SDF and immersed boundary method (IBM) has been developed to create an arbitrary shape wall boundary in the DEM-CFD method [B1-B4]. The SDF/IBM wall boundary model allows us to simulate solid-fluid mixture systems including arbitrary shape walls with simple operation. An implicit algorithm [B5] for the drag force term has been developed, and it significantly improves the stability of the DEM-CFD method. Development of a combined Dual Grid and IBM/SDF technique [B6] makes it possible to perfome the DEM-CFD simulations including thin plates. Adequacy of the above methods has been proved through the several validation tests. My group has also succeeded in deriving the stable calculation conditions for the drag term in the DEM-CFD method. Due to this formula, the time step can be given theoretically, where it was previously determined by trial and error.

[B1] Yao et al., J. Taiwan Inst. Chem. Eng. (2018)
[B2] Mori et al., Powder Technol., 343, 101-112 (2019)
[B3] Li et al., Powder Technol., 377, 89-102 (2021)
[B4] Mori et al., Powder Technol. 372, 48-58 (2020)
[B5] Takabatake & Sakai, Adv. Powder Technol., 31, 1825-1837 (2020)
[B6] Mori & Sakai, Chem. Eng. J., 406, 126841 (2021)

3. Coarse grained DEM

従来のDEMでは、大規模体系のシミュレーションを実質的に実行できないことが課題でした。 酒井研究室で独自開発したDEM粗視化モデルを使用すると、1台のPCであっても莫大な粒子数の体系を計算することができます。 DEM粗視化モデルは、スケーリング則モデルのため、スーパーコンピュータなどの超並列計算の知識がなくても大規模体系の粉体シミュレーションを1台のPCで実行できます。 DEM粗視化モデルは、気流搬送システム[B1]、流動層[B2, B3, B4, B5]などにおいて、実験による検証が行われ、妥当性および信頼性が示されています。 上記の固体-流体連成シミュレーションは、主としてDEMと数値流体力学(Computational Fluid Dynamics (CFD))を連成したDEM-CFD法が使用されます。

Although performing a large-scale DEM simulation on a single PC is always required in industries, it is substantially impossible because of heavy calculation costs. In order to solve this problem, the coarse graining DEM has been proposed in my group. The coarse grain model is one of scaling law models, namely, large-sized computational particle represents group of the original particles. In the coarse grain model, total energy is modeled to agree between the coarse grain particle and the original particles. Hence, the coarse grain model can reduce number of the calculated particles drastically than that in the actual system. The coarse grain model has been applied to pneumatic conveying system [C1] and fluidized beds [C2, C3, C4, C5] so far. Adequacy of the coarse grain model is proven through the verification and validation tests.

[C1] M. Sakai and S. Koshizuka, Chem. Eng. Sci., 64, 533-539 (2009)
[C2] M. Sakai et al., Adv. Powder Technol., 23, 673-681 (2012)
[C3] M. Sakai et al., Int. J. Numer. Meth. Fluids, 64, pp.1319-1335 (2010)
[C4] M. Sakai et al., Chem. Eng. J., 244, 33-43 (2014)
[C5] K. Takabatake et al., Chem. Eng. J., 346, 416-426 (2018)

4. Lagrangian-Lagrangian method

Smoothed Particle Hydrodynamics (SPH)やMoving Particle Semi-implicit (MPS)法は、メッシュフリー法のCFDであり、移流項の数値粘性が考慮されないため、自由液面流れを精度よく計算できると言われています。 酒井研究室では、流体をSPHまたはMPS法で模擬し、固相をDEMで模擬する固液混相流の数値解析手法を独自開発しました。 これらの手法は、DEM-SPH法またはDEM-MPS法と呼ばれています。これらの手法においても、局所体積平均法に基づく流体の支配方程式が採用されています。 DEM-SPH法およびDEM-MPS法は、自由液面を伴う固液混相流を精度よく実行することができます。 DEM-SPH法およびDEM-MPS法は、湿式ボールミル[C1, C2, C3]やビーズミル[C4]において、数値解析結果と実験結果の一致により、妥当な結果が得られることが示されています。

The Smoothed Particle Hydrodynamics (SPH) and the Moving Particle Semi-implicit (MPS) are well known as a mesh-free particle method. These methods have an advantage to simulate a free surface fluid flow precisly because they do not need modeling of the advection term. In my group, the DEM-SPH method and the DEM-MPS method have been developed to simulate a solid-liquid flow involving free surface. Adequacy of these approaches has been illustrated thorough verification and validation tests in wet milling systems [D1, D2, D3, D4].

[D1] X. Sun, M. Sakai, M-T. Sakai, Y. Yamada, Chem. Eng. J., 246, 122-141 (2014)
[D2] M. Sakai et al., Chem. Eng. J., 200-202, 663-672 (2012)
[D3] X. Sun, M. Sakai, Y. Yamada, J. Comput. Phys., 248, 147-176 (2013)
[D4] Y. Yamada and M. Sakai, Powder Technol., 239, 105-114 (2013)

5. DEM-VOF method for a simulation of a gas-solid-liquid flow

自由液面流れの数値シミュレーションについて、SPHやMPS法では実行困難であっても、最先端のVOFを使用すれば容易に実行できることは少なくありません。 酒井研究室では、自由液面流れの世界標準として使用されているVolume-of-Fluid (VOF)を用いて、任意形状壁面内部の固気液三相流れの数値解析技術を独自開発し、DEM-VOF法と名付けました。 このくらいのハイレベルな数値シミュレーションになると、世界トップレベルのグループでもなかなか実行することができません。

Very limited number of groups can develop advanced numerical models for a simulation of a gas-solid-liquid flow, though this is greatly important in science and engineering. Actually, the algorithm is very complex, and the code development is very difficult as well. In my group, the DEM-VOF method [E1] has been developed to simulate the gas-solid-liquid flow in a complex shape domain. In the DEM-VOF, combination of the SDF and the immersed boundary method (IBM) is employed to express the arbitrary shaped wall boundary [E2, E3].

[E1] X. Sun and M. Sakai, Chem. Eng. Sci., 134, 531-548 (2015)
[E2] X. Sun and M. Sakai, Chem. Eng. Sci., 139, 221-240 (2015)
[E3] X. Sun and M. Sakai, J. Chem. Eng. Jpn., 50, 161-169 (2017)

6. MPS method

MPS法は、メッシュフリー法のCFDであり、移流項の数値粘性が考慮されないため、自由液面流れを精度よく計算できると言われています。 MPS法の応用範囲は最先端のVOFに比べて狭いため、産業応用のための研究が必要になります。 酒井研究室では、MPS法を用いて高粘性流体の挙動を効率よく計算するアルゴリズム[B1]を開発しました。 また、MPS法における熱流束境界モデル[B2]の開発も行っています。 本手法は、ガラス固化や炉心溶融の数値シミュレーションに応用されています。

The mesh-free particle method such as the MPS method is well employed in the simulation of a free surface fluid flow. The MPS method has not been established, and hence fundamental studies are still required to improve the applicability. In my group, new numerical models has been developed in the MPS method for efficient calculation of a highly viscous fluid flow [F1] and for setting heat flux at the free surface [F2]. The new methods has been applied in nuclear engineering and resilience engineering.

[F1] X. Sun et al., Nucl. Eng. Des., 248, 14-21 (2012)
[F2] K. Takabatake et al., Int. J. Heat Mass Transf., 103, 635-645 (2016)

7. Parallel computing

Everyone thinks that a calculation should be finished as soon as possible. Recent computational hardware such as CPU and GPU equip multi-cores. The program code should be created corresponding to the latest hardware. In my group, an original multi-thread parallel computation algorithm [G1] has been developed for the DEM simulation, where the OpenMP has been employed for the multi-thread parallel computation.

[G1] Y. Shigeto and M. Sakai, Particuology, 9, 398-405 (2011)