This work is a natural extension of the author’s previous work: “Multiple scattering theory for heterogeneous elastic continua with strong property fluctuation: theoretical fundamentals and applications” (arXiv:1706.09137 [physics.geo-ph]), which established the foundation for developing multiple scattering model for heterogeneous elastic continua with either weak or strong fluctuations in mass density and elastic stiffness. Polycrystalline material is another type of heterogeneous materials that widely exists in nature and extensively used in industry. In this work, the corresponding multiple scattering theory for polycrystalline materials with randomly oriented anisotropic crystallites is developed. To validate the theory, the theoretical results for a series of materials such as OFHC copper, 304 stainless steel, and Inconel 600 are compared to experimental measurements and the numerical results obtained using finite element simulations. Detailed analysis shows that the new theory is capable of predicting the dispersion and attenuation of polycrystals with satisfactory accuracy. The results also show the new model can give an estimate on the average grain size with a relative error equal to or less than ten percent. As applications in ultrasonic nondestructive evaluation, we calculated the dispersion and attenuation coefficient of one of the most important polycrystalline materials in aeronautics engineering: high-temperature titanium alloys. The effects of grain symmetry, grain size, and alloying elements on the dispersion and attenuation behaviors are examined. Key information is obtained which has significant implications for quantitatively evaluating the average grain size, monitoring the phase transition, and even estimating gradual change in chemical composition of titanium components in gas turbine engines. For applications in seismology, the velocities and Q-factors for both hexagonal and cubic polycrystalline iron models for the Earth’s uppermost inner core are obtained in the whole frequency range. Using the realistic material parameters of iron under the high temperature and high-pressure conditions calculated from ab initio simulations, the numerical results show that the Q-factors range from 0.001 to 0.05, which shows good agreement with that inferred from real seismic data. The new model predicts the velocity of longitudinal waves varies between ± 1% to ± 5 % relative to the Voight average velocity, while the velocity of transverse waves varies from ± 10% to ± 20%, which gives promising explanation to the abnormally slow transverse velocity observed in practical measurements. The numerical results support the conjecture that the Earth’s uppermost inner core is a solid polycrystalline medium. The comprehensive numerical examples show the new model is capable of capturing the most important scattering features of both ultrasonic and seismic waves with satisfactory accuracy. This work provides a universal, quantitative model for characterization of a large variety of polycrystalline materials. It also can be extended to incorporate more complicated microstructures, including ellipsoidal grains with or without textures, and even multi-phase polycrystalline materials. The new model demonstrates great potential of applications in ultrasonic nondestructive evaluation and inspection of aerospace and aeronautic structures. It also provides a theoretical framework for quantitative seismic data explanation and inversion for the material composition and structural formations of the Earth’s inner core.