Conference Agenda

Multiphase and Porous Media Flows etc.
Tuesday, 28/June/2022:
2:00pm - 3:20pm

Session Chair: Catherine Choquet, La Rochelle University
Location: Michel Crépeau's Lecture Hall, Pôle Communication, La Rochelle University

Pôle Communication Multimédia Reseaux, La Rochelle University, 44 Avenue Albert Einstein, La Rochelle.

2:00pm - 2:20pm

A lattice Boltzmann study of miscible displacement process with precipitation reaction in porous media

Liu, Gaojie1,2; Shao, Ziyu1,2

1School of Energy and Power Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China, People's Republic of; 2Shanghai Key Laboratory of Multiphase Flow and Heat Transfer in Power Engineering, Shanghai 200093, China, People's Republic of

The miscible displacement process with precipitation reaction in porous media widely exists in many fields such as the formation of submarine hydrate, underground storage of CO2, solidification of metal materials, Liesegang phenomenon and so on. In this process, the coupling of flow, diffusion, reaction and the continuous change of solid structure caused by precipitation, makes the overall process very complex. In this work, the miscible displacement process with precipitation reaction is simulated. The time evolutions of fluid concentration field are observed, and the effect of precipitation reaction on displacement process is compared with the one without precipitation reaction. Then, the changes of porosity and displacement efficiency are studied by changing the relevant dimensionless parameters, the following conclusions are obtained: the precipitation reaction will generate crystal hydrate on the surface of the solid framework of porous media, and the continuous expansion of the solid structure will hinder the displacement process; The greater the Da number, the more hydrate will be generated, and the smaller the displacement efficiency will be. However, when the Da number increases to a certain value, the displacement process cannot be carried out.


This study is supported by the National Natural Science Foundation of China (Grant No. 51806142), Shanghai Sailing Program (Grant No. 18YF1418000).

2:20pm - 2:40pm

A lattice Boltzmann study on displacement processes of thermal miscible fluids with density gradient in porous media

Wang, Yongqiang; Liu, Gaojie

School of Energy and Power, University of Shanghai for Science and Technology, China, People's Republic of

Thermal miscible fluid displacements in porous media widely exist in industry process, such as carbon dioxide storage, fuel cell and chemical industry. In thermal miscible displacement process, the displacement front sharp changes because of temperature difference, viscosity difference and, density gradient. It should be noted that the gravity gradient is also one of the important factors affecting displacement efficiency. However, most of researchers have overlooked it.

In this work, a numerical study of the displacement process of thermal miscible fluids with density gradient in a porous medium is carried out. A lattice Boltzmann method (LBM) is used to simulate the displacement process because of its advantages in simulating flow in porous media. The effects of different Rayleigh number (Ra), inclination angle (θ) and viscosity ratio (M) on interface shape, displacement efficiency and the scalar dissipation of the fluids are quantitatively analyzed. The results show that: under the positive inclination angle, the larger the inclination angle, the stronger the fluid instability and the greater the displacement efficiency, the less sufficient the mixing of two fluids; the fluid flow state and the displacement efficiency are not sensitive to the negative inclination angle. With the increase of the viscosity ratio, the formation of the instability phenomenon of "Kelvin–Helmholtz" and "rolling up" will be accelerated, thereby improving the displacement efficiency. The increase of the Rayleigh number increases the fluid instability and the displacement efficiency.


This study is supported by the National Natural Science Foundation of China (Grant No. 51806142), Shanghai Sailing Program (Grant No. 18YF1418000).

2:40pm - 3:00pm

Geochemical digital twin: Mesoscopic modeling of counter-diffusion experiments for resolving carbonate precipitation mechanisms in porous media

Peng, Haonan1,3; Mokos, Athanasios1; Rajyaguru, Ashish1,2; Curti, Enzo1; Grolimund, Daniel2; Churakov, Sergey1,3; Prasianakis, Nikolaos1

1Laboratory for Waste Management, Paul Scherrer Institut, CH-5232 Villigen PSI, Switzerland; 2Swiss Light Source, Paul Scherrer Institut, CH-5232 Villigen PSI, Switzerland; 3Institute of Geological Sciences, University of Bern, CH-3012 Bern, Switzerland

Mechanistic understanding of calcium carbonate precipitation polymorphs (calcite, vaterite, aragonite) is a challenging research problem due to the complexity of the process which is dictated by the interplay of nucleation, crystallization, precipitation kinetics, species transport, homogeneous and heterogeneous chemical reactions. During experimental investigations and in order to decipher the governing mechanisms, it is necessary to obtain in detail the time evolution of the major diffusion pathways and the local conditions at the solid-fluid interfaces during the crystal growth. For this purpose, the digital twin of a laboratory experiment is developed. The D3Q27 multi-component lattice Boltzmann method (LBM) is used as the basis framework for simulating the reactive transport processes in porous media at the pore scale. It enables to obtain information at very high spatial and temporal resolution and to incorporate all necessary physical and chemical processes [1]. The experimental observations at the synchrotron facilities of the Swiss Light Source (SLS PSI), in combination with numerical modeling results will allow us to correlate the type of polymorph precipitated with the saturation state and local composition of the aqueous solution. The preliminary model is tested by simulating the precipitation process of calcium carbonate with different Damköhler numbers which resulted in different morphology of precipitates. For efficient simulations, two levels of acceleration are needed. First, the computationally expensive geochemical solver can be replaced during the simulations by a neural network based metamodel [2]. Second, the code is accelerated even further by using a hybrid CUDA-MPI programming layout which can be executed on several GPUs in parallel, in this case at the Swiss Super Computing Center (CSCS).

[1] Prasianakis, N.I., Curti, E., Kosakowski, G., Poonoosamy, J, Churakov, S.V., Deciphering pore-level precipitation mechanisms, Scientific Reports, 7(1), 13765 (2017)

[2] Prasianakis, N. I., Haller, R., Mahrous, M., Poonoosamy, J., Pfingsten, W., & Churakov, S. V. (2020) Neural network based process coupling and parameter upscaling in reactive transport simulations. Geochimica et Cosmochimica Acta 291, 126-143

3:00pm - 3:20pm

Lattice Boltzmann Model of fluid flow in porous media: tortuosity and porosity effects

D’Orazio, Annunziata; Karaimipour, Arash; Ranjbarzadeh, Ramin

Sapienza University of Rome, Italy

In this work fluid flow through porous media has been modeled by means of the second-order lattice Boltzmann method (LBM). Pore-Scale (PS) method has been used to simulate a porous channel with different porosity degrees by placing square obstacles into the channel. Different arrangements of the obstacles (in-line, staggered and random configurations) differently sized determine porosity and tortuosity factors of porous media.

The velocity-based method has been used to calculate the tortuosity for four configurations of obstacles, different Reynolds numbers from 30 to 300 and a wide range of porosities 0.65, 0.75, 0.85 and 0.95.

Tortuosity values were computed after modeling the velocity field under the effect of mentioned effective parameters. Finally, we found that a possible forward step in this field of study can be varying the configuration of obstacles with the same porosity factor. Pressure drops are evaluated as a function of porosity and tortuosity factor. In addition, with the same obstacle shape and porosity factor, if we change the angle of the obstacles, we can achieve higher tortuosity values.