BEGIN NEW DATA CASE C BENCHMARK DC-43 C Dynamic arc modeling donated to ATP by Mustafa Kizilcay, Universitat C Hannover Institut fur Elektrische Energieversorgung, Welfengarten 1, C 3000 Hannover 1, WEST GERMANY. This data file consists of two data C cases that were adapted to ATP from Hannover's "M28." CDC by WSM and C Mustafa on October 23rd, 1986. Results were explained by Mustafa to C the 1986 LEC meeting that was held in Leuven on October 27th, 1986. C The 1st of 2 subcases illustrates current interruption in an inductiv C circuit consisting of a voltage source in series with an inductance C and the circuit breaker arc (a Type-91, TACS-controlled, time-varying C resistor). Only thermal behavior of circuit breaker is considered. C The time-step size DELTAT has been increased from .2E-7 sec in order C to speed the simulation. This introduces substantial error, but the C shapes of curves remain valid. See comments preceding plot cards. PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) C .2E-7 .25E-3 { Original card. Note increase of DELTAT to speed simulation. .6E-7 .24E-3 { Decrease TMAX slightly, and triple the step size 1 3 1 1 1 -1 5 5 20 20 100 100 500 500 TACS HYBRID GALT +G 90SCHLT 99STROM = 1.E6 * ABS(SCHLT) 99GSTAT = STROM * STROM / (170000. + 890. * STROM) 99GKOR = 1.003 * GSTAT 99SSGNL = ISTEP - 9. 99G' 58+GSTAT 1.0 .29E-6SSGNL GKOR 88LSGNL = G' .LE. 1.E-6 .OR. GALT .LE. 1.E-6 88G = 1.E-6 * LSGNL + G' * NOT(LSGNL) 98RB = INVRS(G) 33RB G GSTAT STROM SCHLT SSGNL LSGNL G' GALT GKOR 77RB 1. 77G' 1. 77G 1. 77GALT 1. BLANK card ending all TACS cards GEN IND 31.8 SCHLT 1.E-6 1 IND SCHLT 1.E10 91IND SCHLT TACS RB { R(t) controlled by TACS variable "RB" } 2 BLANK card ending branch cards IND SCHLT -1. 1.E20 BLANK card ending switch cards 14GEN 141421.4 1. 179.928 -1. BLANK card ending source cards C Total network loss P-loss by summing injections = 2.504876348524E+05 C GEN -141421.2883382 141421.4 885.90010697039 707796.06505328 C 177.71532574734 179.9280000 707795.51064266 89.9282868 C Step Time IND SCHLT IND GEN SCHLT C SCHLT TERRA C C TACS TACS TACS TACS TACS C STROM SCHLT SSGNL LSGNL G' C *** Phasor I(0) = 8.8590011E+02 Switch "IND " to "SCHLT " closed C 0 0.0 0.0 .8859001E-3 .8859001E-3 -141421.288 885.900107 C 0.0 0.0 0.0 0.0 1.0 C *** Open switch "IND " to "SCHLT " after 6.00000000E-08 sec. C 1 .6E-7 0.0 .8856333E-3 .8856333E-3 -141421.288 885.6332743 C 885.6332743 .8856333E-3 -8. 0.0 .8210062181 C 2 .12E-6 1078.390644 .8853654E-3 1078.391529 -141421.288 885.3654244 C 885.3654244 .8853654E-3 -7. 0.0 .8207138637 1 C 4000 .24E-3 -141421.396 -.141436E-6 -141421.396 -141421.396 -.141435538 C .1414355377 -.141436E-6 3991. 1.0 .1175836E-6 C Variable max : 21435.85382 .8859001E-3 21435.85383 -141421.288 885.900107 C 885.6332743 .8856333E-3 3991. 1.0 1.0 C Times of max : .19614E-3 0.0 .19614E-3 0.0 0.0 C .6E-7 .6E-7 .24E-3 .20172E-3 0.0 C Variable min : -164752.212 -.120238E-5 -164752.212 -141421.4 -1.20237999 C 0.0 -.120238E-5 -8. 0.0 .1175836E-6 C Times of min : .20136E-3 .19752E-3 .20136E-3 .19998E-3 .19752E-3 C 0.0 .19752E-3 .6E-7 0.0 .24E-3 PRINTER PLOT C Original Hannover simulation with DELTAT = .2E-7 sec has limits for this 1st C plot variable (IND, SCHLT) of: minimum = -167207.752 at time 0.19952 msec; C maximum = 23261.09003 at time 0.19618 msec. 185 2.186.204. IND SCHLT { Plot limits : (-1.645, 0.214) C Original Hannover simulation with DELTAT = .2E-7 sec has identical minimum C and maximum for 2nd plot, since R varies between model limits of 1.0 & 1.E6: 195 2.192.204. TACS RB { Plot limits : (0.000, 1.000) BLANK card ending plot cards BEGIN NEW DATA CASE C Second of two subcases. Here the dynamic arc model is imbedded in a more C complicated problem. A short transmission line (modeled by single-phase, C cascaded Pi-circuits) produces arc reignition. Arc model has conductance C dependent parameters P(G) & TAU(G). TACS solution according to A. Gruetz. PRINTED NUMBER WIDTH, 13, 2, C 0.1E-6 400.E-6 .2E-6 400.E-6 1 3 1 1 1 -1 5 5 20 20 100 100 TACS HYBRID GA +G 90IND 90SCHLT 99STRM = 1.E6 * ABS(IND - SCHLT) 99GST = STRM**1.351 / 81686. 99GKOR = 1.0134 * GST 88G" = GST - (GST-GA)*EXP(-DELTAT*INVRS(1.95E-6*GA**(-0.2))) 88G' 60+GKOR +GKOR +G" 19. ISTEP 88LSGNL = G' .LE. 1.E-6 .OR. GA .LE. 1.E-6 88G = 1.E-6*LSGNL + G'*NOT(LSGNL) 98RB = INVRS(G) 33RB STRM GST G GA G' G" LSGNL GKOR 77G 1. 77G' 1. 77RB 1. 77GA 1. BLANK CARD ending TACS data GEN IND 8.0 IND 15. 20. 1LANF LPI2 0.15 0.003 1LPI2 LPI3 LANF LPI2 1LPI3 LPI4 LANF LPI2 1LPI4 LPI5 LANF LPI2 1LPI5 LPI6 LANF LPI2 1LPI6 LANF LPI2 IND SCHLT 1.E-6 1 SCHLT LANF 1.E10 91SCHLT LANF TACS RB { R(t) controlled by TACS variable "RB" } 2 BLANK card ending branch cards SCHLT LANF -1. 1.E20 BLANK card ending switch cards 14GEN 177720. 50. -3.6 -1. BLANK card ending source cards C Total network loss P-loss by summing injections = 9.711630317140E+04 C C GEN 177369.31017627 177720. -3989.277855594 63550.495266844 C -11159.13113075 -3.6000000 -63425.16149646 -93.5990147 GEN IND SCHLT LANF C Step Time SCHLT GEN IND SCHLT LANF C LANF C C TACS TACS TACS TACS TACS C GST G GA G' G" C *** Phasor I(0) = -4.0068784E+03 Switch "SCHLT " to "LANF " closed C 0 0.0 0.0 177369.3102 17964.49305 17964.49706 17964.49706 C 0.0 1.0 1.0 1.0 0.0 C *** Open switch "SCHLT " to "LANF " after 2.00000000E-07 sec. C 1 .2E-6 0.0 177370.011 17964.56422 17964.56822 17964.56822 C .9008627603 .9129343213 1.0 .9129343213 .990336131 C 2 .4E-6 -4320.27872 177370.7111 17143.06414 17143.06808 21463.3468 C .8830444741 .8948772701 .9129343213 .8948772701 .9100706638 C 3 .6E-6 -4491.59677 177371.4105 18332.40835 18332.41237 22824.00913 C .9058958171 .918034821 .8948772701 .918034821 .8959289226 BLANK card ending selective node voltage outputs C 2000 .4E-3 4723.577299 177369.3102 26328.99868 26328.99512 21605.41783 C .7689432041 .7551157337 .7537084043 .7551157337 .7551157337 C Variable max : 34213.83832 177720. 39750.9701 39750.96985 28499.88741 C .9058958171 1.0 1.0 1.0 .990336131 C Times of max : .2624E-3 .2E-3 .282E-3 .282E-3 .3062E-3 C .6E-6 0.0 0.0 0.0 .2E-6 C Variable min : -7935.82748 177369.3102 11591.62848 11591.62859 -3236.73514 C 0.0 .6352189E-3 .6352189E-3 .6352189E-3 0.0 C Times of min : .1686E-3 .4E-3 .1754E-3 .1754E-3 .1916E-3 C 0.0 .2392E-3 .2394E-3 .2392E-3 0.0 PRINTER PLOT C Original Hannover simulation with DELTAT = .1E-6 sec has limits for this C plot that are very close: Plot limits : (-0.328, 3.991) C maximum = at time msec. 14535. 350. SCHLT LANF Plot limits : (-0.323, 3.975) BLANK card ending plot cards BEGIN NEW DATA CASE BLANK CARD