This paper is concerned with the approximation of the effective conductivity σ(A, µ) associated to an elliptic operator ∇ x A(x, η)∇ x where foris a bounded elliptic random symmetric d × d matrix and η takes value in an ergodic probability space (X, µ).We extend this result to non-symmetric operators ∇ x (a + E(x, η))∇ x corresponding to diffusions in ergodic divergence free flows (a is d × d elliptic symmetric matrix and E(x, η) an ergodic skew-symmetric matrix); and to discrete operators corresponding to random walks on Z d with ergodic jump rates. The core of our result is to show that the ergodic Weyl decomposition associated to L 2 (X, µ) can almost surely be approximated by periodic Weyl decompositions with increasing periods, implying that semi-continuous variational formulae associated to L 2 (X, µ) can almost surely be approximated by variational formulae minimizing on periodic potential and solenoidal functions.