A Novel Finite Element-Based Method for Predicting the Permeability of Heterogeneous and Anisotropic Porous Microstructures

Paris Mulye,Elena Syerko,Christophe Binetruy,Adrien Leygue

DOI: https://doi.org/10.3390/ma17122873

IF: 3.4

2024-06-13

Materials

Abstract:Permeability is a fundamental property of porous media. It quantifies the ease with which a fluid can flow under the effect of a pressure gradient in a network of connected pores. Porous materials can be natural, such as soil and rocks, or synthetic, such as a densified network of fibres or open-cell foams. The measurement of permeability is difficult and time-consuming in heterogeneous and anisotropic porous media; thus, a numerical approach based on the calculation of the tensor components on a 3D image of the material can be very advantageous. For this type of microstructure, it is important to perform calculations on large samples using boundary conditions that do not suppress the transverse flows that occur when flow is forced out of the principal directions. Since these are not necessarily known in complex media, the permeability determination method must not introduce bias by generating non-physical flows. A new finite element-based method proposed in this study allows us to solve very high-dimensional flow problems while limiting the biases associated with boundary conditions and the small size of the numerical samples addressed. This method includes a new boundary condition, full permeability tensor identification based on the multiscale homogenization approach, and an optimized solver to handle flow problems with a large number of degrees of freedom. The method is first validated against academic test cases and against the results of a recent permeability benchmark exercise. The results underline the suitability of the proposed approach for heterogeneous and anisotropic microstructures.

materials science, multidisciplinary,metallurgy & metallurgical engineering,physics, applied,chemistry, physical, condensed matter

What problem does this paper attempt to address?

The paper mainly aims to address the following issues:

### Research Background and Objectives

- **Research Background**: Permeability is a fundamental property of porous media, describing the ability of fluids to pass through a porous network under a pressure gradient. For heterogeneous and anisotropic porous materials, measuring permeability is both difficult and time-consuming. Therefore, numerical methods for calculating permeability tensor components based on 3D images have significant advantages.

- **Main Objectives**: Propose a new finite element method to predict the permeability of heterogeneous and anisotropic porous microstructures. This method can handle high-dimensional flow problems and limit biases caused by boundary conditions and numerical sample sizes. Additionally, the method includes a new boundary condition, full permeability tensor identification based on multi-scale homogenization, and an optimized solver for handling large degrees of freedom flow problems.

### Specific Problems Addressed

1. **Developing a New Method**: To accurately predict the permeability of heterogeneous and anisotropic porous materials, a new method capable of solving high-dimensional flow problems is needed, while avoiding non-physical flows introduced by improper boundary conditions.

2. **Improving Computational Efficiency**: For multi-dimensional problems with complex geometries, effective methods are needed to handle surfaces and reduce computational costs. The finite element method excels in this regard as it can naturally handle surfaces without requiring additional boundary conditions.

3. **Determining the Full Permeability Tensor**: To capture the effects of local anisotropy, it is necessary to compute all components of the permeability tensor, including off-diagonal components.

4. **Handling Large Sample Calculations**: To limit biases related to the representativeness of the medium, permeability calculations should be conducted on a sufficiently large representative volume. This requires the efficiency of numerical problem-solving and the amount of memory needed.

5. **Optimizing Boundary Conditions**: The boundary conditions used should not suppress transverse flows, which is particularly important for heterogeneous materials. It is necessary to develop boundary conditions that do not produce non-physical flows.

In summary, this paper aims to propose a new finite element method to overcome the main challenges encountered by existing methods in predicting the permeability of heterogeneous and anisotropic porous media, such as low computational efficiency and inappropriate boundary conditions.