“…),and T i j t = v i,t v j ,t sin(𝜃 i,t − 𝜃 j ,t )are defined as in [20], where 𝜈 i,t , 𝜈 j ,t are the voltage magnitudes of buses i and j at time t, and 𝜃 i,t , 𝜃 j ,t are the phase shift angles of buses i and j at time t, respectively. The power flow through the distribution lines is expressed as a function of the previously defined auxiliary variables in (10) to (11), where variables u i j i,t , u i j j ,t are used to define the voltages at the sending and receiving ends of each line depending on its status, g i j , b i j ,and b i j sh are the conductance, susceptance, and shunt admittance of power line ij, respectively. The distribution line flow is limited to the maximum line capacity S i j in (12).…”