Mesh generation for periodic 3D microstructure models and computation of effective properties
Graphical abstract
Introduction
Porous media can be found in many natural as well as artificial physical, biological and chemical systems. From the composition of soils [1], [2], through which liquids seep into the ground water, to the mechanical stiffness of cements [3], [4], from battery electrodes [5], [6], [7], in which lithium ions are stored, to sponge-based filtration materials [8]: the porous microstructure of the respective system has a crucial impact on the overall behavior [9]. For example, the morphology of electrodes in lithium-ion batteries significantly influences the electrochemical properties [10], [11], [12], [13], which is the main reason why tailored structuring and manufacturing of anodes and cathodes is one promising approach to improve the performance of the cell [14], [15], [16]. Thus, it is a major issue in many research areas to design the microstructure in such a way that the overall performance, e.g. permeability, electrical conductivity, mechanical stiffness, energy density and further quantities, is optimized.
From a mathematical point of view, the impact of the 3D morphology of porous media on their macroscopic behavior, e.g., the flow rate of water through soil or the flux of lithium ions through a battery electrode, can be studied with homogenization techniques. A prominent and mathematically sound tool is periodic homogenization theory [17], which assumes that the porous medium, given as a certain domain Ω, is a periodic repetition of some representative volume element ω, see Fig. 1. This method allows to derive a set of partial differential equations (PDEs) for which the porous microstructure is not spatially resolved anymore. This significantly reduces the numerical complexity of the problem. The method is based on an asymptotic expansion of the balance equation in terms of ε, which is the ratio between a macro-scale length L and the cell-scale length ℓ, i.e., . In the asymptotic limit, where , a set of homogenized balance equations is then obtained, together with some porous media parameters.
Consider the decomposition , where the set is simply connected and corresponds exemplarily to an electrolyte phase, and is multiply connected, denoting exemplarily a solid phase. The interface between and is denoted by . As already mentioned above, Ω is a periodic repetition of the unit cell , and the common interface .
For a scalar balance equation with surface reactions, we have the general model
This problem is formally solved via the introduction of a multi-scale expansion , where , which yields a sequence of PDEs to determine the unknown functions in the orders of the scaling parameter ε. Briefly summarized, with periodic homogenization [17], [18], [19] one obtains the following statements:
Order yields essentially , i.e., the leading order function is independent of the micro-scale y.
Order yields as well as the condition for the (geometrical) corrector function . The typeface emphasizes that is a (numerical) solution of the cell problem CP 1, which depends thus only on the micro-scale y. Since and thus also depend only on the macro-scale x, i.e. , the result for the first successive term is a separation between microscopic geometrical effects and macroscopic gradients of the leading order term . This is a central feature of (periodic) homogenization theory and we show in Section 4 how is related to the tortuosity of a microstructure.
Order yields the PDE for the leading order term and the leading order flux , where the porous media parameters are given by
- 1.
the porosity (or phase fraction) of ,
- 2.
the interfacial area of ,
- 3.
and the (flux) corrector,
For the Stokes problem in a similar manner we obtain where periodic homogenization leads to the Darcy flow The corrector is frequently called a permeability tensor, where and determined from the cell problem [20]
For every PDE problem, e.g., PDE1 or PDE2 described above, periodic homogenization leads to a different cell problem, i.e., CP1 or CP2, which has to be solved in order to determine the effective porous media parameters. However, all of these cell problems do have in common that some stationary PDE system has to be solved on the periodic representative volume element ω. Since this is analytically possible only for a very tiny amount of geometries, the cell problems have in general to be solved numerically. And, in order to so, adequate discretizations of ω are required.
Various approaches for the discretization of the representative volume element ω exist and we briefly review exclusively those which ensure the periodicity of ω. The most simple approach is a voxel based discretization of ω, with equal edge length of the voxel in all three dimensions. This format of discretization is widely used in 3D imaging. However for the purpose of numerical calculations voxel based meshes are inappropriate, (i) because the computational degrees of freedom scale with for N equally sized voxels [3], or (ii) because local refinements (or coarsening), leading to , lead to hanging nodes,1 which are numerically very problematic. To get rid of these problems a mesh for numerical calculations in 3D is typically built by tetrahedra, with which a 3D geometry can be discretized and locally refined without producing hanging nodes. This yields far more efficient numerical calculations and the corresponding meshes are said to be of high quality. Numerical simulations of 3D microstructures allow then for virtual materials testing [21], [22].
The generation of volumetric meshes for realistic microstructures is an interdisciplinary topic and various approaches are found throughout the literature [23], [24], [25]. A central aspect is the assumption regarding the geometrical shape of the particles or inclusions in the microstructure. For example a rather general approach of modeling inclusions in matrix materials assumes ellipsoidal shapes [26], which is explicitly exploited in the periodic mesh generation. However, tomographic imaging methods have shown that realistic 3D microstructures of various functional materials are significantly more complex. A flexible tool to model and simulate the morphology of such particle systems is hence desirable and we propose a complete pipeline for this issue.
In this paper, we propose a robust mesh generation for periodic representative volume elements of realistic microstructures, see Fig. 2, with geometrically more flexible star-shaped particles.
The method is based on a description of the microstructure in terms of spherical harmonics, a subsequent surface mesh generation of and , and finally a volume mesh generation based on TetGen [27]. The proposed method can be applied to a broad spectrum of scenarios arising in different fields of research since numerous scientific problems involve solving a system of differential equations on periodic porous media. Another advantage of the presented approach is that periodic boundary conditions can be easily applied in x-, y- and z-direction as well as to an arbitrary subset of directions. This can be used for example in battery research, where the size of electrodes is typically several orders of magnitudes larger in in-plane direction compared to the thickness of the electrode such that it is reasonable to consider periodic boundary conditions in two directions.
The rest of this paper is organized as follows. In Section 2, we describe the generation of periodic 3D microstructures based on spherical harmonics and a stochastic microstructure model. Then, in Section 3, the generation of a quality volume mesh on the basis of the representation of the particle system via spherical harmonics is explained. In Section 4, some numerical examples are presented. Finally, in Section 5, the paper is concluded by a summary of the main results and an outlook to possible future research is given.
Section snippets
Generation of periodic porous 3D microstructures
To generate a periodic representative volume element ω, we use the stochastic microstructure modeling approach described in [21], which basically consists of three steps. First, a non-overlapping sphere packing is generated, where the volume fraction of the solid phase as well as the particle size distribution can be preset. More precisely, the radii of spheres are drawn from the predefined particle size distribution until the target volume fraction is achieved. Initially, the
Periodic mesh generation
We now discuss the mesh generation for the microstructures described in Section 2, yielding discrete approximations for ω and . Recall that the microstructure to be meshed has the representation , where N is the number of particles. The surface of each particle is parametrized by with being the midpoint and the radius function of . As explained in the previous section, is expanded in terms
Numerical results
In this section we discuss three examples of microstructures to provide an impression of the proposed method. The first example considers simple morphologies of equally sized spherical particles, where the porosity is varied and the resulting porous media parameters are compared to the Bruggeman approximation, see Section 4.1. The second example, considered in Section 4.2, deals with a porous microstructure consisting of N particles, where the particle radii are polydisperse. Such
Conclusion and outlook
In this paper, we presented a novel, robust method for generating high-quality volume meshes based on the spherical harmonics representation of particulate microstructures. For this purpose, a stochastic 3D microstructure model has been used in order to generate virtual, but realistic two-phase microstructures as structural input for the mesh generation procedure. In addition, the presented method is able to handle periodic boundary conditions in a predefined set of directions as well. After
CRediT authorship contribution statement
Manuel Landstorfer: Conceptualization, Funding acquisition, Methodology, Software, Visualization, Writing – original draft. Benedikt Prifling: Methodology, Software, Visualization, Writing – original draft. Volker Schmidt: Funding acquisition, Methodology, Supervision, Writing – review & editing.
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgement
This work was supported by the “Bundesministerium für Bildung und Forschung” (BMBF) through the research grants No. 05M18BCA and 05M18VUA.
References (63)
- et al.
Mechanical damage behaviour of an injected sand by periodic homogenization method
Comput. Mater. Sci.
(2007) An extended homogenized porous electrode model for lithium-ion cell electrodes
J. Power Sources
(2015)- et al.
Derivation of macroscopic filtration law for transient linear viscoelastic fluid flow in porous media
Int. J. Eng. Sci.
(2000) - et al.
Influence of microstructure on the electrochemical performance of spinel cathodes in rechargeable lithium batteries
J. Power Sources
(2004) - et al.
Temperature dependent ageing mechanisms in lithium-ion batteries – a post-mortem study
J. Power Sources
(2014) - et al.
Homogenization of a reaction-diffusion equation with Robin interface conditions
Appl. Math. Lett.
(2006) - et al.
Generation of virtual lithium-ion battery electrode microstructures based on spatial stochastic modeling
Comput. Mater. Sci.
(2018) - et al.
Virtual testing of composites: imposing periodic boundary conditions on general finite element meshes
Compos. Struct.
(2017) - et al.
A mortar based approach for the enforcement of periodic boundary conditions on arbitrarily generated meshes
Comput. Methods Appl. Mech. Eng.
(2014) - et al.
Automatic three-dimensional geometry and mesh generation of periodic representative volume elements for matrix-inclusion composites
Adv. Eng. Softw.
(2016)
Periodic three-dimensional mesh generation for particle reinforced composites with application to metal matrix composites
Int. J. Solids Struct.
Structural characterization of particle systems using spherical harmonics
Mater. Charact.
Stochastic 3D modeling of the microstructure of lithium-ion battery anodes via Gaussian random fields on the sphere
Comput. Mater. Sci.
Parametric microstructure modeling of compressed cathode materials for Li-ion batteries
Comput. Mater. Sci.
On the origin and application of the Bruggeman Correlation for analysing transport phenomena in electrochemical systems
Curr. Opin. Chem. Eng.
Solutions of the flow equation for the drying of soils and other porous media
Soil Sci. Soc. Am. J.
A new model for predicting the hydraulic conductivity of unsaturated porous media
Water Resour. Res.
Numerical homogenization of hardened cement paste
Comput. Mech.
Effective transport properties of porous electrochemical materials – a homogenization approach
J. Electrochem. Soc.
Derivation of micro/macro lithium battery models from homogenization
Transp. Porous Media
Random Heterogeneous Materials: Microstructure and Macroscopic Properties
Influence of microstructure on impedance response in intercalation electrodes
J. Electrochem. Soc.
Morphology effects on the electrochemical performance of
J. Electrochem. Soc.
Increasing the discharge rate capability of lithium-ion cells with laser-structured graphite anodes: modeling and simulation
J. Electrochem. Soc.
Laser-induced breakdown spectroscopy of laser-structured Li(NiMnCo)O2 electrodes for lithium-ion batteries
J. Electrochem. Soc.
Manufacturing process for improved ultra-thick cathodes in high-energy lithium-ion batteries
Energy Technol.
Homogenization and two-scale convergence
SIAM J. Math. Anal.
The periodic unfolding method in domains with holes
SIAM J. Math. Anal.
Homogenization of the Stokes flow in a connected porous medium
Asymptot. Anal.
A mesh adaptation procedure for periodic domains
Int. J. Numer. Methods Eng.
TetGen, a Delaunay-based quality tetrahedral mesh generator
ACM Trans. Math. Softw.
Cited by (7)
A data-driven modeling approach to quantify morphology effects on transport properties in nanostructured NMC particles
2023, International Journal of Solids and StructuresA modelling framework for efficient reduced order simulations of parametrised lithium-ion battery cells
2023, European Journal of Applied MathematicsBridging nano- and microscale X-ray tomography for battery research by leveraging artificial intelligence
2022, Nature Nanotechnology