Topology optimization of unsaturated flows in multi-material porous media: Application to a simple diaper model

Abstract We present a mathematical approach to optimize the material distribution for fluid transport in unsaturated porous media. Our benchmark problem is a simplified diaper model as an exemplary liquid absorber. Our model has up to three materials with vastly different properties, which in the reference configuration are arranged in parallel layers. Here, swelling is neglected and the geometry of a swollen diaper is assumed and treated as a porous medium of high porosity. The imbibition process is then simulated by solving an unsaturated flow problem based on Richards’ equation. Our aim is to maximize the amount of absorbed liquid by redistributing the materials. To this end, a density based multi-material topology optimization (based on the well known SIMP model) is performed. A sensitivity analysis for the nonlinear transient problem is provided, which enables the application of first order optimization solvers. We perform two- and three-material optimization and discuss several variants of the problem setting. We present designs with up to 45% more absorbed liquid compared to the reference configuration.

[1]  Michael A. Saunders,et al.  SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization , 2002, SIAM J. Optim..

[2]  O. Sigmund,et al.  Topology optimization approaches , 2013, Structural and Multidisciplinary Optimization.

[3]  J. Petersson,et al.  Slope constrained topology optimization , 1998 .

[4]  H. Diersch,et al.  Modeling Unsaturated Flow in Absorbent Swelling Porous Media: Part 2. Numerical Simulation , 2011 .

[5]  Zhenyu Liu,et al.  Topology Optimization Theory for Laminar Flow , 2018 .

[6]  A modification of the mixed form of Richards equation and its application in vertically inhomogeneous soils , 2011 .

[7]  J. Petersson,et al.  Numerical instabilities in topology optimization: A survey on procedures dealing with checkerboards, mesh-dependencies and local minima , 1998 .

[8]  H. Diersch,et al.  Modeling Unsaturated Flow in Absorbent Swelling Porous Media: Part 1. Theory , 2010 .

[9]  E. Rohan,et al.  Optimization of the porous material described by the Biot model , 2018, International Journal of Solids and Structures.

[10]  Ole Sigmund,et al.  Topology optimization for transient wave propagation problems in one dimension , 2008 .

[11]  B. Mohanty,et al.  A new convergence criterion for the modified Picard iteration method to solve the variably saturated flow equation , 1996 .

[12]  Jintu Fan,et al.  Structural optimization of porous media for fast and controlled capillary flows. , 2015, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[14]  S. Hassanizadeh,et al.  The Effects of Swelling and Porosity Change on Capillarity: DEM Coupled with a Pore-Unit Assembly Method , 2016, Transport in Porous Media.

[15]  T. E. Bruns,et al.  Topology optimization of non-linear elastic structures and compliant mechanisms , 2001 .

[16]  Xinmei Hao,et al.  A mass-conservative switching method for simulating saturated–unsaturated flow , 2005 .

[17]  B. Schrefler,et al.  The Finite Element Method in the Static and Dynamic Deformation and Consolidation of Porous Media , 1998 .

[18]  Kok-Kwang Phoon,et al.  Numerical simulation of Richards equation in partially saturated porous media: under-relaxation and mass balance , 2007 .

[19]  Xubin Zeng,et al.  Improving the Numerical Solution of Soil Moisture-Based Richards Equation for Land Models with a Deep or Shallow Water Table , 2009 .

[20]  James K. Guest,et al.  Topology optimization of creeping fluid flows using a Darcy–Stokes finite element , 2006 .

[21]  M. Celia,et al.  A General Mass-Conservative Numerical Solution for the Unsaturated Flow Equation , 1990 .

[22]  Erik Lund,et al.  Material interpolation schemes for unified topology and multi-material optimization , 2011 .

[23]  M. Bendsøe,et al.  Topology Optimization: "Theory, Methods, And Applications" , 2011 .

[24]  S. Torquato,et al.  Composites with extremal thermal expansion coefficients , 1996 .

[25]  Nick Croft,et al.  Computational modelling of variably saturated flow in porous media with complex three‐dimensional geometries , 2006 .

[26]  M. Bendsøe.  Optimal shape design as a material distribution problem , 1989 .

[27]  Ole Sigmund,et al.  On the usefulness of non-gradient approaches in topology optimization , 2011 .

[28]  F. Radu,et al.  A study on iterative methods for solving Richards’ equation , 2015, Computational Geosciences.

[29]  L. A. Richards.  Capillary conduction of liquids through porous mediums , 1931 .

[30]  James K. Guest,et al.  Achieving minimum length scale in topology optimization using nodal design variables and projection functions , 2004 .