The Lattice Boltzmann Method with Deformable Boundary for Colonic Flow Due to Segmental Circular Contractions

We extend the 3D Lattice Boltzmann method with a deformable boundary (LBM-DB) for the computations of the full-volume colonic flow of the Newtonian fluid driven by the peristaltic segmented circular contractions which obey the three-step "intestinal law": (i) deflation, (ii) inflation, and (iii) elastic relaxation. The key point is that the LBM-DB accurately prescribes a curved deforming surface on the regular computational grid through precise and compact Dirichlet velocity schemes, without the need to recover for an adaptive boundary mesh or surface remesh, and without constraint of fluid volume conservation. The population "refill" of "fresh" fluid nodes, including sharp corners, is reformulated with the improved reconstruction algorithms by combining bulk and advanced boundary LBM steps with a local sub-iterative collision update. The efficient parallel LBM-DB simulations in silico then extend the physical experiments performed in vitro on the Dynamic Colon Model (DCM, 2020) to highly occlusive contractile waves. The motility scenarios are modeled both in a cylindrical tube and in a new geometry of "parabolic" transverse shape, which mimics the dynamics of realistic triangular lumen aperture. We examine the role of cross-sectional shape, motility pattern, occlusion scenario, peristaltic wave speed, elasticity effect, kinematic viscosity, inlet/outlet conditions and numerical compressibility on the temporal localization of pressure and velocity oscillations, and especially the ratio of retrograde vs antegrade velocity amplitudes, in relation to the major contractile events. The developed numerical approach could contribute to a better understanding of the intestinal physiology and pathology due to a possibility of its straightforward extension to the non-Newtonian chyme rheology and anatomical geometry.