E. Ya. Khruslov: “Homogenized Models of Suspension Dynamics”

The book under review is a research monograph on the homogenization for the motion of suspensions. Starting from a microscopic model describing the dynamics of mixtures of ﬂuid with ﬁne solid particles, the aim is the rigorous derivation of macroscopic models using multi-scale techniques. The homogenization process is performed under the assumption of an a priori known evolution of the particles, and different scenarios are considered in terms of the size of the particles, the distance between neighboring particles, and the interaction between different particles, leading to different macroscopic models. For this several multi-scale techniques are developed which are quite important for homogenization of ﬂuid problems in perforated (non-periodic) domains including small moving rigid particles. While there exists a wide literature on phenomenological ap-proaches for the derivation of macroscopic models for suspension based on physical observations, starting with the pioneering work [5], results on the rigorous mathematical justiﬁcations from ﬁrst principles seem to be rare. Here, the book offers signiﬁcant contributions by using homogenization theory; see the seminal work [2] for an introduction.

the solution of the macroscopic problem. We sketch this procedure for a model problem of diffusion or conductivity with a periodic coefficient: Let ⊂ R n for n ≥ 1 be open and bounded. For 0 < 1 we consider the problem −∇ · a x ∇u = f in , with f ∈ L 2 ( ) and a ∈ L ∞ per (R n ) n×n is symmetric, positive-definite, and Yperiodic, where Y = (0, 1) n is the unit cube in R n . It is easy to check by the Lax-Milgram lemma that (1) admits a unique weak solution u ∈ H 1 0 ( ), and the sequence (u ) >0 is bounded in H 1 ( ). However, even for such simple models the numerical treatment is quite challenging. Starting with the pioneering works [2] and [9], the homogenization theory deals with the derivation of macroscopic models for → 0 with effective/homogenized coefficients, the solution of which gives an approximation of the microscopic solution. In fact, for the solutions u of (1) it can be shown that there exists u ∈ H 1 0 ( ) such that u converges weakly in H 1 ( ) and strongly in L 2 ( ) to u, and u is the unique weak solution of the macroscopic problem −∇ · a * ∇u = f in , Here the homogenized coefficient a * ∈ R n×n is symmetric and positive definite, defined by (i, j = 1, . . . , n) a * ij = Y a(y)(∇ y w i + e i ) · (∇ y w j + e j )dy, with the unit vectors e i , and w i ∈ H 1 per (Y )/R (the space of functions in H 1 loc (R n ) which are Y -periodic and have mean value zero on Y ) are the unique weak solutions of the cell problems There is a huge literature on different methods in the homogenization theory dealing with the derivation of macroscopic models, see for example Tartar's method of oscillating test-functions [10], the two-scale convergence [1,8], or the unfolding method [4]. Another important approach for homogenization, especially for nonlinear problems, using variational methods is the -convergence, see [3,7]. This method is closely related to the method used in the book under review. We illustrate the methodology with the above example: The problem (1) is equivalent to the minimization of the functional where for ξ ∈ R n the function g hom is given by We extend the functionals I and I 0 to the space L 2 ( ) by setting Now, based on the methods of -convergence, see [3], it is possible to show the following results ( -convergence with respect to the L 2 ( )norm) An elementary calculation shows that g hom (ξ ) = a * ξ · ξ . In fact, the minimization problem in g hom for ξ = e i is equivalent to the cell problem (3). Hence, I 0 is the associated energy for the macroscopic problem (2). A crucial point in the homogenization theory is the identification and characterization of the homogenized coefficient a * respectively the homogenized integrand g hom , which also constitutes an important aspect in the book of Khruslov.
Let us describe in more detail the underlying microscopic model in the book under review: Let ⊂ R 3 be a fixed and bounded domain, filled with a viscous incompressible fluid with N = O | | 3 suspended small particles. Here, the parameter 0 < 1 describes the diameter and the distance between nearest neighboring particles. More precisely, is defined as the mean distance between the centers of mass of the i-th particle of the nearest particle in the initial state. The diameter d i of a particle depends on an additional parameter α ∈ [1,3] The domain occupied by the i-th particle for i ∈ {1, . . . , N } at time t ∈ [0, T ] is denoted by Q i (t), and the domain occupied by the fluid at time t is given by Now, the non-cylindrical domain of the i-th particle and the time-varying fluid domain are defined by The evolution of the fluid is described by the incompressible Navier-Stokes equations On the boundaries of the moving particles a stick condition is assumed: with the center of mass x i of the i-th particle, the velocity vector of the mass center u i , and the instantaneous angular velocity ω i . The motion of the i-th particle is described by a system of ordinary differential equations depending on external forces and stresses induced by the fluid. The system is closed with suitable initial and boundary conditions on ∂ , which are not the main focus of the monograph. For existence of a generalized solution for the microscopic problem for fixed the author refers to [6], where a similar problem is considered. However, existence can only be expected locally in time for a small time interval (0, T ), until the first collision of particles or with the boundary ∂ , and T may vanish for → 0.
The numerical treatment of the model (4a)-(4c) leads to high challenges, due to the different scalings like the size of the domain and the small diameter, and simulations in appropriate time scales are almost impossible. A possibility to overcome this problem is the derivation of so called macroscopic models for → 0 with homogenized (effective) coefficients, the solution of which approximates the solution of the microscopic model for fixed . The main part of the monograph is the rigorous derivation of such macro models for various relations between the given data, for example , d , and f . However, the treatment of the full time-dependent problem is rather complicated and perhaps impossible to solve. In fact, due to possible collisions (which can be expected in practical applications even after short times), the maximal existence time T for a solution in general tends to zero for → 0. Furthermore, in contrast to most results in homogenization theory, the evolution of the fluid domain (t) is an unknown of the system and leads to great challenges. Hence, the author uses the so called method of particle fixation to simplify the model. Thereby it is assumed that the evolution of the fluid domain and the particles is known, i.e., the position of the particles, and therefore the varying fluid domain (t), is known for any time t ∈ [0, T ]. Of course, this is a quite strong simplification, but at the moment there seems to be no rigorous mathematical framework for a homogenization of the full model (4a)-(4c) with the unknown evolution of (t). The method of fixing particle positions leads to a system of stationary Stokes equations with suitable boundary conditions on the surface of the particles. This is obtained by introducing a new force (assuming existence and regularity for the microscopic solution)f including the material derivative from the Navier-Stokes equations. The asymptotic analysis for these stationary Stokes systems for → 0 and the derivation of a macroscopic model is a main part of the book, and two geometrical settings for the particles are considered: Then, we are looking for a solution u ∈ • J ( ) of the minimization problem with the symmetric gradient e[v ] = 1 2 ∇v + ∇v T and external moments M i . • The filtering particle mode for d i = O( 3 ) and ρ s = O( −6 ), when particle velocities differ significantly from the mean fluid velocity. In this critical case α = 3 asymptotically the suspension behaves like a multiphase flow with interacting fluids. In this case, it is assumed that beside the centers of mass x i (t) of the i-th particle at time t, also its velocity u i (t) and angular velocity ω i (t) is a priori known.
An additional problem considered in the frozen particle mode is the motion of electrified or magnetized particles in a very strong magnetic or electric field, as well as particles interacting with each other through forces (e.g., electrostatic, elastic, Van der Waals) in both cases (filter and frozen particle mode). The case α > 3 is not considered in detail, as the particles then have almost no influence on the fluid. The principal idea for the derivation of the macroscopic model in the limit → 0 for the problem (5), respectively the problem for the filtering particle mode, is to use a direct approach from the calculus of variations. This could also be stated in the formal language of -convergence, see [7], but the author makes no relation to this method. Perhaps, a formulation in a more abstract framework offers additional opportunities for the treatment of other related problems. Let us exemplarily consider in more detail the homogenization of the problem (5), since the other cases are treated in a similar way. The aim is to show that, under suitable assumptions on the data and the perforated domain (t) for fixed t, the sequence u of minimizers of (5), which we denote by , converges to a divergence free limit function u, which minimizes the functional where : denotes the Frobenius inner product on R 3×3 ,f is the weak limit off , and the fourth order tensor a = (a npqr ) 3 n,p,q,r=1 is a homogenized limit tensor, the characterization of which is an essential part of the book. The proof consists of two main steps: is introduced via a minimization problem over small cubes K(ξ, h) with center ξ and side length h, such that h 1. Hence, the cubes are small compared to the whole domain , but large compared to the size of the particles and include many of them. In other words, with the parameter h an additional mesoscale is introduced, and the mean tensor of suspension a ( , h, ξ ) gives a mesoscopic characteristic of the mixture. Now, under the strong assumption that lim h→0 lim →0 a( , h, x) = a(x) for all x ∈ and a is continuous, it is possible to establish (B).
Altogether, the minimizer of gives an approximation of the minimizer of for every fixed time t ∈ [0, T ], if the evolution of the particles is known. Hence, the homogenized system describes the motion of a carrier fluid perturbed by particles moving in it. In this sense, the macroscopic system is not closed because it does not capture the influence of the fluid on the particles, and it includes unknown coefficients (the limit tensor a) which are depending on the dynamics of the particles and their properties. The closure of the homogenized system is the second main topic of the book. In a first step, the suspension viscosity tensor a is derived for some specific cases: • The particles are distributed locally periodic, leading to an expression of a in terms of solutions of cell problems. • Particles with an axisymmetric shape and arbitrary position, but diameters much less than the distances between neighboring particles (= low concentration). • Formula for the mean value of the viscosity tensor a for random distribution of diameters and orientations of particles.
Finally, the homogenized system is closed by the derivation of evolution equations for the mean orientation vector of • axisymmetric particles in the frozen particle mode, • spherical particles in filtration mode.
For this, the forces and moments acting on a single particle moving in an arbitrary linear unperturbed fluid flow (later the macroscopic fluid flow) are calculated. Combining all these steps, a closed macroscopic model with homogenized coefficients is obtained. Since the closed system is not obtained directly by homogenization from the microscopic model, existence of a weak solution is not guaranteed. This is obtained by some additional assumptions and simplifications on the homogenized coefficients. An additional subject of the book is the treatment of "complex fluids", where the particles are interacting with each other through different forces, in both cases, the frozen and the filtering particle mode. For the homogenization a reduced (linearized) model of (4a)-(4c) with some additional interaction forces between the particles is considered, where especially the microscopic domain is not evolving. In the frozen particle mode, the macroscopic system corresponds to an incompressible viscoelastic medium model including homogenized coefficients, whereas in the filtering particle mode it is shown for spherical particles, that the homogenized system is a two-phase model which can be interpreted as a generalized Brinkman law.
In summary, the monograph deals with an important problem for applications and is also highly interesting from a mathematical point of view. Starting from a general microscopic model, which takes into account the whole complexity of the system, especially the unknown dynamics of the particles, closed macroscopic models are derived by using the method of fixing particles. This seems to be necessary to overcome the difficulties arising due to the evolution of the domain. Nevertheless, the results are still highly nontrivial and state of the art. The monograph is appropriate for readers who are experts in the fields of homogenization and partial differential equations, especially in fluid dynamics. The proofs include methods from the existence and regularity theory for Stokes and Navier-Stokes equations (also for fluids with moving rigid bodies), variational methods, operator and spectral theory, as well as functional analysis. There is no introduction to these methods and the proofs including existing theory from the literature are quite short without many details, which makes it hard to read also for graduated students.
It would be highly interesting to include the results in more general frameworks, like the -convergence, the method of two-scale convergence [1], or the periodic unfolding method [4]. The latter two are developed for locally periodic or stochastic problems. However, for some of the results in the book these methods could be applicable, for example under the assumption of the existence of a transformation to a periodic domain, which should be closely related to the assumption that the suspension velocity tensor a( , h, x) converges to a continuous tensor a. A combination of such results could be beneficial for other applications in this field, for example considering elastic particles instead of rigid ones, leading to a fluid structure interaction between fluid and deformable solids, where it is necessary to take into account elasticity equations.
Funding Note Open Access funding enabled and organized by Projekt DEAL.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/ 4.0/.