균일한 지하수 유속을 가진 불균질한 등방성 대수층 내에서 정류상태로 흐르는 지하수의 흐름과 함께 이동해가는 비반응성 오염물질에 대한 삼차원 몬테카를로 시뮬레이션이 시행되었다. 로그-정규적으로 분포되어 있는 수리전도도 가 임의장으로 설정되었으며 시뮬레이션 동안에 발생할 수 있는 불확실성을 감소하기 위해 여러 가지 방법들이 시도되었다. 3,200개 오염운들에 대한 이차공간적률의 집합적평균 , 그리고 오염운중심분산 이 각기 다른 세가지 불균질도 1.0, 2.5 및 5.0에 대해서 시뮬레이션 되었으며 또한 각기 다른 크기의 평균속도에 수직방향인 선형초기오염원(1.5 및 10)에 대해서 입자추적이 행하여 졌다. 시뮬레이션된 무차원 종적률들은 일차근사법에 의한 비에르고딕 이론적 결과와 비교적 잘 일치하나 시뮬레이션된 무차원 횡적률들은 일차근사법에 의한 이론적 결과들과 잘 일치하지 않으며 특히 불균질도가 큰 대수층에 대해서 그리고 초기 선형오염운의 크기가 작은 무차원 횡이차공간적률에 대해서 뚜렷하게 저평가 했다. 시뮬레이션된 집합적 평균이차적률은 에르고딕 상태에 도달하지 못했으며 일차근사법에 의한 에르고딕 용질이동에 관한 횡이차공간적률은 시뮬레이션 결과를 저평가 했음을 보인다.

three dimensional Monte-Carlo simulations of non-ergodic transport of a non-reactive solute plume by steady-state groundwater flow under a uniform mean velocity in isotropic heterogeneous aquifers were conducted. The log-normally distributed hydraulic conductivity, K(x), is modeled as a random field. Significant efforts are made to reduce the simulation uncertainties. Ensemble averages of the second spatial moments of the plume, and plume centroid variances, were simulated with 3200 Monte Carlo runs for three variances of log K, = 1.0, 2.5, and 5.0, and three dimensionless lengths of line plume sources (1,5 and 10) normal to the mean velocity. The simulated second spatial moment and the plume centroid variance in longitudinal direction fit well to the first order theoretical results while the simulated transverse moments are not fit well with the first order results. The first order theoretical results definitely underestimated the simulated transverse second spatial moments for the aquifers of large and small initial plume sources. The ergodic condition for the second patial moments is far from reaching, and the first order theoretical results of the transverse second spatial moment of the ergodic plume slightly underestimated the simulated moments.