The article presents an analysis of electromagnetic transient processes in long ultra-high voltage transmission lines, taking into account shield wires. It was shown that EMTP and Matlab/Simulink software are currently widely used in the study of transient processes in power lines. The EMTP software package uses the finite element method when integrating the equation that mathematically models the distributed parameter transmission line. The Matlab/Simulink software uses the d’Alembert method. In both cases, it is not known how the boundary conditions for the partial differential equation that mathematically models the transmission line are determined. The power line is analyzed as a distributed parameter system, described by a second-order partial differential equation. The advantage of the proposed method of calculating the boundary conditions for the abovementioned equation is the use of boundary conditions of the second and third types of Neumann and Poincaré, which allowed us to take into account the mutual influence of shield wires and phase conductors of the line in one power system. On this basis, the methodology for obtaining time domain graphs, spatial distributions, and traveling wave distributions of voltages and currents for phase conductors and line shielding wires is presented. The results of a computer simulation of transient processes when switching on the power line, taking into account controlled phase commutation and under single-phase earth fault conditions, are presented. All calculation results of transient processes presented in the article were obtained exclusively using numerical methods.