# Fundamental Solutions for Periodic Media

Authors: Sergey V. Kuznetsov
Journal: Advances in Mathematical Physics
 (2014)

---

## Abstract

Necessity for the periodic fundamental solutions arises when the periodic boundary value problems should be analyzed. The latter are naturally related to problems of finding the homogenized properties of the dispersed composites, porous media, and media with uniformly distributed microcracks or dislocations. Construction of the periodic fundamental solutions is done in terms of the convergent series in harmonic polynomials. An example of the periodic fundamental solution for the anisotropic porous medium is constructed, along with the simplified lower bound estimate.


---

## Body

## 1. Introduction

The closest solutions in mechanics of heterogeneous media with uniformly distributed inclusions can be obtained from application of the two-scale asymptotic analysis [1–5]. In the two-scale asymptotic method, it is assumed that two fields exist: (i) the global field which is described by “slow” variables and (ii) a local field, having high oscillations; these are described by the “fast” variables; see Figure 1. Application of the two-scale asymptotic analysis to porous medium will be considered in more detail later on.Figure 1
“Slow” and “fast” varying processes (horizontal axis can refer to either time or distance).In the two-scale asymptotic method the effective elasticity tensor can be represented by(1)C0=∑p=1NfpCp+K,∑p=1Nfp=1,
where C0 is the effective elasticity tensor, fp is the volume fracture of the pth component, Cp is the elasticity tensor of the pth component, N is the total number of different components of the heterogeneous material, and K is a correcting tensor or “corrector.” It is clear from (1) that the main difficulty in determination of the effective elasticity tensor is in finding the corrector.Remark 1.
It is interesting to note that (1) covers almost all existing methods of homogenization by choosing different expressions for the corrector. (a)
Thus, ifK=0 the well-known Voigt’s homogenization is obtained.(b)
Taking(2)K=∑p=1N-fpCp+fpCp-1-1and assuming that for anyp tensor Cp is invertible along with fpCp-1, the Reuss homogenization for the elasticity tensor comes out. Assumption that Cp is invertible for any p is not valid for media with pores; in this case, the Reuss homogenization produces wrong values for the homogenized elasticity tensor.Determination of the corrector in the two-scale asymptotic method demands the solution of the cell problem, which in turn consists of (i) a boundary value problem on the internal boundaries between inclusions and the matrix material and (ii) a periodic boundary value problem on the outer boundary of a cell. The latter one is nonclassical in the sense that it is formulated on the boundary which, due to periodicity, must have angular points and edges; see Figure2.Figure 2
The typical periodic cell.Along with FEM and finite differences methods, the following other methods for obtaining a solution to the cell problem are known. The methods based on Eshelby’s transformation strain were applied to analysis of isotropic medium with ellipsoidal inclusions [6–8]. The advantage of these methods resides in their principle possibility to analyze media with anisotropic components, while from the computational point of view these are not very convenient since they lead to the three-dimensional integral equations with weakly singular kernels.The media with isotropic components were studied by a method based on the periodic fundamental solution for isotropic medium [9, 10]; such a fundamental solution was originally constructed in [11]. Because of multipolar expansions used for the solution of the inner boundary value problem, this method is confined to inclusions of spherical form. A similar approach was also used in the case of isotropic composites, but it was based on the Galerkin technique for solution of the inner boundary value problem [12].Presumably, for the first time periodic fundamental solutions for media with arbitrary anisotropy were developed in [13]. In combination with the boundary integral equation method (BIEM) these fundamental solutions were applied to the cell problem for composites with anisotropic inhomogeneities and porous media in [14, 15], while analysis of microstructural stresses in the matrix material was considered in [16]. Problems of wave scattering by pores were studied in [17, 18] by the same method. Some of obvious advantages of this method are due to their potential possibility to reduce the solution of the inner boundary value problem to a summation of rapidly convergent series, while periodic boundary conditions on the outer boundary are satisfied automatically due to periodicity of the fundamental solution.The following analyses cover both construction and properties of the periodic fundamental solutions and application of these solutions to analysis of porous anisotropic media with uniformly distributed pores.


## 2. Basic Notations

Initially a homogeneous elastic anisotropic medium is considered. The equations of equilibrium can be written in the following form:(3)A(∂x)u=-divxC0··∇xu=0,
where u is a displacement field. It is assumed that the tensor of elasticity satisfies the condition of positive definiteness, which is generally adopted for problems of mechanics of inhomogeneous media.Applying the Fourier transform(4)f∧(ξ)=∫f(x)exp⁡(2πix·ξ)dx,ξ∈R3,
to (3) gives the following symbol of the operator A:
(5)A∧ξ=2π2ξ·C0·ξ.From the definition of the fundamental solutionE of (3) the following formula for the corresponding symbol can be written:
(6)E∧ξ=A∧ξ-1.
Formula (6) shows that symbol E∧ is also strongly elliptic, positively homogeneous of degree −2 with respect to ξ, and analytical everywhere in R3∖0.Remark 2.
The Fourier inversion of the expression (6) and procedures for construction of the fundamental solution are discussed in [19–21].


## 3. Periodic Fundamental Solution

Consider a homogeneous anisotropic medium, loaded by the periodically distributed force singularities, located in nodesm of a spatial lattice Λ, Figure 3.Figure 3
Periodic lattice.Letai  (i=1,2,3) be linearly independent vectors of the main periods of the lattice, so that each of the nodes can be represented in the following form:
(7)m=∑imiai,
where mi∈Z are integer-valued coordinates of the node m in the basis (ai). The adjoint basis (ai*) is introduced in such a manner that ai*·m=mi. The lattice corresponding to the adjoint basis is denoted by Λ*.Now, periodic delta-function corresponding to singularities disposed in the nodes of the latticeΛ has the following form:
(8)δp(x)=VQ-1∑m*∈Λ*exp⁡(-2πix·m*),
where VQ is the volume of the fundamental region (cell) Q. Formula (8) defines a periodic delta-function uniquely.Substitution of the periodic fundamental solutionEp in (3) should produce
(9)A(∂x)Ep(x)=δp(x)I,
where I is the identity matrix. Looking for Ep also in the form of harmonic series and taking into account representation (8), it is possible to obtain
(10)Ep(x)=VQ-1∑m*∈Λ0*E∧(m*)exp⁡(-2πix·m*),
where Λ0* is an adjoint lattice without the zero node. It should be noted that formula (10) defines a periodic fundamental solution up to an additive (tensorial) constant.Lemma 3 (see [14]).
The series on the right side of (10) is convergent in the L1-topology, defining fundamental solution of the class L1¯(Q,R3⊗R3), where L1¯ is a class of integrable in Q functions with the zero mean value.


## 4. Effective Elasticity Tensor

For clarity and simplicity it will be assumed further that the considered medium has only one kind of uniformly distributed inhomogeneities placed in nodes of spatial latticeΛ (refer to Figure 3). The region occupied by an individual inhomogeneity in a cell Q is denoted by Ω.Two-scale asymptotic analyses being applied to such a medium produce the following expression for the corrector [14]:
(11)K=-VQ-1∫∂ΩC··νY⊗HYdY,
where Y are “fast” variables and H is the third-order tensorial field, being a solution of the following boundary value problem:
(12)A(∂Y)H(Y)=0,Y∈Q∖Ω,T(νY,∂Y)H(Y)∂Ω=-νY·C.
In (11) and (12)  νY represents the field of external unit normal to the boundary ∂Ω, and the elasticity tensor C is defined by
(13)C=C2-C1,
where C2 is referred to the matrix material and C1 to inclusions. Strong ellipticity of the tensor C is also assumed. It should be noted that in (13) tensor C1 can vanish in the case of voids.Lemma 4 (see [14]).
Under assumptions stated above, the boundary value problem (12) admits the unique solution.Solution of the boundary value problem (12) is constructed by the boundary integral equation method, giving the following representation for the desired solution [14]:
(14)12I+SH(Y′)=HcY′∈∂Ω,
where Hc is a constant tensor and S is a singular integral operator resulting from a restriction of the double-layer potential on the surface ∂Ω. Some of the relevant properties of operator S are discussed in [15].Substitution of (10) for periodic fundamental solutions in the expression for the operator S allows obtaining a lower (in energy) bound for the corrector; that is,
(15)Kl=-8π2VQ-2∑m*∈Λ0*χ∧Ω2C··m*⊗E∧(m*)⊗m*··C,
where χ∧Ω is the Fourier image of the characteristic function of the region Ω. An expression for the upper bound can be obtained similarly [14, 15].Theorem 5 (see [14, 15]).
Series appearing on the right side of (15) is absolutely convergent, provided Ω is a proper open region in Q.Remark 6.
Proof of convergence for very thin inclusions, including cracks, is to be studied separately, as the special asymptotic analysis is needed.


## 5. Microstructural Stresses and Scattering Cross-Sections

As was shown in [15, 16], the energy level Wosc of microstructural highly oscillating stresses for porous medium is defined by
(16)Wosc=12ε0··K··ε0,
where ε0 represents the uniform deformation field and K is the corrector.Similarly, having applied terminology used in quantum mechanics, scattering cross-sectionS for the porous medium has the following form:
(17)S=1-f-1ε0··K··ε0ε0··C··ε0,
where f is the porous ratio and C=C2 is the elasticity tensor of the matrix.Remark 7.
Problems for obtaining expressions for microstructural stresses and scattering cross-sections related to media with inclusions of nonzero stiffness are not studied yet by the two-scale asymptotic analysis and periodic fundamental solutions.


## 6. Conclusions

The developed methodology allowed us to construct the spatially periodic fundamental solution (10) for periodic media with anisotropic components. On the basis of thus constructed periodic fundamental solution, the closed form expression for the corrector tensor for porous medium with anisotropic matrix was obtained (11), along with the approximate formula (15).An estimate for the scattering cross-section for porous medium with anisotropic matrix was derived; see (17).



---

*Source: 473068-2014-12-14.xml*