added data analysis
This commit is contained in:
parent
dd2a0f93b3
commit
1456d0f0c4
@ -19,9 +19,9 @@ def Q_heat(flow,d_T):
|
||||
def main():
|
||||
|
||||
|
||||
d_T = 1.5
|
||||
d_T = 2
|
||||
|
||||
flow = 0.3/5#/5 #m/s
|
||||
flow = 1#/5 #m/s
|
||||
#flow *=2
|
||||
print(f"flow = {flow}m/s")
|
||||
V_t = 4.8 * 3.2 * 1e-6 * flow
|
||||
|
||||
64
Data_Coils/Data.py
Normal file
64
Data_Coils/Data.py
Normal file
@ -0,0 +1,64 @@
|
||||
# %%
|
||||
import matplotlib.pyplot as plt
|
||||
import RigolWFM.wfm as rigol
|
||||
|
||||
# %%
|
||||
|
||||
hhpion = 'C:/Users/Joschka/Desktop/Coil_Data/New/hhpion.wfm'
|
||||
hhpioff = 'C:/Users/Joschka/Desktop/Coil_Data/New/hhpioff.wfm'
|
||||
hhon = 'C:/Users/Joschka/Desktop/Coil_Data/New/hhon.wfm'
|
||||
hhoff = 'C:/Users/Joschka/Desktop/Coil_Data/New/hhoff.wfm'
|
||||
|
||||
scope = 'DS1104Z-S'
|
||||
|
||||
# %%
|
||||
hhpion = rigol.Wfm.from_file(hhpion, scope)
|
||||
hhpion.plot()
|
||||
plt.show()
|
||||
# %%
|
||||
hhpioff = rigol.Wfm.from_file(hhpioff, scope)
|
||||
hhpioff.plot()
|
||||
plt.show()
|
||||
|
||||
# %%
|
||||
hhoff = rigol.Wfm.from_file(hhoff, scope)
|
||||
hhoff.plot()
|
||||
plt.show()
|
||||
|
||||
# %%
|
||||
hhon = rigol.Wfm.from_file(hhon, scope)
|
||||
hhon.plot()
|
||||
plt.show()
|
||||
|
||||
|
||||
# %%
|
||||
print(w.channels[0].times)
|
||||
print(w.channels[0].volts)
|
||||
|
||||
# %%
|
||||
|
||||
ahhpion = 'C:/Users/Joschka/Desktop/Coil_Data/New/ahhpion.wfm'
|
||||
ahhpioff = 'C:/Users/Joschka/Desktop/Coil_Data/New/ahhpioff.wfm'
|
||||
ahhon = 'C:/Users/Joschka/Desktop/Coil_Data/New/ahhon.wfm'
|
||||
ahhoff = 'C:/Users/Joschka/Desktop/Coil_Data/New/ahhoff.wfm'
|
||||
|
||||
scope = 'DS1104Z-S'
|
||||
|
||||
# %%
|
||||
ahhpion = rigol.Wfm.from_file(ahhpion, scope)
|
||||
ahhpion.plot()
|
||||
plt.show()
|
||||
# %%
|
||||
ahhpioff = rigol.Wfm.from_file(ahhpioff, scope)
|
||||
ahhpioff.plot()
|
||||
plt.show()
|
||||
|
||||
# %%
|
||||
ahhoff = rigol.Wfm.from_file(ahhoff, scope)
|
||||
ahhoff.plot()
|
||||
plt.show()
|
||||
|
||||
# %%
|
||||
ahhon = rigol.Wfm.from_file(ahhon, scope)
|
||||
ahhon.plot()
|
||||
plt.show()
|
||||
181
Test_class.py
181
Test_class.py
@ -10,177 +10,16 @@ from src import B_field_calculation as bf
|
||||
|
||||
from src import coil_class as BC
|
||||
|
||||
from IPython import get_ipython
|
||||
get_ipython().run_line_magic('matplotlib', 'qt')
|
||||
#get_ipython().run_line_magic('matplotlib', 'inline')
|
||||
|
||||
#set up axis
|
||||
x_m = np.linspace(-0.05, 0.05, 101)
|
||||
z_m = np.linspace(-0.05, 0.05, 101)
|
||||
# %%
|
||||
|
||||
res = BC.BCoil.resistivity_copper(20)
|
||||
res2 = BC.BCoil.resistivity_copper(50)
|
||||
l = 64 * 2*np.pi *48e-3
|
||||
A = (0.25e-3)**2 * np.pi
|
||||
R1 = res * l/A
|
||||
print(R1)
|
||||
R2 = res2 *l/A
|
||||
print(R2)
|
||||
|
||||
z = z_m*1e3
|
||||
x = x_m*1e3 #for plotting in mm
|
||||
|
||||
#Import Values from simulation
|
||||
|
||||
################# My simulation #########################
|
||||
I = 5
|
||||
HH = 1
|
||||
d_coils = 44
|
||||
R_mid = 44
|
||||
R_inner = 44-3*1.7
|
||||
|
||||
layers = 6
|
||||
windings = 2
|
||||
wire_width = 1.7
|
||||
wire_height = 2.6
|
||||
|
||||
|
||||
HH_Coil1 = BC.BCoil(HH, d_coils ,R_mid, layers, windings, wire_width, wire_height)
|
||||
#HH_Coil1.print_info()
|
||||
#B_z_sim, B_x_sim = HH_Coil1.B_field(5, x, z)
|
||||
|
||||
#Bz, B_x = bf.B_multiple_raster(I,HH,R_inner,d_coils,layers,windings,wire_width, wire_height, x_m,z_m)
|
||||
|
||||
#B_test = B_field_ideal_AHH(layers*windings,I,R_inner*1e-3,d_coils*1e-3,z_m)
|
||||
|
||||
#B_x = np.concatenate((-np.flip(B_r),B_r[1:len(B_r)]))
|
||||
|
||||
|
||||
HH_Coil1.B_quick_plot(I, x, z)
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
#Calculate gradients/curvature
|
||||
B_z_sim_grad = np.gradient(np.gradient(B_z_sim,z_m),z_m)/1e4
|
||||
B_x_sim_grad = np.gradient(B_x_sim,x_m)/100
|
||||
#B_z_grad = np.gradient(np.gradient(Bz,z_m),z_m)/1e4
|
||||
B_z_grad = np.gradient(B_z,z_m)/100
|
||||
|
||||
B_z_sim_grad = np.gradient(B_z_grad,z_m)/100
|
||||
|
||||
B_x_grad = np.gradient(B_x,x_m)/100
|
||||
|
||||
|
||||
#Calculate relative differences in permille
|
||||
rel_diff_Bz = (B_z-B_z_sim)/np.mean(B_z)
|
||||
#rel_diff_Bx = (B_x-B_x_sim)/np.mean(B_x)
|
||||
rel_diff_Bz_grad = (B_z_grad-B_z_sim_grad)/np.mean(B_z_grad)
|
||||
|
||||
rel_diff_Bz_grad_mean = (B_z_grad-B_z_sim_grad)/np.mean(B_z_grad)
|
||||
#rel_diff_Bx_grad = (B_x_grad-B_x_sim_grad)/np.mean(B_x_grad)
|
||||
|
||||
#Plotting
|
||||
plt.figure(1,figsize=(20,18))
|
||||
|
||||
plt.rcParams.update({'font.size': 15})
|
||||
plt.suptitle("Helmholtz coil field Bz along z-axis, comparison of simulations", fontsize=30)
|
||||
|
||||
|
||||
#Field plot
|
||||
##########################
|
||||
plt.subplot(3,2,1)
|
||||
plt.plot(z,B_z,linestyle = "solid", label = r"$Bz$: Result via elliptic integrals")
|
||||
plt.plot(z,B_z_sim,linestyle = "dashdot", label = r"$B_{z, sim}$: Numerical Matlab simulation")
|
||||
plt.plot(z,(B_z-B_z_sim), label = r"$Bz - B_{z, sim}$")
|
||||
#plt.xlim(-0.01,0.01)
|
||||
plt.title("B-field" ,fontsize = 30)
|
||||
|
||||
plt.ylabel(r"$Bz$ [G]")
|
||||
plt.xlabel("z-axis [mm]")
|
||||
plt.legend()
|
||||
|
||||
|
||||
#############################
|
||||
plt.subplot(3,2,3)
|
||||
plt.plot(z,(B_z-B_z_sim), label = r"$Bz - B_{z, sim}$")
|
||||
plt.ylabel("absolute deviation [G]")
|
||||
plt.xlabel("z-axis [mm]")
|
||||
plt.legend()
|
||||
|
||||
#############################
|
||||
plt.subplot(3,2,5)
|
||||
plt.plot(z,1000*rel_diff_Bz, label = "$(Bz - B_{z, sim}) / Bz$")
|
||||
plt.ylabel("relative deviation [‰]")
|
||||
plt.xlabel("z-axis [mm]")
|
||||
plt.legend()
|
||||
|
||||
|
||||
######################Gradient plot############################
|
||||
|
||||
################
|
||||
plt.subplot(3,2,2)
|
||||
plt.plot(z,B_z_grad,linestyle = "solid", label = r"$\nabla_z^2 Bz$: Result via elliptic integrals")
|
||||
plt.plot(z,B_z_sim_grad,linestyle = "dashdot", label = r"$\nabla_z^2 B_{z, sim}$: Numerical Matlab sim.")
|
||||
plt.plot(z,(B_z_grad-B_z_sim_grad), label = r"$\nabla_z^2 Bz - \nabla_z^2 B_{z, sim}$")
|
||||
|
||||
plt.ylabel(r"$\nabla_z^2 Bz [G/cm^2]$")
|
||||
plt.xlabel("z-axis [mm]")
|
||||
plt.title("Curvature of B-field",fontsize = 30)
|
||||
plt.legend(loc='lower right')
|
||||
|
||||
|
||||
#################
|
||||
|
||||
plt.subplot(3,2,4)
|
||||
plt.plot(z,(B_z_grad-B_z_sim_grad), label = r"$\nabla_z^2 Bz - \nabla_z^2 B_{z, sim}$")
|
||||
plt.ylabel(r"absolute deviation $[G/cm^2]$")
|
||||
plt.xlabel("z-axis [mm]")
|
||||
plt.legend()
|
||||
|
||||
#####################
|
||||
plt.subplot(3,2,6)
|
||||
plt.plot(z,1000*rel_diff_Bz_grad, label = r"$(\nabla_z^2 Bz - \nabla_z^2 B_{z, sim}) / \nabla_z^2 Bz$")
|
||||
#plt.ylim(-57,10)
|
||||
plt.ylabel("relative deviation [‰]")
|
||||
plt.xlabel("z-axis [mm]")
|
||||
plt.legend()
|
||||
|
||||
|
||||
plt.savefig("output/HH_benchmark_5A_6x2.pdf")
|
||||
plt.show()
|
||||
|
||||
############### relative deviation with averaging by the mean not the individual value ########################################
|
||||
plt.figure(2)
|
||||
|
||||
plt.plot(z,1000*rel_diff_Bz_grad_mean, label = r"$(\nabla_z^2 Bz - \nabla_z^2 B_{z, sim}) / mean(\nabla_z^2 Bz)$")
|
||||
#plt.ylim(-57,10)
|
||||
plt.ylabel("relative deviation [‰]")
|
||||
plt.xlabel("z-axis [mm]")
|
||||
plt.legend()
|
||||
plt.savefig("output/HH_benchmark_5A_6x2_rel_deviation_via_mean.pdf")
|
||||
plt.show()
|
||||
|
||||
##################### x-Axis #########################################################
|
||||
|
||||
plt.figure(3)
|
||||
|
||||
plt.rcParams.update({'font.size': 15})
|
||||
plt.suptitle("Helmholtz coil field B_x along x-axis, comparison of simulations", fontsize=30)
|
||||
|
||||
|
||||
#Field plot
|
||||
##########################
|
||||
|
||||
plt.plot(x,B_x,linestyle = "solid", label = r"$B_x$: Result via elliptic integrals")
|
||||
plt.plot(x,B_x_sim,linestyle = "dashdot", label = r"$B_{x, sim}$: Numerical Matlab simulation")
|
||||
plt.plot(x,(B_x-B_x_sim), label = r"$B_x - B_{x, sim}$")
|
||||
#plt.xlim(-0.01,0.01)
|
||||
plt.title("B-field" ,fontsize = 30)
|
||||
|
||||
plt.ylabel(r"$B_x$ [G]")
|
||||
plt.xlabel("x-axis [mm]")
|
||||
plt.legend()
|
||||
|
||||
|
||||
|
||||
|
||||
#################
|
||||
|
||||
|
||||
plt.savefig("output/HH_benchmark_5A_6x2_x-axis.pdf")
|
||||
plt.show()
|
||||
|
||||
print((R2-R1)/R1)
|
||||
@ -26,7 +26,7 @@ f_626 = 478.839e12
|
||||
lambda_626 = 626.082e-9
|
||||
|
||||
#specific constants
|
||||
rho_copper_20 = 1.68e-8
|
||||
rho_copper_20 = 1.72e-8
|
||||
density_copper = 8.940e3
|
||||
|
||||
# water
|
||||
|
||||
Loading…
Reference in New Issue
Block a user