Infrastructure/Stabilization/Description/BasicIdeas.py


# coding: utf-8

# In[25]:

from scipy.constants import pi,mu_0
from IPython.display import display, Math, Markdown
from math import sqrt

Rho_Co=0.0718;display(Markdown(r'Copper resistivity $\rho_{{Co}}$={} Ohm*mm²/m'.format(Rho_Co)))


# # Geometricals

# In[26]:

R_HFS=0.25;display(Markdown(r'Radius of the HFS stabilization coil $R_{{HFS}}$={} m'.format(R_HFS)))
R_LFS=0.65;display(Markdown(r'Radius of the LFS stabilization coil $R_{{LFS}}$={} m'.format(R_LFS)))
S_HFS=pi*R_HFS**2;display(Markdown(r'Surface of the HFS stabilization coil $S_{{HFS}}$={} m$^2$'.format(round(float(S_HFS),3)))) 
S_LFS=pi*R_LFS**2;display(Markdown(r'Surface of the LFS stabilization coil $S_{{LFS}}$={} m$^2$'.format(round(float(S_LFS),3))))
l=0.04;display(Markdown(r'Lenght of the coil could be $l$={} m'.format(round(float(l),3))))


# In[27]:

def CoilInductance(CoilTurns,Surface, length):
    L=mu_0*CoilTurns**2*Surface/length
    return L


# ## Horizontal stabilization
# <img src="HorizontStab.jpg" width="25%">

# In[28]:

L_HFS_UP=CoilInductance(4,S_HFS,l)
display(Markdown(r'Inductance of the individual HFS upper coil is $L_{{HFS}}$={} mH'.format(round(float(L_HFS_UP*1000),3))))
L_HFS_BOT=CoilInductance(4,S_HFS,l)
display(Markdown(r'Inductance of the individual HFS bottom coil is $L_{{HFS}}$={} mH'.format(round(float(L_HFS_BOT*1000),3))))
L_LFS_UP=CoilInductance(4,S_LFS,l)
display(Markdown(r'Inductance of the individual LFS upper coil is $L_{{LFS}}$={} mH'.format(round(float(L_LFS_UP*1000),3))))
L_LFS_BOT=CoilInductance(4,S_LFS,l)
display(Markdown(r'Inductance of the individual LFS bottom coil is $L_{{LFS}}$={} mH'.format(round(float(L_LFS_BOT*1000),3))))


# In[31]:

L_HFS=L_HFS_UP+L_HFS_BOT+2*0.75*sqrt(L_HFS_UP*L_HFS_BOT) # 75% koeficient
display(Markdown(r'Inductance of the HFS coil is $L_{{HFS}}$={} mH'.format(round(float(L_HFS*1000),3))))
L_LFS=L_LFS_UP+L_LFS_BOT+2*0.75*sqrt(L_LFS_UP*L_LFS_BOT) # 75% koeficient
display(Markdown(r'Inductance of the LFS coil is $L_{{LFS}}$={} mH'.format(round(float(L_LFS*1000),3))))


# In[ ]:




# ## Reference
# * [Mutual inductance](https://phys.libretexts.org/Bookshelves/University_Physics/Book%3A_University_Physics_(OpenStax)/Map%3A_University_Physics_II_-_Thermodynamics%2C_Electricity%2C_and_Magnetism_(OpenStax)/14%3A_Inductance/14.02%3A_Mutual_Inductance)

# In[ ]: