From b18347ffc9db9641e215995edea1c04c363b2bdf Mon Sep 17 00:00:00 2001 From: Angelo Rossi Date: Wed, 21 Jun 2023 12:04:16 +0000 Subject: Initial commit. --- benchmarks/dc43.dat | 166 ++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 166 insertions(+) create mode 100644 benchmarks/dc43.dat (limited to 'benchmarks/dc43.dat') diff --git a/benchmarks/dc43.dat b/benchmarks/dc43.dat new file mode 100644 index 0000000..63fcea5 --- /dev/null +++ b/benchmarks/dc43.dat @@ -0,0 +1,166 @@ +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 + -- cgit v1.2.3