A survey on modeling of microgrids—from fundamental
physics to phasors and voltage sources
Johannes Schiffera,∗, Daniele Zonettib , Romeo Ortegab , Aleksandar Stankovićc , Tevfik Sezid , Jörg Raische
a School

of Electronic and Electrical Engineering, University of Leeds, Leeds LS2 9JT, UK
des Signaux et Systémes, École Supérieure d0 Electricité (SUPELEC), Gif-sur-Yvette 91192, France
c Tufts University, Medford, MA 02155, USA
d Siemens AG, Smart Grid Division, Energy Automation, Humboldtstr. 59, 90459 Nuremberg, Germany
e Technische Universität Berlin, Einsteinufer 11, 10587 Berlin, Germany and Max-Planck-Institut für Dynamik komplexer technischer
Systeme, Sandtorstr. 1, 39106 Magdeburg, Germany

arXiv:1505.00136v2 [] 17 May 2016

b Laboratoire

Abstract
Microgrids have been identified as key components of modern electrical systems to facilitate the integration of renewable
distributed generation units. Their analysis and controller design requires the development of advanced (typically
model-based) techniques naturally posing an interesting challenge to the control community. Although there are widely
accepted reduced order models to describe the dynamic behavior of microgrids, they are typically presented without
details about the reduction procedure—hampering the understanding of the physical phenomena behind them. Preceded
by an introduction to basic notions and definitions in power systems, the present survey reviews key characteristics and
main components of a microgrid. We introduce the reader to the basic functionality of DC/AC inverters, as well as
to standard operating modes and control schemes of inverter-interfaced power sources in microgrid applications. Based
on this exposition and starting from fundamental physics, we present detailed dynamical models of the main microgrid
components. Furthermore, we clearly state the underlying assumptions which lead to the standard reduced model with
inverters represented by controllable voltage sources, as well as static network and load representations, hence, providing
a complete modular model derivation of a three-phase inverter-based microgrid.
Keywords: Microgrid modeling, microgrid analysis, smart grid applications, inverters
1. Introduction
1.1. Motivation
It is a widely accepted fact that fossil-fueled thermal
power generation highly contributes to greenhouse gas emissions [66, 69, 67]. In addition, a growing stream of scientific results [51, 45, 97] has substantiated claims that these
emissions are a key driver for climate change and global
warming. As a consequence, many countries have agreed
to reduce their greenhouse gas emissions.
Apart from a reduction of energy consumption, e.g.,
through an increase in efficiency, one possibility to reduce
greenhouse gas emissions is to shift the energy production
from fossil-fueled plants towards renewable sources [66, 67,
19]. Therefore, the worldwide use of renewable energies
has increased significantly in recent years [101].
∗ Corresponding author J. Schiffer. Tel. +44 (0)113 343 9719.
Fax +44 (0)113 343 2032.
Email addresses: j.schiffer@leeds.ac.uk (Johannes Schiffer),
zonetti@lss.supelec.fr (Daniele Zonetti),
ortega@lss.supelec.fr (Romeo Ortega), astankov@ece.tufts.edu
(Aleksandar Stanković), tevfik.sezi@arcor.de (Tevfik Sezi),
raisch@control.tu-berlin.de (Jörg Raisch)

Preprint submitted to Elsevier

Unlike fossil-fueled thermal power plants, the majority of renewable power plants are relatively small in terms
of their generation power. An important consequence of
this smaller size is that most of them are connected to the
low voltage (LV) and medium voltage (MV) levels. Such
generation units are commonly denoted as distributed generation (DG) units [1]. In addition, most renewable DG
units are interfaced to the network via DC/AC inverters.
The physical characteristics of such power electronic devices largely differ from the characteristics of synchronous
generators (SGs), which are the standard generating units
in existing power systems. Hence, different control and
operation strategies are needed in networks with a large
amount of renewable DG units [40, 104, 101].
1.2. The microgrid concept
One potential solution to facilitate the integration of
large shares of renewable DG units are microgrids [61, 46,
40, 55, 19, 99]. A microgrid gathers a combination of generation units, loads and energy storage elements at distribution or sub-transmission level into a locally controllable
system, which can be operated either in grid-connected
mode or in islanded mode, i.e., in a completely isolated
May 18, 2016

manner from the main transmission system. The microgrid concept has been identified as a key component in
future electrical networks [19, 32, 62, 99].
Many new control problems arise for this type of networks. Their satisfactory solution requires the development of advanced model-based controller design techniques
that often go beyond the classical linearization-based nestedloop proportional-integral (PI) schemes. This situation
has, naturally, attracted the attention of the control community as it is confronted with some new challenging control problems of great practical interest.
It is clear that to carry out this task it is necessary
to develop a procedure for assembling mathematical models of a microgrid that reliably capture the fundamental
aspects of the problem. Such models have been developed by the power systems and electronics communities
and their pertinence has been widely validated in simulations and applications [20, 53, 81, 74]. However, these
are reduced or simplified, i.e., linearized, models that are
typically presented without any reference to the reduction
procedure—hampering the understanding of the physical
phenomena behind them.
1.3. Existing literature
For the purposes of this survey, previous work on microgrid modeling can be broadly categorized into two classes.
The first class focusses on modeling and control of inverterinterfaced DG units in microgrid applications, but the
model derivation is restricted to individual DG units and
the current and power flows between different units are
not considered explicitly [40, 81, 53, 112, 84, 41, 10, 11].
The second class discusses models of microgrids including electrical network interactions, but the model derivation is based on linearization (i.e., the so-called smallsignal model) [54, 81, 74]. Furthermore, this class of modeling is often tied to specific network control schemes,
such as droop control [20, 81] or to specific test networks
[53, 54, 74, 73]. Building on this previous work and in a
survey-like manner, the present paper brings both aforementioned classes together to formulate a generic modular
model of a microgrid.
Going beyond a mere review of existing microgrid models, we employ model reduction via a time-scale separation
together with the subsequent derivation of the well-known
power flow equations, which is a standard procedure in
SG-based networks [105]. A similar approach has also
been employed for microgrids in [72, 71, 83, 68]. However, neither reference provides a detailed model derivation
for inverter-interfaced units. Also, the analysis in [72, 71]
is restricted to an AC microgrid consisting of two inverters connected via a resistive-inductive line and two local
resistive-inductive loads, while the modeling procedure of
the present paper applies to networks with generic topology and arbitrary number of units.

2

1.4. About the survey
The present survey is an attempt to provide a guideline for control engineers attracted by this fundamental
application for Smart Grids to assess the importance of
the main dynamical components of a three-phase inverterbased microgrid as well as the validity of different models
used in the power literature. To this end, we at first review some fundamental concepts and definitions in power
systems, including a survey on the notion of instantaneous power. Subsequently, we introduce the reader to
the microgrid concept and discuss its main components.
We illustrate that inverter-interfaced units are the main
new elements in future power networks, detail the basic functionality of inverters and review the most common operation modes of inverter-interfaced units together
with their corresponding control schemes. This paves the
path for—starting from fundamental physics—presenting
detailed dynamical models of the individual microgrid components. Subsequently, we clearly state the underlying
assumptions which lead to the standard reduced model
with inverters represented by controllable voltage sources,
as well as static network and load representations. This
reduced model is used in most of the available work on
microgrid control design and analysis [96, 90, 2, 76, 29].
We focus on purely inverter-based networks, since inverter-interfaced units are the main new elements in microgrids compared to traditional power systems. However,
we remark that the employed modeling and model reduction techniques can equivalently be applied to standard
bulk power system models as well as to power systems
with mixed generation pool. For modeling of traditional
electro-mechanical SG-based units, the reader is referred
to standard textbooks on power systems [59, 69, 5].
The main contributions of the present survey paper are
summarized as follows.
• Provide a detailed comprehensive model derivation
of a microgrid based on fundamental physics and
combined with detailed reviews of the microgrid concept, its components and their main operation modes.
• Answer the question, when an inverter can be modeled as a controllable AC voltage source and depict
the necessary underlying model assumptions.
• Show that the usual power flow equations can be
obtained from a network with dynamic line models
via a suitable coordinate transformation (called dqtransformation) together with a singular perturbation argument.
• By combining the two latter contributions, recover
the reduced-order microgrid model currently widely
used in the literature.
We emphasize that the aim of the present survey is not
to give an overarching justification for the final (simplified) model, but to provide a comprehensive overview of

Definition 2.2. A signal x : R≥0 → R3 is said to be a
three-phase AC signal if it is of the form

the modeling procedure for main microgrid components
together with their dynamics, as well as of the main necessary assumptions, which allow the reduction of model
complexity. Which of the presented models (if any) is appropriate for a specific control design and analysis cannot
be established in general, but has to be decided by the
user. Any model used in simulation and analysis necessarily involves certain assumptions. Therefore, it is of great
importance that the user is aware of the pertinence of the
employed model to appropriately assess the implications
of a model-based analysis.
The remainder of this survey is structured as follows.
Basic preliminaries, such as definitions and transformations in power systems are given in Section 2. The microgrid concept is reviewed in Section 3. A detailed dynamical
model of a microgrid is derived in Section 4. In particular,
common operation modes of inverter-interfaced units are
discussed therein. The model reduction yielding models of
inverters as AC voltage sources and phasorial power flow
equations is conducted in Section 5. The survey paper
is wrapped-up with some conclusions and a discussion of
topics of future research in Section 6.

xABC = col(xA , xB , xC ),
where xA : R≥0 → R, xB : R≥0 → R and xC : R≥0 → R
are AC signals.
A special kind of three-phase AC signals are symmetric
AC three-phase signals, defined below.
Definition 2.3. [4, Chapter 2] A three-phase AC signal
xabc : R≥0 → R3 is said to be symmetric if it can be described by




sin(δ(t))
xa (t)

xabc (t) =  xb (t)  = A(t) sin(δ(t) − 2π
3 ) ,
xc (t)
)
sin(δ(t) + 2π
3

where A : R≥0 → R≥0 is called the amplitude and δ :
R≥0 → S is called the phase angle of the signal.
Clearly, from the preceding definition, a symmetric
three-phase AC signal xabc can be described completely
by two signals: its angle δ and its amplitude1 A.

Notation. We define the sets R≥0 := {x ∈ R|x ≥ 0},
R>0 := {x ∈ R|x > 0} and S := [0, 2π). For a set V,
let |V| denote its cardinality and [V]k denote the set of all
subsets of V that contain k elements. For a set of, possibly
unordered, positive natural numbers V = {l, k, . . . , n}, the
short-hand i ∼ V denotes i = l, k, . . . , n. Given a positive
integer n, we use 0n to denote the vector of all zeros, 1n
the vector with all ones and In the n × n identity matrix.
Let x = col(x1 , . . . , xn ) ∈ Cn denote a column vector with
entries xi ∈ C. Whenever clear from the context, we simply
n
write x = col(xi ) ∈ Cp
. The 2-norm of a vector x ∈ Cn
is denoted by kxk2 = |x1 |2 + . . . + |xn |2 . Let diag(ai ) ∈
Cn×n denote a diagonal matrix with entries ai ∈ C. Let
j denote the imaginary unit. The conjugate transpose of
a matrix A ∈ Cn×n is denoted by A∗ . For a function f :
Rn → R, ∇f denotes the transpose of its gradient. The
operator ⊗ denotes the Kronecker product.

Definition 2.4. [4, Chapter 2] A three-phase AC signal
is said to be asymmetric if it is not symmetric.
Definition 2.5. [48, Chapter 3] A three-phase AC electrical system is said to be symmetrically configured if a
symmetrical feeding voltage yields a symmetrical current
and vice versa.
Definition 2.6. [48, Chapter 3] A three-phase AC power
system is said to be operated under symmetric conditions
if it is symmetrically configured and symmetrically fed.
Examples of symmetric and asymmetric three-phase
AC signals2 are given in Fig. 1. Note that the signal in
Fig. 1b satisfies Definition 2.3 and only differs from the
signal in Fig. 1a in that it possesses a time-varying periodic
amplitude.

2. Preliminaries and basic definitions

Remark 2.7. The terms “balanced” and “unbalanced”
are frequently used as synonyms of “symmetric”, respectively “asymmetric” in the literature [4, 37].

2.1. Symmetric AC three-phase signals
Definition 2.1. [27] A signal x : R≥0 → R is said to be
an AC signal if it satisfies the following conditions

Remark 2.8. Three-phase electrical power systems consist of three main conductors in parallel. Each of these
conductors carries an AC current. A three-phase system
can be arranged in ∆- or Y-configuration, see Fig. 2. The

1. it is periodic with period T ∈ R>0 , i.e.,
x(t) = x(t + nT ),

∀n ∈ N,

∀t ≥ 0,

2. its arithmetic mean is zero, i.e.,
Z

t

1 To simplify notation the time argument of all signals is omitted,
whenever clear from the context. The same applies to the definition
of signals, i.e., a signal x : R≥0 → R, is defined equally as x ∈ R,
whenever clear from the context.
2 In this paper only AC systems and signals are considered. Therefore, the qualifier AC is dropped from now on.

t+T

x(τ )dτ = 0

∀t ≥ 0.

3

0

−1

S → R3×3 ,

1

xabc

xabc

1

-π

-2π

π

0

−1

2π

-2π

-π

δ

xabc

xabc

0

-π

π

π

0

va

fdq0 (x(t), %(t)) = Tdq0 (%(t))x(t)

−1
>
Note that the mapping (2.1) is unitary, i.e., Tdq0
= Tdq0
.
From a geometrical point of view, the dq0-transformation
is a concatenation of two rotational transformations, see
[77] for further details. The variables in the transformed
coordinates are often denoted by dq0-variables.
The dq0-transformation offers various advantages when
analyzing and working with power systems and is therefore
widely used in applications [77, 5, 4, 101, 112]. For example, the dq0-transformation permits, through appropriate
choice of %, to map three-phase AC signals to constant signals, i.e., to transform periodic orbits into constant equilibria. This simplifies the control design and analysis in
power systems, which is the main reason why the transformation (2.2) is introduced in the present case. In addition, the transformation (2.2) exploits the fact that, in
a power system operated under symmetric conditions, a
three-phase signal can be represented by two quantities.
To see this, let xabc : R≥0 → R3 be a symmetric threephase signal with amplitude A : R≥0 → R≥0 and phase
angle θ : R≥0 → S, as in Definition 2.3. Applying the
mapping (2.1) with some angle % : R≥0 → S to xabc yields


 
r
sin(θ − %)
xd
3
A cos(θ − %) .
xdq0 = xq  = Tdq0 (%)xabc =
2
0
x0

-2π

-π

0

π

2π

(d) Asymmetric three-phase AC
signal resulting of an asymmetric superposition of a symmetric
signal with signals oscillating at
higher frequencies

∆-configuration

vb
vc
ZY

ZY

ZY

(2.2)

0

δ

Y-configuration

(2.1)

Then, fdq0 : R3 × S → R3 ,

Figure 1: Symmetric and asymmetric AC three-phase signals. The
lines correspond to xa ’—’, xb ’- -’, xc ’· · · ’.
Generation

2
2

is called dq0-transformation.

δ

(c) Asymmetric three-phase AC
signal with phases not shifted
equally by 2π
3

2
2


cos(% + 23 π)
sin(%√+ 32 π)  .

1

−1

2π

cos(% − 23 π)
sin(%(t)
− 32 π)
√

2

2π

(b) Symmetric three-phase AC
signal with time-varying amplitude

1

-2π

0

δ

(a) Symmetric three-phase AC
signal with constant amplitude

−1

r cos(%)
2
sin(%)
Tdq0 (%) :=
√
3
2

0

Z∆

Z∆

Z∆

Neutral point

Figure 2: Standard Y- and ∆-configurations of three-phase AC
power systems based on [48, Chapter 3]. ZY ∈ C denotes the
impedance in Y-configuration, while Z∆ ∈ C denotes the impedance
in ∆-configuration.

Hence, x0 = 0 for all t ≥ 0. Therefore and as in this work
only symmetric three-phase signals are considered, it is
convenient to introduce the mapping Tdq : S → R2×3 ,
r 

2 cos(%) cos(% − 23 π) cos(% + 23 π)
Tdq (%) :=
,
3 sin(%) sin(% − 23 π) sin(% + 23 π)
(2.3)
which, when applied to the symmetric three-phase signal
xabc defined above, yields
r


 
3
sin(θ − %)
xd
A
.
xdq =
= Tdq (%)xabc =
cos(θ − %)
xq
2

latter is also called wye-configuration. Frequently, in a system with Y-configuration an additional fourth (grounded)
neutral conductor is used to reduce transient overvoltages
and to carry asymmetric currents [37, Chapter 2]. Such
systems are typically called three-phase four-wire systems.
Most three-phase power systems are four-wire Y-connected
systems with grounded neutral conductor [37, Chapter 2].
However, it can be shown that, under symmetric operating conditions, this fourth wire does not carry any current
and can therefore be neglected [37, Chapter 2].
2.2. The dq0-transformation

In the following, xdq are referred to as the dq-coordinates
of xabc . Note that xabc = Tdq (%)> xdq .

An important coordinate transformation known as dq0transformation in the literature [78, 77, 5, 69, 101, 6, 112]
is introduced.

Remark 2.10. There are several variants of the mapping
(2.1) available in the literature. They may differ from the
mapping (2.1) in the order of the rows and the sign of the
entries in the second row of the matrix given in (2.1), see

Definition 2.9. [5, Chapter 4], [69, Chapter 11] Let x :
R≥0 → R3 and % : R≥0 → S. Consider the mapping Tdq0 :
4

[5, 6, 112]. However, all representations are equivalent in
the sense that they can all be represented by Tdq0 as given
in (2.1) by choosing an appropriate angle % and, possibly,
rearranging the row order of the matrix Tdq0 . The same
applies to the mapping Tdq given in (2.3). Since—with a
slightly different scaling factor—this transformation was
first introduced by Robert H. Park in 1929 [78] it is also
often called Park transformation [101, Appendix A].

Based on the preceding discussion, the following definitions of instantaneous active, reactive and apparent power
under symmetric, but not necessarily steady-state, conditions are used in this work. The definitions are based on [3,
4], in which they are given in αβ-coordinates. For the purpose of the present paper, it is more convenient to equivalently define the instantaneous powers in dq-coordinates.
Definition 2.11. Let vdq (t) and idq (t) be given by (2.5).
The instantaneous three-phase active power is defined as

2.3. Instantaneous power
Power is one of the most important quantities in control, monitoring and operation of electrical networks. The
first theoretical contributions to the definition of the power
flows in an AC network date back to the early 20th century. However, these first definitions are restricted to sinusoidal steady-state conditions and based on the root-meansquare values of currents and voltages. As a consequence,
these definitions of electric power are not well-suited for
the purposes of network control under time-varying operating conditions [4].
The extension of the definition of electrical power to
time-varying operating conditions is called “instantaneous
power theory” in the power system and power electronics
community [4, 101]. The development of this theory already begun in the 1930s with the study of active and nonactive components of currents and voltages [35]. Among
others, relevant contributions are [15, 22, 3, 109, 23, 79,
57].
Today, it is widely agreed by researchers and practitioners [109, 79, 101] that the definitions of instantaneous power proposed in [3] and contained in [4] are wellsuited for describing the power flows in three-phase threewire systems and symmetric three-phase four-wire systems. However, a proper definition of instantaneous power
in asymmetric three-phase four-wire systems with nonzero
neutral current and voltage is still an open (and controversial) field of research [24, 4, 7, 101]. A good overview
of the research history on instantaneous power theory is
given in [101, Appendix B].
Consider a symmetric three-phase voltage, respectively
current, given by




sin(θ)
sin(ϕ)
√
√


iabc = 2I sin(ϕ − 2π
vabc = 2V sin(θ − 2π
3 ) ,
3 ) ,
2π
2π
sin(θ + 3 )
sin(ϕ + 3 )
(2.4)
where
θ
:
R
→
S
and
ϕ
:
R
→
S
are
the
phase
angles
≥0
≥0
√
√
and 2V : R≥0 → R≥0 , respectively 2I : R≥0 → R≥0 ,
the amplitudes of the respective three-phase signal. As
shown in Section 2.2, applying the transformation (2.3) to
the signals given in (2.4) yields
 


√
Vd
sin(θ − %)
vdq =
= 3V
,
cos(θ − %)
Vq
 


(2.5)
√
Id
sin(ϕ − %)
idq =
= 3I
.
Iq
cos(ϕ − %)

>
P (t) := vdq
(t)idq (t) = Vd (t)Id (t) + Vq (t)Iq (t).

The instantaneous three-phase reactive power is defined as


0 1
>
Q(t) := vdq (t)
i (t) = Vd (t)Iq (t) − Vq (t)Id (t).
−1 0 dq
Finally, the instantaneous three-phase (complex) apparent
power is defined as
S(t) := P (t) + jQ(t).

From the above definition, straightforward calculations
together with standard trigonometric identities yield
P (t) = 3V (t)I(t) cos(θ(t) − ϕ(t)),

Q(t) = 3V (t)I(t) sin(θ(t) − ϕ(t)).

It follows that whenever vabc and iabc given in (2.4) possess constant amplitudes, as well as the same frequency,
i.e., θ̇ = ϕ̇, all quantities P, Q and S are constant. Moreover, then the given definitions of power are in accordance
with the conventional definitions of power in a symmetric
steady-state [39, 4, 37]. For further information on definitions and physical interpretations of instantaneous power,
also under asymmetric conditions, the reader is referred to
[109, 79, 30, 110, 4, 101].
Since this work is mainly concerned with dynamics of
generation units, all powers are expressed in “Generator
Convention” [37, Chapter 2]. That is, delivered active
power is positive, while absorbed active power is negative.
Furthermore, capacitive reactive power is counted positively and inductive reactive power is counted negatively.
2.4. Algebraic graph theory
An undirected graph of order n is a tuple G := (V, E),
where V := {n1 , . . . , nn } is the set of nodes and E ⊆ [V]2 ,
E := {e1 , . . . , em }, is the set of undirected edges, i.e., the
elements of E are subsets of V that contain two elements.
In the case of multi-agent systems, each node in the graph
typically represents an individual agent. For the purpose
of the present work, an agent represents a DG or storage
unit, respectively a load. The l-th edge connecting nodes i
and k is denoted by el = {i, k}. By associating an arbitrary
ordering to the edges, the node-edge incidence matrix B ∈
R|V|×|E| is defined element wise as bil = 1, if node i is
5

Main electrical network

the source of the l-th edge el , bil = −1, if i is the sink
of el and bil = 0 otherwise. For further information on
graph theory, the reader is referred to, e.g., [26, 38] and
references therein.

PCC
Transformer
PV

PV
1

Load

3. The microgrid concept

2

4

Load

PV

Microgrids have attracted a wide interest in different
research and application communities over the last decade
[95, 46, 81, 99]. However, the term “microgrid” is not
uniformly defined in the literature [61, 46, 40, 55, 19, 37,
99]. Based on [40, 46, 99], the following definition of an
AC microgrid is employed in this survey paper.

8

Load

3

PV

Load

11
Load
7
PV
Storage

PV

FC
5

Storage

1. It is a connected subset of the LV or MV distribution
system of an AC electrical power system.
2. It possesses a single point of connection to the remaining electrical power system. This point of connection is called point of common coupling (PCC).
3. It gathers a combination of generation units, loads
and energy storage elements.
4. It possesses enough generation and storage capacity
to supply most of its loads autonomously during at
least some period of time.
5. It can be operated either connected to the remaining electrical network or as an independent island
network. The first operation mode is called gridconnected mode and the second operation mode is
called islanded, stand-alone or autonomous mode.
6. In grid-connected mode, it behaves as a single controllable generator or load from the viewpoint of the
remaining electrical system.
7. In islanded mode, frequency, voltage and power can
be actively controlled within the microgrid.

FC

Wind

PV

Load

FC
Load

Definition 3.1. An AC electrical network is said to be an
AC microgrid if it satisfies the following conditions.

9

10

Load
SG
6

Load

PV
Load

Figure 3: Schematic representation of a microgrid. The microgrid
is composed of several DG units, loads and storage devices. The
DG units are inverter-interfaced photovoltaic (PV), fuel cell (FC)
and wind power plants. In addition, a power generation unit is connected to the network via a synchronous generator (SG). The point
of connection of the microgrid to the main network is called point of
common coupling (PCC).

elements is also an important assumption for the inverter
models derived in this paper. An illustration of an exemplary microgrid is given in Fig. 3.
Most of the named DG and storage units are either
DC sources (PV, FC, batteries) or are often operated at
variable or high-speed frequency (wind turbines, microturbines, flywheels). Therefore, they have to be connected to
an AC network via AC or DC/AC inverters [40, 101]. For
ease of notation, such devices are simply called “inverters” in the following. Overviews on existing test-sites and
experimental microgrids around the globe are provided in
the survey papers [46, 8, 63, 80, 43].

According to Definition 3.1, the main components in a
microgrid are DG units, loads and energy storage elements.
Typical DG units in microgrids are renewable DG units,
such as photovoltaic (PV) units, wind turbines, fuel cells
(FCs), as well as microturbines or reciprocating engines
in combination with SGs. The latter two can either be
powered with biofuels or fossil fuels [63, 37].
Typical loads in a microgrid are residential, commercial
and industrial loads [61, 55, 63]. It is also foreseen to
categorize the loads in a microgrid with respect to their
priorities, e.g., critical and non-critical loads. This enables
load shedding as a possible operation option in islanded
mode [61, 63].
Finally, storage elements play a key-role in microgrid
operation [63, 37]. They are especially useful in balancing
the power fluctuations of intermittent renewable sources
and, hence, to contribute to network control. Possible
storage elements are, e.g., batteries, flywheels or supercapacitors. The combination of renewable DGs and storage

Remark 3.2. While not comprised in Definition 3.1, true
island power systems are sometimes also called microgrids
in the literature [46]. This can be justified by the fact
that island power systems operating with a large share of
renewable energy sources face similar technical challenges
as microgrids. Nevertheless, an island power system differs
from a microgrid in that it cannot be frequently connected
to and disconnected from a larger electrical network [40].
Remark 3.3. Microgrids can also be implemented as DC
systems [85, 86, 60]. Definition 3.1 can easily be adapted to
this scenario by removing the property “frequency control”
in point 7. Recent reviews of the main differences and
challenges for AC and DC microgrids are given in [107,
80, 52].

6

vAC,1

vDC

4. Modeling of inverter-based microgrids
As outlined in the previous section, an AC microgrid
is a spatially small power system the main components of
which are renewable DG units, loads and energy storage elements interconnected through a network of AC transmission lines and transformers. Furthermore, most renewable
DG and storage units are interfaced to the network via
inverters. As a consequence, fundamental network control actions, such as frequency or voltage control, have
to be performed by inverter-interfaced units. This fact
represents a fundamental difference to the operation of
conventional power systems, where mainly SG units are
responsible for network control. Therefore and since the
modeling of SGs is a well-covered topic in the literature
[59, 5, 69], we focus in the following on microgrids with
purely inverter-interfaced DG and storage units. Also, it
is straightforward to incorporate SG-based units into the
microgrid model presented hereafter.
In line with these considerations, an inverter-based microgrid can be represented by an undirected graph G =
(N , E), in which—similarly to [34]—nodes represent voltage buses, edges represent dynamic power lines and the
topology of the network is fully described by the incidence
matrix B. The set of neighbors of node i is denoted by
Ni and contains all k for which el = {i, k} ∈ E. Please
see Section 2.4 for a brief introduction to algebraic graph
theory. In the present case, we further assume that the
set of nodes N can be partitioned into two subsets NI and
NR , associated to inverter and load nodes respectively. We
next proceed as follows. First, we provide a description of
the basic functionality and common operation modes of
inverters that are instrumental for the modeling. Thereafter, we present models of inverters—depending on their
mode of operation—loads, power lines and transformers.
The section is concluded by combining the individual models to an overall representation of an inverter-based microgrid. To enhance readability, the subindex i, preceded by a
comma when necessary, denotes in the sequel the elements
corresponding to the i-th subsystem.

vAC,2

t

t

t
L

vDC

vAC,1

C1

vAC,2

C2

Power
electronics

∼
=
Inverter

Figure 4: Schematic representation of a DC-AC voltage conversion
by a DC/AC inverter. The DC signal vDC : R≥0 → R on the left side
is converted into an AC signal via power semiconductor devices. The
generated AC signal vAC,1 : R≥0 → R3 at the output of the power
electronics is not sinusoidal. Therefore, an LC filter is connected in
series with the power electronics to obtain an approximately sinusoidal ouput voltage vAC,2 : R≥0 → R3 with low harmonic content.

vDC
2

Rf 1

Lf

vDC
vDC
2

va

Rg

Lg v
Ga

vb

vG b

vc

vGc

Cf
Rf 2

Figure 5: Typical circuit of a two-level three-phase inverter with
LC output filter to convert a DC into a three-phase AC voltage.
The inverter is constructed with insulated-gate bipolar transistors
(IGBTs) and antiparallel diodes. The DC voltage is denoted by vDC :
R≥0 → R, the three-phase AC voltage generated by the inverter
with vabc : R≥0 → R3 , vabc = col(va , vb , vc ) and the three-phase
grid-side AC voltage by vG : R≥0 → R3 , vG = col(vGa vGb , vGc ).
The components of the output filter are an inductance Lf ∈ R>0 , a
capacitance Cf ∈ R>0 and two resistances Rf1 ∈ R>0 , respectively
Rf2 ∈ R>0 . Typically, the resistance Rg ∈ R>0 and the inductance
Lg ∈ R>0 represent a transformer or an output impedance. At
the open connectors denoted by “o“ the circuit can be grounded if
desired.

generated AC signal is typically processed through a lowpass filter constructed with LC(L) elements. Further information on the hardware design of inverters and related
controls is given in [31, 75, 112].
In microgrids, two main operation modes for inverters can be distinguished [107, 84]: grid-forming and gridfeeding mode. The latter is sometimes also called gridfollowing mode [55] or PQ control [65], whereas the first
is also referred to as voltage source inverter (VSI) control [65]. The main characteristics of these two different
operation modes are as follows [65, 55, 107, 84].

4.1. Basic functionality and common operation modes of
inverters
Recall that inverters are key components of microgrids.
Therefore, this section is dedicated to the model derivation
of an inverter in a microgrid. The basic functionality of
an inverter is illustrated in Fig. 4. The main elements of
inverters are power semiconductor devices [31, 75]. An exemplary basic hardware topology of the electric circuit of a
two-level three-phase inverter constructed with insulatedgate bipolar transistors (IGBTs) and antiparallel diodes
is shown in Fig. 5. The conversion process from DC to
AC is usually achieved by adjusting the on- and off-times
of the transistors. These on- and off-time sequences are
typically determined via a modulation technique, such as
pulse-width-modulation [31, 75]. To improve the quality
of the AC waveform, e.g., to reduce the harmonics, the

1. Grid-forming mode (also: VSI control).
The inverter is controlled in such way that its output voltage can be specified by the designer. This
is typically achieved via a cascaded control scheme
consisting of an inner current control and an outer
7

voltage control as shown in Fig. 6a, based on [84].
The feedback signal of the current control loop is
the current through the filter inductance, while the
feedback signal of the voltage control loop is the inverter output voltage vabc : R≥0 → R3 . The inner
loop of the control cascade is not necessary to control the output voltage of the inverter and can hence
also be omitted. Nevertheless, it is often included to
improve the control performance.
2. Grid-feeding mode (also: grid-following mode, PQ
control).
The inverter is operated as power source, i.e., it provides a pre-specified amount of active and reactive
power to the grid. The active and reactive power setpoints are typically provided by a higher-level control or energy management system, see [84, 13, 44].
Also in this case, a cascaded control scheme is usually implemented to achieve the desired closed-loop
behavior of the inverter, as illustrated in Fig. 6b. As
in the case of a grid-forming inverter, the inner control loop is a current control the feedback signal of
which is the current through the filter inductance.
However, the outer control loop is not a voltage, but
rather a power (or, sometimes, a current) control.
The feedback signals of the power control are the
active and reactive power provided by the inverter.

However, to the best of the authors’ knowledge, these approaches are not as commonly used as the control schemes
shown in Fig. 6a and Fig. 6b.
4.2. Modeling of grid-forming inverters
A suitable model of a grid-forming inverter for the purpose of control design and stability analysis of microgrids
is derived. There are many control schemes available to
operate an inverter in grid-forming mode, such as PI control in dq-coordinates [81], proportional resonant control
[36, 100] or repetitive control [108, 50] among others. An
overview of the most common control schemes with an
emphasis on H∞ repetitive control is given in [112]. For
a comparison of different control schemes, the reader is
referred to [64]. The assumption below is key for the subsequent model derivation.
Assumption 4.3. Whenever an inverter operated in gridforming mode connects a fluctuating renewable generation
source, it is equipped with a fast-reacting storage.
Assumption 4.3 implies that the inverter can increase
and decrease its power output within a certain range. This
is necessary if the inverter should be capable of providing
a fully controllable voltage also when interfacing a fluctuating renewable DG unit to the network. Furthermore,
since the storage element is assumed to be fast-reacting,
the DC-side dynamics can be neglected in the model. The
capacity of the required DC storage element depends on
the specific source at hand. Generally, the standard capacitive elements of an inverter don’t provide sufficient
energy storage capacity and an additional storage component, e.g., a battery or flywheel, is required if an inverter
is operated in grid-forming mode [101] See [16] for a survey of energy storage technologies in the context of power
electronic systems and renewable energy sources.
Due to the large variety of available control schemes,
it is difficult to determine a standard closed-loop model of
an inverter operated in grid-forming mode together with
its inner control and output filter. Therefore, the approach taken in this work is to represent such a system
as a generic dynamical system. Note that the operation
of the IGBTs of an inverter occurs typically at very high
switching frequencies (2-20 kHz) compared to the network
frequency (45-65 Hz). It is therefore common practice
[65, 40, 81, 74, 18] to model an inverter in network studies
with continuous dynamics by using the so-called averaged
switch modeling technique [31, 18], i.e., by averaging the
internal inverter voltage and current over a suitably chosen
time interval such as one switching period.
It is convenient to partition the set NI into two subsets,
i.e., NI = N1 ∪ N2 , such that N1 contains all nodes associated to grid-forming inverters and N2 contains those associated to grid-feeding inverters. Consider a grid-forming
inverter located at the i-th node of a given microgrid, i.e.,
i ∈ N1 . Denote its three-phase symmetric output voltage

In both aforementioned operation modes, the current
and voltage control loops are, in general, designed with
the objectives of rejecting high frequency disturbances, enhancing the damping of the output LC(L) filter and providing harmonic compensation [82, 12, 74, 81]. Furthermore, nowadays, most inverter-based DG units, such as
PV or wind plants, are operated in grid-feeding mode [84].
However, grid-forming units are essential components in
AC power systems, since they are responsible for frequency
and voltage regulation in the network. Therefore, in microgrids with a large share of renewable inverter-based DG
units, grid-forming capabilities often also have to be provided by inverter-interfaced sources [65, 55].
Remark 4.1. Some authors [107, 84] also introduce a third
operation mode for inverters called grid-supporting mode.
Nevertheless, this last category is not necessary to classify
typical operation modes of inverters in microgrids in the
context of this work, since grid-supporting inverters are
grid-forming inverters equipped with an additional outer
control-loop to determine the reference output voltage.
Such outer control-loops are, e.g., the usual droop controls
[17, 41]. Therefore, the term “grid-supporting inverter” is
not used in the following.
Remark 4.2. In addition to the two control schemes introduced above, there also exist other approaches to operate inverters in microgrid applications. For example,
[9, 106, 111] propose to design the inverter control based
on the model of an SG with the aim of making the inverter mimic as closely as possible the behavior of an SG.
8

vDC
2

Rf 1

Lf

vDC

Rg

vDC
2

Lg v
Ga

vb

vG b

vc

vGc
vDC
2

Cf
Rf 2
Modulator

Current controller

if,abc

iref
Voltage controller

Rf 1

Lf

vDC

Digital signal processor
(DSP)

Digital signal processor
(DSP)

vDC
2

va

vabc

(a) Schematic representation of an inverter operated in grid-forming
mode based on [84]. Bold lines represent electrical connections, while
dashed lines represent signal connections. The current through the
filter inductance is denoted by if,abc : R≥0 → R3 and the inverter
output voltage by vabc : R≥0 → R3 . Both quantities are fed back
to a cascaded control consisting of an outer voltage and an inner
current control. The reference signal vref : R≥0 → R3 for the voltage
controller is set by the designer, respectively a higher-level control.
The IGBTs of the inverter are then controlled via signals generated
by a modulator. The control structure can also be reduced to a pure
voltage control.

Rg

Lg v
Ga

vb

vG b

vc

vGc

Cf
Rf 2
Modulator

Current controller

if,abc

iref
vabc

Power controller
Pref

vref

va

Qref

(b) Schematic representation of an inverter operated in grid-feeding
mode based on [84]. Bold lines represent electrical connections, while
dashed lines represent signal connections. As in Fig. 6a, the current
through the filter inductance is denoted by if,abc : R≥0 → R3 and the
inverter output voltage by vabc : R≥0 → R3 . In grid-feeding mode,
both quantities are fed back to a cascaded control consisting of an
outer power and an inner current controller. The reference active
and reactive powers Pref ∈ R, respectively Qref ∈ R, are set by the
designer or a higher-level control.

Figure 6: Schematic representation of the two main operation modes of inverters in microgrid applications: grid-forming and grid-feeding
mode.

3
by vabc,i : R
q≥0 → R with phase angle αi : R≥0 → S and
amplitude 23 Vi : R≥0 → R≥0 , i.e.,

vabc,i =

r

control scheme discussed, e.g., in [42, 41] this next higher
control level corresponds to the primary control layer of a
microgrid. Let zi : R≥0 → Rp denote the state signal of
this higher-level control system, ui : R≥0 → Rq its input
signal and vref,i its output signal. Furthermore, let gi :
Rp × Rq → Rp and wi : Rp × Rq → R3 be continuously
differentiable functions. Then, the outer control system of
the inverter can be described by




sin(αi )
2 
Vi sin(αi − 23 π) .
3
sin(αi + 23 π)

Furthermore, denote by ωi := α̇i the frequency of the voltage vabc,i . Denote the state signal of the inverter with its
inner control and output filter by xi : R≥0 → Rm , its input signal by vref,i : R≥0 → R3 and its interconnection
port signals by vabc,i and iabc,i : R≥0 → R3 , see Fig. 6a.
Let fi : Rm × R3 × R3 → Rm and hi : Rm × R3 → R3 denote continuously differentiable functions and νi denote a
nonnegative real constant. Then, the closed-loop inverter
dynamics with inner control and output filter can be represented in a generic manner as
νi γi ẋi = fi (xi , vref,i , iabc,i ),
vabc,i = hi (xi , vref,i ),

γi żi = gi (zi , ui ),
vref,i = wi (zi , ui ).

(4.2)

Combining (4.1) and (4.2) yields the overall inverter dynamics for the i-th node, i ∈ N1 ,
γi żi = gi (zi , ui ),
νi γi ẋi = fi (xi , wi (zi , ui ), iabc,i ),

(4.3)

vabc,i = hi (xi , wi (zi , ui )).

(4.1)

4.3. Modeling of grid-feeding inverters and loads
As discussed in Section 4.1, grid-feeding inverters are
typically operated as current or power sources. In order
to achieve such behavior, the control methods employed to
design the inner control loops of grid-forming inverters (see
Section 4.2) can equivalently be applied to operate inverters in grid-feeding mode. The current or power reference
values are typically provided by a higher-level control, e.g.,
a maximum power point tracker (MPPT) [84].

where the positive real constant γi denotes the time-drift
due to the clock drift of the processor used to operate
the inverter, see [91] for further details. Note that iabc,i
represents a disturbance for the inner control system.
One key objective in microgrid applications is to design
suitable higher-level controls to provide a reference voltage
vref,i for the system (4.1) [41]. Within the hierarchical
9

We define the set N` := N2 ∪ NR that contains the
nodes associated to grid-feeding inverters and loads. As
done for the model of a grid-forming inverter in (4.3), let
x`,i : R≥0 → Rr denote the state signal, vabc,i : R≥0 → R3
and iabc,i : R≥0 → R3 denote the interconnection port
signals, f`,i : Rr × R3 → Rr and h`,i : Rr × R3 → R3
denote continuously differentiable functions and κi denote
a nonnegative real constant. We assume then a generic
dynamic model of the form
κi ẋ`,i = f`,i (x`,i , iabc,i ),
vabc,i = h`,i (x`,i , iabc,i ),

For a compact derivation of the network dynamics, it is
convenient to define the aggregated nodal voltages and
currents
vabc := col(vabc,i ) ∈ R3|N | , iabc := col(iabc,i ) ∈ R3|N | ,
the aggregated line voltages and currents
vL := col(vL,l ) ∈ R3|E| ,
as well as the matrices
L := diag(Ll ) ∈ R|E|×|E| ,

(4.4)

Hence, by combining (4.5) with (4.6) the dynamical system
representing the network is given by

4.4. Modeling of power lines and transformers

L ⊗ I3 ẋL = −R ⊗ I3 xL + B > ⊗ I3 vabc ,

The main purpose of the present paper is to provide a
structured modeling procedure for microgrids. For ease of
presentation, we make the following assumption.

iabc = B ⊗ I3 xL .

(4.7)

We next transform the model (4.7) into dq-coordinates
by means of the transformation Tdq introduced in (2.3).
This coordinate transformation is instrumental for the model
reduction carried out in Section 5. Let

Assumption 4.4. All power lines and transformers can
be represented by symmetric three-phase RL elements.
In light of Assumption 4.4 and to ease presentation,
we solely use the term power lines to refer to the network
interconnections in the following sections. Also, note that
it is straightforward to extent the modeling approach presented hereafter to more detailed power line or transformer
models, as well as to DG units interfaced to the network
via SGs.
Recall that the topology of a microgrid can be conveniently described by an undirected graph G = (N , E),
where E denotes the set of power lines interconnecting
the different network nodes i ∼ N . We associate an arbitrary ordering to the power lines el ∼ E. Likewise, we
assign to each power line el ∈ E a three-phase line current
iL,l : R≥0 → R3 . With Assumption 4.4, the power line
el ∈ E connecting a pair of nodes {i, k} ∈ [N ]2 is symmetric, i.e., each phase of the power line el is composed
of a constant ohmic resistance Rl ∈ R>0 in series with a
constant inductance Ll ∈ R>0 and Rl as well as Ll have
the same value for each phase. Furthermore, the voltage
drop across the line is given by

φ := mod2π (ω com t) ∈ S,

(4.8)

where the operator3 mod2π is added to respect the topology of the torus. Applying the transformation Tdq with
transformation angle φ to the signals vabc,i and iabc,i , i ∼
N , gives


 
V̂
Iˆ
v̂dq,i := Tdq (φ)vabc,i = d,i , îdq,i := Tdq (φ)iabc,i = ˆd,i ,
V̂q,i
Iq,i
where the superscript ”ˆ· ” is introduced to denote signals
in dq-coordinates with respect to the angle φ. This notation is used in the subsequent section, where a reducedorder model of a microgrid is derived by using several dqtransformation angles. Furthermore, following standard
notation in power systems, the constant φ̇ = ω com is referred to as the rotational speed of the common reference
frame. Likewise, the signal xL,l in (4.5) becomes


X̂L,d,l
x̂L,dq,l := Tdq (φ)xL,abc,l =
.
X̂L,q,l
Note that

vL,l := vabc,i − vabc,k .

ˆL,dq,l = Ṫdq (φ)xL,l + Tdq (φ)ẋL,l
ẋ


com −X̂L,q,l
=ω
+ Tdq (φ)ẋL,l .
X̂L,d,l

We denote the state of the l-th line by xL,l := iL,l and its
interconnection port variables by vL,l and iL,l . Then, the
model of the l-th line el ∈ E is given by
iL,l = xL,l .

R := diag(Rl ) ∈ R|E|×|E| .

Then, the three-phase interconnection laws can be obtained by following the approach used in [34], where Kirchhoff’s current and voltage laws (KCL and KVL) are expressed in relation to the node-edge incidence matrix B,
i.e.,
iabc = B ⊗ I3 iL ,
B > ⊗ I3 vabc = vL .
(4.6)

for any node i ∈ N` . In addition to grid-feeding inverters, the model (4.4) can equivalently represent impedance
(e.g., R parallel to L), current- or power-controlled loads.
Furthermore, a large variety of other load behaviors can
be modeled by (4.4). We refer the reader to [59, 103] for
further details on load modeling.

Ll ẋL,l = −Rl xL,l + vL,l ,

iL := col(iL,l ) ∈ R3|E| ,

3 The operator mod
2π : R → [0, 2π), is defined as follows: y =
mod2π {x} yields y = x−k2π for some integer k, such that y ∈ [0, 2π).

(4.5)
10

Network

Hence, (4.5) reads in dq-coordinates as




com −X̂L,q,l
ˆ
Ll ẋL,dq,l = Ll ω
+ Tdq (φ)ẋL,l
X̂L,d,l


(4.9)
−X̂L,q,l
= −Rl x̂L,dq,l + Ll ω com
+ v̂dq,l ,
X̂L,d,l

L ⊗ I2 x̂˙ L,dq = (−R ⊗ I2 + X )x̂L,dq + B > ⊗ I2 v̂dq ,
v̂dq = I|N | ⊗ Tdq (φ)vabc ,
iabc = B ⊗ Tdq (φ)> x̂L,dq

vabc,i

îL,dq,l = x̂L,dq,l .

Grid-forming inverters

the aggregated line voltage and currents in dq-coordinates

vabc,k

k ∼ N`

γi żi = gi (zi , ui ),
νi γi ẋi = fi (xi , wi (zi , ui ), iabc,i ),
vabc,i = hi (xi , wi (zi , ui ))

v̂dq := col(v̂dq,i ) ∈ R2|N | , îdq := col(îdq,i ) ∈ R2|N | ,

κk ẋ`,k = f`,k (x`,k , iabc,k ),
vabc,k = h`,k (x`,k , iabc,k )
Grid-feeding inverters and loads

Figure 7: Port-based representation of an inverter based-microgrid
modeled by (4.11), (4.12).

x̂L,dq := col(îL,dq,l ) ∈ R2|E| ,

as well as the matrix


com 0
X := diag Ll ω
1

5. To phasors and voltage sources via time-scale
separation


−1
∈ R2|E|×2|E| ,
0

For the purpose of deriving an interconnected network
model suitable for network control design and stability
analysis, it is customary to make the following assumptions on (4.11), (4.12), where  stands for a generic small
positive real constant.

(4.7) becomes in dq-coordinates
L ⊗ I2 x̂˙ L,dq = (−R ⊗ I2 + X )x̂L,dq + B > ⊗ I2 v̂dq ,
îdq = B ⊗ I2 x̂L,dq .

Assumption 5.1. νi <  in (4.3), i ∼ N1 . Therefore,
ẋi (t) = 0m for all t ≥ 0. Furthermore, vabc,i = wi (zi , ui ),
i ∼ N1 .

(4.10)

4.5. Overall model

Assumption 5.2. κk <  in (4.4), k ∼ N` . Therefore,
ẋ`,k (t) = 0r for all t ≥ 0. Furthermore, the power balance
at each node k ∈ N` can be described by a ZIP model [59],
i.e.,

Pk (v̂dq,k , îdq,k ) = − aP,k kv̂dq,k k22 + bP,k kv̂dq,k k2 + cP,k

By defining the state vectors z ∈ Rp|N1 | , x ∈ Rm|N1 | ,
x` ∈ Rr|N` | , the input u ∈ Rq|N1 | , the matrices
Γ : = diag(γi ) ∈ R|N1 |×|N1 | , V := diag(νi ) ∈ R|N1 |×|N1 | ,

K : = diag(κi ) ∈ R|N2 |×|N2 | ,

:= Pk∗ (kv̂dq,k k2 ),

Qk (v̂dq,k , îdq,k ) = − aQ,k kv̂dq,k k22 + bQ,k kv̂dq,k k2 + cQ,k

and combining (4.3), (4.4) and (4.10), the overall microgrid
model (see Fig. 7) is given by the differential equations


Γ ⊗ Ip|N1 | ż = g(z, u),


VΓ ⊗ Im|N1 | ẋ = f (x, w(z, u), iabc ),


K ⊗ Ir|N` | ẋ` = f` (x` , iabc ),
[L ⊗ I2 ] x̂˙ L,dq = (−R ⊗ I2 + X )x̂L,dq + B > ⊗ I2 v̂dq ,

:= Q∗k (kv̂dq,k k2 ),

Assumption 5.3. L < I|E| in (4.10). Therefore,
x̂˙ L,dq (t) = 02|E| for all t ≥ 0.

together with the algebraic relations

Assumption 5.1 is equivalent to the assumption that
the inner current and voltage controllers track the voltage
and current references instantaneously and exactly. Usually, the current and voltage controllers in (4.1) (see also
Fig. 6a) are designed such that the resulting closed-loop
system (4.1) has a very large bandwidth compared to the
control system located at the next higher control level represented by (4.2) [65, 20, 74]. If this time-scale separation
is followed in the design of the system (4.3), the first part

iabc = B ⊗ I3 xL = B ⊗ Tdq (φ)> x̂L,dq ,

vabc,i = hi (xi , wi (zi , ui )),

(4.12)

vabc,k = h`,k (x`,k , iabc,k ),
i ∼ N1 ,



where aP,k , bP,k , cP,k , aQ,k , bQ,k and cQ,k are real constants and Pk (v̂dq,k , îdq,k ) and Qk (v̂dq,k , îdq,k ) are calculated as given in Definition 2.11.

(4.11)

v̂dq = I|N | ⊗ Tdq (φ)vabc ,

iabc,k
i ∼ N1

By defining the aggregated nodal voltages and currents in
dq-coordinates

v̂L,dq := col(v̂L,dq,l ) ∈ R2|E| ,

iabc,i

k ∼ N` .

11

of Assumption 5.1 can be mathematically formalized by invoking singular perturbation theory [56, Chapter 11], [58].
The second part of Assumption 5.1 expresses the fact that
the inner control system (4.1) is assumed to track the reference vref,i = wi (zi , ui ) exactly, independently of the disturbance iabc,i . Typical values for the bandwidth of (4.1)
reported in [74, 81] are in the range of 400 − 600 Hz, while
those of (4.2) are in the range of 2 − 10 Hz.
Assumption 5.2 implies that the dynamics of loads and
grid-feeding units can be neglected. This assumption is
also frequently employed in microgrid and power system
stability studies, where loads are often modeled as either
constant impedance (Z), constant current (I) or constant
power loads (P) or a combination of them (ZIP) [59, 103].
Similarly, grid-feeding units with positive active power injection are represented by setting aP,k = aQk = 0 and
bP,k or cP,k to negative values. The values for bQ,k and
cQ,k should be chosen in dependency of the reactive power
contribution of the unit.
Assumption 5.3 is standard in power system analysis
[59, 39, 87, 5, 69, 37]. The usual justification of Assumption 5.3 is that the line dynamics evolve on a much faster
time-scale than the dynamics of the generation sources. In
the present case, Assumption 5.3 is justified whenever Assumption 5.1 is employed, since the line dynamics (4.10)
are typically at least as fast as those of the internal inverter
controls (4.1), see, e.g., [81]. Again, Assumption 5.3 can
be mathematically formalized by invoking singular perturbation arguments [56, Chapter 11], [58].
Under Assumption 5.1, the model of each grid-forming
inverter (4.3) reduces to
γi żi = gi (zi , ui ),
vabc,i = wi (zi , ui ),

i ∼ N1 .

vref

vref

va

Inverter
with LC
filter and
inner
control
loops

va

Rg

Lg

vGa

vb

vGb

vc

vG c

Rg

Lg

vG a

vb

vGb

vc

vG c

Figure 8: Simplified representation of an inverter operated in gridforming mode as ideal controllable voltage source. Bold lines represent electrical connections, while dashed lines represent signal connections.

through a filter to obtain the power components corresponding to the fundamental frequency [81, 20, 74]
γi τPi Ṗim = −Pim + Pi ,

m
γi τPi Q̇m
i = −Qi + Qi .

(5.3)

Here, Pi and Qi are the active and reactive power injections of the inverter, Pim : R≥0 → R and Qm
i : R≥0 → R
their measured values and τPi ∈ R>0 is the time constant
of the low pass filter.
Note that whenever the particular form (5.2), (5.3) of
(4.2) is considered and the measured and filtered power
signals are used as feedback signals in the controls uδi , respectively uV
i , then the bandwidth of the overall control
system is limited by the bandwidth of the measurement
filter. Hence, if τPi  νi , then Assumption 5.1 is justified.
With Assumption 5.2, (4.4) can be represented by the
algebraic relation

(5.1)

The model (5.1) represents the inverter as an AC voltage
source, the amplitude and frequency of which can be defined by the designer. The system (5.1) is a very commonly
used model of a grid-forming inverter in microgrid control
design and analysis [65, 40, 55, 90].
Furthermore, often a particular structure of (5.1) is
used in the literature [95, 96, 88, 90, 2, 76]. As discussed
in Section 2.1, a symmetric three-phase voltage can be
completely described by its phase angle and its amplitude.
In addition, it is usually preferred to control the frequency
of the inverter output voltage, instead of the phase angle.
Hence, a suitable model of the inverter at the i-th node is
given by [90, 88]

Pk (v̂dq,k , îdq,k ) = Pk∗ (kv̂dq,k k2 ),

Qk (v̂dq,k , îdq,k ) = Q∗k (kv̂dq,k k2 ),

k ∼ N` .

(5.4)

Finally, under Assumption 5.3, the network model (4.10)
is also static and given by
−1

îdq = B ⊗ I2 (R ⊗ I2 − X )

γi α̇i = ωi = uδi ,
Vi = uV
i ,

Inverter
with LC
filter and
inner
control
loops

(5.2)

vabc,i = vabc,i (αi , Vi ),
where uδi : R≥0 → R and uVi : R≥0 → R are control signals.
Usually, it is also assumed that the active and reactive
power output of the inverter is measured and processed

(5.5)

The reduced-order model (5.2) - (5.5) is still rather
complex to handle, as the variables of (5.2) are expressed
in abc-coordinates, while those of (5.4) and (5.5) are expressed in common dq-coordinates. Therefore, a more
compact representation of (5.2) - (5.5) is derived in the
following. To this end, it is convenient to recall that αi
is the angle of the voltage at the i-th node with initial
condition α0i , i ∼ N and to define
Z t
δi := α0i +
(α̇i − ω com )dτ ∈ S, i ∼ N .
(5.6)
0

12

B > ⊗ I2 v̂dq .

vabc,i =

q

2
3 Vi




sin(αi )
2π
sin(αi − )
3
sin(αi + 2π
3 )

Tdq (αi )

vdq,i =

Note that the reactances Xl = Ll ω com are calculated at
the frequency ω com , which, under the made assumptions,
should be chosen as the (constant) synchronous frequency
of the network—denoted by ω s ∈ R in the following4 . Typically, ω s ∈ 2π[45, 65] rad/s.



 
Vd,i
0
= Vi
Vq,i
1

Tδ (δi )

Tdq (φ)

v̂dq,i

Remark 5.4. The form (5.12) is a very popular representation and these complex quantities are often denoted as
phasors [5, 112]. Furthermore, by using Euler’s formula
[47], (5.12) can also be rewritten in polar form. Note,
however, that, unlike, e.g., [5, 112], other authors define a
phasor as a complex sinusoidal quantity with a constant
frequency [37].





V̂
sin(αi − φ)
= d,i = Vi
cos(αi − φ)
V̂q,i

Figure 9: Illustration of the different coordinate frames used to
derive the model of an electrical network given in (5.15). The signal
vabc,i : R≥0 → R3 denotes the three-phase voltage at the i-th bus
with phase angle αi : R≥0 → S and amplitude Vi : R≥0 → R≥0 ,
i ∼ N . The mappings Tdq and Tδ are given in (2.3), respectively
(5.7). The angle δi : R≥0 → S is defined in (5.6). Note that, by
construction, αi − δi = mod2π (ω com t) = φ, where the real constant
ω com denotes the speed of the common dq-reference frame.

Define the admittance matrix of the electrical network
by

Gii :=<(Yii ),

i ∼ N,

Tdq (φ) = Tδ (δi )Tdq (αi ).

and

Hence, by construction,

Gii + jBii =

 
0
= Tδ (δi )Vi
,
1
(5.8)

i ∼ N.

(5.9)

Recall that Vqd and Iqd defined in (5.12) are expressed in
local dq-coordinates. By making use of (5.9) and (5.14),
(5.15) can be written component-wise as
Iqd,i = Iq,i + jId,i ,
X
Iq,i = Gii Vi −
(Gik cos(δik ) + Bik sin(δik )) Vk ,
k∼Ni

Id,i = Bii Vi −

(5.11)

and define
V̂qd :=V̂q + j V̂d ∈ C

Vqd :=Vq + jVd ∈ C

,

|N |

,

Iˆqd := Iˆq + j Iˆd ∈ C|N | ,
Iqd := Iq + jId ∈ C|N | .

X

k∼Ni

(Bik cos(δik ) − Gik sin(δik )) Vk ,
(5.16)

i ∼ N , where, for ease of notation, angle differences are
written as δik := δi − δk . Furthermore, from Definition 2.11

(5.12)

4 Under the made assumptions, (5.5) is the equilibrium of the
”fast” line dynamics (4.9) [56, Chapter 11]. Hence, in order for the
currents îdq and voltages v̂dq to be constant in steady-state, ω com
has to be chosen identically to the synchronous steady-state network
frequency.

Then, with X := diag(Xl ) = diag(Ll ω com ) ∈ R|E|×|E| ,
we can rewrite (5.5) as
−1
Iˆqd = B (R + jX) B > V̂qd .

(Rl + jXl )−1 ,

where Ei denotes the set of edges associated to node i.
Inserting (5.10) and (5.11) into (5.13) yields


Iqd = diag e−jδi Ydiag ejδi Vqd .
(5.15)

V̂qd,i := V̂q,i + j V̂d,i = (cos(δi ) + j sin(δi ))Vqd,i = ejδi Vqd,i ,
(5.10)
where Vqd,i = Vq,i + jVd,i , i ∼ N . Equivalently, let

|N |

X

l∼Ei

The variables vdq,i are referred to as local dq-coordinates
of vabc,i in the following. The relation between vabc,i , v̂dq,i
and vdq,i is illustrated in Fig. 9.
It is convenient to represent (5.8) in the complex plane

Iˆqd,i := Iˆq,i + j Iˆd,i = ejδi Iqd,i

i 6= k.

Moreover, it follows immediately that

if nodes i and k are



0
not connected
Yik =

if nodes i and k are

−(Rl + jXl )−1

connected by line l

and that straightforward algebraic manipulations yield

which makes it convenient to define


 
Vd,i
0
vdq,i :=
= Vi
,
Vq,i
1

(5.14)

Bii := =(Yii ),

Yik :=Gik + jBik := −Yik ,

Note that, with δi defined in (5.6),

v̂dq,i = Tdq (φ)vabc,i = Tδ (δi )Tdq (αi )vabc,i

B > ∈ C|N |×|N |

and

Let $ : R≥0 → S and consider the mapping Tδ : S → R2×2 ,


cos($) sin($)
Tδ ($) :=
.
(5.7)
− sin($) cos($)
αi − δi = mod2π (ω com t) = φ,

−1

Y := B (R + jX)

(5.13)
13

Network

together with (5.9) and (5.16), the power flows in the network are given by
Pi =Gii Vi2 −

Pi = Vi Iq,i =
X
Gii Vi2 −
(Gik cos(δik ) + Bik sin(δik )) Vk Vi ,

X

=

−Bii Vi2

+

X

Qi = − Bii Vi2 +

k∼Ni

(δi , Vi )

X

k∼Ni

(Bik cos(δik ) − Gik sin(δik )) Vk Vi

(Pi , Qi )
i ∼ N1

(Bik cos(δik ) − Gik sin(δik )) Vk Vi .

(Pk , Qk )

0 = Pk − Pk∗ (Vk ),

γi τPi Ṗim = −Pim + Pi ,
Vi =

The equations (5.17) are the standard power flow equations used in most recent work on microgrid control design
and stability analysis, e.g., [95, 90, 2, 76].

0 = Qk − Q∗k (Vk )

uV
i ,

m
γi τPi Q̇m
i = −Qi + Qi

Grid-forming inverters

(δk , Vk )

k ∼ N`

γi δ̇i = uδi − ω com ,

(5.17)

Grid-feeding inverters and loads

Figure 10: Reduced microgrid model in standard representation
with grid-forming inverters modeled by (5.18), as well as loads and
grid-feeding inverters given by (5.4) with kv̂dq,k k2 = Vk . The network
is represented by the power flow equations (5.17). This compact
model is obtained from the model (5.2) - (5.5) by expressing the
current and power flows in local dq-coordinates vdq,i and idq,i , i ∼ N ,
see (5.15) and Fig. 9.

Remark 5.5. Note that for any other choices of the transformation angle in local dq-coordinates Vdi 6= 0. This is
usually the case when modeling SGs, since the angle of
the internal machine electromagnetic force (EMF) is in
general not known. Then, the equations (5.17) become
slightly more involved, see [5, Chapter 9].

Then the system (5.2) - (5.5) can be written equivalently
by means of (5.4), (5.17), (5.18) as

Furthermore, in local dq-coordinates, the particular inverter model (5.2), (5.3), is given by

Γδ˙I = uδ − ω com 1|N1 | ,

γi δ˙i = ωi − ω com = uδi − ω com ,

γi τPi Ṗim = −Pim + Pi ,

ΓT Ṗ m = −P m + PI ,

(5.18)

Vi = uVi ,

VI = uV ,

ΓT Q̇m = −Qm + QI ,

m
γi τPi Q̇m
i = −Qi + Qi ,

(5.19)

0|N` | = P` − P`∗ ,

with Vqd,i = Vi (see (5.9)) and Pi and Qi given by (5.17).
Finally, recall (5.4) and note that
kv̂dq,k k2 = kV̂qd,k k2 = kVqd,k k2 = Vk ,

(Gik cos(δik ) + Bik sin(δik )) Vk Vi ,

k∼Ni

k∼Ni

Qi = −Vi Id,i

i∼N

0|N` | = Q` − Q∗` ,
where the last 2|N` | algebraic equations correspond to the
power balances at nodes k ∼ N` .
This section has illustrated the main modeling steps
and assumptions, which lead from the detailed microgrid
model (4.11), (4.12) to the model (5.19), (5.17), respectively (5.4), (5.17), (5.18). The model (5.19), (5.17) is
frequently used in the analysis and control design of microgrids [95, 96, 14, 89, 90, 2, 76, 92, 93]. Some of the mentioned work is conducted under additional assumptions
such as instantaneous power measurements [95, 2, 76], constant voltage amplitudes [95, 14, 2, 89] or small phase angle
differences [96, 92, 93]. In addition, ideal clocks are usually
assumed, i.e., Γ = I|N1 | . Furthermore, whenever constant
impedance or constant current loads are assumed, the algebraic equations in (5.19), (5.17) can be eliminated by an
appropriate network reduction. This process is commonly
known as Kron reduction and frequently employed in microgrid and power system studies. For further details on
Kron reduction, the reader is referred to [59, 28].

k ∈ N` .

This completes the reformulation of the model (5.2) - (5.5).
The final overall microgrid model is given by (5.4), (5.17),
(5.18) and shown in Fig. 10. This is the standard model
employed throughout the literature.
The section is concluded by deriving a vector-based
formulation of the microgrid model (5.4), (5.17), (5.18).
To this end, we define the vectors
|N |

δI :=col(δi ) ∈ S|N1 | , VI := col(Vi ) ∈ R>01 ,

uδ :=col(uδi ) ∈ R|N1 | , uV := col(uVi ) ∈ R|N1 | ,

PI :=col(Pi ) ∈ R|N1 | , QI := col(Qi ) ∈ R|N1 | ,

P` :=col(Pk ) ∈ R|N` | , Q` := col(Qk ) ∈ R|N` | ,

P`∗ :=col(Pk∗ (Vk )) ∈ R|N` | , Q∗` := col(Q∗k (Vk )) ∈ R|N` | ,
with Pi , Pk , Qi , and Qk given by (5.17), as well as the
matrix
T := diag(τPi ) ∈ R|N1 |×|N1 | .
14

voltage stability [95, 96, 90, 102, 25, 2, 76, 92, 93, 29],
secondary control [10, 95, 11, 29, 94] or optimal dispatch
[29, 13, 44]. Compared to the model derived in the present
paper, most of the aforementioned work is conducted under certain additional simplifying assumptions, such as
constant voltage amplitudes, identical DG unit dynamics,
as well as lossless or identical lines.
On the modeling side, one direct extension of the presented modeling framework is to investigate the suitability of the use of more refined modeling techniques such
as dynamic phasors [98, 70, 21] or symmetric components
[77, 37], to describe the dynamics of a microgrid in asymmetric operating conditions. In this context, one main
challenge is to accurately consider these phenomena while
maintaining an analytically tractable model. Another subject of future research is the derivation of more detailed
load models for microgrid applications. As in conventional
power systems, accurate load modeling is a very important, but also very difficult task [59, 103]. The main reason for this is that there are typically many different kinds
of loads connected within one power system or microgrid,
see, e.g., [59, Chapter 7]. As a consequence, it is difficult
to obtain suitable generically valid abstractions.
In conclusion, there are many challenging open research
questions regarding a reliable, safe and efficient operation of microgrids. Therefore, the authors hope that the
present survey on modeling of microgrids may serve as a
stimulating base for a large variety of future research on
both the theoretical and the application side.

6. Conclusions and topics of future research
6.1. Summary
The present survey paper has introduced the reader to
the microgrid concept with the main focus of providing
a detailed procedure for the model derivation of a threephase inverter-based microgrid. In particular, it has been
shown how—and under which assumptions—the microgrid
models usually used in the literature can be obtained from
a significantly more complex model derived from fundamental physical laws. The assumptions invoked in the
reduction process are often satisfied in standard applications. Therefore, the reduced model represents a valid approximation and may, hence, be useful for control design
and system analysis. In addition, the employed model reduction techniques can equivalently be applied to standard
bulk power system models.
Nevertheless, it is important to note that the model
derived in the present paper neglects effects such as asymmetric operation, DC-side dynamics of DG units or line
capacitances. These facts have to be kept in mind, when
performing microgrid analysis based on the derived model
and assessing the results.
Also, it is worth mentioning that numerical simulation
of the introduced microgrid models (4.11), (4.12), respectively (5.19), (5.17), requires careful selection of the numerical integration method to be employed. The main
reason for this is that the model (4.11), (4.12) contains
dynamics evolving at a wide range of time-scales, i.e., it is
a stiff model [49, Chapter 7], [33, Chapter 8]. As is wellknown, certain (standard) numerical integration methods
will lead to numerical instability, when applied to stiff
models—unless an extremely small step size is employed
[49, Chapter 7], [33, Chapter 8]. On the contrary, in the
reduced model (5.19), (5.17) the fast dynamics have been
eliminated and replaced by their corresponding steadystate equations. Hence, this model is not stiff and simpler
integration methods can be used. Very similar situations
are usually encountered in simulation of large conventional
power systems [69, Chapter 13], [21].

References
[1] Ackermann, T., Andersson, G., & Söder, L. (2001). Distributed generation: a definition. Electric power systems research, 57 , 195–204.
[2] Ainsworth, N., & Grijalva, S. (2013). A structure-preserving
model and sufficient condition for frequency synchronization of
lossless droop inverter-based AC networks. IEEE Transactions
on Power Systems, 28 , 4310–4319.
[3] Akagi, H., Kanazawa, Y., & Nabae, A. (1983). Generalized
theory of the instantaneous reactive power in three-phase circuits. In International Power Electronic Conference (IPEC)
(pp. 1375–1386). volume 83.
[4] Akagi, H., Watanabe, E. H., & Aredes, M. (2007). Instantaneous power theory and applications to power conditioning.
John Wiley & Sons.
[5] Anderson, P., & Fouad, A. (2002). Power system control and
stability. J.Wiley & Sons.
[6] Andersson, G. (2012). Dynamics and control of electric power
systems. Lecture notes, EEH - Power Systems Laboratory,
ETH Zürich, .
[7] Aredes, M., Akagi, H., Watanabe, E. H., Vergara Salgado, E.,
& Encarnacao, L. F. (2009). Comparisons between the p–q
and p–q–r theories in three-phase four-wire systems. IEEE
Transactions on Power Electronic, 24 , 924–933.
[8] Barnes, M., Kondoh, J., Asano, H., Oyarzabal, J., Ventakaramanan, G., Lasseter, R., Hatziargyriou, N., & Green, T.
(2007). Real-world microgrids-an overview. In IEEE Int. Conf.
on System of Systems Engineering, 2007. SoSE ’07 (pp. 1 –8).
[9] Beck, H.-P., & Hesse, R. (2007). Virtual synchronous machine.
In 9th Int. Conf. on Electr. Power Quality and Utilisation.
(pp. 1 –6).

6.2. Future research
To conclude this survey, we very briefly highlight some
topics of future research. Following up on the discussion
in the introduction, there are numerous challenges related
to system and control theory in microgrid applications.
To further motivate these, we briefly review control goals
in microgrids. At the present, the following are considered to be among the most relevant control objectives in
microgrids [32, 61, 40, 46, 55, 37, 41]: frequency stability, voltage stability, operational compatibility of inverterinterfaced and SG-interfaced units, desired power sharing
in steady-state, seamless switching from grid-connected to
islanded-mode and vice-versa, robustness with respect to
uncertainties and optimal dispatch.
Some of these problems have been addressed in recent
work within the control community, e.g., frequency and
15

[10] Bidram, A., Davoudi, A., Lewis, F. L., & Qu, Z. (2013). Secondary control of microgrids based on distributed cooperative
control of multi-agent systems. IET Generation, Transmission
& Distribution, 7 , 822–831.
[11] Bidram, A., Lewis, F., & Davoudi, A. (2014). Distributed control systems for small-scale power networks: Using multiagent
cooperative control theory. IEEE Control Systems Magazine,
34 , 56–77.
[12] Blaabjerg, F., Teodorescu, R., Liserre, M., & Timbus, A. V.
(2006). Overview of control and grid synchronization for distributed power generation systems. IEEE Transactions on Industrial Electronics, 53 , 1398–1409.
[13] Bolognani, S., & Zampieri, S. (2013). A distributed control
strategy for reactive power compensation in smart microgrids.
IEEE Transactions on Automatic Control, 58 , 2818–2833.
[14] Bouattour, H., Simpson-Porco, J. W., Dörfler, F., & Bullo,
F. (2013). Further results on distributed secondary control in
microgrids. In 52nd Conference on Decision and Control (pp.
1514–1519).
[15] Buchholz, F. (1950).
Das Begriffssystem Rechtleistung,
Wirkleistung, totale Blindleistung. Selbstverl.;[Lachner in
Komm].
[16] Carrasco, J. M., Franquelo, L. G., Bialasiewicz, J. T., Galván,
E., Guisado, R. C. P., Prats, M. Á. M., León, J. I., & MorenoAlfonso, N. (2006). Power-electronic systems for the grid integration of renewable energy sources: A survey. IEEE Transactions on Industrial Electronics, 53 , 1002–1016.
[17] Chandorkar, M., Divan, D., & Adapa, R. (1993). Control of
parallel connected inverters in standalone AC supply systems.
IEEE Transactions on Industry Applications, 29 , 136 –143.
[18] Chiniforoosh, S., Jatskevich, J., Yazdani, A., Sood, V., Dinavahi, V., Martinez, J., & Ramirez, A. (2010). Definitions
and applications of dynamic average models for analysis of
power systems. IEEE Transactions on Power Delivery, 25 ,
2655–2669.
[19] Chowdhury, S., & Crossley, P. (2009). Microgrids and active distribution networks. The Institution of Engineering and
Technology.
[20] Coelho, E., Cortizo, P., & Garcia, P. (2002). Small-signal stability for parallel-connected inverters in stand-alone AC supply
systems. IEEE Transactions on Industry Applications, 38 , 533
–542.
[21] Demiray, T. H. (2008). Simulation of power system dynamics
using dynamic phasor models. Ph.D. thesis ETH Zurich.
[22] Depenbrock, M. (1962). Untersuchungen über die Spannungsund Leistungsverhältnisse bei Umrichtern ohne Energiespeicher . Ph.D. thesis Mikrokopie GmbH.
[23] Depenbrock, M. (1993). The FBD-method, a generally applicable tool for analyzing power relations. IEEE Transactions
on Power Systems, 8 , 381–387.
[24] Depenbrock, M., Staudt, V., & Wrede, H. (2004). Concerning”
instantaneous power compensation in three-phase systems by
using pqr theory”. IEEE Transactions on Power Electronics,
19 , 1151–1152.
[25] Dhople, S. V., Johnson, B. B., Dörfler, F., & Hamadeh, A. O.
(2014). Synchronization of nonlinear circuits in dynamic electrical networks with general topologies. IEEE Transactions on
Circuits and Systems I: Regular Papers, 61 , 2677–2690.
[26] Diestel, R. (2000). Graduate texts in mathematics: Graph
theory. Springer.
[27] DIN 40110-1:1994 (1994).
DIN 40110-1:1994 - Wechselstromgrößen; Zweileiter-Stromkreise, (engl: ”Quantities used
in alternating current theory; two-line circuits).
[28] Dörfler, F., & Bullo, F. (2013). Kron reduction of graphs with
applications to electrical networks. IEEE Transactions on Circuits and Systems I: Regular Papers, 60 , 150–163.
[29] Dörfler, F., Bullo, F. et al. (2016). Breaking the hierarchy: Distributed control & economic optimality in microgrids. IEEE
Transactions on Control of Network Systems, to appear., .
[30] Emanuel, A. E. (2004). Summary of IEEE standard 1459: definitions for the measurement of electric power quantities under

[31]
[32]
[33]
[34]

[35]

[36]

[37]
[38]
[39]
[40]

[41]

[42]

[43]

[44]

[45]

[46]

[47]
[48]

[49]
[50]

[51]

[52]

[53]

16

sinusoidal, nonsinusoidal, balanced, or unbalanced conditions.
IEEE Transactions on Industry Applications, 40 , 869–876.
Erickson, R. W., & Maksimovic, D. (2001). Fundamentals of
power electronics. Springer.
Farhangi, H. (2010). The path of the smart grid. IEEE Power
and Energy Magazine, 8 , 18 –28.
Fatunla, S. O. (2014). Numerical methods for initial value
problems in ordinary differential equations. Academic Press.
Fiaz, S., Zonetti, D., Ortega, R., Scherpen, J., & van der
Schaft, A. (2013). A port-hamiltonian approach to power network modeling and analysis. European Journal of Control, 19 ,
477–485.
Fryze, S. (1932). Wirk-, Blind-und Scheinleistung in elektrischen Stromkreisen mit nichtsinusförmigem Verlauf von
Strom und Spannung. Elektrotechnische Zeitschrift, 25 , 569–
599.
Fukuda, S., & Yoda, T. (2001). A novel current-tracking
method for active filters based on a sinusoidal internal model
[for pwm invertors]. IEEE Transactions on Industry Applications, 37 , 888–895.
Glover, J. D., Sarma, M. S., & Overbye, T. J. (2011). Power
system analysis and design. Cengage Learning.
Godsil, C., & Royle, G. (2001). Algebraic graph theory.
Springer.
Grainger, J. J., & Stevenson, W. D. (1994). Power system
analysis volume 621. McGraw-Hill New York.
Green, T., & Prodanovic, M. (2007). Control of inverter-based
micro-grids. Electric Power Systems Research, Vol. 77 , 1204–
1213.
Guerrero, J., Loh, P., Chandorkar, M., & Lee, T. (2013). Advanced control architectures for intelligent microgrids – part I:
Decentralized and hierarchical control. IEEE Transactions on
Industrial Electronics, 60 , 1254–1262.
Guerrero, J., Vasquez, J., Matas, J., de Vicuna, L., & Castilla,
M. (2011). Hierarchical control of droop-controlled AC and DC
microgrids; a general approach toward standardization. IEEE
Transactions on Industrial Electronics, 58 , 158 –172.
Guo, Y., & Gawlik, W. (2014). A survey of control strategies
applied in worldwide microgrid projects. Tagungsband ComForEn 2014 , (p. 47).
Hans, C. A., Nenchev, V., Raisch, J., & Reincke-Collon, C.
(2014). Minimax model predictive operation control of microgrids. In 19th IFAC World Congress (p. 1028710292). Cape
Town, South Africa.
Hansen, J., Sato, M., Ruedy, R., Nazarenko, L., Lacis, A.,
Schmidt, G., Russell, G., Aleinov, I., Bauer, M., Bauer, S. et al.
(2005). Efficacy of climate forcings. Journal of Geophysical
Research: Atmospheres (1984–2012), 110 .
Hatziargyriou, N., Asano, H., Iravani, R., & Marnay, C.
(2007). Microgrids. IEEE Power and Energy Magazine, 5 ,
78–94.
Hazewinkel, M. (1993). Encyclopaedia of Mathematics (9) volume 9. Springer.
Heuck, K., Dettmann, K.-D., & Schulz, D. (2013). Elektrische
Energieversorgung: Erzeugung, Übertragung und Verteilung
elektrischer Energie für Studium und Praxis. Springer DE.
Hoffman, J. D., & Frankel, S. (2001). Numerical methods for
engineers and scientists. CRC press.
Hornik, T., & Zhong, Q.-C. (2011). A current-control strategy
for voltage-source inverters in microgrids based on and repetitive control. IEEE Transactions on Power Electronics, 26 ,
943–952.
Houghton, J. T. (1996). Climate change 1995: The science of
climate change: contribution of working group I to the second
assessment report of the Intergovernmental Panel on Climate
Change volume 2. Cambridge University Press.
Justo, J. J., Mwasilu, F., Lee, J., & Jung, J.-W. (2013). Acmicrogrids versus dc-microgrids with distributed energy resources: A review. Renewable and Sustainable Energy Reviews, 24 , 387–405.
Katiraei, F., & Iravani, M. (2006). Power management strate-

[54]

[55]

[56]
[57]

[58]

[59]
[60]

[61]
[62]
[63]

[64]

[65]

[66]
[67]
[68]

[69]
[70]

[71]

[72]

[73]

[74]

[75]
[76]

gies for a microgrid with multiple distributed generation units.
IEEE Transactions on Power Systems, 21 , 1821 –1831.
Katiraei, F., Iravani, M., & Lehn, P. (2007). Small-signal dynamic model of a micro-grid including conventional and electronically interfaced distributed resources. IET Generation,
Transmission Distribution, 1 , 369–378.
Katiraei, F., Iravani, R., Hatziargyriou, N., & Dimeas, A.
(2008). Microgrids management. IEEE Power and Energy
Magazine, 6 , 54–65.
Khalil, H. K. (2002). Nonlinear systems volume 3. Prentice
Hall.
Kim, H., & Akagi, H. (1999). The instantaneous power theory
on the rotating pqr reference frames. In International Conference on Power Electronics and Drive Systems (PEDS) (pp.
422–427). IEEE volume 1.
Kokotovic, P., Khali, H. K., & O’reilly, J. (1999). Singular perturbation methods in control: analysis and design volume 25.
Siam.
Kundur, P. (1994). Power system stability and control.
McGraw-Hill.
Kwasinski, A., & Onwuchekwa, C. N. (2011). Dynamic behavior and stabilization of DC microgrids with instantaneous
constant-power loads. IEEE Transactions on Power Electronics, 26 , 822–834.
Lasseter, R. (2002). Microgrids. In IEEE Power Engineering
Society Winter Meeting, 2002 (pp. 305 – 308 vol.1). volume 1.
Lasseter, R. H. (2011). Smart distribution: Coupled microgrids. Proceedings of the IEEE , 99 , 1074–1082.
Lidula, N., & Rajapakse, A. (2011). Microgrids research: A
review of experimental microgrids and test systems. Renewable
and Sustainable Energy Reviews, 15 , 186–202.
Loh, P. C., & Holmes, D. G. (2005). Analysis of multiloop
control strategies for LC/CL/LCL-filtered voltage-source and
current-source inverters. IEEE Transactions on Industry Applications, 41 , 644–654.
Lopes, J., Moreira, C., & Madureira, A. (2006). Defining control strategies for microgrids islanded operation. IEEE Transactions on Power Systems, 21 , 916 – 924.
Lund, H. (2007). Renewable energy strategies for sustainable
development. Energy, 32 , 912–919.
Lund, H. (2009). Renewable energy systems: the choice and
modeling of 100% renewable solutions. Academic Press.
Luo, L., & Dhople, S. (2014). Spatiotemporal model reduction
of inverter-based islanded microgrids. IEEE Transactions on
Energy Conversion, 29 , 823–832.
Machowski, J., Bialek, J., & Bumby, J. (2008). Power system
dynamics: stability and control. J.Wiley & Sons.
Maksimović, D., Stanković, A. M., Thottuvelil, V. J., & Verghese, G. C. (2001). Modeling and simulation of power electronic
converters. Proceedings of the IEEE , 89 , 898–912.
Mariani, V., Vasca, F., & Guerrero, J. M. (2014). Analysis of
droop controlled parallel inverters in islanded microgrids. In
IEEE International Energy Conference (ENERGYCON) (pp.
1304–1309). IEEE.
Mariani, V., Vasca, F., & Guerrero, J. M. (2014). Dynamicphasor-based nonlinear modelling of AC islanded microgrids
under droop control. In 11th International Multi-Conference
on Systems, Signals & Devices (SSD) (pp. 1–6). IEEE.
Miao, Z., Domijan, A., & Fan, L. (2011). Investigation of microgrids with both inverter interfaced and direct AC-connected
distributed energy resources. IEEE Transactions on Power
Delivery, 26 , 1634 –1642.
Mohamed, Y., & El-Saadany, E. (2008). Adaptive decentralized droop controller to preserve power sharing stability of paralleled inverters in distributed generation microgrids. IEEE
Transactions on Power Electronics, 23 , 2806 –2816.
Mohan, N., & Undeland, T. M. (2007). Power electronics:
converters, applications, and design. John Wiley & Sons.
Münz, U., & Metzger, M. (2014). Voltage and angle stability
reserve of power systems with renewable generation. In 19th
IFAC World Congress (pp. 9075–9080). Cape Town, South

Africa.
[77] Paap, G. C. (2000). Symmetrical components in the time
domain and their application to power network calculations.
IEEE Transactions on Power Systems, 15 , 522–528.
[78] Park, R. H. (1929). Two-reaction theory of synchronous machines generalized method of analysis-part I. Transactions of
the American Institute of Electrical Engineers, 48 , 716–727.
[79] Peng, F. Z., & Lai, J.-S. (1996). Generalized instantaneous
reactive power theory for three-phase power systems. IEEE
Transactions on Instrumentation and Measurement, 45 , 293–
297.
[80] Planas, E., Gil-de Muro, A., Andreu, J., Kortabarria, I., &
Martı́nez de Alegrı́a, I. (2013). General aspects, hierarchical
controls and droop methods in microgrids: A review. Renewable and Sustainable Energy Reviews, 17 , 147–159.
[81] Pogaku, N., Prodanovic, M., & Green, T. (2007). Modeling,
analysis and testing of autonomous operation of an inverterbased microgrid. IEEE Transactions on Power Electronics,
22 , 613 –625.
[82] Prodanovic, M., & Green, T. C. (2003). Control and filter design of three-phase inverters for high power quality grid connection. IEEE Transactions on Power Electronics, 18 , 373–380.
[83] Riverso, S., Sarzo, F., & Ferrari-Trecate, G. (2015). Plugand-play voltage and frequency control of islanded microgrids
with meshed topology. IEEE Transactions on Smart Grid, 6 ,
1176–1184.
[84] Rocabert, J., Luna, A., Blaabjerg, F., & Rodriguez, P. (2012).
Control of power converters in AC microgrids. IEEE Transactions on Power Electronics, 27 , 4734–4749.
[85] Salomonsson, D. (2008). Modeling, Control and Protection
of Low-Voltage DC Microgrids. Ph.D. thesis KTH, Electric
Power Systems.
[86] Salomonsson, D., Söder, L., & Sannino, A. (2009). Protection
of low-voltage dc microgrids. IEEE Transactions on Power
Delivery, 24 , 1045–1053.
[87] Sauer, P., & Pai, M. (1998). Power system dynamics and
stability. Prentice Hall.
[88] Schiffer, J., Anta, A., Trung, T. D., Raisch, J., & Sezi, T.
(2012). On power sharing and stability in autonomous inverterbased microgrids. In 51st Conference on Decision and Control
(pp. 1105–1110). Maui, HI, USA.
[89] Schiffer, J., Goldin, D., Raisch, J., & Sezi, T. (2013). Synchronization of droop-controlled microgrids with distributed
rotational and electronic generation. In 52nd Conference on
Decision and Control (pp. 2334–2339). Florence, Italy.
[90] Schiffer, J., Ortega, R., Astolfi, A., Raisch, J., & Sezi, T.
(2014). Conditions for stability of droop-controlled inverterbased microgrids. Automatica, 50 , 2457–2469.
[91] Schiffer, J., Ortega, R., Hans, C., & Raisch, J. (2015). Droopcontrolled inverter-based microgrids are robust to clock drifts.
In American Control Conference (pp. 2341–2346).
[92] Schiffer, J., Seel, T., Raisch, J., & Sezi, T. (2014). A consensusbased distributed voltage control for reactive power sharing in
microgrids. In 13th European Control Conference (pp. 1299–
1305). Strasbourg, France.
[93] Schiffer, J., Seel, T., Raisch, J., & Sezi, T. (2016). Voltage
stability and reactive power sharing in inverter-based microgrids with consensus-based distributed voltage control. IEEE
Transactions on Control Systems Technology, 24 , 96–109.
[94] Shafiee, Q., Guerrero, J., & Vasquez, J. (2014). Distributed
secondary control for islanded microgrids – a novel approach.
IEEE Transactions on Power Electronics, 29 , 1018–1031.
[95] Simpson-Porco, J. W., Dörfler, F., & Bullo, F. (2013). Synchronization and power sharing for droop-controlled inverters
in islanded microgrids. Automatica, 49 , 2603 – 2611.
[96] Simpson-Porco, J. W., Dörfler, F., & Bullo, F. (2013). Voltage
stabilization in microgrids using quadratic droop control. In
52nd Conference on Decision and Control (pp. 7582–7589).
Florence, Italy.
[97] Solomon, S. (2007). Climate change 2007-the physical science
basis: Working group I contribution to the fourth assessment

17

report of the IPCC volume 4. Cambridge University Press.
[98] Stanković, A. M., & Aydin, T. (2000). Analysis of asymmetrical faults in power systems using dynamic phasors. IEEE
Transactions on Power Systems, 15 , 1062–1068.
[99] Strbac, G., Hatziargyriou, N., Lopes, J. P., Moreira, C.,
Dimeas, A., & Papadaskalopoulos, D. (2015). Microgrids: Enhancing the resilience of the European megagrid. IEEE Power
and Energy Magazine, 13 , 35–43.
[100] Teodorescu, R., Blaabjerg, F., Liserre, M., & Loh, P. C.
(2006). Proportional-resonant controllers and filters for gridconnected voltage-source converters. IEE Proceedings-Electric
Power Applications, 153 , 750–762.
[101] Teodorescu, R., Liserre, M., & Rodriguez, P. (2011). Grid converters for photovoltaic and wind power systems volume 29.
John Wiley & Sons.
[102] Tôrres, L. A., Hespanha, J. P., & Moehlis, J. (2015). Synchronization of identical oscillators coupled through a symmetric
network with dynamics: A constructive approach with applications to parallel operation of inverters. IEEE Transactions
on Automatic Control, 60 , 3226–3241.
[103] Van Cutsem, T., & Vournas, C. (1998). Voltage stability of
electric power systems volume 441. Springer.
[104] Varaiya, P. P., Wu, F. F., & Bialek, J. W. (2011). Smart
operation of smart grid: Risk-limiting dispatch. Proceedings
of the IEEE , 99 , 40–57.
[105] Venkatasubramanian, V., Schattler, H., & Zaborszky, J.
(1995). Fast time-varying phasor analysis in the balanced
three-phase large electric power system. IEEE Transactions
on Automatic Control, 40 , 1975–1982.
[106] Visscher, K., & De Haan, S. (2008). Virtual synchronous machines for frequency stabilisation in future grids with a significant share of decentralized generation. In SmartGrids for
Distribution, 2008. IET-CIRED. CIRED Seminar (pp. 1 –4).
[107] Wang, X., Guerrero, J. M., Blaabjerg, F., & Chen, Z. (2012).
A review of power electronics based microgrids. Journal of
Power Electronics, 12 , 181–192.
[108] Weiss, G., Zhong, Q.-C., Green, T. C., & Liang, J. (2004). H∞
repetitive control of DC-AC converters in microgrids. IEEE
Transactions on Power Electronics, 19 , 219–230.
[109] Willems, J. L. (1992). A new interpretation of the akaginabae power components for nonsinusoidal three-phase situations. IEEE Transactions on Instrumentation and Measurement, 41 , 523–527.
[110] Willems, J. L., Ghijselen, J. A., & Emanuel, A. E. (2005). The
apparent power concept and the IEEE standard 1459-2000.
IEEE Transactions on Power Delivery, 20 , 876–884.
[111] Zhong, Q., & Weiss, G. (2011). Synchronverters: Inverters
that mimic synchronous generators. IEEE Transactions on
Industrial Electronics, 58 , 1259 –1267.
[112] Zhong, Q.-C., & Hornik, T. (2012). Control of power inverters
in renewable energy and smart grid integration. John Wiley
& Sons.

18

