Tomasz N. Koltunowicz

The paper presents studies of the site percolation phenomenon for square matrixes with dimensions L = 55, 101 and 151 using the Monte Carlo computer simulation method. The number of iterations for each matrix was 5 × 106. An in-depth analysis of the test results using the metrological approach consisting of determining the uncertainty of estimating the results of iterations with statistical methods was performed. It was established that the statistical distribution of the percolation threshold value is a normal distribution. The coefficients of determination for the simulation results in approximations of the percolation threshold using the normal distribution for the number of iterations 5 × 106 are 0.9984, 0.9990 and 0.9993 for matrixes with dimensions 55, 101 and 151, respectively. The average value of the percolation threshold for relatively small numbers of iterations varies in a small range. For large numbers of iterations, this value stabilises and practically does not depend on the dimensions of the matrix. The value of the standard deviation of the percolation threshold for small numbers of iterations also fluctuates to a small extent. For a large number of iterations, the standard deviation values reach a steady state. Along with the increase in the dimensions of the matrix, there is a clear decrease in the value of the standard deviation. Its value is about 0.0243, about 0.01 and about 0.012 for matrixes with dimensions 55, 101 and 151 for the number of iterations 5 × 106. The mean values of the percolation threshold and the uncertainty of its determination are (0.5927046 ± 1.1 × 10−5), (0.5927072 ± 7.13 × 10−6) and (0.5927135 ± 5.33 × 10−6), respectively. It was found that the application of the metrological approach to the analysis of the percolation phenomenon simulation results allowed for the development of a new method of optimizing the determination and reducing the uncertainty of the percolation threshold estimation. It consists in selecting the dimensions of the matrix and the number of iterations in order to obtain the assumed uncertainty in determining the percolation threshold. Such a procedure can be used to simulate the percolation phenomenon and to estimate the value of the percolation threshold and its uncertainty in matrices with other matrix shapes than square ones.