Energy conservation law in industrial architecture: an approach through geometric algebra

symmetry S S Article Energy Conservation Law in Industrial Architecture: An Approach through Geometric Algebra Juan C. Bravo and Manuel V. Castilla * Higher Polytechnic School, University of Seville, Seville 41011, Spain; [email protected] *Correspondence: [email protected]; Tel.: +34-661-753-960 Academic Editor: Ka Lok Man Received: 5 May 2016; Accepted: 29 August 2016; Published: 6 September 2016 Abstract: Since 1892, the electrical engineering scientific community has been seeking a power theory for interpreting the power flow within electric networks under non-sinusoidal conditions. Although many power theories have been proposed regarding non-sinusoidal operation, an adequate solution is yet to be found. Using the framework based on complex algebra in non-sinusoidal circuit analysis (frequency domain), the verification of the energy conservation law is only possible in sinusoidal situations. In this case, reactive energy turns out to be proportional to the energy difference between the average electric and magnetic energies stored in the loads and its cancellation is mathematically trivial. However, in industrial architecture, apparent power definition of electric loads (non-sinusoidal conditions) is inconsistent with the energy conservation law. Up until now, in the classical complex algebra approach, this goal is only valid in the case of purely resistive loads. Thus, in this paper, a new circuit analysis approach using geometric algebra is used to develop the most general proof of energy conservation in industrial building loads. In terms of geometric objects, this powerful tool calculates the voltage, current, and apparent power in electrical systems in non-sinusoidal, linear/nonlinear situations. In contrast to the traditional method developed by Steinmetz, the suggested powerful tool extends the concept of phasor to multivector-phasors and is performed in a new Generalized Complex Geometric Algebra structure ( CGn ), where Gn is the Clifford algebra in n-dimensional real space and C is the complex vector space. To conclude, a numerical example illustrates the clear advantages of the approach suggested in this paper. Keywords: energy conservation law; building loads; harmonics; Geometric Algebra 1. Introduction Most of the harmonic problems affecting electrical distribution networks in industrial architecture are generated within the building. This is partly due to the proliferation of linear and nonlinear loads connected to the circuits in the building; air conditioning, computers, CCTV, servers, adjustable speed drive (ASD), and other electronic equipment, are the main sources of problems. The result of using such a high number of nonlinear loads is that the current waveform is distorted, causing excessive harmonic voltages to be generated. Also, the close proximity of many of these buildings with similar industrial activities will definitely contribute to the distortion of the electric power quality of feed supplying these constructions. These harmonics can cause serious problems in power systems, for example, excessive heat in appliances, premature aging of electrical equipment, derating of equipment, fault of protection, and lower power factor. Consequently, harmonic energy losses reduce the efficiency of power systems and cause significant costs in distribution networks. The disadvantages of harmonics have been discussed extensively, including such factors as the cost effect [1,2]. The challenge of developing a power theory suitable for harmonic analysis has continued over last century. Since the appearance of the first power theory developed by Budeanu in 1927, significant research efforts [ 3 – 11 ] have been made to find out a valid definition of the apparent Symmetry 2016,8, 92; doi:10.3390/sym8090092 www.mdpi.com/journal/symmetry Symmetry 2016,8, 92 2 of 13 power in n-sinusoidal situations by analysing the harmonic effects on linear and nonlinear loads. However, these works have been directed towards particular decompositions of the apparent power that only address the problem partially, focusing on different objectives such as mathematical meaning, physical meaning, power factor improvement, distortionless conditions, power quality improvement, etc. It is our conviction that the apparent power concept should be global and represented by a set of adaptable orthogonal terms for any association criteria. Nevertheless, no power theory has achieved general acceptance yet. Interestingly, although the formulation of every power theory relies on the classic apparent power, S=U rms I rms , the principle of conservation of energy does not apply to this quantity [ 12 ]. Then, the question that arises immediately is: why?. The correct answer is based on the fact that, in sinusoidal conditions and with linear and/or nonlinear loads, the traditional apparent power definition is erratic, except for resistive loads. This is a direct consequence of having only magnitudes for currents and voltages in a circuit branch instead of an expression composed by signed quantities; with this limitation, network analysis involving all the harmonics simultaneously cannot even be performed. Moreover, Kirchhoff’s circuit laws become simply inapplicable because the addition of quantities that represent the time signal of different frequencies has been not defined until now. Thus, the principle of conservation of energy cannot be corroborated for unified systems where sources and loads work simultaneously. This principle states: The instantaneous rate of instantaneous volt-amperes at the input terminal is equal to the sum of the instantaneous volt-amperes at each load component. The limitations of complex algebra and the impossibility of applying the principle of conservation of energy to the apparent power quantity have encouraged the development of an alternative method founded on geometric algebra (GA) concepts. It has been concluded that the typical linear/nonlinear behaviour of power systems require, for its complete analysis, a new mathematical structure that could guarantee the particular character of different components. In this sense, our work considers a new representation of power quantities based on multivectors supported by a Generalized Complex Geometric Algebra (GCGA) [ 11 ]. This transformed frequency domain permits the definition of the apparent power multivector in terms of oriented geometric objects full of power information (magnitude, direction, and sense) from voltage and current vector-phasors, (geometric phasors). The importance of this representation in developing solutions for power theory analysis has been recognized in some recent papers [ 13 – 16 ]. In contrast to classical versions, our representation of apparent power considers the net flow of all the power components representing source–load interactions. This new vector space seems suitable for developing a theory more useful to generalize and to interpret energy conservation in the linear/nonlinear loads of industrial buildings, and also to solve the problem of quantification of losses in the complex case of industrial architecture. 2. Geometric Algebra Foundations William K. Clifford (1845–1879) introduced geometric algebra in a work entitled “On the classification of geometric algebras” [ 17 ]. He unified into the same algebra both Grassmann’s exterior algebra and Hamilton’s quaternions by means of a special product named geometric product. The quickest way to approach the construction of geometric algebra (GA) is through its familiarity with the concept of vector space Vn spanned by orthonormal basis vectors {σ1,σ2,σ3...σn} . In this space a geometric product is defined. This product satisfies the associative, distributive, and contraction rules. It has also the particularity to verify the anticommutative rule. Then, by multiplying vectors Vn , the cited GA is generated as a higher linear space Gn=G(Vn). In GA, the geometric product ab of vectors is defined in terms of the inner product (dot product) and outer product (wedge product): ab =a·b+a∧b(1) Symmetry 2016,8, 92 3 of 13 and reciprocally the inner product and the outer product can be defined in terms of the geometric product as follows: a·b=1 2(ab +ba) a∧b=1 2(ab −ba) The dot product of two vectors is symmetric, while the wedge product is antisymmetric. To better understand the geometric product let us consider a three-dimensional space with basis {σ1,σ2,σ3}satisfying, σ2 1=σ2 2=σ2 3=1, σ1·σ2=σ2·σ3=σ1·σ3=0 For orthogonal vectors the geometric product is a pure bivector (or 2-vector), σ1σ2=σ1·σ2+σ1∧σ2=σ1∧σ2 and that orthogonal vectors anticommute: σ1σ2=σ1∧σ2=−σ2∧σ1=−σ2σ1 so bivectors square to minus 1, (σ1σ2)2=σ1σ2σ1σ2=−σ1σ2σ2σ1=−1 Then, from the basis vectors, another three independent bivectors are generated {σ1σ2,σ2σ3,σ3σ1} and they represent independent directed plane segments, in much the same way as a vector represents directed line segments. The final entity present in the G3algebra is the product of all three vectors: σ1σ2σ3=σ1∧σ2∧σ3 This corresponds to a grade-3 object, called a trivector, and represents a volume element in three dimensions. It also squares to minus 1: (σ1σ2σ3)2=σ1σ2σ3σ1σ2σ3=−σ1σ2σ3σ1σ3σ2=σ1σ2σ3σ3σ1σ2=σ1σ2σ1σ2=−1 In the G3 algebra the trivector corresponds to the highest grade element, usually named pseudoscalar (J), and this top grade coincides with the dimension of the underlying vector space. The full algebra named G3is spanned by the basis set: 1, {σ1,σ2,σ3},{σ1σ2,σ2σ3,σ3σ1},{σ1σ2σ3} and their respective representations are shown in Figure 1. Symmetry 2016,8, 92 4 of 13 Symmetry2016,9,923of12 andreciprocallytheinnerproductandtheouterproductcanbedefinedintermsofthegeometric productasfollows: 1 = ( ) 2 ab ab ba  1 =( ) 2 ab abba  Thedotproductoftwovectorsissymmetric,whilethewedgeproductisantisymmetric. Tobetterunderstandthegeometricproductletusconsiderathree‐dimensionalspacewithbasis  123 ,,σσ σ satisfying, 222 123 1σσσ , 12 23 13 0σσ σσ σσ  Fororthogonalvectorsthegeometricproductisapurebivector(or2‐vector), 12 1 2 1 2 1 2 σσ σσσσ σσ  andthatorthogonalvectorsanticommute: 12 1 2 2 1 21 σσ σ σ σ σ σσ  sobivectorssquaretominus1,  2 12 1212 1221 1σσ σσσσ σσσσ  Then,fromthebasisvectors,anotherthreeindependentbivectorsaregenerated   12 23 31 ,,σσ σσ σσ  andtheyrepresentindependentdirectedplanesegments,inmuchthesamewayasavector representsdirectedlinesegments. ThefinalentitypresentintheG 3 algebraistheproductofallthreevectors: 123 1 2 3 σσσ σ σ σ  Thiscorrespondstoagrade‐3object,calledatrivector,andrepresentsavolumeelementinthree dimensions.Italsosquarestominus1:  2 123 123123 123132 123312 1212 1σσσ σσσσσσ σσσσσσ σσσσσσ σσσσ  IntheG 3 algebrathetrivectorcorrespondstothehighestgradeelement,usuallynamedpseudoscalar (J),andthistopgradecoincideswiththedimensionoftheunderlyingvectorspace. ThefullalgebranamedG 3 isspannedbythebasisset: 1,   123 ,,σσσ ,   12 23 31 ,,σσ σσ σσ ,   123 σσσ  andtheirrespectiverepresentationsareshowninFigure1.  Figure1.BasissetsofG 3 :scalaraspoint,vectorsasdirectedlinesegments,bivectorsasorientedplanes, andtrivectorasorientedvolume. Regardingtheaboverulesformultiplyingvectors,thegeometricproductoftwogenericvectors withinG 3 system, 11 2 2 3 3 =λλ λaσσσ and 11 2 2 33 μμ μbσσσ ,isgivenby   11 22 33 11 22 33 11 22 33 12 12 12 23 23 23 31 31 31 λλ λμμ μ λμ λμ λμ λμ μλ λμ μ λ λμ μλ       ab σσσσσ σ σσ σσ σσ ,  wherethedotproductandthewedgeproductresultinascalartermandasetofbivectorterms, respectively:  11 2 2 3 3 λμ λμ λμ  ab  Figure 1. Basis sets of G3 :scalar as point, vectors as directed line segments, bivectors as oriented planes, and trivector as oriented volume. Regarding the above rules for multiplying vectors, the geometric product of two generic vectors within G3system, a=λ1σ1+λ2σ2+λ3σ3and b=µ1σ1+µ2σ2+µ3σ3, is given by ab =(λ1σ1+λ2σ2+λ3σ3) (µ1σ1+µ2σ2+µ3σ3)= =(λ1µ1+λ2µ2+λ3µ3)+(λ1µ2−µ1λ2)σ1σ2+(λ2µ3−µ2λ3)σ2σ3+(λ3µ1−µ3λ1)σ3σ1 where the dot product and the wedge product result in a scalar term and a set of bivector terms, respectively: a·b=(λ1µ1+λ2µ2+λ3µ3) a∧b=(λ1µ2−µ1λ2)σ1σ2+(λ2µ3−µ2λ3)σ2σ3+(λ3µ1−µ3λ1)σ3σ1 In GA framework ( Gn ), higher dimensional oriented subspaces are also called blades. Therefore, the k-blade term is used for a k-dimensional homogeneous subspace. Thus, a vector is a 1-blade, a bivector is a 2-blade, and so on, up to the pseudoscalar n-blade. The combination of any of these objects configures a multivector. The k-grade part of a multivector is obtained from the grade operatorc, [ 18 ]. Scalar, vectors, bivectors, trivectors, and, in general, multivectors are called simply “geometric objects”. A general multivector Min G3can be written in the expanded form [18], M=α+a+B+Jβ where αand βare scalars, ais a vector, Bis a bivector, and Jis the pseudoscalar. To facilitate algebraic manipulations, it is convenient to introduce an important operation in GA: the reversion (“ † ”), which reverses the order of vectors in any product. Scalars and vectors are invariant under reversion, but bivectors and pseudoscalar change sign: (σ1σ2)†=σ2σ1=−σ1σ2 J†=(σ1σ2σ3)†=σ3σ2σ1=−σ3σ1σ2=σ1σ3σ2=−σ1σ2σ3=−J Thus, the reverse of M,M†, is M†=α+a−B−Jβ Based on the contraction rule, the norm, magnitude, modulus, or length of M is defined by the obviously positive definite expression, ||M||2=|M|2=M2=DMM†E0=α2+a2+B2+β2 M=|M|=sn ∑ k=0|hMik|2 where a2=aa†=aa =a2=λ2 1+λ2 2+λ2 3and B2=BB†=µσpσqµσpσq†=σpσqσqσp=µ2. Symmetry 2016,8, 92 5 of 13 Regarding the GA structure, we define a Generalized Complex Geometric Algebra {CGn,⊗} (GCGA), where C is the complex vector space, Gn is the GA associated to the n-dimensional real space Vn, and ⊗is a new outer product. In brief, the coefficients of multivectors presented above (λ,µ) are replaced by complex numbers in this new framework. The second slight but significative modification over the classic Clifford algebra is made on the geometric product definition. In particular, the wedge product (∧) is replaced by a rotated version (⊗:(<◦∧)) that only applies to the reverse complex bivectors; the dot product remains the same. A detailed description of this new structure (GCGA) is given in reference [11]. In GCGA, typical complex multivectors are in the form of complex vectors V=Vejϕσp=(α+jβ)σp and joint complex scalar-complex bivector object S=S0ejϕσ0+S2ejδB=(α+jβ)σ0+(λ+jµ)B whose respectives magnitudes can be achieved, considering that (“*”) operates over complex numbers and (“†”) operates over multivectors as follows: |V|2=VV∗†=VejϕσpVe−jϕσp=V2 S2=DSS∗†E0=DhS0ejϕσ0+S2ejδBihS0e−jϕσ0−S2e−jδBiE0=S2 0+S2 2 3. Industrial Building Loads: Two Types of Harmonic Sources Nonlinear loads used in industrial buildings introduce current harmonics at the utility. This causes malfunctioning of the sensitive loads connected at the point of common coupling (PCC). Hence, harmonic analysis of the nonlinear loads is essential [ 19 ]. An important aspect is that some nonlinear loads should be modeled by harmonic current sources or harmonic voltage sources [ 20 ]. This option depends on what is being imposed by nonlinear loads. 3.1. Current-Source Load Type (Harmonic Current Source) It is widely accepted that electromagnetic devices (transformers and motors) and consequently the three-phase thyristor converters behave as thyristor rectifiers. The harmonic currents result from the switching operation. Figure 2a shows a simplified three-phase thyristor rectifier and Figure 2b shows their corresponding equivalent circuit, where ui is the voltage source and iN is the harmonic current source, Zi and Zo are input and output impedances, respectively. In Figure 2b, in order not to lose generality, the harmonic current source is represented as a Norton’s equivalent, where the Norton impedance ZN=∞. Symmetry2016,9,925of12 Hence,harmonicanalysisofthenonlinearloadsisessential[19].Animportantaspectisthatsome nonlinearloadsshouldbemodeledbyharmoniccurrentsourcesorharmonicvoltagesources[20]. Thisoptiondependsonwhatisbeingimposedbynonlinearloads. 3.1.Current‐SourceLoadType(HarmonicCurrentSource) Itiswidelyacceptedthatelectromagneticdevices(transformersandmotors)andconsequently thethree‐phasethyristorconvertersbehaveasthyristorrectifiers.Theharmoniccurrentsresultfrom theswitchingoperation.Figure2ashowsasimplifiedthree‐phasethyristorrectifierandFigure2b showstheircorrespondingequivalentcircuit,where i u isthevoltagesourceand N i isthe harmoniccurrentsource, i Z and o Z areinputandoutputimpedances,respectively.InFigure2b, inordernottolosegenerality,theharmoniccurrentsourceisrepresentedasaNorton’sequivalent, wheretheNortonimpedance N Z .  (a)(b) Figure2.(a)Simplifiedthree‐phasethyristorrectifier.(b)Equivalentcircuit. 3.2.Voltage‐SourceLoadType(HarmonicCurrentSource) Ontheotherhand,anothertypeofharmonicsourceisthedioderectifierwithanoutputfilter. Thissystemisanexampleofnonlinearloadwithharmonicvoltage‐sourcebehaviour.Thiskindof systemispresentinalmostallindustrialbuildingnonlinearloads:residential,commercial,and industrialnonlinearloads.Figure3ashowsatypicalperphaseequivalentdioderectifierandFigure 3btheircorrespondingequivalentcircuit,where i u isthevoltagesourceand Th u istheharmonic voltagesource, i Z and o Z areinputandouputimpedances,respectively.InFigure3b,the harmonicvoltagesourceisrepresentedasaThevenin’sequivalent,wheretheTheveninimpedance 0 Th Z .  (a)(b) Figure3.(a)Simplifiedthreephasedioderectifierand.(b)Equivalentcircuit. 4.MultivectorEnergyConservationLawinIndustrialArchitecture Ingeneralform,itcanbeassumedthatanon‐sinusoidalvoltage () i ut isappliedtoageneric linear/nonlinearindustrialbuildingload,asinFigure4. Figure 2. (a) Simplified three-phase thyristor rectifier. (b) Equivalent circuit. Symmetry 2016,8, 92 6 of 13 3.2. Voltage-Source Load Type (Harmonic Current Source) On the other hand, another type of harmonic source is the diode rectifier with an output filter. This system is an example of nonlinear load with harmonic voltage-source behaviour. This kind of system is present in almost all industrial building nonlinear loads: residential, commercial, and industrial nonlinear loads. Figure 3a shows a typical per phase equivalent diode rectifier and Figure 3b their corresponding equivalent circuit, where ui is the voltage source and uTh is the harmonic voltage source, Zi and Zo are input and ouput impedances, respectively. In Figure 3b, the harmonic voltage source is represented as a Thevenin’s equivalent, where the Thevenin impedance ZTh =0. Symmetry2016,9,925of12 Hence,harmonicanalysisofthenonlinearloadsisessential[19].Animportantaspectisthatsome nonlinearloadsshouldbemodeledbyharmoniccurrentsourcesorharmonicvoltagesources[20]. Thisoptiondependsonwhatisbeingimposedbynonlinearloads. 3.1.Current‐SourceLoadType(HarmonicCurrentSource) Itiswidelyacceptedthatelectromagneticdevices(transformersandmotors)andconsequently thethree‐phasethyristorconvertersbehaveasthyristorrectifiers.Theharmoniccurrentsresultfrom theswitchingoperation.Figure2ashowsasimplifiedthree‐phasethyristorrectifierandFigure2b showstheircorrespondingequivalentcircuit,where i u isthevoltagesourceand N i isthe harmoniccurrentsource, i Z and o Z areinputandoutputimpedances,respectively.InFigure2b, inordernottolosegenerality,theharmoniccurrentsourceisrepresentedasaNorton’sequivalent, wheretheNortonimpedance N Z .  (a)(b) Figure2.(a)Simplifiedthree‐phasethyristorrectifier.(b)Equivalentcircuit. 3.2.Voltage‐SourceLoadType(HarmonicCurrentSource) Ontheotherhand,anothertypeofharmonicsourceisthedioderectifierwithanoutputfilter. Thissystemisanexampleofnonlinearloadwithharmonicvoltage‐sourcebehaviour.Thiskindof systemispresentinalmostallindustrialbuildingnonlinearloads:residential,commercial,and industrialnonlinearloads.Figure3ashowsatypicalperphaseequivalentdioderectifierandFigure 3btheircorrespondingequivalentcircuit,where i u isthevoltagesourceand Th u istheharmonic voltagesource, i Z and o Z areinputandouputimpedances,respectively.InFigure3b,the harmonicvoltagesourceisrepresentedasaThevenin’sequivalent,wheretheTheveninimpedance 0 Th Z .  (a)(b) Figure3.(a)Simplifiedthreephasedioderectifierand.(b)Equivalentcircuit. 4.MultivectorEnergyConservationLawinIndustrialArchitecture Ingeneralform,itcanbeassumedthatanon‐sinusoidalvoltage () i ut isappliedtoageneric linear/nonlinearindustrialbuildingload,asinFigure4. Figure 3. (a) Simplified three phase diode rectifier and (b) Equivalent circuit. 4. Multivector Energy Conservation Law in Industrial Architecture In general form, it can be assumed that a non-sinusoidal voltage ui(t) is applied to a generic linear/nonlinear industrial building load, as in Figure 4. ui(t)=√2∑ p∈L∪N Upsin pωt+αp(2) where pis the harmonic order of ui(t). The resulting current has an instantaneous value given by ii(t)=√2∑ q∈N Iqsin qωt+βq(3) where qis the harmonic order of i(t), βq=αq−ϕq in linear case, and ϕq is the impedance phase angle. On the other hand, Nis the set of harmonic indexes present in the voltage and current, and Land M are sets of harmonic indexes only present in the voltage or current, respectively. The capital Up and Iq represent the Root Mean Square (RMS) values of up(t) and iq(t) . In the GCGA structure, we assume that the associated p-th harmonic voltage multivector-phasor (m-phasor) and q-th harmonic current m-phasor are, respectively: Up=Upejαpσpand Iq=Iqejβqσq(4) where Up=Up, Iq=Iq. Then U=∑ p∈L∪N Upand I=∑ q∈N∪M Iq. Symmetry2016,9,926of12  Figure4.Phaseofagenericindustrialbuildingload.    2 ipp pL N ut Usinpt    (2) wherepistheharmonicorderof() i ut.Theresultingcurrenthasaninstantaneousvaluegivenby    2     iqq qN it Isinqt (3) whereqistheharmonicorderofi(t),    qqq inlinearcase,and  qistheimpedancephase angle.Ontheotherhand,Nisthesetofharmonicindexespresentinthevoltageandcurrent,andL andMaresetsofharmonicindexesonlypresentinthevoltageorcurrent,respectively.Thecapital p Uandq IrepresenttheRootMeanSquare(RMS)valuesof   p ut and   q it.IntheGCGA structure,weassumethattheassociatedp‐thharmonicvoltagemultivector‐phasor(m‐phasor)andqthharmoniccurrentm‐phasorare,respectively: p jα p pp UeUσandq jβ qq q IeIσ(4) where p UpU,qq II.Then p pL N   UU andq qN M   II . Applyingthenewgeometricproductconceptand(“*”)beingthecomplexconjugatedoperation, themultivectorapparentpowerattheanonlinearload[11–14],isnowgivenby   qp pq * pq p p p q pp p 0 pp p 0 jφjφj(αβ) pq qp pq pq pq , ,, UIcosφUI sinφj UI e UIe UI e pN L pqN pq qN M pq pq pq pL NqM pq N p Lq N                   *** SUI UI UI U I σσ σσ (5) where    qp pq pp p 0 pp p 0 0 0 jφjφj(αβ) p qqppq pq pq , ,, UIcosφP UI sinφjQj = UIe UIe UIe pq pq pq pL NqM pq N p Lq N            Pσ Qσσ σ Dσσ (6) aretheactive,reactiveanddistortionpowermultivectors,respectively.Notethatactiveandreactive powersaretherealandimaginarypartofthescalarterm,respectively;also,inthedistortionpower, thefirstsummationisduetothelinearelementsandthesecondoneisaconsequenceofa nonlinearbehaviour. Moreover,thesquaredvalue2 Sverifies, Figure 4. Phase of a generic industrial building load. Symmetry 2016,8, 92 7 of 13 Applying the new geometric product concept and (“*”) being the complex conjugated operation, the multivector apparent power at the a nonlinear load [11–14], is now given by S=U I∗=∑ p∈N∪L q∈N∪M UpI∗ q=∑ p=q∈N Up·I∗ p+∑ p6=q Up⊗I∗ q= =∑ p=qUpIpcos ϕpσ0+∑ p=qUpIpsin ϕpjσ0+ +∑ p<q p,q∈NUpIqejϕq−UqIpejϕpσpq +∑ p∈L∪N,q∈M p∈L,q∈N UpIqej(αp−βq)σpq (5) where P=∑ p=qUpIpcos ϕpσ0=Pσ0 Q=∑ p=qUpIpsin ϕpjσ0=Qjσ0 D=∑ p<q p,q∈NUpIqejϕq−UqIpejϕpσpq +∑ p∈L∪N,q∈M p∈L,q∈N UpIqej(αp−βq)σpq (6) are the active, reactive and distortion power multivectors, respectively. Note that active and reactive powers are the real and imaginary part of the scalar term, respectively; also, in the distortion power, the first summation is due to the linear elements and the second one is a consequence of a nonlinear behaviour. Moreover, the squared value |S|2verifies, |S|2=|U I∗|2=∑ p∈N∪L q∈N∪M U2 p∑ p∈N∪L q∈N∪M I2 q=|U|2|I|2=S2(7) The magnitude of multivector S in Equation (7) equals the well-known apparent power S in classical theory. However, the formal product of the voltage and current RMS values is only a simple value without any other attributes. On the other hand, S is full of algebraic, geometric, and physical meaning. This is indeed the most important difference between S =|S|2and S. To probe the energy conservation law [ 21 ] let consider a generic z-element, Figure 5, in which the associated p-th voltage and the q-th current harmonic m-phasors are the complex multivectors Up=Upejαpσp and Iq=Iqej(αq−ϕq)σq , respectively. The pq-th component of the multivector apparent power S, delivered in this element, is given by the geometric product, Sz pq =Uz pIz∗ q(8) Symmetry2016,9,927of12 2 222 22 2 p UI q pN L pN L qN M qN M S         * SUI UI (7) ThemagnitudeofmultivectorSinEquation(7)equalsthewell‐knownapparentpowerSin classicaltheory.However,theformalproductofthevoltageandcurrentRMSvaluesisonlyasimple valuewithoutanyotherattributes.Ontheotherhand,Sisfullofalgebraic,geometric,andphysical meaning.Thisisindeedthemostimportantdifferencebetween2 SSandS. Toprobetheenergyconservationlaw[21]letconsideragenericz‐element,Figure5,inwhich theassociatedp‐thvoltageandtheq‐thcurrentharmonicm‐phasorsarethecomplexmultivectors p jα p pp UeUσandqq j(αφ) qq q Ie  Iσ,respectively.Thepq‐thcomponentofthemultivectorapparent powerS,deliveredinthiselement,isgivenbythegeometricproduct, zzz* p qpq SUI (8) Thissuggestedapproachcanalsobeexpressedthroughpotentialsofthez‐elementnodes1and 2.Itispossibletodesignatetwocurrents,1z z qq  II and2z z qq  II ,whichhavethesamevaluebut oppositesenses. 1z q I 2z q I z p U  Figure5.z‐elementsignscriteria. Thus(8)canbetransformedasfollow:   z z z* 1z 2z z* 1z 1z* 2z 2z* p qpq ppq pq pq   SUI UUIUIUI (9) where   z1z2z ppp UUU isthep‐thharmonicvoltagem‐phasorofthez‐element,n p Uisthen‐nodepthharmonicvoltagem‐phasor,nz q Iistheq‐thharmoniccurrentm‐phasorinz‐thelementthatleaves then‐nodeandnz* q Iistheconjugateofnz q I.Thedefinitionofharmonicmultivectorapparentpower (9),iseasilygeneralizedtoanygenericsystemwithn+1nodes,n=1…Nandz‐elementswithz= 1...Z.Foreachpq‐thmultivectorcomponent: * pq 11 NZ nnz pq nz  SUI (10) In(10),thesummationextendstonnodesthatagreewithexternalterminalsofthez‐thnetwork elementandalsotozelementswhichareconnectedton‐thnode.Oneofthemostimportant propertiesofthismultivectorialpowerrepresentationisitsvalidityifKirchhoff´slawsarefulfilled. Inthiscontext,(10)interpretsthatthetotalapparentpowertransferredbythesourcesofthesystem mustequalthatonedeliveredbytheloads.Thus,itisclearfrom(10)thattheprincipleofenergy conservationissatisfied. 5.NumericalExample Inordertoshowthevalidityoftheenergyconservationlawstatedintheabovesection,a numericalexampleisillustratednext.Thistheoreticalapplicatoncaseconsidersanexampleofa typicalresidential,commercial,andofficeloadwithharmonicvoltagesourcebehaviour,(seeSection 3.2).Asmentionedbefore,thiskindofcircuitispresentinalmostallresidential,commercial,and industrialnonlinearloads,suchascomputers,videomonitors,CCTV,electroniclampballasts,etc. Thesimplifiedequivalentcircuitoftheseloadsisassumedtobeadioderectifierwithanoutputfilter, asseeninFigure6a,b. Figure 5. z-element signs criteria. This suggested approach can also be expressed through potentials of the z-element nodes 1 and 2. It is possible to designate two currents, I1z q=Iz q and I2z q=−Iz q , which have the same value but opposite senses. Thus (8) can be transformed as follow: Sz pq =Uz pIz∗ q=U1z p−U2z pIz∗ q=U1z pI1z∗ q+U2z pI2z∗ q(9) Symmetry 2016,8, 92 8 of 13 where Uz p=U1z p−U2z p is the p-th harmonic voltage m-phasor of the z-element, Un p is the n-node p-th harmonic voltage m-phasor, Inz q is the q-th harmonic current m-phasor in z-th element that leaves the n-node and Inz∗ q is the conjugate of Inz q . The definition of harmonic multivector apparent power (9), is easily generalized to any generic system with n+ 1 nodes, n= 1 . . . N and z-elements with z= 1...Z. For each pq-th multivector component: Spq = N ∑ n=1 Z ∑ z=1 Un pInz∗ q(10) In (10), the summation extends to nnodes that agree with external terminals of the z-th network element and also to zelements which are connected to n-th node. One of the most important properties of this multivectorial power representation is its validity if Kirchhoff ´ s laws are fulfilled. In this context, (10) interprets that the total apparent power transferred by the sources of the system must equal that one delivered by the loads. Thus, it is clear from (10) that the principle of energy conservation is satisfied. 5. Numerical Example In order to show the validity of the energy conservation law stated in the above section, a numerical example is illustrated next. This theoretical applicaton case considers an example of a typical residential, commercial, and office load with harmonic voltage source behaviour, (see Section 3.2). As mentioned before, this kind of circuit is present in almost all residential, commercial, and industrial nonlinear loads, such as computers, video monitors, CCTV, electronic lamp ballasts, etc. The simplified equivalent circuit of these loads is assumed to be a diode rectifier with an output filter, as seen in Figure 6a,b. Symmetry2016,9,928of12  (a)(b) Figure6.(a,b)Loadfedbyrectifierandharmonicvoltagesourceequivalent,bothwithshuntfilter. Figure7showstheinstantaneouswaveformsofthevoltageandtheloadcurrent,obtainedby simulationatthePCC.  Figure7.PCCdistortedvoltageandloadcurrent. ConditioningthecircuitoftheFigure6b,theequivalentcircuitobtained,Figure8,containsz (1,2…7)‐elementsandn=1,2,...5nodes. z=2 ++ 1 U 3 11 I 0 1 2 5 4 z=6 z=3 z=4 z=5 5 U PCC 12 I 22 I 23 I 33 I 34 I 35 I 45 I 46 I 57 I 56 I  Figure8.Conditionedequivalentcircuit. Intheproposedexample,theinstantaneousvoltagesourcesatnodes1and5aregivenby 1 5 u (t) 220 2sin( t ) 5 2sin(3 t ) (V) 63 u ( t ) 100 2 sin( t ) 25 2 sin( 3 t ) (V ) 36        InthenewGCGAfrequencydomainframeworktheircorrespondingm‐phasor  nn p 1,3p  UU are 0 5 10 50 0 50 upcc t() iload t() t Figure 6. (a,b) Load fed by rectifier and harmonic voltage source equivalent, both with shunt filter. Figure 7shows the instantaneous waveforms of the voltage and the load current, obtained by simulation at the PCC. Symmetry2016,9,928of12  (a)(b) Figure6.(a,b)Loadfedbyrectifierandharmonicvoltagesourceequivalent,bothwithshuntfilter. Figure7showstheinstantaneouswaveformsofthevoltageandtheloadcurrent,obtainedby simulationatthePCC.  Figure7.PCCdistortedvoltageandloadcurrent. ConditioningthecircuitoftheFigure6b,theequivalentcircuitobtained,Figure8,containsz (1,2…7)‐elementsandn=1,2,...5nodes. z=2 ++ 1 U 3 11 I 0 1 2 5 4 z=6 z=3 z=4 z=5 5 U PCC 12 I 22 I 23 I 33 I 34 I 35 I 45 I 46 I 57 I 56 I  Figure8.Conditionedequivalentcircuit. Intheproposedexample,theinstantaneousvoltagesourcesatnodes1and5aregivenby 1 5 u (t) 220 2sin( t ) 5 2sin(3 t ) (V) 63 u ( t ) 100 2 sin( t ) 25 2 sin( 3 t ) (V ) 36        InthenewGCGAfrequencydomainframeworktheircorrespondingm‐phasor  nn p 1,3p  UU are 0 5 10 50 0 50 upcc t() iload t() t Figure 7. PCC distorted voltage and load current. Symmetry 2016,8, 92 9 of 13 Conditioning the circuit of the Figure 6b, the equivalent circuit obtained, Figure 8, contains z (1,2 . . . 7)-elements and n = 1,2,...5 nodes. Symmetry2016,9,928of12  (a)(b) Figure6.(a,b)Loadfedbyrectifierandharmonicvoltagesourceequivalent,bothwithshuntfilter. Figure7showstheinstantaneouswaveformsofthevoltageandtheloadcurrent,obtainedby simulationatthePCC.  Figure7.PCCdistortedvoltageandloadcurrent. ConditioningthecircuitoftheFigure6b,theequivalentcircuitobtained,Figure8,containsz (1,2…7)‐elementsandn=1,2,...5nodes. z=2 ++ 1 U 3 11 I 0 1 2 5 4 z=6 z=3 z=4 z=5 5 U PCC 12 I 22 I 23 I 33 I 34 I 35 I 45 I 46 I 57 I 56 I  Figure8.Conditionedequivalentcircuit. Intheproposedexample,theinstantaneousvoltagesourcesatnodes1and5aregivenby 1 5 u (t) 220 2sin( t ) 5 2sin(3 t ) (V) 63 u ( t ) 100 2 sin( t ) 25 2 sin( 3 t ) (V ) 36        InthenewGCGAfrequencydomainframeworktheircorrespondingm‐phasor  nn p 1,3p  UU are 0 5 10 50 0 50 upcc t() iload t() t Figure 8. Conditioned equivalent circuit. In the proposed example, the instantaneous voltage sources at nodes 1 and 5 are given by u1(t) = 220√2sin(ωt+π 6) + 5√2sin(3ωt+π 3) (V) u5(t) = 100√2sin(ωt+π 3) + 25√2sin(3ωt−π 6) (V) In the new GCGA frequency domain framework their corresponding m-phasor Un=∑ p∈{1,3} Un pare U1=U1 1+U1 3=220ej30oσ1+5ej60oσ3(V) U5=U5 1+U5 3=100ej60oσ1+25e−j30oσ3(V))(11) and the resulting current m-phasors Inz =∑ q∈{1,3} Inz qare given by I11 =I11 1+I11 3=100.04e−j133.90oσ1+1.19e−j30.96oσ3=−I12 =I22 =−I23 =I33 (A) I34 =I34 1+I34 3=119.68ej43.40oσ1+2.25e−j40.73oσ3(A) I35 =I35 1+I35 3=20.31e−j150.03oσ1+3.43e−j37.36oσ3=I46 =I57 =−I45 =−I56 (A)     (12) For each n-th node, and q-th harmonic, the equation of first Kirchhoff ´ s law can be rewritten in the following form: Z ∑ z=1 Inz∗ q=0 (13) Multiplying (13) by the p-th voltage harmonic m-phasor in n-th node, Un p , the resulting equation still satisfies the null value, Z ∑ z=1 Un pInz∗ q=0 (14) where the sum is applied to all z-components connecting the n-th node. The set (14) remains unaltered when all n-th nodes are considered, N ∑ n=1 Z ∑ z=1 Un pInz∗ q=0 (15) Equation (15) is the proposed generic expression of the power conservation law.