1 SUBROUTINE tsw_model (ht, time, parms, rvar, currentr2)
3 c january 26, 1993; hnw
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
16 c define passed variables
18 REAL ht,
time, parms(*), rvar, currentr2
20 c define internal variables
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
27 parameter(pie=3.14159265)
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
39 c this
model is based on the screamer
module of rls and modifies the
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;
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.
51 c this
model was developed by t h martin.
64 go to(10,20,30,40,50,60), diel1
65 c setting up for the various dielectrics
103 c calculate values for
use
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
111 IF(
time.LT.tbd1)
THEN
116 iamp1 = abs(currentr2)
118 irtb1 = (ht*(igaus1**(.66667)))+irta1
120 rchn1 = cona1*irtb1**(.5)+1.0e-5
121 resc1 = dcm1/rchn1**2/mhos1/pie/sf
124 c providing outputs for screamer
127 c Save values for next
time step
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
c *****************************************************************************c Various format statements for read_screamer_data output c To get these into use
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
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
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
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
subroutine tsw_model(ht, time, parms, rvar, currentr2)
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
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