Permeability is a parameter that measures the resistance that fluid faces when flowing through a porous medium. Usually, this parameter is determined in routine laboratory tests by applying Darcy’s law. Those tests can be complex and time-demanding, and they do not offer a deep understanding of the material internal microstructure. Currently, with the development of new computational technologies, it is possible to simulate fluid flow experiments in computational labs. Determining permeability with this strategy implies solving a homogenization problem, where the determination of the macro parameter relies on the simulation of a fluid flowing through channels created by connected pores present in the material’s internal microstructure. This is a powerful example of the application of fluid mechanics to solve important industrial problems (e.g., material characterization), in which the students can learn basic concepts of fluid flow while practicing the implementation of computer simulations. In addition, it gives the students a concrete opportunity to work with a problem that associates two different scales. In this work, we present an educational code to compute absolute permeability of heterogeneous materials. The program simulates a Stokes flow in the porous media modeled with periodic boundary conditions using finite elements. Lastly, the permeability of a real sample of sandstone, modeled by microcomputed tomography (micro-CT), is obtained.

This paper formulates and solves a time harmonic inverse scattering problem to reconstruct the effective stiffness distribution of an adhesive bond in a layered elastic plate. The motivation is based on the assumption that localized adhesion flaws that diminish bond stiffness also tend to diminish bond strength. The formulation is based on the invariant imbedding method, applies to isotropic and anisotropic elastic layers, and is essentially that of identifying embedded acoustic sources in elastic layered structures. This paper presents two solutions for the inverse problem: the Born approximation and the exact solution. The example calculations compare the two solutions and show that when imperfections are too large in either magnitude or extent the accuracy of the Born approximation breaks down. The impact of noise and uncertainties in the background properties in the inversion is also investigated. A regularization strategy is introduced in the exact solution that controls solution sensitivity in regions with low signal to noise ratio.

Herein, we describe in detail a methodology to estimate the effective elastic parameters of nodular cast iron, using micro-tomography in conjunction with multiscale finite elements. We discuss the adjustment of the image acquisition parameters, address the issue of the representative-volume choice, and present a brief discussion on image segmentation. In addition, the finite-element computational implementation developed to estimate the effective elastic parameters from segmented microstructural images is described, indicating the corresponding computational costs. We applied the proposed methodology to a nodular cast iron, and estimated the graphite elastic parameters through a comparison between the numerical and experimental results.