**Abstract**

The computation of thermalt properties of dusty plasmas is substantial task in the area of science and technology. The thermal conductivity (λ) has been computed by applying polarization effect through molecular dynamics (MD) simulations of two dimensional (2D) strongly coupled complex dusty plasmas (SCCDPs). The effects of polarization on thermal conductivity have been measured for a wide range of Coulomb coupling (Γ) and Debye screening (κ) parameters using homogeneous non-equilibrium molecular dynamics (HNEMD) method for suitable system sizes. The HNEMD simulation method is employed at constant external force field strength (*F\** ) and varying polarization effects. The algorithm provides precise results with rapid convergence and minute dimension effects. The outcomes have been compared with earlier available simulation results of molecular dynamics, theoretical predictions and experimental results of complex dusty plasma liquids. The calculations show that the kinetic energy of SCCDPS depends upon the system temperature ( 1/Г) and it is independent of higher screening parameter. Furthermore, it has shown that the presented HNEMD method has more reliable results than those obtained through earlier known numerical methods.

**Keywords:** Plasma thermal conductivity, complex dusty plasma, Homogenous non-equilibrium molecular dynamics, force field strength, system size, plasma parameters etc.
