Electrorheological complex (dusty) plasmas (ER-CDPs) are type of plasmas with rheological behavior under external electric field (EEF). In this chapter, the nonequilibrium molecular dynamics simulations method is used to tune thermal conductivity for two-dimensional (2D) complex plasma liquids using EEF in different directions. Anisotropic thermal conductivity is investigated in three different cases, (i) Ex(x-axis), (ii) Ey(y-axis), and (iii) Exy(xy-axis) and with constant external perturbation force (Px = 0.02). The thermal conductivity under the influence of EEF is different in different directions. Obtained results are compared and discussed with previous known theoretical, simulation, and experimental data for 2D systems in constant EEF. The appropriate normalized thermal conductivity with Einstein frequency at constant EEF follows universal temperature scaling law. These comparisons and discussions show that algorithms of EEF of different cases with Yukawa potential have accuracy and consistency. These comparisons validated a new numerical model that can be used for variations of EEF along with different system sizes and plasma parameters.