sheild = [0.85, 0.1] # neutron sheilding
Icoil = 17.5e6 # coil amp-turns
Jmax = 15e6 # max current density
Acoil = Icoil / Jmax
dRcoil = np.sqrt(Acoil)
dRsteel = 0.15 * dRcoil
flip = 1
graze = 1.5 * np.pi / 180
with open('./plot_data/' + config + '_sol.pkl', 'rb') as input:
sf = pickle.load(input)
plot = pickle.load(input)
geom = pickle.load(input)
sol = solCalc(sf, geom, config, flip=flip)
rb = RB(geom, sf, config, Np=400, flip=flip)
#rb.internal_coils(['P6','P5','PS2','P6B','P5B','PS2','PS5'])
pl.figure(figsize=(18, 12))
pl.axis('equal')
pl.xlim([-15.5, 15.5]), pl.ylim([-15, 9])
#pl.xlim([6,9]),pl.ylim([-5,-3])
#pl.xlim([6,14]),pl.ylim([-10,-6])
#pl.xlim([2,17.5]),pl.ylim([-11,5])
rb.fill(dR=gap, dt=tfw, alpha=0.7, color='r', label='first wall')
rb.fill(dt=BB, alpha=0.7, color='b', label='blanket')
rb.fill(dt=tBBsupport, alpha=0.4, color='g', label='blanket support')
P6support = {}
import pickle
from radial_build import RB
from feild_calc import solCalc
import matplotlib
font = {'family': 'serif', 'serif': 'Arial', 'weight': 'normal', 'size': 16}
matplotlib.rc('font', **font)
fig_path = '../Figs/'
config = 'DD3' # DD1/DD3
with open('./plot_data/' + config + '_sol.pkl', 'rb') as input:
sf = pickle.load(input)
plot = pickle.load(input)
geom = pickle.load(input)
sol = solCalc(sf, geom, config)
rb = RB(geom, sf, config, Np=400)
divertor_coils = ['P6', 'P5', 'PS2', 'P6B', 'P5B', 'PS2', 'PS5']
for name in divertor_coils:
geom.coil[name]['Fr'], geom.coil[name]['Fz'] = Fcoil(geom.coil, name)
fig = pl.figure(figsize=(10, 6))
plot = layout(geom)
plot.coil_keys = divertor_coils
#plot.sol() # scrape-off layer
#plot.plates() # divertor plates
#plot.TF() # TF coil
#plot.FW()
#plot.plasma()
import pickle
import pylab as pl
import numpy as np
from feild_calc import solCalc
with open('./plot_data/DD3_sol.pkl', 'rb') as input:
sf = pickle.load(input)
plot = pickle.load(input)
geom = pickle.load(input)
sol = solCalc(sf, geom)
fig = pl.figure(figsize=(14, 12))
pl.axis('equal')
pl.grid()
pl.xlim([4.5, 13.5])
pl.ylim([-10, -5])
#plot.sol() # scrape-off layer
plot.plates() # divertor plates
plot.TF() # TF coil
plot.FW()
plot.plasma()
plot.coil_fill()
#plot.coil_current([0,20],[-11,-2])
plot.coil_sheild()
#sf.sol()
p3arc = [[10, -9.53], [12, -7.8], [13, -5.8]] # three-point strike arc
sol.strike_arc(p3arc, 0.25, plot=True)
for leg in ['inner', 'outer']: