DE-R Four Module + PRS Load Simulation
!
! Rick Spielman 2018-06-18
!
Time-step 2.0e-12
Resolution-time 100e-12
End-time 1.0e-6
Number-prints 5
Execute-cycles all
Grids no
Echo-setup no
Detail-prints full
Max-points 5001
!
!Start circuit definition
!
! DE-R_rev6.txt is based on DE-R_rev5final.txt but we are changing the sub-Marxes to 20 stages
! and charging to 85 kV.
!
! DE-R_rev5final.txt is based on DE-R_rev5b.txt but the water line details are
! changed to reflect the final value of the conceptual design. Also adding a length scaling factor
! to correct for the overestimate of the water switch resistance of 40%. Decrease length by .4.
!
! DE-R_rev5b.txt is based on DE-R_rev5a.txt. In this run deck we are splitting the two halves
! of the water lines AND Marxes in half so that we can drive each side of the vacuum disk feed 
! separately. This allows more accurate treatment of the individual MITLs and the convolute.
!  MITL and Z_flow models will be incorporated into the deck. We will divide C by 2X and multiply L
! and R by 2X. Remember that the IOUT/IIN calls are backward in the second branch.
! In this run each Marx has a capacitance of 70.5 nF and the L1 with ~ 42.85 nF
!
! DE-R Rev5a.txt is based on DE-R Rev4a.txt. We are including a PRS load using the general 
!  parameters of the Z51 deck used in Z15_rev2_Z51.txt. The exact inductance numbers have
!  been tweaked a bit to be self consistent. We output V, I before the MITL TR and V,I after the MITL.
!
! DE-R Rev4a is based on DE-R Rev4. I am removing all of the extra branches to speed iteration. 
! This forces me to divide all marx and water line impedances and inductances by four
! and to increase the Marx C by 4.
!
! DE-R Rev 3a is based on DE-R Rev 2a. We are lowering the WL impedance to 0.4 Ω (1.6/4) and
! explicitly putting in the water flare impedance of the stack. In addition, the insulator, vacuum flare, 
! and MITL impedances will be rechecked. Goal - explicit identification of all impedance and 
! inductance components.
!
! DE-R Rev 2a uses the DE-R Rev 2 deck with 4 parallel modules that are driving an insulator stack,
! MITL, and e-beam load. We will use 4, separate modules joined to the insulator stack with
! End Branches. This way we can keep the Rev1b structure exactly. We also allow jitter in the 
! Future. Note: the module end branches are in reverse order away from the branch.
!
! DE-R Rev 1 uses the DE Rev 4 as the baseline. We will be changing all the line impedances to
! 2 Ω and will be changing the line lengths. The goal here is to have a driver that delivers a clean 100-ns FWHM pulse to a constant impedance load.
!
! Rev 4 changes the CPL water switch gap to the actual gap 3.75”
! also changes the output line to constant 1.3 Ω impedance with a matched load.
!
! Rev 3 Marx switch losses and more accurate water line parameters
! Change marx to 30 stages of 300 nF/stage
!
! Gas switch dimensions from drawings
!
! The sub Marx capacitance is 290 nF/40= 7.25 nF. We operate at a charge voltage of 85 kV.
! We have a 7.25-nF sub Marx with 40 cap pairs in series and 20 switches. Vch = 3.4 MV
! We have 8 sub Marxes in parallel per module.
! The total Marx module capacitance is 58 nF.
!
! Include 130 mΩ individual cap ESR or 65 mΩ per pair of caps in parallel in a half stage
!   X 20 stages per sub-Marx = 2.6 Ω per sub-Marx
!
! With 167 nH every stage. L = 167 nH * 20 = 3.34 µH.
! The total Marx module inductance (caps and switches) would be 3.34 µH/8 = 417 nH
!
! The total Marx capacitance is 58 nF. This implies
!   that there will a ringing gain on the CPL. This allows a faster rise time than otherwise to 
!   get to the desired CPL voltage.
!
! CPL capacitance used is 42.85 nF. This is distributed in a TL that is ~60 ns long and
!   has a constant impedance of 1.4 Ω. (C = t/Z = 60/1.4 nF = 42.85 nF) 
!  Optimum C ratio is 1.5X to Cmarx = 64.275 nF (we are a bit low in the Marx but OK)
!
! Energy per Marx module at 85 kV = 335.2 kJ, 2 marxes = 670.5 kJ, 4 marxes = 1.34 MJ
!
! The matched impedance of a sub-Marx = sqrt(L/C) = 21.46 Ω.
!   The total Marx impedance is 21.46 Ω/8 = 2.68 Ω
!   Note CPL impedance is 1.4 Ω.
!   The ESR amounts to ~ 12.1% of the matched impedance.
!
! The switch model will have 20x the length of one switch per sub Marx and then will
! have 8 arc channels (sub Marxes) in parallel.
! Then there will be two Marxes in parallel representing half of the machine
! 
! We will assume a TOTAL inductance of 417 + 100 = 517 nH for one Marx and one L1.
! The total inductance for Branch A (and also Branch B) two modules in parallel is 258.5 nH.
!
! 
! Main Branch - Branch #1
!
BRANCH
!
! B LEVEL
!
! Marx capacitance and charge voltage
! CX2 for two modules (58 nF each Marx module) = 116 nF
! +/-85 kV charge x 40 = 3.4 MV
!
RCGround 1e+12 116e-9
Initial VC1 3.4e6

TXT VC1 SCAle 1.0e-06
$V_marx_B (MV)
TXT EC1 SCAle 1.0e-03
$E_marx_B (kJ)

!
! Cap inductance, case & parasitic inductance, and ESR
! ESR = 0.325 for the total Marx module (2.6 Ω sub Marx/8)
! Inductance is 417 nH per module 
! Divide R and L by 2 for 2 modules (driving one MITL)
!

RLSeries 0.1625 208.5e-9

!
! Switch resistance - Martin Model and switch inductance - the total gap for a single
! switch from DWGs is 0.53” (1.346 cm) gap and 93 psig air X 20 = 26.92 cm
!
! Initial switch resistance was chosen to be higher than that printed in the first
! few time steps so as to appear monotonically decreasing in R plots
! A careful few runs showed that 10 GΩ looked best. In any event by 1 ns all of the 
! resistance values are the same.
!  20 switches for 8 sub-Marxes, 2 channels for 2 modules
!  Divide parallel switch inductance for one Marx (~50 nH) by 2 for 2 modules
! Not right - switch inductance is already included in inductance above.
!

RLSeries 10e9 25e-9
Var R2 Switch
!dielectric switchtime gap        pressure     nswitch  nchannels
AIR           0.0               0.2692  7.3                8             2

TXT IR2 SCAle 1.0e-06
$I_marx_B (MA)
TXT R2
$R_switch_B (Ω)
TXT ER2 SCAle 1.0e-03
$E_switch_B (kJ)
TXT QR2
$Q_switch_B (C)

!
! Output Marx header inductance
! Divide 100 nH per interface by 2 for 2 Marxes in parallel per module
! QE will have two modules per level (divide inductance by 2) 
!

RLSeries 0.00 25e-9

!
! Line 1 (CPL), 60 ns, 1.4 Ω, 60 ns/1.4= 42.85 nF per L1
! Divide Z by 2 for one half of DE-R driving one insulator stack
! Include the radius transition of the L1 
!

TRLine EXP 2.29e-09  2  .7
TRLine LIN 56.0e-09   .7 .7
TRLine EXP 2.29e-09 .7  2

TXT VOUT SCAle 1.0e-06
$V_L1_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_L1_out_B (MA)
TXT EOUT SCAle 1.0e-03
$E_L1_out_B (kJ)
TXT POUT SCAle 1.0e-12
$P_L1_out_B (TW)

!
! Line 1 water switches - gap is scaled by 0.6
! 5 output switches, 4” gap x .6 = 2.4"  (DE today) + switch inductance
!  Increase number of channels x2 for the two modules
!

RLSeries 10e9 5e-9
Var R2 Switch
!dielectric switchtime gap         pressure   nswitch  nchannels
H2O          310e-09     0.06096   1.0              5             2

TXT IR2 SCAle 1.0e-06
$I_L1_sw_B (MA)
TXT R2 (Ω)
$R_L1_sw_B
TXT ER2 SCAle 1.0e-03
$E_L1_sw_B (kJ)
TXT QR2
$Q_L1_sw_B (C)

!
! Line 2, 72 ns, 1.4 Ω constant impedance
! Divide by 2 for 2 modules
!

TRLine EXP 2.29e-09 2 .7

TXT VIN SCAle 1.0e-06
$V_PFL_in_B (MV)
TXT IIN SCAle 1.0e-06
$I_PFL_in_B (MA)

TRLine LIN 68e-09      .7 .7
TRLine EXP 2.29e-09 .7 2

TXT VOUT SCAle 1.0e-06
$V_PFL_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_PFL_out_B (MA)
TXT POUT SCAle 1.0e-12
$P_PFL_out_B (TW)
TXT EOUT SCAle 1.0e-03
$E_PFL_out_B (kJ)

!
! Line 2 water switches - gap is scaled by 0.6
! 7 output switches, 0.5" x .6 = 0.3” gap + switch inductance
! Increase number of channels to two for two modules
!

RLSeries 10e9 5e-9
Var R2 Switch
!dielectric switchtime gap      pressure     nswitch  nchannels
H2O         385e-09    0.00762   1.0                    7               2

TXT IR2 SCAle 1.0e-06
$I_PFL_sw_B (MA)
TXT R2
$R_PFL_sw_B (Ω)
TXT ER2 SCAle 1.0e-03
$E_PFL_sw_B kJ)

!
! Output Line (OL),  72 ns, 1.4 Ω constant impedance
!  Divide Z by 2 for two modules
!

TRLine EXP 2.29e-09 2 .7

TXT VIN SCAle 1.0e-06
$V_OL_in_B (MV)
TXT IIN SCAle 1.0e-06
$I_OL_in_B (MA)

TRLine LIN 70e-09 .7 .7

TXT VOUT SCAle 1.0e-06
$V_OL_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_OL_out_B (MA)
TXT POUT SCAle 1.0e-12
$P_OL_out_B (TW)
TXT EOUT SCAle 1.0e-03
$E_OL_out_B (kJ)

!
!  ********* Water flare transition external to the insulator stack ****************
!  Router = 1.165 m , stack outer = 1.1 m, constant 13.35-cm gap, l = 1.96 ns
!  —— Divide by 2X for half QE Zouter = 0.76 Ω, Zinner = 0.81 Ω
!

TRLine LIN 1.96e-9 0.76 0.81

!
!   ***********************  Insulator stack  *****************************
! The water/vacuum interface Router=1.1 m, Rinner=1.0 m, height=13.35 cm (0.64-cm grading rings) 
!  Plastic height is 11.43 cm. L per side = 2.1788 nH
!

RLSeries  0.0    2.1788E-9

!
TXT VOUT SCAle 1.0e-06
$V_stack_B (MV)
TXT IOUT SCAle 1.0e-06
$I_stack_B (MV)
TXT POUT SCAle 1.0e-12
$P_stack_B (TW)
TXT EOUT SCAle 1.0e-03
$E_stack_B (kJ)

!
!   **************************  First vacuum piece inside insulator  *****************************
!  Constant 13.35 cm height, Router = 1 m, Rinner = 0.975 m
! Inductance = 0.676 nH
!

RLSeries  0.0    0.676E-9

!
!
!   **************************  Vacuum flares  *****************************
! The flare includes the flare itself and the rectangular piece that extends down to the cathode.
! L1=0.639 nH, L2=0.924 nH, L1+L2=1.563 nH
! This section would be better modeled as MITLs as the lower portion of the feed could be 
!  emissive.
!

RLSeries  0.0    1.563E-9

!
! Extra stray L in the flare region 0.5 nH, just a WAG
!

RLSeries  0.0    0.5E-9

!
TXT VOUT SCAle 1.0e-06
$V_vf_B (MV)
TXT IOUT SCAle 1.0e-06
$I_vf_B (MA)
TXT POUT SCAle 1.0e-12
$P_vf_B TW)
TXT EOUT SCAle 1.0e-03
$E_vf_B (kJ)

!
!   *****************************   Constant Z MITLs   ********************************
! Z= 3.7 Ω, length = 2.5 ns (Router = 89.65 cm gap=5.53 cm, Rinner=16.2 cm gap=1 cm)
! L=Zt - 9.25 nH
! I will divide the spacing between Router and Rinner into 10 segments of equal distance
!  and use the mid point of each segment as the circumference for the MITL model.
! 7.345 cm long segments,
! MITL midpoints (cm) - 85.98, 78.63, 71.29, 63.94, 56.60, 49.25, 41.91, 34.56, 27.22, 19.87
! circumferences (m) - 5.40, 4.94, 4.48, 4.017, 3.556, 3.094, 2.633, 2.171, 1.71, 1.248
! vac gap for 3.7 Ω (m) - 5.306e-2, 4.852e-2, 4.4e-2, 3.946e-2, 3.492e-2, 3.039e-2, 2.586e-2, 2.133e-2, 1.68e-2, 1.226e-2
!
!MITL      Circum            Gap            Length-s       Impedance        Resolution       E-Turnon
MITL         5.40            5.306e-02     0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl1_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl1_out_B (MA)
TXT CLOS
$I_mitl1_loss_B (A)
TXT ALOS
$J_mitl1_loss_B (A/m2)

!MZFlow 3.7
!TXT ICA
!$Icath_mitl1_B
!TXT IPL
!$Iflow_mitl1_B
!TXT ZOT
!$Zflow_mitl1_B
               
MITL         4.94          4.852e-02        0.245E-09          3.7                      0.05e-9               

TXT VOUT SCAle 1.0e-06
$V_mitl2_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl2_out_B (MA)
TXT CLOS
$I_mitl2_loss_B (A)
TXT ALOS
$J_mitl2_loss_B (A/m2)

MITL         4.48          4.4e-02        0.245E-09          3.7                      0.05e-9               

TXT VOUT SCAle 1.0e-06
$V_mitl3_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl3_out_B (MA)
TXT CLOS
$I_mitl3_loss_B (A)
TXT ALOS
$J_mitl3_loss_B (A/m2)

MITL         4.017        3.946e-02       0.245E-09          3.7                      0.05e-9               

TXT VOUT SCAle 1.0e-06
$V_mitl4_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl4_out_B (MA)
TXT CLOS
$I_mitl4_loss_B (A)
TXT ALOS
$J_mitl4_loss_B (A/m2)

MITL         3.556         3.492e-02        0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl5_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl5_out_B (MA)
TXT CLOS
$I_mitl5_loss_B (A)
TXT ALOS
$J_mitl5_loss_B (A/m2)

MITL         3.094         3.039e-02        0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl6_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl6_out_B (MA)
TXT CLOS
$I_mitl6_loss_B (A)
TXT ALOS
$J_mitl6_loss_B (A/m2)

MITL         2.633         2.586e-02        0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl7_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl7_out_B (MA)
TXT CLOS
$I_mitl7_loss_B (A)
TXT ALOS
$J_mitl7_loss_B (A/m2)

MITL        2.171         2.133e-02        0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl8_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl8_out_B (MA)
TXT CLOS
$I_mitl8_loss_B (A)
TXT ALOS
$J_mitl8_loss_B (A/m2)

MITL         1.71         1.678e-02        0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl9_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl9_out_B (MA)
TXT CLOS
$I_mitl9_loss_B (A)
TXT ALOS
$J_mitl9_loss_B (A/m2)

MITL         1.248         1.226e-02        0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl10_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl10_out_B (MA)
TXT CLOS
$I_mitl10_loss_B (A)
TXT ALOS
$J_mitl10_loss_B (A/m2)
TXT POUT SCAle 1.0e-12
$P_mitl10_out_B (TW)
TXT EOUT SCAle 1.0e-03
$E_mitl10_out_B (kJ)

!
! Inductance of the disk - plate transition at 16.2 cm to 10.4 cm
!    Gap = 1 cm at this location, top and bottom 0.886 nH
!    Will need MITL calls here at some point.
!

RLSeries 0.0       0.886E-9

!
! Stray L in the PHC region 0.5 nH - just a WAG
! Should be Router = 10 cm, Rinner 6.5 cm, 1-cm gap, L = .862 nH
! This could be a MITL piece but I am treating this in the Zflow losses
!

RLSeries  1.0E-12    0.5E-9

TXT VOUT SCAle 1.0e-06
$V_disk_out_B (MV)
TXT IOUT SCAle 1.0e-06
$I_MITL_B (MA)

!
!   *****************************  Convolute  ******************************
! Convolute diameter 6 inches = 15.24 cm diameter, 7.62 cm radius
!  Assume a 1-cm inner gap/2-cm outer gap for all convoute spacings, Post diameter - 1.59 cm
!  Convolute height = 1 cm + 1.25 cm + 1 cm + 1.25 cm + 1 cm + 1.25 cm + 1 cm = 7.75 cm
!  Divide by 12 posts = 1.37 nH (1.5 cm gap coax estimate)
! The Z51 inductance for B level was 1.89 nH. A 3-D E&M code would give the correct information
!
! In series only with B level

RLSeries 0.0       1.37E-9

!
! Attach A level
! Branch #2
!

ENDbranch

!   *****  Z-flow current loss immediately downstream of the convolute  ******
!
!   We assume the MITL system has a Z-flow impedance of 1.20 ohms for two levels:
!  The 1.2 Ω here seems to be arbitrary. For shot Z51, the value was adusted to match data.
!  The losses could be half that of Z (single PHC vs. double PHC), in which case the
!   Zflow would be 0.6 Ω
!

RCG     0.0001      1.00E-12
VARIABLE R1      POS-MODEL 
!  TSW   CURSW   TOPEN  ZFLOW  GSWMIN  GSWMAX  CBFLAG
   2E-9   1E2            2E-9      1.20        0.0001      10000         0

TXT IR1 SCAle 1.0e-06
$I_zfloss (MA)
TXT R1
$R_zfloss (Ω)

!
!   *******************  Inner-MITL, coax, and pinch at t = 0  ********************
!   THE INNER MITL INDUCTANCE Disk from PHC
!   Router = 7.62 - .625 -1 = 6 cm, Rinner = 2.5 cm, gap = 0.5 cm
! We are treating this as a lumped L but it could be a MITL.
!

RLSeries   0.0       1.04E-9

!
! Extra stray L in the inner MITL region - just a WAG
!

RLSeries  0.0      0.25E-9

TXT VOUT SCAle 1.0e-06
$V_in_coax (MV)

!
!   The lower coax INDUCTANCE Router = 2.5 cm, Rinner = 2 cm, height = 0.5 + 1.25 cm
!      this is to the base of the load
!

RLSeries 0.0       0.781E-9

TXT VOUT SCAle 1.0e-06
$V_out_coax (MV)

!
!   Initial PRS Load Inductance 2 cm height, 2 cm wire array radius, 0.5 cm AK gap
!

RLSeries 0.0  0.89E-9

!     
!    *********   Time-dependent z-pinch model   ****************************
!
!   We choose a 2-cm initial radius, a 2-cm length, 1-mg mass, and a 20:1 convergence ratio.
! 
!                INITIAL R     LENGTH   TOTAL MASS    FINAL R
CYLFOIL      0.02             0.02             1.0e-6            0.001
!
!

TXT FRAD SCAle 1.0e-06
$Rad_PRS (TW)
TXT FVEL
$Vel_PRS (m/s)
TXT FKE SCAle 1.0e-03
$E_kinetic_PRS (kJ)
TXT VIN SCAle 1.0e-06
$V_load (MV)
TXT IIN SCAle 1.0e-06
$I_load (MA)
TXT EIN SCAle 1.0e-03
$E_load (kJ)
TXT PIN SCAle 1.0e-12
$P_load (TW)
TXT L2 SCAle 1.0e+09
$L_load (nH)
TXT R2
$R-liner (Ω)
TXT VR2 SCAle 1.0e-06
$V_liner (MV)
TXT PR2 SCAle 1.0e-12
$P_liner (TW)

!
!   ********   Tie it all back to ground   ********************************
!

RCG  1.0E-12  1.0E-12

!
! End Main Branch (Branch #1)
!
! Level 2 Branches
!
! Branch #2
! Note: all elements are in reverse order inside to outside
!

Branch

!
! L through the PHC region to the post center 0.5 nH -  just a WAG now
! Should be Router = 10 cm, Rinner 6.5 cm, 1-cm gap, L = .862 nH
! This could be a MITL piece but I am treating this in the Zflow losses
!

RLSeries  0.0    0.5E-9

TXT IIN SCAle 1.0e-06
$I_mitl_A (MA)

!
! Inductance of the disk - plate transition at 16.2 cm to 10.4 cm
!    Gap = 1 cm at this location, top and bottom 0.886 nH
! Need to replace this with several MITL elements.
!

RLSeries 0.0       0.886E-9


!   *****************************   Constant Z MITLs   ********************************
! Z= 3.7 Ω, length = 2.5 ns (Router = 89.65 cm gap=5.53 cm, Rinner=16.2 cm gap=1 cm)
! L=Zt - 9.25 nH
! I will divide the spacing between Router and Rinner into 10 segments of equal distance
!  and use the mid point of each segment as the circumference for the MITL model.
! 7.345 cm long segments,
! MITL midpoints (cm) - 85.98, 78.63, 71.29, 63.94, 56.60, 49.25, 41.91, 34.56, 27.22, 19.87
! circumferences (m) - 5.40, 4.94, 4.48, 4.017, 3.556, 3.094, 2.633, 2.171, 1.71, 1.248
! vac gap for 3.7 Ω (m) - 5.306e-2, 4.852e-2, 4.4e-2, 3.946e-2, 3.492e-2, 3.039e-2, 2.586e-2, 2.133e-2, 1.68e-2, 1.226e-2
!
! Remember in the branch outputs and inputs are in reverse order
!
!MITL      Circum            Gap            Length-s       Impedance        Resolution       E-Turnon
MITL         1.248         1.226e-02        0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl10_out_A (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl10_out_A (MA)
TXT CLOS
$I_mitl10_loss_A (A)
TXT ALOS
$J_mitl10_loss_A (A/m2)
TXT POUT SCAle 1.0e-12
$P_mitl10_out_A (TW)
TXT EOUT SCAle 1.0e-03
$E_mitl10_out_A (kJ)

MITL         1.71         1.678e-02        0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl9_out_A (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl9_out_A (MA)
TXT CLOS
$I_mitl9_loss_A (A)
TXT ALOS
$J_mitl9_loss_A (A/m2)

MITL        2.171         2.133e-02        0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl8_out_A (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl8_out_A (MA)
TXT CLOS
$I_mitl8_loss_A (A)
TXT ALOS
$J_mitl8_loss_A (A/m2)

MITL         2.633         2.586e-02        0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl7_out_A (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl7_out_A (MA)
TXT CLOS
$I_mitl7_loss_A (A)
TXT ALOS
$J_mitl7_loss_A (A/m2)

MITL         3.094         3.039e-02        0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl6_out_A (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl6_out_A (MA)
TXT CLOS
$I_mitl6_loss_A (A)
TXT ALOS
$J_mitl6_loss_A (A/m2)

MITL         3.556         3.492e-02        0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl5_out_A (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl5_out_A (MA)
TXT CLOS
$I_mitl5_loss_A (A)
TXT ALOS
$J_mitl5_loss_A (A/m2)

MITL         4.017        3.946e-02       0.245E-09          3.7                      0.05e-9               

TXT VOUT SCAle 1.0e-06
$V_mitl4_out_A (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl4_out_A (MA)
TXT CLOS
$I_mitl4_loss_A (A)
TXT ALOS
$J_mitl4_loss_A (A/m2)

MITL         4.48          4.4e-02        0.245E-09          3.7                      0.05e-9               

TXT VOUT SCAle 1.0e-06
$V_mitl3_out_A (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl3_out_A (MA)
TXT CLOS
$I_mitl3_loss_A (A)
TXT ALOS
$J_mitl3_loss_A (A/m2)

MITL         4.94          4.852e-02        0.245E-09          3.7                      0.05e-9               

TXT VOUT SCAle 1.0e-06
$V_mitl2_out_A (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl2_out_A (MA)
TXT CLOS
$I_mitl2_loss_A (A)
TXT ALOS
$J_mitl2_loss_A (A/m2)

MITL         5.40            5.306e-02     0.245E-09          3.7                      0.05e-9

TXT VOUT SCAle 1.0e-06
$V_mitl1_out_A (MV)
TXT IOUT SCAle 1.0e-06
$I_mitl1_out_A (MA)
TXT CLOS
$I_mitl1_loss_A (A)
TXT ALOS
$J_mitl1_loss_A (A/m2)
               
!
! Stray L in the flare region 0.5 nH - divide by two, just a WAG
!

RLSeries  1.0E-12    0.5E-9

TXT VIN SCAle 1.0e-06
$V_vf_A (MV)
TXT IIN SCAle 1.0e-06
$I_vf_A (MA)
TXT PIN SCAle 1.0e-12
$P_vf_A (TW)
TXT EIN SCAle 1.0e-03
$E_vf_A (kJ)

!
!
!   **************************  Vacuum flares  *****************************
! The flare includes the flare itself and the piece that extends down to the cathode.
! L1=0.639 nH, L2=0.924 nH, L1+L2=1.563 nH
! Two vacuum flares in parallel - divide by 2
!

RLSeries  0.0    1.563E-9

!
!   **************************  First vacuum piece inside insulator  *****************************
!  Constant 13.35 cm height, Router = 1 m, Rinner = 0.975 m
! Inductance = 0.676 nH - divide by 2
!

RLSeries  0.0    0.676E-9

!
!   ***********************  Insulator stack  *****************************
! The water/vacuum interface Router=1.1 m, Rinner=1.0 m, height=13.35 cm (0.64-cm grading rings) 
!  Plastic height is 11.43 cm. L per side = 2.1788 nH
! Two stacks in parallel - divide L by 2
!

RLS  1.0E-12    2.1788E-9

TXT VIN SCAle 1.0e-06
$V_stack_A (MV)
TXT IR2 SCAle 1.0e-06
$I_stack_A (MV)
TXT PIN SCAle 1.0e-12
$P_stack_A (TW)
TXT EIN SCAle 1.0e-03
$E_stack_A (kJ)

!
!  ********* Water flare transition external to the insulator stack ****************
!  Router = 1.165 m , stack outer = 1.1 m, constant 13.35-cm gap, l = 1.96 ns
!  Zouter = 0.76 Ω, Zinner = 0.81 Ω 
!

TRLine LIN 1.96e-9 0.81 0.76

!
! Output Line (OL),  72 ns, 1.4 Ω constant impedance
!  Divide by 2 for 2 modules
!

TRLine LIN 70e-09 .7 .7

TXT VIN SCAle 1.0e-06
$V_OL_out_A (MV)
TXT IIN SCAle 1.0e-06
$I_OL_out_A (MA)
TXT PIN SCAle 1.0e-12
$P_OL_out_A (TW)
TXT EIN SCAle 1.0e-03
$E_OL_out_A (kJ)

TRLine EXP 2.29e-09 .7 2

TXT VOUT SCAle 1.0e-06
$V_OL_in_A (MV)
TXT IOUT SCAle 1.0e-06
$I_OL_in_A (MA)

!
! Line 2 water switches - gap is scaled by 0.6
! 7 output switches, 0.5" x .6 = 0.3” gap + switch inductance
! Increase number of channels to two for two modules
!

RLSeries 10e9 5e-9
Var R2 Switch
!dielectric switchtime gap      pressure     nswitch  nchannels
H2O         385e-09    0.00762   1.0                    7               2

TXT IR2 SCAle 1.0e-06
$I_PFL_sw_A (MA)
TXT R2
$R_PFL_sw_A (Ω)
TXT ER2 SCAle 1.0e-03
$E_PFL_sw_A (kJ)
TXT QR2
$Q_PFL_sw_A (C)

!
! Line 2, 72 ns, 1.4 Ω constant impedance
! Divide by 2 for 2 modules
!

TRLine EXP 2.29e-09 2 .7

TXT VIN SCAle 1.0e-06
$V_PFL_out_A (MV)
TXT IIN SCAle 1.0e-06
$I_PFL_out_A (MA)
TXT PIN SCAle 1.0e-12
$P_PFL_out_A (TW)
TXT EIN SCAle 1.0e-03
$E_PFL_out_A (kJ)

TRLine LIN 68e-09 .7 .7
TRLine EXP 2.29e-09 .7 2

TXT VOUT SCAle 1.0e-06
$V_PFL_in_A (MV)
TXT IOUT SCAle 1.0e-06
$I_PFL_in_A (MA)

!
! CPL water switches
! 5 output switches, 4” gap + switch inductance
!  Increase number of channels for the two modules
!

RLSeries 10e9 5e-9
Var R2 Switch
!dielectric switchtime gap         pressure   nswitch  nchannels
H2O          310e-09     0.06096   1.0              5             2

TXT IR2 SCAle 1.0e-06
$I_L1_sw_A (MA)
TXT R2
$R_L1_sw_A (Ω)
TXT ER2 SCAle 1.0e-03
$E_L1_sw_A (kJ)

!
! Line 1 (CPL), 60 ns, 1.4 Ω, 60 ns/1.4= 42.85 nF per L1
! Divide Z by 2 for one half of DE-R driving one insulator stack
! Include the radius transition of the L1 
!

TRLine EXP 2.29e-09    2  .7

TXT VIN SCAle 1.0e-06
$V_L1_out_A (MV)
TXT IIN SCAle 1.0e-06
$I_L1_out_A (MA)
TXT EIN SCAle 1.0e-03
$E_L1_out_A (kJ)
TXT PIN SCAle 1.0e-12
$P_L1_out_A (TW)

TRLine LIN   56.0e-09   .7  .7
TRLine EXP 2.29e-09    .7   2

!
! Output Marx header inductance
! Divide 100 nH per interface by 2 for 2 Marxes in parallel per module
! QE will have two modules per level (divide inductance by 2) 
!

RLSeries 0.00 25e-9

!
! Switch resistance - Martin Model and switch inductance - the total gap for a single
! switch from DWGs is 0.53” (1.346 cm) gap and 93 psig air X 20 = 26.92 cm
!
! Initial switch resistance was chosen to be higher than that printed in the first
! few time steps so as to appear monotonically decreasing in R plots
! A careful few runs showed that 10 GΩ looked best. In any event by 1 ns all of the 
! resistance values are the same.
!  20 switches for 8 sub-Marxes, 2 channels for 2 modules
!  Divide parallel switch inductance for one Marx (~50 nH) by 2 for 2 modules
! Not right - switch inductance is already included in inductance above.
!

RLSeries 10e9 25e-9
Var R2 Switch
!dielectric switchtime gap        pressure     nswitch  nchannels
AIR           0.0               0.2692  7.3                8             2

TXT IR2 SCAle 1.0e-06
$I_marx_B (MA)
TXT R2
$R_switch_B (Ω)
TXT ER2 SCAle 1.0e-03
$E_switch_B (kJ)
TXT QR2
$Q_switch_B (C)

!
! Cap inductance, case & parasitic inductance, and ESR
! ESR = 0.325 for the total Marx module (2.6 Ω sub Marx/8)
! Inductance is 417 nH per module 
! Divide R and L by 2 for 2 modules (driving one MITL)
!

RLSeries 0.1625 208.5e-9

!
! Marx capacitance and charge voltage
! CX2 for two modules (58 nF each Marx module) = 116 nF
! +/-85 kV charge x 40 = 3.4 MV
!
RCGround 1e+12 116e-9
Initial VC1 3.4e6

TXT VC1 SCAle 1.0e-06
$V_marx_A (MV)
TXT EC1 SCAle 1.0e-03
$E_marx_A (kJ)

!
! End branch #2
