Japan Geoscience Union Meeting 2025

Presentation information

[E] Oral

M (Multidisciplinary and Interdisciplinary) » M-IS Intersection

[M-IS01] Particulate Gravity Current

Fri. May 30, 2025 3:30 PM - 5:00 PM Exhibition Hall Special Setting (6) (Exhibition Hall 7&8, Makuhari Messe)

convener:Hajime Naruse(Department of Geology and Mineralogy, Graduate School of Science, Kyoto University), Yuichi Sakai(Faculty of Agriculture, Utsunomiya University), Hiroyuki A. Shimizu(Sabo and Landslide Technical Center), Takahiro Tanabe(National Research Institute for Earth Science and Disaster Resilience), Chairperson:Yuichi Sakai(Faculty of Agriculture, Utsunomiya University), Hiroyuki A. Shimizu(Sabo and Landslide Technical Center), Takahiro Tanabe(National Research Institute for Earth Science and Disaster Resilience)

4:00 PM - 4:15 PM

[MIS01-02] Numerical simulation of lava flow with plastic skin using a particle-based method

*Shingo Tomita1, Takuma Sato1, Satoshi Murakami1, Joe Yoshikawa3, Makoto Sugimoto1, Hisaya Komen2, Masaya Shigeta1 (1.Tohoku University, 2.Joining and Welding Research Institute, Osaka University, 3.Industrial Technology Institute Miyagi Prefectural Government)


Keywords:Pahoehoe lava, Plastic skin, Non-Newtonian fluid, Solid-liquid phase change, Numerical simulation, Smoothed Particle Hydrodynamics

Characteristic phenomenon in lava flows, the formation of thin plastic skin by surface cooling, was numerically simulated successfully using the Smoothed Particle Hydrodynamics (SPH) method. This specific phenomenon may affect the extent where lava flows. Therefore, it is indispensable to elucidate and incorporate the phenomenon into models to accurately predict areas which lava strikes.
In order to simulate the formation of thin plastic skin by surface cooling, the model considers for temperature-dependent viscosity. The skin was treated as a low-temperature and high-viscosity region. In addition, an apparent viscosity model was introduced to simulate Bingham fluid behavior, in which apparent viscosity changes depended on strain rates. Using this model, the process of lava flowing on an inclined plate with formation of surface skin was simulated. Due to computational constraints, the viscosity coefficient was adjusted to 1/100 of its original value, with upper and lower limits set in the simulation. The following outcomes were observed in the simulations:
1. Rapid deceleration was observed by the introduction of the apparent viscosity model. In addition, the process was identified in which the apparent viscosity increased as the flow velocity decreased, and then the flow velocity further decreased. Finally, the process led the flow to stop. This process replicated the behavior of Bingham fluids qualitatively. Comparisons with Newtonian fluids confirmed that Bingham behavior significantly influenced lava flow velocity and distance.
2. The model successfully simulated the formation of thin skin (low-temperature and high-viscosity regions) observed on the surface of actual lava flow, by considering the temperature-dependent viscosity. Comparison with simulations without skin formations revealed that the skin suppressed the flow velocity. This is a qualitative replication of the flow suppression effect of the skin in actual lava flow. The simulation results confirmed that the formation of the skin affects the velocity and distance of lava flow. The suppression effect of the skin was also confirmed to influence the appearance of the Bingham behavior described above.
In summary, this model successfully simulated the non-Newtonian behavior of lava as a Bingham fluid, as well as the flow suppression effect of surface skin. The results indicated that Bingham behavior and the suppression effect of the skin have a significant impact on the flow distance. These findings emphasize the necessity of consideration of both non-Newtonian behavior and surface skin formation for accurate predictions of lava flow behavior.