<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">CS</journal-id><journal-title-group><journal-title>Circuits and Systems</journal-title></journal-title-group><issn pub-type="epub">2153-1285</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/cs.2012.33028</article-id><article-id pub-id-type="publisher-id">CS-20196</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject><subject> Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Reactive Compensator Synthesis in Time-Domain as an Alternative to Harmonic Method
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>aciej</surname><given-names>Siwczyński</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Marcin</surname><given-names>Jaraczewski</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Electrical and Computer Engineering, Cracow University of Technology, Poland</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>e-3@pk.edu.pl(AS)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>25</day><month>06</month><year>2012</year></pub-date><volume>03</volume><issue>03</issue><fpage>205</fpage><lpage>209</lpage><history><date date-type="received"><day>July</day>	<month>6,</month>	<year>2011</year></date><date date-type="rev-recd"><day>March</day>	<month>1,</month>	<year>2012</year>	</date><date date-type="accepted"><day>March</day>	<month>8,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The source reactive-current compensation is crucial in energy transmission efficiency. The compensator design in frequency-domain was already widely discussed and examined. This paper presents results of a study on how to design reactive compensators in time-domain. It’s the first time the reactive compensator have been designed in time domain. The example of compensator design was presented.
 
</p></abstract><kwd-group><kwd>Power and Energy Theory; Reactive Power Filters; Reactive Current; Optimization; Harmonic Filter Design</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>This article is a discussion on issue raised in the article of L. S. Czarnecki [<xref ref-type="bibr" rid="scirp.20196-ref1">1</xref>] where the author consider if it is possible to make current decomposition into active, reactive and unbalanced current in time domain and basing on it build reactive compensators. The current decomposition in time domain was presented in the previous articles [2, 3] and in this article is presented the reactive compensators design in time-domain.</p></sec><sec id="s2"><title>2. Reactive Current Compensation in Time-Domain</title><p>In the article [<xref ref-type="bibr" rid="scirp.20196-ref3">3</xref>] was shown that source-receiver current can be decomposed into active and reactive current in “s” domain i.e. for Laplace transform of signals. Reactive current can be compensated with the reactive compensator. The source-receiver current decomposition is given below:</p><disp-formula id="scirp.20196-formula12507"><label>(1)</label><graphic position="anchor" xlink:href="1-7600095\bf05bfd7-0662-440c-bd99-ab1d0dd38a45.jpg"  xlink:type="simple"/></disp-formula><p>where:</p><disp-formula id="scirp.20196-formula12508"><label>(2)</label><graphic position="anchor" xlink:href="1-7600095\6a544a36-7115-414f-8979-8dac692ac7a5.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.20196-formula12509"><label>(3)</label><graphic position="anchor" xlink:href="1-7600095\3d97fac0-bdaa-4d46-9397-a9475a1589e4.jpg"  xlink:type="simple"/></disp-formula><p>stand for the active and reactive parts of receiver admittance operator Y<sup>o</sup>(s).</p><p>The I(s) transform for T-periodic signals is derived using the following relation between non-periodic and periodic signals transform [<xref ref-type="bibr" rid="scirp.20196-ref3">3</xref>].</p><p><img src="1-7600095\2646756c-aefb-4ca4-9863-d9dd835a3bc3.jpg" /> <img src="1-7600095\d83a9587-96eb-4ee0-b717-0d6a48dee709.jpg" /> (4)</p><p>where t &#206; [0, T), Re(σ) &gt; 0, T-time period.</p><p>It can be also calculated directly in time-domain as the T-periodic convolution:</p><disp-formula id="scirp.20196-formula12510"><label>(5)</label><graphic position="anchor" xlink:href="1-7600095\8c31a352-843c-4495-b958-2d8140d76958.jpg"  xlink:type="simple"/></disp-formula><p>where g<sup>o</sup>(t), b<sup>o</sup>(t) stands for T-periodic impulse response of admittance active and reactive part.</p><p>Connecting in parallel the compensator (<xref ref-type="fig" rid="fig1">Figure 1</xref>) the reactive current balance in time-domain states that:</p><disp-formula id="scirp.20196-formula12511"><label>(6)</label><graphic position="anchor" xlink:href="1-7600095\d670026d-88ff-402a-822b-c597109ff124.jpg"  xlink:type="simple"/></disp-formula></sec><sec id="s3"><title>3. T-Periodic Impulse Response of Compensator and Receiver</title><p>Admittance of the elementary RLC compensator branch</p><p>is:</p><p><img src="1-7600095\31d24730-f05a-4b25-949b-7b1ba1be8b46.jpg" />;<img src="1-7600095\ce567913-67e2-42af-a54e-a60a4ede2ba3.jpg" />; <img src="1-7600095\e2dd0f38-7fe2-4498-b231-b169d7d46082.jpg" /></p><p>(We assume later that the compensator is composed of almost lossless elementary branches). And its reactive part is then (3)</p><p><img src="1-7600095\3d9d24e9-b05f-414a-93c0-13be78031e3c.jpg" /></p><p>which leads to general form</p><disp-formula id="scirp.20196-formula12512"><label>(7)</label><graphic position="anchor" xlink:href="1-7600095\6a014948-3dbe-48df-b0be-6238540f32a0.jpg"  xlink:type="simple"/></disp-formula><p>where L(s), M(s)—odd and even polynomials.</p><p>The residues meet the relations: if</p><p><img src="1-7600095\5fc80e27-7ead-4265-83f1-49780bfb3510.jpg" /></p><p>then:</p><p><img src="1-7600095\3fb9c01a-09b5-48bc-aec6-eb9d34e34d56.jpg" />;<img src="1-7600095\5417d5bb-899c-49dd-b558-8e71a5d37d0e.jpg" />; <img src="1-7600095\87932d22-65ab-4c43-b292-70deee9acfc7.jpg" /></p><p>where:</p><p>M′(s) is the derivative of M(s) with respect to s, d— real number.</p><p>Thus (7) reduces to</p><disp-formula id="scirp.20196-formula12513"><label>(8)</label><graphic position="anchor" xlink:href="1-7600095\42c8ff87-b2e3-4912-abc6-db4b37a0f25b.jpg"  xlink:type="simple"/></disp-formula><p>and under (4) and trigonometric identity</p><p><img src="1-7600095\27831d7f-b0c9-4959-86bc-233de2708a8d.jpg" /><img src="1-7600095\c1bd144f-4e84-43b5-b350-5d3a4b20b99a.jpg" /></p><p>where aT = α + jβ.</p><p>We get</p><disp-formula id="scirp.20196-formula12514"><label>(9)</label><graphic position="anchor" xlink:href="1-7600095\7708a7c3-a258-4488-9a7b-1c4322889e55.jpg"  xlink:type="simple"/></disp-formula><p>In the case of almost lossless compensator i.e. for α → 0.</p><p><img src="1-7600095\4c9f9404-eace-49ba-993b-02416330ecd3.jpg" /> <img src="1-7600095\4b9bafdb-1216-490c-bffb-85539ee318b2.jpg" /> (10)</p><p>where:</p><p><img src="1-7600095\48d7eb11-82c7-48a0-ba1e-895a9ad4ad6d.jpg" /></p><p><img src="1-7600095\8f39b9b6-38d9-40e7-bbff-bfd4e40a44ea.jpg" />,<img src="1-7600095\eb44538e-34d8-4dcb-98fe-ff43eb9d9fac.jpg" />—capacitive and inductive reactance for the main frequency<img src="1-7600095\bdb3015c-98e3-460c-bf58-005e97522a6b.jpg" />.</p><p>Residue for B<sup>k</sup>(s) in <img src="1-7600095\f56783b8-16fc-4867-9a13-6bde1d40598a.jpg" /> can be calculated as</p><p><img src="1-7600095\4d960ad7-06ef-4293-a92a-f71e9baf7053.jpg" /></p><p>Thus the T-periodic impulse response of reactive part of the elementary RLC branch (without R) has form</p><p><img src="1-7600095\7243d6b0-7bc5-4978-89b7-ecce0fb70713.jpg" /> <img src="1-7600095\d1c4f332-2e5e-47db-9b70-9b89ce11d95a.jpg" /> (11)</p><p>where:</p><p><img src="1-7600095\46992ac1-2bdc-4f7f-9093-f1d15c2dbdbb.jpg" />,<img src="1-7600095\f47f1161-6fd4-40ea-b372-1d5ef3bf24af.jpg" />—angular and relative resonance frequency of m-th branch, and is depicted in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><p>Later, in the article, it assumes that the reactive part of receiver Y<sup>o</sup>(s) has only real poles, so the receiver is nonoscillatory circuit not as the compensator.</p><p>For the single pole receiver e.g. R<sub>L</sub> or R<sub>C</sub> type the operational admittance is</p><p><img src="1-7600095\4308a705-67df-48e5-a5c0-64de93b4e94e.jpg" />&#160;&#160; <img src="1-7600095\0646c589-f2e2-40ca-b973-d820b32dbe7f.jpg" /></p><p>thus its impulse response is</p><disp-formula id="scirp.20196-formula12515"><label>(12)</label><graphic position="anchor" xlink:href="1-7600095\ccd52262-bbe6-489f-bcd7-17a511c7b874.jpg"  xlink:type="simple"/></disp-formula><p>where A = aT, <img src="1-7600095\8de95f16-2584-4534-a8d8-3854739627d7.jpg" />&#206; [0, 1).</p><p>The coefficient b can be both positive and negative what is shown below for the RRLC receiver (see <xref ref-type="fig" rid="fig3">Figure 3</xref>).</p><p>Its operator admittance is</p><p><img src="1-7600095\a910c391-9a14-4e24-bd05-67fa932d7e7f.jpg" /></p><p>where<img src="1-7600095\ba45a3c3-b8bd-46c2-b5c8-0617ec7ecdf4.jpg" />, <img src="1-7600095\e9829b04-ca8a-491a-a2f5-51e2a189cdea.jpg" />, <img src="1-7600095\68e59767-70a8-4a0c-a687-eb5f89c478bc.jpg" />, <img src="1-7600095\f351827c-b61f-4593-b2bb-f223cc801443.jpg" /></p><p>Then the reactive part of Y<sup>o</sup>(s) is</p><p><img src="1-7600095\2a44f0f6-46db-4624-8d73-16d1f3bceece.jpg" /></p><p>and its T-periodic impulse response is</p><p><img src="1-7600095\d6bcf2e6-92a8-479f-91e1-efaacc23ce2d.jpg" /></p><p>where A<sub>L</sub> = a<sub>L</sub>T, A<sub>C</sub> = a<sub>C</sub>T.</p><p>For the receiver shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>, the operational admittance is</p><p><img src="1-7600095\085d9fd7-a91e-471b-b466-b378e812cf69.jpg" /></p><p>where<img src="1-7600095\0ecbb937-a7a5-4a2f-a2ba-f89dbbaf7735.jpg" />, <img src="1-7600095\218a6cd2-84a8-4f07-bb3e-c9e6282166e5.jpg" />and for the positive poles</p><p><img src="1-7600095\ce61f834-4350-4906-be24-014a385e7460.jpg" />, Y<sup>o</sup>(s) takes form</p><p><img src="1-7600095\3b37ec50-6218-4b77-803a-13aef2ee3433.jpg" /></p><p>where<img src="1-7600095\c225be12-bedd-4a6e-b1a8-4b524434b5bf.jpg" />, <img src="1-7600095\4f7cb134-b66a-4efa-8540-a7af8645bab1.jpg" />and the coefficients b<sub>1</sub>, b<sub>2</sub> are then</p><p><img src="1-7600095\b5295380-1a15-421d-83ca-77621a8a9aa9.jpg" /><img src="1-7600095\7b80ebac-e4d4-4ae1-8b4e-ca56be59c234.jpg" /></p><p>thus</p><p><img src="1-7600095\5df533c4-7f2e-43bb-b4dd-85d7ef58a778.jpg" /></p><p>The reactive function <img src="1-7600095\1701c0b5-42ab-4f48-86d2-112d4eb3ebfa.jpg" /> (12) is shown in <xref ref-type="fig" rid="fig5">Figure 5</xref>. It is necessary to distinguish the two cases when b &gt; 0 and b &lt; 0. When <img src="1-7600095\fbb67877-60b1-4294-bd7b-45e695c9757e.jpg" /> → 0 the curve become straight line.</p><p><img src="1-7600095\6fef1e72-dd8c-457b-aa2d-c3879c8c26e3.jpg" /> <img src="1-7600095\3ea95879-bf98-4c60-aebe-aeebbfc29b20.jpg" /> (13)</p><p>because sh(x) → 0 for x→ 0.</p><p>Thus we arrive to the compensatory balance Equation (6) in a new form</p><disp-formula id="scirp.20196-formula12516"><label>(14)</label><graphic position="anchor" xlink:href="1-7600095\547eb1f7-4bec-4a32-a40d-7d18e1a6b955.jpg"  xlink:type="simple"/></disp-formula><p>where M—total number of compensator branches.</p><p>The solution of (14) (see Figures 6 and 7) for the unknowns L<sub>m</sub> and C<sub>m</sub>, can by find with optimization method. The (14) can be rewritten in respect to D<sub>m</sub> and w<sub>m</sub></p><disp-formula id="scirp.20196-formula12517"><label>(15)</label><graphic position="anchor" xlink:href="1-7600095\f078a765-4355-4fb8-ac0b-4e80fa0af731.jpg"  xlink:type="simple"/></disp-formula><p>where</p><p><img src="1-7600095\e033edea-f73b-4c9c-999e-f025411d2e55.jpg" />, <img src="1-7600095\95360168-4247-4b8c-bd3c-740ce3779efb.jpg" /></p><p>The relative frequencies of L<sub>C</sub> compensator branches have to meet the condition</p><p><img src="1-7600095\349c3d20-b8aa-40b4-a875-9ed2ab8b1235.jpg" /></p><p>Thus we must choose relative resonance frequency w<sub>m</sub> of compensator branches as not the even numbers <img src="1-7600095\0dffc7a0-36c9-41da-80af-926df0b385b4.jpg" />; p—even number.</p><p>The set of Equations (14) can be then solved for D<sub>m</sub> by minimizing<img src="1-7600095\8b9d74d0-20b9-4e55-9a23-1e740916bc2d.jpg" />.</p><disp-formula id="scirp.20196-formula12518"><label>(16)</label><graphic position="anchor" xlink:href="1-7600095\b3f182a5-40df-4123-b6c3-0ae5f1655163.jpg"  xlink:type="simple"/></disp-formula><p>where</p><p><img src="1-7600095\56e10726-1205-46ac-9219-d613d7b6c6b0.jpg" /></p><p>After equate to zero appropriate partial derivatives</p><p><img src="1-7600095\cf56c04c-9982-4579-a801-91bd3df84e5c.jpg" /></p><p>and assuming that</p><p><img src="1-7600095\a3186924-c4c9-4eca-ba3c-ad76d00985c4.jpg" />,<img src="1-7600095\30bd15ef-3a3f-4ea5-8f71-868f49bdebbd.jpg" />;<img src="1-7600095\02e07d95-aa59-4a85-9d84-0d79c5f14359.jpg" /> (17)</p><p>we get necessary minimum condition in form of the set of M linear equations for D<sub>m</sub></p><disp-formula id="scirp.20196-formula12519"><label>(18)</label><graphic position="anchor" xlink:href="1-7600095\50e470b6-94cc-4ecb-9c8e-78797bf226b1.jpg"  xlink:type="simple"/></disp-formula><p>Then we use relation <img src="1-7600095\28c040a8-2a06-4908-8186-0a4163e1f3eb.jpg" /> as to calculate L<sub>m</sub> of compensator branches.<sub></sub></p><p>The offset α in (17) must be less then 0.5 as to assure L<sub>p</sub> positive.</p></sec><sec id="s4"><title>4. Frequency-Domain Compensator Design</title><p>The frequency-domain approach is a well known method (see M. Pasko [4-6]).</p><p>The counterpart of (6) in frequency-domain is</p><disp-formula id="scirp.20196-formula12520"><label>(19)</label><graphic position="anchor" xlink:href="1-7600095\a6216e20-1d3c-4a7e-b990-787a94cacf8d.jpg"  xlink:type="simple"/></disp-formula><p>where:</p><p><img src="1-7600095\e8df648b-ce22-47c4-b1ac-071a49c184d2.jpg" />,<img src="1-7600095\7212a903-588c-4f34-bb74-6d5bc701c6b3.jpg" />—frequency response of compensator and receiver susceptances,</p><p><img src="1-7600095\c92f6fd3-418d-4664-bdb8-5c31928120ae.jpg" />—harmonic number.</p><p>For the elementary compensator branch (L<sub>C</sub> in series) the branch elementary susceptance is</p><p><img src="1-7600095\7292a9dc-fafe-4dbf-a2af-86a5eeae374a.jpg" /></p><p>Thus the formula of reactances balance (19) (see Figures 6 and 7) takes form</p><disp-formula id="scirp.20196-formula12521"><label>(20)</label><graphic position="anchor" xlink:href="1-7600095\b7eb50ec-cb31-4622-9122-91063158b8e5.jpg"  xlink:type="simple"/></disp-formula><p>The Equation (20) is the counterpart of (11) transformed to optimization task (18).</p><p>Then in the particular case of the R<sub>L</sub> in series receiver we get the set of M linear equations for all L<sub>m</sub> of compensator branches</p><disp-formula id="scirp.20196-formula12522"><label>(21)</label><graphic position="anchor" xlink:href="1-7600095\3f93ce63-e478-4d6c-ba36-1ebffe45996d.jpg"  xlink:type="simple"/></disp-formula><p>where n—harmonic numberm—branch number.</p><p>Comparing (18) with (21) we can see that both formulas are the system of linear equations, but in (18) we have</p><p>the impulse response instead of the frequency response of the receiver.</p></sec><sec id="s5"><title>5. Calculation Example</title><p>Let consider the RL load in series for which: P = 500 [W], T = 0.02 [s], ω<sub>1</sub> = 314 [rad/s], AL = T/τ = 10, L<sub>o</sub> = 25.7 [mH] and M = 10, τ—time-constant of the load.</p><p>Effective compensation is up to 5-th harmonic (see Figures 6 and 7).</p><p>The integral in the right side of (18) was calculated numerically using 21 samples and time samples was shifted by T/21/2 due to singularity problem.</p></sec><sec id="s6"><title>6. Conclusion</title><p>The frequency response method used until now to synthesis L<sub>C</sub> compensators and considered the only one [<xref ref-type="bibr" rid="scirp.20196-ref1">1</xref>], has its counterpart in time-domain. In both approaches the L<sub>C</sub> parameters can be found with simple optimization techniques for linear system. The only difference is that in (18) we can use directly the impulse response of the load (differentiated step response) instead of harmonic analysis. Moreover it is the first time in literature that the time-domain reactive compensator design is presented.</p></sec><sec id="s7"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.20196-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">L. S. Czarnecki, “Discussion on ‘a Uniform Concept of Reactive Power of Nonsinusoidal Currents in a TimeDomain’,” Electrical Review, Vol. 85 No. 6, 2009, CDROM.</mixed-citation></ref><ref id="scirp.20196-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">M. Siwczyński and M. Jaraczewski, “The L1-Impulse Method as an Alternative to the Fourier Series in the Power Theory of Continues Time Systems,” Bulletin of the Polish Academy of Science, Technical Sciences, Vol. 57, No 1, 2009, pp. 79-86.</mixed-citation></ref><ref id="scirp.20196-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">M. Siwczyński, “Decompositions: Active Current, Scattered Current, Reactive Current in Time-Domain—Single-Phase Circuits,” Electrical Review, Vol. 6, 2010, pp. 11-17. </mixed-citation></ref><ref id="scirp.20196-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">M. Pasko, “Choice of the Two-Terminal Components, Compensating Reactive Current Component of the Linear Load Fed with Distorted Voltage,” Seminarium z Podstaw Elektrotechniki i Teorii Obwodow, Gliwice-Wis?a, 20-23 May 1992, CD-ROM.</mixed-citation></ref><ref id="scirp.20196-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">M. Pasko, “Improvement of the Working Conditions of the Real Distorted Voltage Sources Using Two-Terminal LC components,” Scientific Notebook of the Silesium University of Technology Electrotechnics, Vol. 117, 1991, pp. 45-62.</mixed-citation></ref><ref id="scirp.20196-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">M. Pasko and J. Walczak, “Optimization of the Power Quality Coefficients of the Electrical Circuits with Nonsinusoidal Periodic Signals,” Wydawnictwo Politechniki ?l?skiej, Gliwice, 1996.</mixed-citation></ref></ref-list></back></article>