Implementation of a Small Type DC Microgrid Based on Fuzzy Control and Dynamic Programming

Abstract: A DC microgrid (DC-MG) is a novel power system that uses DC distribution in order to provide high quality power. The study system is made by a photovoltaic array (PV), a wind generator (WG), a fuel cell (FC), and an energy storage system (ESS) to establish a small type DC microgrid, with the bus being established by DC/DC converters with fuzzy controllers. An overall power dispatch was designed for the proposed system to distribute the power flows among the different energy sources and the storage unit in the system in order to satisfy the load requirements throughout an entire 24-h period. The structure of a power supervisor based on an optimal power dispatch algorithm is here proposed. Optimization was performed using dynamic programming (DP). In this paper, a system configuration of a DC microgrid is analyzed in different scenarios to show the efficacy of the control for all devices for the variable weather conditions with different DC loads. Thus, the voltage level and the power flow of the system are shown for different load conditions.


Introduction
Global environmental concerns and the ever-increasing need for energy, coupled with steady progress in renewable energy technologies, are opening up new opportunities for the utilization of renewable energy resources.Developments of wind and photovoltaic array (PV) energy technologies have increased their use in hybrid wind/PV configurations.Further, integrating of PV and wind energy sources with fuel cells leads to a non-polluting reliable energy source [1].
This paper investigates power flow management for a DC microgrid (DC-MG) with an energy storage system focused on optimal scheduling.The method for performing power management optimization is chosen according to the nature of the problem (components, constraints, and performance index).The main purpose of this work is to find the power flow dispatch that controls for all units through specific converters, especially in case of variable weather conditions with different loads over the studied period.A fuzzy control system tuned by a dynamic programming approach is used to manage the microgrid system.
For ease of implementation, fuzzy controllers with microprocessors have been successfully applied to power converters [2][3][4] and inverters [5][6][7][8].The design of a fuzzy controller is mainly related to knowledge acquisition, tuning, and the improvement of the optimality and stability of the control system.Tuning the scaling factors, the membership function, and the set of control rules have been proposed recently for tuning the fuzzy controller [9].The membership function parameters of input variables and output variables in a fuzzy control system are commonly used to conduct proper transformations between the real input data and the universe of the discourses of the fuzzy variables in the system.Sometimes, membership function parameters can be used to fine-tune the performance of the designed system in a similar way to tune a proportional-integral-derivative (PID) controller for a better dynamic response [10][11][12].Thus, fuzzy controllers are used for voltage regulation of the system for all the converters used in the microgrid.A dynamic programming (DP) approach to a problem discretizes the problem's state space and uses a backwards recursion to obtain the optimal value in every state.The state variable describing the problem at a particular stage is defined that it completely describes the process.Given the state of the process at the beginning of a stage, we made a decision that transforms the process to the final state at the end of the stage.The objective was to maximize (or minimize) the expected objective over all stages.DP is a graph-based technique corresponding to the shorter path algorithms.The performance index and the constraints can hold all the natures, and no specific mathematical solver is needed are the advantages of the method.Its high memory needs when the studied period is long and discretized within a small time frame is the weakness of this technique.However, it is not problematic if the computation parameters are well chosen, and the computation time can be reduced by appropriate modifications [13,14].

Microgrid Architecture
The studied DC-MG consists of a cluster of generators, a storage system, a load, and power electronic interfaces all connected to the DC grid as shown in Figure 1.The regarded distributed generators (DGs) are a 1.2-W wind generator (WG), a 1.8-W PV, and a 1.0-W fuel cell (FC).A stationary lithium iron phosphate battery (3.2 V, 1100 mAh) composes the energy storage system (ESS) and is used for storing and supplying electrical power.The specifications of these devices are shown in Table 1.
The voltage of the DC bus is set to 10 V in this study.A BUCK converter is used for the PV to decrease the voltage level.A SEPIC converter is used for the WG through a rectifier to regulate the DC voltage.A BOOST converter is used for the FC to increase the voltage level.A bidirectional DC/DC converter is used for the ESS to charge and discharge the energy from the DC bus to balance the power flow and stabilize the voltage level.The deviation of power of each DG is obtained by the measured current (i) and the deviation of voltage ( V).Each of the four fuzzy logic controllers (FLCs) has two inputs, the deviation of load power (∆P load ), and the deviation of converter power (∆P DG ).The output of each FLC is the control signal of duty cycle (∆duty).The experiment was executed on a personal computer (PC).The A/D converter was used to convert the analog output voltage of the DC/DC converter into digital data, and a D/A converter was used to convert the digital results calculated by the PC into an analog control signal.The analog control signal was fed into ICs TL494 for a pulse width modulation (PWM) to generate the appropriate duty cycle.

Fuzzy Logic Controller (FLC)
There are four steps for the design of FLC: Step 1: Defining the Model's Functional and Operational Characteristics Defining exactly where the fuzzy subsystem fits into the total system architecture, which helps the designer to estimate the numbers and ranges of inputs and outputs that will be required is the essential step.
Step 2: Designing the Fuzzifier Mapping from an observed input space to fuzzy sets in certain input universes of discourse is the fuzzification interface performed by the fuzzifier.This strategy interprets a crisp value x as a fuzzy set A with membership function ( ) μ A x belonging to [0, 1].Each variable in the fuzzy model is divided into a set of fuzzy regions, which are given unique names, with labels such as "positive big", "negative small", "cool", "slow", etc.These labels are related to the physical states of the fuzzy variable.

Fuzzy Logic Controller (FLC)
There are four steps for the design of FLC: Step 1: Defining the Model's Functional and Operational Characteristics Defining exactly where the fuzzy subsystem fits into the total system architecture, which helps the designer to estimate the numbers and ranges of inputs and outputs that will be required is the essential step.
Step 2: Designing the Fuzzifier Mapping from an observed input space to fuzzy sets in certain input universes of discourse is the fuzzification interface performed by the fuzzifier.This strategy interprets a crisp value x as a fuzzy set A with membership function ( ) μ A x belonging to [0, 1].Each variable in the fuzzy model is divided into a set of fuzzy regions, which are given unique names, with labels such as "positive big", "negative small", "cool", "slow", etc.These labels are related to the physical states of the fuzzy variable.

Fuzzy Logic Controller (FLC)
There are four steps for the design of FLC: Step

1: Defining the Model's Functional and Operational Characteristics
Defining exactly where the fuzzy subsystem fits into the total system architecture, which helps the designer to estimate the numbers and ranges of inputs and outputs that will be required is the essential step.
Step 2: Designing the Fuzzifier Mapping from an observed input space to fuzzy sets in certain input universes of discourse is the fuzzification interface performed by the fuzzifier.This strategy interprets a crisp value x as a fuzzy set A with membership function ( ) μ A x belonging to [0, 1].Each variable in the fuzzy model is divided into a set of fuzzy regions, which are given unique names, with labels such as "positive big", "negative small", "cool", "slow", etc.These labels are related to the physical states of the fuzzy variable.

Fuzzy Logic Controller (FLC)
There are four steps for the design of FLC: Step 1: Defining the Model's Functional and Operational Characteristics Defining exactly where the fuzzy subsystem fits into the total system architecture, which helps the designer to estimate the numbers and ranges of inputs and outputs that will be required is the essential step.
Step 2: Designing the Fuzzifier Mapping from an observed input space to fuzzy sets in certain input universes of discourse is the fuzzification interface performed by the fuzzifier.This strategy interprets a crisp value x as a fuzzy set A with membership function ( ) μ A x belonging to [0, 1].Each variable in the fuzzy model is divided into a set of fuzzy regions, which are given unique names, with labels such as "positive big", "negative small", "cool", "slow", etc.These labels are related to the physical states of the fuzzy variable.

Fuzzy Logic Controller (FLC)
There are four steps for the design of FLC: Step 1: Defining the Model's Functional and Operational Characteristics Defining exactly where the fuzzy subsystem fits into the total system architecture, which helps the designer to estimate the numbers and ranges of inputs and outputs that will be required is the essential step.
Step 2: Designing the Fuzzifier Mapping from an observed input space to fuzzy sets in certain input universes of discourse is the fuzzification interface performed by the fuzzifier.This strategy interprets a crisp value x as a fuzzy set A with membership function µ A(x) belonging to [0, 1].Each variable in the fuzzy model is divided into a set of fuzzy regions, which are given unique names, with labels such as "positive big", "negative small", "cool", "slow", etc.These labels are related to the physical states of the fuzzy variable.

Step 3: Designing the Inference Engine
This step involves writing the rules with syntax, e.g., IF <fuzzy proposition> THEN <fuzzy proposition>.The fuzzy proposition might be originally written in English, e.g., IF the change of load power is decreasing very much and the change of DG power is also decreasing very much, THEN the change of duty of the system shall be decreased greatly.Thus, the statement of the rule will be IF ∆P load is NB and ∆P DG is NB THEN ∆duty is NB.There are five methods to generate fuzzy control rules.These methods are not mutually exclusive, and it seems likely that a combination of them may be necessary to construct an effective method for the derivation of fuzzy control rules.These methods are (1) expert experience and knowledge; (2) modeling the operator's control actions; (3) modeling the process; (4) self-organization; and (5) learning from numerical examples.

Step 4: Selecting a Defuzzification Method
The final step of creating the FLC is the selection of the defuzzification interface.There are several ways to convert an output fuzzy set into a crisp variable, but the most commonly used strategies include max criterion, mean of maximum method (MOM), and center of area method (COA).The COA method is used in this paper to generate the center of gravity of the possibility distribution of a control action.The method considers both membership function value and its shape.In the case of two-inputs (x 1 ,x 2 ) and one-output (y), , where N is the total number of fuzzy rules, y i denotes the center value of the output set in their domain, which is fired at the ith fuzzy rule, and µ Ai (µ Bi ) is the grade of membership at which x 1 and x 2 belong to the sets Ai (Bi).
All membership functions (MFs) of each of the FLC inputs and outputs are defined on the common normalized domain [0, 1], as shown in Figures 2-4.The acronyms NB, NS, ZE, PS, and PB stand for negative big, negative small, zero, positive small, and positive big, respectively.Here, triangular MFs are chosen for all fuzzy sets.The fuzzy rules for ∆P load and ∆P DG are shown in Table 2.

Step 3: Designing the Inference Engine
This step involves writing the rules with syntax, e.g., IF <fuzzy proposition> THEN <fuzzy proposition>.The fuzzy proposition might be originally written in English, e.g., IF the change of load power is decreasing very much and the change of DG power is also decreasing very much, THEN the change of duty of the system shall be decreased greatly.Thus, the statement of the rule will be IF  load P is NB and  DG P is NB THEN duty is NB.There are five methods to generate fuzzy control rules.These methods are not mutually exclusive, and it seems likely that a combination of them may be necessary to construct an effective method for the derivation of fuzzy control rules.These methods are (1) expert experience and knowledge; (2) modeling the operatorʹs control actions; (3) modeling the process; (4) self-organization; and (5) learning from numerical examples.
Step 4: Selecting a Defuzzification Method The final step of creating the FLC is the selection of the defuzzification interface.There are several ways to convert an output fuzzy set into a crisp variable, but the most commonly used strategies include max criterion, mean of maximum method (MOM), and center of area method (COA).The COA method is used in this paper to generate the center of gravity of the possibility distribution of a control action.The method considers both membership function value and its shape.In the case of two-inputs ( 1x , 2 x ) and one-output (y), ( , ).
( , ) where N is the total number of fuzzy rules, i y denotes the center value of the output set in their domain, which is fired at the ith fuzzy rule, and  Ai (  Bi ) is the grade of membership at which 1 x and 2 x belong to the sets Ai (Bi).
All membership functions (MFs) of each of the FLC inputs and outputs are defined on the common normalized domain [0, 1], as shown in Figures 2-4.The acronyms NB, NS, ZE, PS, and PB stand for negative big, negative small, zero, positive small, and positive big, respectively.Here, triangular MFs are chosen for all fuzzy sets.The fuzzy rules for  load P and  DG P are shown in Table 2.

Dynamic Programming (DP)
DP was used to solve the optimization problem.The problem was formulated as a system evolution divided into a multi-stage decision process.To illustrate the type of functional equation that arises from an application of the principle of optimality, begin with the simplest case of a deterministic process where the system is described at any time by an M-dimensional vector Assume that an N-stage process are considered to be carried out to maximize some scalar function, ( )  S p of the final state.This function is called the N-stage return.A policy consists of a selection of N transformations, where

Dynamic Programming (DP)
DP was used to solve the optimization problem.The problem was formulated as a system evolution divided into a multi-stage decision process.To illustrate the type of functional equation that arises from an application of the principle of optimality, begin with the simplest case of a deterministic process where the system is described at any time by an M-dimensional vector Assume that an N-stage process are considered to be carried out to maximize some scalar function, ( )  S p of the final state.This function is called the N-stage return.A policy consists of a selection of N transformations, where

Dynamic Programming (DP)
DP was used to solve the optimization problem.The problem was formulated as a system evolution divided into a multi-stage decision process.To illustrate the type of functional equation that arises from an application of the principle of optimality, begin with the simplest case of a deterministic process where the system is described at any time by an M-dimensional vector p = (p 1 , p 2 , ...p M ), constrained to lie within some region D. Let F = {F k }, where k runs over a set which may be finite, enumerable, or continuous, be a set of transformations with the property that p ∈ D implies that F k (p) ∈ D for all k.
Assume that an N-stage process are considered to be carried out to maximize some scalar function, S(p) of the final state.This function is called the N-stage return.A policy consists of a selection of N transformations, where p = (F 1 , F 2 , ..., F N ), yielding successively the states If D is a finite region, if each F k (p) is continuous in p, and if S {p} is a continuous function of p for p ∈ D, it is clear that an optimal power exists.The maximum value of S(p N ), determined by an optimal power, will be a function only of the initial vector p and the number of stages N. Let us then define This is the N-stage return obtained using an optimal power starting from the initial state p.
To derive a functional equation for P f N (p), we employ the principle cited above.Assuming that we choose some transformation F k as a result of our first decision, thereby obtaining a new state F k (p), the maximum return from the following (N-1) stages is, by definition, P f N−1 (F k (p)).It follows that k must now be chosen so as to maximize this result, the basic functional equation being It is clear that a knowledge of any particular optimal power will yield P f N (p), which is unique.
Conversely, given the sequence P f N (p) , all optimal powers may be determined [15].The constraints of the generated power of DGs are as follows: where P PV is the generated power by PV, P WG is the generated power by WG, and P FC is the generated power by FC.
There are i (i = 0-4) stages between time (t − 1) and time t.The stages relation table is shown in Table 3.The DP models are described by the following equations.
where P f 3 (s) is only determined by FC.P f 2 (s) is determined by WG and P f 3 (s).P f 1 (s) is determined by PV and P f 2 (s).
The detailed computing procedures are as follows: P f 3 (4) = P FC (x 4 ) .( 14)  Case B: At 0:00, there is no wind, the load demand is supplied by the FC, and the extra energy is stored in the ESS.At 8:00, the load demand increases, causing the PV and the WG to supply the load demand simultaneously.At 12:00, the load demand decreases and the ESS is not full, and some of energy is stored in the ESS.At 14:00, the load demand increases, and the PV and the WG are out of service, causing the output of the FC increases.At 19:00, the load demand increases, the FC cannot satisfy the load demand, causing the ESS to provide the output power.At 20:00, the load demand decreases, causing the FC to supply the load demand and the ESS to stop providing the output power.The voltage level and the power flow in the MG are shown in Figures 7 and 8.   Case B: At 0:00, there is no wind, the load demand is supplied by the FC, and the extra energy is stored in the ESS.At 8:00, the load demand increases, causing the PV and the WG to supply the load demand simultaneously.At 12:00, the load demand decreases and the ESS is not full, and some of energy is stored in the ESS.At 14:00, the load demand increases, and the PV and the WG are out of service, causing the output of the FC increases.At 19:00, the load demand increases, the FC cannot satisfy the load demand, causing the ESS to provide the output power.At 20:00, the load demand decreases, causing the FC to supply the load demand and the ESS to stop providing the output power.The voltage level and the power flow in the MG are shown in Figures 7 and 8. Case B: At 0:00, there is no wind, the load demand is supplied by the FC, and the extra energy is stored in the ESS.At 8:00, the load demand increases, causing the PV and the WG to supply the load demand simultaneously.At 12:00, the load demand decreases and the ESS is not full, and some of energy is stored in the ESS.At 14:00, the load demand increases, and the PV and the WG are out of service, causing the output of the FC increases.At 19:00, the load demand increases, the FC cannot satisfy the load demand, causing the ESS to provide the output power.At 20:00, the load demand decreases, causing the FC to supply the load demand and the ESS to stop providing the output power.The voltage level and the power flow in the MG are shown in Figures 7 and 8.  Case C: At 0:00, due to the lower load demand, the load demand is sufficiently supplied by the FC.At 8:00, the PV becomes the main source of energy due to the sun light.At time 9:00, the WG also supplies the energy.Since the PV is the main source, the WG is the slave source.Extra energy is used to charge the ESS.At 12:00, the load demand decreases, causing the charge rate of the ESS to increase.At 14:00, the load demand increases, causing the used power to charge, the ESS is switched to the load to balance the power.At 15:00, the ESS is full.At 18:00, the PV cannot supply the power due to a lack of sunlight, and the WG becomes the main power.At 20:00, the WG cannot supply the power due to a lack of wind.The FC becomes the main power.The shortage of power, causes the stored power in the ESS to switch to discharge power to supply load demand.The voltage level and the power flow in the MG are shown in Figures 9 and 10    Case C: At 0:00, due to the lower load demand, the load demand is sufficiently supplied by the FC.At 8:00, the PV becomes the main source of energy due to the sun light.At time 9:00, the WG also supplies the energy.Since the PV is the main source, the WG is the slave source.Extra energy is used to charge the ESS.At 12:00, the load demand decreases, causing the charge rate of the ESS to increase.At 14:00, the load demand increases, causing the used power to charge, the ESS is switched to the load to balance the power.At 15:00, the ESS is full.At 18:00, the PV cannot supply the power due to a lack of sunlight, and the WG becomes the main power.At 20:00, the WG cannot supply the power due to a lack of wind.The FC becomes the main power.The shortage of power, causes the stored power in the ESS to switch to discharge power to supply load demand.The voltage level and the power flow in the MG are shown in Figures 9 and 10.Case C: At 0:00, due to the lower load demand, the load demand is sufficiently supplied by the FC.At 8:00, the PV becomes the main source of energy due to the sun light.At time 9:00, the WG also supplies the energy.Since the PV is the main source, the WG is the slave source.Extra energy is used to charge the ESS.At 12:00, the load demand decreases, causing the charge rate of the ESS to increase.At 14:00, the load demand increases, causing the used power to charge, the ESS is switched to the load to balance the power.At 15:00, the ESS is full.At 18:00, the PV cannot supply the power due to a lack of sunlight, and the WG becomes the main power.At 20:00, the WG cannot supply the power due to a lack of wind.The FC becomes the main power.The shortage of power, causes the stored power in the ESS to switch to discharge power to supply load demand.The voltage level and the power flow in the MG are shown in Figures 9 and 10.

Conclusions
In this paper, we propose the concept of a small type DC microgrid and show a system configuration and control methods for power converters and generators.A DC-MG-based power generation system with a PV, WT, GE, and ESS generation that is connected to the grid is proposed.During islanded operation, the power and voltage of a microgrid may change rapidly due to a power imbalance between supply and demand.Fuzzy controllers and a dynamic programming method are suggested in order to increase the injected power and balance the power flow and the voltage level of the distributed generators.The energy storage system is also considered to overcome the power flow issue and stabilize the voltage level.The experimental results show that the proposed control method provides good operation performance for different load conditions.

Conclusions
In this paper, we propose the concept of a small type DC microgrid and show a system configuration and control methods for power converters and generators.A DC-MG-based power generation system with a PV, WT, GE, and ESS generation that is connected to the grid is proposed.During islanded operation, the power and voltage of a microgrid may change rapidly due to a power imbalance between supply and demand.Fuzzy controllers and a dynamic programming method are suggested in order to increase the injected power and balance the power flow and the voltage level of the distributed generators.The energy storage system is also considered to overcome the power flow issue and stabilize the voltage level.The experimental results show that the proposed control method provides good operation performance for different load conditions.

PV BP365, 12
-V, 65-W PV Module (When exposed to sunlight, it will generate 12 V of DC power into a battery, up to a maximum of 65 W according to the strength of solar radiation.)WG Rutland WG913 Windcharger (1 × Main Generator Assembly, 6 × Aerofoil Blades).The 12 V of DC power generated relates directly to the wind speed at the turbine.FC/cell Size (L × W × H: 65 × 85× 38 mm), maximum power of 0.25 W, maximum current of 1 A, maximum voltage of 0.9 V. ESS IRF18650P, Lithium iron phosphate battery (3.2 V, 1100 mAh).

BP365, 12 Table 1 .
-V, 65-W PV Module (When exposed to sunlight, it will generate 12 V of DC power into a battery, up to a maximum of 65 W according to the strength of solar radiation.)WG Specifications of the devices.PV BP365, 12-V, 65-W PV Module (When exposed to sunlight, it will generate 12 V of DC power into a battery, up to a maximum of 65 W according to the strength of solar radiation.)WG Rutland WG913 Windcharger (1 × Main Generator Assembly, 6 × Aerofoil Blades).The 12 V of DC power generated relates directly to the wind speed at the turbine.FC/cell Size (L × W × H: 65 × 85× 38 mm), maximum power of 0.25 W, maximum current of 1 A, maximum voltage of 0.9 V. ESS IRF18650P, Lithium iron phosphate battery (3.2 V, 1100 mAh).

Rutland WG913 Windcharger ( 1 × 12 Table 1 .
Main Generator Assembly, 6 × Aerofoil Blades).The 12 V of DC power generated relates directly to the wind speed at the turbine.FC/cell Energies 2016, 9, 781 4 of Specifications of the devices.PV BP365, 12-V, 65-W PV Module (When exposed to sunlight, it will generate 12 V of DC power into a battery, up to a maximum of 65 W according to the strength of solar radiation.)WG Rutland WG913 Windcharger (1 × Main Generator Assembly, 6 × Aerofoil Blades).The 12 V of DC power generated relates directly to the wind speed at the turbine.FC/cell Size (L × W × H: 65 × 85× 38 mm), maximum power of 0.25 W, maximum current of 1 A, maximum voltage of 0.9 V. ESS IRF18650P, Lithium iron phosphate battery (3.2 V, 1100 mAh).

12 Table 1 .
Size (L × W × H: 65 × 85× 38 mm), maximum power of 0.25 W, maximum current of 1 A, maximum voltage of 0.9 V.ESSEnergies 2016, 9, 781 4 of Specifications of the devices.PV BP365, 12-V, 65-W PV Module (When exposed to sunlight, it will generate 12 V of DC power into a battery, up to a maximum of 65 W according to the strength of solar radiation.)WG Rutland WG913 Windcharger (1 × Main Generator Assembly, 6 × Aerofoil Blades).The 12 V of DC power generated relates directly to the wind speed at the turbine.FC/cell Size (L × W × H: 65 × 85× 38 mm), maximum power of 0.25 W, maximum current of 1 A, maximum voltage of 0.9 V. ESS IRF18650P, Lithium iron phosphate battery (3.2 V, 1100 mAh).

Figure 2 .
Figure 2. Normalized membership function of the  load P .

Figure 2 .
Figure 2. Normalized membership function of the ∆P load .

Figure 3 .
Figure 3. Normalized membership function of the  DG P .

Figure 4 .
Figure 4. Normalized membership function of the control signal.
p , constrained to lie within some region D. Let    k F F , where k runs over a set which may be finite, enumerable, or continuous, be a set of transformations with the property that  p D implies that ( ) k F p D for all k.

Figure 3 .
Figure 3. Normalized membership function of the ∆P DG .

Figure 3 .
Figure 3. Normalized membership function of the  DG P .

Figure 4 .
Figure 4. Normalized membership function of the control signal.
p , constrained to lie within some region D. Let    k F F , where k runs over a set which may be finite, enumerable, or continuous, be a set of transformations with the property that  p D implies that ( ) k F p D for all k.

Figure 4 .
Figure 4. Normalized membership function of the control signal.

Figure 5 .
Figure 5.The voltage levels in the microgrid (MG) for case A.

Figure 6 .
Figure 6.The power flows in the MG for case A.

Figure 5 .
Figure 5.The voltage levels in the microgrid (MG) for case A.

Figure 5 .
Figure 5.The voltage levels in the microgrid (MG) for case A.

Figure 6 .
Figure 6.The power flows in the MG for case A.

Figure 6 .
Figure 6.The power flows in the MG for case A.

Figure 7 .
Figure 7.The voltage levels in the MG for case B.

Figure 8 .
Figure 8.The power flows in the MG for case B. .

Figure 7 .
Figure 7.The voltage levels in the MG for case B.

Figure 7 .
Figure 7.The voltage levels in the MG for case B.

Figure 8 .
Figure 8.The power flows in the MG for case B.

Figure 8 .
Figure 8.The power flows in the MG for case B.

Figure 9 .
Figure 9.The voltage levels in the MG for case C.

Figure 10 .
Figure 10.The power flows in the MG for case C.

Figure 9 .
Figure 9.The voltage levels in the MG for case C.

Figure 10 .
Figure 10.The power flows in the MG for case C.

Table 1 .
Specifications of the devices. PV

Table 1 .
Specifications of the devices.

Table 2 .
Fuzzy rules for  load P and  DG P .

Table 2 .
Fuzzy rules for  load P and  DG P .

Table 2 .
Fuzzy rules for ∆P load and ∆P DG .