Geometric Objects: A Quality Index to Electromagnetic Energy Transfer Performance in Sustainable Smart Buildings

symmetry S S Article Geometric Objects: A Quality Index to Electromagnetic Energy Transfer Performance in Sustainable Smart Buildings Juan C. Bravo and Manuel V. Castilla * Higher Polytechnic School, University of Seville, Seville 41011, Spain; [email protected] *Correspondence: [email protected]; Tel.: +34-954-556-440 Received: 29 October 2018; Accepted: 28 November 2018; Published: 29 November 2018   Abstract: Sustainable smart buildings play an essential role in terms of more efficient energy. However, these buildings as electric loads are affected by an important distortion in the current and voltage waveforms caused by the increasing proliferation of nonlinear electronic devices. Overall, buildings all around the world consume a significant amount of energy, which is about one-third of the total primary energy resources. Optimization of the power transfer process of such amount of energy is a crucial issue that needs specific tools to integrate energy-efficient behaviour throughout the grid. When nonlinear loads are present, new capable ways of thinking are needed to consider the effects of harmonics and related power components. In this manner, technology innovations are necessary to update the power factor concept to a generalized total or a true one, where different power components involved in it calculation, properly reflect each harmonic interaction. This work addresses an innovative theory that applies the Poynting Vector philosophy via Geometric Algebra to the electromagnetic energy transfer process providing a physical foundation. In this framework, it is possible to analyse and detect the nature of disturbing loads in the exponential growth of new globalized buildings and architectures in our era. This new insight is based on the concept of geometric objects with different dimension: vector, bivector, trivector, multivector. Within this paper, these objects are correlated with the electromagnetic quantities responsible for the energy flow supplied to the most common loads in sustainable smart buildings. Besides, it must be considered that these phenomena are characterized by a quality index multivector appropriate even for detecting harmonic sources. A numerical example is used to illustrate the clear capabilities of the suggested index when it applies to industrial loads for optimization of energy control systems and enhance comfort management in smart sustainable buildings. Keywords: smart building; harmonics; geometric algebra; Poynting Multivector 1. Introduction 1.1. Motivation Nowadays, professional and academic experts have started to consider the term “Smart Sustainable Cities” [ 1 ] so as to incorporate the different aspects of sustainability in the classical “smart cities” new concept. In fact, literature tends to consider a sustainable city as a whole place of sustainable smart buildings that have a strong environmental focus with a balance within the buildings and the city between infrastructures, information and communications technology (ICT), smart technologies, and urban metabolism, focusing mainly on consumption and energy saving [2]. Buildings all around the world consume a significant amount of energy, which is about one-third of the total primary energy resources [ 3 ]. For this reason, building energy efficiency has turned out Symmetry 2018,10, 676; doi:10.3390/sym10120676 www.mdpi.com/journal/symmetry Symmetry 2018,10, 676 2 of 17 to be a multi-faceted problem and the majority of harmonic problems affecting sustainable smart buildings are generated within new applications and grids of global architectures. This is due, in part, to a proliferation of non-linear loads connected to the networks of the building. These technologies as CCTV recognition systems, automatic smart air conditioning equipment, artificial intelligence (AI), the latest generation of computers, all types of smart detectors and warning systems, and other power electronic equipment are the main sources of such problems. The result of using such highly non-linear load is that the current waveform is distorted. Thus, causing an excessive harmonic of current and voltage. Besides the proximity of many of this new smart building category (industrial and residential smart constructions) will contribute to the distortion of the electric power quality of the feeder, which supplies these constructions and new architectures. These harmonics can cause serious problems in power systems, excessive heat of appliances and machines, aging of electronic component and a decreased capacity, failures of the safety devices and measures of protection, lower power factor and consequently, a reducing power system efficiency due to increasing losses. All these effects are some of the main results of harmonics in power distribution systems. Note, that Harmonic distortion can cause significant costs in distribution networks. Harmonic cost consists of harmonic energy losses, premature aging and de-rating of electrical equipment. The difference between the known generation and the estimated consumption is considered as the energy loss. Other causes of energy loss and, consequently, an increase in the cost of it, are due to the lack of control of energy efficiency in the thermal performance of buildings and energy balance. Most of the building heat losses occur through the building envelope. In recent years, an important number of papers on quantification and optimization of energy efficiency in buildings has been published referring to the standard for buildings. Many of these works have been developed in different areas of application such energy losses in the building envelope as HVAC systems (heating, ventilation and air conditioning), windows, etc. All these works have in common the aim of making efficient buildings from an energy point of view to be sustainable. Particularly, in [ 4 ], an original approach for the U-value evaluation (analogies with coeval buildings, the calculation method, the in-situ measurements and the laboratory tests) is taken. In [ 5 ], measures of energy efficiency and optimization in the building sector are also evaluated. In this article, it is proposed an energy quality multivector index (EQI) based on the Generalized Poynting Multivector (GPM) theory, that possesses clear advantages from the viewpoint of harmonic sources detection and minimization in sustainable smart buildings. 1.2. Literature Review Valuable contributions in electromagnetic field applied to the electric power theory analysis under multi-sinusoidal conditions have been appeared so far in the literature. Despite they have different interpretations, most of them share the common denominator of dealing with the suitability of the Poynting Vector to explain the electromagnetic energy flow in electric systems [6]. In sinusoidal systems, Complex Algebra provides an appropriate framework to analyse the relationship between the Complex Poynting Vector and the energy flow [7]. In one-ports under periodic multi-sinusoidal linear/non-linear operation this issue has still some fundamental unsolved aspects. Nevertheless, some progress has unquestionably been made from numerous valuable contributions in the literature [ 8 – 21 ], each one of them trying to clarify different aspects of the problem by applying the classic Poynting vector (PV) concept. Among them [ 8 , 9 , 14 , 15 ] masterfully explain the power factor concept and the physical mechanism of energy propagation in electric power systems; Ferrero et al. in [ 17 ] reconsiders the physical meaning assigned to the non-active components of the Park instantaneous power; Balci et al. in [ 20 ] describes the transition between PV and instantaneous active and reactive powers; and Faria et al. in [ 11 ] computes the instantaneous power directly from Maxwell’s equations together with the evaluation of the PV flux. On the other hand, several applications are given by means of PV: Lundin et al. in [ 10 ] analyses synchronous generators using field simulations; De León et al. in [ 12 ] identifies terms in nonlinear-switched circuits; Symmetry 2018,10, 676 3 of 17 Cheng et al. in [ 13 ] calculates the reactive power of iced transmission line; Todeschini et al. in [ 16 ] detects and explains the process of compensation and restoration of symmetry in an unbalanced system and Stahlhut et al. in [ 21 ] examines critically the PV possibilities in the area of instrumentation of losses. However, critics of PV calculations [ 18 , 19 ] argue that electromagnetic theory is useless for practical applications of electric power theory. Nevertheless, the large number of papers published on the physical and/or mathematical nature of electromagnetic energy transfer suggest that the work has not been finished. The multidimensionality of power equation and energy quality indexes in the multi-sinusoidal case is the underlying obstacle that considerably complicates the issue at hand. Instead, reference [ 13 ] is a pioneering contribution to the role of a Poynting Multivector to the interpretation of the energy flows based on an original Clifford Vector space. Accordingly, this work adds a new representation of electromagnetic power theory deduced from a Generalized Poynting Multivector [22]. 1.3. Contribution An introduction to “geometric objects” in Geometric Algebra (GA) and the associated phenomena within the electrical systems is developed. In addition, these entities permit a unified treatment of the energy flow concept. By means of the GA framework, the classic explanations of the energy flow process based on interactions between electric and magnetic fields of like-frequency is overcome. The proposed generalization adds the cross-fields interactions in a natural manner, as well. By this way, this paper is concerned to the application of the Generalized Poynting Multivector (GPM) concept [ 22 ] to provide both of them, a physical foundation to the non-active electromagnetic geometric objects as well as a new multivector index to assess the efficiency of the complete energy process in sustainable and smart building loads supplied from a transmission line. 2. Geometric Objects Geometric Algebra [ 23 – 26 ] is based on the concept of objects with different geometrical dimension that result from the geometric product of distinct graded basis elements, e.g., scalars, vectors, 2-vectors and so on. Thus, starting from vectors within an n-dimensional linear space over the real numbers Vn , the geometric product of vectors ab if a,b∈ Vnis formed by a symmetric inner product a·b=1 2(ab +ba)(1) and an antisymmetric outer product a∧b=1 2(ab −ba)(2) Therefore, ab has the canonical decomposition e M=ab =a·b+a∧b(3) The resulting multivector e M is the sum of a scalar ( a·b ) and a bivector or 2-vector ( a∧b ) object. Despite this sum of two distinct objects might seem strange at first sight and it is against the common rule that only same objects should be added, this Clifford ´ s brilliant idea [ 23 ] allows to generalize easily the product to arbitrary higher dimensions and incorporates geometric interpretations of objects and operators. Thus, a bivector can be viewed as directed plane segment, in much the same way as a vector represents a directed line segment, Figure 1. The bivector a∧b has a magnitude |a∧b| equal to the usual scalar area of the circle in Figure 1, with the direction of the plane in which the circle lies, and with sense, which can be assigned to the circle in the plane. Then, just as a vector arepresents (or is represented by) a directed line segment and a bivector a∧b represents a directed plane segment, the trivector (3-vector) a∧b∧cis a grade-3 object that represents a volume (the sphere). Symmetry 2018,10, 676 4 of 17 In G3, the geometric algebra of three-dimensional space, a general multivector can be written as e M=α+λ+B+µJ(4) where ais a vector, Bis a bivector, Iis a trivector (pseudoscalar) and λ and µ are both scalars. The three orthogonal basis vectors {σ1,σ2,σ3} , the three bivectors {σ1σ2,σ2σ3,σ3σ1} , the unit scalar, and the trivector σ1σ2σ3define a graded linear space of total dimension 8=23and are shown in Figure 1. The unit right-handed pseudoscalar for the space Jsquares to − 1, J2=− 1. The pseudoscalar, as well as bivectors, change sign under reversion, but vectors do not. This reversion operation reverses the order of vectors in any product and the convention. So, the reverse of multivector e M derived only from vectors e M†=ba =a·b−a∧b(5) or in the general form resulting from all the possible basis product in G3 e M†=α+λ−B−µJ(6) For all the formulae presented below the followed notation uses the tilde ~ for multivectors, Euclidean vectors are written in bold font, the dagger symbol † denotes the reverse operation and the upper asterisk * represents the complex conjugate operation. For a complete understanding of notation, a list of symbols is summarized in the glossary. Symmetry 2018, 10, 676 4 of 18 where a is a vector, B is a bivector, I is a trivector (pseudoscalar) and λ and μ are both scalars. The three orthogonal basis vectors { } 123 ,, σσσ , the three bivectors { } 12 23 31 ,, σσ σσ σσ , the unit scalar, and the trivector 123 σσσ define a graded linear space of total dimension 8=23 and are shown in Figure 1. The unit right-handed pseudoscalar for the space J squares to −1, 21J=− . The pseudoscalar, as well as bivectors, change sign under reversion, but vectors do not. This reversion operation reverses the order of vectors in any product and the convention. So, the reverse of multivector M  derived only from vectors † M ba a b a b==⋅−∧  (5) or in the general form resulting from all the possible basis product in 3  † M BJ αλ μ =+−−  (6) For all the formulae presented below the followed notation uses the tilde ~ for multivectors, Euclidean vectors are written in bold font, the dagger symbol † denotes the reverse operation and the upper asterisk * represents the complex conjugate operation. For a complete understanding of notation, a list of symbols is summarized in the glossary. Figure 1. Geometric objects. 3. The System Model An “electrical system” is shown in Figure 2 and it can be considered as a space confined by a closed surface to/from which electric power is supplied/received by conductors carrying an electric current. The total instantaneous power transmitted is kk k s (t) u i= (7) Two main electric power processes occur in the system: - Electric power (energy) is dissipated. - Electric power (energy) in the system derives from an electric field and magnetic field interactions. In a general form, the electromagnetic phenomena occurring in the system under study can be roughly modelled by electric circuits. Despite the application of GA to circuit analysis and power theory has a very short history some relevant advances have been made in this area from different perspectives [27–35]. Figure 1. Geometric objects. 3. The System Model An “electrical system” is shown in Figure 2and it can be considered as a space confined by a closed surface to/from which electric power is supplied/received by conductors carrying an electric current. The total instantaneous power transmitted is s(t)=∑ k ukik(7) Two main electric power processes occur in the system: - Electric power (energy) is dissipated. - Electric power (energy) in the system derives from an electric field and magnetic field interactions. In a general form, the electromagnetic phenomena occurring in the system under study can be roughly modelled by electric circuits. Despite the application of GA to circuit analysis and power theory has a very short history some relevant advances have been made in this area from different perspectives [27–35]. Symmetry 2018,10, 676 5 of 17 Symmetry 2018, 10, 676 5 of 18 Figure 2. Electrical system. 4. Generalized Poynting Multivector In this work, space coordinates and time (frequency) coordinates have a differentiated treatment. Space domain remains represented by the classical Euclidean approach, the harmonic field in both, time and frequency domains, are characterized by complex geometric algebra, n  The complete space-time framework results in the hybrid 3 − t n  structure named Generalized Euclidean Space ( 3 − n  for space-frequency domain) [22]. In this 3 − n  algebraic approach, the Generalized Poynting Multivector (GPM) is defined as () , =  * pq pq ΠEH (8) where  p E and q H vectors are called ‘spatial geometric phasors’ of the electric and magnetic harmonic fields, * q H vector is the conjugate of the q-th harmonic spatial geometric phasor q H and ‘’ is the generalized geometric product (GGP) [27]. From this multivectorial field theory, the flux of the GPM quantity for a volume υ enclosed by a surface ς , (GPM theorem), can be expanded into two terms. nn nS ςς ς ςς ς ≠ ⋅= ⋅+ ⋅ =      ppq ppq Πdd d (9) where n is the unitary vector orthogonal to the infinitesimal surface ς d,   is the Poynting Multivector (PM) and   is the Complementary Poynting Multivector (CPM) and * SU I=   is the Power Multivector . In Equation (10) are present three “electromagnetic geometric objects”: the complex vectors  p E and q H, the complex scalar =  * pp p E H (10) and the complex bivector () < =+    ** pqqp pq E HEH (11) From n  structure, Equation (10) explains the total power energy flow because contains the component  related with the power contribution due to frequency-like products and the component   associated to the mutual influence of the harmonic components of the fields. Observe that in classical approach this last term vanishes due to the orthogonality of harmonics basis. Figure 2. Electrical system. 4. Generalized Poynting Multivector In this work, space coordinates and time (frequency) coordinates have a differentiated treatment. Space domain remains represented by the classical Euclidean approach, the harmonic field in both, time and frequency domains, are characterized by complex geometric algebra, CGn . The complete space-time framework results in the hybrid CGt n−R3 structure named Generalized Euclidean Space (CGn−R3for space-frequency domain) [22]. In this CGn−R3 algebraic approach, the Generalized Poynting Multivector (GPM) is defined as e Π=∑ p,qe Epe H∗ q(8) where e Epand e Hqvectors are called ‘spatial geometric phasors’ of the electric and magnetic harmonic fields, e H∗ q vector is the conjugate of the q-th harmonic spatial geometric phasor e Hq and ‘  ’ is the generalized geometric product (GGP) [27]. From this multivectorial field theory, the flux of the GPM quantity for a volume υ enclosed by a surface ς, (GPM theorem), can be expanded into two terms. x ς n·e Πdς=x ς ∑ p n·e Ppdς+x ς ∑ p6=q n·e Dpq dς=e S(9) where n is the unitary vector orthogonal to the infinitesimal surface dς , e P is the Poynting Multivector (PM) and e D is the Complementary Poynting Multivector (CPM) and e S=e UeI∗ is the Power Multivector. In Equation (10) are present three “electromagnetic geometric objects”: the complex vectors e Ep and e Hq, the complex scalar e P=∑ pe Epe H∗ p(10) and the complex bivector e D=∑ p<qe Epe H∗ q+e Eqe H∗ p(11) From CGn structure, Equation (10) explains the total power energy flow because contains the component e P related with the power contribution due to frequency-like products and the component e D associated to the mutual influence of the harmonic components of the fields. Observe that in classical approach this last term vanishes due to the orthogonality of harmonics basis. Nevertheless, a null average value of power without any net energy transferred must be also considered to understand the Symmetry 2018,10, 676 6 of 17 real power flux processes and to fully evaluate the energetic efficiency of electric systems in presence of harmonics. The above mentioned electromagnetic geometric objects are shown in Figure 3and must be considered to explain the main aspects of energy transfer quality. Symmetry 2018, 10, 676 6 of 18 Nevertheless, a null average value of power without any net energy transferred must be also considered to understand the real power flux processes and to fully evaluate the energetic efficiency of electric systems in presence of harmonics. The above mentioned electromagnetic geometric objects are shown in Figure 3 and must be considered to explain the main aspects of energy transfer quality. Figure 3. Complex electromagnetic geometric objects. 5. Formulation of Electromagnetic Energy Quality Index 5.1. Energy Flow on Electrical Systems Consider the case of an electrical system in the form of a circuit, Figure 4, where the load is a linear/nonlinear system. The supply system voltage is periodic but multi-sinusoidal and can be represented by () () p p jωt+α p p ut= 2Im Ue  (12) where p is the harmonic order of u(t). The instantaneous supply current responsible of the generation of the magnetic-field is given by () () qq jωt+β q q it= 2Im Ie  (13) where q is 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, and that components M of current exist without corresponding voltages. In linear operation, β αϕ =− qqq , ϕ q is the impedance phase angle of the consumer electrical system. Figure 4. Flux of the generalized Poynting multivector components of electrical systems. From Equation (10), the energy balance can be expressed as a multivector S  in { } n  , generated by the geometric product “ ” of the voltage and conjugate current geometric phasors [14] given by the following set ** * SU I UI U I==⋅+∧     (14) where  U is the voltage geometric phasor and ∗  I is the conjugated current geometric phasor. The electromagnetic object   can be separated into active and non-active components Figure 3. Complex electromagnetic geometric objects. 5. Formulation of Electromagnetic Energy Quality Index 5.1. Energy Flow on Electrical Systems Consider the case of an electrical system in the form of a circuit, Figure 4, where the load is a linear/nonlinear system. The supply system voltage is periodic but multi-sinusoidal and can be represented by u(t)=√2Im∑ p Upej(ωpt+αp)(12) where pis the harmonic order of u(t). The instantaneous supply current responsible of the generation of the magnetic-field is given by i(t)=√2Im∑ q Iqej(ωqt+βq)(13) 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, and that components M of current exist without corresponding voltages. In linear operation, βq=αq−ϕq , ϕq is the impedance phase angle of the consumer electrical system. Symmetry 2018, 10, 676 6 of 18 Nevertheless, a null average value of power without any net energy transferred must be also considered to understand the real power flux processes and to fully evaluate the energetic efficiency of electric systems in presence of harmonics. The above mentioned electromagnetic geometric objects are shown in Figure 3 and must be considered to explain the main aspects of energy transfer quality. Figure 3. Complex electromagnetic geometric objects. 5. Formulation of Electromagnetic Energy Quality Index 5.1. Energy Flow on Electrical Systems Consider the case of an electrical system in the form of a circuit, Figure 4, where the load is a linear/nonlinear system. The supply system voltage is periodic but multi-sinusoidal and can be represented by () () p p jωt+α p p ut= 2Im Ue  (12) where p is the harmonic order of u(t). The instantaneous supply current responsible of the generation of the magnetic-field is given by () () qq jωt+β q q it= 2Im Ie  (13) where q is 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, and that components M of current exist without corresponding voltages. In linear operation, β αϕ =− qqq , ϕ q is the impedance phase angle of the consumer electrical system. Figure 4. Flux of the generalized Poynting multivector components of electrical systems. From Equation (10), the energy balance can be expressed as a multivector S  in { } n  , generated by the geometric product “ ” of the voltage and conjugate current geometric phasors [14] given by the following set ** * SU I UI U I==⋅+∧     (14) where  U is the voltage geometric phasor and ∗  I is the conjugated current geometric phasor. The electromagnetic object   can be separated into active and non-active components Figure 4. Flux of the generalized Poynting multivector components of electrical systems. From Equation (10), the energy balance can be expressed as a multivector e S in {CGn} , generated by the geometric product “  ” of the voltage and conjugate current geometric phasors [ 14 ] given by the following set e S=e UeI∗=e U·eI∗+e U∧eI∗(14) where e Uis the voltage geometric phasor and eI∗is the conjugated current geometric phasor. The electromagnetic object e Pcan be separated into active and non-active components Symmetry 2018,10, 676 7 of 17 e P=∑ pe Ppact +∑ pe Ppnon−act =Re  ∑ pe Pp   | {z } e Pact +jIm  ∑ pe Pp   | {z } e Pnon−act (15) where the set e Pact transfers the harmonic active power and the set e Pnon−act transport the scalar non-active power (classical reactive power). The only component of e P that transfer useful energy is e Pact, Figure 4. Hence, by virtue of (10), it is obtains x ς ∑ p → 1Z·e Ppdς=x ς ∑ p → 1Z·e Epe H∗ pdς=∑ P UpIpejφpσ0(16) The PM ( e P ) is associated to the energy density at a point on the surface given in terms of the harmonic spatial geometric phasor of electric and magnetic fields at that point. Observe that the PM coincides with the classic Complex Poynting Vector only in pure sinusoidal case. The real part of e P in (16), P=Re(s ς∑ p → 1Z·e Ppdς) , permits a direct interpretation in terms of average energy flow, (i.e. active power P). The imaginary part of e P in (16), Q=Im(s ς∑ p → 1Z·e Ppdς) , is the scalar non-active power. In general, it can be verified that PM e P considers direction and sense, not only of active power components e Ppact but also of scalar non-active components e Ppnon−act Thus, the possible reverse sense of any harmonic e Ppact is very important for a correct identification of harmonic source locations and for determining the responsibility of the utility electrical system and the associated load. Thus, if e Ppact ≥0Pp≥0 the energy flow is unidirectional toward the load. On the other way, if e Ppact ≤0Pp≤0the flow is generated in the nonlinearity of the load. A more detailed develop of the CPM ( e D) is e D=∑ p<qe Epe H∗ q+e Eqe H∗ p=∑ p<q p,q∈Ne Dpq +∑ p∈N,q∈Me Dpq =Re             ∑ p<q p,q∈Ne Dpq +∑ p∈N,q∈Me Dpq              | {z } e Dact +jIm             ∑ p<q p,q∈Ne Dpq +∑ p∈N,q∈Me Dpq              | {z } e Dnon−act (17) Moreover, the flux of the Complementary Poynting Multivector is given by x ς ∑ p<q → 1Z·e Dpq dς=e D(18) where e D can be associated to the polluting nature components, e Dact and e Dnon−act , that do not help in transferring useful energy and e D=∑ p<q linear {UpIqejϕq−UqIpejϕpej(αp−αq)oσpq +∑ p<q non−linear UpIqej(αp−βq)σpq (19) Symmetry 2018,10, 676 8 of 17 is the non-active power bivector called distortion power. Through (9), (16) and (18), it can be observed that the energy flow that it is originated from the surface of the source in the electrical system equals the flow that enters in the load surface. 5.2. Electromagnetic Quality Index (EQI) The aim of this paper is to propose a multivectorial index e ξ , based on electromagnetic geometric objects, capable to characterize the performance of the power transfer or the efficiency of the transmission equipment in single-phase circuits with multi-harmonic signals. This figure of merit can be an important piece of information for the interpretation of the electromagnetic energy flow in presence of disturbing loads. On the contrary, the classic THD index (Total Harmonic Distortion) for the measurement of the harmonic distortion level is not able to give any information about the disturbance nature. To this end, a new index that integrally reflects the nature of the different energy quality situations is defined. This is expressed as e ξ=s ς → 1Ze Πdς Re(s ς∑ p∈N → 1Z·e Ppdς)=1+j Im(s ς∑ p∈N → 1Z·e Ppdς) Re(s ς∑ p∈N → 1Z·e Ppdς)+ s ς∑ p6=q → 1Z·e Dpq dς Re(s ς∑ p∈N → 1Z·e Ppdς)(20) or also, e ξ= Re(s ς ∑ p=1 → 1Z·e P1dς)+jIm(s ς ∑ p=1 → 1Z·e P1dς) Re(s ς ∑ p → 1Z·e Ppdς)+ + Re(s ς ∑ p6=1 → 1Z·e Ppdς)+jIm(s ς ∑ p6=1 → 1Z·e Ppdς) Re(s ς ∑ p → 1Z·e Ppdς)+ + s ς ∑ p6=q → 1Z·e Dpq dς Re(s ς ∑ p → 1Z·e Ppdς) (21) The power factor PF can be written as PF =1 e ξ =Re(s s∑ p → 1Z·e Ppds) s s → 1Z·e Πds (22) It is noteworthy that the power factor is not an exhaustive index for energy quality. In fact, combining (20) and (22) the power factor could be brought to unity, but the electrical systems would still operate in multi-sinusoidal mode. Equations (21), (22) and (23) show that the index e ξ contains terms that keep direction and sense, allowing harmonic source detection. This property is very important to achieve an appropriated compromise between the energy quality index and the power factor. Moreover, the dominant harmonic source should be based on an evaluation of non-active scalar and bivector electromagnetic geometric object components of the Generalized Poynting Multivector. This subject is discussed in the next subsection. Symmetry 2018,10, 676 9 of 17 5.3. Illustrative Comparison between Different Non-Active Electromagnetic Geometric Objects It is seen from (12) and (17) that e U1 and eI1 are the sinusoidal voltage and current geometric phasor and e Up and eIq are the harmonic voltage and current geometric phasor when p6= 1, q6= 1 respectively. Obviously, the voltage geometric phasor e U1 produces the sinusoidal electric field spatial geometric phasor e E1 and the harmonic voltage e Up generates the harmonic field e Ep . Similarly the sinusoidal current eI1 geometric phasor produces the sinusoidal magnetic field spatial geometric phasor e H1 and the harmonic current eIq generates the harmonic field e Hq . The interaction among these fields produces the cited characteristic electromagnetic geometric objects e Π , e P and e D . The objects e P and e D are separated into an active and a non-active component as in Equations (15) and (17). In these equations, the different non-active scalar and bivector electromagnetic objects assume different expressions depending on the electrical system nature. Starting from these considerations, a new philosophy is proposed for the detection of the dominant harmonic source that is based on the comparison of these non-active components to explain the power quality concept. For more detail, the following situations are considered: •Sinusoidal case: p=q=1⇒N={1} e Πsin =e P1=Rene P1o | {z } e P1act +jImne P1o | {z } e P1non−act (23) and e ξsin =1+j Im(s ς∑ p∈1 → 1Z·e P1dς) Re(s ς∑ p∈1 → 1Z·e P1dς)(24) In this case, the quantity Im{e P1}=e P1non−act can be considered as minimum reference value to improve the energy quality, since it is the only non-active electromagnetic geometric object. It is the well established fundamental reactive power (Q) that can be reduced to zero by shunts capacitors Im  x ς ∑ p∈1 → 1Z·e P1dς  =0 (25) and consequently, e ξsin=1⇒PF1=1 (26) where PF1is the fundamental power factor, also known as the displacement power factor PF1=cos(ϕ1)(27) The quantity Rene P1o=e P1act transfers the fundamental active power (useful energy) and is associated to the instantaneous active power. •Multi-sinusoidal linear case: p∈N,q∈N Symmetry 2018,10, 676 16 of 17 References 1. Ibrahim, M.; El-Zaart, A.; Adams, C. Smart Sustainable Cities roadmap: Readiness for transformation towards urban sustainability. Sustain. Cities Soc. 2018,37, 530–540. [CrossRef] 2. Höjer, M.; Wangel, J. Smart sustainable cities: Definition and challenges. In ICT Innovations for Sustainability; Hilty, L., Aebischer, B., Eds.; Springer: Cham, Switzerland, 2015; pp. 333–349. 3. Pervez, H.S.; Nursyarizal, B.M.N.; Perumal, N.; Irraivan, E.; Taib, I. A review on optimized control systems for building energy and comfort management of smart sustainable buildings. Renew. Sustain. Energy Rev. 2014,34, 409–429. [CrossRef] 4. Nardi, I.; Lucchi, E.; Rubeis, T.; Ambrosini, D. Quantification of heat energy losses through the building envelope: A state-of-the-art analysis with critical and comprehensive review on infrared thermography. Build. Environ. 2018,146, 190–205. [CrossRef] 5. Karmellos, M.; Kiprakis, A.; Mavrotas, G. A multi-objective approach for optimal prioritization of energy efficiency measures in buildings: Model, software and case studies. Appl. Energy. 2015 ,139, 131–150. [CrossRef] 6. Poynting, J.H. On the transfer of energy in the electromagnetic field. Philos. Trans. R. Soc. 1884 ,175, 343–361. [CrossRef] 7. Steinmetz, C.P. Theory and Calculation of Alternating Current Phenomena; McGraw Publishing Company: New York, NY, USA, 1900. 8. Emanuel, A.E. Powers in nonsinusoidal situations a review of definitions and physical meaning. IEEE Trans. Power Deliv. 1990,5, 1377–1389. [CrossRef] 9. Emanuel, A.E. Poynting vector and the physical meaning of nonactive powers. IEEE Trans. Instrum. Meas. 2005,54, 1457–1462. [CrossRef] 10. Lundin, U.; Bolund, B.; Leijon, M. Poynting vector analysis of synchronous generators using field simulations. IEEE Trans. Magn. 2007,43, 3601–3606. [CrossRef] 11. Faria, J.A.B. The role of Poynting’s vector in polyphase power calculations. Eur. Trans. Electr. Power 2009 ,19, 683–688. [CrossRef] 12. De León, F.; Cohen, J. AC power theory from Poynting Theorem: Accurate identification of instantaneous power components in nonlinear-switched circuits. IEEE Trans. Power Deliv. 2010,25, 2104–2112. [CrossRef] 13. Cheng, P.; Yang, F.; Luo, H.; Guo, H.; Ran, W.; Yang, Y.; Ullah, I. A method to calculate the reactive power of iced transmission line based on Poynting vector and FDFD. Int. J. Appl. Electromagn. Mech. 2016 ,50, 417–433. [CrossRef] 14. Emanuel, A.E. About the Rejection of Poynting Vector in Power Systems Analysis. Electr. Power Qual. Util. J. 2007,13, 43–49. 15. Emanuel, A.E. Power Definitions and the Physical Mechanism of Power Flow; Wiley: Hoboken, NJ, USA, 2010. 16. Todeschini, G.; Emanuel, A.E.; Ferrero, A.; Morando, A.P. A Poynting vector approach to the study of the Steinmetz compensator. IEEE Trans. Power Deliv. 2007,22, 1830–1833. [CrossRef] 17. Ferrero, A.; Leva, S.; Morando, A.P. An approach to the non-active power concept in terms of the poynting-park vector. Eur. Trans. Electr. Power 2001,11, 291–299. [CrossRef] 18. Czarnecki, L.S. Considerations on the Concept of Poynting Vector Contribution to Power Theory Development. In Proceedings of the 6th International Workshop on Power Definitions and Measurements Under Nonsinusoidal Conditions, Milano, Italy, 13–15 October 2003. 19. Czarnecki, L.S. Could Power Properties of Three-Phase Systems Be Described in Terms of the Poynting Vector? IEEE Trans. Power Deliv. 2006,21, 339–344. [CrossRef] 20. Balci, M.; Hocaoglu, M.; Power, S.A. Transition from Poynting vector to instantaneous power. In Proceedings of the 12th International Conference on Harmonics and Quality of Power, Cascais, Portugal, 1–5 October 2006. 21. Stahlhut, J.W.; Browne, T.J.; Heydt, G.T. The Assessment of the Measurement of the Poynting Vector for Power System Instrumentation. In Proceedings of the 39th North American Power Symposium, Las Cruces, NM, USA, 30 September–2 October 2007; pp. 237–241. 22. Castilla, M.; Bravo, J.C.; Ordonez, M.; Montano, J.C. An approach to the multivectorial apparent power in terms of a generalized poynting multivector. Prog. Electromagn. Res. B 2009,15, 401–422. [CrossRef] 23. Clifford, W.K. Applications of Grassmann’s Extensive Algebra. Am. J. Math. 1878,1, 350. [CrossRef] Symmetry 2018,10, 676 17 of 17 24. Hestenes, D. A unified language for mathematics and physics. their Appl. Math. Phys. 1986. [CrossRef] 25. Chappell, J.M.; Drake, S.P.; Seidel, C.L.; Gunn, L.J.; Iqbal, A.; Allison, A.; Abbott, D. Geometric algebra for electrical and electronic engineers. Proc. IEEE 2014. [CrossRef] 26. Doran, C.; Lasenby, A. Geometric Algebra for Physicists; Cambridge University Press: Cambridge, UK, 2009. 27. Castilla, M.; Bravo, J.C.; Ordóñez, M.; Montaño, J.C. Clifford theory: A geometrical interpretation of multivectorial apparent power. IEEE Trans. Circuits Syst. I Regul. Pap. 2008,55, 3358–3367. [CrossRef] 28. Castilla, M.; Bravo, J.C.; Ordonez, M.; Montaño, J.C.; Lopez, A.; Borras, D.; Gutierrez, J. The geometric algebra as a power theory analysis tool. In Proceedings of the 2008 International School on Nonsinusoidal Currents and Compensation, Lagow, Poland, 10–13 June 2008. 29. Castilla, M.; Bravo, J.C.; Ordóñez, M. Geometric algebra: A multivectorial proof of Tellegen’s theorem in multiterminal networks. IET Circuits, Devices Syst. 2008,2, 383–390. [CrossRef] 30. Bravo, J.C.; Castilla, M.V. Energy Conservation Law in Industrial Architecture: An Approach through Geometric Algebra. Symmetry 2016,8, 92. [CrossRef] 31. Castro-Núñez, M.; Londoño-Monsalve, D.; Castro-Puche, R.M. the conservative power quantity based on the flow of energy. J. Eng. 2016,2016, 269–276. [CrossRef] 32. Petroianu, A.I. A geometric algebra reformulation and interpretation of Steinmetz’s symbolic method and his power expression in alternating current electrical circuits. Electr. Eng. 2015,97, 175–180. [CrossRef] 33. Menti, A.; Zacharias, T.; Milias-Argitis, J. Geometric Algebra: A Powerful Tool for Representing Power Under Nonsinusoidal Conditions. IEEE Trans. Circuits Syst. I Regul. Pap. 2007,54, 601–609. [CrossRef] 34. Castro-Nunez, 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 Regul. Pap. 2012,59, 2056–2064. [CrossRef] 35. Herrera, R.S.S.; Salmeron, P.; Vazquez, J.R.R.; Litran, S.P.P.; Perez, A. Generalized instantaneous reactive power theory in poly-phase power systems. In Proceedings of the 13th European Conference on Power Electronics and Applications, Barcelona, Spain, 8–10 September 2009; pp. 1–10. © 2018 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 (http://creativecommons.org/licenses/by/4.0/).