FFT-based solver for upscaling the thermoporomechanical behavior of infinitely contrasted porous materials

We propose a Fast Fourier Transform based solver for estimating the thermoporomechanical behavior of porous materials, aiming at the periodic homogenization of microstructures with infinitely contrasted properties. Grounded on a thermoporomechanical framework, we first present the complete set of homogenized coupled properties. Subsequently, we isolate each physics, considering an elastoplastic formulation for the mechanical part, and fluid flow and heat conduction being governed by Darcy and Fourier laws before deriving the corresponding homogenized coupled operators. We employ, across the entire numerical framework, the Adaptive Eyre–Milton (AEM) algorithm as an enhanced iterative algorithm that guarantees convergence under infinitely contrasted material properties. Practical implementation aspects and illustrative applications in two- and three-dimensional settings are validated using analytical micromechanics solutions and a FEM-based solver. We emphasize the capability of the proposed framework to compute all the thermohydromechanical (THM) operators within a single tool while maintaining robustness even for materials with infinitely high contrasts in material properties for the three considered physics.