Abstract
Electrical utilities apply condition monitoring on power transformers (PTs) to prevent unplanned outages and detect incipient faults. This monitoring is often done using dissolved gas analysis (DGA) coupled with engineering methods to interpret the data, however the obtained results lack accuracy and reproducibility. In order to improve accuracy, various advanced analytical methods have been proposed in the literature. Nonetheless, these methods are often hard to interpret by the decision-maker and require a substantial amount of failure records to be trained. In the context of the PTs, failure data quality is recurrently questionable, and failure records are scarce when compared to nonfailure records. This work tackles these challenges by proposing a novel unsupervised methodology for diagnosing PT condition. Differently from the supervised approaches in the literature, our method does not require the labeling of DGA records and incorporates a visual representation of the results in a 2D scatter plot to assist in interpretation. A modified clustering technique is used to classify the condition of different PTs using historical DGA data. Finally, well-known engineering methods are applied to interpret each of the obtained clusters. The approach was validated using data from two different real-world data sets provided by a generation company and a distribution system operator. The results highlight the advantages of the proposed approach and outperformed engineering methods (from IEC and IEEE standards) and companies legacy method. The approach was also validated on the public IEC TC10 database, showing the capability to achieve comparable accuracy with supervised learning methods from the literature. As a result of the methodology performance, both companies are currently using it in their daily DGA diagnosis.

Abstract
Modern power grids have high levels of distributed energy resources, automation, and inherent flexibility. Those characteristics have been proven to be favorable from an environmental, social and economic perspective. Despite the increased versatility, modern grids are becoming more vulnerable to high-impact low-probability (HILP) threats, particularly for the distribution networks. On one hand, this is due to the increasing frequency and severity of weather events and natural disasters. On the other hand, it is aggravated by the increased complexity of smart grids. Resilience is broadly defined as the capability of a system to mitigate the effects of and recover from HILP events, which is often confused with reliability that is concerned with low-impact high-probability (LIHP) ones. In this paper, a distribution system in Portugal is simulated to showcase how the utilization of flexibility and mobile energy resources (MERs) should be considered differently relative to HILP vs LIHP threats.

Abstract
This paper presents an approach to optimize the placement of fault indicator devices in distribution systems using the cross-entropy method and results from traffic simulations. The problem formulation takes into account the impact of the devices on restoration times and costs due to fines related to service interruption reliability indices. Candidate solutions to the problem are evaluated using sequential Monte Carlo simulations, where travel times of maintenance crews are sampled according to data acquired from mobility traffic simulations. Results show the applicability of the approach in different simulation scenarios and the benefits of installing the devices in distribution networks.

Abstract
In this paper, a multi-temporal optimal power flow (OPF) model for radial networks is proposed. The mathematical problem formulation is presented as a mixedinteger quadratically constrained programming (MIQCP) problem. The main core of the developed OPF problem is benefiting from the second-order conic programming (SOCP) approach while the quadratic constraints of the power flow equations have been efficiently handled. In the developed model, the dynamic behaviour of the electrical energy storage (EES) has been addressed for the day-ahead operation problem. In addition, the developed model is tested and verified for both normal and contingent events and the obtained results are satisfactory in terms of feasibility and optimality. In the islanded operation, a grid-forming unit is the main responsible for maintaining the voltage reference while other units behave as slave. The model is tested on the modified IEEE 33-bus network to verify the performance of the developed tool.

Abstract
This paper presents a methodology to estimate flexibility existing on TSO-DSO borderline, for the cases where multiple TSO-DSO connections exist (meshed grids). To do so, the work conducted exploits previous developments regarding flexibility representation through the adoption of active and reactive power flexibility maps and extends the concept for the cases where multiple TSO-DSO connection exists, using data-driven approach to determine the equivalent impedance between TSO nodes, preserving the anonymity regarding sensitive grid information, such as the topology. This paper also provides numerical validation followed by real-world demonstration of the methodology proposed.

Abstract
Today's power grid is in a transitional stage to cater to the needs of energy efficiency, climate change, and environmental targets. In the process of designing the future power grid, one of the most fundamental models to be utilized is AC optimal power flow (AC-OPF). Since the feasible space of AC-OPF is non-convex, the optimization models developed using it often result in multiple local minima. To avoid such computational challenges in solving optimization models, various relaxation methods have been developed in the past. In the literature, these relaxation methods are mainly tested on specific networks. However, the scalability of relaxation techniques on branch-flow-based AC-OPF is yet to be explored. In this context, this paper compares the performance of different relaxation methods with the well-established MATPOWER AC-OPF solver in terms of the mean square error (MSE), maximum squared error, minimum and maximum values of voltage magnitude, and the average simulation time. In addition, the scalability of these models is tested on various radial and mesh networks with nodes ranging from 33 to 6655 nodes and 9 to 6515 nodes, respectively. In this manner, the trade-off between computational complexity and solution accuracy is demonstrated and analyzed in depth. This provides an enhanced understanding of the suitability and efficiency of the compared relaxation methods, helping, in turn, the efficiency of optimization models for varying sizes and types (i.e., radial or meshed) of networks.