Object-Oriented Usability Indices for Multi-Objective Demand Side Management Using Teaching-Learning Based Optimization

: This paper proposes Object-Oriented Usability Indices (OOUI) for multi-objective Demand Side Management (DSM). These indices quantify the achievements of multi-objective DSM in a power network. DSM can be considered as a method adopted by utilities to shed some load during peak load hours. Usually, there are service contracts, and the curtailments or dimming of load are automatically done by service providers based on contract provisions. This paper formulates three indices, namely peak power shaving, renewable energy integration, and an overall usability index. The ﬁrst two indices indicate the amount of peak load shaving and integration of renewable energy, while the third one combines the impact of both indices and quantiﬁes the overall beneﬁt achieved through DSM. The application of the proposed indices is presented through simulation performed in a grid-tied microgrid environment for a multi-objective DSM formulation. The adopted microgrid structure consists of three units of diesel generators and two renewable energy sources. Simulation has been done using MATLAB software. Teaching-Learning-Based Optimization (TLBO) is adopted as the optimization tool due to its simplicity and independency of algorithm-speciﬁc control parameters. Five different cases of renewable energy availability with results validate the efﬁciency of the proposed approach. The results indicate the usefulness in determining the suitable condition regarding


Motivation
DSM can be considered as the coordinated reduction in load during a specific time so as to maximize the energy usage from renewable energy and minimize the generation cost. Most of the research on DSM have tried to explore the technical benefit of DSM in design and operation area, recently some studies have tried to quantify the economic benefits of DSM. Thus motivation of the work is to quantify both the technical and economic benefits of DSM using OOUI.

Literature Review
Currently, the power industry is facing numerous challenges because of the fast-changing structure of the power network. The integration of small-size renewable energy sources with the conventional grid has become a significant challenge. In addition to this, the uncertainty of power availability from these renewable energy sources is also a serious concern. The search for a suitable scheme which provides a solution for technical and commercial challenges associated with the above-discussed

Microgrid Modeling
The schematic diagram of the microgrid is given in Figure 1. This microgrid consists of three units of diesel generators (DG1, DG2, DG3) and two renewable energy sources: one wind and one solar generator. The switching logic block in Figure 1 takes care of scheduling all six energy resources. Three different types of customer are assumed to be connected in this system using advanced metering infrastructure (AMI). AMI provides all the required information of the customer load to the utility. A scheduling interval of 24 h is used in this paper. Renewable energy sources are coupled with DGs and the main grid through a solid state interface converter which allows power flow from the grid and distributed energy resources (diesel, wind, solar energy resources) to loads connected at the point of common coupling. The maximum interchangeable power allowed between the main grid and microgrid is 4 kW. The schematic diagram of the microgrid is given in Figure 1. This microgrid consists of three units of diesel generators (DG1, DG2, DG3) and two renewable energy sources: one wind and one solar generator. The switching logic block in Figure 1 takes care of scheduling all six energy resources. Three different types of customer are assumed to be connected in this system using advanced metering infrastructure (AMI). AMI provides all the required information of the customer load to the utility. A scheduling interval of 24 h is used in this paper. Renewable energy sources are coupled with DGs and the main grid through a solid state interface converter which allows power flow from the grid and distributed energy resources (diesel, wind, solar energy resources) to loads connected at the point of common coupling. The maximum interchangeable power allowed between the main grid and microgrid is 4 kW.   Figure 4 shows the hourly power interruptibility cost( , ) in dollars. Descriptions of the dieselgeneration fuel cost function with its parameters are as follow: where , is the power output and , , are the fuel cost coefficients related to the lth DG at any time instant t; for example, if the cost function of DG2 is , = , + , then here and are parameters related to DG2 and , is the output power of DG2 during the 5th hour. Information about the conventional DG parameters are tabulated in Table 1, providing all the cost parameter data of the diesel generators. The solar and wind modeling presented here is based on a complex and detailed model. For the detailed model please referred to [12,31].   Figure 3 shows the initial hourly demand of three customers. Figure 4 shows the hourly power interruptibility cost(λ j,t ) in dollars. Descriptions of the dieselgeneration fuel cost function with its parameters are as follow: C l (P l,t ) = a l P 2 l,t + b l P l,t where P l,t is the power output and a l , b l , are the fuel cost coefficients related to the lth DG at any time instant t; for example, if the cost function of DG2 is C 2 (P 2,t ) = a 2 P 2 2,5 + b 2 P 2,5 then here a 2 and b 2 are parameters related to DG2 and P 2,5 is the output power of DG2 during the 5th hour. Information about the conventional DG parameters are tabulated in Table 1, providing all the cost parameter data of the diesel generators. individual load types (i.e., air conditioners, refrigerators, etc.), or presented as characteristic impedance equations (i.e., constant impedance loads, non-linear loads, etc.), this paper considers loads as cumulative consumer units. The flexibility of the load is presented through a term which indicates the willingness of individual customers to participate in the rescheduling process. The willingness value varies from 0 to 1, which indicates zero percent willingness to one-hundred percent willingness. The details of customer-related parameters are discussed in Section 3.2 and tabulated in Table 2.  The microgrid conceptually utilizes the potential of solar and wind energy generation units. The small-scale generation units based on renewable energy sources are normally interfaced into the conventional grid at the low-voltage distribution level. Such microgrid structures lower the cost, negating traditional economies [32]. The studied system adopted in this paper is a microgrid structure in grid-tied mode. The adopted microgrid structure consists of a solar-based generator, wind-based generator, and three conventional diesel generators at the source end.  The mathematical model describing the structure and behavior of individual components and the complete system is prepared and presented in the next sections.

Microgrid DSM Problem Formulation
This section presents the microgrid-DSM formulation. Optimization results of this formulation are used in the upcoming section to show the application and usability of the proposed indices.
Adopted microgrid DSM problem formulation has two objectives, namely, minimization of fuel cost function G(X) for diesel generators by maximum utilization of renewable energy sources, and maximization of utility benefit function using DSM. The formulation of the multi-objective function involves these two objectives, as shown in Figure 5.   The mathematical model describing the structure and behavior of individual components and the complete system is prepared and presented in the next sections.

Microgrid DSM Problem Formulation
This section presents the microgrid-DSM formulation. Optimization results of this formulation are used in the upcoming section to show the application and usability of the proposed indices.
Adopted microgrid DSM problem formulation has two objectives, namely, minimization of fuel cost function G(X) for diesel generators by maximum utilization of renewable energy sources, and maximization of utility benefit function using DSM. The formulation of the multi-objective function involves these two objectives, as shown in Figure 5.  The solar and wind modeling presented here is based on a complex and detailed model. For the detailed model please referred to [12,31].
Wind energy conversion units develop electrical energy from wind speed using the mathematical relation of Equation (1a), referred to in [12,31]: where η is the efficiency of the wind generator, A is the swept area of the wind turbine rotor; V is the velocity of the wind at the hub height, C p is the power coefficient of the wind turbine, and ρ is the air density. The PV generator output per hour is given below [12,31]: where η is the efficiency of the solar pv generator, A c is the area of the solar pv array, and I pvt is the per hour incident of solar irradiation (kWh/m 2 ) on the solar array. Load modeling is an important part of the demand response program as the main target of the program is to reschedule the consumer loads. The loads in this paper are taken as aggregate loads as individual customers. Unlike some other load modeling where loads are presented in terms of individual load types (i.e., air conditioners, refrigerators, etc.), or presented as characteristic impedance equations (i.e., constant impedance loads, non-linear loads, etc.), this paper considers loads as cumulative consumer units. The flexibility of the load is presented through a term θ i which indicates the willingness of individual customers to participate in the rescheduling process. The willingness value varies from 0 to 1, which indicates zero percent willingness to one-hundred percent willingness. The details of customer-related parameters are discussed in Section 3.2 and tabulated in Table 2. The microgrid conceptually utilizes the potential of solar and wind energy generation units. The small-scale generation units based on renewable energy sources are normally interfaced into the conventional grid at the low-voltage distribution level. Such microgrid structures lower the cost, negating traditional economies [32]. The studied system adopted in this paper is a microgrid structure in grid-tied mode. The adopted microgrid structure consists of a solar-based generator, wind-based generator, and three conventional diesel generators at the source end.
The mathematical model describing the structure and behavior of individual components and the complete system is prepared and presented in the next sections.

Microgrid DSM Problem Formulation
This section presents the microgrid-DSM formulation. Optimization results of this formulation are used in the upcoming section to show the application and usability of the proposed indices.
Adopted microgrid DSM problem formulation has two objectives, namely, minimization of fuel cost function G(X) for diesel generators by maximum utilization of renewable energy sources, and maximization of utility benefit function using DSM. The formulation of the multi-objective function involves these two objectives, as shown in Figure 5.   The mathematical model describing the structure and behavior of individual components and the complete system is prepared and presented in the next sections.

Microgrid DSM Problem Formulation
This section presents the microgrid-DSM formulation. Optimization results of this formulation are used in the upcoming section to show the application and usability of the proposed indices.
Adopted microgrid DSM problem formulation has two objectives, namely, minimization of fuel cost function G(X) for diesel generators by maximum utilization of renewable energy sources, and maximization of utility benefit function using DSM. The formulation of the multi-objective function involves these two objectives, as shown in Figure 5.

Objective1: Minimization of the Fuel Cost Function, G(X)
In this mode the reduction of the fuel cost of microgrid generators is the main objective. In this work, a power exchange scheme is assumed to exist between the main grid and microgrid. This  In this mode the reduction of the fuel cost of microgrid generators is the main objective. In this work, a power exchange scheme is assumed to exist between the main grid and microgrid. This scheme is needed to cater to the intermittent nature of renewables. When the microgrid system is not able to complete the load demand then power has to be taken from the main grid, and if the microgrid supply is a surplus with respect to its need then the excess power is given back to the main grid. Thus, this power exchange cost needs to be considered during modeling of the function to minimize the fuel cost of conventional generators of the microgrid. Equation (3) represent the required mathematical model of this cost function where C r (Pr t ) is the cost of transferable power. "γ t ", (dollars per kWh) is the rate of power exchange at a specific bus [12], whereas (P rt ) > 0 and (P rt ) < 0 reflect the power transfer between the microgrid and main grid, while the condition (P rt ) = 0 reflects there is no power exchange between the grid and microgrid. Positive and negative values reflect monetary gain and loss. The function related to minimization of fuel cost in this mode is as follows: where l, t, T, L hold integer values only: l: Variable to represent conventional generating units. Its value ranges from 1 ≤ l ≤ L; t: Dispatch interval, 1 ≤ t ≤ T; in present work, is expanded over T = 24 time horizon. G(X): Operating cost function. Cost t : Total cost at time t to deliver power; P l,t : Conventional generator, number of units; C l (P l,t ): Fuel cost function for conventional generators. C r (Pr t ): Transferable power cost; P rt : Transferable power; γ: Location marginal prices [33]; L: Total number of the conventional generating unit.
Equation (2) subjected to following constraint: Equation (4) signifies the power balance. Collective generation of a conventional microgrid must fulfill all demand inside the microgrid. A small three-feeder radial low voltage scheme is taken in the present study. Losses for the model under consideration are assumed as having very low numeric value. The rest of the constraints are for power generation capacity of conventional and renewable power generators. The constraint, defined by Equation (5), is for generation limits of the conventional generators. This ensures that limits are not violated. Constraints defined using Equations (6) and (7) stand for renewable energy generators to ensure that the optimal values for renewable energy generators are within the forecast range. The fourth constraint, Equation (8), is to restrict the exchange of power. The sixth constraint, Equation (9), takes care of the maximum allowed ramp up and down rates.

Objective2: Maximization of Utility Benefit Function, H(X), Using Demand Response
In general, demand response is viewed as the involvement of the end-user in the power system. This involvement is in the form of cost variation with time and obtained incentive. Demand response can be divided into two classes: the first one based on incentive, and the second based on a time-based rate program [34,35]. In this work we have adopted the first one.
The main interest of the utility is to maximize its benefit. The mathematical model of the function to represent this interest may be modeled as Equation (10). In the present work, we categorized electricity consumers in three different categories, and the following equations demonstrate the modeling for their behavior.
In this work, the quadratic outage cost function for more than one customer is used [36]. Equation (12) defines the customer cost function, where the "Z" term sorts users according to θ, whereas θ is giving information about the user type and this information is used to categorize the customers [36][37][38][39]. Details of parameters related to the customer cost function are tabulated in Table 2. Customer benefit is defined as the difference between the incentive received by him and the cost incurred by him. Utility benefit can be stated as the difference between the cost of the load curtailed by user and the incentive obtained by the user. The mathematical model of the benefit function is as follows: Benefit function of the utility: Equation (10) denotes the benefit function of the utility where term "λx" is the total hourly benefit to the utility when it does not deliver "x" kW power to a specific customer at the rate of "λ " dollars per kWh [12]. The term y is compensation in terms of incentive which is paid to customer for reducing its load by "x" kW.
Benefit function of the customer: where: V 1 : Customer benefit in $/kWh. It must be V 1 ≥ 0 for user participation. V 2 : Utility benefit in $/kWh. Equation (11) denotes the benefit function (V 1 ) of the customer, where c is the monetary loss to the user who reduces his power requirement by "x" kW when DSM is applied and "y" is monetary incentive value that the customer receives. The customer would only participate if V 1 ≥ 0, meaning he gets some benefit.
The customer cost function for the ith customers can be modeled as: Energies 2019, 12, 370 9 of 25 and K 1,i , K 2,i are coefficients related to the costincurred by customer i.
The main interest of the utility is to maximize its benefit (V 2 ). If the number of users is "i" then it can be mathematically expressed as max(λ i x i − y i ). A block diagram representation of this expression is shown by Figure 6.
The main interest of the utility is to maximize its benefit ( ). If the number of users is " "then it can be mathematically expressed as max( − ).A block diagram representation of this expression is shown by Figure 6. The above demand response model maximizes the expected benefit for the utility when implemented over the total optimization period of twenty-four hours. With considering maximum powertarget and the total budget as practical constraintsit can be modeled as: subject to: , ≤ , ≤ where , , , hold integer values only, and i ranges between1 ≤ ≤ , in Equation (15) and 2 ≤ ≤ , in Equation (16), and " " denotesthe total budget of the utility. The daily limit of power interruptibility for customer is denoted by . Equation (15) denotes the constraint which restricts the incentive of an end-user on a daily basis. The constraint of Equation (16) is to restrict the customer benefit, which must be greater for greater curtailed customer power. The constraint denoted by Equation (17) limits the total utility payment to end-users as an incentive within its budget. The constraint represented by Equation (18) controls the individual curtailed customer power (on total time horizon).

Multi-Objective Formulation
The DSM problem in this paper involves two objectives, fuel cost minimization and utility benefit maximization. Putting Equations (2) and (14) together represents the DSM problem as a multi- The above demand response model maximizes the expected benefit for the utility when implemented over the total optimization period of twenty-four hours. With considering maximum powertarget and the total budget as practical constraintsit can be modeled as: subject to: where i, n, t, T hold integer values only, and i ranges between 1 ≤ i ≤ n, in Equation (15) and 2 ≤ i ≤ n, in Equation (16), and "TB" denotesthe total budget of the utility. The daily limit of power interruptibility for customer i is denoted by CP i . Equation (15) denotes the constraint which restricts the incentive of an end-user on a daily basis. The constraint of Equation (16) is to restrict the customer benefit, which must be greater for greater curtailed customer power. The constraint denoted by Equation (17) limits the total utility payment to end-users as an incentive within its budget. The constraint represented by Equation (18) controls the individual curtailed customer power (on total time horizon).

Multi-Objective Formulation
The DSM problem in this paper involves two objectives, fuel cost minimization and utility benefit maximization. Putting Equations (2) and (14) together represents the DSM problem as a multi-objective optimization problem. However, the optimization tool adopted in this paper is TLBO, which is a single-objective optimization algorithm.
Thus, to implement TLBO both of the functions need to be converted into the same type. We converted the objective represented by Equation (14) from a maximization problem to a minimization problem by multiplying it with (−1). This transformation is represented in Equation (19). The two objectives represented by Equations (2) and (19) are scalarized by pre-multiplying each of them with a weighing factor "w" and "(1 − w)" such as the total weight: w + (1 − w) = The overall single-objective optimization problem is represented by Equation (20): where x i,t , y i,t , P w,t , P s,t , P r,t , P l,t are decision variables.
In this work we transformed the constrained objective function into the unconstrained objective by adding penalty terms for each constraint violation (Equations (4)- (9) and (15)-(18)). In the present work we a have given equal weighting to both objectives.
There are many methods to handle the constraint of an objective function in an optimization problem. In the present work, we designed the constraints of our objective function using an exterior penalty approach in which infeasible design variables are penalized, but feasible design variables are not penalized [40].

Teaching-Learing-Based Optimization (TLBO) Algorithm
Under this section TLBO is first introduced, and then its adaptation for multi-objective formulation has been presented.
In [29] six different meta-heuristic algorithms including TLBO are compared. These algorithms are, namely, Whale Optimization, Fire Fly, PSO, Differential Evaluation, and GA. This compression is done mainly in terms of selecting the best suitable optimization method for cost minimization of the micro-grid. In terms of convergence characteristics analysis among these methods TLBO reported the fastest convergence. Search agents of TLBO are found more effective to explore the feasible search space. It is reported in the [29] that TLBO search agents require the least effort. Among the above-mentioned six meta-heuristic techniques TLBO is reported as the best in most of the conditions.
The microgrid DSM is a multi-objective problem, and TLBO is a promising method for optimization. TLBO is proposed by [27,28]. The algorithm is divided in two phases, namely, teacher and learner. This phase of TLBO is explained below.
A. Teacher Phase: In this phase a teacher tries to improve the average result of the class in his subject. In this identification of the best solution from the population is done on the based on the objective function value.
If X k,i is the fittest solution at any iteration i then it can be represented as X j,kbest,i . Calculation of the mean result of the learner, in a subject j will be done in the next step. The effort of the teacher to improve the mean result of the class in his subject depends on the capability of the teacher. This improvement in the existing mean result is modeled as: where: X j,kbest,i : Result of the best learner in subject j; r j,i : Random number in the range [0, 1]; and T f : Teaching factor. It decides the value of mean to be changed the value of T f is selected in random way with equal probability as: updated solution in the teacher phase: where: X j,k,i : A value in the solution, X j,k,i : Updated value of X j,k,i . Accept X j,k,i if it gives better value of the function j: jth design variable. Denotes subject chosen by the learners. j = 1, 2, . . . , m; k: kth member of population. Denotes learner. k = 1, 2, . . . , n; i: ith iteration, i = 1, 2, . . . , G max , G max : Denotes maximum iterations.
All the values of the function accepted in this phase are reserved and treated as the learner phase input.
B. Learner phase: In this phase learners improve their knowledge by interaction. A population size of n is considered. At any iteration, random compression between learners is done. Random selection of two learners is conducted, namely, P and Q. For these two randomly-selected learners X P,i = X Q,i . Updated values of the first phase: X j,P,i is accepted if it provides a better function value. Accepted function values at the end of the learner phase are reserved and treated as the input of the next iteration.
The steps involved in programming for TLBO implementation in presented paper and its flow chart are given in below section.

TLBO Steps Adopted to Optimize Multi-Objective DSM Function F(X)
In this section step needed to adopt TLBO to optimize F(X) is presented first and followed by its flow chart illustration as given in Figure 7.
Step1: Initialize population, design variables (X) of DSM program, and termination criterion: Step2: Calculate the mean of each design variable: In the above equation X(s) mean the result while "r": random number, ∈ [0, 1]. The value of T f determined randomly as: Step3: Determine the fittest solution.
Step4: Improve the solution according to the fittest solution: Step5: Determine the better solution between X i and X i .
Step5(a): If the above condition is not satisfied then reject the new result and keep the previous one and go to the next step.
Step5(b): If the condition is satisfied then accept the new one by replacing the old and go to step6.
Step6: Choose two random solution X i and X k .
Step7: Determine fittest solution between X i and X k . Step7(a): If X i is better than X k then perform the update below and go to Step8: Step7(b): If X i not better than X k than perform the below update and go to the next step.
Step8: Is solution X i better than X i ? Step8(a): If the condition is not satisfied then reject X i , keep X i and go to step9.
Step8(b): If the condition is satisfied then accept X i , replace X i and go to the next step.
Step9: Is the termination criteria satisfied? If the answer is yes then go to the next step; otherwise start the whole process again from step2.
Step10: Report the optimum solution. Step5(a) : If the above condition is not satisfied then reject the new result and keep the previous one and go to the next step.
Step5(b): If the condition is satisfied then accept the new one by replacing the old and go to step6.
Step7(a): If ′ is better than ′ then perform the update below and go to Step8: Step7(b): If ′ not better than ′ than perform the below update and go to the next step.
Step8(a): If the condition is not satisfied then reject " , keep ′ and go to step9.
Step8(b): If the condition is satisfied then accept " , replace ′ and go to the next step.
Step9: Is the termination criteria satisfied? If the answer is yes then go to the next step; otherwise start the whole process again from step2.
Step10: Report the optimum solution.

Object-Oriented Usability Indices (OOUI)
Under this section the proposed indices are defined and discussed.

Object-Oriented Usability Indices (OOUI)
Under this section the proposed indices are defined and discussed. Peak load shaving and increased integration of renewable sources were identified as the main objectives of a demand-side management program. A usability index is expected to present the fulfillment of these two objectives in any specific case study via quantified numerical factors. The outcome of any demand response program does not contain any such indices directly based on which the relative fulfillment of these objectives could be identified. In the absence of such indices, the policy decision of going with demand-side management may not produce the optimized result all the time. Therefore, OOUI are proposed in this section to indicate the relative fulfillment of the two main objectives.
In a power network, the electricity demand varies with time which is presented through the load profile. The power operator faces a challenge of balancing load and generation when peak load conditions appear. Demand-side management, when appropriately performed, reduces some load in peak hours and relieves the operator's stress. Load factor, which is the ratio of the average load on the generator over a period of time to the peak load in the same time interval, indicates the variability of the load. A low load factor indicates that the load is highly variable. A high load factor is desirable for the economic feasibility of plant.
If a demand response program performs as per the requirements, it should reduce the load in peak hours and, therefore, the load factor should improve after allying DSM. A peak power shaving factor is defined as: where LF DSM is the load factor of the power network when DSM is applied and LF WDSM is the load factor of the power network when it operates without demand response. If DSM successfully shaves the load in peak hours, then the load factor will improve and the peak power shaving index will be greater. For better peak power shaving operation, m pps should be as high as possible.
Renewable energy generators are often small in size and unpredictable in terms of availability. Therefore, the generation support from these sources may or may not always be helpful for the power operator. A good demand response program, in terms of renewable energy integration, is one which ensures a good amount of load satisfied with renewable energy sources. A renewable energy integration index indicates the percentage utilization of electrical energy from renewable-based sources. This index is given as: where P gr active power is supplied by renewable energy sources and P gt is the total active power supplied. An overall usability index which indicates the fulfillment of both the objectives of DSM presented by OOUI. Depending upon the specific system architecture and policy-related decisions, the relative weighting of these two indices may vary. Therefore, the overall usability index also includes the individual weighting and is defined as: m dsmu f = w pps m pps + w rei m rei ; where w pps + w rei = 1; Above, w pps and w rei are the relative weighting of the peak power shaving index and renewable energy integration index, respectively.

Theoretical Analysis of the Proposed Indices
This section presents the theoretical analysis of proposed indices. First of all, the analysis of the peak power shaving index is presented, and then the renewable energy integration index is analyzed. Finally, the overall usability index is analyzed. Let load demand to utility over a time interval T is given as: where d 1 , d 2 , d 3 , . . . . . . , d T , are loads during 1st, 2nd, 3rd . . . and Tth, dispatch intervals (in the presented work hourly load variation with dispatch intervals of one hour are considered (Figure 3a). Average load demand to utility over time period T is calculated using the expression: Now the load factor with DSM (LF DSM ) and without DSM (LF WDSM ) is calculated using the following expression: where: t: dispatch interval. d t : load demand to utility when no DSM is applied. d max : maximum load over the time period T when no DSM is applied d t : load demand to utility when DSM is applied. d max : maximum load over the time period T when DSM is applied m pps is defined as the ratio of LF DSM /LF WDSM . This can be represented as: Any index proposed as an indicator of DSM usability should quantify the obtained benefit. One of the main purposes of DSM is to reduce the stress on the operator. The sudden high demand from the load in some hours of operation puts stress over the system. Therefore, the index which shows the benefit in terms of reduction of load should be an indicator.
Peak power index analysis can be done under the following points: 1. Peak power shaving index greater than one that is m pps > 1.

2.
Peak power shaving index less than one that is m pps < 1.

3.
Peak power shaving index equal one that is m pps = 1.

4.
Comparative result of peak power shaving index obtained for DSM solutions.
The first condition, m pps > 1, will arise if the load factor with DSM (LF DSM ) is higher than the load factor without DSM (LF WDSM ). In this condition the value of the peak power shaving index will be higher. Here, the value of m pps indicates how much improvement is achieved in terms of load factor after implementing DSM.
The second condition, m pps < 1, will arise if the load factor with DSM (LF DSM ) is lower than the load factor without DSM (LF WDSM ). In this condition the value of the peak power shaving index will be lower than 1. The value of m pps in this condition indicates a decrease in the load factor value with respect to the without-DSM case.
The third condition, m pps = 1, will arise if the load factor with DSM (LF DSM ) is equal to the load factor without DSM (LF WDSM ). This also indicates that the load factor without DSM is equal.
If the peak power shaving index is calculated and presented for different solutions of DSM, then its relative higher value indicates a better load factor compared to the DSM outcome in other cases. Coincidently, if all DSM solutions are equally good considering load reduction, then the higher value of the peak power shaving index in a specific DSM program can give a ranking of achieved relative peak shaving.
The following example gives an illustrative view about the application of the peak power shaving index. Let us assume that there is a DSM scheme, F(X), for a grid-tied microgrid. This DSM is adopted for five different availability conditions of renewable energy. These conditions are named as case 1, case 2, case 3, case 4, and case 5. The proposed peak power shaving index is calculated for each case and their values are presented as 'm pps '.
If the operator wants to know that, among all five cases, in which case maximum peak shaving has been achieved, the answer can be provided by the values of m pps . To demonstrate the working of m pps , a bar chart is plotted. Indicator values as given in Figure 8 represent that in case 3 the peak power shaving index hold maximum value which indicate that the percentage contribution of peak power shaving is the maximum in case 3 relative to other DSM cases. Additionally, the higher value of case 3 indicates that the load factor in this case with respect to other cases is much better. is adopted for five different availability conditions of renewable energy. These conditions are named as case1, case2, case3, case4, and case5. The proposed peak power shaving index is calculated for each case and their values are presented as ' '. If the operator wants to know that, among all five cases, in which case maximum peak shaving has been achieved, the answer can be provided by the values of . To demonstrate the working of , a bar chart is plotted. Indicator values as given in Figure 8 represent that in case 3 the peak power shaving index hold maximum value which indicate that the percentage contribution of peak power shaving is the maximum in case 3relative to other DSM cases. Additionally, the higher value of case 3 indicates that the load factor in this case with respect to other cases is much better. The ratio of active power from renewables ( ) to total active power ( ) quantifies the involvement of renewable energy. Suppose if the value of this index is 0.26, then it means renewable energy contribution in this case to fulfill total demand is 26%. The lowest possible value for this index can be zero, which will indicate the condition of no-contribution of renewable energy in fulfilling the demand. If the value of is equal to one than it mean 100% renewable energy is used to fulfill all the demand.
If the operator wants to know that, among all cases, which one achieves relatively major integration of renewable energy sources while implementing DSM, can help here by providing a numerical value in the form of an index. The DSM case with a higher value of involves a greater renewable energy percentage to fulfill load demand.
An overall usability index indicates DSM usability combined with peak load shaving and renewable energy integration. Depending upon the specific system architecture and policy-related decisions, the individual weighting of and may vary. In extreme cases, if the weight associated with is set at 0and the weight associated is set at 1, then the overall usability index (represented by ) becomes equal to the peak power shaving index. Thus, it indicates that the desired aim of DSM is only to take care of peak power shaving. Similarly, if this extreme setting is interchanged for DSM, it can be concluded that the DSM scheme is only focusing on renewable energy integration and has zero interest in achieving load factor improvement or peak power shaving.

Case Studies, Simulation Results, and Discussion
Under this section we first present the different renewable energy availability cases. Simulation results of the multi-objective formulation of the DSM using the TLBO method in the MATLAB environment is presented next. Finally, determination of OOUI, followed by a discussion, is conducted. Figure 9 represents the approach adopted in the presented work with respect to the application and usability of the proposed indices. First of all, a microgrid multi-objective DSM formulation ( ( )) is conducted. In the next step the TLBO algorithm is adopted to solve ( ), and the optimized results are obtained. From these solutions the proposed OOUI are calculated and made available to the The ratio of active power from renewables (P gr ) to total active power P gt quantifies the involvement of renewable energy. Suppose if the value of this index is 0.26, then it means renewable energy contribution in this case to fulfill total demand is 26%. The lowest possible value for this index can be zero, which will indicate the condition of no-contribution of renewable energy in fulfilling the demand. If the value of m rei is equal to one than it mean 100% renewable energy is used to fulfill all the demand.
If the operator wants to know that, among all cases, which one achieves relatively major integration of renewable energy sources while implementing DSM, m rei can help here by providing a numerical value in the form of an index. The DSM case with a higher value of m rei involves a greater renewable energy percentage to fulfill load demand.
An overall usability index indicates DSM usability combined with peak load shaving and renewable energy integration. Depending upon the specific system architecture and policy-related decisions, the individual weighting of m pps and m rei may vary.
In extreme cases, if the weight associated with m rei is set at 0 and the weight associated m pps is set at 1, then the overall usability index (represented by m dsmu f ) becomes equal to the peak power shaving index. Thus, it indicates that the desired aim of DSM is only to take care of peak power shaving. Similarly, if this extreme setting is interchanged for DSM, it can be concluded that the DSM scheme is only focusing on renewable energy integration and has zero interest in achieving load factor improvement or peak power shaving.

Case Studies, Simulation Results, and Discussion
Under this section we first present the different renewable energy availability cases. Simulation results of the multi-objective formulation of the DSM using the TLBO method in the MATLAB environment is presented next. Finally, determination of OOUI, followed by a discussion, is conducted. Figure 9 represents the approach adopted in the presented work with respect to the application and usability of the proposed indices. First of all, a microgrid multi-objective DSM formulation (F(X)) is conducted. In the next step the TLBO algorithm is adopted to solve F(X), and the optimized results are obtained. From these solutions the proposed OOUI are calculated and made available to the operator. A MATLAB-based simulation has been performed to obtain all of the simulation results.
It is considered that the participation of all loads is equally shared, and each type of load has equal right to make a decision to get involved in the demand response program. The microgrid operator has the information about the interruptible energy limit of the user on per day bases and this information is used to index the end users in order of increasing interest to diminish their energy need. Additionally, the outage cost function coefficients (K 1,i and K 2,i ) of involved users are considered to be known to the operator. The daily utility budget is $500. It is considered that the participation of all loads is equally shared, and each type of load has equal right to make a decision to get involved in the demand response program. The microgrid operator has the information about the interruptible energy limit of the user on per day bases and this information is used to index the end users in order of increasing interest to diminish their energy need. Additionally, the outage cost function coefficients ( , and , ) of involved users are considered to be known to the operator. The daily utility budget is $500.

Case Studies, Simulation Results.
Solution of the multi-objective formulation of DSM is obtained using TLBO for five different scenarios of renewable energy availability in the microgrid. These scenarios are presented as case 1-5. Case 1 deals with the condition of average wind and average solar availability; in this case total availability of renewable power is 346.24 kW. Case 2 considers the maximum wind and average solar availability and total renewable power availability in this case is 400.5 kW. Case 3 is when wind availability is at a minimum, and solar availability is average and overall availability of renewable sources is 275.03 kW. Average wind and maximum solar availability are termed as case 4; in this case renewable energy availability is 362.76 kW. Finally, case 5 deals with the generation scenario when the wind is average and solar has minimum availability, and the total availability of renewable energy is 296.01 kW. The solar power availability and wind power availability under different cases are represented in Figures 10 and 11; these figures present the hourly variation of solar and wind energy availability in different cases.
Under all of these conditions the DSM program is run, and the iteration vs. fitness value plots are given in Figure 12. The results indicate that the optimization target is achieved around 2000 iterations in all cases. As a result of the optimization new load profiles are obtained for each case. These new load profiles are in accordance with the two objectives of F(X). These load curves for different cases are plotted and presented in Figure 13, with and without the DSM load curve. Values of parameters obtained from these load curves are listed in Table 3. Total, average, and peak demand when DSM is not implemented is 864.9 kW, 36.0375 kW, and 42.09 kW, and the load factor is 0.84. These parameters are needed to calculate OOUI values.

Case Studies, Simulation Results
Solution of the multi-objective formulation of DSM is obtained using TLBO for five different scenarios of renewable energy availability in the microgrid. These scenarios are presented as cases 1-5. Case 1 deals with the condition of average wind and average solar availability; in this case total availability of renewable power is 346.24 kW. Case 2 considers the maximum wind and average solar availability and total renewable power availability in this case is 400.5 kW. Case 3 is when wind availability is at a minimum, and solar availability is average and overall availability of renewable sources is 275.03 kW. Average wind and maximum solar availability are termed as case 4; in this case renewable energy availability is 362.76 kW. Finally, case 5 deals with the generation scenario when the wind is average and solar has minimum availability, and the total availability of renewable energy is 296.01 kW. The solar power availability and wind power availability under different cases are represented in Figures 10 and 11; these figures present the hourly variation of solar and wind energy availability in different cases.
Under all of these conditions the DSM program is run, and the iteration vs. fitness value plots are given in Figure 12. The results indicate that the optimization target is achieved around 2000 iterations in all cases. As a result of the optimization new load profiles are obtained for each case. These new load profiles are in accordance with the two objectives of F(X). These load curves for different cases are plotted and presented in Figure 13, with and without the DSM load curve. Values of parameters obtained from these load curves are listed in Table 3. Total, average, and peak demand when DSM is not implemented is 864.9 kW, 36.0375 kW, and 42.09 kW, and the load factor is 0.84. These parameters are needed to calculate OOUI values.

Calaulation of Proposed Indices
The calculations of the proposed indices are based on the mathematical relation as given in Section 5. Table 4 presents the calculated values of these indices under all five cases studied. The renewable energy integration index presents how much the participation of renewable sources has been achieved. The peak power shaving index indicates how much peak load shaving has been achieved. The overall usability index is calculated with the help of the weighted sum of these two indices.
This work considers three different weighting scenarios named as scenario 1, scenario 2, and scenario 3. In scenario 1, equal weighting w pps = w rei = 0.5 is given to both objectives. In the second scenario, the renewable energy integration has been given higher weighting (w rei = 0.9). In the third scenario, the peak power shaving remains the preferred objective and has higher weighting (w pps = 0.9). Table 4 presents the calculated indices and overall index under all five cases of energy availability and weighting scenarios. The peak power shaving targets are achieved with the maximum amount in case 1 when wind and solar energy availabilities are average.  The next best case is case 4 when wind availability is average and solar energy is at maximum availability. This indicates that a significant amount of peak power shaving can be achieved when at least an average solar availability is ensured. The reason can be attributed to the solar availability in peak load duration. With average or maximum availability, solar provides additional support to conventional generators and the grid in easing the peak load conditions. Renewable energy utilization is another objective of DSM. The renewable energy integration index indicates the fractional quantity of renewable generation concerning a conventional generator. The calculated values indicate that maximum integration of renewable energy is achieved with DSM in case number 4, where wind availability is average and solar availability is the maximum. In general, the condition where both have maximum availability should be the best scenario for renewable energy utilization.
For such condition to occur, a combination of suitable weather conditions is required. Therefore, the next best situation is when one source generates the maximum and the other generates average energy. These two conditions occur in case 2 and case 4. However, as the values of the renewable energy integration index indicate, solar energy has more impact on fulfilling this target. When both the indices are combined to obtain a final index, which presents the object-oriented DSM usability, the condition for the maximum benefit of DSM is obtained. As can be seen from Table 3, the overall maximum benefit from DSM can be achieved in case 3, when wind provides minimum energy and solar provides average energy. This conclusion remains valid for any weighting given to either objective of DSM. Whether either objectives have equal weighting, or one of them has higher weighting, the best result can be obtained only in case number 3. This conclusion can help in policy-making while deciding to adopt the DSM for any specific system. Figures 14 and 15 present the different indices in plots where the conclusion can be easily identified and established. weighting, the best result can be obtained only in case number 3. This conclusion can help in policymaking while deciding to adopt the DSM for any specific system. Figures 14 and 15 present the different indices in plots where the conclusion can be easily identified and established.   Table 4.

Discussion of the Results with Comparative Analysis
This section is divided into the following two parts:(A) The findings of the present work are discussed keeping the focus on the obtained values of the proposed indices for DSM; and (B) comparative discussion with similar studies published by other authors to show the uniqueness of the present work.
(A) Observation of Tables 3 and 4, focusing on the numeric values of peak power shaving index, reveals two things: (1) regarding peak shaving and (2) the load factor. A close look of Tables 3 and 4 shows that the peak shaving actor holds adirect relation with the load factor and a reciprocal relation with peak shaving. Observation of Table5 reveals that case 1 holds rank 1 when ranking is done on the basis of the peak power shaving index. This indicates that when DSM is applied to case1 then it achieves maximum peak shaving in obtaining its objectives. Case 2 is ranked 5th on the basis of the peak shaving index value, which means that peak reduction is the minimum in this case. In case 2 the maximum peaks are 41.67 kW (Table 3) which is very close to without DSM load curve peak (42.09 kW). Figure 13 (case2) shows that in this case DSM also reduces the overall demand and achieves its objectives, but is poor in terms of peak shaving. In terms of higher peak reduction, cases 1-5 can be arranged as follows: case1>case5>case3>case4>case2. weighting, the best result can be obtained only in case number 3. This conclusion can help in policymaking while deciding to adopt the DSM for any specific system. Figures 14 and 15 present the different indices in plots where the conclusion can be easily identified and established.   Table 4.

Discussion of the Results with Comparative Analysis
This section is divided into the following two parts:(A) The findings of the present work are discussed keeping the focus on the obtained values of the proposed indices for DSM; and (B) comparative discussion with similar studies published by other authors to show the uniqueness of the present work.
(A) Observation of Tables 3 and 4, focusing on the numeric values of peak power shaving index, reveals two things: (1) regarding peak shaving and (2) the load factor. A close look of Tables 3 and 4 shows that the peak shaving actor holds adirect relation with the load factor and a reciprocal relation with peak shaving. Observation of Table5 reveals that case 1 holds rank 1 when ranking is done on the basis of the peak power shaving index. This indicates that when DSM is applied to case1 then it achieves maximum peak shaving in obtaining its objectives. Case 2 is ranked 5th on the basis of the peak shaving index value, which means that peak reduction is the minimum in this case. In case 2 the maximum peaks are 41.67 kW (Table 3) which is very close to without DSM load curve peak (42.09 kW). Figure 13 (case2) shows that in this case DSM also reduces the overall demand and achieves its objectives, but is poor in terms of peak shaving. In terms of higher peak reduction, cases 1-5 can be arranged as follows: case1>case5>case3>case4>case2.  Table 4.

Discussion of the Results with Comparative Analysis
This section is divided into the following two parts: (A) The findings of the present work are discussed keeping the focus on the obtained values of the proposed indices for DSM; and (B) comparative discussion with similar studies published by other authors to show the uniqueness of the present work.
(A) Observation of Tables 3 and 4, focusing on the numeric values of peak power shaving index, reveals two things: (1) regarding peak shaving and (2) the load factor. A close look of Tables 3 and 4 shows that the peak shaving actor holds adirect relation with the load factor and a reciprocal relation with peak shaving. Observation of Table 5 reveals that case 1 holds rank 1 when ranking is done on the basis of the peak power shaving index. This indicates that when DSM is applied to case 1 then it achieves maximum peak shaving in obtaining its objectives. Case 2 is ranked 5th on the basis of the peak shaving index value, which means that peak reduction is the minimum in this case. In case 2 the maximum peaks are 41.67 kW (Table 3) which is very close to without DSM load curve peak (42.09 kW). Figure 13 (case 2) shows that in this case DSM also reduces the overall demand and achieves its objectives, but is poor in terms of peak shaving. In terms of higher peak reduction, cases 1-5 can be arranged as follows: case 1 > case 5 > case 3 > case 4 > case 2.
Observation of Tables 3 and 4, while focusing on the numeric values of renewable energy index, reveals that its value is the maximum in case 2 and the minimum in case 3. This indicates that the percentage contribution of renewable energy in fulfilling the total load demand is the maximum with respect to other cases. In case 2 total power requirements for the users is 793.92 kW and the renewable energy integration index value is 0.504, which indicates that the contribution of renewable energy in the total requirement is 50.4%. Cases 1-5 can be arranged according to renewable energy contribution in their total demand as follows: case 2 > case 4 > case 1 > case 5 > case 3.
Observation of Table 4, while focusing on the numeric values obtained for the overall usability index, giving equal weighting (0.5) to peak power shaving and the renewable energy integration index shows that it holds a maximum value for case 2 and a minimum for case 3.
The overall usability index indicates the relative effectiveness of DSM in reducing peak power and utilizing renewable energy simultaneously while achieving its objectives. Observation of Tables 4  and 5 for m s1 dsmu f shows that case 2 has rank 1 as in this case the combined impact of peak power shaving and renewable contribution is the highest (Tables 3 and 4, verified this). Table 5 shows that ranking of the overall usability index followed the pattern of m rei ranking when high weighting (w rei = 0.9) is associated with the renewable energy integration index ( Table 4). The overall usability index tends to follow the pattern of m pps ranking when weighting (w pps = 0.9) shifted toward m pps . Table 5 shows the ranking based on different criteria.
(B) Previous works have also tried to present a framework or factor for DSM. In a study, an impact analysis of DSM is presented through the Inverted Load Duration Curve (ILDC) [41]. This study remodels the power network by taking the Inverted Load Duration Curve. DSM is applied for the renewed model and the DSM impact is presented through a model, which is named as "VPI model". This presents the benefit of DSM in a quantified manner in corresponding to only the load, whereas this work presents the benefit quantification corresponding to both renewable generation and loads.
As discussed in the introduction section, similar indices are proposed in [30], wherein the economic and technical benefits are enumerated and quantified and two indices are proposed based on current reduction and cost of one kWh. However, the technical benefits take in account only the reduction in loads and do not consider the effect of different types of generation.
The VPI model as proposed in [41] is further extended and analyzed for a big power network in [42]. The influence of DSM is calculated and evaluated based on load reduction only. The impact of renewable energy generation and small DG's are not discussed. Additionally, the quantified value is not presented in terms of any general index.
In [43], a set of performance indicators are proposed to design a new energy management system. The focus of these indicators is consumption allocation to individual load, total consumption, cost allocation, total cost, maximum consumption, energy usage, and hourly distribution of energy. The systems under study in this paper were Supervisory Control and Data Acquisition (SCADA)-based connected loads and generation units. The formulation of DSM is not discussed, and the mathematical formulation of indicators is not presented through any simulation results. In another work, a decision support framework is presented for selection of the demand response method [44]. The proposed framework of [44] is based on cost minimization and load scheduling. The work, however, does not consider any role of renewable and focuses only on energy management. Table 6 shows the uniqueness of the present work with respect to studies done by other authors [30,[41][42][43][44] in a comparative manner.

Conclusions
The main finding of this paper is to provide the OOUI, namely, peak power shaving index, renewable energy integration index, and overall usability index. OOUI allows quantifying of the technical and economic benefits from DSM. The formulated indices are helpful for the operator in identifying the relatively more suitable operating condition, and relatively more beneficial, for DSM, which can help in policy-related decisions. The economic benefit of demand response can also be quantified using the proposed index mechanism.
The concluding numeric values of the indices obtained are listed in Table 4. The ranking of DSM, in Table 5, is done on the basis of these values. Observation of Table 5 gives clear information that the peak power shaving index value is the maximum in case 1. This indicates that the DSM application in case 1 results in the maximum peak shaving. Similarly, the renewable energy integration indicator value is the maximum in case 2. This indicates that the renewable energy contribution to fulfilling load demand is the maximum in case 2. Based on the two indices as described above, an overall usability index is designed. It is observed from results that the overall usability index has the maximum value in case 2 when other two indices are given equal weight.
The overall importance of this paper is that it provides a mathematical tool for OOUI. These indices provide additional aid to the operator in identifying the conditions in which DSM provides more benefit in terms of peak power shaving and renewable energy integration. Thus, it can help the operator in policy-related decisions. In other existing works, the quantification of the DSM benefit is made either through a load curve or through generation reduction. However, these existing works have not considered both load change and renewable energy integration simultaneously. The presented simulation-based analysis establishes that the proposed OOUI can successfully quantify the two important benefits, curtailments of peak load and harnessing the renewable energy. Additionally, the two benefits can be clubbed together and presented as a single overall usability index.
The proposed approach shows that OOUI have the potential to determine a DSM solution, which can provide relief to both the utility and customers by promoting higher renewable energy penetration and simultaneously shaving the peak load. A limitation of the current approach is that OOUI are obtained and analyzed for the DSM solutions obtained for multi-objective formulation, F(X), using single objective optimization algorithm TLBO. In future work, this index-based approach can be analyzed for different multiobjective DSM models incorporating a multi-objective optimization algorithm. These DSM models can be based on different demand response programs, such as an emergency demand response program, and critical peak pricing.