Metamodeling and porous media


Flow through porous media represents an area in fluid dynamics that has wide applications, in both natural and technological fields. To name but a few, it is the process governing the transport of oil, gas and water in oil reservoirs, the flow of non-aqueous phase liquids in contaminated aquifers, or the storage of nuclear or other hazardous wastes. In modern technology, it has been successfully applied in many physical, chemical, thermal and biological processes. This flow problem has been studied empirically for a long time since its formulation in 1856 by Henry Darcy who, first, related the flow rate through the porous bed to the pressure drop across the bed’s sides. Since then, many corrections have the proposed, initially on a pure empirical basis. When the flow through the pores is not too slow, the currently accepted formulation employs a permeability coefficient and a second, inertia-related coefficient, called the Forchheimer coefficient. These two coefficients can be incorporated into a single one, termed the apparent permeability. The possibly anisotropic microstructural features of the porous bed render the apparent permeability a tensorial quantity, in which case the mean flow is not necessarily aligned with the driving pressure gradient across the porous bed.

Recently, scientists at the Institute of Fluid Mechanics of Toulouse, Dr. Nicola Luminari and Professor Christophe Airiau in collaboration with Professor Alessandro Bottaro at the University of Genova, investigated how the direction of the macroscopic pressure gradient, the porosity and the Reynolds number could modify the Darcy and Forchheimer closures arising from a volume-averaged model of a fibrous porous medium. In particular, they focused on a transversely isotropic material composed by parallel fibers of circular cross-section, with one axis of symmetry, a model for example of fluid flow through canopies. Their work is published in the International Journal of Multiphase Flow

In brief, the work by the research team entailed the consideration of a three-dimensional unit cell for the microscopic model, with a generic non-linear forcing whose direction was defined by two Euler angles. Specifically, they first simulated motion in periodic cells, and then solved successive closure problems leading – after applying an intrinsic averaging procedure – to the components of the apparent permeability tensor. Finally, they identified response parameters in a large parametric space by the use of a metamodel based on kriging interpolation.

The authors observed that the medium, composed of bundles of parallel fibers, favored the deviation of the mean flow towards the fibers’ axis when the driving pressure gradient had even a small component along it, and this was seen to be enhanced by a decreasing porosity. In addition, they noted that, for the over one hundred cases which they simulated, the apparent permeability tensor remained, to a very good approximation, diagonal, a fact mainly related to the transversely isotropic nature of the medium. “This is a nice property”. “since it permits to approximate and simulate complex fibrous media by employing a relatively small number of permeability components.” Said Airiau in a statement to Advances in Engineering.

In summary, Luminari, Airiau and Bottaro solved for the microscopic velocity fields through pores in representative volume elements, by varying forcing amplitude and direction, treating over one hundred different cases of flows through arrangements of parallel fibers. This way, they were able to obtain the components of the apparent permeability tensor. The data point computed directly represented the backbone of a very large data base which was subsequently created by the use of a metamodel. “The data base is freely available on an on-line repository and will be useful for validation purposes by other groups”, states Bottaro, adding that “the computed response surfaces should prove invaluable in determining the force caused by the presence of solid, filamentous inclusions in macroscopic simulations of the flow through bundles of fibers whose orientations and dimensions can vary in space and/or time. The possible applications are countless.

Metamodeling and porous media - Advances in Engineering


About the author

Christophe Airiau studied mechanical engineering at Ecole Normale Supérieure de Cachan, near Paris, and aeronautics at ISAE-SupAero, Toulouse. He obtained a PhD in aerodynamics after carrying out research work on the nonlocal stability of boundary layers at ONERA, Toulouse (1994). His “habilitation thesis”, defended in 2004, dealt with receptivity and flow control in shear flows. 

He has been a faculty member at Université Paul Sabatier, Toulouse, since 1995, and full professor since 2008. He is also a member of the Institut de Mécanique des Fluides de Toulouse, where he has led a research group on transition and turbulence.

He has recently published a book on the fundamentals of aerodynamics (in French), for students and engineers, with the goal of transmitting his passion for anything which has to do with aeronautics and aircrafts.

About the author

A graduate of the University of Genova (Laurea in Mechanical Engineering, 1983), Alessandro Bottaro holds a PhD in Mechanical and Aerospace Engineering from Rutgers University, NJ (1988), where he spent four years as a Fulbright Fellow.  He has been Researcher and Lecturer at the Swiss Federal Institute of Technology, Lausanne (1988-1997) and Professor of Mechanics at University Paul Sabatier, Toulouse (1997-2003), before joining the University of Genova as Professor of Fluid Dynamics.

He has been awarded the Edoardo Kramer prize in 2013 from the Istituto Lombardo, Accademia di Scienze e Lettere, as the “Italian engineer with significant contributions to fluid dynamics and turbulence”. In 2014 he was granted a Chaire d’Attractivité from the Institut National Polytechnique de Toulouse, to conduct research in fluid mechanics and biomimetics at the Institut de Mécanique des Fluides de Toulouse for the period 2015-2019.

About the author

Nicola Luminari graduated cum laude in Marine Engineering and Naval Architecture at the University of Genova, with a thesis on the numerical modelling of the fluid behavior around aerodynamic appendices used by small yacht at anchor. He then joined the Institut de Mécanique des Fluides de Toulouse , France, where he obtained a PhD in 2018.  His research was centered on modelling fibrous porous media, meant to be employed as passive flow control devices. During this research period, Nicola Luminari has refined his skills in theoretical and numerical techniques, ranging from hydrodynamic stability analysis, to computational fluid dynamics and meta-modelling. 

Since July 2018, he is employed at Delair, a drone company based in Toulouse, as Head of the Data Science R&D group, developing algorithms for image manipulation and recognition. His group specializes in various computer vision techniques, including machine learning and deep learning.


Nicola Luminari, Christophe Airiau, Alessandro Bottaro. Effects of porosity and inertia on the apparent permeability tensor in fibrous media. International Journal of Multiphase Flow, volume 106 (2018) page 60–74.

Go To International Journal of Multiphase Flow

Check Also

Featured Video Play Icon

A Shocking Development — Researchers Engineer Thinner, Lighter and More Efficient Shock Absorbers