summaryrefslogtreecommitdiffstats
path: root/benchmarks/dc19.dat
blob: d803ae7c531119679f8a9c4c418715a31d60f1e5 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
BEGIN NEW DATA CASE
C     BENCHMARK DC-19                                                           
C      Test  of the control of TACS sources by node voltages and switch currents  
C      of the electric network.  The electric network has no dynamics, with both
C      the  node voltage and switch current equal  to  one  half,  always.  TACS
C      supplemental variable  GENT  adds these to give unity.  The output signal
C      is  then  passed through a function block  1/(1+S)  to  give  ELEMT,  and  
C      through a  S/(1+S)  block to give  GT/ELT.   There also is a disconnected
C      validation of the  Type-17  electric network source,  which provides  for
C      modulation of the following source function by a TACS variable.   For use
C      here,  the TACS variable is the constant  TWO = 2.0  ---  easy to verify.       
     .05     2.0                                                               
       1       1       1       1       1                                        
TACS HYBRID
 1ELEMT   +GENT                                                                 
       1.0                                                                      
       1.0       1.0                                                            
 1GT/ELT  +GENT                                                                 
                 1.0                                                            
       1.0       1.0                                                            
90BUS2                                                                          
91BUS1 
11TWO            2.0       { Constant for modulation of electric network  COSINE                                                                       
99GENT    = BUS1 + BUS2
33ELEMT GENT  GT/ELTBUS1  BUS2                                                  
BLANK card ending all TACS data
  BUS3  BUS2                 1.0                                                
  BUS2  BUS1                 1.0                                                
  COSINE                     1.0      { Dummy branch connects source of interest
  SURGE                      1.0  { Dummy branch connects Type-15 surge function
  TAKUSG                     1.0    { Dummy branch for Taku Noda's Type-15 surge
BLANK card ending electric network branches
  BUS1                                                MEASURING                1
BLANK card ending switches
11BUS3           1.0                                                            
17TWO         { The constant 2.0 of TACS multiplies following 1/2 to give unity:
14COSINE         0.5       0.5  { Cosine starts at peak 0.5, and also ends there
15SURGE          2.0      -1.0      -3.0  { Surge function is a smooth pulse > 0
15TAKUSG         2.0      -1.0      -3.0    -6666. { Same except 2.0 = peak valu
BLANK card ending electric network source cards.
C  Step    Time      BUS1       BUS2       BUS3       COSINE     SURGE      TAKUSG     BUS1       TACS       TACS       TACS
C                                                                                     TERRA      ELEMT      GENT       GT/ELT
C ***        Switch  "BUS1  "  to  "      "  closed before  0.00000000E+00  sec.
C     0     0.0       0.0        0.0        0.0        0.0        0.0        0.0        0.0        0.0        0.0        0.0
C     1     .05       0.0        0.5        1.0        0.0 .181042896 .470363242        0.5 .024390244        1.0 .975609756
 1                                                                              
C    39      1.95         0.0        0.5        1.0 .987688341 .278788345 .724313367        0.5 .854137166        1.0 .145862834
C    40       2.0         0.0        0.5        1.0        1.0 .265713062 .690342786        0.5 .861252427        1.0 .138747573
C maxima :    0.0         0.5        1.0        1.0 .769799804 1.99999856        0.5 .861252427        1.0 .975609756
C maxima :    0.0         .05        .05        2.0        .55        .55        .05        2.0        .05        .05
  PRINTER PLOT
 143 .4 0.0 2.0         SURGE COSINE              { Axis limits: (-1.000, 1.000)
 193 .4 0.0 2.0         TACS  GT/ELTTACS  ELEMT    { Axis limits: (0.000, 9.756)
BLANK card ending plot cards
BEGIN NEW DATA CASE         
C     2nd of 5 subcases of DC-19 tests important supplemental variable functions
C     Unlike 1st,  there is no electric network (in theory).  In practice, ATP
C     will automatically supply a dummy electric network.  This is since around
C     1993 when  TACS STAND ALONE  was internally converted to  TACS HYBRID.
    .100     6.0
       1       1       0       1       0       0       0       0                
TACS STAND ALONE                                                        
  X2      -X3     +TQ                             6.6667 -100. 0.2                
 1X3      +X2                                             0.0   1.05                       
       1.0                                                    
       0.0   1.                                              
11TQ        1.0                                                    0.0          
98SUPVAR  = 0.4 * ( ABS(X3) + 1.E-6 ) ** 3.0 * 2.0
26T            { Type-26 is given value equal to ATP variable in columns 3-8
33X2    X3    T     SUPVAR
C   Step      Time      TACS       TACS       TACS       TACS
C                       X2         X3         T          SUPVAR
C      0       0.0         0.0        0.0        0.0        0.0
C      1       0.1         0.2        .01        0.1  .80024E-6
C      2       0.2         0.2        .03        0.2 .216022E-4
BLANK card ending all TACS data cards          
C     58       5.8  .203206E-4 .999997968        5.8 .799997523
C     59       5.9  .677346E-5 .999999323        5.9 .800000774
C     60       6.0   .22578E-5 .999999774        6.0 .800001858
  PRINTER PLOT                                                                  
 143 1. 0.0 6.0         X2    SUPVAR              { Axis limits: (0.000,  8.000)
BLANK card ending plot cards
BEGIN NEW DATA CASE         
C      3rd of 5 subcases consists only of Type-15 lightning functions, including
C      both old (the 2-exponential surge function)  and new (Bernd Stein of FGH)
C      models.  Bernd's FGH improvement is documented in a letter dated  30.7.86
C      Smoothing of plot data points after 50 ups and downs is also illustrated.
  1.0E-8   3.E-6
       1       1       1       1       1      -1               1
       5       5      20      20                                            
  LIGHTN                     1.0                                                
  ALGHTN                     1.0                                                
  LIGHT1                           1.0 
  ALGHT1                           1.0 
  ALGHT2                     1.0                                                
  ALGHT3                           1.0 
  STAN1                      1.0  { Resistor for Standler surge added 29 Sept 00 
  STAN2                      1.0  { Resistor for 2nd Standler surge added 27 Oct 00 
  TWOE1                      1.0  { Resistor for TWO EXP surge added 19 Jan 01 
  TWOE2                      1.0  { Resistor for 2nd TWO EXP surge added 19 Jan 01 
  HEID1                      1.0  { Resistor for Heidler surge added 19 Jan 01 
  HEID2                      1.0  { Resistor for 2nd Heidler surge added 19 Jan 01 
C      Orlando Hevia mentions legality and danger of negative L of series R-L-C.
C      Add a branch to illustrate as mentioned in newsletter dated January, 2001
  HEVIA                    10.E3  -3.0       { Series R-L with negative L }    1
BLANK card ending electric network branches
BLANK card ending switches
15ALGHTN-1       1.0     2.E-6    50.E-6       5.0
15LIGHTN-1       1.0   -20000.  -500000.  
15ALGHT2-1       1.0     2.E-6    50.E-6      10.0 
15ALGHT1-1       1.0     2.E-6    50.E-6       5.0
15LIGHT1-1       1.0   -20000.  -500000.  
15ALGHT3-1       1.0     2.E-6    50.E-6      10.0 
C 345678901234567890123456789012345678901234567890
C         <--- Off-line Standler 3F10.0 ->
15STAN1 -1       0.6    20.E-6      0.04  Standler { "Standler" is same as -7777
C                 C        A          B      { Identify of 3 Standler parameters
C         <--- In-line Standler 4F8.0 ---><-Request word->  <T-start >< T-stop >
C          T1      T2         Perc  peak
15STAN2 -1  1.0E-6 20.0E-6     0.0  1.000 Standler in-line         0.0       1.0
C
C         <--- Off-line Two exp  3F10.0 ->
15TWOE1 -1    1.20   -36000.0    -5.0E+6          { it is an old two exp source
C                 C        A          B      { Identify of 3 Two exp parameters
C         <--- In-line Two exp  4F8.0 ---><-Request word->  <T-start >< T-stop >
C          T1      T2         Perc  peak
15TWOE2 -1  1.0E-6 20.0E-6     0.0  1.0   Two exp  in-line         0.0       1.0
C
C         <--- Off-line Heidler 3F10.0 -->
15HEID1 -1       0.6   1.0E-06    3.0E-5       5.0  { it is an old Heidler source
C                 C        A          B        N   { Identify of 4 Heidler parameters
C         <-Heidler in line     4F8.0 ---><-key----------><><T-start >< T-stop >
C          T1      T2         Perc  peak                   N
15HEID2 -1  1.0E-6 20.0E-6     0.0  1.000 Heidler  in-line 5       0.0       1.0
C
14HEVIA          1.0       60.    { Voltage source is shorted at 1 us }    1.E-6
BLANK card ending electric network source cards.
  STAN1 STAN2 ALGHT2ALGHT1LIGHT1ALGHTNLIGHTNALGHT3HEID1 HEID2 TWOE1 TWOE2  
C  First  12     output variables are electric-network voltage differences (upper voltage minus lower voltage);
C  Next   1     output variables are branch currents (flowing from the upper node to the lower node);
C   Step      Time      STAN1      STAN2      ALGHT2     ALGHT1     LIGHT1     ALGHTN     LIGHTN     ALGHT3     HEID1      HEID2
C 
C 
C                       TWOE1      TWOE2      HEVIA
C                                             TERRA
C      0       0.0         0.0        0.0        0.0        0.0        0.0        0.0        0.0        0.0        0.0        0.0
C                          0.0        0.0        0.0
C      1     .1E-7  .442479244 .855069065 .25405E-20 .812971E-4 957.508161 .406485E-9 .004787541 .50811E-15 .930195E-8 .140651E-7
C                   .058092768  .04993186 -.16949E-5
C      2     .2E-7  .454691626  .88074615  .2601E-17 .002438392 -4.9670744 .130049E-7 .009550246 .51918E-12 .297563E-6 .449919E-6
C                   .113331409 .097441938 -.51422E-5
C      3     .3E-7  .461895211 .895961187 .14996E-15  .01470791 952.565899 .987364E-7  .01428824 .28952E-10 .225886E-5  .34153E-5
C                   .165855128 .142646809 -.87063E-5
C      4     .4E-7  .467007506 .906801471 .26623E-14 .048743046 -9.8846491 .415991E-6 .019001647 .473526E-9 .951554E-5 .143866E-4
C                    .21579634 .185657436 -.12391E-4
C      5     .5E-7  .470959015  .91521116  .2479E-13 .121908619 947.672889 .126925E-5 .023690588 .395202E-8 .290285E-4  .43887E-4
C                   .263281003 .226579442 -.16201E-4
BLANK card ending output requests                                                                              
C    260    .26E-5   .48556872 .972472279 .991750958 -15.735499 -362.71992 .993881314 .676797074 -18.724315 .572623725 .949657569
C                   1.09277363  .94980214 -.55047746
C    280    .28E-5  .482164403 .967184298 .987942471 -17.163302 -375.61152 .990578487 .698942172 -19.269102 .568828878 .942698718
C                   1.08493561 .942824605 -1.0722499
C    300     .3E-5  .478686003 .961711993 .984066679 -17.995991  -387.2694 .987055255 .718634373 -19.452928 .565055727 .935784527
C                   1.07715275 .935897248 -2.0885866
C Variable maxima :  .506829483 1.00000002        1.0 1757.33053 957.508161        1.0 .718634373 2231.07839        0.6        1.0
C                    1.14948506 .999999945        0.0
C Times of maxima :       .8E-6      .1E-5      .2E-5      .7E-6      .1E-7      .2E-5      .3E-5    .113E-5      .1E-5      .1E-5
C                        .99E-6      .1E-5        0.0
C Variable minima :         0.0        0.0        0.0 -17.995991  -387.2694        0.0        0.0 -19.452928        0.0        0.0
C                           0.0        0.0 -2.0885866
C Times of minima :         0.0        0.0        0.0      .3E-5      .3E-5        0.0        0.0      .3E-5        0.0        0.0
C                           0.0        0.0      .3E-5
  PRINTER PLOT
 145 .4 0.0 3.0         LIGHTNALGHTNALGHT2       { Plot limits: (0.000,   1.000)
 145 .4 0.0 3.0         LIGHT1ALGHT1ALGHT3       { Plot limits: (-0.110,  2.231)
C 145 .4 0.0 3.0         STAN1 STAN2       
C 145 .4 0.0 3.0         HEID1 HEID2       
C 145 .4 0.0 3.0         TWOE1 TWOE2        
BLANK card ending plot cards
BEGIN NEW DATA CASE
C      4th of 5 subcases was added  17 September 1999  to illustrate delayed
C      opening of a Type-13 TACS-controlled switch.  There are two parallel,
C      identical circuits on the electrical side.  A sinusoidal source drives
C      current through a series R-L branch.  The old logic is illustrated by
C      the  QUICK  alternative whereas the new is the  DELAY  alternative.
C      The new optional  DO NOT OPEN UNTIL CURRENT 0  is illustrated by the
C      latter of these two:  opening will be delayed until current passes
C      through zero.  Both circuits oscillate following opening because both
C      have a hanging inductor.  But  QUICK  opens first,  and it oscillates
C      more.  The graph is distinctive.
PRINTED NUMBER WIDTH, 10, 2,  { Request maximum precision (for 8 output columns)
   .0004    .010                   
       1       1    
TACS HYBRID
14CLAMP         100.       50.       30.        0.                 -1.
33CLAMP 
BLANK card ending TACS
  GEN   QUICK                10.    4.   { Series R-L will oscillate when opened                                   
  GEN   DELAY                10.    4.                                     
BLANK card ending branches
13DELAY                 Do not open until current 0   CLOSED          CLAMP   11
$DISABLE
C     The preceding involves a high-level, English-language declaration.  It
C     is equivalent to the following low-level, numeric declation.
13RESIS                    -77333.                    CLOSED          CLAMP   11
C     I.e., the data field of otherwise-unused columns 25-34 can be used to
C     hold the special numeric value.  About the English,  this must begin
C     to the right of column 24.  If 55-60 are used as here for CLOSED, the
C     beginning of  "do not open until current 0"  can be delayed at most 3
C     columns (note 3 blanks between ending 0 and start of CLOSED).  Case is
C     arbitrary (note 1st letter has been made upper case on switch card).
$ENABLE
13QUICK                                               CLOSED          CLAMP   11
BLANK card ending switches
14GEN           100.       50.       0.0        0.                 -1.
BLANK card ending sources
  GEN   QUICK DELAY
BLANK card ending voltage printout
  CALCOMP PLOT  { Needed to cancel  PRINTER PLOT  of preceding subcase
  SUPERIMPOSE                  2    { Superimpose plots of the following 2 cards
 144 1. 0.0 10.         GEN   DELAY QUICK      { First, three node voltages
 194 1. 0.0 10.-200.200.TACS  CLAMP            { 2nd of 2 in overlay is I-branch
BLANK card ending plot
BEGIN NEW DATA CASE         
C      5th of 5 subcases was added  28 July 2002  to illustrate 3 new source
C      types:  CIGRE,  CESI,  and  USRFUN.   These are all special sub-types
C      within the Type-15 source.  The supporting code comes from Orlando P.
C      Hevia of UTN in  Santa Fe, Argentina.  For an alternate and different
C      illustration of the  USRFUN  alternative,  see  DC-5.   For comments
C      about required precision, see note following this subcase.
  .05E-6  10.E-6     { Hevia used dT = 1.0E-9, but such detail is not needed
       1       1       1               1      -1
       5       5      20      20
  CESI1                      1.0
  CESI2                      1.0
  CIGRE1                     1.0
  CIGRE4                     1.0
  USER2                      1.0
  USER8                      1.0   { Add another USRFUN test on 10 December 2002 
BLANK card ending branch cards
BLANK card ending switch cards (none for this data)
C CORRENTE DI CRESTA -0.001 [kA]
C TEMPO DI CRESTA 1 [us]
C TEMPO EMIVALORE 10 [us]
C GENERATORI EQUIVALENTI FRONTE
15CESI1 -1       1.5    1.0E-6    10.0-6  CESI
15CESI2 -1      -1.2    1.0E-6    10.0-6  CESI        
15CIGRE1   1.0          1.0E-6   10.0E-6  CIGRE      20.0E+5    0.5E-6    9.5E-6
15CIGRE4   1.0          1.0E-6   10.0E-6  CIGRE      50.0E+5    0.5E-6    9.5E-6
15USER2    2.0                            USRFUN 
15USER8    8.0                            USRFUN 
BLANK card ending electric network sources
 1 { Request for all possible node voltage outputs (here, just 4)
  CALCOMP PLOT
 145 .5 0.0 5.0         CESI1 CESI2 USER2 CIGRE1
BLANK card ending plot cards
BEGIN NEW DATA CASE
BLANK
EOF

Warning about the preceding.  The CIGRE1 output depends greatly on arithmetic
precision.  Even 64 bits seems inadequate since turning the Salford debugger
on and off (minor update using SUBR5) results in the following differences.
Note that the 3rd row is what will be published in a future newsletter.  The
CIGRE1 value on  step 40  rises from  .661  to  .668  when the debugger is
turned on.

FC version 2.0 - Copyright (c) 1990 Mike Albert        Wed Aug 21 14:08:13 2002
Options are: -c1 -ds5 -t8 

Compare files:
   dc19.lis          71885   8-21-102   2:07p
   dc19.sal          71885   8-21-102   1:48p

Changed lines 859-868
  >     15    .75E-6   1.000594  -.800475  .0648095  .0662431  .0783034
  >     20     .1E-5   1.477403  -1.18192  .1299899  .1324862  .1044063
  >     40     .2E-5   1.338601  -1.07088   .668485  .4280678  .2093254
  >     60     .3E-5    1.23813  -.990504  .9529333  .9505637  .3267037
  >     80     .4E-5   1.167656  -.934125  .8811803  .8800206  .5541328
  >    100     .5E-5   1.104946  -.883957  .8148301  .8147126  .9983771
  >    120     .6E-5   1.044245  -.835396  .7534758  .7542513  .9931261
  >    140     .7E-5   .9840634  -.787251  .6967413  .6982769  .9835613
  >    160     .8E-5   .9240164  -.739213  .6442788  .6464564  .9740871
  >    180     .9E-5   .8640042  -.691203  .5957665  .5984817  .9647042
To
  >     15    .75E-6   1.000594  -.800475  .0644145  .0662431  .0783034
  >     20     .1E-5   1.477403  -1.18192  .1291901  .1324862  .1044063
  >     40     .2E-5   1.338601  -1.07088   .661331  .4280678  .2093254
  >     60     .3E-5    1.23813  -.990504  .9533223  .9505637  .3267037
  >     80     .4E-5   1.167656  -.934125    .88154  .8800206  .5541328
  >    100     .5E-5   1.104946  -.883957  .8151627  .8147126  .9983771
  >    120     .6E-5   1.044245  -.835396  .7537834  .7542513  .9931261
  >    140     .7E-5   .9840634  -.787251  .6970257  .6982769  .9835613
  >    160     .8E-5   .9240164  -.739213  .6445418  .6464564  .9740871
  >    180     .9E-5   .8640042  -.691203  .5960097  .5984817  .9647042

Changed line 872
  >Variable maxima :   1.497392       0.0   .999673  .9989941  .9999945
To
  >Variable maxima :   1.497392       0.0  .9997968  .9989941  .9999945

Comparison complete

21 August 2002, remove from before the 5th, which was not being executed:
BLANK        -----------------    Protect following 4th data subcase from execution 
C          Note:  Following is copy of 3rd subcase between October, 1990,  and
C                 July, 1991.   It was wrong during 9 months.   The solution
C                 prior to October, 1990,  was kept below here,  and discovered
C                 to be correct once again on July 22, 1991.   So,  it would 
C                 seem that Guido's "correction" of Oct, 1990, was wrong,  and
C                 during the spring or summer of 1991 he merely corrected the
C                 correction.   Well,  as history,  we keep the wrong data:
C      3rd of 3 subcases consists only of Type-15 lightning functions, including
C      both old (the 2-exponential surge function)  and new (Bernd Stein of FGH)
C      models.  Bernd's FGH improvement is documented in a letter dated  30.7.86
C      Smoothing of plot data points after 50 ups and downs is also illustrated.
C         Comment cards altered 23 Oct 1990.   Original contents preserved below
C         as 4th data subcase until Bernd Stein approves of changes.
   1.E-8 300.E-8
       1       1       1       1       1      -1
       5       5      20      20                                            
  LIGHTN                     1.0                                                
  ALGHTN                     1.0                                                
  LIGHT1                           1.0                                                
  ALGHT1                           1.0                                                
  ALGHT2                     1.0                                                
  ALGHT3                           1.0                                                
BLANK card ending electric network branches
BLANK card ending switches
15ALGHTN-1       1.0     2.E-6    50.E-6       5.0
15LIGHTN-1       1.0   -20000.  -500000.  
15ALGHT2-1       1.0     2.E-6    50.E-6      10.0 
15ALGHT1-1       1.0     2.E-6    50.E-6       5.0
15LIGHT1-1       1.0   -20000.  -500000.  
15ALGHT3-1       1.0     2.E-6    50.E-6      10.0 
BLANK card ending electric network source cards.
C  Step   Time     ALGHT2     ALGHT1     LIGHT1     ALGHTN     LIGHTN     ALGHT3                                            
C   0     0.0         0.0        0.0        0.0        0.0        0.0        0.0
C   1   .1E-7  .26503E-20 .853443E-4 957.508161 .426722E-9 .004787541 .53005E-15
C   2   .2E-7  .27133E-17 .002559784 -4.9670744 .136524E-7 .009550246 .54161E-12
C   3   .3E-7  .15643E-15 .015440124 952.565899 .103652E-6  .01428824 .30202E-10
C   4   .4E-7  .27773E-14 .051169655 -9.8846491 .436701E-6 .019001647 .493979E-9
 1                                                                              
C 300   .3E-5  1.02657121 -18.891898  -387.2694 1.03619451 .718634373 -20.293153
C  maxima :    1.04319273 1844.81694 957.508161  1.0497837 .718634373 2327.44477
C Times of max :    .2E-5      .7E-6      .1E-7      .2E-5      .3E-5    .113E-5
C  minima :           0.0 -18.891898  -387.2694        0.0        0.0 -20.293153
C Times of min :      0.0      .3E-5      .3E-5        0.0        0.0      .3E-5
  PRINTER PLOT
 145 .4 0.0 3.0         LIGHTNALGHTNALGHT2       { Plot limits: (0.000,   1.050)
 145 .4 0.0 3.0         LIGHT1ALGHT1ALGHT3       { Plot limits: (-0.110,  2.327)
BLANK card ending plot cards
BEGIN NEW DATA CASE