diff options
| author | Angelo Rossi <angelo.rossi.homelab@gmail.com> | 2023-06-21 12:04:16 +0000 |
|---|---|---|
| committer | Angelo Rossi <angelo.rossi.homelab@gmail.com> | 2023-06-21 12:04:16 +0000 |
| commit | b18347ffc9db9641e215995edea1c04c363b2bdf (patch) | |
| tree | f3908dc911399f1a21e17d950355ee56dc0919ee /benchmarks/dc43.dat | |
Initial commit.
Diffstat (limited to 'benchmarks/dc43.dat')
| -rw-r--r-- | benchmarks/dc43.dat | 166 |
1 files changed, 166 insertions, 0 deletions
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
+
|