In modern high power communication distortion problems, the nonlinear mixing distortion effect in passive devices is always the most difficult to investigate, and it has attracted great interests in the modern high power and high density wireless communication system. Facing the typical nonlinear effect on passive device, Passive Intermodulation(PIM), this work concludes the common modeling methods and the important steps in PIM mechanism researches in recent years. In particular, targeting to the typical microscopic contact equivalence, PIM numerical transformation, and PIM prediction method in the engineered cases, the available methods are concluded for present, while a statistics based PIM prediction method for engineered device with dynamic range is proposed for the first time. This method can model the specific contact sources by using statistic strategy. The PIM prediction can be expressed as a statistical interval with the probability of each PIM value. This method can provide a big sample data for the engineered PIM prediction, and can guide the optimization of the PIM stability by controlling the detailed parameters.