Abstract
Following the previous work on a simulation of powder compaction on two-dimensional model, where 400 discs were used and the cohesion probability P' and insertion probability P" were introduced, 1600 discs (40×40) were used and it was assumed that discs did not shift to the hole present in the lower layer when the coordination number of a disc or the probability variable (R) related to the coordination number was greater than a given value (q). The elimination or jumping probability P'" were newly introduced. Almost all the data of porosity and tapping number followed the first-order rate equation and the rate constant k increased with decreasing the probability P' and with increasing the probability P". As the tapping number increased, the porosity approached a definite value, which varied according to the magnitude of the values of q and P"'.
