Screamer Project  V3.3.1
Screamer Structure
 All Files Functions Variables
lsyswmdl.f
Go to the documentation of this file.
1  SUBROUTINE tsw_model (ht, time, parms, rvar, currentr2)
2 c
3 c january 26, 1993; hnw
4 c
5 c modifications:
6 c 1994-03-07 kws: put values that are needed from one use of
7 c the routine to the next in variable arrays. these arrays are also available
8 c for output. deleted versions 2, 3, 4, so that only this routine is used when
9 c the model is called. note that the arrays are specific to each call of the
10 c model, so that there is no cross-talk between different uses.
11 c 2014-05-02 rbs: change real*4 to real
12 c
13  include 'zdemmax.h' ! parameters
14  include 'zdemout.h' ! common blocks
15 c
16 c define passed variables
17 c
18  REAL ht, time, parms(*), rvar, currentr2
19 c
20 c define internal variables
21 c
22  REAL iamp1, mhog1, igaus1, cona1, irta1, irtb1, rchn1, resc1, pie
23  REAL zeta1, mhos1, rho1, rhob1, sf
24  real tbd1, dcm1, psia1, dcm1n, xsw1, xch1
25 c
26  INTEGER diel1
27  parameter(pie=3.14159265)
28 c
29 c switch loss model subroutine for screamer. this model uses the braginskii
30 c formulations as modified by t h martin to determine the switch loss in
31 c h2o(1),oil(2),sf6(3),air(4), he(5), or h2(6). the model calculates
32 c the channel radius and then provides the channel resistance by assuming the
33 c appropriate constant plasma conductivity between 150 and 600 mhos-cm.
34 c the input of switch current is taken to the 2/3 power then integrated
35 c and multiplied by the appropriate constant to provide the channel radius
36 c squared.
37 c
38 c
39 c this model is based on the screamer module of rls and modifies the
40 c resistor r2.
41 c
42 c the input data consists of the type of dielectric,
43 c diel1 = h2o(1), oil(2), sf6(3), air(4), he(5), or h2(6);
44 c tbd1 = the breakdown time in seconds;
45 c dm1 = the total gap length d in m;
46 c atm1 = the gas pressure in atmospheres. Use 1.0 for h2o or oil.
47 c xsw1 = no of sw modules
48 c xch1 = no of channels per sw
49 c the switch is active for time greater than or equal to tbd1.
50 c
51 c this model was developed by t h martin.
52 c
53  diel1 = parms(1)
54  tbd1 = parms(2)
55  dm1 = parms(3)
56  atm1 = parms(4)
57  xsw1 = parms(5)
58  xch1 = parms(6)
59  irta1 = parms(7)
60 c
61  dcm1 = dm1 * 100.0
62  psia1 = atm1 * 14.70
63 c
64  go to(10,20,30,40,50,60), diel1
65 c setting up for the various dielectrics
66 c
67  10 CONTINUE ! H2O
68  zeta1=4.5
69  mhos1=600.
70  rhob1=1.
71  go to 100
72 c
73  20 CONTINUE ! OIL
74  zeta1=4.5
75  mhos1=600.
76  rhob1=.9
77  go to 100
78 c
79  30 CONTINUE ! SF6
80  zeta1=4.5
81  mhos1=160.
82  rhob1=6.16e-3
83  go to 100
84 c
85  40 CONTINUE ! AIR
86  zeta1=4.5
87  mhos1=200.
88  rhob1=1.2929e-3
89  go to 100
90 c
91  50 CONTINUE ! HE
92  zeta1=4.5
93  mhos1=140.
94  rhob1=1.78e-4
95  go to 100
96 c
97  60 CONTINUE ! H2
98  zeta1=4.5
99  mhos1=300.
100  rhob1=8.99e-5
101  go to 100
102 c
103 c calculate values for use
104 c
105  100 rho1 = rhob1/14.7*psia1
106  mhog1 = mhos1*.8985e12
107  cona1 = (4./mhog1/rho1/(pie*pie)/zeta1)**.1667
108  sf = (xsw1*xch1)**.3333333
109 c
110 c putting in the breakdown time for this switch
111  IF(time.LT.tbd1) THEN
112  rvar = 1.0e6
113  rchn1 = 0.0
114  goto 200
115  ELSE
116  iamp1 = abs(currentr2)
117  igaus1= iamp1*3.0e9
118  irtb1 = (ht*(igaus1**(.66667)))+irta1
119  irta1 = irtb1
120  rchn1 = cona1*irtb1**(.5)+1.0e-5
121  resc1 = dcm1/rchn1**2/mhos1/pie/sf
122  end if
123 c
124 c providing outputs for screamer
125  rvar = resc1
126 c
127 c Save values for next time step
128 c
129  200 continue
130  parms(7)=irta1
131 c
132 c Save switch channel radius in meters
133 c
134  parms(8)=rchn1/100.
135 c
136  RETURN
137  END
c *****************************************************************************c Various format statements for read_screamer_data output c To get these into made format to be characters c for each line corrected spelling errors in format added statement for Zflow Plasma Loss Model added format for CSV output type fixed more lines longer than characters added format for Measure Zflow Block and forward c reverse current directions in Zflow plasma loss c and Zflow POS models added format for SFC output type c removed from all code calls c c c c c a80 c i10 c No grids on plots c Do not write files containing the plotted points c Execute only one cycle c Do not echo the setup parameters and indicies c c &exitting c a13 c c c c c c c c102 c &described as a function of time c c c &described as a function of time c c c &function of time c c shell c &min A K gap
Definition: zdemfmt.h:84
c *****************************************************************************c Various format statements for read_screamer_data output c To get these into use
Definition: zdemfmt.h:7
c *****************************************************************************c Various format statements for read_screamer_data output c To get these into made format to be characters c for each line corrected spelling errors in format added statement for Zflow Plasma Loss Model added format for CSV output type fixed more lines longer than characters added format for Measure Zflow Block and forward c reverse current directions in Zflow plasma loss c and Zflow POS models added format for SFC output type c removed from all code calls c c c c c a80 c i10 c No grids on plots c Do not write files containing the plotted points c Execute only one cycle c Do not echo the setup parameters and indicies c c &exitting c a13 c c c c c c c c102 c &described as a function of time c c c &described as a function of time c c c &function of time c c shell c &min A K c &trapped field time
Definition: zdemfmt.h:85
c *****************************************************************************c Various format statements for read_screamer_data output c To get these into made format to be characters c for each line corrected spelling errors in format added statement for Zflow Plasma Loss Model added format for CSV output type fixed more lines longer than characters added format for Measure Zflow Block and forward c reverse current directions in Zflow plasma loss c and Zflow POS models added format for SFC output type c removed from all code calls c c c c c a80 c i10 c No grids on plots c Do not write files containing the plotted points c Execute only one cycle c Do not echo the setup parameters and indicies c c &exitting c a13 c c c c c c c c102 c &described as a function of time c c c &described as a function of time c c c &function of time c c shell c &min A K c &trapped field c153 c c c Sin c c c remaining calls in rdscrelem c Tabular a10 Cond c c c &used with this model c Exponential model of a resistive c c c Rise model of a resistive c c Z FLOW POS model
Definition: zdemfmt.h:131
c *****************************************************************************c Various format statements for read_screamer_data output c To get these into made format to be characters c for each line corrected spelling errors in format added statement for Zflow Plasma Loss Model added format for CSV output type fixed more lines longer than characters added format for Measure Zflow Block and forward c reverse current directions in Zflow plasma loss c and Zflow POS models added format for SFC output type c removed from all code calls c c c c c a80 c i10 c No grids on plots c Do not write files containing the plotted points c Execute only one cycle c Do not echo the setup parameters and indicies c c &exitting c a13 c c c c c c c c102 c &described as a function of time c c c &described as a function of time c c c &function of time c c shell length
Definition: zdemfmt.h:83
c *****************************************************************************c Various format statements for read_screamer_data output c To get these into made format to be characters c for each line corrected spelling errors in format added statement for Zflow Plasma Loss Model added format for CSV output type fixed more lines longer than characters added format for Measure Zflow Block and forward c reverse current directions in Zflow plasma loss c and Zflow POS models added format for SFC output type c removed from all code calls c c c c c a80 c i10 c No grids on plots c Do not write files containing the plotted points c Execute only one cycle c Do not echo the setup parameters and indicies c c &exitting c a13 c c c c c c c c102 c &described as a function of time c c c &described as a function of time c c c &function of time c c shell c &min radius
Definition: zdemfmt.h:84
subroutine tsw_model(ht, time, parms, rvar, currentr2)
Definition: lsyswmdl.f:1
c *****************************************************************************c Various format statements for read_screamer_data output c To get these into made format to be characters c for each line c
Definition: zdemfmt.h:7
c *****************************************************************************c Various format statements for read_screamer_data output c To get these into made format to be characters c for each line corrected spelling errors in format added statement for Zflow Plasma Loss Model added format for CSV output type fixed more lines longer than characters added format for Measure Zflow Block and forward c reverse current directions in Zflow plasma loss c and Zflow POS models added format for SFC output type c removed from all code calls c c c c c a80 c i10 c No grids on plots c Do not write files containing the plotted points c Execute only one cycle c Do not echo the setup parameters and indicies c c &exitting c a13 c c c c c c c c102 c &described as a function of time c c c &described as a function of time c c c &function of time c c shell c &min A K c &trapped field c153 c c c Sin c c c remaining calls in rdscrelem c Tabular a10 Cond c c c &used with this model c Exponential model of a resistive switch
Definition: zdemfmt.h:118