The focus of our research is the migration of fluids in porous media. Results of this work have included the following:


1

Development of a numerical simulation for gas-water, non-aqueous-phase-liquid (NAPL), three-phase flow with a dusty gas model.

We developed, in particular, a method for numerically simulating gas-water-NAPL, three-phase flow and the migration of a chemical substance in a porous medium. The gas-water-NAPL, three-phase flow simulation employs a modified Picard iteration and finite element methodology. We have applied the dusty gas model in the case of multiple components to the migration of a chemical substance in the gas phase of a porous medium. Finally, in our laboratory we have integrated the dusty gas model with the gas-water-NAPL, three-phase flow simulation and developed the following simulation program (anwflowDGM written in fortran):
Simulator name anwflowDGM
sourceSave SourceCode(.zip)
Example of input data　Manual (in Japanese)

References

  • Hibi, Y., Nakata, S., Sugiyama, A., 2010. Estimation of effects on NAPL residual saturation for air-water-NAPL phase flow in porous media. Journal of Geotechnical Engineering Japan Society and Civil Engineers. 66(2), pp. 418 E29 (in Japanese).
  • Hibi, Y., Fujinawa, K., Nishizaki, S., Okamura K., Tasaki, M., 2009. Multi-component migration in the gas phase of soil:　Comparison between results of experiments and simulation by Dusty Gas Model. Soils and Foundations 49(4), 569 E82.
  • Hibi, Y., 2008. Formulation of a dusty gas model for multi-component diffusion in the phase of soil. Soils and Foundations 48(3), 419 E32.
  • Hibi, Y., Fujinawa, K., 2005. A comparison of numerical model for air-water-NAPL three-phase flow in porous media. Journal of Geotechnical Engineering Japan Society and Civil Engineers. 797/VII, pp. 81 E4 (in Japanese).
  • Hibi, Y., Fujinawa, K., Fujiwara, Y., 2001. A comparison of finite element solution for pressure based and mixed type equations of two-phase flow in porous media, In International Association for Research Committee on Groundwater Hydraulics 1st Groundwater Seminar between China, Korea and Japan International Association for Hydraulic Research, Fukuoka, pp. 41 E2.
2

Development of a method for obtaining Knudsen coefficients and diffusion coefficients with tortuosity and mechanical dispersion.

Determination of the Knudsen coefficient and diffusion coefficient with tortuosity and mechanical dispersion requires simulation with the dusty gas model. We developed a method for obtaining these parameters by using column experiments and an inverse simulation program as follows:
Inverse simulator name invgasflow
sourceSave SourceCode(.zip)
Example of input data

References

  • Hibi, Y., Kanou, Y., Ohira, Y., 2012. Estimation of mechanical dispersion and dispersivity in a soil-gas system by column experiments and dusty gas model. Journal of Contaminant Hydrology 131, pp. 39 E3.
  • Kanou, Y., Hibi, Y., Ohira, Y., 2012. Valuation of parameters for migration of components in the gas system obtained by using column experiments, Journal of Geotechnical Engineering Japan Society and Civil Engineers 68(1), pp. 57 E7 (in Japanese).
  • Hibi, Y., Taguchi, T., 2011. Development of a method to estimate the dispersion and Knudsen diffusion coefficients in soil by using the dusty gas model. Journal of Geotechnical Engineering Japan Society and Civil Engineers 67(2), 198 E04 (in Japanese).
3

Atmosphere-Surface water-Groundwater Multiphase flow numerical simulation method (ASGMF numerical simulation method)

ASGMF numerical simulation method couplies fluid flow in a surface system consisting of atmosphere and surface water and gas-groundwater two phase flow in soil (a porous medium) and is able to simulate a interaction between the fluid in surface system and the gas-groundwater two phase flow in soil (a porous medium). This numerical simulation employs a one-field model for multiple immiscible fluids to simulate the surface system and takes into account effects of circulation in atmosphere and surface water. The gas-groundwater two phase flow in soil (a porous medium) is modeled with a water-saturation formulation including total velocity obtained from global pressure.
　　　　　　　　　　　　

References of Atmosphere-Surface water-Groundwater Multiphase flow numerical simulation method (ASGMF numerical simulation method)

　　　　　　　　　　　
  • Hibi, Y., Tomigashi, A., 2015. Evaluation of a coupled model for numerical simulation of a multiphase flow system in a porous medium and a surface fluid. Journal of Contaminant Hydrology 180, 34 - 55.
  • Hibi, Y., Tomigashi, A., Hirose, M., 2015. Evaluation of a numerical simulation model for a system coupling atmospheric gas surface water and unsaturated or saturated porous medium. Journal of Contaminant Hydrology 183, 121-134
  • 日比義彦，冨樫聡：多孔質体内の流れの飽和度型支配方程式を用いた堤体の越流数値解析，混相流，Vol.29 No.4(2015)，pp.326-334，2015.
  • 日比義彦，冨樫聡：複雑な形状の堤体の越流問題への大気−表面水−多孔質体連成数値解析手法の適応性の検討，混相流，Vol.31 No.1，pp.29-36，2017.
  • Hibi, Y., Tomigashi, A., 2018. A numerical simulation model for a coupled porous medium and surface fluid system with multiphase flow. Journal of Groundwater Hydrology 60(4), 409-434.
　　　　　　　　　　
4

大気-表面水-地下水連成解析プログラムの講習会のお知らせ

11月15日14時から17時まで地盤工学会地下会議室にて大気-表面水-地下水連成解析プログラムの講習を無料で開催します。講習会では解析プログラムのソースを無料で配布し、解析プログラムについて詳細に説明します。参加希望の方はhibiy@meijo-u.ac.jpまでご連絡ください。２日前まで受け付けます。
We lecture about Atmosphere-Surface water-Groundwater Multiphase flow numerical simulation method (ASGMF numerical simulation method) and the source code wrote with FORTRAN at 2 to 5 p.m. on the 15th November, at a basement lecture room in Japanese geotechnical society office (4-38-2 Sengoku , Bunkyo-ku ,Tokyo ,112-0011 , Japan). The lecture fee is free, and the source code is distributed in the lecture. If you will participate in the lecture, you sent us e-meal (hibiy@meijo-u.ac.jp) until the 13th November. By the way, the method and the source code in the lecture will be explained in Japanese.