fix['value'] = np.append(fix['value'],0) fix['BC'] = np.append(fix['BC'],'Bz') ''' del sf.coil['Coil9'] #del sf.coil['Coil10'] #del sf.coil['Coil11'] #del sf.coil['Coil12'] #del sf.coil['Coil13'] eq.coils(delta=0.25) # multi-filiment coils conf.TF(sf) rb.TFopp(False,objF=conf.TFopp) # L==length, V==volume rb.TFfill() eq.coil_currents(fix) sf.plot_coils(Color,coils=sf.coil,label=False,plasma=False,current=True) sf.plot_coils(Color,coils=eq.coil,label=False,plasma=False) eq.update_psi() # levels=sf.cs.levels sf.contour(Nstd=1.5) #sf.sol(plot=True) eq = EQ([3,14],[-10,10],8e3,sf) #eq = EQ([2,16],[-10,10],2e3,sf)
sf.plot_coils(Color) conf.TF(sf) rb = RB(conf, sf, Np=200) rb.divertor_outline(True) rb.fill(dt=conf.tfw, alpha=0.7, color=next(Color)) rb.fill(dt=conf.BB, alpha=0.7, color=next(Color)) rb.fill(dt=conf.tBBsupport, alpha=0.7, color=next(Color)) rb.BBsheild_fill(dt=conf.sheild, ref_o=2 / 8 * np.pi, alpha=0.7, color=next(Color)) rb.fill(dt=conf.VV, alpha=0.7, color=next(Color), loop=True) rb.set_TFbound() # TF boundary conditions #rb.TFbound['ro_min'] -= 0.25 #rb.plot_TFbounds() rb.TFopp(False, objF='L') # L==length, V==volume rb.TFfill() pl.text(16, 9.5, config, fontsize=36, ha='right', va='top') ''' Color = cycle(sns.color_palette('Set2')) for leg in conf.targets.keys(): R,Z = rb.sol.legs(leg) pl.plot(R,Z,color=next(Color)) pl.tight_layout() col_labels=['volume','length'] row_labels=['TF','LCFS','ratio'] table_vals=[[r'{:1.0f}m$^3$'.format(rb.TFvol),r'{:1.1f}m'.format(rb.TFlength)], [r'{:1.0f}m$^3$'.format(rb.Pvol),r'{:1.1f}m'.format(rb.Plength)], ['{:1.2f}'.format(rb.Rvol),r'{:1.2f}'.format(rb.Rlength)]]
P6support['bottom'] = {} P5support = {} P5support['top'] = {} P5support['bottom'] = {} if config is 'DD3': rb.VVfill(connect=[0.36, 0.802], r=[13.56, 13.5, 10, 6.8, 5.2], z=[-1.5, -5, -10, -9.5, -3.5], w=[0.9, 2, 2, 3.5, 2], dt=VV, ref_o=3 / 8 * np.pi, alpha=0.3, color='b', label='vacuum vessel') rb.TFopp(TFopp) rb.rzGet() # VV ports rb.fill(trim=[0.795, 0.795 + 0.04], dt=VV, ref_o=3 / 8 * np.pi, alpha=1, color='r') rb.rzGet() rb.fill(trim=[0.745, 0.745 + 0.04], dt=VV, ref_o=3 / 8 * np.pi, alpha=1, color='r', label='ports')
r95,z95 = sf.get_boundary(alpha=0.95) ru = r95[np.argmax(z95)] # triangularity rl = r95[np.argmin(z95)] del_u = (R-ru)/a del_l = (R-rl)/a kappa = (np.max(z95)-np.min(z95))/(2*a) print('R',R,'a',a,'AR',AR,'del_u',del_u,'del_l',del_l,'kappa',kappa) Pvol = rb.loop_vol(rb.Rp,rb.Zp,plot=False) print(Pvol) ''' rb.TFopp(False,config='SXex') # load TF outline rb.TFfill() # construct TF geometory ''' conf.TFopp = 'L' rb.set_TFbound() # TF boundary conditions #rb.TFbound['ro_min'] -= 0.5 rb.plot_TFbounds() rb.TFopp(True,objF=conf.TFopp) # L==length, V==volume rb.TFfill() # construct TF geometory ''' fig.tight_layout(rect=[-0.3,-0.54,1.1,1.5]) pl.savefig('../Figs/SXcoils_const.png',dpi=300)