Homogenisation method based on energy conservation and independent of boundary conditions

The foundation of homogenisation methods rests on the postulate of Hill–Mandel, describing energy consistency throughout the transition of scales. The consideration of this principle is therefore crucial in the discipline of Digital Rock Physics which focuses on the upscaling of rock properties. For this reason, numerous studies have developed numerical schemes for porous media to enforce the Hill–Mandel condition to be respected. The most common method is to impose specific boundary conditions, such as periodic ones. However, these boundary conditions influence both the effective property and the size of the REV. The recent study of Thovert and Mourzenko (2020) has shown that most boundary conditions still result in the same intrinsic effective physical property if the averaging is applied outside the range of the boundary layer. From this discovery, it becomes logical to question the status of Hill–Mandel postulate in porous media when homogenising away from the boundary. In this contribution, we simulate Stokes flow through random packings of spheres and a range of rock microstructures. For each, we plot the evolution of the ratio micro- vs macro-scale of the energy of the fluid transport outside the boundary layer, for a growing subsample size of porous media. Here, we prove that we naturally find energy consistency across scales when reaching the size of the Representative Elementary Volume (REV), which is a known condition for rigorous upscaling. Furthermore, we show that this index for the energy consistency is a more accurate indicator of REV convergence since the mean value is already known to be unitary.