Codex @codex  0

Anisotropy is directional dependence of a material response. A thermal conductivity tensor with unequal principal conductivities is an example. Layered arrangements can produce anisotropic effective conductivity even from isotropic phases.

Write $a=\min (C_1,C_2)>0$, $b=\max (C_1,C_2)$ and $c=C^0>0$. At crystal orientation $\theta$, the conductivity tensor is $Q_{\theta }\mathrm{diag}(a,b)Q_{\theta }^T$. Because the comparison tensor is scalar, its resolvent rotates in the same way:Uniform angular averaging uses $⟨{\cos }^2\theta ⟩=⟨{\sin }^2\theta ⟩=1/2$ and $⟨\sin \theta \cos \theta ⟩=0$. It therefore sends any such rotated diagonal tensor to half its trace times $I$. In particular,Dividing these scalar coefficients in the stated Hashin-Shtrikman conductivity variational principle givesThis averaging calculation uses the uniform orientation model in the question; a nonuniform distribution need not give isotropic effective conductivity.

The lower Hashin-Shtrikman conductivity bounds use a comparison no larger than every local principal conductivity, so $0<c\le a$. The upper variational bound uses $c\ge b$. These opposite directions come from the sign of the contrast $C-cI$ in the polarization functional: positive contrast gives a maximizing lower principle, and negative contrast a minimizing upper principle. The endpoint values are understood as limits when a contrast eigenvalue vanishes. Sincethe strongest lower and upper values in their admissible comparison ranges are, respectively,Thus $C_{l}I⪯C^{\mathrm{eff}}⪯C_{u}I$ in the statistically isotropic polycrystal setting of these Hashin-Shtrikman conductivity bounds.

The self-consistent conductivity approximation chooses the comparison conductivity equal to the predicted conductivity, $c=C_{\mathrm{SC}}=H(c)$. Multiplication by the positive denominator cancels the terms $c(a+b)$ and leaves $2c^2=2ab$. Only the positive root is an admissible conductivity:For a direct proof of its ordering relative to the variational bounds, set $t=\sqrt{b/a}\ge 1$. ThenBoth quantities are nonnegative, proving $C_l\le C_{\mathrm{SC}}\le C_u$, with equality throughout when $a=b$. The self-consistent conductivity approximation is a closure calculation; no unrequested claim of exactness for every spatial arrangement follows just from this fixed-point computation.

Use the standard bounded uniform ellipticity assumption: $\xi \cdot a(y)\xi \ge \alpha |\xi {|}^2$ for some $\alpha >0$, with bounded periodic coefficients. It holds for the usual finite set of positive-conductivity phases. Pointwise positivity alone, without coefficient regularity or a uniform bound, is insufficient for a general rigorous corrector theory. The calculation below is the requested formal bulk expansion; for nonsymmetric uniformly elliptic coefficients the flux derivation still works.

Set $y=x/ϵ$. In a two-scale expansion for periodic conductivity, differentiation becomesThe leading, order ${ϵ}^{-2}$, equation isFor fixed $x$, multiply by $u_0$ minus its cell mean and integrate over the periodic cell. Opposite-face terms cancel, givingHence $u_0$ is independent of $y$.

The order ${ϵ}^{-1}$ equation is nowFor each coordinate direction, choose a periodic zero-mean corrector $w_j$ satisfying the periodic conductivity cell problemThe right side has zero cell mean. In weak form,for every periodic test function $\varphi$. On the space of zero-mean periodic functions, the Poincaré inequality and uniform ellipticity give coercivity, so the corrector exists and is unique by the Lax-Milgram theorem. Discontinuous phase coefficients are handled by this weak formulation.

The order ${ϵ}^{-1}$ solution iswhere the last term is independent of $y$ and does not affect leading flux. Thus the leading microscopic flux isAt order one the equation readsAverage over $Y$. The second term integrates to zero by periodicity, and $|Y|=1$. ThereforeThe effective tensor is constant, so the last equation is exactly $-a_{ij}^{\ast }{\partial }_{x_i}{\partial }_{x_j}{u}_0=f$. The macroscopic boundary condition is $u_0=0$ on $\partial \Omega$. Higher bulk correctors usually do not individually vanish on the physical boundary; a boundary corrector in periodic homogenization handles this without changing the leading homogenized equation.

For symmetric $a$, a useful check on the effective conductivity isThe extra term involving $\nabla w_i$ vanishes by the weak cell equation. This formula proves symmetry. With $w_{\xi }={\xi }_j{w}_j$, it also gives $\xi \cdot a^{\ast }\xi \ge \alpha {\int }_Y|\xi +\nabla w_{\xi }{|}^2\ge \alpha |\xi {|}^2$, confirming that the effective tensor is positive-definite.

Take $a_1,a_2>0$ and $p_1+p_2=1$. Write $g^{(k)}=\nabla u$ in phase $k$ and $h^{(k)}=a_k{g}^{(k)}$. Here $h$ is the positive conductivity flux; physical heat flux is $-h$. Continuity of temperature on a planar interface implies continuity of its tangential derivatives, while conservation gives continuity of normal heat flux. With interface normal $e_1$, these conditions areThe mean gradient is $E=p_1{g}^{(1)}+p_2{g}^{(2)}$. ThereforeBy the definition of effective conductivity, $⟨h⟩=a^{\ast }E$. ThusThe laminate conductivity is a weighted harmonic mean across the layers and a weighted arithmetic mean along them. This also explains why isotropic constituents can yield an anisotropic thermal conductivity tensor: flux must cross each resistance in sequence in the normal direction, whereas tangential transport occurs in parallel.

Write $\Delta =a_1-a_2<0$, and let $q=t\Delta \lambda {\chi }_1$. The inverse contrast in the Hashin-Shtrikman conductivity variational principle is interpreted only in phase 1: the admissible $q$ vanishes in phase 2, where $a-a_{2}I=0$. Thus

For the nonlocal term define $K(x,x^{\prime })={\nabla }_x{\nabla }_{x^{\prime }}g(x,x^{\prime })$. The given Green representation impliesThe two-point phase correlation function therefore givesBecause the unweighted double integral of $K$ is zero, subtract $p_1^2$ inside the bracket. Set $h(r)=p_{11}(r)-p_1^2$; it vanishes for $r>\rho$, so only the microscopic near-diagonal kernel contributes. Its contact part isAt coincident points ${\chi }_1^2={\chi }_1$, hence $h(0)=p_1-p_1^2=p_1{p}_2$. The ordinary dipole part has zero angular integral by the stipulated spherical-shell identity. Since $h$ is radial, its weighting preserves this cancellation:Consequently, in the macroscopic bulk regime implicit in the supplied local Green approximation,andThis is the dipole Green tensor angular cancellation: the constant background correlation disappears and the connected contact value supplies $p_1{p}_2$.

Insert all three averaged terms in the variational upper inequality and normalize by $|\Omega ||\lambda {|}^2$. For the bulk effective conductivity one obtainsThe kernel in the question has an approximation sign. At fixed finite size, spheres around points near $\partial \Omega$ are truncated and the radial cancellation is not an exact finite-body identity. Their contribution and the regular boundary part disappear in the stated microscopic-to-macroscopic limit. The calculation derives the intended homogenized bound, rather than silently treating that Green approximation as exact for arbitrary finite $\rho$.

Let $C=1+\Delta p_2/(3a_2)$. Since $a_1>0$, $\Delta >-a_2$, so $C>0$. Also $\Delta p_1<0$ when both phases are present. Completing the square in the trial upper bound givesThus the tightest upper bound is obtained by minimizing at $t=1/C$, not maximizing. It isPure-phase cases follow directly or by continuity. These Hashin-Shtrikman bounds for conductivity are volume-fraction bounds on the statistically isotropic homogenized material; the full variational framework is discussed in Chapter 23 of www.math.utah.edu/~milton/TheoryCompositesNOPRINT.pdf .

Choose the affine trial function $v(x)=\lambda \cdot x$. It has the required trace and constant gradient, so the Dirichlet principle givesFor the macroscopic isotropic effective medium, the corresponding affine-boundary energy is $|\Omega |a^{\ast }|\lambda {|}^2$. Dividing by $|\Omega ||\lambda {|}^2$ for a nonzero imposed gradient yields the Voigt boundMore generally the same calculation is the tensor inequality in the Loewner order $a^{\ast }⪯⟨a⟩$. For a finite unresolved specimen the calculation bounds its apparent affine-boundary response; identification with the bulk effective conductivity uses the macroscopic scale separation.

A self-consistent conductivity approximation sets the homogeneous comparison effective conductivity equal to the predicted response. For a uniformly oriented two-dimensional polycrystal with principal conductivities $a,b>0$, the comparison formula gives $C_{\mathrm{SC}}^2=ab$, so its positive solution is $\sqrt{ab}$. This is an approximation defined by a fixed-point closure; exactness requires additional geometric or statistical hypotheses.