A Powerful Tool for Optimal Control of Energy Systems in Sustainable Buildings: Distortion Power Bivector

energies Article A Powerful Tool for Optimal Control of Energy Systems in Sustainable Buildings: Distortion Power Bivector Castilla Manuel V. 1,* and Martin Francisco 2   Citation: Manuel V., C.; Francisco, M. A Powerful Tool for Optimal Control of Energy Systems in Sustainable Buildings: Distortion Power Bivector. Energies 2021,14, 2177. https://doi.org/10.3390/en14082177 Academic Editor: Junemo Koo Received: 18 March 2021 Accepted: 8 April 2021 Published: 13 April 2021 Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/). 1Higher Polytechnic School, University of Seville, 41011 Sevilla, Spain 2Higher School of Industrial Engineering, University of Malaga, 29071 Málaga, Spain; [email protected] *Correspondence: [email protected] Abstract: In the field of building constructions, there is undeniably a growing need to optimize the energy systems which are a key target in new modern constructions and industrial buildings. In this sense, energy systems are being traced for the development of energy distribution networks that are increasingly smart, efficient, and sustainable. Modern generation and distribution energy systems, such as microgrids control systems, are being affected by the presence of linear and nonlinear loads, resulting a distorted voltage and current waveforms. Thus, it is stated that industrial and residential building heating and cooling loads behave essentially like sources of harmonics. This paper presents a new framework based on geometric algebra (GA) to the definition of a multivectorial distortion power concept, which is represented by a bivector that is geometrically interpreted to distinguish the rotated distortion and distortion power bivectors in these kinds of loads. Both bivectors, and their relations to the phase angles of distorted voltage are the main subject of this paper to interpret an optimal control of building energy. Numerical examples are used to illustrate of the suggested distortion power concept, as well as the information it provides for energy control in new buildings in a more sustainable way. Keywords: energy systems; building energy consumption; harmonics; geometric algebra; efficient constructions; heating and cooling loads 1. Introduction Nowadays, the electric power system has become one of the fundamental installations in construction. A new architecture is not only sustainable by the control of its energy, but also by the optimization of the installation and design itself. One of the most controversial issues in electrical engineering applications is the universal representation of the power equation of electrical systems in harmonic pollution situations. Moreover, the confusion extends to the concepts of reactive, not active and distortion powers. This is partly due to the proliferation of linear and non-linear loads connected to the circuits in the buildings. The ever-growing proliferation of power switching devices in modern industrial applications has increased the eventuality of unbalanced currents, unacceptable harmonic levels, and poor power factor in distribution systems. A clear example is seen in distributed generation systems, smart buildings or control systems, where many receivers are installed such as cycloconverters, speed drives, household appliances, battery power converters, power inverters and more. Thus, “modern home appliances are increasingly adopting power electronic frontends nowadays due to energy efficiency considerations. This trend has led to increased harmonic distortions in building power distribution systems and its supply transmission systems” [ 1 ]. Some works have been published to analyse this issue. For example, [ 2 ] investigates the current and voltage unbalance of nonlinear consumers in the residential and educational sectors. In [ 3 ], the objective is to present a systematic, and versatile technique to study the harmonic impact of residential load. Reference [ 4 ] presents a study of the use of active filters placed at Energies 2021,14, 2177. https://doi.org/10.3390/en14082177 https://www.mdpi.com/journal/energies Energies 2021,14, 2177 2 of 17 different locations for harmonic mitigation methods in wind plants. Reference [ 5 ] proposed the harmonic compensation method on a residential distribution system and in [ 6 ], an advanced model is used to retrieve data from publications related to energy quality. The result of using such a high number of non-linear loads in buildings is that the current waveform is distorted, causing excessive harmonic voltages to be generated. This situation feeds back and causes a progressive degradation to the power quality of the supply [ 7 ]. In [ 8 ] a new indicator designated called a power quality factor (PQF) is suggested to integrally reflect the power transfer quality of a general three phase network feeding non-linear loads. Modern energy generation and distribution systems are a great advance for today’s society, since they allow better management and integration of energy with new sources such as renewable energy, wind energy, solar energy, nuclear energy etc. These systems, such as multi-energy micro grids (MEMGs) are composed by distributed generators with adequate placement for the system energy efficiency [9]. Also, the proximity of many of these buildings with similar industrial activities will contribute to the distortion of the electric power quality of feed supplying these constructions. These harmonics can cause serious problems in the power system of industrial architecture, for example excessive heat in industrial appliances, premature aging of electrical equipment, de-rating of equipment, fault of protection and lower power factor or inefficiencies in the transmission of energy. In different countries buildings are responsible for the 21–25% of total energy consumption and the amount of this energy used for heating and cooling systems is about 55% in the residential sector [ 10 ]. It is a situation that causes a progressive degradation to the power quality of the supply, [ 11 ] give a description of state of the art for the distributed power generation systems (DPGS) based on renewable energy; in [ 12 ] an extensive survey of loads is conducted in various residential buildings and [ 13 ] aimed to analyse the individual influences, which aggravate the power quality of the distribution grid. In this sense, other lines of research related to heating and cooling loads of industrial and residential buildings analyse the phase change materials as cold thermal energy storage media for district cooling [14]. These harmful and costly effects of harmonics have been discussed extensively and spurred stringent requirements by international institutions regarding the allowed levels of harmonics at the point of connection to the power supply [ 15 ], (IEEE Recommended Practice and Requirements for Harmonic Control in Electric Power Systems). On the other hand, as is well known, a phase displacement between corresponding voltages and currents [ 16 ] indicates both a low utilization of the generation and distribution equipment and increased line losses for the same power consumption level. Significant research efforts have been focussed to the analysis of this area, especially definition and interpretation of reactive and distortion powers and the challenge of developing a power theory suitable for harmonic analysis has continued forever last century In sinusoidal situations, the mathematical tools used have been based on Steinmetz theory [ 17 ]. In non-sinusoidal and linear conditions, Budeanu [ 18 ] decomposes the apparent power in active, reactive and distortion components and Shepherd et al. [ 19 ] suggest a new concept of reactive power. Likewise, in [ 20 ] a reactive power was defined that permits the derivation of the value of passive linear devices to achieve a relative optimum power factor and in [21], Fryze published his power equation in time domain. Other formulations developed in frequency domain [ 22 – 26 ], affect the search for equations that interpret certain physical phenomena or analyse the distortion component in non-sinusoidal regimes. Among them [ 25 ], provides a meticulous discussion of the distortion power concept and [ 26 ] based on the Fourier series, assumed that the distortion power may consist of active and reactive components. Unfortunately, the mentioned work [ 26 ], needs to know the load parameters to define these components. Moreover, it is not possible to generalize this approach in vector form to any load model, because the basic quantities, voltage and current, are not ordered within a suitable mathematical structure. Energies 2021,14, 2177 3 of 17 Foundations: Literature Review In sinusoidal conditions, Steinmetz established his theory of power, which with slight modifications remains in force to-day. However, it is essential to know the electrical energy balances in any microgrid to be able to make the right decisions for its control and optimization. In these conditions, the apparent power is defined by S=P+jQ , where P is the active power, Qis the reactive power and jis the imaginary unit. In non-sinusoidal conditions, Budeanu and Fryze formulated the first power equations for n-sinusoidal operation between 1920 and 1930. The first author [ 18 ], addressed the question in the frequency domain. He defined the distortion power as: D2 B=SB2−P2+Q2 B(1) where P=∑ nUnIncos φn , QB=∑ nUnInsin φn and SB are the active, reactive and apparent powers, respectively. The second author [ 21 ], developed the current decomposition into active ( ia ), and reactive (ir) currents in the time domain. The active current is given by: ia=P ku(t)k2u(t)=Geu(t)(2) where Ge is associated to an equivalent conductance of the arbitrary linear load. The residual part of the total current, is named reactive or non-active current: ir=i−ia(3) The rms values of the currents fulfil the equality: kirk2=kik2−kiak2(4) and the scalar product of these currents is given by: hia,iri=1 TZT 0iairdt =1 TZT 0ia(i−ia)dt =0 (5) If (4) is multiplied by ku(t)k2it is seen that: Qr2=S2−P2(6) where Qr6=QB,P,S, are the nonactive, active and apparent powers, respectively. Motivated by these theories, numerous valuable works have appeared as mentioned above, but not all of them consider the multivectorial character of apparent power components. In this context, in reference [ 27 ] the authors use geometric algebra to define a multivector power based on decomposition of the instantaneous current into the active and reactive components. It should be noted that this approach does not distinguish between reactive and distortion power from a mathematical viewpoint. They state that the scalar part of the apparent power represents the average power (active power). Furthermore, the bivectorial part of the apparent power contains power components involving like-frequency and cross-frequency products which are not suitable for reactive and/or distortion compensator design. Moreover, a serious limitation of this theory is based on the difficulty to consider voltage phase angles. This restriction will probably lead to erroneous results. In particular [ 28 ] is an important and original contribution to the power analysis in non-sinusoidal networks. This theory was extended by the author for linear circuits in n-sinusoidal conditions [ 29 ]. In [ 30 ], Montoya et al. have analysed power flows under n-sinusoidal conditions via GA a as new advance in power theory. The large number of papers published on the physical and/or mathematical nature of distortion power and their applications, suggest that work in this field is not finished. For Energies 2021,14, 2177 4 of 17 this reason, in this paper a complete solution to the power theory analysis and particularly the distortion power concept for single-phase periodic n-sinusoidal linear and nonlinear operation and its relation to the optimization of energy control is presented. This power concept is based on a new frequency-domain approach, which uses geometric algebra [31–34], to define and calculate any component of the distortion power as a bivector from a generalized complex geometric algebra (GCGA). Thus, the distortion power bivector incorporates in a single equation all required information (magnitude, direction, and sense) where the experimental phenomena are implicit. In this sense, the geometrical interpretation of the distortion power bivector is quite effective to clarify its nature from a mathematical and physical nature viewpoint and optimize energy flow in distributed generation, control systems or intelligent buildings. 2. Mathematical Preliminaries Geometric algebra is a mathematical structure developed over the last 40 years based on Clifford Algebras. A Geometric algebra can be defined simply by specifying appropriate rules for multiplying vectors. Thus, let Vn an n-dimensional linear space over the real numbers. The geometric product of vectors a⊗b if a , b∈ Vn , [ 28 ], can be decomposed into a symmetric inner product: a·b=1 2(a⊗b+b⊗a)(7) and an antisymmetric outer product: a∧b=1 2(a⊗b−b⊗a)(8) Therefore, ab has the canonical decomposition: a⊗b=a·b+a∧b(9) The inner product a·b is a scalar and the outer product a∧b is called bivector (or 2-vector). Geometrically, it represents a directed plane, in much the same way as a vector represents a directed line segment (Figure 1). Energies 2021, 14, 2177 4 of 17 For this reason, in this paper a complete solution to the power theory analysis and particularly the distortion power concept for single-phase periodic n-sinusoidal linear and nonlinear operation and its relation to the optimization of energy control is presented. This power concept is based on a new frequency-domain approach, which uses geometric algebra [31–34], to define and calculate any component of the distortion power as a bivector from a generalized complex geometric algebra (GCGA). Thus, the distortion power bivector incorporates in a single equation all required information (magnitude, direction, and sense) where the experimental phenomena are implicit. In this sense, the geometrical interpretation of the distortion power bivector is quite effective to clarify its nature from a mathematical and physical nature viewpoint and optimize energy flow in distributed generation, control systems or intelligent buildings. 2. Mathematical Preliminaries Geometric algebra is a mathematical structure developed over the last 40 years based on Clifford Algebras. A Geometric algebra can be defined simply by specifying appropriate rules for multiplying vectors. Thus, let 𝒱 an n-dimensional linear space over the real numbers. The geometric product of vectors 𝑎⊗𝑏 if 𝑎,𝑏∈𝒱, [28], can be decomposed into a symmetric inner product: 𝑎⋅𝑏= 12(𝑎⊗𝑏+𝑏⊗𝑎) (7) and an antisymmetric outer product: 𝑎∧𝑏=12(𝑎⊗𝑏−𝑏⊗𝑎) (8) Therefore, 𝑎𝑏 has the canonical decomposition: 𝑎⊗𝑏=𝑎⋅𝑏+𝑎∧𝑏 (9) The inner product 𝑎⋅𝑏 is a scalar and the outer product 𝑎∧𝑏 is called bivector (or 2-vector). Geometrically, it represents a directed plane, in much the same way as a vector represents a directed line segment (Figure 1). Figure 1. Vector, bivector and trivector representation. We can regard 𝑎∧𝑏 as a directed area with a norm ‖𝑎∧𝑏‖ equal to the usual scalar area of each parallelogram in Figure 1, with the direction of the plane in which the parallelogram lies, and sense that can be assigned to the parallelogram in the plane. Then, just as a vector 𝑎 represents (or is represented by) a directed line segment and a bivector 𝑎∧ Figure 1. Vector, bivector and trivector representation. We can regard a∧b as a directed area with a norm ka∧bk equal to the usual scalar area of each parallelogram in Figure 1, with the direction of the plane in which the parallelogram lies, and sense that can be assigned to the parallelogram in the plane. Then, just as a vector a represents (or is represented by) a directed line segment and a bivector a∧b represents Energies 2021,14, 2177 5 of 17 a directed plane segment, the trivector (3-vector) a∧b∧c represents a directed space segment (the parallelepiped with edges a, b, c). However the electrical quantities voltage and current do not have an easy interpretation in classic Clifford algebra. For a clear representation of these waveforms, a new geometric algebra has been constructed—a generalization of the classic Clifford algebra [ 32 ] which we have termed “generalized complex geometric algebra” (GCGA) and denoted by {CGn,} . In this framework, C is the complex vector space, Gn is the Clifford algebra or geometric algebra associated to the n-dimensional real space Vnand: :(<◦⊗)(10) is the new generalized geometric product. A detailed description of this new structure (GCGA), of the geometric product “  ” and their properties, are given in references [ 32 , 33 ]. 3. Geometric Phasor Representation of Periodic Voltage and Current Waveforms We start by considering periodic voltage and current waveforms of an arbitrary singlephase system. Both signals can be considered elements of the vector space with inner and outer product. Suppose a non-sinusoidal voltage u(t)=√2∑p∈L∪NUpsinpωt+αp where pis the harmonic order of u(t) and a current i(t)=√2∑q∈N∪MIqsinqωt+βq where qis the harmonic order of i(t) . It is assumed that a group of voltage harmonics N exist that have corresponding current harmonics of the same frequencies, that components L of the supply voltage exist without corresponding currents, and that components M of current exist without corresponding voltages. In linear operation, βq=αq−φq , φq is the harmonic impedance phase angle and L={ϕ} , M={ϕ} . The capital Up and Iq represent rms values of up(t) and iq(t) . In the {CGn,} structure spanned by an orthonormal basis {σ1,σ2,σ3, . . . σn} , the associated p-th harmonic voltage and q-th harmonic current can be represented by the geometric-phasors: e Up=Upejαpσp(11) eIq=Ipejβqσq(12) where Up=ke Upk,Iq=keIqk. Then: e U=∑p∈L∪NUpejαpσp(13) eI=∑q∈N∪MIqej(αq−φq)σqq∈N+Iqejβqq∈M(14) In particular, on linear operation the Equation (14) is given by eI=∑q∈NIqej(αq−φq)σq geometric -phasor, where eIqharmonic current can be decomposed as follows: eIq=Iqejαqe−jφqσq=eIqe−jφqσq=eIqcos φqσq−jeIqsinφqσq=eIqk−jeIq⊥(15) The subscripts “ k ” and “ ⊥ ” indicates “in phase” and “in quadrature” respectively. The eI∗ qconjugate harmonic current is given by: eI∗ q=eI∗ qk+jeI∗ q⊥(16) where (∗)is the standard “conjugate complex”operation. Energies 2021,14, 2177 6 of 17 4. Distortion Power Bivector 4.1. Power Multivector Consider an arbitrary nonlinear single-phase system supplied by the voltage u(t). The resulting current has an instantaneous value given by i(t). The instantaneous voltage and current signals have been defined in the section above. In order to represent the concept of power multivector in GCGA, we introduce a new rule for multiplying geometric-phasors (complex-vectors). This rule is the new geometric product (10), where < [ 32 ] is an application over complex planes and “ ⊗ ” is the standard geometric product. Apparent power can thus be expressed as a multivector e S in {CGn} , generated by “” of the voltage and conjugate current geometric phasors: e S=  e U·eI∗ |{z} e Ω• +e U∧eI∗ |{z} e Ω∧   (17) and which obeys the usual conservation law [ 33 ]. From Equation (17), one easily defines that power multivector is given by: e S=∑ p∈N∪L q∈N∪Me UpeI∗ q ∑ p=qUpIpcos φp+j∑ p=qUpIpsin φp!σ0+∑ p<q p,q∈Nnej(αp−αq)UpIqejφq−UqIpejφpσpqo | {z } e Dpq,Lin +∑ p∈L∪N,q∈M p∈L,q∈N ejαpUpIqe−jβqσpq | {z } e Dpq,Nonlin =e P+je Q+e ∆ (18) which consist of a complex-scalar part e Ω•=e P+je Q and a complex-bivector part e Ω∧=e ∆ . In Equation (18), “  ” is the new “generalized complex geometric product” [ 32 ]. The components e P,e Q,e ∆in Equation (18), can be represented as: e P=Rehe U·eI∗i=∑ p∈N UpIpcosφpσ0(19) e Q=Imhe U·eI∗i=∑ p∈N UpIpsinφpσ0(20) e ∆=∑ p<q p,q∈N <e Up∧eI∗ qσpq +∑ p∈L∪N,q∈M p∈L,q∈N ejαpUpIqe−jβq(21) The generic term e Dpq (18) is a bivector named distortion power term and it is related with the generic term e ∆pq named rotated distortion power term: e ∆pq =ej(αp−αq)e Dpq ⇒ ke ∆pqk=ke Dpqk(22) The e ∆pq terms in Equations (18) and (21), can be collected into four groups and written in a short notation related with the sets L, N and Mas: e ∆N,N=∑ p,q∈Ne ∆pq ⇒e ∆pq =ej(αp−αq)e Dpq,Lin (23) e ∆L,N=∑ p∈L q∈Ne ∆pq ⇒e ∆pq =ej(αp−αq)e Dpq,Nonlin (24) Energies 2021,14, 2177 7 of 17 e ∆L,M=∑ p∈L q∈Me ∆pq ⇒e ∆pq =ejαpe Dpq,Nonlin (25) e ∆N,M=∑ p∈N q∈Me ∆pq ⇒e ∆pq =ejαpe Dpq,Nonlin (26) Consequently, the complex-bivector part, Equation (21), named rotated distortion power, is given by: e ∆=e ∆N,N+e ∆L,N+e ∆L,M+e ∆N,M(27) Thus, in the cross products (21), each e ∆pq term represents the interaction of voltage and current harmonics, with different frequencies. This term is associated to e Dpq linear and/or e Dpq nonlinear distortion power bivector term. More exactly, in linear operation the generic e ∆pq term represents to a e Dpq term (distortion power bivector) associated to a rotation defined by αp−αq angle. In non-linear operation e Dpq term (distortion power bivector) is associated to a rotation αp−αqand/or αpangle. The new quantity, rotated distortion power proposed in Equation (21), have three basic attributes: magnitude, direction and sense. In the numerous situations direction and sense are not required; then the magnitude ke ∆k=ke Dk was established to manage a concrete problem. However, there are situations where direction and sense are necessary, as those that depend on the origin and nature of distortion. In these cases, the distortion power bivector e ∆ is the appropriate quantity to solve the problem, since it incorporates in a succinct expression all required information. The convenience of one or the other is dictated by the necessities of the situation. In this sense, ke ∆k=ke Dk is a consequence of e ∆ and both, magnitude and bivector are complementary. It is seen from Equation (27), that due to the independence of the basis σ1...n , the following equality is true: ke ∆k2=ke Dk2=ke ∆N,Nk2+ke ∆L,Nk2+ke ∆L,Mk2+ke ∆N,Mk2(28) 4.2. Distortion Power Bivector: Admittance Matrix Representation When the load impedance is linear, the only current harmonics that will flow must correspond to the supply voltage harmonics. Thus, in linear operation, the harmonic indexes of voltage and current, p,q∈Ncan define the quantities: ˆ Y=Yppn p=1=Yp⇒Ypp =Gp−jBp(29) Gp=Ip Up cosφp,Bp=Ip Up sinφp(30) where ˆ Y is called the n-frequency admittance matrix, Gp is the p-th conductance and Bp is the p-th susceptance. The multivectorial Ohm’s law is now given by: eI=e Uˆ Y=∑ qeIq⇒eI∗=e U∗ˆ Y∗=∑ qeI∗ q(31) where: eIq=e UqGq−jBqσq=eIk−jeI⊥(32) eI∗ q=e U∗ qGq+jBqσq=eI∗ k+jeI∗ ⊥(33) that coincide with Equations (15) and (16), respectively. In this case, e S multivector may be written in the form: e S=e UeI∗=e UeI∗ k+je UeI∗ ⊥(34) Energies 2021,14, 2177 8 of 17 Now it is seen from Equations (A2) and (34) that: e UeI∗ k=e U·eI∗ k+e U∧eI∗ k(35) and the following magnitude expressions of the active and active rotated distortion power components are obtained: e P=e U·eI∗ k=∑ Pe Upe U∗ pGp⇒P=∑ P U2 pGp(36) e ∆Act =e U∧eI∗ k=∑p<q p,q∈Ne Up∧eI∗ qkke ∆Actk2=∑p<qUp2Iq2Gp−Gq2(37) where from Appendix A, Equation (A4): e Up∧eI∗ qk=Upσp∧U∗ qGqσq+Uqσq∧U∗ pGpσp=ej(αp−αq)e Dpq,Act =e ∆pq,Act (38) On the other hand, from Equations (A5) and (34) can be obtained that: e UeI∗ ⊥=e U·eI∗ ⊥+e U∧eI∗ ⊥(39) Similarly, another two magnitude expressions for reactive and reactive rotated distortion powers are obtained: e Q=e U·eI∗ ⊥=∑pUpU∗ pBp⇒Q=∑pU2 pBp(40) e ∆Re act =e U∧eI∗ ⊥=∑p<q p,q∈Ne Up∧eI∗ q⊥ke ∆Re actk2=∑p<qUp2Iq2Bp−Bq2(41) where from Appendix A, Equation (A7): e Up∧eI∗ q⊥=Upσp∧U∗ qBqσq+Uqσq∧U∗ pBpσp=ej(αp−αq)e Dpq,Re act =e ∆pq,Re act (42) Combination of Equations (36)–(40) gives the apparent power multivector: e S=e UeI∗=e UeI∗ k+je UeI∗ ⊥=(P+jQ)σ0+∑p<q p,q∈N ej(αp−αq)(e Dpq,Act +je Dpq,Re act) = e P+je Q+e ∆(43) Note that the magnitudes of Equations (36)–(41) agree with the ones defined by a linear operation in [ 23 , 27 ]. These magnitudes are not simply calculated from algebraic operations, but rather they are a result of the multivectorial analysis. It should be noted that the introduction of a multivector representation became possible by employing the GCGA. In general, it can be said that today’s accepted approaches cannot explain the results obtained here. 4.3. Geometric Interpretation of Distortion Power Bivectors In the (CGn,) algebraic structure, the bivectorial part e Ω∧ of Equations (18) and (21), can be represented as a linear combination of basis bivectors, each encoding a distinct independent pq-plane involving cross-frequency products ( p6=q ) in n-dimensional space. Thus, Figure 2shows the geometric interpretation into generic pq-planes of e ∆pq rotated distortion power and the corresponding e Dpq distortion power term when p , q∈N , q6=p and/or p∈L , q∈N . Therefore, it is seen that Figure 3shows, into generic pq-planes, the geometric interpretation of e ∆pq rotated distortion power and the corresponding e Dpq distortion power term when p∈L∪N , q∈M . It can be seen from Figures 2and 3, that the geometric interpretation of Equation (21) is based on the association of complex planes (C) to each multivectorial element of Clifford basis Gn. Energies 2021,14, 2177 9 of 17 Energies 2021, 14, 2177 9 of 17 Figure 2. 𝜎 𝑝lane. Representation of the 𝛥󰆻and 𝐷  distortion bivector term in linear operation (𝑝,𝑞∈𝑁,𝑞≠𝑝)and/or non-linear operation, (𝑝∈𝐿, 𝑞∈𝑁) . Figure 3. 𝜎 𝑝lane. Representation of the 𝛥󰆻and 𝐷  distortion bivector term in non-linear operation (𝑝∈𝐿∪𝑁, 𝑞∈𝑀). Note in Figure 2 that: 𝑅𝑒𝛥󰆻=𝑅𝑒𝑒()𝐷 ≠𝛥 󰆻, (44) 𝐼𝑚𝛥󰆻=𝐼𝑚𝑒()𝐷 ≠𝛥 󰆻, (45) Similarly, in Figure 3: 𝑅𝑒𝛥󰆻=𝑅𝑒𝑒𝐷 ≠𝛥 󰆻, (46) 𝐼𝑚𝛥󰆻=𝐼𝑚𝑒𝐷 ≠𝛥 󰆻, (47) It is certain that the value (magnitude),𝛥󰆻=𝐷 , of the distortion power does not provide information to the waveform distortion, but each 𝛥󰆻 distortion power bivector term, has a direct dependence of the harmonic voltage phase angles. Notice that the choice Figure 2. σpq p lane. Representation of the e ∆pq and e Dpq distortion bivector term in linear operation (p,q∈N,q6=p)and/or non-linear operation, (p∈L,q∈N). Energies 2021, 14, 2177 9 of 17 Figure 2. 𝜎 𝑝lane. Representation of the 𝛥󰆻and 𝐷  distortion bivector term in linear operation (𝑝,𝑞∈𝑁,𝑞≠𝑝)and/or non-linear operation, (𝑝∈𝐿, 𝑞∈𝑁) . Figure 3. 𝜎 𝑝lane. Representation of the 𝛥󰆻and 𝐷  distortion bivector term in non-linear operation (𝑝∈𝐿∪𝑁, 𝑞∈𝑀). Note in Figure 2 that: 𝑅𝑒𝛥󰆻=𝑅𝑒𝑒()𝐷 ≠𝛥 󰆻, (44) 𝐼𝑚𝛥󰆻=𝐼𝑚𝑒()𝐷 ≠𝛥 󰆻, (45) Similarly, in Figure 3: 𝑅𝑒𝛥󰆻=𝑅𝑒𝑒𝐷 ≠𝛥 󰆻, (46) 𝐼𝑚𝛥󰆻=𝐼𝑚𝑒𝐷 ≠𝛥 󰆻, (47) It is certain that the value (magnitude),𝛥󰆻=𝐷 , of the distortion power does not provide information to the waveform distortion, but each 𝛥󰆻 distortion power bivector term, has a direct dependence of the harmonic voltage phase angles. Notice that the choice Figure 3. σpq p lane. Representation of the e ∆pq and e Dpq distortion bivector term in non-linear operation (p∈L∪N,q∈M). Note in Figure 2that: Rene ∆pqo=Renej(αp−αq)e Dpqo6=e ∆pq,Act (44) Imne ∆pqo=Imnej(αp−αq)e Dpqo6=e ∆pq,Re act (45) Similarly, in Figure 3: Rene ∆pqo=Renejαpe Dpqo6=e ∆pq,Act (46) Imne ∆pqo=Imnejαpe Dpqo6=e ∆pq,Re act (47) Energies 2021,14, 2177 16 of 17 e ∆L,N= ∑ p∈Le Up! ∑ q∈NeI∗ q!=∑ p∈L,q∈N ej(αp−αq)UpIqejφqσpq | {z } e Dpq,Nonlin                e ∆L,M Act =∑ p,qejαpRe UpIqejβqσpq | {z } e Dpq,Act e ∆L,M Re act =∑ p,qejαpImUpIqejβqσpq | {z } e Dpq,Re act (A13) e ∆N,M= ∑ p∈Ne Up! ∑ q∈MeI∗ q!=∑ p,qejαpUpIqejβqσpq | {z } e Dpq,Nonlin                e ∆N,M Act =∑ p,qejαpRe UpIqejβqσpq | {z } e Dpq,Act e ∆N,M Re act =∑ p,qejαpImUpIqejβqσpq | {z } e Dpq,Re act (A14) References 1. Wang, Y.; Yong, J.; Sun, Y.; Xu, W.; Wong, D. Characteristics of Harmonic Distortion in Residential Distribution Systems. IEEE Trans. Power Deliv. 2017,32, 1495–1504. [CrossRef] 2. Popa, G.N.; Iagär, A.; Dinis, C.M. Considerations on Current and tage Unbalance of Nonlinear Loads in Residential and Educational Sectors. Energies 2021,14, 102. [CrossRef] 3. Salles, D.; Jiang, C.; Xu, W.; Freitas, W.; Mazin, H.E. Assessing the collective harmonic impact of modern residential loads—Part I: Methodology, IEEE Trans.Power Deliv. 2012,27, 1937–1943. [CrossRef] 4. Schwanz, D.; Bollen, M.; Larsson, A.; Kocewiak, Ł.H. Harmonic mitigation in wind power plants: Active filter solutions. In Proceedings of the IEEE 2016 17th International Conference on the Harmonics and Quality of Power (ICHQP), Belo Horizonte, Brazil, 16–19 October 2016; pp. 220–225. 5. Munir, S.; Li, Y.W. Residential distribution system harmonic compensation using PV interfacing inverter. IEEE Trans. Smart Grid 2013,4, 816–827. [CrossRef] 6. Montoya, F.G.; Baños, R.; Alcayde, A.; Montoya Maria, G.; Manzano-Agugliaru, F. Power Quality: Scientific Collaboration Networks and Research Trends. Energies 2018,11, 2067. [CrossRef] 7. Bollen, M.H.; Das, R.; Djokic, S.; Ciufo, P.; Meyer, J.; Rönnberg, S.K.; Zavodam, F. Power quality concerns in implementing smart distribution-grid applications. IEEE Trans. Smart Grid 2017,8, 391–399. [CrossRef] 8. Sharon, D.; Montaño, J.C.; Castilla, M.; Borrás, D.; López, A.; Gutierrez, J. Power Quality Factor for Networks Supplying Unbalanced Nonlinear Loads. IEEE Trans. Instrum. Meas. 2008,57, 1268–1274. [CrossRef] 9. Li, Z.; Xu, Y.; Fang, S.; Mazzoni, S. Optimal placement of hetherogeneous distributed generators in a grid-connected multi-energy microgrid under uncertainties. IET Renew. Power Gener. 2019,13, 2623–2633. [CrossRef] 10. IEA International Energy Agency. Energy Efficiency Indicators Highlights; OECD/IEA: Paris, France, 2016. 11. Bouzid, A.M.; Guerrero, J.M.; Cheriti, A.; Bouhamida, M.; Sicard, P.; Benghanem, M. A survey on control of electric power distributed generation systems for microgrid applications. Renew. Sustain. Energy Rev. 2015,44, 751–766. [CrossRef] 12. Cherian, E.; Bindu, G.; Nair, P.C. Pollution impact of residential loads on distribution system and prospects of DC distribution. Eng. Sci. Technol. Int. J. 2016,19, 1655–1660. [CrossRef] 13. Otcenasova, A.; Bolf, A.; Altus, J.; Regula, M. The Influence of Power Quality Indices on Active Power Losses in a Local Distribution Grid. Energies 2019,12, 1389. [CrossRef] 14. Mazzoni, S.; Yin Sze, J.; Nastasi, B.; Ooi, S.; Desideri, U.; Romagnoli, A. A techo-economic assessment on the adoption of latent heat thermal energy storage systems for district cooling optimal dispatch & operations. Appl. Energy 2021,289, 116646. 15. IEE Standars Association. 519-2014—IEEE Recommended Practice and Requirements for Harmonic Control in Electric Power Systems; IEEE Power and Energy Society: Piscataway, NJ, USA, 2014; pp. 1–17. [CrossRef] 16. Cheng, Q.; Wang, C.; Wang, J. Analysis on Displacement Angle of Phase-Shifted Carrie PWM for Modular Multilevel Converter. Energies 2020,13, 6743. [CrossRef] 17. Steinmetz, C.P. Theory and Calculation of Alternating Current Phenomena; McGraw-Hill Book Company, Incorporated: New York, NY, USA, 1916. 18. Budeanu, C.I. Puisances Reactives et Fictives; Instytut Romain de l´Energie: Bucharest, Romania, 1927. Energies 2021,14, 2177 17 of 17 19. Shepherd, W.; Zhakikhani, P. Suggested definition of reactive power for nonsinusoidal systems. Proc. Inst. Elect. Eng. 1972 ,119, 1361–1362. [CrossRef] 20. Sharon, D. Reactive power definitions and power factor improvement in nonlinear systems. Proc. IEE 1973 ,120, 704–706. [CrossRef] 21. Fryze, S. Wik-,blind-, unscheinleistung in elektrischens tromkreisen mit nicht-sinusoidalen verlauf von stromund spanning. Elekt. Z. 1932,53, 596–599, 625–627, 700–702. 22. Czarnecki, L.S. Considerations on the reactive power in non-sinusoidal situations. IEEE Trans. Instr. Meas. 1985 ,IM-34, 399–404. [CrossRef] 23. Ghassemi, F. New Apparent Power and Power Factor with Non-Sinusoidal waveforms. In Proceedings of the Power Engineering Society Winter Meeting, Singapore, 23–27 January 2000; pp. 2852–2857. [CrossRef] 24. Sommariva, A.M. Power Analysis of One-Ports Under Periodic Multi-Sinusoidal Operation. IEEE Trans. Circuits Syst. I Regul. Pap. 2006,53, 2068–2074. [CrossRef] 25. Czarnecki, L.S. Distortion power in systems with nonsinusoidal tage. IEE Proc. B 1992,139, 276–280. 26. Slonim, M.A. Distortion power in linear and non-linear systems Int.J. Electron. 1990,68, 769–778. [CrossRef] 27. Menti, A.; Zacharias, T.; Milias-Argitis, J. Geometric Algebra: A powerful tool for representing power under nonsinusoidal conditions. IEEE Trans. Circuits Syst. I Reg. Papers 2007,54, 601–609. [CrossRef] 28. Castro-Núñez, M. The Use of Geometric Algebra in the Analysis of Non-sinusoidal Networks and the Construction of a Unified Power Theory for Single Phase Systems-A Paradigm Shift. Ph.D. Thesis, University of Calgary, Calgary, AB, USA, 2013. 29. Castro-Nuñez, M.; Castro-Puche, R. Advantages of geometric algebra over complex numbers in the analysis of networks with nonsinusoidal sources and linear loads. IEEE Trans. Circuits Syst. I Regular Pap. 2012,59, 2056–2064. [CrossRef] 30. Montoya Francisco, G.; Baños, R.; Alcayde, A.; Arrabal-Campos, F.M. A new approach to single phase systems under sinusoidal and non-sinusoidal supply using geometric algebra. Electric Power Syst. Res. 2020,189, 106605. [CrossRef] 31. Doran, C.; Lasenby, A. Geometric Algebra for Physicists; Cambridge University Press: Cambridge, UK, 2005. 32. Castilla, M.; Bravo, J.C.; Ordoñez, M.; Montaño, J.C. Clifford Theory: A Geometrical Interpretation of Multivectorial Apparent Power. IEEE Transactions Circuits Syst. Regular Pap. 2008,55, 3359–3367. [CrossRef] 33. Castilla, M.; Bravo, J.C.; Ordoñez, M. Geometric Algebra: A Multivectorial Prof. of Tellegen’s Theorem in Multiterminal networks. IET Circuits Devices Syst. 2008,2, 383–390. [CrossRef] 34. Bravo, J.C.; Castilla, M.V. Energy Conservation Law in Industrial Architecture: An Approach through Geometric Álgebra. Symmetry 2016,8, 92. [CrossRef]