Effect of macropore tortuosity and morphology on preferential flow through saturated soil: A Lattice Boltzmann study

Unlike micropores where water moves upward or downward based on hydraulic gradient, in macropores, water flows predominantly downward due to the gravity. Therefore, models based on capillary flow are not capable of simulating macropore flow. There are attempts to model the macropore flow using two domains, one for capillary flow and another one for macropores. These models use Richard’s equation for capillary flow and Poiseuille’s law for macropores in which the macropore is approximated to be cylindrical or planar. This study quantifies the magnitudes of the errors induced by this assumption. Influence of macropore shapes and tortuosity was quantified by using a 3D Lattice Boltzmann model, which is capable of simulating fluid flow in micropores as well as macropores of cracked clays. Artificial macropores of constant sectional area and volume, but different shapes were generated in 3D and the influence of macropore shapes, shape related parameters, and tortuosity were systematically investigated. Macropore flow rate decreases with different shapes compared to cylindrical macropores and increase in aspect ratio of sectional shape leads to decrease in macropore flow rate. The maximum effect of bends/turnings along the tortuous macropore was about 25% on overall decrease of flow rate due to tortuosity. However, more detailed study is required on the influence of bends on macropore flow rate. The macropore flow rate reduces by about 70% for tortuosity of 1.41. A prediction equation is verified to predict the flow rate of different shapes and tortuous macropores based on straight cylindrical macropore using aspect ratio and tortuosity factor.