An efficient percolation aggregation solver methodology captures media connectivity and continuity to reliably incorporate relevant flow solution trends in subterranean formation models. The approach allows introduction of meaningful physical information that is generally overlooked by state-of-the-art algebraic algorithms in the solution process. Percolation aggregation extends the efficiency and robustness of solution methods used to solve scientific and engineering problems.