VIRTUAL OSCILLATOR CONTROL

Abstract

Virtual oscillator control systems, devices, and techniques are provided. One example device includes a processor configured to implement a virtual oscillator circuit and output an oscillating waveform based on the virtual oscillator circuit and power electronics operatively coupled to the processor and configured to convert, based on the oscillating waveform, direct current (DC) electricity to alternating current (AC) electricity. The processor may be further configured to extract, from the virtual oscillator circuit, a virtual current based on an output current of the AC electricity, and output the oscillating waveform further based on an input voltage of the DC electricity.

Claims

1. A device comprising: a processor configured to implement a virtual oscillator circuit and output an oscillating waveform based on the virtual oscillator circuit; and power electronics operatively coupled to the processor and configured to convert, based on the oscillating waveform, direct current (DC) electricity to alternating current (AC) electricity, wherein the processor is further configured to: extract, from the virtual oscillator circuit, a virtual current based on an output current of the AC electricity, and output the oscillating waveform further based on an input voltage of the DC electricity.

2. The device of claim 1, wherein the virtual oscillator circuit comprises a virtual negative-conductance element.

3. The device of claim 2, wherein: the virtual oscillator circuit further comprises a virtual inductor and a virtual capacitor, and the virtual negative-conductance element, the virtual inductor, and the virtual capacitor are connected in parallel.

4. The device of claim 1, wherein the virtual oscillator circuit comprises a virtual cubic voltage-dependent current source.

5. The device of claim 4, wherein: the virtual oscillator circuit further comprises a virtual inductor and a virtual capacitor, and the virtual cubic voltage-dependent current source, the virtual inductor, and the virtual capacitor are connected in parallel.

6. The device of claim 5, wherein the virtual cubic voltage-dependent current source consumes current proportional to a virtual voltage across the virtual capacitor.

7. The device of claim 5, wherein the virtual oscillator circuit further comprises a virtual negative-conductance element that is connected in parallel with the virtual cubic voltage-dependent current source, the virtual inductor, and the virtual capacitor.

8. The device of claim 1, wherein the processor is configured to output the oscillating waveform by: receiving, from the virtual oscillator circuit, a virtual output voltage; scaling the virtual output voltage based on the input voltage of the DC electricity to produce a scaled virtual output voltage; and comparing the scaled virtual output voltage to a carrier wave to produce the oscillating waveform.

9. The device of claim 1, wherein: the virtual oscillator circuit comprises a virtual inductor and a virtual capacitor, and the processor is further configured to: compute a linear combination of a virtual voltage across the virtual capacitor and a virtual current through the virtual inductor to produce an alternate virtual output voltage; and output the oscillating waveform further based on the alternate virtual output voltage.

10. The device of claim 9, wherein the linear combination is computed as κ.sub.v.sub.c cos φ+κ.sub.vεi.sub.L sin φ, wherein: κ.sub.v represents a voltage scaling factor, .sub.C represents the virtual voltage across the virtual capacitor, φ represents a rotation angle, ε represents a current scaling factor, and i.sub.L represents the virtual current through the virtual inductor.

11. The device of claim 1, further comprising an output filter that includes a filter inductor and that has a resistance, wherein: the virtual oscillator circuit comprises a virtual capacitor, and the processor is configured to output the oscillating waveform by: receiving, from the virtual oscillator circuit, a virtual output voltage; scaling the virtual output voltage based on a value of the resistance and an inductance value of the filter inductor to produce a first scaled virtual output voltage; combining the first scaled virtual output voltage with a virtual capacitor current flowing through the virtual capacitor to produce a combined virtual output voltage; and outputting the oscillating waveform further based on the combined virtual output voltage.

12. The device of claim 1, wherein: the virtual oscillator circuit comprises a virtual inductor and a virtual cubic voltage-dependent current source, and the processor is further configured to implement a power tracking control unit that modifies at least one of a virtual inductance of the virtual inductor or a scaling factor of the virtual cubic voltage-dependent current source based on a specified output power.

13. The device of claim 12, wherein the power tracking control unit modifies the at least one of the virtual inductance or the scaling factor based further on a scaled virtual output voltage of the virtual oscillator circuit and the output current of the AC electricity.

14. The device of claim 1, wherein the power electronics are configured to convert the DC electricity to AC electricity by alternatively engaging and disengaging power electronics switches based on the oscillating waveform.

15. A device comprising: a processor configured to: implement a virtual oscillator circuit comprising a virtual capacitor, a virtual inductor, a virtual negative-conductance element, and a virtual cubic voltage dependent current source; generate, based on the virtual oscillator circuit, an oscillating waveform; and output, based on the oscillating waveform, at least one control signal to cause power electronics to convert direct current (DC) electricity to alternating current (AC) electricity; receive an indication of an output current of the AC electricity; extract, from the virtual oscillator circuit, a virtual current based on the output current of the AC electricity; and receive an indication of an input voltage of the DC electricity, wherein the oscillating waveform is generated further based on the input voltage of the DC electricity.

16. The device of claim 15, wherein the processor is configured to: compute a linear combination of a virtual voltage across the virtual capacitor and a virtual current through the virtual inductor to produce an alternate virtual output voltage, and generate the oscillating waveform based on the alternate virtual output voltage.

17. The device of claim 15, wherein the processor is further configured to modify at least one of a virtual inductance of the virtual inductor or a scaling factor of the virtual cubic voltage-dependent current source based on a specified output power and at least one of: a scaled virtual output voltage of the virtual oscillator circuit or the output current of the AC electricity.

18. A method comprising: determining, for a processor operatively coupled to power electronics configured to convert direct current (DC) electricity to alternating current (AC) electricity, based on a specified open-circuit voltage value for the power electronics, a voltage scaling factor value, wherein the processor is configured to implement a virtual oscillator circuit comprising a virtual capacitor, a virtual inductor, a virtual negative-conductance element and a virtual cubic voltage-dependent current source; configuring the processor to scale, using the voltage scaling factor value, a virtual oscillator output voltage of the virtual oscillator circuit thereby producing a scaled virtual oscillator output voltage for use in controlling the power electronics; determining, for the processor, based on a specified rated power output value for the power electronics and a specified rated-power voltage value for the power electronics that corresponds to the specified rated power output value, a current scaling factor value; configuring the processor to generate a virtual oscillator input current by scaling, using the current scaling factor value, an output current of the power electronics; determining, for the processor, based on the specified open-circuit voltage value for the power electronics and the specified rated-power voltage value for the power electronics, a conductance value; configuring the processor to implement the virtual negative-conductance element based on the conductance value; determining, for the processor, based on the conductance value, a cubic coefficient value; configuring the processor to implement the virtual cubic voltage-dependent current source based on the cubic coefficient value; determining, for the processor, a capacitance value based on at least one of: (i) a specified rated reactive-power output value for the power electronics, (ii) a specified maximum frequency offset value, (iii) the specified rated power output value for the power electronics, (iv) the specified open-circuit voltage value for the power electronics, (v) the specified rated-power voltage value for the power electronics, (vi) a specified rise time value, (vii) a specified nominal system frequency value, or (viii) a specified ratio of an amplitude of a third harmonic to an amplitude of a fundamental; configuring the processor to implement the virtual capacitor based on the capacitance value; determining, for the processor, based on the capacitance value and the specified nominal system frequency value, an inductance value; and configuring the processor to implement the virtual inductor based on the inductance value.

19. The method of claim 18, wherein: the specified rated power output value for the power electronics comprises a specified rated real power output value for the power electronics, and the specified rated-power voltage value for the power electronics comprises a specified rated real-power voltage value for the power electronics.

20. The method of claim 18, wherein: the processor is further configured to implement a transformation unit that computes a linear combination of a voltage across the virtual capacitor and a current through the virtual inductor to produce an alternate virtual output voltage, the specified rated power output value for the power electronics comprises a specified rated reactive power output value for the power electronics, and the specified rated-power voltage value for the power electronics comprises a specified rated reactive-power voltage value for the power electronics.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0009] FIG. 1 is a conceptual diagram illustrating one example of a microgrid system (system 100) using VOC, in accordance with one or more aspects of the present disclosure.

[0010] FIG. 2 is a circuit diagram illustrating one example of a single-phase power-electronic inverter implemented using VOC, in accordance with one or more aspects of the present disclosure.

[0011] FIG. 3 is a flow diagram illustrating example operations for designing a VOC system, in accordance with one or more aspects of the present disclosure.

[0012] FIG. 4 is a graphical plot illustrating an example voltage-regulation curve resulting from example design specifications and voltage-regulation parameters determined therefrom, in accordance with one or more aspects of the present disclosure.

[0013] FIG. 5 is a graphical plot illustrating example synchronization times of droop control and virtual oscillator-based control, in accordance with one or more aspects of the present disclosure.

[0014] FIGS. 6A and 6B are graphical plots illustrating active- and reactive-power sharing among two inverters for virtual oscillator-based control and droop control, in accordance with one or more aspects of the present disclosure.

[0015] FIG. 7 is a circuit diagram illustrating another example of a single-phase power-electronic inverter implemented using VOC, in accordance with one or more aspects of the present disclosure.

[0016] FIG. 8 is a graphical plot illustrating VOC system dynamics on a polar plane, in accordance with one or more aspects of the present disclosure.

[0017] FIG. 9 is a flow diagram illustrating additional example operations for designing a VOC, system, in accordance with one or more aspects of the present disclosure.

[0018] FIG. 10 is a circuit diagram illustrating another example of a single-phase power-electronic inverter implemented using VOC, in accordance with one or more aspects of the present disclosure.

[0019] FIG. 11 is a circuit diagram illustrating another example of a single-phase power-electronic inverter implemented using VOC, in accordance with one or more aspects of the present disclosure.

[0020] FIG. 12 is a block diagram illustrating one example of a dynamic model of a power-tracking loop, in accordance with one or more aspects of the present disclosure.

[0021] FIGS. 13A and 13B are graphical plots showing real and reactive power references and measured values for a prototype VOC-based power-tracking inverter, in accordance with one or more aspects of the present disclosure.

DETAILED DESCRIPTION

[0022] The present disclosure provides systems, devices, and methods for rapidly stabilizing arbitrary initial conditions and load transients to a stable limit cycle using virtual oscillator control (VOC). For example, the techniques described herein may allow a plurality of oscillating signal generators, such as power-electronic inverters connected to an electrical network (e.g., a microgrid), to synchronize with one another without centralized management or communication among the signal generators. Instead, each signal generator may use local measurements and local control actions such that system-wide synchronization is obtained.

[0023] In another aspect, the present disclosure provides a formalized design approach for implementing VOC to meet required design specifications. That is, the present disclosure also provides for the design of the virtual oscillators that underpin the control techniques described herein. This may be referred to herein as the oscillator synthesis problem.

[0024] In yet another aspect, the present disclosure provides a suite of extensions to VOC, including the ability to trade off real power for frequency and/or trade off reactive power for voltage amplitude, as may be desirable in inductive networks. Additionally, the suite of extensions to VOC described herein also provides improvements to such formalized design approaches, which may pave the way for multi-mode (grid-connected and islanded) operation.

[0025] The techniques of the present disclosure are described herein with respect to power-electronic inverters, such as those operating in a microgrid. However, these techniques may also be applicable to various other areas such as standalone systems for remote power applications, uninterruptible AC power supplies, forward operating base power systems, or in any other systems that would benefit from fast, self-referential synchronization of multiple oscillating signals. Furthermore, Van der Pol oscillators are primarily discussed herein. However, the techniques of the present disclosure may be applied using any of a family of weakly nonlinear, limit-cycle oscillators, such as the Liénard oscillator and others. Hereinafter, inverters controlled in accordance with the present disclosure are termed virtual-oscillator controlled (VO-controlled) inverters.

[0026] FIG. 1 is a conceptual diagram illustrating one example of a microgrid system (system 100) using VOC, in accordance with one or more aspects of the present disclosure. FIG. 1 represents only one example of a system configured to use VOC, and various other systems and devices may be configured to use VOC in accordance with the techniques described herein.

[0027] Microgrids often include a collection of heterogeneous energy sources (e.g., photovoltaic arrays, fuel cells, and energy-storage devices). Microgrid 100, for instance, includes energy sources 102A-102C (collectively “energy sources 102”). Energy sources may be interfaced to an AC electric distribution network, such as distribution network 104.

[0028] Distribution network 104 is connected to bulk AC electrical system 106. For instance, bulk AC electrical system 106 may represent the “grid,” generally. In some examples, network 104 may be “islanded” and operated independently from the bulk AC system. This is shown in the example of FIG. 1 by the dotted line between distribution network 104 and bulk AC electrical system 106.

[0029] In microgrids and other electrical networks, energy conversion is often performed by power-electronic inverters. In the example of FIG. 1, for instance, energy from energy sources 102 is converted by inverters 110A-110C (collectively “inverters 110”), respectively.

[0030] In islanded settings, the control challenge is to regulate the amplitude and frequency of the inverters' terminal voltage such that high power quality can be guaranteed to the loads in the network. The related-art control strategy of droop control is relatively ubiquitous in this domain. Droop control linearly trades off the inverter-voltage amplitude and frequency with real- and reactive-power output.

[0031] The techniques of the present disclosure, however, depart from droop control to pioneer time-domain control methods. In particular, the present disclosure provides a nonlinear control strategy where inverters 110 are controlled to emulate the dynamics of weakly nonlinear limit-cycle oscillators. This strategy is termed “virtual oscillator control” or “VOC” herein, since such nonlinear oscillators may be virtual (e.g., programmed on a digital controller) instead of physical.

[0032] VOC presents appealing circuit- (e.g., inverter) and system- (e.g., microgrid electrical network) level advantages. From a system-level perspective, synchronization emerges in connected electrical networks of inverters with VOC without any communication, and primary-level voltage- and frequency-regulation objectives are ensured in a decentralized fashion. At the circuit level, each of inverters 110 is able to rapidly stabilize arbitrary initial conditions and load transients to a stable limit cycle. As such, VOC is fundamentally different than droop control. While VOC acts on instantaneous time-domain signals, droop control is based on phasor electrical quantities and the notion of an electrical frequency that are only well defined on slow AC-cycle time scales. As further described herein, however, the sinusoidal state behavior of VOC can be engineered to correspond to droop laws as well.

[0033] AC performance requirements for inverters 110 (e.g., voltage and frequency regulation, harmonics, dynamic response, etc.) may typically be specified with the aid of phasor quantities that are only valid in the quasi-stationary sinusoidal steady state. As such, given the intractability of obtaining closed-form solutions to the oscillator dynamic trajectories, from the outset it may be unclear how to design the nonlinear oscillators such that the controlled inverters meet prescribed specifications. The present disclosure provides one solution to the oscillator synthesis problem by using averaging- and perturbation-based nonlinear-systems analysis methods as detailed below. This process focuses on an averaged dynamical model for the nonlinear oscillators that couples the real- and reactive-power outputs to the terminal-voltage dynamics of an inverter. Analyzing this averaged model in the sinusoidal steady state uncovers the voltage- and frequency-regulation characteristics of virtual-oscillator-controlled inverters. In addition, the present disclosure provides perturbation-based methods that, in general, may allow for approximating solutions to periodic nonlinear dynamical systems when analytical closed-form solutions cannot be found. This analysis parameterizes the higher-order harmonic content in the inverter output as a function of the oscillator parameters, further aiding in oscillator design.

[0034] FIG. 2 is a circuit diagram illustrating one example of a single-phase power-electronic inverter (e.g., inverter 200) implemented using VOC, in accordance with one or more aspects of the present disclosure. Specifically, FIG. 2 illustrates a Van der Pol VOC implementation for a single phase inverter, inverter 200 represents only one example of an inverter implemented using VOC, and various other structures and configurations may use one or more of the techniques described herein. For example, other devices may include additional or alternative components not shown in the example of FIG. 2.

[0035] In the example of FIG. 2, inverter 200 includes power electronics switches 202. Switches 202 implement a single-phase H-bridge. As such, alternative sets of switches 202 may be opened and closed in sequence to convert a DC signal (shown as .sub.dc in FIG. 2) into an AC signal. As shown in FIG. 2, the opening and closing of power electronics switches 202 is controlled by signals received from control unit 210, described below.

[0036] Inverter 200, in the example of FIG. 2, also includes LCL output filter 204. LCL output filter 204 may reduce high-order harmonics in the inverter terminal voltage. LCL output filter 204 includes inverter-side inductors (shown as L.sub.f in FIG. 2), an AC filter capacitor (shown as C.sub.f in FIG. 2), and grid-side inductors (shown as L.sub.g in FIG. 2).

[0037] In the example of FIG. 2, inverter 200 also includes control unit 210. Control unit 210 may be a processor, such as a microcontroller, a digital signal processor (DSP), a field-programmable gate array (FPGA) or other suitable device configured to implement the various modules described herein. Control unit 210 implements a discrete realization of the nonlinear dynamics of a Van der Pol oscillator (e.g., virtual oscillator 212).

[0038] Virtual oscillator 212, in the example of FIG. 2, is composed of the parallel connection of virtual capacitor 220, virtual inductor 222, virtual negative-conductance element 224 and virtual cubic voltage-dependent current source 226. Virtual capacitor 220 and virtual inductor 222 represent a harmonic oscillator with inductance, L, and capacitance, C, yielding a resonant frequency, ω*=1/√{square root over (LC)}). Virtual negative-conductance element 224 may have a conductance value of −σ. A voltage across virtual capacitor 220 is denoted by .sub.C and a current through virtual inductor 222 is denoted by i.sub.L. A current consumed by virtual cubic voltage-dependent current source 226 is given by α.sub.C.sup.3, where α is a positive constant.

[0039] In the example of FIG. 2, virtual oscillator 212 is coupled to physical electrical signals in inverter 200 through current-scaling unit 213 and voltage-scaling unit 214 (using scaling factor values κ.sub.i and κ.sub.v, respectively). Specifically, an output current of inverter 200, i, may be processed by analog-to-digital converter (ADC) 215 and provided to current-scaling unit 213. Current-scaling unit 213 may multiply the received signal by κ.sub.i. The resulting scaled current may then be extracted from virtual oscillator circuit 212. The resulting value of .sub.C (the voltage across virtual capacitor 220) is provided to voltage-scaling unit 214. Voltage-scaling unit 214 may multiply the received signal by κ.sub.V to produce the signal .

[0040] Virtual oscillator 212, in the example of FIG. 2, also includes ADC 216 and modulation unit 217. ADC 216 converts the DC bus voltage V.sub.dc to a digital value, which is provided to modulation unit 217. Modulation unit 217 receives the digitized DC bus voltage and the signal from voltage-scaling unit 214 and divides the latter by the former. That is, modulation unit 217 divides by V.sub.dc to produce a modulation signal m.

[0041] For the single-phase inverter topology depicted in FIG. 2, the PWM signal, m(t), may be constructed as follows

[00001]$\begin{array}{cc}m\left(t\right)=\frac{v\left(t\right)}{v_{\mathrm{dc}}},& \left(1\right)\end{array}$

where .sub.dc is the DC-bus voltage.

[0042] In the example of FIG. 2, virtual oscillator 212 includes pulse width modulation (PWM) unit 218. PWM receives the modulation signal m and compares its instantaneous value to a carrier waveform to generate the control signals for power-electronics switches 202. If the modulation signal value is higher than the carrier waveform value, one set of power-electronics switches 202 will be “on” and the others will be “off.” If the modulation falls below the carrier, another set of switches turns “on” and the first turns “off.” In the example of FIG. 2, the switching frequency may be equal to the carrier frequency.

[0043] The carrier waveform, in the example of FIG. 2, is a free-running triangular waveform. In practice, the carrier waveform may be generated by a digital counter (not shown) that counts up and/or down between a minimum and maximum value to create the triangular waveform, or by a suitable analog signal generator. The modulation signal (which, for example, may look like a sinusoid in steady-state), is typically scaled such that its max and min lie within the max and min of the carrier waveform. Since the carrier has a much higher frequency (e.g., tens of kHz) than the modulation signal (e.g., 60 Hz), there will be many transition points where the carrier will rise above and fall below the modulation signal.

[0044] The switching period, T.sub.sw, is much smaller than the period of the modulation signal, which in this setting is approximately 2π/ω*, where ω* is the resonant frequency of the LC harmonic oscillator. Consequently, the switch-cycle average of the instantaneous inverter-terminal voltage—denoted by in FIG. 2—is approximately equal to the scaled virtual capacitor voltage, :

[00002]$\begin{array}{cc}\frac{1}{T_{\mathrm{sw}}}.Math.{\int }_{s=t-T_{\mathrm{sw}}}^t.Math.v\left(s\right).Math.\mathrm{ds}=m\left(t\right).Math.v_{\mathrm{dc}}=v\left(t\right).& \left(2\right)\end{array}$

It is assumed hereafter that the switch-cycle-average inverter terminal voltage is equal to and this is referred to as the inverter terminal voltage.

[0045] The following is an analysis of the terminal-voltage dynamics of the VO-controlled inverter 200 as shown in FIG. 2. The resulting dynamic model may then be leveraged to synthesize the underlying virtual oscillator such that inverter 200 meets AC performance specifications. The description herein is focused on real- and reactive-power ratings, voltage limits, output-current ratings, maximum frequency deviation, dynamic response, and higher-order harmonics. These performance specifications are connected below to the inherently nonlinear dynamics of the Van der Pol oscillator. In other examples, additional and/or other performance specifications may also be utilized.

[0046] To begin, the voltage- and frequency-regulation characteristics of a VO-controlled inverter (e.g., inverter 200 of FIG. 2) may be derived. The derivation is based on an averaging analysis of the Van der Pol oscillator dynamics.

[0047] The dynamics of the virtual-oscillator inductor current, i.sub.L, and inverter terminal voltage, , are given by

[00003]$\begin{array}{cc}L.Math.\frac{{\mathrm{di}}_L}{\mathrm{dt}}=\frac{v}{{\kappa }_v},.Math.C.Math.\frac{\mathrm{dv}}{\mathrm{dt}}=-\alpha .Math.\frac{v^3}{{\kappa }_v^2}+\sigma .Math..Math.v-{\kappa }_v.Math.i_L-{\kappa }_v.Math.{\kappa }_i.Math.i.& \left(3\right)\end{array}$

The inverter terminal voltage is parameterized in one of the two forms below

(t)=√{square root over (2)}V(t)cos(ωt+θ(t))=√{square root over (2)}V(t)cos(φ(t)),  (4)

where ω is the electrical frequency, θ(t) represents the phase offset with respect to ω, and φ(t) is the instantaneous phase angle. In order to completely specify the terminal voltage at any instant, it is necessary to obtain the dynamical equations that govern the evolution of the RMS-voltage, V(t), and phase offset, θ(t), (or equivalently, the instantaneous phase angle, φ(t)). To this end, the following definitions will help simplify notation:

[00004]$\begin{array}{cc}\mathrm{.Math.}:=\sqrt{\frac{L}{C}},g\left(y\right):=y-\frac{\beta }{3}.Math.y^3,\beta :=\frac{3.Math.\alpha }{{\kappa }_v^2.Math.\sigma }.& \left(5\right)\end{array}$

A state-space model of the VO-controlled inverter in Cartesian coordinates aimed at recovering the RMS-voltage amplitude and instantaneous phase dynamics can be formulated with a scaled version of the inductor current and the inverter terminal voltage selected as states, x :=κ.sub.vεi.sub.L, and y :=. With the aid of the g(•) and ε defined in (5), (3) may be rewritten in the time coordinates τ=ω*t=(1/√{square root over (LC)})t as follows:

[00005]$\begin{array}{cc}\stackrel{.}{x}=\frac{\mathrm{dx}}{d.Math..Math.\tau }=y,& \left(6\right)\\ \stackrel{.}{y}=\frac{\mathrm{dy}}{d.Math..Math.\tau }=-x+\mathrm{.Math.\sigma }.Math..Math.g\left(y\right)-{\mathrm{.Math.\kappa }}_v.Math.{\kappa }_i.Math.i.& \left(7\right)\end{array}$

Next, with the coordinate transformation

[00006]$\begin{array}{cc}\sqrt{2}.Math.V=\sqrt{x^2+y^2},\phi ={\tan }^{-1}\left(\frac{x}{y}\right),& \left(8\right)\end{array}$

applied to (6)-(7), the following dynamical model for the RMS terminal-voltage amplitude, V, and the instantaneous phase angle, φ, are recovered:

[00007]$\begin{array}{cc}\stackrel{.}{V}=\frac{d.Math.V}{d.Math..Math.\tau }=\frac{\mathrm{.Math.}}{\sqrt{2}}.Math.\left(\sigma .Math..Math.g\left(\sqrt{2}.Math.V.Math..Math.\cos \left(\phi \right)\right)-{\kappa }_v.Math.{\kappa }_i.Math.i\right).Math.\cos \left(\phi \right),.Math.\stackrel{.}{\phi }=\frac{d.Math..Math.\phi }{d.Math..Math.\tau }=1-\frac{\mathrm{.Math.}}{\sqrt{2}.Math.V}.Math.\left(\sigma .Math..Math.g\left(\sqrt{2}.Math.V.Math..Math.\cos \left(\phi \right)\right)-{\kappa }_v.Math.{\kappa }_i.Math.i\right).Math.\sin \left(\phi \right).& \left(9\right)\end{array}$

As ε0, the so-called quasi-harmonic limit is transitioned, where the (unloaded) oscillator exhibits near-sinusoidal oscillations at the resonant frequency of the LC harmonic oscillator:

[00008]$\begin{array}{cc}{\omega }^{*}=\frac{1}{\sqrt{\mathrm{LC}}}.& \left(10\right)\end{array}$

[0048] A demonstration of ε as a design parameter that has bearing on the dynamic response and the harmonics of the system can be found below with respect to design procedure. In particular, a small value of ε may ensure near-sinusoidal oscillations, potentially at the expense of a sluggish dynamic response.

[0049] By focusing on the parametric regime characterized by ε0, notations of periodic averaging may be leveraged to further simplify and analyze the weakly nonlinear periodic dynamics in (9). As a primer, consider a time-varying dynamical system {dot over (x)}=εf(x,τ,ε), with time-periodic vector field

f(x,τ,ε)=f(x,τ+T,ε),  (11)

with period T>0 and a small parameter ε>0. The associated time-averaged dynamical system may be defined as

[00009]$\begin{array}{cc}\stackrel{.}{\stackrel{\_}{x}}=\mathrm{.Math.}.Math..Math.f_{\mathrm{avg}}\left(\stackrel{\_}{x}\right)=\frac{\mathrm{.Math.}}{T}.Math.{\int }_{\tau =0}^T.Math.f\left(\stackrel{\_}{x},\tau ,0\right).Math.d.Math..Math.\tau .& \left(12\right)\end{array}$

The averaged system, f.sub.avg(x), is autonomous and in general, more amenable to analysis compared to the original non-averaged system f(x,τ,ε). Furthermore, for ε0, the accuracy is not compromised by analyzing the averaged system, since it follows that the difference in the state variables corresponding to the original and averaged systems is of (ε). In particular,

x(τ,ε)−x(ετ)=(ε).  (13)

Subsequently, the time average of a periodic signal x(t) with period T>0 is denoted by x and defined as follows:

[00010]$\begin{array}{cc}\stackrel{\_}{x}:=\frac{1}{T}.Math.{\int }_{s=0}^T.Math.x\left(s\right).Math.\mathrm{ds}.& \left(14\right)\end{array}$

[0050] With these preliminaries in place, the dynamics in (9) averaged over one AC cycle, 2π/ω*, under the implicit assumption that ε=√{square root over (L/C)}0 are given by the following set of coupled nonlinear differential equations:

[00011]$\begin{array}{cc}\frac{d}{\mathrm{dt}}.Math.\stackrel{\_}{V}=\frac{\sigma }{2.Math..Math.C}.Math.\left(\stackrel{\_}{V}-\frac{\beta }{2}.Math.{\stackrel{\_}{V}}^3\right)-\frac{{\kappa }_v.Math.{\kappa }_i}{2.Math..Math.C.Math.\stackrel{\_}{V}}.Math.\stackrel{\_}{P},& \left(15\right)\\ \frac{d}{\mathrm{dt}}.Math.\stackrel{\_}{\theta }={\omega }^{*}-\omega +\frac{{\kappa }_v.Math.{\kappa }_i}{2.Math..Math.C.Math.{\stackrel{\_}{V}}^2}.Math.\stackrel{\_}{Q},& \left(16\right)\end{array}$

where P and Q are the average real- and reactive-power outputs of the VO-controlled inverter (measured with respect to the nominal frequency, ω*), respectively.

[0051] Using this nonlinear control strategy, close to the sinusoidal steady state (recovered in the quasi-harmonic limit ε0), the voltage-amplitude and phase dynamics are directly linked to the average real- and reactive-power outputs of the inverter, respectively. Consequently, these averaged dynamics are leveraged herein for designing virtual oscillator controlled inverters that satisfy voltage- and frequency-regulation specifications in sinusoidal steady state. To this end, the equilibrium solutions corresponding to (15)-(16) are useful, as they establish the voltage- and frequency-regulation characteristics of the VO-controlled inverter.

[0052] For the voltage-regulation characteristic, the equilibria of (15) can be recovered from the solutions of the nonlinear equation:

[00012]$\begin{array}{cc}0=\frac{\sigma }{2.Math..Math.C}.Math.\left({\stackrel{\_}{V}}_{\mathrm{eq}}-\frac{\beta }{2}.Math.{\stackrel{\_}{V}}_{\mathrm{eq}}^3\right)-\frac{{\kappa }_v.Math.{\kappa }_i}{2.Math..Math.C.Math.{\stackrel{\_}{V}}_{\mathrm{eq}}}.Math.{\stackrel{\_}{P}}_{\mathrm{eq}},& \left(17\right)\end{array}$

where V.sub.eq and P.sub.eq represent the equilibrium steady-state RN/IS-voltage amplitude and average real-power output, respectively. Rearranging terms in (17) provides the following power-balance condition for the VO-controlled inverter:

[00013]$\begin{array}{cc}\frac{\mathrm{\sigma \beta }}{2}.Math.{\stackrel{\_}{V}}_{\mathrm{eq}}^4-\sigma .Math.{\stackrel{\_}{V}}_{\mathrm{eq}}^2+{\kappa }_v.Math.{\kappa }_i.Math.{\stackrel{\_}{P}}_{\mathrm{eq}}=0.& \left(18\right)\end{array}$

The positive roots of (18) are given by

[00014]$\begin{array}{cc}{\stackrel{\_}{V}}_{\mathrm{eq}}={{\kappa }_v\left(\frac{\sigma \pm \sqrt{{\sigma }^2-6.Math.\alpha \left({\kappa }_i/{\kappa }_v\right).Math.{\stackrel{\_}{P}}_{\mathrm{eq}}}}{3.Math.\alpha }\right)}^{\frac{1}{2}},& \left(19\right)\end{array}$

where the fact that σβ=3α/κ.sub.v.sup.2 was used (see (5)). Notice that (19) has two roots. Both roots are real valued if the equilibrium real-power output satisfies

[00015]$\begin{array}{cc}0<{\stackrel{\_}{P}}_{\mathrm{eq}}<{\stackrel{\_}{P}}_{\mathrm{cr}}:=\frac{{\sigma }^2}{6.Math.\alpha \left({\kappa }_i/{\kappa }_v\right)},& \left(20\right)\end{array}$

where P.sub.cr is referred to as the critical value for real power. The corresponding critical value of the terminal voltage is given by

[00016]$\begin{array}{cc}{\stackrel{\_}{V}}_{\mathrm{cr}}:={\kappa }_v.Math.\sqrt{\frac{\sigma }{3.Math.\alpha }}.& \left(21\right)\end{array}$

[0053] Under a set of mild power-flow decoupling approximations, the high-voltage solution in (19) is locally asymptotically stable. Hereinafter, the high-voltage solution of (19) is referred to with a slight abuse of notation as V.sub.eq. It is worth mentioning that the high-voltage root is a decreasing function of P.sub.eq over the range 0≦P.sub.eq≦P.sub.cr. That is, steady-state voltage “droops” with increasing real power output. Finally, also note that the open-circuit voltage of the VO-controlled inverter, V.sub.oc, can be Obtained by substituting P.sub.eq=0 into the high-voltage root in (19):

[00017]$\begin{array}{cc}{\stackrel{\_}{V}}_{\mathrm{oc}}={\kappa }_v.Math.\sqrt{\frac{2.Math.\sigma }{3.Math.\alpha }}.& \left(22\right)\end{array}$

[0054] For the frequency-regulation characteristic, consider the phasor-angle dynamics in (16). The equilibrium of (16) returns the frequency of the VO-controlled inverter:

[00018]$\begin{array}{cc}{\omega }_{\mathrm{eq}}={\omega }^{*}+\frac{{\kappa }_v.Math.{\kappa }_i}{2.Math..Math.C.Math.{\stackrel{\_}{V}}_{\mathrm{eq}}^2}.Math.{\stackrel{\_}{Q}}_{\mathrm{eq}},& \left(23\right)\end{array}$

where V.sub.eq is the stable high-voltage equilibrium obtained from (19), and Q.sub.eq is the average reactive-power output of the VO-controlled inverter.

[0055] The next focus is the dynamic response of the VO-controlled inverter (e.g., inverter 200 shown in FIG. 2). Specifically of interest is how quickly the terminal voltage of an unloaded inverter builds up to the open-circuit voltage. The voltage dynamics of interest are recovered from (15) by setting P=0:

[00019]$\begin{array}{cc}\frac{d}{\mathrm{dt}}.Math.\stackrel{\_}{V}=\frac{\sigma }{2.Math..Math.C}.Math.\left(V-\beta .Math.\frac{1}{2}.Math.{\stackrel{\_}{V}}^3\right).& \left(24\right)\end{array}$

Since (24) is a variable-separable ordinary differential equation, both sides can be integrated, setting the limits from 0.1V.sub.oc to 0.9V.sub.oc (without loss of generality). The time taken for this excursion is defined as the rise time. It is denoted by t.sub.rise, and given by the solution of:

[00020]$t_{\mathrm{rise}}={\frac{2}{\mathrm{.Math.\sigma }}\left[\log .Math..Math.\stackrel{\_}{V}-\frac{1}{2}.Math.\log .Math..Math.1-\frac{\beta }{2}.Math.{\stackrel{\_}{V}}^2.Math.\right]}_{0.1.Math..Math.{\stackrel{\_}{V}}_{\mathrm{oc}}}^{0.9.Math..Math.{\stackrel{\_}{V}}_{\mathrm{oc}}}.$

Evaluating the limits above results in the recovery of

[00021]$\begin{array}{cc}{t}_{\mathrm{rise}}\approx \frac{6}{\omega *\mathrm{.Math.\sigma }}.& \left(25\right)\end{array}$

The approximation in (25) indicates that the rise time, t.sub.rise, is inversely proportional to E. This aspect will be leveraged in the design procedure below, specifically with regard to the virtual oscillator capacitance, C.

[0056] The above analysis pertains to an unloaded inverter. For an inverter loaded to its rated real power rating with a resistive load, R.sub.rated, the rise time is given by

[00022]$t_{\mathrm{rise}}\approx \frac{6}{\omega *{\mathrm{.Math.\sigma }}^{\prime }},$

where

[00023]${\sigma }^{\prime }=\left(\sigma -\frac{{\kappa }_v.Math.{\kappa }_i}{R_{\mathrm{rated}}}\right).$

The analysis on inverter design may additionally or alternatively be performed with this specification of rise time, if need be.

[0057] The following portion derives a closed-form analytical expression for the amplitude of the third harmonic of an unloaded VO-controlled inverter (e.g., inverter 200 of FIG. 2). In particular, the effect of the nonlinear forcing term in (3), as executed within control unit 204, on the low-order harmonic content of the AC output is investigated. This analysis provides parameter selection with the aim of bounding the ratio of the amplitude of the third harmonic to the fundamental. To this end, perturbation methods and the method of multiple scales are used. These methods seek approximate analytical solutions to nonlinear dynamical systems where exact solutions cannot be found.

[0058] Consider the non-averaged dynamics of the terminal-voltage magnitude in an unloaded VO-controlled inverter:

−εσ(1−β.sup.2)+=0,  (26)

where υ and β are defined in (5), and as before, while operating in the quasi-harmonic limit ε0. This model follows from expressing (6) and (7) as a second-order system with the input current i=0. The objective is an approximate solution to (26) that can be expressed as:

(τ,ε)≈.sub.0(τ,{tilde over (τ)})+ε.sub.1(τ,{tilde over (τ)}).  (27)

The solution is written with respect to two time scales: the original time scale τ, and a slower time scale, {tilde over (τ)}:=ετ. While higher-order time scales (e.g., ε.sup.2τ, ε.sup.3τ, etc.) can be analyzed in a similar fashion to obtain approximate solutions correct to higher-order terms, the analysis performed herein is: i) valid up to (ε), ii) yields an approximate amplitude for the third harmonic, and iii) provides error terms of (ε.sup.2). Substituting (27) in (26), and retaining only (ε) terms results in:

[00024]$\begin{array}{cc}\left(\frac{{\partial }^2.Math.v_0}{\partial {\tau }^2}+v_0\right)+\mathrm{.Math.}\left(\frac{{\partial }^2.Math.v_0}{\partial {\tau }^2}+v_1+2.Math.\frac{{\partial }^2.Math.v_0}{\partial \tau .Math.\partial \tilde{\tau }}-\sigma \left(1-\beta .Math..Math.v_0^2\right).Math.\frac{\partial v_0}{\partial \tau }\right)=0.& \left(28\right)\end{array}$

Note that (28) must hold for any small parameter ε. This can be ensured if:

[00025]$\begin{array}{cc}\frac{{\partial }^2.Math.v_0}{\partial {\tau }^2}+v_0=0,& \left(29\right)\\ \frac{{\partial }^2.Math.v_1}{\partial {\tau }^2}+v_1+2.Math.\frac{{\partial }^2.Math.v_0}{\partial \tau .Math.\partial \tilde{\tau }}-\sigma \left(1-\beta .Math..Math.v_0^2\right).Math.\frac{\partial v_0}{\partial \tau }=0.& \left(30\right)\end{array}$

[0059] Since (29) represents the dynamics of a simple harmonic oscillator, the corresponding closed-form solution can be expressed as:

.sub.0(τ,{tilde over (τ)}):=a.sub.0({tilde over (τ)})cos(τ+ρ.sub.0({tilde over (τ)})),  (31)

where a.sub.0({tilde over (τ)}) and ρ.sub.0({tilde over (τ)}) are amplitude and phase terms that vary in the slow time scale specified by {tilde over (τ)}. For notational convenience, define the orthogonal signal .sub.0.sup.⊥(τ,{tilde over (τ)}) associated with .sub.0(τ,{tilde over (τ)}) in (31) as follows:

.sub.0.sup.⊥(τ,{tilde over (τ)}):=a.sub.0({tilde over (τ)})sin(τ+ρ.sub.0({tilde over (τ)})).  (32)

Substituting for .sub.0 from (31) into (30):

[00026]$\begin{array}{cc}\begin{array}{c}\frac{\partial v_1}{\partial {\tau }^2}+v_1=.Math.-2.Math.\frac{{\partial }^2.Math.v_0}{\partial \tau .Math.\partial \tilde{\tau }}+\sigma \left(1-\beta .Math..Math.v_0^2\right).Math.\frac{\partial v_0}{\partial \tau }\\ =.Math.2.Math.\frac{\partial v_0^{\perp }}{\partial \tilde{\tau }}-\sigma .Math..Math.v_0^{\perp }+\mathrm{\sigma \beta }.Math..Math.v_0^2.Math.v_0^{\perp }\\ =.Math.2.Math.\frac{\partial a_0}{\partial \tilde{\tau }}.Math.\sin \left(\tau +{\rho }_0\right)+2.Math..Math.v_0.Math.\frac{\partial {\rho }_0}{\partial \tilde{\tau }}-\\ .Math.\sigma .Math..Math.v_0^{\perp }+\mathrm{\sigma \beta }.Math..Math.a_0^3\left(\sin \left(\tau +{\rho }_0\right)-{\sin }^3\left(\tau +{\rho }_0\right)\right)\\ =.Math.2.Math.\frac{\partial a_0}{\partial \tilde{\tau }}.Math.\sin \left(\tau +{\rho }_0\right)+2.Math..Math.v_0.Math.\frac{\partial {\rho }_0}{\partial \tilde{\tau }}-\\ .Math.\sigma .Math..Math.v_0^{\perp }+\frac{\mathrm{\sigma \beta }.Math.}{4}.Math.a_0^3\left(\sin \left(3.Math.\tau +3.Math.{\rho }_0\right)+\sin \left(\tau +{\rho }_0\right)\right),\end{array}& \left(33\right)\end{array}$

where, in the last line of (33), the trigonometric identity sin 3θ=3 sin θ−4 sin.sup.3 θ was used. Grouping together the coefficients that multiply the sin(τ+ρ.sub.0) and cos(τ+ρ.sub.0) terms, the last line of (33) can be rewritten as follows:

[00027]$\begin{array}{cc}\frac{{\partial }^2.Math.v_1}{\partial {\tau }^2}+v_1=\left(2.Math.\frac{\partial a_0}{\partial \tilde{\tau }}-\sigma .Math..Math.a_0+\frac{\mathrm{\sigma \beta }}{4}.Math.a_0^3\right).Math.\sin \left(\tau +{\rho }_0\right)+\left(2.Math..Math.a_0.Math.\frac{\partial {\rho }_0}{\partial \tilde{\tau }}\right).Math.\cos \left(\tau +{\rho }_0\right)+\frac{\mathrm{\sigma \beta }}{4}.Math.a_0^3.Math.\sin \left(3.Math.\tau +3.Math.{\rho }_0\right).& \left(34\right)\end{array}$

[0060] The coefficients that multiply the sin(τ+ρ.sub.0) and cos(τ+ρ.sub.0) terms have to be forced to zero to ensure that unbounded terms of the form τ sin(τ+ρ.sub.0) and τ cos(τ+ρ.sub.0) do not appear in the solution for .sub.1. Consequently, the following can be recovered:

[00028]$\begin{array}{cc}2.Math.\frac{\partial a_0\left(\tilde{\tau }\right)}{\partial \tilde{\tau }}-\sigma .Math..Math.a_0\left(\tilde{\tau }\right)+\frac{\mathrm{\sigma \beta }}{4}.Math.{a_0\left(\tilde{\tau }\right)}^3=0,& \left(35\right)\\ 2.Math..Math.a_0\left(\tilde{\tau }\right).Math.\frac{\partial {\rho }_0\left(\tilde{\tau }\right)}{\partial \tilde{\tau }}=0.& \left(36\right)\end{array}$

Solving (36) with initial condition a.sub.0 (0) results in the following expression for a.sub.0({tilde over (τ)}):

[00029]$\begin{array}{cc}{a}_0\left(\tilde{\tau }\right)={\left(\frac{\beta }{4}+e^{-\eta -\sigma .Math.\tilde{\tau }}\right)}^{-\frac{1}{2}},e^{\eta }=\frac{a_0^2\left(0\right)}{1-\frac{\beta }{4}.Math.a_0^2\left(0\right)}.& \left(37\right)\end{array}$

It follows that the peak amplitude of the first harmonic in sinusoidal steady state is given by:

[00030]$\begin{array}{cc}{\lim }_{\tilde{\tau }.fwdarw.\infty }.Math.a_0\left(\tilde{\tau }\right)=:a_0=\frac{2}{\sqrt{\beta }}.& \left(38\right)\end{array}$

Note that the RMS value corresponding to the peak amplitude in (38) matches the expression for the open-circuit voltage in (22). From (37), it can also be inferred that a.sub.0({tilde over (τ)})≠0 if a.sub.0 (0)≠0. Therefore, it can be seen from (35) that ρ.sub.0({tilde over (τ)})=ρ.sub.0. That is, ρ.sub.0 is independent of {tilde over (τ)}.

[0061] With these observations in place, the following equation that governs the evolution of .sub.1(τ,{tilde over (τ)}) can be recovered from (34)

[00031]$\begin{array}{cc}\frac{{\partial }^2.Math.v_1}{\partial {\tau }^2}+v_1=\frac{2.Math.\sigma }{\sqrt{\beta }}.Math.\sin \left(3.Math.\tau +3.Math.{\rho }_0\right).& \left(39\right)\end{array}$

The particular solution to (39) is given by the general form:

[00032]$\begin{array}{cc}{v}_1\left(\tau \right)=-\frac{\sigma }{4.Math.\sqrt{\beta }}.Math.\sin \left(3.Math.\tau +3.Math.{\rho }_0\right).& \left(40\right)\end{array}$

From (31), (38), (40), and (27), it can be seen that the ratio of the amplitude of the third harmonic to the fundamental, a quantity denoted herein by δ.sub.3:1, is given by

[00033]$\begin{array}{cc}{\delta }_{3.Math.\text{:}.Math.1}=\frac{\mathrm{.Math.\sigma }}{8}.& \left(41\right)\end{array}$

If initial conditions for .sub.1 are taken into account while solving (39), (41) is correct up to (ε). Moreover, the expression in (41) indicates that the undesirable third-order harmonic is directly proportional to ε. This aspect is also leveraged in the design procedure below with respect to the virtual oscillator capacitance, C.

[0062] FIG. 3 is a flow diagram illustrating example operations for designing a VOC system, in accordance with one or more aspects of the present disclosure. More generally, FIG. 3 describes a systematic design strategy to synthesize nonlinear oscillators for inverter control. The design strategy shown in the example of FIG. 3 may aid in sculpting a desirable sinusoidal limit cycle that meets performance requirements specified in the lexicon of AC power systems. The disclosed strategy can be conveniently applied to design controllers for inverters with different power, voltage, and current ratings. In fact, for a system of parallel-connected inverters, the disclosed design strategy natively ensures power sharing in proportion to the inverter power ratings.

[0063] Tangentially related to the present disclosure are related-art efforts in the design of droop-controlled inverters. Particularly, there exists a wide body of literature that attempts to identify droop-control parameters that ensure inverters satisfy steady-state performance metrics and constraints such as voltage regulation, maximum frequency deviation, and proportional power sharing. Previous efforts in realizing VOC relied on an iterative design procedure involving simulation-based, open-circuit and full-rated-load tests to design the virtual oscillators. However, such ad-hoc design methods rely on repeated time-domain simulations to tune parameters and have not been affirmed by a rigorous nonlinear-systems analysis approach. The unified and formal design methodology for VOC described herein may particularly benefit practicing engineers and others interested in specific implementation aspects.

[0064] The operations shown in FIG. 3 are described with reference to inverter 200 as shown in FIG. 2. FIG. 3 represents only one example method for determining Van der Pol oscillator parameter values such that a VO-controlled inverter (e.g., inverter 200 of FIG. 2) satisfies a set of AC performance specifications. In other examples, the design procedure may include additional or different operations. The virtual oscillator parameter values to be determined in the example of FIG. 3 are summarized in Table I below. The parameters include i) scaling factors κ.sub.v and κ.sub.i, ii) voltage-regulation parameters σ and α, and iii) harmonic-oscillator parameters L and C.

TABLE-US-00001 TABLE I VOC PARAMETERS Symbol Description Units κ.sub.v Voltage scaling factor V/V κ.sub.i Current scaling factor A/A σ Conductance Ω.sup.−1 α Cubic coefficient A/V.sup.3 C Harmonic-oscillator capacitance F L Harmonic-oscillator inductance H

[0065] The performance specifications which the parameters are designed to satisfy include the open-circuit voltage V.sub.oc; rated real-power output and corresponding voltage P.sub.rated and V.sub.min, respectively; rated reactive-power output |Q.sub.rated|; maximum-permissible frequency deviation, rise time, and ratio of the amplitude of the third harmonic to the fundamental |Δω|.sub.max, t.sub.rise.sup.max, and δ.sub.3:1.sup.max, respectively.

[0066] Solely for ease of understanding, a running example implementing the operations of FIG. 3 is included herein using the set of AC performance specifications shown in Table II below. These example specifications result in a voltage regulation of ±5% around a nominal voltage of 120 V, and frequency-regulation of ±0.5 Hz around a nominal frequency of 60 Hz. A hardware prototype implementing this example design is described below.

TABLE-US-00002 TABLE II EXAMPLE AC PERFORMANCE SPECIFICATIONS Symbol Description Value Units V.sub.oc Open-circuit voltage 126 V (RMS) P.sub.rated Rated real power 750 W V.sub.min Voltage at rated power 114 V (RMS) |Q.sub.rated| Rated reactive power 750 VARs ω* Nominal system frequency 2π60  rad/s |Δω|.sub.max Maximum frequency offset 2π0.5 rad/s t.sub.rise.sup.max Rise time (25)    0.2 s δ.sub.3:1.sup.max Ratio of third-to-first harmonic (41)  2 %

[0067] Notice from FIG. 2 that the parameters κ.sub.v and κ.sub.i, respectively, determine the voltage and current scaling between the physical inverter terminal voltage and output current, and those of virtual oscillator 212. To standardize design, κ.sub.v is chosen herein such that when the voltage of virtual capacitor 220 is 1 V RMS, the inverter-terminal voltage is equal to the open-circuit voltage V.sub.oc. Furthermore, κ.sub.i is chosen herein such that when the output of virtual oscillator 212 is 1 A, inverter 200 is loaded to full rated capacity P.sub.rated.

[0068] As a result, in the example of FIG. 3, a value of the voltage scaling factor, κ.sub.v, may be determined based on the open-circuit voltage of the power electronics (operation 300). Specifically,

κ.sub.v:=V.sub.oc.  (42a)

A processor implementing VOC as described herein may be configured to scale an output voltage of the virtual oscillator using the voltage scaling factor value (operation 302). For instance, voltage-scaling unit 214 of control unit 210 may be configured to scale V.sub.C using the voltage scaling factor value to generate as shown in FIG. 2.

[0069] In the example of FIG. 3, a value of the current scaling factor, κ.sub.i, may be determined based on the rated real-power output and corresponding output voltage of the power electronics (operation 304). Specifically,

[00034]$\begin{array}{cc}{\kappa }_i:=\frac{{\stackrel{\_}{V}}_{\min }}{{\stackrel{\_}{P}}_{\mathrm{rated}}}.& \left(42.Math.b\right)\end{array}$

The processor implementing VOC as described herein may be configured to generate a virtual oscillator input current by scaling the power electronics output current using the current scaling factor value (operation 306). For instance, current-scaling unit 213 of control unit 210 may be configured to scale i using the current scaling factor value to generate κ.sub.ii as shown in FIG. 2.

[0070] A system of inverters with different power ratings, connected in parallel, share the load power in proportion to their ratings if the current gains are chosen as suggested by (42). This directly follows as a consequence of (19) since κ.sub.iP.sub.eq is constant in the parallel configuration, and therefore, P.sub.eq/P.sub.rated is the same for each inverter with identical voltage drops across each output impedance.

[0071] The example AC design specifications shown in Table II require an open-circuit voltage V.sub.oc=126 V. This translates to a voltage scaling factor value κ.sub.v=126 V/V. A rated power and corresponding voltage P.sub.rated=750 W and V.sub.min=114 V results in a current scaling factor κ.sub.i=114/750=0.152 A/A.

[0072] To determine a value of the VO conductance σ and the cubic coefficient α, the closed-form expression for the voltage-regulation characteristic in (19) may be used. Effectively, this design strategy will ensure that the equilibrium RMS terminal voltage of the inverter V.sub.eq is bounded between the limits V.sub.oc≧V.sub.eq≧V.sub.min as the average real-power output P.sub.eq is varied between the limits 0≦P.sub.eq≦P.sub.rated.

[0073] First, notice from (22) that the choice of κ.sub.v in (42a) implies that α is related to σ through

[00035]$\begin{array}{cc}\alpha =\frac{2.Math.\sigma }{3}.& \left(43\right)\end{array}$

Next, substituting P.sub.eq=P.sub.rated and V.sub.eq=V.sub.min into the high-voltage solution of (19) results in

[00036]$\begin{array}{cc}{\stackrel{\_}{V}}_{m.Math..Math.i.Math..Math.n}={{\kappa }_v\left(\frac{\sigma +\sqrt{{\sigma }^2-6.Math.\alpha \left({\kappa }_i/{\kappa }_v\right).Math.{\stackrel{\_}{P}}_{\mathrm{rated}}}}{3.Math.\alpha }\right)}^{\frac{1}{2}}.& \left(44\right)\end{array}$

Substituting κ.sub.v and κ.sub.i from (42) and for α from (43)

[00037]$\begin{array}{cc}{\stackrel{\_}{V}}_{m.Math..Math.i.Math..Math.n}={{\stackrel{\_}{V}}_{\mathrm{oc}}\left(\frac{\sigma +\sqrt{{\sigma }^2-4.Math.\sigma \left({\stackrel{\_}{V}}_{m.Math..Math.i.Math..Math.n}/{\stackrel{\_}{V}}_{\mathrm{oc}}\right)}}{2.Math.\sigma }\right)}^{\frac{1}{2}}.& \left(45\right)\end{array}$

Solving for σ above results in

[00038]$\begin{array}{cc}\sigma =\frac{{\stackrel{\_}{V}}_{\mathrm{oc}}}{{\stackrel{\_}{V}}_{m.Math..Math.i.Math..Math.n}}.Math.\frac{{\stackrel{\_}{V}}_{\mathrm{oc}}^2}{{\stackrel{\_}{V}}_{\mathrm{oc}}^2-{\stackrel{\_}{V}}_{m.Math..Math.i.Math..Math.n}^2}.& \left(46\right)\end{array}$

[0074] Thus, in the example of FIG. 3, the conductance value σ may be determined based on the open circuit voltage and the output voltage of the power electronics that corresponds to the rated real-power output (operation 308). The processor implementing VOC as described herein may be configured to implement a virtual negative conductance element based on the conductance value (operation 310). For example, Control unit 210 may be configured to implement virtual negative conductance element 224 of virtual oscillator 212 with a conductance of 1−/σ.

[0075] Using (43) above, in the example of FIG. 3, the cubic coefficient value α may be determined based on the conductance value (operation 312). The processor implementing VOC as described herein may be configured to implement a cubic voltage-dependent current source based on the cubic coefficient value (operation 314). For example, control unit 210 may be configured to implement virtual cubic voltage-dependent current source 226 of virtual oscillator 212 to multiply the cubed voltage by α.

[0076] The choice of α and σ in (46), respectively, inherently establishes the critical power value P.sub.cr in (20). The design may be iterated if the margin or difference between the rated and critical power values is insufficient.

[0077] Using the example AC design specifications shown in Table II (e.g., the RMS open-circuit and rated-voltage values V.sub.oc, =126 V and V.sub.min=114 V) and applying (46) results in σ=6.09 Ω.sup.−1. From (43), α=4.06 A/V.sup.3.

[0078] FIG. 4 is a graphical plot illustrating an example voltage-regulation curve resulting from the example design specifications of Table II and the voltage-regulation parameters determined therefrom, in accordance with one or more aspects of the present disclosure. Specifically, FIG. 4 shows the equilibrium terminal voltage V.sub.eq as a function of the real-power output P.sub.eq. The inverter is designed to operate in the voltage regime V.sub.oc≧V.sub.eq≧V.sub.min.

[0079] Continuing the design process illustrated in FIG. 3, expressions for i) the frequency-regulation characteristic (23), ii) the rise time (25), and iii) the ratio of amplitudes of the third harmonic to the fundamental (41) may be leveraged to obtain a set of design constraints for the harmonic-oscillator parameter values (e.g., the capacitance C and inductance L).

[0080] In the example of FIG. 3, the capacitance value may be determined (operation 316) For instance, the capacitance value for an inverter may be determined based on (i) a specified rated reactive-power output value for the power electronics, (ii) a specified maximum frequency offset value, (iii) the specified rated real-power output value for the power electronics, (iv) the specified open-circuit voltage value for the power electronics, (v) the specified rated-real-power voltage value for the power electronics, (vi) a specified rise time value, (vii) a specified nominal system frequency value, and (viii) a specified ratio of an amplitude of a third harmonic to an amplitude of a fundamental. Begin with the equilibrium frequency analysis detailed above and the frequency-regulation characteristic in (23). The maximum permissible frequency deviation, denoted by |Δω|.sub.max, is a design input. Substituting for κ.sub.v and κ.sub.i from (42) into (23), and considering the worst-case operating condition for the terminal voltage (which corresponds to consuming or sourcing the maximum reactive power at the minimum permissible terminal voltage, V.sub.min, which is defined in (45)) results in the following lower bound on the capacitance C:

[00039]$\begin{array}{cc}C\ge \frac{1}{2.Math.{.Math.\Delta .Math..Math.\omega .Math.}_{m.Math..Math.a.Math..Math.x.Math.}}.Math.\frac{{\stackrel{\_}{V}}_{\mathrm{oc}}}{{\stackrel{\_}{V}}_{m.Math..Math.i.Math..Math.n.Math.}}.Math.\frac{.Math.{\stackrel{\_}{Q}}_{\mathrm{rated}}.Math.}{{\stackrel{\_}{P}}_{\mathrm{rated}}}=\text{:}.Math.C_{{.Math.\Delta .Math..Math.\omega .Math.}_{\mathrm{ma}.Math..Math.x}}^{m.Math..Math.i.Math..Math.n}& \left(47\right)\end{array}$

where Q.sub.rated is the maximum average reactive power that can be sourced or consumed by the VO-controlled inverter.

[0081] Next, consider the analysis of the (open-circuit) voltage amplitude dynamics detailed above, and the expression for the rise time in (25). With the maximum permissible rise time t.sub.rise.sup.max, serving as design input, the following upper bound for the capacitance C comes from (25) and (46):

[00040]$\begin{array}{cc}C\le \frac{t_{\mathrm{rise}}^{\mathrm{ma}.Math..Math.x}}{6}.Math.\frac{{\stackrel{\_}{V}}_{\mathrm{oc}}}{{\stackrel{\_}{V}}_{m.Math..Math.i.Math..Math.n.Math.}}.Math.\frac{{\stackrel{\_}{V}}_{\mathrm{oc}}^2}{{\stackrel{\_}{V}}_{\mathrm{oc}}^2-{\stackrel{\_}{V}}_{m.Math..Math.i.Math..Math.n}^2}=\text{:}.Math.C_{t_{\mathrm{rise}}}^{\mathrm{ma}.Math..Math.x}.& \left(48\right)\end{array}$

[0082] Finally, consider the harmonics analysis detailed above, and the expression for the ratio of the amplitudes of the third harmonic to the fundamental in (41). With the maximum permissible ratio δ.sub.3:1.sup.max, serving as a design input, an additional lower bound on the capacitance C results from (41) and (46)

[00041]$\begin{array}{cc}C\ge \left(\frac{1}{8.Math.\omega *{\delta }_{3:1}^{\mathrm{ma}.Math..Math.x}}\right).Math.\frac{{\stackrel{\_}{V}}_{\mathrm{oc}}}{{\stackrel{\_}{V}}_{m.Math..Math.i.Math..Math.n}}.Math.\frac{{\stackrel{\_}{V}}_{\mathrm{oc}}^2}{{\stackrel{\_}{V}}_{\mathrm{oc}}^2-{\stackrel{\_}{V}}_{m.Math..Math.i.Math..Math.n}^2}=\text{:}.Math.C_{{\delta }_{3:1}}^{m.Math..Math.i.Math..Math.n}.& \left(49\right)\end{array}$

[0083] A capacitance value C may be determined that satisfies (47)-(49). The processor implementing VOC as described herein may be configured to implement a virtual capacitor based on the capacitance value (operation 318). For example, control unit 210 may be configured to implement virtual capacitor 220 of virtual oscillator 212 using the capacitance value.

[0084] In the example of FIG. 3, the inductance value L may be determined based on the capacitance value and a nominal system frequency (operation 320). That is, the inductance value may follow from rearranging terms in (10)

[00042]$\begin{array}{cc}L=\frac{1}{{C\left({\omega }^{*}\right)}^2}.& \left(50\right)\end{array}$

[0085] The processor implementing VOC as described herein may be configured to implement a virtual inductor based on the inductance value (operation 322). For example, control unit 210 may be configured to implement virtual inductor 222 of virtual oscillator 212 using the inductance value.

[0086] Combining (47)-(49) results in the following range in which C must be selected to meet the performance specifications of frequency regulation, rise time, and harmonics

max{C.sub.|Δω|.sub.max.sup.min,C.sub.δ.sub.3:1.sup.min}≦C≦C.sub.t.sub.rise.sup.max.  (51)

[0087] If C.sub.t.sub.rise.sup.max≦max{C.sub.|Δω|.sub.max.sup.min,C.sub.δ.sub.3:1.sup.min}, then it is not possible to simultaneously meet the specifications of frequency regulation, rise time, and harmonics. Therefore, (51) specifies a fundamental tradeoff in specifying performance requirements and designing VO-controlled inverters. In particular, a VO-controlled inverter that offers a short rise time will have a larger frequency offset and harmonic distortion, while a tightly regulated VO-controlled inverter (e.g., smaller frequency offset and harmonic distortion) will have a longer rise time.

[0088] The example AC design specifications of Table II include |Δω|.sub.max=2π0.5 rad/s, t.sub.rise.sup.max=0.2 s, and δ.sub.3:1.sup.max=2%. Substituting these into (47)-(49) results in C.sub.|Δω|.sub.max.sup.min=0.1.759 F, C.sub.t.sub.rise.sup.max=0.2031 F, and C.sub.δ.sub.3:1.sup.min=0.1010 F. Therefore, to meet the performance specifications, a harmonic-oscillator capacitance C must be selected in the range 0.1759 F≦C≦0.2031 F. In this specific example, the rise time specification is prioritized and a value of C=0.1759 F is selected at the lower bound of the specified range. In various examples other AC design specifications may be prioritized. Since ω*=2π60 rad/s, it follows from (50) that L=39.99 μH.

[0089] Note that the voltage- and frequency-regulation specifications for the example discussed with respect to FIG. 3 are provided herein in terms of worst-case limits. Given the ubiquity of droop control in this domain, such specifications could be specified in terms of the active- and reactive-power droop coefficients m.sub.P and m.sub.Q, respectively. Additionally, by leveraging the correspondences in (54) and (55), the design procedure shown in FIG. 3 may be modified to ensure that the VO-controlled inverter mimics a droop-controlled inverter with the specified m.sub.P and m.sub.Q.

[0090] The example operations of FIG. 3 constitute a design procedure through which a VO-controlled inverter may be configured to meet a specified set of AC performance specifications. Such procedure may, in some examples, be performed by design engineers as part of a design process, or may be performed by one or more systems or machines, such as those used for fabrication.

[0091] As one example, a design system may be configured to receive, as input, AC design specifications for a particular device (e.g., an inverter). The design system may perform the operations described with respect to FIG. 3 above, and thereby configure a digital controller or other processor to emulate a virtual oscillator as described herein to control the corresponding device. In this way, the design procedure disclosed herein may facilitate more efficient production of multiple devices having different AC design specifications.

[0092] For resistive distribution lines, droop control linearly trades off the inverter terminal-voltage amplitude versus active power; and inverter frequency versus reactive power. In the context of the notation established above, these linear laws can be expressed as:

V.sub.eq=V.sub.oc+m.sub.PP.sub.eq,  (52)

ω.sub.eq=ω*+m.sub.PQ.sub.eq,  (53)

where m.sub.P<0 is the active-power droop coefficient and m.sub.Q>0 is the reactive-power droop coefficient.

[0093] It has been shown in the literature that the relations in (52) and (53) provide robust performance for various types of line impedances and are thus referred to as universal droop laws. The equilibria of the averaged VOC dynamics in (15) and (16) can be engineered to be in close correspondence with the droop laws in (52) and (53). For instance, a first-order expansion of V.sub.eq (as a function of P.sub.eq) around the open-circuit voltage, V.sub.oc, is of the form (52) with the following choice of m.sub.P:

[00043]$\begin{array}{cc}{m}_p=\frac{{\kappa }_i.Math.{\kappa }_v}{2.Math.\sigma }.Math.{\left({\stackrel{\_}{V}}_{\mathrm{oc}}-\beta .Math..Math.{\stackrel{\_}{V}}_{\mathrm{oc}}^3\right)}^{-1}.& \left(54\right)\end{array}$

This expression can be derived by evaluating dV.sub.eq/dP.sub.eq from (18) at the open-circuit voltage, V.sub.oc. Similarly, by inspecting (23), it can be seen that ω.sub.eq as a function of Q.sub.eq around the open-circuit voltage, V.sub.oc, is of the form (53) with the following choice of m.sub.Q:

[00044]$\begin{array}{cc}{m}_Q=\frac{{\kappa }_i.Math.{\kappa }_v}{2.Math.C.Math.{\stackrel{\_}{V}}_{\mathrm{oc}}^2}.& \left(55\right)\end{array}$

[0094] With the design strategy described above with respect to FIGS. 3 and 4 for parameters C, κ.sub.v, κ.sub.i, α, and σ, it emerges that the voltage-regulation characteristic in (19) and the frequency regulation characteristic in (2 are close to linear over a wide load range. Experimental results using a hardware prototype described below validate this claim, conclusively demonstrating that droop laws are embedded within the equilibria of the nonlinear VOC dynamics. In other words, the techniques described herein are backward compatible with existing droop control as the techniques described herein subsume droop control in sinusoidal steady state.

[0095] The converse scenario is also considered, wherein droop coefficients, m.sub.P and m.sub.Q, are translated into VOC parameters. The choice of κ.sub.v and κ.sub.i (as given by (42)) and the choice of α (as given by (43)) would remain unchanged. With regard to σ, (54) and β from (5) result in

[00045]$\begin{array}{cc}\sigma =-m_Q^{-1}.Math.\frac{{\kappa }_i}{2}.& \left(56\right)\end{array}$

Furthermore, from (55), it can be seen that the choice of capacitance, C would be given by

[00046]$\begin{array}{cc}C=m_Q^{-1}.Math.\frac{{\kappa }_i}{2.Math.{\stackrel{\_}{V}}_{\mathrm{oc}}},& \left(57\right)\end{array}$

while the inductance, L, would still be specified by (50). Limits on C can be considered in a similar fashion as before, if the specification on m.sub.Q is in terms of an upper bound.

[0096] Although correspondences between the quasi-steady-state behavior of the VOC techniques disclosed here and droop control exist as outlined above, their time-domain performance is markedly different. The main advantage of the techniques of the present disclosure is that they provide time-domain control which acts directly on unprocessed AC measurements when controlling the inverter terminal voltage, as evident in (3). This is unlike droop control which processes AC measurements to compute phasor-based quantities, namely real and reactive power, which are then used to update the inverter voltage amplitude and frequency setpoints. Since phasor quantities are not well-defined in real-time, droop control must necessarily employ a combination of low-pass filters, cycle averaging, coordinate transformations, or π/2 delays to compute P.sub.eq and Q.sub.eq in (52) and (53). These filters, which typically have a cutoff frequency in the range of 1 Hz to 15 Hz, act as a bottleneck to control responsiveness which in turn cause a sluggish response. In contrast, the VO-controlled inverters described herein may operate on real-time measurements and respond to disturbances as they occur. To illustrate these concepts, a simple case study is presented below that demonstrates the time-domain performance of VOC and its power-sharing capabilities.

[0097] Consider two identical single-phase inverters connected in parallel through resistors to a parallel R-L load. The simulated time-domain behavior of VOC in such a scenario is shown in FIGS. 5, 6A, and 6B. Also shown, for comparison, is the performance of droop control. For the droop control implementation, the droop laws in resistive networks (52) and (53) are leveraged and the control architecture is adopted from the literature. The virtual oscillator that emulates the regulation characteristics is then derived by using the aforementioned analysis (equations (56) and (57)).

[0098] FIG. 5 is a graphical plot illustrating example synchronization times of droop control and virtual oscillator-based control, in accordance with one or more aspects of the present disclosure. More specifically, FIG. 5 depicts the time it takes for two inverters to synchronize starting from arbitrary initial conditions. Synchronization error, captured from the deviation of the inverter terminal voltages from the average, as a function of time, for VOC is shown by 500 and for droop control by 502. The waveforms shown in FIG. 5 may be obtained from switching-level simulations.

[0099] In the example of FIG. 5, a metric ∥∥.sub.2 is used, where =[.sub.1,.sub.2].sup.T collects terminal voltages at the inverter. The matrix π:=I.sub.2−½1.sub.21.sub.2.sup.T (1.sub.2×2 is the 2×2 identity, and 1.sub.2×1 is the 2×1 vector with all entries equal to one) is the so-called projector matrix, and by construction, it can be seen that π returns a vector where the entries capture deviations from the average of the vector . From FIG. 5, it can be seen that the inverters configured in accordance with the techniques described herein synchronize by around t=0.1 s, while with droop control, the inverters synchronize by t=0.6 s.

[0100] FIGS. 6A and 6B are graphical plots illustrating active- and reactive-power sharing among two inverters for virtual oscillator-based control and droop control, in accordance with one or more aspects of the present disclosure. In FIG. 6A, solid lines 600 and 601 represent the first and second VO controlled inverters, respectively and dotted lines 602 and 603 represent the first and second droop-controlled inverters, respectively. In FIG. 6B, solid lines 610 and 611 represent the first and second VO controlled inverters, respectively and dotted lines 612 and 613 represent the first and second droop-controlled inverters, respectively. FIGS. 6A and 6B show that identical active- and reactive-power sharing is achieved with both control strategies. It is worth noting, however, that VOC reaches steady-state faster than droop control. The R-L load considered in this particular setup has values of R.sub.load=20Ω and L.sub.load=0.1 H with interconnecting conductances g.sub.1,line=5 Ω.sup.−1 and g.sub.2,line=4.sup.Ω−1 respectively for Inverters 1 and 2.

[0101] A laboratory-scale hardware prototype using the virtual oscillator control techniques described herein with respect to FIGS. 2-6 was also created using the design specifications in Table II. The scaling, voltage-regulation, and harmonic-oscillator parameters (e.g., the parameters that completely characterize the virtual oscillator and ensure the VO-controlled inverter satisfies the design specifications) were computed in the running example discussed with respect to FIG. 3, above. For convenience, these values are listed in Table III below.

TABLE-US-00003 TABLE III VOC PARAMETER VALUES Symbol Description Value Units κ.sub.v Voltage scaling factor 126 V/V κ.sub.i Current scaling factor 0.15 A/A σ Conductance 6.09 Ω.sup.−1 α Cubic coefficient 4.06 A/V.sup.3 C Harmonic-oscillator capacitance 0.18 F L Harmonic-oscillator inductance 3.99 × 10.sup.−5 H

[0102] The component values of the inverter LCL filter for the prototype (e.g., corresponding to LCL output filter 204) were L.sub.f=600 μH, C.sub.f=24 μF, and L.sub.g=44 μH, where L.sub.f, C.sub.f, and L.sub.g are the inverter-side inductor, AC-filter capacitor, and grid-side inductor, respectively. The switching frequency of the prototype inverter was T.sub.sw.sup.−1=15 kHz, the dead time was 200 ns, and the three-level unipolar sine-triangle PWM was utilized. In the specific example of the prototype, the nonlinear dynamics of the virtual-oscillator circuit were programmed on a TMS320F28335 microcontroller available from Texas instruments headquartered in Dallas, Tex., USA.

[0103] To implement the virtual oscillator control on a digital control unit, the continuous functions may need to be discretized. To that end, denote the sampling time utilized in the numerical integration by T.sub.s. In this particular implementation, T.sub.s.sup.−1=15 kHz is chosen. To discretize the virtual-oscillator dynamics (3), the trapezoidal rule of integration may be adopted and the following difference equations can be recovered:

[00047]$\begin{array}{cc}v\left[k\right]={\left(1-\frac{T_s.Math.\sigma }{2.Math.C}+\frac{T_s^2}{4.Math.\mathrm{LC}}\right)}^{-1}\left[\left(1+\frac{T_s.Math.\sigma }{2.Math.C}-\frac{T_s^2}{4.Math.\mathrm{LC}}\right).Math.v\left[k-1\right]-\frac{T_s}{C}.Math.{\kappa }_v.Math.i_L\left[k-1\right]-\frac{T_s}{2.Math.C}.Math.{\kappa }_v.Math.{\kappa }_i\left(i\left[k\right]+i\left[k-1\right]\right)-\frac{\alpha .Math..Math.T_s}{2.Math.C.Math..Math.{\kappa }_v^2}.Math.\left(v^3\left[k\right]+v^3\left[k-1\right]\right)\right],.Math..Math.i_L\left[k\right]=i_L\left[k-1\right]+\frac{T_s}{2.Math.L.Math..Math.{\kappa }_v}.Math.\left(v\left[k\right]+v\left[k-1\right]\right).& \left(58\right)\end{array}$

where kε.sub.≧0 denotes the kth sampling instance, i[k] is the sampled inverter-output current, i.sub.L[k] is the sampled Van der Poll oscillator inductor current, and [k] is the sampled inverter-terminal voltage. In some examples, the difference equations (58) may not be directly implemented on a digital controller, since they contain an algebraic loop through the cubic term .sup.3 [k]. Therefore, it may be necessary to make a simplifying assumption to eliminate the algebraic loop. There are many approaches to accomplish this task. In the specific example of the prototype inverter described here, the assumption .sup.3 [k]≈.sup.3[k−1] was made, allowing (58) to be approximated as:

[00048]$\begin{array}{cc}v\left[k\right]={\left(1-\frac{T_s.Math.\sigma }{2.Math.C}+\frac{T_s^2}{4.Math.\mathrm{LC}}\right)}^{-1}\left[\left(1+\frac{T_s.Math.\sigma }{2.Math.C}-\frac{T_s^2}{4.Math.\mathrm{LC}}\right).Math.v\left[k-1\right]-\frac{T_s}{C}.Math.{\kappa }_v.Math.i_L\left[k-1\right]-\frac{T_s}{2.Math.C}.Math.{\kappa }_v.Math.{\kappa }_i\left(i\left[k\right]+i\left[k-1\right]\right)-\frac{\alpha .Math..Math.T_s}{2.Math.C.Math..Math.{\kappa }_V^2}.Math.v^3\left[k-1\right]\right],.Math..Math.i_L\left[k\right]=i_L\left[k-1\right]+\frac{T_s}{2.Math.L.Math..Math.{\kappa }_v}.Math.\left(v\left[k\right]+v\left[k-1\right]\right).& \left(59\right)\end{array}$

[0104] The difference equations (59) yield reasonable (albeit approximate) dynamics of the virtual oscillator circuit (e.g., as shown in FIG. 2) and they can be implemented directly on a digital controller. The inverter PWM modulation signal (denoted m in the example of FIG. 2) can be constructed as:

[00049]$\begin{array}{cc}m\left[k\right]:=\frac{v\left[k\right]}{v_{d.Math..Math.c}\left[k\right]},& \left(60\right)\end{array}$

where .sub.dc[k] is the measured value of the inverter DC-bus voltage at the kth sampling instance.

[0105] FIG. 7 is a circuit diagram illustrating another example of a single-phase power-electronic inverter (e.g., inverter 700) implemented using VOC, in accordance with one or more aspects of the present disclosure. Inverter 700 builds upon inverter 200, shown in the example of FIG. 2, and provides an extension that allows recovery of a family of droop relations (e.g., real power can be traded off for frequency and reactive power for voltage, as would be desired in inductive networks).

[0106] In the example of FIG. 7, inverter 700 includes power electronics switches 702, LCL output filter 704, and control unit 710. Control unit 710 is configured to implement virtual oscillator 712, current-scaling unit 713, voltage-scaling unit 714, ADCs 715 and 716, modulation unit 717, and PWM unit 718. Virtual oscillator 712 includes virtual capacitor 720, virtual inductor 722, virtual negative-conductance element 724, and virtual cubic voltage-dependent current source 726. In the example of FIG. 7, components 702, 704, 710, 712, 713, 714, 715, 716, 717, 718, 720, 722, 724, and 726 may have functionality similar to that of components 202, 204, 210, 212, 213, 214, 215, 216, 217, 218, 220, 222, 224, and 226 as shown FIG. 2, respectively.

[0107] In the example of FIG. 7, however, control unit 710 is further configured to implement current-scaling unit 730 and transformation unit 732. Current-scaling unit 730 scales the inductor current i.sub.L by Σκ.sub.v. Transformation unit 732 receives, as inputs, the voltage across virtual capacitor 720 and the current through virtual inductor 722, which have been scaled by κ.sub.v (by voltage-scaling unit 714) and εκ.sub.v (by current-scaling unit 730), respectively. Transformation unit 732 then performs an angular displacement of the virtual oscillator dynamic states and outputs for use in constructing the PWM switching signals for inverter 700. As the angular displacement is varied, the VOC dynamics in sinusoidal steady state vary from mimicking droop laws in resistive networks to those in inductive networks. That is, in essence, the coordinate transformation performed by transformation unit 732 allows for the recovery of a suite of nonlinear supersets to droop control.

[0108] The circuit equations that describe the operation of virtual oscillator 712 follow. Specifically, the dynamics of the current through virtual inductor 722, i.sub.L, and the voltage across virtual capacitor 720, .sub.C, are:

[00050]$\begin{array}{cc}L.Math..Math.\frac{{\mathrm{di}}_L}{\mathrm{dt}}=v_c,C.Math..Math.\frac{{\mathrm{dv}}_c}{\mathrm{dt}}=-\alpha .Math..Math.\frac{v_c^3}{{\kappa }_v^3}+\sigma .Math..Math.v_c-i_L-{\kappa }_i.Math.i.& \left(61\right)\end{array}$

Note that (61) is essentially the same as (3) except that (61) makes reference to .sub.C whereas (3) referenced /κ.sub.v. The definitions of (5) still hold to help simplify notation and thus the defined states x and y are the same as discussed with respect to FIG. 2:

[0109] In the example of FIG. 7, the system dynamics can now be represented on the polar plane as shown in FIG. 8, where √{square root over (2)}V=(x.sup.2+y.sup.2).sup.1/2. The x and y coordinates are orthogonal and can be expressed as

x(t)=√{square root over (2)}V(t)sin(ωt+θ(t))=√{square root over (2)}V(t)sin(φ(τ)),

y(t)=√{square root over (2)}V(t)cos(ωt+θ(t))=√{square root over (2)}V(t)cos(φ(τ)),  (62)

where ω, θ(t), and φ(t) still represent the electrical frequency, the phase offset with respect to ω, and the instantaneous phase angle, respectively.

[0110] With reference to FIG. 8, transformation unit 732 may perform a 2-D rotation to define another set of orthogonal signals, v(t) and .sup.⊥(t), where (t) will subsequently be used to control the inverter terminal voltage. The values of and .sup.⊥ can be written in terms of the states=κ.sub.v.sub.C and y=κ.sub.vεi.sub.L using the rotation matrix, Ξ, as:

[00051]$\begin{array}{cc}\left[\begin{array}{c}v\\ v^{\perp }\end{array}\right]=\left[\begin{array}{cc}\cos .Math..Math.\varphi & \sin .Math..Math.\varphi \\ -\sin .Math..Math.\varphi & \cos .Math..Math.\varphi \end{array}\right]\left[\begin{array}{c}y\\ x\end{array}\right]=\text{:}.Math.\Xi \left[\begin{array}{c}y\\ x\end{array}\right].& \left(63\right)\end{array}$

In effect, the output voltage command is given by the first row of (63) and is now given by =κ.sub.v .sub.C cos φ+κ.sub.vεi.sub.L sin φ. In other words, the output voltage command is now a linear combination of the virtual oscillator capacitor voltage and inductor current.

[0111] To extract amplitude and phase information, the voltage signal, , and its corresponding orthogonal signal, .sup.⊥, are transformed to polar coordinates as: =√{square root over (2)}V cos(φ+ω), and .sup.⊥=√{square root over (2)}V sin(φ+φ). Differentiating the identities:

[00052]$\sqrt{2}.Math.V=\sqrt{x^2+y^2}.Math..Math.\mathrm{and}.Math..Math.\phi ={\tan }^{-1}.Math.\frac{x}{y},$

the following dynamical equations for V and φ are obtained:

[00053]$\begin{array}{cc}\stackrel{.}{V}=\frac{\mathrm{.Math.}}{\sqrt{2}}.Math.\left(\sigma .Math..Math.g\left(\sqrt{2}.Math.V.Math..Math.\cos \left(\phi \right)\right)-{\kappa }_v.Math.{\kappa }_i.Math.i\right).Math.\cos .Math..Math.\left(\phi \right),.Math.\stackrel{.}{\phi }=1-\frac{\mathrm{.Math.}}{\sqrt{2}.Math.V}.Math.\left(\sigma .Math..Math.g\left(\sqrt{2}.Math.V.Math..Math.\cos \left(\phi \right)\right)-{\kappa }_v.Math.{\kappa }_i.Math.i\right).Math.\sin \left(\phi \right).& \left(64\right)\end{array}$

[0112] Again, the parametric regime characterized by ε0 is focused on to ensure that dynamics mimic those of a harmonic oscillator. Furthermore, the dynamics in (64) are averaged to focus on AC-cycle time scales. The resulting autonomous system is far easier to analyze than the time-varying system (without compromising much on the accuracy as the analysis is still correct up to (ε)). V and θ are denoted to be 2π-averaged values of V and θ. Then, recognizing that ε0, the averaged dynamics can be written as:

[00054]$\begin{array}{cc}\left[\begin{array}{c}\stackrel{\stackrel{.}{\_}}{V}\\ \stackrel{\stackrel{.}{\_}}{\theta }\end{array}\right]=\frac{\mathrm{.Math.}.Math..Math.\sigma }{2}\left[\begin{array}{c}\stackrel{\_}{V}-\frac{\beta }{2}.Math.{\stackrel{\_}{V}}^3\\ 0\end{array}\right]-\frac{\mathrm{.Math.}.Math..Math.{\kappa }_v-{\kappa }_i}{2.Math.\pi .Math.\sqrt{2}}.Math.{\int }_0^{2.Math.\pi }.Math.i\left[\begin{array}{c}\cos \left(\tau +\stackrel{\_}{\theta }\right)-\\ \frac{1}{\stackrel{\_}{V}}.Math.\sin \left(\tau +\stackrel{\_}{\theta }\right)\end{array}\right].Math.d.Math..Math.\tau .& \left(65\right)\end{array}$

[0113] The relationship between terminal voltage and power is of interest. To that end, the instantaneous active- and reactive-power injections are defined as:

[00055]$\begin{array}{cc}P\left(t\right)=v\left(t\right).Math.i\left(t\right),Q\left(t\right)=v\left(t-\frac{\pi }{2}\right).Math.i\left(t\right)=\text{:}.Math.v^{\perp }\left(t\right).Math.i\left(t\right),& \left(66\right)\end{array}$

where i(t) is the output current. The average real and reactive power over an AC cycle of period 2π/ω are:

[00056]$\begin{array}{cc}\stackrel{\_}{P}=\frac{\omega }{2.Math.\pi }.Math.{\int }_{s=0}^{2.Math.\pi /\omega }.Math.P\left(s\right).Math.\mathrm{ds},\stackrel{\_}{Q}=\frac{\omega }{2.Math.\pi }.Math.{\int }_{s=0}^{2.Math.\pi /\omega }.Math.Q\left(s\right).Math.\mathrm{ds}.& \left(67\right)\end{array}$

[0114] Transitioning (65) from τ to t coordinates, retaining only (ε) terms, and using the definitions of active and reactive power in (66) results in:

[00057]$\begin{array}{cc}\frac{d}{\mathrm{dt}}.Math.\stackrel{\_}{V}=\frac{\sigma }{2.Math.C}.Math.\left(\stackrel{\_}{V}-\frac{\beta }{2}.Math.{\stackrel{\_}{V}}^3\right)-\frac{{\kappa }_V.Math.{\kappa }_i}{2.Math.C}.Math.\left(\cos .Math..Math.\phi .Math..Math.\frac{\stackrel{\_}{P}}{V}+\sin .Math..Math.\phi .Math..Math.\frac{\stackrel{\_}{P}}{V}\right),.Math.\frac{d}{\mathrm{dt}}.Math.\stackrel{\_}{\theta }=\frac{{\kappa }_v.Math.{\kappa }_i}{2.Math.C}.Math.\left(-\sin .Math..Math.\phi .Math..Math.\frac{\stackrel{\_}{P}}{{\stackrel{\_}{V}}^2}+\cos .Math..Math.\phi .Math..Math.\frac{\stackrel{\_}{P}}{{\stackrel{\_}{V}}^2}\right)& \left(68\right)\end{array}$

Thus, the terminal-voltage amplitude and frequency dynamics that are linked to active- and reactive-power outputs at the inverter terminals through the 2-D rotation matrix, Ξ, are obtained. The steady-state voltage- and frequency-regulation equations are equilibrium points of the system in (68):

[00058]$\begin{array}{cc}\sigma .Math..Math.{\stackrel{\_}{V}}_{\mathrm{eq}}-\frac{\sigma .Math..Math.\beta }{2}.Math.{\stackrel{\_}{V}}_{\mathrm{eq}}^3-{\kappa }_v.Math.{{\kappa }_i\left[\begin{array}{c}\cos .Math..Math.\varphi \\ \sin .Math..Math.\varphi \end{array}\right]}^T\left[\begin{array}{c}\frac{{\stackrel{\_}{P}}_{\mathrm{eq}}}{{\stackrel{\_}{V}}_{\mathrm{eq}}}\\ \frac{{\stackrel{\_}{Q}}_{\mathrm{eq}}}{{\stackrel{\_}{V}}_{\mathrm{eq}}}\end{array}\right]=0,.Math.{\omega }_{\mathrm{eq}}-\omega --.Math.{\frac{{\kappa }_v.Math.{\kappa }_i}{2.Math.C}\left[\begin{array}{c}-\sin .Math..Math.\varphi \\ \cos .Math..Math.\varphi \end{array}\right]}^T\left[\begin{array}{c}\frac{{\stackrel{\_}{P}}_{\mathrm{eq}}}{{\stackrel{\_}{V}}_{\mathrm{eq}}^2}\\ \frac{{\stackrel{\_}{Q}}_{\mathrm{eq}}}{{\stackrel{\_}{V}}_{\mathrm{eq}}^2}\end{array}\right]=0.& \left(69\right)\end{array}$

with V.sub.eq, P.sub.eq, and Q.sub.eq denoting RMS voltage amplitude, average real- and reactive-power output respectively, and ω.sub.eq denoting the steady-state frequency. With appropriate decoupling assumptions, various droop-like laws can be recovered from (69): with φ=π/2, terminal voltage is traded off for reactive power and frequency is traded off for active power; similarly, with φ=0, the terminal voltage is traded off for active power and frequency is traded off for reactive power. These are reminiscent of the ubiquitous droop control laws for inductive (and resistive, respectively) line impedences.

[0115] FIG. 9 is a flow diagram illustrating additional example operations for designing a VOC system, in accordance with one or more aspects of the present disclosure. The operations shown in FIG. 9 are described with reference to inverter 700 as shown in FIG. 7. FIG. 9 represents only one example method for determining virtual oscillator parameter values such that a VO-controlled inverter (e.g., inverter 700 of FIG. 7) satisfies a set of AC performance specifications. In other examples, the design procedure may include additional or different operations.

[0116] The averaged models described with respect to FIG. 7 above reveal conclusive links between real- and reactive-power outputs and the terminal-voltage dynamics. Therefore, they suggest avenues to design the virtual oscillators so that the controlled inverter (e.g., inverter 700) satisfies performance requirements. With reference to the model in (69), the operations of FIG. 3 provide an example design procedure for the case where φ=0. The example operations of FIG. 9 explore an analogous design strategy for the case φ=π/2, which is applicable to inductive networks.

[0117] Accompanying the design strategy to pick system parameters is a running example for VOC design corresponding to the following set of AC performance specifications: rated real power 750 W, rated reactive power 750 VAR, voltage-regulation of ±5% around a nominal voltage of 120 V, frequency-regulation of ±0.5 Hz around a nominal frequency of 60 Hz, rise-time of 0.1 s, ratio of third-to-first harmonic 2% and load sharing in a parallel setting.

[0118] With the choice φ=π/2, the expressions in (69) boil down to the following:

[00059]$\begin{array}{cc}\sigma .Math..Math.{\stackrel{\_}{V}}_{\mathrm{eq}}-\frac{\sigma .Math..Math.\beta }{2}.Math.{\stackrel{\_}{V}}_{\mathrm{eq}}^3-{\kappa }_v.Math.{\kappa }_i.Math.\frac{Q}{{\stackrel{\_}{V}}_{\mathrm{eq}}}=0,& \left(70\right)\\ {\omega }_{\mathrm{eq}}-\omega +\frac{{\kappa }_v.Math.{\kappa }_i}{2.Math.C}.Math.\frac{\stackrel{\_}{P}}{{\stackrel{\_}{V}}_{\mathrm{eq}}^2}=0& \left(71\right)\end{array}$

[0119] In the example of FIG. 9, a value of the voltage scaling factor, κ.sub.v, may be determined based on the open-circuit voltage of the power electronics (operation 900). To standardize design, κ.sub.v is again chosen so that when the voltage across virtual capacitor 720 is 1 V RMS, the output voltage is the open-circuit voltage, V.sub.oc. A processor implementing VOC as described herein may be configured to scale an output voltage of the virtual oscillator using the voltage scaling factor value (operation 902).

[0120] A value of the current scaling factor, κ.sub.i, may be determined based on the rated reactive-power output and corresponding output voltage of the power electronics (operation 904). In this example, κ.sub.i is chosen so that the output current is 1 A when rated reactive power, Q.sub.rated, is drawn from it. Thus, the choice of κ.sub.v and κ.sub.i are taken to be:

[00060]$\begin{array}{cc}{\kappa }_v={\stackrel{\_}{V}}_{\mathrm{oc}},{\kappa }_i=\frac{V_{m.Math..Math.i.Math..Math.n}}{.Math.{\stackrel{\_}{Q}}_{\mathrm{rated}}.Math.}.& \left(72\right)\end{array}$

The processor implementing VOC as described herein may be configured to generate a virtual oscillator input current by scaling the power electronics output current using the current scaling factor value (operation 906). For instance, current-scaling unit 713 of control unit 710 may be configured to scale i using the current scaling factor value to generate κ.sub.ii as shown in FIG. 7.

[0121] For the specified voltage and current scaling factor values, σ can be calculated in terms of V.sub.min, which is the inverter terminal voltage when rated reactive power is drawn from it. Rearranging the terms and substituting values of κ.sub.v and κ.sub.i in voltage regulation equation gives

[00061]$\begin{array}{cc}\sigma =\frac{{\stackrel{\_}{V}}_{\mathrm{oc}}^3}{{\stackrel{\_}{V}}_{m.Math..Math.i.Math..Math.n}\left({\stackrel{\_}{V}}_{\mathrm{oc}}^2-{\stackrel{\_}{V}}_{m.Math..Math.i.Math..Math.n}^2\right)}.& \left(73\right)\end{array}$

Thus, in the example of FIG. 9, the conductance value σ may be determined based on the open circuit voltage and the output voltage of the power electronics that corresponds to the rated reactive-power output (operation 908). The processor implementing VOC as described herein may be configured to implement a virtual negative conductance element based on the conductance value (operation 910). For example, Control unit 710 may be configured to implement virtual negative conductance element 724 of virtual oscillator 712 with a conductance of −1/σ. In the running example described here, upon substituting V.sub.oc=120 V and V.sub.min=114 V, the result is σ=10.79.

[0122] Similar to operation 312 of FIG. 3, and again using (43) above, the cubic coefficient value α may, in the example of FIG. 9, be determined based on the conductance value (operation 912). The processor implementing VOC as described herein may be configured to implement a cubic voltage-dependent current source based on the cubic coefficient value (operation 914). For example, control unit 710 may be configured to implement virtual cubic voltage-dependent current source 726 of virtual oscillator 712 to multiply the cubed voltage by α.

[0123] In the example of FIG. 9, the capacitance value may be determined (operation 916). To quantify the dynamic response of the VO-controlled inverter, the terminal voltage rise time to the open-circuit voltage can be analyzed. As before, the rise time, t.sub.rise is defined as the time taken for the inverter terminal voltage to rise from 0.1V.sub.oc to 0.9V.sub.oc. Furthermore, the harmonic content of the terminal voltage is captured and quantified through the ratio of the amplitude of the third harmonic to the fundamental, δ.sub.3:1. Following the analysis with respect to FIGS. 2 and 3 above for rise-time and harmonic content in the terminal voltage for an unloaded inverter,

[00062]$\begin{array}{cc}{t}_{\mathrm{rise}}\approx \frac{6.Math.C}{\sigma },{\delta }_{3:1}\approx \frac{\sigma }{8.Math.\omega .Math..Math.C}& \left(74\right)\end{array}$

where ω is the nominal frequency. Thus, upon imposing the design requirements (rise-time of 0.1 s, ratio of third-to-first harmonic 2%) results in the following from (74):

[00063]$\begin{array}{cc}\frac{10}{192.Math.\pi }\le \frac{C}{\sigma }\le \frac{1}{60}.& \left(75\right)\end{array}$

[0124] Also, to ensure frequency-regulation of ±0.5 Hz, the following must hold:

[00064]$\begin{array}{cc}C\ge \frac{1}{4.Math.\mathrm{\pi 0}.Math.\mathrm{.5}}.Math.\frac{{\stackrel{\_}{V}}_{\mathrm{oc}}}{{\stackrel{\_}{V}}_{m.Math..Math.i.Math..Math.n.Math.}}.Math.\frac{{\stackrel{\_}{P}}_{\mathrm{rated}}}{.Math.{\stackrel{\_}{Q}}_{\mathrm{rated}}.Math.},& \left(76\right)\end{array}$

where κ.sub.v and κ.sub.i from (72) have been substituted in (71) and where the worst-case operating condition is assumed (where maximum average power, P.sub.rated, is sourced at the minimum permissible terminal voltage, V.sub.min).

[0125] The processor implementing VOC as described herein may be configured to implement a virtual capacitor based on the capacitance value (operation 918). For example, control unit 710 may be configured to implement virtual capacitor 720 of virtual oscillator 712 using the capacitance value. For the specific example discussed here, a value of C=1019/60 F may be selected to satisfy the lower bound specified by (76).

[0126] Similar to operation 320 of FIG. 3, and again using (50) above, the inductance value L may be determined based on the capacitance value and a nominal system frequency (operation 920). The processor implementing VOC as described herein may be configured to implement a virtual inductor based on the inductance value (operation 922). For example, control unit 710 may be configured to implement virtual inductor 722 of virtual oscillator 712 using the inductance value.

[0127] The example operations of FIG. 9 constitute a design procedure through which a VO-controlled inverter may be configured to meet a specified set of AC performance specifications for inductive networks. Such procedure may, in some examples, be performed by design engineers as part of a design process, or may be performed by one or more systems or machines, such as those used for fabrication.

[0128] As one example, a design system may be configured to receive, as input, AC design specifications for a particular device (e.g., an inverter). The design system may perform the operations described with respect to FIG. 9 above, and thereby configure a digital controller or other processor to emulate a virtual oscillator as described herein to control the corresponding device. In this way, the design procedure disclosed herein may facilitate more efficient production of multiple devices having different AC design specifications.

[0129] FIG. 10 is a circuit diagram illustrating another example of a single-phase power-electronic inverter (e.g., inverter 1000) implemented using VOC, in accordance with one or more aspects of the present disclosure. Inverter 1000 builds upon inverter 200, shown in the example of FIG. 2, and provides an extension that allows for global asymptotic synchronization of inverters in uniform linear, time-invariant (LTI) electrical networks with arbitrary topologies. Inverter 1000 uses a particular feedback scheme. This implementation thereby provides a closed-loop system that is equivalent to a static diffusive interconnection. This allows for tailoring incremental-stability arguments to the particular setting of uniform electrical networks.

[0130] In the example of FIG. 10, inverter 1000 includes power electronics switches 1002, LCL output filter 2004, and control unit 1010. Control unit 1010 is configured to implement virtual oscillator 1012, current-scaling unit 1013, voltage-scaling unit 1014, ADCs 1015 and 1016, modulation unit 1017, and PWM unit 1018. Virtual oscillator 1012 includes virtual capacitor 1020, virtual inductor 1022, virtual negative-conductance element 1024, and virtual cubic voltage-dependent current source 1026. In the example of FIG. 10, components 1002, 1004, 1010, 1012, 1013, 1014, 1015, 1016, 1017, 1018, 1020, 1022, 1024, and 1026 may have functionality similar to that of components 202, 204, 210, 212, 213, 214, 215, 216, 217, 218, 220, 222, 224, and 226 as shown FIG. 2, respectively.

[0131] Of note, in the example of FIG. 10, LCL output filter 1004 includes two resistors, each having a resistance of r.sub.f/2. These resistances represent the intrinsic line resistance of the respective branch. These resistances may also be present in one or more other example inverters described herein but, for brevity, are not shown.

[0132] In the example of FIG. 10, control unit 1010 is further configured to implement voltage-scaling unit 1040 and summing unit 1042. Voltage-scaling unit 1040 is configured to scale an output voltage .sub.C of virtual oscillator 1012 by a voltage scaling factor, γ. This voltage scaling factor is a constant that is determined based on the inductance value L.sub.f and resistance value r.sub.f of the inverter-side branch for inverter 1000 (e.g., output filter 1004). Specifically,

[00065]$\begin{array}{cc}\gamma \ge \frac{r_f}{L_f}.& \left(77\right)\end{array}$

[0133] Summing unit 1042 adds the resulting scaled voltage, γ.sub.C, to the current through virtual capacitor 1020, i.sub.C and provides the resulting sum, γ.sub.C+i.sub.C to voltage-scaling unit 1014.

[0134] With the additional implementation of voltage-scaling unit 1040 and summing unit 1042, inverter 1000 can be shown to be globally asymptotically stable. That is, inverter 1000 may achieve global asymptotic synchronization.

[0135] FIG. 11 is a circuit diagram illustrating another example of a single-phase power-electronic inverter (e.g., inverter 1100) implemented using VOC, in accordance with one or more aspects of the present disclosure. Inverter 1100 builds upon inverter 200, shown in the example of FIG. 2, and provides an extension that allows for power-tracking capability. That is, inverter 1100 may be able to regulate the power output to match that available from the primary DC source. Inverter 1100 may adaptively adjust the parameters of the underlying VOC to regulate its output power via feedforward linearization.

[0136] In the example of FIG. 11, inverter 1100 includes power electronics switches 1102, LCL output filter 1104, and control unit 1110. Control unit 1110 is configured to implement virtual oscillator 1112, current-scaling unit 1113, voltage-scaling unit 1114, ADCs 1115 and 1116, modulation unit 1117, and PWM unit 1118. Virtual oscillator 1112 includes virtual capacitor 1120, virtual variable inductor 1122, virtual negative-conductance element 1124, and virtual variable cubic voltage-dependent current source 1126. Components 1102, 1104, 1110, 1112, 1113, 1114, 1115, 1116, 1117, 1118, 1120, and 1124 may have functionality similar to that of components 202, 204, 210, 212, 213, 214, 215, 216, 217, 218, 220, and 224 as shown FIG. 2, respectively.

[0137] In the example of FIG. 11, however, control unit 1110 is configured to implement the virtual inductor and virtual cubic voltage-dependent current source within the virtual oscillator 1112 in a variable fashion (e.g., virtual variable inductor 1127 and virtual variable cubic voltage-dependent current source 1128). Furthermore, control unit 1110 is configured to also implement power tracking control unit 1150 to control the variable nature of virtual variable inductor 1127 and virtual variable cubic voltage-dependent current source 1128. This may allow inverter 1100 to track a specified real and reactive power reference while maintaining its “grid-forming” behaviors. These additional features may be necessary, in some cases, to match the capabilities of the primary DC energy source (such as a photovoltaic array or battery).

[0138] The addition of power-tracking control unit 1150 provides an outer control loop that adaptively modifies the parameters of the VOC control techniques to adjust the output power of inverter 1100. In addition, a feedforward linearization is used to linearize the power-tracking dynamics and to allow power-tracking control to be designed with specified dynamic behavior. A closed-form procedure to design a power-tracking VOC inverter such as inverter 1100 to meet specified stability margins is further described below with respect to FIG. 12.

[0139] As a corollary to the analysis of FIG. 2, and with a few slight changes in notation, it can be shown by singular-perturbation and averaging techniques that the AC voltage magnitude and phase of = in inverter 200 follow the dynamics:

[00066]$\begin{array}{cc}\frac{d}{\mathrm{dt}}.Math.\stackrel{\_}{V}=\frac{\sigma }{2.Math.C}.Math.\left(\stackrel{\_}{V}-\frac{\sigma }{V_{\mathrm{oc}}^2}.Math.{\stackrel{\_}{V}}^3\right)-\frac{{\kappa }_v.Math.{\kappa }_i}{2.Math.C.Math..Math.\stackrel{\_}{V}}.Math.\stackrel{\_}{P}& \left(78\right)\\ \frac{d}{\mathrm{dt}}.Math.\stackrel{\_}{\theta }={\omega }_{\mathrm{oc}}-\omega +\frac{{\kappa }_v.Math.{\kappa }_i}{2.Math.C.Math..Math.{\stackrel{\_}{V}}^{2.Math.}}.Math.\stackrel{\_}{Q},& \left(79\right)\end{array}$

where V and θ are the AC-cycle-average inverter voltage magnitude and synchronous phase, P and Q are the inverter real and reactive output power, and ω is the system frequency. The open-circuit parameters V.sub.oc and ω.sub.oc are found by:

[00067]$\begin{array}{cc}{V}_{\mathrm{oc}}={\kappa }_v.Math.\sqrt{\frac{2.Math.\sigma }{3.Math.\alpha }},& \left(80\right)\\ {\omega }_{\mathrm{oc}}=\frac{1}{\sqrt{\mathrm{LC}}}.& \left(81\right)\end{array}$

[0140] The equilibria of (78) and (79) were also compared above to the standard universal droop

V.sub.eq=V.sub.oc+m.sub.PP.sub.eq,  (82)

ω.sub.eq=ω.sub.oc+m.sub.QQ.sub.eq,  (83)

where V.sub.eq and ω.sub.eq are the equilibrium voltage and frequency, m.sub.P<0 is the active-power droop coefficient, and m.sub.Q>0 is the reactive-power droop coefficient. It can further be shown that the equilibria of (78) and (79) track (82) and (83) with less than 1% error for the following selections of m.sub.P and m.sub.Q:

[00068]$\begin{array}{cc}{m}_P=\frac{{\kappa }_v.Math.{\kappa }_i}{2.Math.\sigma }.Math.{V_{\mathrm{Nom}}\left(1-2.Math.\sigma \right)}^{-1},& \left(84\right)\\ m_Q=\frac{{\kappa }_v.Math.{\kappa }_i}{2.Math.{\mathrm{CV}}_{\mathrm{Nom}}^2},& \left(85\right)\end{array}$

where V.sub.Nom is the system nominal AC voltage.

[0141] In summary, while the relationship between the fundamental VOC parameters σ, α, L, and C) to the equilibrium states (V.sub.eq, ω.sub.eq, P.sub.eq, Q.sub.eq) are highly nonlinear, the open-circuit parameters V.sub.oc and ω.sub.oc provide a method of establishing a simple linear relationship between the VOC characteristics and the equilibrium states. These relationships form the basis of the VOC-based power-tracking control techniques described with respect to FIG. 11 herein.

[0142] In some examples, a VOC-based power-tracking inverter such as inverter 1100 may be based on a feedforward-linearization approach. One purpose of these techniques may be to allow a VOC-controlled inverter to track a reference of real and reactive power while maintaining the “grid-forming” nature of VOC. This can be accomplished by a real-time modulation of the VOC fundamental parameters α and L in order to regulate inverter real and reactive output power.

[0143] Because the relationship of the fundamental parameters with real and reactive power P.sub.eq and Q.sub.eq is highly nonlinear (as shown in (78) and (79)), it may be difficult to control real and reactive power with them directly. However, the open-circuit parameters V.sub.oc and ω.sub.oc are linearly related to P.sub.eq and Q.sub.eq by (82) and (83). These relationships allow the design of linear control techniques with feedforward to determine V.sub.oc and ω.sub.oc to achieve the desired P.sub.eq and Q.sub.eq. The VOC fundamental parameters σ, α, L, and C can then be calculated from V.sub.oc and ω.sub.oc by (80) and (81).

[0144] The feedforward-linearization approach used by inverter 1100 is based on the relationships described above. The control signals ΔV.sub.oc and Δω.sub.oc may be defined as follows:

ΔV.sub.oc:=V.sub.oc−V.sub.Nom,  (86)

Δω.sub.oc:=ω.sub.oc−ω.sub.Nom,  (86)

where V.sub.Nom and ω.sub.Nom are the nominal voltage and frequency at the inverter terminals. Substituting (86) and (87) into (82) and (83) and solving for P.sub.eq and Q.sub.eq results in the following equilibrium relationships:

P.sub.eq:=−m.sub.P.sup.−1ΔV.sub.oc+m.sub.P.sup.−1(V.sub.eq−V.sub.Nom),  (88)

Q.sub.eq:=−m.sub.Q.sup.−1Δω.sub.oc+m.sub.Q.sup.−1(ω.sub.eq−ω.sub.Nom).  (89)

[0145] The feedfoward linearization is accomplished by application of (80) and (81) to transform the open-circuit parameters V.sub.oc and ω.sub.oc (which are linearly related to P.sub.eq and Q.sub.eq) into the VOC fundamental parameters σ, α, L, and C. To completely define the four VOC fundamental parameters (80) and (81) must be supplemented with two additional algebraic constraints. In addition, from (80) and (81) it is clear that adjusting the values of σ and C will affect the droop slopes m.sub.P and m.sub.Q, which are desired to be fixed (since these linear relationships will form the basis of the linearized control).

[0146] By performing the VOC design procedure detailed with respect to FIG. 3 above, a “nominal” set of fixed VOC parameters σ.sub.Nom, α.sub.Nom, L.sub.Nom, and C.sub.Nom are obtained that achieve a set of specified AC performance targets. The value of σ may be fixed at σ.sub.Nom in order to fix the droop slope m.sub.P, and similarly the value of C may be fixed at C.sub.Nom in order to fix m.sub.Q. That is, in FIG. 11, the conductance of virtual negative conductance element 1124 and the capacitance of virtual capacitance 1120 may be fixed at −1/σ.sub.Nom and C.sub.Nom, respectively. Therefore, solving (80) and (81) for α and L respectively for fixed σ and C, results in:

[00069]$\begin{array}{cc}\alpha =\frac{2}{3}.Math.\frac{{\kappa }_v^2}{{\stackrel{\_}{V}}_{\mathrm{oc}}^2.Math.{\sigma }_{\mathrm{Nom}}}=\text{:}.Math.f_{\alpha }\left({\stackrel{\_}{V}}_{\mathrm{oc}}\right),& \left(90\right)\\ L=\frac{1}{{\omega }_{\mathrm{oc}}^2.Math.C_{\mathrm{Nom}}}=\text{:}.Math.f_L\left({\omega }_{\mathrm{oc}}\right).& \left(91\right)\end{array}$

The nonlinear functions f.sub.α(V.sub.oc) and f.sub.L(ω.sub.oc) above allow the VOC fundamental parameters α and L to be calculated from the control variables V.sub.oc and ω.sub.oc. In addition, since V.sub.oc and ω.sub.oc are linearly related to P.sub.eq and Q.sub.eq, (90) and (91) constitute a feedforward linearization of the VOC real and reactive power dynamics.

[0147] In the example of FIG. 11, control unit 1110 operates using techniques similar to the VOC control techniques described with respect to FIG. 2 above, but with dynamically-defined parameters α and L for virtual variable inductance 1127 and virtual variable cubic voltage-dependent current source 1128, respectively. σ and C are fixed at their nominal values.

[0148] Power tracking control unit 1150, in the example of FIG. 11, includes power calculation unit 1152, power feedback filters 1153 and 1154, feedforward units 1155 and 1156, linear compensators 1157 and 1158, and nonlinear operators 1159 and 1160. Power tracking unit 1150 acts as an outer loop to the VOC and adaptively modulates a and L to regulate P and Q to track the specified P.sub.ref and Q.sub.ref.

[0149] Power calculation unit 1152 receives the inverter output current and the scaled output voltage of virtual oscillator 1112 and determines the current real and reactive powers being output by inverter 1100. For a single-phase inverter, the measured real and reactive power P and Q are calculated from and i of the VOC as follows:

P=i−Hil()Hil(i)  (92)

P=Hil()i+Hil(i)  (93)

where Hil(•) represents the Hilbert transform for generation of an orthogonal signal. For a three-phase inverter, (92) and (93) can be replaced by a three-phase power calculation method.

[0150] In the example of FIG. 11, power feedback filters 1153 and 1154, implementing H.sub.FB,P and H.sub.FB,Q, respectively, may receive the respective power measurements from power calculation unit 1152 and filter out harmonic content therefrom. As one example, H.sub.FB,P=H.sub.FB,Q may be selected as a moving average filter with a window equal to one nominal AC cycle period (T.sub.AC=2π/ω.sub.Nom). Linear compensators 1157 and 1158 (G.sub.P and G.sub.Q, respectively) may operate on the real and reactive power error signals (e.g., P.sub.ref−{circumflex over (P)} and Q.sub.ref−{circumflex over (Q)}, respectively) to produce, after summation with a feedforward term, the control signals ΔV.sub.oc and Δω.sub.oc. Further details of linear compensators 1157 and 1158 are provided below with respect to FIG. 12. The control signals ΔV.sub.oc and Δω.sub.oc are then summed with their respective nominal values to generate V.sub.oc and ω.sub.oc.

[0151] Nonlinear operators 1159 and 1160 respectively implement functions f.sub.α(V.sub.oc) and f.sub.L(ω.sub.oc) as defined in (90) and (91) to provide feedforward linearization and generate the VOC fundamental parameters α and L. These are then dynamically fed into virtual variable inductance 1127 and virtual variable cubic voltage-dependent current source 1128, respectively.

[0152] Feedforward units 1157 and 1158 multiply their respective inputs by the terms −m.sub.P and −m.sub.Q, respectively. −m.sub.P and −m.sub.Q represent the nominal relationship of voltage to real power and frequency to reactive power respectively. Rearranging (82) and (83) substituting P=P.sub.ref, Q=Q.sub.ref, V.sub.eq=V.sub.Nom, and ω.sub.eq=ω.sub.Nom results in

ΔV.sub.ff:=V.sub.oc−V.sub.Nom=−m.sub.PP.sub.ref,  (94)

Δω.sub.ff:=ω.sub.oc−ω.sub.Nom=−m.sub.QQ.sub.ref.  (95)

[0153] The feedforward linearization technique described herein decouples the nonlinearity of the relationship between the VOC fundamental parameters and its equilibrium states, thereby allowing power-tracking control to be designed based on the linear relationships in (82) and (83). However, (82) and (83) are algebraic equations describing the VOC equilibrium response, not dynamic equations describing its transient response. In order to design power-tracking control to meet desired time-domain response specifications (such as rise time, overshoot, etc.), it may be necessary to model the (feedforward-linearized) VOC real and reactive power dynamics.

[0154] It has been shown with respect to FIGS. 2 and 3 above that a VOC-controlled inverter exhibits a voltage rise time t.sub.rise that can be approximated as

[00070]$\begin{array}{cc}{t}_{\mathrm{rise}}\approx 6.Math..Math.\frac{C}{\sigma }.& \left(96\right)\end{array}$

It was further observed from simulation and experiment that the VOC power dynamics exhibit a first-order response with similar rise time. Thus the dynamics of the relationship between the control variables ΔV.sub.oc and Δω.sub.oc and the inverter real and reactive power P and Q may be approximated by substituting (86) and (87) into (82) and (83) and assuming a first order response with rise time t.sub.rise as follows:

[00071]$\begin{array}{cc}\stackrel{\_}{P}\left(s\right)=-m_P^{-1}.Math.\frac{{\omega }_c}{s+{\omega }_c}.Math.\left(\Delta .Math..Math.{\stackrel{\_}{V}}_{\mathrm{oc}}\left(s\right)+\stackrel{\_}{V}\left(s\right)-{\stackrel{\_}{V}}_{\mathrm{Nom}}\right),& \left(97\right)\\ \stackrel{\_}{Q}\left(s\right)=-m_Q^{-1}.Math.\frac{{\omega }_c}{s+{\omega }_c}.Math.\left(\Delta .Math..Math.{\omega }_{\mathrm{oc}}\left(s\right)+\omega \left(s\right)-{\omega }_{\mathrm{Nom}}\right),& \left(98\right)\end{array}$

where ω.sub.c=2π/(3t.sub.rise) is the cutoff frequency of the first-order response and s is the Laplace variable. In the calculation of ω.sub.c, the rise time of the first-order filter has been assumed to be equal to three times its time constant.

[0155] FIG. 12 is a block diagram illustrating one example of a dynamic model of a power-tracking loop (e.g., based on the feedforward-linearized model of the VOC power response in (97) and (98)), in accordance with one or more aspects of the present disclosure. FIG. 12 includes power feedback filters 1253 and 1254, feedforward units 1255 and 1256, linear compensators 1257 and 1258, and nonlinear operators 1259 and 1260. Components 1253-1260 represent operations similar to the functionality of components 1153-1160 shown in FIG. 11, respectively.

[0156] A first step in one example design procedure for VOC-based power-tracking control aiming to achieve a set of specified AC and dynamic performance targets is the selection of a set of “nominal” VOC fundamental parameters (σ.sub.Nom, α.sub.Nom, L.sub.Nom, and C.sub.Nom) to achieve the specified AC performance criteria. This set of parameters represents the value of the VOC parameters when ΔV.sub.oc=Δω.sub.oc=0. Table IV defines an example set of AC and dynamic performance targets.

TABLE-US-00004 TABLE IV EXAMPLE AC AND DYNAMIC PERFORMANCE SPECIFICATIONS Symbol Description Value Units P.sub.rated Inverter Rated Real Power 750 W |Q.sub.rated| Inverter Rated Reactive Power 750 VAr V.sub.Nom Nominal Voltage Magnitude 120 Vrms V.sub.min Minimum Voltage Magnitude 114 Vrms ω.sub.Nom Nominal AC Frequency 2π60  rad/sec |Δω.sub.max| Maximum Frequency Offset 2π0.5 rad/sec t.sub.rise VOC Rise Time    0.05 sec φ.sub.PM Power-Tracking Phase Margin  65 degrees

[0157] Based on the example AC and dynamic performance targets, a set of nominal VOC fundamental parameters may be determined by following the VOC design procedure detailed with respect to FIG. 3, above. These values will inform the design of the VOC-based power-tracking control described here. The nominal VOC fundamental parameters (determined by the method presented in FIG. 3) are shown in Table V below.

[0158] Next, the design of the linear compensators (e.g., 1155 and 1156 of FIG. 11), implementing G.sub.P(s) and G.sub.Q(s), respectively, may be determined based on the feedforward-linearized model in FIG. 12. The open-loop transfer function of the real-power-tracking loop can be written as

[00072]$\begin{array}{cc}{G}_{\mathrm{OL},P}\left(s\right)=-G_P\left(s\right).Math.m_P^{-1}.Math.\frac{{\omega }_c}{s+{\omega }_c}.Math.H_{\mathrm{FB},P}\left(s\right).& \left(99\right)\end{array}$

The compensator G.sub.P(s) is selected as a series Proportional-Integral controller of the form

[00073]$\begin{array}{cc}{G}_P\left(s\right)={\kappa }_{p,P}\left(1+{\kappa }_{i,P}.Math.\frac{1}{s}\right)={\kappa }_{p,P}.Math.\frac{s+{\kappa }_i,P}{s}.& \left(100\right)\end{array}$

where κ.sub.p,P and κ.sub.i,P are the proportional and integral gains Whose values are to be designed.

[0159] For brevity, the transfer function of the power feedback filter H.sub.FB,P is approximated as a delay equal to half the AC period T.sub.AC/2:

H.sub.FB,P(s)=H.sub.FB,Q(s)≈e.sup.−sT.sup.AC.sup./2.  (101)

This approximation is accurate for frequencies below half the AC frequency, which is the range of interest for power-tracking control.

[0160] Substituting (100) and (101) into (99) yields

[00074]$\begin{array}{cc}{G}_{\mathrm{OL},P}\left(s\right)=-{\kappa }_{p,P}.Math.m_P^{-1}.Math.\frac{{\omega }_c\left(s+{\kappa }_{i,P}\right)}{s\left(s+{\omega }_c\right)}.Math.e^{-{\mathrm{sT}}_{A.Math..Math.C}/2}.& \left(102\right)\end{array}$

[0161] The assignment κ.sub.i,P=ω.sub.C simplifies (102) to

[00075]$\begin{array}{cc}{G}_{\mathrm{OL},P}\left(s\right)=-k_{p,P}.Math.m_P^{-1}.Math.\frac{{\omega }_c}{s}.Math.e^{-{\mathrm{sT}}_{A.Math..Math.C}/2}.& \left(103\right)\end{array}$

[0162] To determine the proportional gain κ.sub.p,P, a target phase margin

[00076]${\phi }_{\mathrm{PM}}\in \left[\frac{\pi }{4},\pi \right]$

must first be selected. Evaluating (103) at s=jω and separating into magnitude and phase results in

[00077]$\begin{array}{cc}.Math.G_{\mathrm{OL},P}\left(j.Math..Math.\omega \right).Math.=.Math.k_{p,P}.Math.m_P^{-1}.Math.\frac{{\omega }_c}{\omega }.Math.,& \left(104\right)\\ \angle .Math..Math.G_{\mathrm{OL},P}\left(j.Math..Math.\omega \right)={\tan }^{-1}\left(\frac{\cos \left(\omega .Math..Math.T_{A.Math..Math.C}/2\right)}{\sin \left(\omega .Math..Math.T_{A.Math..Math.C}/2\right)}\right)=\left(\frac{\omega .Math..Math.T_{A.Math..Math.C}}{2}+\frac{\pi }{2}\right).& \left(105\right)\end{array}$

[0163] Setting (105) equal to −π+φ.sub.PM and solving for ω yields the target gain crossover frequency ω.sub.GC:

[00078]$\begin{array}{cc}{\omega }_{G.Math..Math.C}=\frac{\pi }{T}-\frac{2}{T}.Math.{\phi }_{\mathrm{PM}}.& \left(106\right)\end{array}$

[0164] Finally, (104) is evaluated at ω=ω.sub.GC and set equal to unity to find the value of κ.sub.p,P to yield the target phase margin:

[00079]$\begin{array}{cc}{\kappa }_{p,P}=.Math.m_P.Math.\frac{{\omega }_{\mathrm{GC}}}{{\omega }_c}.Math..& \left(107\right)\end{array}$

Applying the same procedure to the reactive power-tracking loop results in the value of κ.sub.p,Q that yields the same phase margin (assuming κ.sub.i,Q=κ.sub.i,P=ω.sub.c):

[00080]$\begin{array}{cc}{\kappa }_{p,Q}=.Math.m_Q.Math.\frac{{\omega }_{\mathrm{GC}}}{{\omega }_c}.Math..& \left(108\right)\end{array}$

[0165] The power-tracking parameters for the candidate design that result from the above example design procedure are also shown in Table V. These parameters were calculated based on the AC and dynamic performance targets shown in Table IV above.

TABLE-US-00005 TABLE V DESIGN VOC AND POWER-TRACKING PARAMETERS Symbol Description Value Units σ.sub.Nom VOC Conductance 10.7962 Ohms.sup.−1 α.sub.Nom VOC Cubic Coefficient 7.1975 A/V.sup.−3 L.sub.Nom VOC Inductance 39.899 μH C.sub.Nom VOC Capacitance 176.35 mF m.sub.P P-V Droop Slope −7.410e−3  V/W m.sub.Q Q-ω Droop Slope 3.780e−3 rad/sec/VAr ω.sub.c Cutoff Frequency 2π6.583 rad/sec ω.sub.GC Gain Crossover Frequency 2π8.33  rad/sec κ.sub.p, P P Proportional Gain 9.380e−3 V/W κ.sub.i, P P Integral Gain 2π6.583 rad/sec κ.sub.p, Q Q Proportional Gain 4.786e−3 rad/sec/VAr κ.sub.i, Q Q Integral Gain 2π6.583 rad/sec

[0166] As validation of the proposed VOC-based power-tracking designs disclosed herein, control unit 1110 of FIG. 11 was implemented in a digital microcontroller for a single-phase test inverter with the parameters presented in Table V. The test inverter was connected to a stiff 120 V AC voltage source (e.g., representing a mains utility connection) with a series impedance Z.sub.g=R.sub.g+jωL.sub.g. The value of Z.sub.g was selected so that it represents a 5% impedance on the basis of the inverter specifications with an X/R ratio of unity. This may represent a typical case in which an inverter is located at the end of a long underground distribution feeder line.

[0167] The test system was constructed in the Energy Systems Integration Facility (ESIF) at the National Renewable Energy Laboratory in Golden, Colo. The grid was operated at nominal voltage and frequency (V.sub.g=120 Vrms and ω.sub.eq=2π60 rad/sec) and the inverter (with VOC and power-tracking enabled) was synchronized to the grid and connected. A sequence of inverter power reference steps was then commanded, and the resulting power-tracking performance of the inverter recorded.

[0168] FIGS. 13A and 13B are graphical plots showing real and reactive power references and measured values for a prototype VOC-based power-tracking inverter, in accordance with one or more aspects of the present disclosure. In FIGS. 6A and 6B, solid lines 1302 and 1306 represent the actual real power (P) and reactive power (Q), respectively. Dotted lines 1304 and 1308 represent the reference real power (P.sub.ref) and reactive power (Q.sub.ref), respectively. Table VI shows the power reference steps and the resulting settling time, overshoot, and regulation error.

TABLE-US-00006 TABLE VI QUANTITATIVE EXPERIMENT RESULTS Worst Quantity Unit Step 1 Step 2 Step 3 Case Step time sec 1.0 1.5 2.0 N/A P.sub.ref W 750 −750 −500 N/A Q.sub.ref VAr 0 0 500 N/A P final W 749.99 −749.91 −499.94 N/A Q final VAr −0.013 0.015 499.99 N/A P regulation error % 1.09e−3 1.14e−3 0.812e−3 1.14e−3 Q regulation error % 1.17e−3 2.03e−3  1.49e−3 2.03e−3 P settling time sec 0.1248 0.203 0.1483 0.2031 Q settling time sec 0.1661 0.166 0.1142 0.1661 P overshoot % 10.013 2.220 1.1963 10.013 Q overshoot % 24.204 18.41 51.360 51.360

[0169] Observe in FIGS. 13A and 13B that when P.sub.ref and/or Q.sub.ref are stepped, P and Q both track their respective references with a slight overshoot and ringing. The maximum settling time is on the order of 200 msec, and the steady-state regulation error is near zero (see Table VI).

[0170] The overshoot of P is less than 10%, which is as expected from the design procedure (the selected φ.sub.PM=65 deg phase margin should yield an overshoot of about 10%). The overshoot of Q is as much as 51%. This may be a result of inverter terminal voltage rise due to the unmodeled impedances of the inverter filter and the network.

[0171] The response of Q experiences significant ringing in response to the step of P.sub.ref. This cross-coupling effect may be explained by inclusion of the effects of filter and network impedances in the model.

[0172] In one aspect, the techniques of the present disclosure provide enhanced VOC with power-tracking capability. The techniques described herein may adaptively adjust the parameters of the underlying VOC to regulate real and reactive power. The control behavior may be linearized by use of a feedforward linearization. A linearized model of the control dynamics is provided herein, along with a method for designing the control to meet specified AC and dynamic performance targets. Furthermore, experimental validation of one example control implementation was performed in a single-phase, grid-tied test inverter system.

[0173] The example control setup successfully regulated the inverter's real and reactive output power to track the specified references with near-zero steady-state error. Dynamic performance of the control prototype was similar to that predicted by the proposed model. The overshoot of reactive power was larger than expected. In addition, some cross-coupling was observed between real and reactive power. The control prototype successfully met target performance criteria with settling time on the order of 200 msec.

[0174] In some examples, the control techniques described herein may be modified to ensure that they maintain VOC's “grid-forming” behaviors such as synchronization and maintaining voltage and frequency within specified limits. In addition, the control techniques of the present disclosure may, in some examples, be used as the basis for the design of VOC-based control for specific types of distributed generation, such as photovoltaics or battery energy storage systems. This may result in a systemic solution to the challenges of power-electronic-dominated distribution feeders.

[0175] While described herein within the context of inverters (e.g., for PV system outputs), the techniques of the present disclosure may be utilized in various other contexts where virtual-oscillator-based power-tracking control may be beneficial. These additional contexts may include use of the techniques of the present disclosure in any power-electronic systems, such as those of power generation sources, power consumption sources, power storage sources, or any combination thereof.

[0176] The following examples may additionally or alternatively describe the techniques of the present disclosure

Example 1

[0177] A device comprising: a processor configured to implement a virtual oscillator circuit and output an oscillating waveform based on the virtual oscillator circuit; and power electronics operatively coupled to the processor and configured to convert, based on the oscillating waveform, direct current (DC) electricity to alternating current (AC) electricity, wherein the processor is further configured to: extract, from the virtual oscillator circuit, a virtual current based on an output current of the AC electricity, and output the oscillating waveform further based on an input voltage of the DC electricity.

Example 2

[0178] The device of example 1, wherein the virtual oscillator circuit comprises a virtual negative-conductance element.

Example 3

[0179] The device of example 2, wherein: the virtual oscillator circuit further comprises a virtual inductor and a virtual capacitor, and the virtual negative-conductance element, the virtual inductor, and the virtual capacitor are connected in parallel.

Example 4

[0180] The device of any of examples 1-3, wherein the virtual oscillator circuit comprises a virtual cubic voltage-dependent current source.

Example 5

[0181] The device of example 4, wherein: the virtual oscillator circuit further comprises a virtual inductor and a virtual capacitor, and the virtual cubic voltage-dependent current source, the virtual inductor, and the virtual capacitor are connected in parallel.

Example 6

[0182] The device of example 5, wherein the virtual cubic voltage-dependent current source consumes current proportional to a virtual voltage across the virtual capacitor.

Example 7

[0183] The device of any of examples 5-6, wherein the virtual oscillator circuit further comprises a virtual negative-conductance element that is connected in parallel with the virtual cubic voltage-dependent current source, the virtual inductor, and the virtual capacitor.

Example 8

[0184] The device of any of examples 1-7, wherein the processor is configured to output the oscillating waveform by: receiving, from the virtual oscillator circuit, a virtual output voltage; scaling the virtual output voltage based on the input voltage of the DC electricity to produce a scaled virtual output voltage; and comparing the scaled virtual output voltage to a carrier wave to produce the oscillating waveform.

Example 9

[0185] The device of any of examples 18, wherein: the virtual oscillator circuit comprises a virtual inductor and a virtual capacitor, and the processor is further configured to: compute a linear combination of a virtual voltage across the virtual capacitor and a virtual current through the virtual inductor to produce an alternate virtual output voltage; and output the oscillating waveform further based on the alternate virtual output voltage.

Example 10

[0186] The device of example 9, wherein the linear combination is computed as κ.sub.v.sub.c cos ω+κ.sub.vεi.sub.L sin φ wherein: κ.sub.v represents a voltage scaling factor, .sub.C represents the virtual voltage across the virtual capacitor, φ represents a rotation angle, ε represents a current scaling factor, and i.sub.L represents the virtual current through the virtual inductor.

Example 11

[0187] The device of any of examples 1-10, further comprising an output filter that includes a filter inductor and has a resistance, wherein: the virtual oscillator circuit comprises a virtual capacitor, and the processor is configured to output the oscillating waveform by: receiving, from the virtual oscillator circuit, a virtual output voltage; scaling the virtual output voltage based on a value of the resistance and an inductance value of the filter inductor to produce a first scaled virtual output voltage; combining the first scaled virtual output voltage with a virtual capacitor current flowing through the virtual capacitor to produce a combined virtual output voltage; and outputting the oscillating waveform further based on the combined virtual output voltage.

Example 12

[0188] The device of any of examples 1-11, wherein: the virtual oscillator circuit comprises a virtual inductor and a virtual cubic voltage-dependent current source, and the processor is further configured to implement a power tracking control unit that modifies at least one of a virtual inductance of the virtual inductor or a scaling factor of the virtual cubic voltage-dependent current source based on a specified output power.

Example 13

[0189] The device of example 12, wherein the power tracking control unit modifies the at least one of the virtual inductance or the scaling factor based further on a scaled virtual output voltage of the virtual oscillator circuit and the output current of the AC electricity.

Example 14

[0190] The device of any of examples 1-13, wherein the power electronics are configured to convert the DC electricity to AC electricity by alternatively engaging and disengaging power electronics switches based on the oscillating waveform.

Example 15

[0191] A device comprising: a processor configured to: implement a virtual oscillator circuit comprising a virtual capacitor, a virtual inductor, a virtual negative-conductance element, and a virtual cubic voltage dependent current source; generate, based on the virtual oscillator circuit, an oscillating waveform; and output, based on the oscillating waveform, at least one control signal to cause power electronics to convert direct current (DC) electricity to alternating current (AC) electricity; receive an indication of an output current of the AC electricity; extract, from the virtual oscillator circuit, a virtual current based on the output current of the AC electricity; and receive an indication of an input voltage of the DC electricity, wherein the oscillating waveform is generated further based on the input voltage of the DC electricity.

Example 16

[0192] The device of example 15, wherein the processor is configured to: compute a linear combination of a virtual voltage across the virtual capacitor and a virtual current through the virtual inductor to produce an alternate virtual output voltage, and generate the oscillating waveform based on the alternate virtual output voltage.

Example 17

[0193] The device of any of examples 15-16, wherein the processor is further configured to modify at least one of a virtual inductance of the virtual inductor or a scaling factor of the virtual cubic voltage-dependent current source based on a specified output power and at least one of: a scaled virtual output voltage of the virtual oscillator circuit or the output current of the AC electricity.

Example 18

[0194] A method comprising: determining, for a processor operatively coupled to power electronics configured to convert direct current (DC) electricity to alternating current (AC) electricity, based on a specified open-circuit voltage value for the power electronics, a voltage scaling factor value, wherein the processor is configured to implement a virtual oscillator circuit comprising a virtual capacitor, a virtual inductor, a virtual negative-conductance element and a virtual cubic voltage-dependent current source; configuring the processor to scale, using the voltage scaling factor value, a virtual oscillator output voltage of the virtual oscillator circuit thereby producing a scaled virtual oscillator output voltage for use in controlling the power electronics; determining, for the processor, based on a specified rated power output value for the power electronics and a specified rated-power voltage value for the power electronics that corresponds to the specified rated power output value, a current scaling factor value; configuring the processor to generate a virtual oscillator input current by scaling, using the current scaling factor value, an output current of the power electronics; determining, for the processor, based on the specified open-circuit voltage value for the power electronics and the specified rated-power voltage value for the power electronics, a conductance value; configuring the processor to implement the virtual negative-conductance element based on the conductance value; determining, for the processor, based on the conductance value, a cubic coefficient value; configuring the processor to implement the virtual cubic voltage-dependent current source based on the cubic coefficient value; determining, for the processor, a capacitance value based on at least one of: (i) a specified rated reactive-power output value for the power electronics, (ii) a specified maximum frequency offset value, (iii) the specified rated power output value for the power electronics, (iv) the specified open-circuit voltage value for the power electronics, (v) the specified rated-power voltage value for the power electronics, (vi) a specified rise time value, (vii) a specified nominal system frequency value, or (viii) a specified ratio of an amplitude of a third harmonic to an amplitude of a fundamental; configuring the processor to implement the virtual capacitor based on the capacitance value; determining, for the processor, based on the capacitance value and the specified nominal system frequency value, an inductance value; and configuring the processor to implement the virtual inductor based on the inductance value.

Example 19

[0195] The method of example 18, wherein: the specified rated power output value for the power electronics comprises a specified rated real power output value for the power electronics, and the specified rated-power voltage value for the power electronics comprises a specified rated real-power voltage value for the power electronics.

Example 20

[0196] The method of example 18, wherein: the processor is further configured to implement a transformation unit that computes a linear combination of a voltage across the virtual capacitor and a current through the virtual inductor to produce an alternate virtual output voltage, the specified rated power output value for the power electronics comprises a specified rated reactive power output value for the power electronics, and the specified rated-power voltage value for the power electronics comprises a specified rated reactive-power voltage value for the power electronics.

[0197] In one or more examples, the functions described may be implemented in hardware, software, firmware, or any combination thereof. If implemented in software, the functions may be stored on or transmitted over, as one or more instructions or code, a computer-readable medium and executed by a hardware-based processing unit. Computer-readable media may include computer-readable storage media, which corresponds to a tangible medium such as data storage media, or communication media, which includes any medium that facilitates transfer of a computer program from one place to another, e.g., according to a communication protocol. In this manner, computer-readable media generally may correspond to (1) tangible computer-readable storage media, which is non-transitory or (2) a communication medium such as a signal or carrier wave. Data storage media may be any available media that can be accessed by one or more computers or one or more processors to retrieve instructions, code and/or data structures for implementation of the techniques described in this disclosure. A computer program product may include a computer-readable storage medium.

[0198] By way of example, and not limitation, such computer-readable storage media can comprise RAM, ROM, EEPROM, CD-ROM or other optical disk storage, magnetic disk storage, or other magnetic storage devices, flash memory, or any other medium that can be used to store desired program code in the form of instructions or data structures and that can be accessed by a computer. Also, any connection is properly termed a computer-readable medium. For example, if instructions are transmitted from a website, server, or other remote source using a coaxial cable, fiber optic cable, twisted pair, digital subscriber line (DSL), or wireless technologies such as infrared, radio, and microwave, then the coaxial cable, fiber optic cable, twisted pair, DSL, or wireless technologies such as infrared, radio, and microwave are included in the definition of medium. It should be understood, however, that computer-readable storage media and data storage media do not include connections, carrier waves, signals, or other transient media, but are instead directed to non-transient, tangible storage media. Disk and disc, as used herein, includes compact disc (CD), laser disc, optical disc, digital versatile disc (DVD), floppy disk and Blu-ray disc, where disks usually reproduce data magnetically, while discs reproduce data optically with lasers. Combinations of the above should also be included within the scope of computer-readable media.

[0199] Instructions may be executed by one or more processors, such as one or more digital signal processors (DSPs), general purpose microprocessors, application specific integrated circuits (ASICs), field programmable logic arrays (FPGAs), or other equivalent integrated or discrete logic circuitry. Accordingly, the term “processor,” as used herein may refer to any of the foregoing structure or any other structure suitable for implementation of the techniques described herein. In addition, in some aspects, the functionality described herein may be provided within dedicated hardware and/or software modules. Also, the techniques could be fully implemented in one or more circuits or logic elements.

[0200] The techniques of this disclosure may be implemented in a wide variety of devices or apparatuses, including a wireless handset, an integrated circuit (IC) or a set of ICs (e.g., a chip set). Various components, modules, or units are described in this disclosure to emphasize functional aspects of devices configured to perform the disclosed techniques, but do not necessarily require realization by different hardware units. Rather, as described above, various units may be combined in a hardware unit or provided by a collection of inter-operative hardware units, including one or more processors as described above, in conjunction with suitable software and/or firmware.

[0201] The foregoing disclosure includes various examples set forth merely as illustration. The disclosed examples are not intended to be limiting. Modifications incorporating the spirit and substance of the described examples may occur to persons skilled in the art. These and other examples are within the scope of this disclosure and the following claims.