Abstract
An original graph-based model and algorithm for optimal industrial task scheduling is proposed in this article. The innovative algorithm designed, dubbed "Dijkstra optimal tasking" (DOT), is suitable for fully distributed task scheduling of autonomous industrial agents for optimal resource allocation, including energy use. The algorithm was designed starting from graph theory fundamentals, from the ground up, to guarantee a generic nature, making it applicable on a plethora of tasking problems and not case-specific. For any industrial setting in which mobile agents are responsible for accomplishing tasks across a site, the objective is to determine the optimal task schedule for each agent, which maximizes the speed of task achievement while minimizing the movement, thereby minimizing energy consumption cost. The DOT algorithm is presented in detail in this manuscript, starting from the conceptualization to the mathematical formulation based on graph theory, having a thorough computational implementation and a detailed algorithm benchmarking analysis. The choice of Dijkstra as opposed to other shortest path methods (namely, A* Search and Bellman-Ford) in the proposed graph-based model and algorithm was investigated and justified. An example of a real-world application based on a refinery site is modeled and simulated and the proposed algorithm's effectiveness and computational efficiency is duly evaluated. A dynamic obstacle course was incorporated to effectively demonstrate the proposed algorithm's applicability to real-world applications.

Abstract
This paper proposes a data-driven chance-constrained optimal gas-power flow (OGPF) calculation method without any prior assumption on the distribution of uncertainties of wind power generation. The Gaussian mixture model is employed to fit the uncertainty distribution, where the Bayesian nonparametric Dirichlet process is adopted to tune the component number. To facilitate the online application of the proposed methods, an online-offline double-track distribution construction approach is established, where the frequency of training the relatively time-consuming Dirichlet process Gaussian mixture model can be reduced. On account of the quadratic gas consumption expression of gas-fired generators as well as the linear decision rule based uncertainty mitigation mechanism, the chance constraints would become quadratic ones with quadratic terms of uncertainties, which makes the proposed model more intractable. An equivalent linear separable counterpart is then provided for the quadratic chance constraints, after which the intractable chance constraints could be converted into traditional linear ones. The convex-concave procedure is used to crack the nonconvex Weymouth equation in the gas network and the auxiliary quadratic equalities. Simulation results on two test systems validate the effectiveness of the proposed methods.

Abstract
Unpredictable faults always reduce the stability and reliability of the electrical system. The increasing use of renewable energy sources (RES) in recent decades has exacerbated power system problems. Micro grids (MG) participation in Ancillary Services (AS) market is a suitable solution for the optimal performance of power systems in these conditions. MGs can also maximize their profits by participating in the AS market. In this paper, the optimal stochastic bidding strategy in joint energy and AS (regulation up and regulation down, spinning reserve and non-spinning reserve) market is modeled. Uncertainties of wind speed and solar radiation are modeled using Weibull and Beta probability distribution functions (PDFs) and probability of call AS is computed for all available AS. Therefore, the risk of the bidding strategy is controlled using conditional value at risk (CVaR). ERCOT market simulation has been carried out in order to determine the participation of each generator in all of the mentioned markets for different prices of energy and also to present the bidding curve, based on real-world data.

Abstract
Energy hubs (EHs) are units in which multiple energy carriers are converted, conditioned and stored to simultaneously supply different forms of energy demands. In this research, the objective is to develop a new stochastic model for unit commitment in EHs including an intelligent electric vehicle (EV) parking lot, boiler, photovoltaic (PV) module, fuel cell, absorption chiller, electric heat pump, electric/thermal/ cooling storage systems, with electricity and natural gas (NG) as inputs and electricity, heat, cooling and NG as demands. The uncertainties of demands, PV power and initial energy of EV batteries are modeled with Monte Carlo Simulation. The effect of demand response and demand participation factors as well as effect of EVs and storage systems on EH operation are investigated. The results indicate that thermal demand response is more effective than electric and cooling demand response; as it decreases EH operation cost by 12%, while electric demand response and cooling demand response decrease it respectively by 9.3% and 4.2%. The results show that at low electric/thermal/cooling demand participation factors, an increase in participation factor sharply decreases EH operation cost, while the same amount of increase at higher participation factors leads to a smaller decrease in operation cost. The results also indicate that thermal storage system and cooling storage system have significant effect on reduction of EH operation cost, while the effect of electric storage system is trivial.

Abstract
High penetration of renewable energy sources and energy storage systems has considerably increased the flexibility in power distribution networks operation. However, employing converter-interfaced energy and storage sources may significantly reduce the mechanical inertia and as a result, the power grids may confront serious stability challenges during transient conditions. This article introduces a strategy for enhancing transient stability margin of active distribution networks with high penetration of electric vehicles (EVs). The proposed optimization strategy intends to control EVs contributions during transient stability conditions. The EVs contributions are controlled through a new index proposed based on the system's total corrected critical kinetic energy (TCCKE). The proposed procedure for TCCKE calculation is driven by a hybrid algorithm taking into account the equal area criterion and sensitivity analysis. The suggested procedure for TCCKE only depends on the during fault data and as a result, the proposed optimization strategy is useful to prevent transient instability in the case of first swing instability. The proposed optimization is applied and evaluated on the IEEE test systems. The results clearly demonstrate the applicability and efficacy during a multitude of fault and emergency conditions.

Abstract
Increasing the penetration of renewables on prosumers' side brings about operational challenges in the distribution grid due to their variable and uncertain behavior. In fact, these resources have increased the distribution grid net load fluctuation during recent years. In this article, the flexibility-oriented stochastic scheduling of a microgrid is suggested to capture the net load variability at the distribution grid level. In this scheduling, the flexibility limits are set to manage the net load fluctuation at a desirable level for the main grid operator. The uncertainties of load and renewables are considered, and their uncertainties are under control by the risk-averse strategy. Moreover, multiperiod islanding constraints are added to the problem, preparing the microgrid for a resilient response to disturbances. The model is examined on a typical distribution feeder consisting of prosumers and a microgrid. The numerical results are compared for both flexibility-oriented and traditional scheduling of a microgrid at the distribution level. The proposed model reduces the net load ramping of the distribution grid using an efficient dispatch of resources in the microgrid. A sensitivity analysis is also carried out to show the effectiveness of the model.