Quantitative homogenization and fluctuations
Distinguish qualitative averaging, quantitative dependence assumptions, estimator bias, and central-limit fluctuations.
On this page
The question for today
What you’ll understand
- Explain why ergodicity alone provides no prescribed algebraic rate.
- Identify which dependence and observable assumptions enter quantitative statements.
- Compute exact resistance fluctuations and an asymptotic conductivity fluctuation law.
Knowing that an effective coefficient exists does not tell an engineer how large a specimen must be or how uncertain its estimated response remains. Quantitative homogenization addresses those questions under additional assumptions on the random environment. A rate must name the observable, norm, dimension, geometry, and dependence structure. We use an exactly solvable one-dimensional estimator to distinguish variance, nonlinear bias, and an asymptotic central limit law. Its explicit constants belong to that estimator, not to every random elliptic equation.
Dependence assumptions determine available rates
Stationarity and ergodicity allow arbitrarily slow averaging, so they do not imply a universal algebraic error bound. A much stronger example is finite-range dependence: the sigma fields generated by coefficients in two regions are independent whenever their separation exceeds a fixed range. Together with uniform ellipticity and boundedness, such hypotheses support quantitative stochastic homogenization and large-scale regularity results. Other theories use specified functional inequalities or quantified mixing assumptions. These conditions are alternatives with their own consequences, not interchangeable consequences of ergodicity.
Even after a dependence assumption is fixed, the target matters. An averaged flux, a finite-volume effective tensor, a solution in L², and a pointwise gradient can have different rates and dimensional corrections. Boundary conditions and coefficient regularity also affect estimates. Decay of one pairwise covariance does not automatically control dependence of nonlinear local functions or prove the functional inequalities a theorem uses. Before importing a rate, state exactly which random observable is estimated and which assumptions of the cited theorem apply to it.
The nonlinear response has both bias and fluctuations
For independent cells with conductivities one and four, let ri be the cell resistance and let r̄N be its sample mean. The finite-specimen conductivity with a unit macroscopic slope is aN=1/r̄N. Resistance averages linearly, so its expectation and variance are easy to compute. Conductivity is a nonlinear transformation: its expectation is generally not the reciprocal of expected resistance. Convexity of the reciprocal on positive numbers gives an upward finite-sample bias. Centering a confidence calculation at the limiting effective value therefore requires attention to bias as well as variance.
The classical independent central limit theorem applies to resistance because it is bounded with nonzero variance. Expanding the reciprocal at the limiting mean yields the conductivity central limit law by the delta method. This asymptotic statement describes the distribution after multiplication by √N as N tends to infinity. It is not an exact Gaussian distribution for a small specimen, and its asymptotic variance is not the exact finite-N variance. Quantitative remainder estimates need their own bounds, although the one-dimensional variables here are exceptionally well controlled.
Sampling error is only one component of a homogenization calculation
In computations, the sample volume, boundary treatment, discretization, and number of independent realizations all influence uncertainty. Repeating a biased finite-volume estimator reduces its sampling variance around its own mean; it does not automatically remove its boundary bias. Enlarging a domain can reduce both bias and fluctuations, but with different scalings. A useful error budget keeps the discretization error, finite-volume systematic error, and random statistical error separate, rather than treating a single empirical standard deviation as the full error certificate.
The exact one-dimensional calculation teaches a transferable discipline rather than a universal exponent. First identify the random quantity that averages linearly. Next understand how the desired response transforms it. Finally check whether the dependence assumptions justify a law of large numbers, a central limit theorem, or a quantitative tail estimate. In multidimensional random PDEs the corresponding response can involve global correctors and long-range correlations even when the underlying coefficients have short-range dependence. That is why fluctuation theory requires more than counting independent microscopic cells.
Work it through
Propagate resistance fluctuations through a reciprocal
Derive the limiting variance of √N(aN−8/5) for independent equal-probability conductivities one and four.
Write ri=1−(3/4)Bi with Bi Bernoulli of parameter one half. Then Eri=5/8 and Var(ri)=(9/16)(1/4)=9/64, giving Var(r̄N)=9/(64N).
The independent central limit theorem gives √N(r̄N−5/8) converging to a centered normal variable of variance 9/64. The reciprocal function g(r)=1/r is differentiable near 5/8 and has derivative g′(5/8)=−64/25.
The delta method multiplies the limiting variance by the squared derivative. Thus (64/25)²(9/64)=576/625. The asymptotic standard deviation of aN is therefore 24/(25√N), while its exact small-sample distribution remains discrete.
The conductivity fluctuation constant is 576/625. It follows from independent resistance averaging plus a nonlinear transformation, not from a general PDE rate.
The living notebook
Explore the model
Interactive lab · sampling
Check an exact variance before transforming it
Reset, set success probability p=0.5 and sample size n=25; predict the sample-proportion variance 0.01. For a fresh comparison Reset, set p=0.5 and n=100; predict 0.0025. Since resistance mean equals 1−0.75 times the sample proportion, its variances are 0.005625 and 0.00140625. The lab gives the Bernoulli sampling law; the conductivity variance requires the reciprocal transformation and is only asymptotically given by the delta-method formula.
- Mean sample proportion
- 0.5
- Sample-mean variance
- 0.0125
- Uniform CDF error bound
- ≤ 0.1062
Exact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forQuadrupling independent sample size quarters the exact variance of the linear resistance average.
Pause & reason
Check your understanding
Your turn
Try it yourself
For N=2 cells, enumerate the possible conductivities a₂ and calculate its expectation. Compare it with the infinite-volume value 8/5.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
There can be zero, one, or two cells of conductivity four, with probabilities one quarter, one half, and one quarter.
Compare with the solution
With no high-conductivity cells the mean resistance is one and a₂=1. With one such cell, mean resistance is 5/8 and a₂=8/5. With two, mean resistance is 1/4 and a₂=4. Hence Ea₂=(1/4)·1+(1/2)·(8/5)+(1/4)·4=41/20. Its upward bias relative to 8/5 is 9/20. This exact finite-sample bias illustrates the convex reciprocal transformation; it does not contradict almost-sure convergence or the centered asymptotic central limit law as N grows.
Keep the idea
Take it with you
- Quantitative rates need stronger dependence hypotheses than ergodicity.
- Nonlinear effective responses can have finite-sample bias as well as fluctuations.
- An asymptotic fluctuation law is distinct from an exact finite-volume error formula.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.