Sommaire As shown in Fig. Traffic of nication sans fil. La performation des techniques de GRD a un admitted calls is then controlled by other RRM techniques such as impact direct sur la performance de chaque usager et celle du scheduling, power control and transmission rate control.

Packet Arrival ticated due to the unique features of wireless networks such as channel multiple access interference, channel impairments, handoff require- ments, and limited bandwidth. Call Arrival CAC schemes developed for cellular wireless networks including sec- ond and third generation systems are extensively studied in the literature.

Reference [2] includes a comprehensive survey on CAC on wireless cellular networks. Call 2. Satellite systems are also considered for overlapping Channel with terrestrial wireless networks to provide service for high mobile Conditions RRM Controller users. However, two main challenges have to be tackled in satellite sys- Power and Rate Control tems. The first challenge is the large propagation delay that limits the adaptation capability of RRM techniques including CAC schemes.

The second challenge is the spectrum partitioning between terrestrial and Figure 1: Radio Resource Management Model satellite systems. CAC is employed to manage the assignment of arriving and resource availability for burst admission in ad-hoc wireless net- calls new and handoff to one of the two layers depending on the call works.

For instance, services with lower delay tolerance are admitted type. The admission decision is probabilistic where the admission first, then services with higher delay tolerance, which can be queued probability is chosen to minimizes the blocking probability subject to until resources become available. Arriving bursts send their requests to constraints on the dropping probability and average percentage of calls the cluster head that manage the resource availability and prioritization assigned to the satellite coverage given a certain bandwidth partition- scheme.

Results show that the proposed scheme outperforms classical ing plan. The second constraint is used to represent the consideration of non-prioritized burst admission schemes such as first-come-first-serve the large propagation delay in the satellite connection.

When a new call known policies, namely, cellular first CF and satellite first SF. The delay incurred by the probing packets is used to determine the service curve, which quantifies A threshold-based CAC scheme has been proposed in [5]. Call admis- the network loading status. The measured service curve is compared by sion is based on resource availability for constant bit rate CBR , a pre-specified service curve corresponding to the QoS requirements.

The threshold the universal service curve; otherwise the call is rejected. In order to avoid excess delay in the resource Three CAC schemes for ad-hoc wireless local are networks LANs allocation, the CAC and other resource management are processed on have been proposed in [14]. The master device node decides whether board the satellite. The three schemes CAC for variable bit rate VBR real-time and non-real-time and differ mainly in the estimation technique of the aggregate link utiliza- CBR services has been proposed in [6] using a probabilistic measure of tion taking into account the burst nature of the traffic.

The first scheme the QoS guarantee by estimating the excess demand probability which uses the sum of the peak rates of different users as an estimate of the measures the probability of the resource unavailability of all admitted aggregate link utilization. Although this scheme is very simple and can calls in ATM-Satellite network.

Unspecified bit rate UBR is also guarantee a low packet loss rate, the conservative estimate leads to a considered but without any CAC, i. The second and third schemes use description of the signaling and this CAC implementation is provided the effective bandwidth technique to estimate the link utilization. The in [7]. A CAC scheme for voice and data services over low earth orbit probability of the aggregate link utilization is approximated using the satellite LEOS system has been proposed in [8].

The admission deci- Hoeffding bound [15] and Gaussian distribution in the second and third sion is based on the resource availability with a higher priority to the schemes respectively.

Results show that when a low packet loss rate is voice service. A similar strategy is used in [16] for mobile ad-hoc networks compared with classical wireless networks. However, it should be noticed that the utili- networks. These CAC schemes have to consider the lack of infrastruc- zation factor value is sensitive to many systems parameters and it has to ture for ad-hoc networks , network connectivity, new interference be determined for each particular network configuration. This framework tries to strike a balance between the network High altitude aeronautical platform HAAP has been proposed to com- connectivity, which is enhanced by admitting more users and the sig- bine the advantages of terrestrial and satellite systems while avoiding nal quality in terms of the interference level that increases by admitting most of the disadvantages of both systems [17].

As shown in Fig. The CAC concept classifies the incoming user HAAP provides the advantage of covering large areas with minimum as class 1 if by admitting this user the number of links will equal one infrastructure and having centralized system control and global informa- of the critical values otherwise it is classified as class 2.

The critical tion. Nevertheless, the transmitted power in the downlink is significantly numbers of links, determined by the graph theory, are the ones that limited compared with terrestrial wireless networks. The first algorithm restricts the maxi- can reach other nodes using one or more hops, where n is the number mum transmitted power per base station BS while the second scheme of existing nodes.

Class 1 users are admitted if the advantage of restricts the maximum total transmitted power. It is clear that the second increasing the connectivity by admitting those users compensates the algorithm is more efficient and causes less blocking due to the statisti- degradation in the signal quality due to the potential increase in the cal multiplexing.

In both algorithms the call is only admitted if the SIR interference level while class 2 users are only admitted if the interfer- constraints of all users in all cells can be satisfied without violating the ence level after admitting the incoming users is acceptable. The maximum power constraint. The CAC scheme admits the network has been proposed. The SIR is vation at included nodes are employed to explore the possibility of checked by calculating the total received power at all BSs.

If no routes could be found such that all nodes in that route can be allocated the required resources in terms of the number of tome slots, the call is rejected. Time slot reallocation 5. Before admitting a direct the incoming calls to the proper layer. The first one, called Uniform Call Admission pair is checked whether it is less than a threshold value. Also, the sum UCA , directs all calls voice or data, new or handoff to micro-cells of the number in each pair of circuits intersecting at any node is first.

If no channels are available in the micro-cells, the incoming call is checked to ensure that it is less than another threshold value. The then redirected to the corresponding macro-cells. Unlike UCA, the sec- threshold values are chosen to minimize the blocking probability using ond algorithm, called Non-uniform Call Admission NCA , directs data the ordinal optimization techniques.

There- fore, novel RRM in general and CAC in particular are needed to deal with the anticipated new composite radio wireless environment. Perros; K. Cell j [3]. ElBatt; A. Wieselthier; C. Barnhart; A. A CAC scheme has been proposed in [21] to maximize the system [5].

Koraitim; S. Modi- fied linear programming techniques used to solve the optimization [6]. Iera; A. Molinaro, S. This is because the [7]. Molinaro; G. Aloi; S. Manhattan-model assumed in the micro-cells.

Ween; A. Qureshi; M. Kraetz; M. It has been shown that CAC schemes play a [9]. Chiang, M. In future wireless networks, RRM is becoming more challenging because of the anticipated heterogeneous environment taking into It is anticipated that different networks and access technologies will coexist JSAC , vol.

Ayyagari; A. Valaee; B. Razzano; A. Letters to the Editor [15]. Subject: Erratum , pp. Dear Editor, [16]. Dong; D. Makrakis; T. The article should 2, April , pp. Djuknic; J.

This additional reference is: [18]. Foo; W. Lim; R. Lauzon, P. Hubbard, C. Rabbath, E. Gagnon, B. Kim, P. Rabbath, Dr. Gagnon and Mr. Lauzon Sept. Quebec [20]. Pandey; D. Ghosal; B. Dear Editor, [21]. Ho; C. Lea; G. Wang; G. Vehicular Technology, vol. Ying; Z. Jingmei; W. Weidong; Z. Using different techniques and improvements, Can you point me to a good article or paper on the light bulb invention story. I would like to understand why there is an apparent lack of agree- ment on who invented the light bulb.

From Mar. Newfoundland, as an assistant professor. This year, , is the first that the Mobile Computing. He won the Ontario Graduate Scholarship for Prize has been awarded.

He can be reached at Congratulations! Shoaib A. The authors presented the FACTS models in simple and effective way their research work from and also added Matlab code in the book which is an added dimension. November 25, - Published on Amazon. It seems that the authors could have saved pages giving a CD with the codes. Before buying, it is a good idea to check other books. Go to Amazon. Discover the best of shopping and entertainment with Amazon Prime.

Prime members enjoy FREE Delivery on millions of eligible domestic and international items, in addition to exclusive access to movies, TV shows, and more. Google Classroom. Teaching Technology. Educational Technology.

Blended Learning. Science Classroom. Find this Pin and more on ideas by Marianna Weimann. Multimedia Technology. Medical Technology. Energy Technology.

Computer Technology. Technology Gadgets. Hardware Software. Cloud Computing. Computer Architecture. Find this Pin and more on titopomi by Mihaya. Find this Pin and more on Budgeting by Budgeting Tips.

Making A Budget. Create A Budget. Diy On A Budget. Monthly Budget Planner. Savings Planner. Living On A Budget. Family Budget. Frugal Living. Financial Budget. Money Tips. Online Savings Account. Savings Calculator. Money Challenge. Managing Your Money. Selected type: E-Book. Added to Your Shopping Cart. View on Wiley Online Library. Elgerd, O. Grainger, J.

Harker, B. Hesse, M. Laughton, M. Wasley, R. Wedephol, L. Weedy, B. This is a problem of great engineering complexity where the following operational policies must be observed: 1 nodal voltage magnitudes and system frequency must be kept within narrow boundaries; 2 the alternating current AC voltage and current waveforms must remain largely sinusoidal; 3 transmission lines must be operated well below their thermal and stability limits; and 4 even short-term interruptions must be kept to a minimum.

Moreover, because of the very competitive nature of the electricity supply business in an era of deregulation and open access, transmission costs must be kept as low as possible. If voltage magnitudes are outside bounds in one or more points of the network, then appropriate action is taken in order to regulate such voltage magnitudes.

The set of equations representing the power system are nonlinear. In its most basic form, these equations are derived by assuming that a perfect symmetry exists between the phases of the three-phase power system Arrillaga and Arnold, This applies to both active power and reactive power. PGk and QGk represent, respectively, the active and reactive powers injected by the generator at bus k.

PLk and QLk represent the active and reactive powers drawn by the load at bus k, respectively. Provided the nodal voltages throughout the power network are known to a good degree of accuracy then the transmitted powers are easily and accurately calculated. However, if the nodal voltages are not known precisely then the calculated transmitted powers will have only approximated values and the corresponding mismatch powers are not zero.

Upon convergence, the nodal voltage magnitudes and angles yield useful information about the steady-state operating conditions of the power system and are known as state variables. Based on Figure 4. It should be remarked that Equations 4. However, a practical power system will consist of many buses and many transmission elements. This calls for Equations 4.

This is illustrated in Figures 4. In the broadest sense, one can think of voltage magnitudes and phase angles as state variables, and active and reactive powers as control variables. Furthermore, the active and reactive powers drawn by the load PL and QL are known from available measurements. Constant voltage operation is possible only if the generator reactive power design limits are not violated, that is, QG min 4.

Early approaches were based on loop equations and numerical methods using Gauss-type solutions. Nevertheless, reliability towards convergence was still the main concern. Further developments led to the introduction of the Gauss—Seidel method with acceleration factors.

The drawback is that these algorithms exhibit poor convergence characteristics when applied to the solution of networks of realistic size Elgerd, Based on this premise, all high-order derivative terms in Equation 4. Also, the rows and columns corresponding to reactive power and voltage magnitude for PV buses are discarded. Furthermore, when buses k and m are not directly linked by a transmission element, the corresponding k—m entry in the Jacobian is null.

An additional characteristic is that they are symmetric in structure but not in value Zollenkoff, Consider the lth element connected between buses k and m in Figure 4. Also, QG in PV buses are calculated in each iteration in order to check if the generators are within reactive power limits.

Details of this computation are given in the next section. One of the main strengths of the Newton—Raphson method is its reliability towards convergence. Contrary to non-Newton—Raphson solutions, such a characteristic is independent of the size of the network being solved and the number and kinds of control equipment present in the power system.

Aspects that may dent its quadratic convergence performance are reactive power limit violations in generator PV buses and extreme loading conditions. The initial voltage phase angles are selected to be 0 at all buses. The nodal voltage magnitude at bus k is allowed to vary and Vk becomes a state variable. It should be remarked that bus k may revert to being a generator PV bus at some point during the iterative process if better estimates of Qcal k , calculated with more accurate nodal voltages, indicate that the reactive power requirements at bus k can, after all, be met by the generator connected at bus k.

Hence, reactive power limit checking is carried out at each iteration. The switching of buses from PV to PQ and vice versa impose additional numerical demands on the iterative solution and retard convergence.

Bus 1 is selected to be the slack bus and bus 2 is a generator bus. Bus 3 contains no generation and becomes a load bus. A transformer and a transmission line link buses 1 and 2 and buses 2 and 3, respectively.

One shunt element and one load are connected at bus 3. At a given bus, the power balance is obtained by adding the contribution of each plant component connected to that bus. Redrawn, with permission, from C. Acha, S. Tan, and J. The contribution of all three buses is shown in this example for completeness, but it should be remembered that in actual calculations active and reactive power mismatch entries are not required for the slack bus.

Likewise, the reactive power mismatch entry is not required for the generator PV bus. Reproduced, with permission, from C. Nevertheless, the basic procedure illustrated above, based on superposition, will also apply to the formation of the Jacobian. For each plant component, relevant Jacobian equations are chosen based on the type of buses to which the plant component is connected.

These buses determine the location of the individual Jacobian terms in the overall Jacobian structure. The contributions of the line, transformer, and shunt components to the Jacobian are shown in Figure 4.

It should be noted that entries for the slack bus and the reactive power entry of the generator bus are not considered in the Jacobian structure. The program is general, as far as the topology of the network is concerned, and caters for any number of PV and PQ buses.

Moreover, any bus in the network may be designated to be the slack bus. No transformers are represented in this base program and no sparsity techniques Zollenkoff, are incorporated. These assumptions are based on physical properties exhibited by electrical power systems, in particular in highvoltage transmission systems. The power mismatch equations of both methods are identical but their Jacobians are quite different; the Jacobian elements of the Newton—Raphson method are voltage-dependent whereas those of the fast decoupled method are voltage-independent i.

Moreover, the number of Jacobian entries used in the fast decoupled method is only half of those used in the Newton—Raphson method.

The trade-off lies in the weakening of the strong convergence characteristic exhibited by the Newton—Raphson method; the convergence characteristics of the fast decoupled method are linear as opposed to quadratic.

However, an asset of the fast decoupled method is that one of its iterations only takes a fraction of the time required by one of the Newton—Raphson method iterations.

Incorporating these assumptions in the Jacobian elements of Equations 4. Matrices B0 and B00 are identical if no generator buses exist in the system. However, in the more general case, when generator buses do exist in the system then the row and column corresponding to each generator bus are removed from matrix B Equations 4. Matrices B0 and B00 are symmetric in structure and, provided no phase-shifting transformers are present in the system, they are also symmetric in value.

This is in contrast to with the Newton—Raphson method, where the Jacobian is evaluated and inverted factorised; Zollenkoff, at each iteration. As shown in Figure 4. The data are given in function PowerFlowsData, suitable for use with either the Newton—Raphson or the fast decoupled Matlab1 programs.

The former method takes 6 iterations to converge whereas the latter takes 27 iterations. However, it should be mentioned that one iteration of the fast decoupled method executes much faster than one iteration of the Newton—Raphson method; no inversion refactorisation of the Jacobian is required at each iterative step of the fast decoupled. It can be observed from the results presented in Table 4. From G.

Stagg and A. Phase angle deg North 1. This is also the transmission line that incurs higher active power loss i. The active power system loss is 6. The operating conditions demand a large amount of reactive power generation by the generator connected at North i.

This amount is well in excess of the reactive power drawn by the system loads i. The generator at South draws the excess of reactive power in the network i.

This amount includes the net reactive power produced by several of the transmission lines. The model makes provisions for complex taps on both the primary and the secondary windings, and the magnetising branch of the transformer is included to account for core losses.

However, the LTC model does not require complex taps, and Equation 3. This is with a view to developing LTC models aimed at systems applications. Also, the subscript sc is dropped in the transformer admittance terms. Comprehensive bus power injection equations for the LTC transformer may be derived based on Equation 4. Simpler expressions may be derived if a number of practical assumptions are introduced in this equation.

Incorporating these simplifying assumptions in Equation 4. For this mode of operation Vk is maintained constant at the target value. The Jacobian elements in matrix Equation 4. It resembles a generator PV bus but here the voltage control is exerted by an LTC as opposed to a generator. It should be remarked that a more comprehensive set of nodal power equations may be derived for the two-winding transformer by basing the power equation derivations on Equation 4.

There is no need to assume that the transformer impedance is all placed on the primary side. Also, the effect of the magnetising admittance may be included in the nodal power equations of the LTC transformer. The status of LTC taps is checked at each iterative step to assess whether or not the LTC is still operating within limits and capable of regulating voltage magnitude.

The nodal voltage magnitude at bus k is allowed to vary and Vk replaces Tk as the state variable. Similar criteria would apply if the LTC tapping facilities were on the secondary winding, with Um and Tk changing roles in Equation 4. Moreover, relevant power equations and Jacobian elements, equivalent to Equations 4.

The linearised Equation 4. The initial condition of the tap is set to a nominal value i. The winding impedance contains no resistance, and an inductive reactance of 0. The data given in function PowerFlowsData in Section 4. The transmission line originally connected between Lake and Main is now connected between Lakefa bus 6 and Main bus 4. The nodal voltages are given in Table 4. It should be noted that the LTC upholds the target value of 1 p. Table 4. Phase angle deg 1. It is also interesting to note that the LTC achieves its voltage regulation objective at the expense of consuming reactive power; it draws 10 MVAR from the system.

The system active power loss is 6. To show the prowess of the Newton—Raphson method towards convergence, in Table 4. If the generator hits one of its reactive limits then the master LTC tap becomes active and the bus is converted to PVT; the bus becomes controlled by the LTC as opposed to the generator. The control of nodal voltage magnitude by the generator has higher priority.

If the set of LTCs associated with a given generator are controlling buses different from the generator bus and the generator reaches one of its reactive limits then the LTC is switched to control the generator bus so that it changes to a PVT bus.

These control actions are shown schematically in Figure 4. The LTC tap, located on the primary winding, is used to control voltage magnitude when the generator violates its minimum reactive power limit. The winding impedance contains no resistance and an inductive reactance of 0. Once the generator violates reactive limits the LTC becomes active. For the condition when the target voltage magnitude at South is 1 p. Convergence is obtained in 7 iterations. The nodal voltages are very similar to the base case presented in Section 4.

The value of LTC tap required to achieve 1 p. As expected, the LTC achieves its operating point at the expense of consuming reactive power. However, in this case it draws only 2. The system active power loss increases to 6. It is derived from the two-winding, single-phase transformer model presented by Section 3.

Comprehensive bus power injection equations for the phase shifter may be derived with reference to Equation 3. However, simpler expressions may be derived if some practical assumptions are introduced at this stage. For instance, it is reasonable to assume that the phase-changing facility is only on the primary side, i.

Based on Equation 4. For instance, the effect of the magnetising admittance may be included in the nodal power equations of the transformer. Checking of phase-shifter tap limits starts from iteration 1. This bus is added to enable connection of the phase shifter. The initial value of the complex tap is set to the nominal value i.

The winding contains no resistance, and an inductive reactance of 0. The phase shifter upholds its target value. The maximum absolute power mismatches of the system buses and phase shifter are shown in Table 4. As expected, the nodal voltage magnitudes do not change compared with the base case presented in Section 4. This is slightly more than a twofold increase in transmitted power, and the phase angle difference changes from —0.

The phase shifter achieves its operating point at the expense of consuming 1. Figure 4. For instance, simulations are presented in Table 4. The size of the feasible active power control region is a function of the phase angle controller range; as the range increases so too does the sizes of the regions. Note: ptr, pointer. However, this programming philosophy incurs excessive cpu overheads. Pointers are used to move from one structure to another. This is illustrated in Figure 4.

An array of pointers is created, the size of which equals the number of rows in the matrix. Each element points to the address of the start of a list. Moreover, one list is created for each row. Each list consists of an array of structures used to store information associated with off-diagonal Jacobian elements. The information associated with diagonal elements is stored in a separate array of structures. The result may be poor convergence, or more seriously, divergent solutions. Such unwanted problems can be avoided quite effectively by limiting the size of correction, with the actual computed adjustments being replaced by truncated adjustments.

This is a straightforward software solution to a common problem when dealing with utilitysize power systems. The network contains two generators and four synchronous condensers.

Transformers connected between buses 4—12, 6—10, and 27—28 are taken to be LTC transformers. The nodal voltage magnitudes at buses 4, 6, and 12 are controlled at 1, 1, and 1. Mismatch p. However, the algorithm fails to reach convergence if the state variable increments are not truncated. It is assumed in the study that none of the LTCs violates tap limits.

The active and reactive powers generated by the two synchronous generators GE and four synchronous condensers CO are shown in Table 4. Note: GE, generator; CO, condenser.

This situation is shown in Figure 4. The parallel condition occurs when bus k is regulated by two or more LTCs, as shown in Figure 4. It must be noted that buses m and n may not necessarily be electrically connected. When two or more LTCs are controlling one nodal voltage magnitude multiple solutions become a possibility because the number of unknown variables is greater than the number of equations. An entire group of parallel LTCs may be treated as a single control criterion if they are started from the same tapping initial condition.

This Equation is linearised with respect to the common tap and incorporated in the overall Jacobian Equation 4. From the LTC set, the LTC that draws less reactive power is selected to be the master, and its tapping position becomes the master tapping position.

Moreover, a new master is selected from the remaining active LTCs. The tap is adjusted by using Equation 4. The nodal voltage magnitude at bus 6 is kept at 1. The voltage magnitude at buses 4 and 27 are controlled at 1.

The transformers reactance and offnominal tap values given in Freris and Sasson, are taken to be on the secondary and primary windings, respectively. The primary windings of the four transformers are assumed connected to buses 6, 4, and 27, respectively. A comparison is made for the various cases given in Table 4.

However, the use of sensitivity factors guarantees better results in terms of the number of iterations required to get to the solution, compared with the case in which identical tapping updates is carried out. Owing to the idiosyncrasies of the electrical power network, the mathematical model that describes its operation during steady-state is nonlinear.

Furthermore, for most practical situations, the power network is a very large-scale system. The additional burden imposed on the numerical solution by the many constraint actions resulting from the various power system controllers in the network does not impair the ability of the Newton—Raphson method to converge in a quadratic fashion.

Looks like you are currently in Finland but have requested a page in the United States site. Would you modellong to change to the United Facts modelling and simulation in power networks free download site? Enrique ModelingClaudio R. Claudio Free download windows 8.1 iso 64 bit. Undetected location. NO YES. Selected type: E-Book. Added to Your Shopping Cart. View on Wiley Online Library. This is a dummy description. Table of contents Preface. Modelling of Facts modelling and simulation in power networks free download Power Plant. Conventional Power Flow. Three-phase Power Flow. Optimal Power Flow. Power Flow Tracing. Reviews "I certainly recommend this book to all power system planning engineers and students who wish to follow careers in this area. FACTS: Modelling and Simulation in Power Networks [Acha, Enrique, Fuerte-Esquivel, Get your Kindle here, or download a FREE Kindle Reading App. FACTS: Modelling and Simulation in Power Networks - Kindle edition by Acha, Download it once and read it on your Kindle device, PC, phones or tablets. Read with the free Kindle apps (available on iOS, Android, PC & Mac), Kindle. Request PDF | FACTS: Modelling and Simulation in Power Networks | Background Flexible Alternating Current Transmission SystemsInherent Limitations of. FACTS: modelling and simulation in power networks Newsletter Editor - Jeffery Mackinnon Free of charge to all IEEE members in Canada. For IEEE members. This book is printed on acid-free paper responsibly manufactured from FACTS: Modelling and Simulation in Power Networks, is to provide a. FACTS FACTS Modelling and Simulation in Power Networks Enrique Acha University of Glasgow, UK Claudio downloads Views 5MB Size Report. FACTS: Modelling and Simulation in Power Networks Book Free Download. Find this Pin and more on Electrical engineer by James Lee. Tags. Esquivel. May 24, - FACTS: Modelling and Simulation in Power Networks Book Free Download. FACTS: Modelling and Simulation in Power Networks Enrique Acha, Claudio R. Fuerte-Esquivel, Hugo Ambriz-Pérez, César Angeles-Camacho. John Wiley. Defend I. Provides a thorough grounding in the mathematical representation of FACTS controllers and power plant components. Lichen Biology. Whether you've loved the book or not, if you give your honest and detailed thoughts then people will find new books that are right for them. Igneous and Metamorphic Petrology. Maya 5 Fundamentals. English Choose a language for shopping. Report Writing for Architects and Project Managers. Tell the Publisher! Offering a comprehensive understanding of FACTS controllers with the power network, this book is an asset to electrical and electronic engineers involved in the planning, design and operation of power systems.- duty free charles de gaulle terminal 2a, free online photo resizer in mm, free live streaming psg vs rennes, free mobile phone contacts transfer software, crew neck soccer t shirt mockup front view free, free angel card reading new age, forfait free mobile illimit? sans engagement FACTS: Modelling and Simulation in Power Networks | WileyDescription