'savefig.dpi': 150, #*12/16
'savefig.jpeg_quality': 100,
'savefig.pad_inches': 0.1,
'lines.linewidth': 0.75
}
sns.set(context='paper',
style='white',
font='sans-serif',
palette='Set2',
font_scale=7 / 8,
rc=rc)
Color = cycle(sns.color_palette('Set2'))
nTF = 18 #,nPF,nCS = 18,6,5
config = {'TF': 'dtt', 'eq': 'SFp'}
config, setup = select(config, nTF=nTF, update=False) # nPF=nPF,nCS=nCS
sf = SF(setup.filename)
rb = RB(setup, sf)
pf = PF(sf.eqdsk)
tf = TF(
Profile(config['TF'], family='S', part='TF', nTF=nTF, obj='L', load=True))
rb.firstwall(symetric=True, DN=True, plot=True)
#pf.plot(coils=pf.coil,label=False,plasma=False,current=True)
tf.fill()
#sf.cpasma *= 1.1
eq = EQ(sf, pf, dCoil=1.5, sigma=0, boundary=sf.get_sep(expand=1.1), n=3e3)
eq.gen_bal()
#eq.gen_opp()
alpha=0.3,color=sns.color_palette('Set2',5)[3])
#self.cage.plot_loops(sticks=False)
if len(ax) > 1:
pl.sca(ax[1])
pl.cla()
plot_oppvar(shp.loop.xo,shp.loop.oppvar)
if __name__ is '__main__':
nTF = 16
family='S'
ripple = False
config = {'TF':'demo','eq':'SN'}
config,setup = select(config,nTF=nTF)
demo = DEMO()
profile = Profile(config['TF'],family=family,part='TF',nTF=nTF) # ,load=False
shp = Shape(profile,nTF=nTF,obj='L',eqconf=config['eq'],ny=1)
shp.add_vessel(demo.parts['Vessel']['out'])
shp.minimise(ripple=ripple,verbose=False)
cage = shp.cage
#shp.update()
#shp.tf.fill()
#shp.loop.plot({'flat':0.3,'tilt':13})
#shp.loop.plot()
#demo.fill_part('TF_Coil',alpha=0.8)
from nova.loops import Profile
from nova.inverse import SWING
from nova.firstwall import main_chamber
from nova.shape import Shape
import seaborn as sns
sns.set(context='talk', style='white', font='sans-serif', font_scale=7 / 8)
Color = cycle(sns.color_palette('Set2'))
pkl = PKL('SXdev', directory='../../Movies/')
name = 'SX'
nTF, ripple = 18, True
base = {'TF': 'dtt', 'eq': 'DTT_SN', 'name': name}
config, setup = select(base, nTF=nTF, update=False)
sf = SF(setup.filename)
pf = PF(sf.eqdsk)
for coil in pf.index['CS']['name']:
pf.coil[coil]['r'] = 2.7
pf.coil[coil]['dr'] = 0.8
profile = Profile(config['TF'], family='S', part='TF', nTF=nTF, obj='L')
tf = TF(profile=profile, sf=sf)
#tf.adjust_xo('dz',value=-2)
#tf.fill()
eq = EQ(sf,
pf,
'savefig.jpeg_quality': 100,
'savefig.pad_inches': 0.1,
'lines.linewidth': 0.75
}
sns.set(context='paper',
style='white',
font='sans-serif',
palette='Set2',
font_scale=7 / 8,
rc=rc)
Color = cycle(sns.color_palette('Set2'))
#base = {'TF':'dtt','eq':'SN'}
base = {'TF': 'dtt', 'eq': 'SN2014_EOF'}
#config,setup = select(base={'TF':'dtt','eq':'SN2014_EOF'},nTF=18)
config, setup = select(base=base, update=False, nTF=18, nPF=5, nCS=3)
sf = SF(setup.filename)
rb = RB(setup, sf)
pf = PF(sf.eqdsk)
tf = TF(
Profile(config['TF'],
family='S',
part='TF',
nTF=config['nTF'],
obj='L',
load=True))
tf.fill()
eq = EQ(sf,
pf,
dCoil=1.0,
sigma=0,
nsel,r,node,,wp ! ensure all nodes from winding pack
*get,nnd,node,0,count ! count nodes
*do,j,1,3 ! Fx,Fy,Fz - all nodes attached to element
F,all,F%xyz(j)%,Fbody_%xyz(j)%(el_l(i),el_dr(i),el_o(i))*el_v(i)/nnd
*enddo
*enddo
allsel
'''
ans.f.write(apdlstr)
ans.close()
if __name__ is '__main__':
nTF, ripple = 13, True
base = {'TF': 'demo_nTF', 'eq': 'DEMO_SN_SOF'}
config, setup = select(base, nTF=nTF, update=True)
occ = OCC(config, setup, nTF=nTF)
occ.OIS()
occ.write()
occ.plot()
#occ.ansys(plot=True,nl=50,nr=3,ny=1)
pl.plot(0, -12, 'o', alpha=0)
pl.plot(17, 10, 'o', alpha=0)
pl.tight_layout()
#pl.savefig('../Figs/CoilForces.png')
rc = {
'figure.figsize': [7, 7 * 16 / 12],
'savefig.dpi': 125, #
'savefig.jpeg_quality': 200,
'savefig.pad_inches': 0.1,
'lines.linewidth': 1.5
}
sns.set(context='talk',
style='white',
font='sans-serif',
palette='Set2',
font_scale=7 / 8,
rc=rc)
color = sns.color_palette('Set2', 5)
config, setup = select(base={'TF': 'dtt', 'eq': 'DEMO_FW_SOF'}, nTF=18)
sf = SF(setup.filename)
pf = PF(sf.eqdsk)
pf.plot(coils=pf.coil, label=True, plasma=False, current=True)
levels = sf.contour(plot_vac=False)
rb = RB(setup, sf)
rb.firstwall(plot=True, debug=False, color=color[1])
rb.trim_sol()
filename = trim_dir('../../Data/') + 'CATIA_FW.xlsx'
wb = load_workbook(filename=filename, read_only=True, data_only=True)
ws = wb[wb.get_sheet_names()[0]]
FW = {}
for col, var in zip([5, 6], ['r', 'z']):
row = ws.columns[col]
from nova.force import force_feild
import json
from amigo.IO import trim_dir
from nova.shelf import PKL
import seaborn as sns
rc = {'figure.figsize':[5,5*16/12],'savefig.dpi':120, #
'savefig.jpeg_quality':200,'savefig.pad_inches':0.1,
'lines.linewidth':1.5}
sns.set(context='talk',style='white',font='sans-serif',palette='Set2',
font_scale=7/8,rc=rc)
#config,setup = select(base={'TF':'dtt','eq':'SN'},nTF=18,nPF=5,nCS=3)
#config,setup = select(base={'TF':'dtt','eq':'SX'},nTF=18,nPF=5,nCS=3)
#config,setup = select(base={'TF':'dtt','eq':'DEMO_FW_SOF'},nTF=18)
config,setup = select(base={'TF':'dtt','eq':'SFp'},nTF=18)
sf = SF(setup.filename)
pf = PF(sf.eqdsk)
'''
eq = EQ(sf,pf,dCoil=1.5,sigma=0,boundary=sf.get_sep(expand=1.05),
zmin=-8,rmin=5.5,n=1e4)
eq.gen_opp()
'''
#pf.plot(coils=pf.coil,label=True,plasma=False,current=True)
sf.contour()
rb = RB(setup,sf)
