פרק 12 - שדות מגנטיים בחומר: Difference between revisions

From EM Fields - TAU
Jump to navigation Jump to search
No edit summary
RonelM (talk | contribs)
No edit summary
Line 55: Line 55:
\Rightarrow
\Rightarrow
\begin{cases}  
\begin{cases}  
\nabla \cdot (\mu_0\vec H) = \rho_{mf} +\rho_m = \rho_m\\  
\nabla \cdot (\mu_0\vec H) = \underbrace{\rho_{mf}}_{=0 }  +\rho_m = \rho_m\\  
\nabla \times\vec E = -\frac{\partial}{\partial t}(\mu_0\vec H) -\underbrace{\frac{\partial}{\partial t}(\mu_0\vec M)}_{J_m}-J_{mf}\\
\nabla \times\vec E = -\frac{\partial}{\partial t}(\mu_0\vec H) -\underbrace{\frac{\partial}{\partial t}(\mu_0\vec M)}_{J_m}-
\hat n\cdot(\mu_0\vec H_2 - \mu_0\vec H_1) = \eta_{mf} + \eta_m = -\hat n\cdot(\mu_0\vec M_2 - \mu_0\vec M_1)
\underbrace{J_{mf}}_{=0}
\\
\hat n\cdot(\mu_0\vec H_2 - \mu_0\vec H_1) = \underbrace{\eta_{mf}}_{=0} + \eta_m = -\hat n\cdot(\mu_0\vec M_2 - \mu_0\vec M_1)
\end{cases}</math>
\end{cases}</math>


Line 64: Line 66:
\nabla \times \vec E = -\frac{\partial(\mu_0\vec H_a)}{\partial t}\\
\nabla \times \vec E = -\frac{\partial(\mu_0\vec H_a)}{\partial t}\\
\nabla \cdot \vec D = \rho _f\\  
\nabla \cdot \vec D = \rho _f\\  
\nabla \times \vec H_a-\vec M = \frac{\partial \vec D}{\partial t} + \vec J_f\\
\nabla \times (\vec H_a-\vec M) = \frac{\partial \vec D}{\partial t} + \vec J_f\\
\nabla \cdot (\mu_0\vec H_a) = 0\\
\nabla \cdot (\mu_0\vec H_a) = 0\\
\hat n \times (\vec E_2-\vec E_1) = 0\\
\hat n \times (\vec E_2-\vec E_1) = 0\\
Line 74: Line 76:
\nabla \cdot \vec D = \rho _f\\  
\nabla \cdot \vec D = \rho _f\\  
\nabla \times \vec H = \frac{\partial \vec D}{\partial t} + \vec J_f\\
\nabla \times \vec H = \frac{\partial \vec D}{\partial t} + \vec J_f\\
\nabla \cdot (\mu_0\vec [\vec H+\vec M]) = 0\\
\nabla \cdot (\mu_0 [\vec H+\vec M]) = 0\\
\hat n \times (\vec E_2-\vec E_1) = 0\\
\hat n \times (\vec E_2-\vec E_1) = 0\\
\hat n \cdot (\vec D_2-\vec D_1) = \eta_f\\
\hat n \cdot (\vec D_2-\vec D_1) = \eta_f\\
Line 105: Line 107:
כדור בעל מגנטיזציה אחידה. מהו <math>\vec B</math> בכל המרחב?
כדור בעל מגנטיזציה אחידה. מהו <math>\vec B</math> בכל המרחב?


נשתמש במודל המטען:<math display="block">\eta_m = -\hat n\cdot(\vec M_{out}-\vec M_{in})\mu_0 = -\hat r\cdot(0-M\hat z\mu_0) = M\hat r\cdot\hat z = M\cos\theta\mu_0</math>צפיפות המטען:<math display="block">\rho_m = -\nabla\cdot(\mu_0\vec M) = 0</math>נפתור באמצעות פוטנציאל סקלרי:<math display="block">\nabla\times\vec H = \vec J_f + \frac{\partial \vec D}{\partial t} + \vec J_a = 0 \Rightarrow \vec H = -\nabla\phi_m</math>נציב ונקבל ממקסוול:<math display="block">\nabla\cdot(\mu_0\vec H) = \rho_m = 0 \Rightarrow \nabla\cdot(\mu_0\cdot(-\nabla\phi_m)) = 0</math>קיבלנו את משוואת לפלס:<math display="block">\nabla^2\phi_m = 0</math>נפתור את משוואת לפלס עם מקורות משטחיים בלבד:<math display="block">\begin{cases}  
נשתמש במודל המטען:<math display="block">\eta_m = -\hat n\cdot(\vec M_{out}-\vec M_{in})\mu_0 = -\hat r\cdot(0-M\hat z\mu_0) = M\hat r\cdot\hat z = M\cos\theta\mu_0</math>צפיפות המטען:<math display="block">\rho_m = -\nabla\cdot(\mu_0\vec M) = 0</math>נפתור באמצעות פוטנציאל סקלרי:<math display="block">\nabla\times\vec H = \underbrace{\vec J_f}_{=0} +  
\underbrace{\frac{\partial \vec D}{\partial t}}_{=0 \text{ static}} +  
\underbrace{\vec J_a}_{=0 \text{ Not using this field} }
= 0 \Rightarrow \vec H  
= -\nabla\phi_m</math>נציב ונקבל ממקסוול:<math display="block">\nabla\cdot(\mu_0\vec H) = \rho_m = 0 \Rightarrow \nabla\cdot(\mu_0\cdot(-\nabla\phi_m)) = 0</math>קיבלנו את משוואת לפלס:<math display="block">\nabla^2\phi_m = 0</math>נפתור את משוואת לפלס עם מקורות משטחיים בלבד:<math display="block">\begin{cases}  
\phi_m(r>>a)\rightarrow0\\
\phi_m(r>>a)\rightarrow0\\
\phi_m(r\rightarrow0)<\infty\\  
\phi_m(r\rightarrow0)<\infty\\  
\hat n \times (\vec H_2-\vec H_1) = \vec K_f = 0\\
\hat n \times (\vec H_2-\vec H_1) = \vec K_f = 0\\
\hat n \cdot (\mu_0\vec H_2 - \mu_0\vec H_1) = \eta_m = \mu_0M\cos\theta
\hat n \cdot (\mu_0\vec H_2 - \mu_0\vec H_1) = \eta_m = \mu_0M\cos\theta
\end{cases}</math>נבחר פתרון כללי <math>(l=0, n=1)</math>:<math display="block">\phi = (c_1r+\frac{c_2}{r^2})\cos\theta</math><math display="block">\phi_{m_1} =Ar\cos\theta \quad , \quad \phi_{m_2} =\frac{C}{r^2}\cos\theta</math>נציב בתנאי השפה:<math display="block">Aa\cos\theta = \frac{C}{a^2}\cos\theta \Rightarrow a^3A = C</math>מתנאי השפה האחרון:<math display="block">\hat r \cdot [\mu_0\cdot(-\nabla\phi_{m_2}) - \mu_0(-\nabla\phi_{m_1})] = \mu_0M\cos\theta </math><math display="block">-\frac{\partial \phi_{m_2}}{\partial r} + \frac{\partial \phi_{m_1}}{\partial r} = M\cos\theta \Rightarrow -[\frac{-2C}{a^3}\cos\theta]+A\cos\theta=M\cos\theta \Rightarrow \frac{2C}{a^3}+A=M  </math>נקבל את המקדמים:<math display="block">A=\frac{M}{3} \quad, \quad C = a^3\frac{M}{3}</math>נציב את המקדם חזרה בפוטנציאל הראשון:<math display="block">\phi_{m_1} =\frac{M}{3}r\cos\theta \quad \Rightarrow \vec H_1 = -\nabla\phi_{m_1} = -\frac{M}{3}\hat z</math>נמצא את השדה המגנטי:<math display="block">\Rightarrow \vec B_1 = \mu_0\cdot(\vec H_1 +\vec M) = \mu_0\cdot(-\frac{M}{3}\hat z+M\hat z)=\frac{2}{3}\mu_0M\hat z</math>כעת נציב את המקדם בפוטנציאל השני:<math display="block">\phi_{m_1} =\frac{M}{3}\frac{a^3}{r^2}\cos\theta \quad \Rightarrow \vec H_2 = -\nabla\phi_{m_2} = \frac{Ma^3}{3r^3}[2\cos\theta\hat r + \sin\theta\hat\theta] \quad , \vec B_2 = \mu_0\vec H_2 </math>תזכורת - שדה מגנטי של דיפול:<math display="block">\vec H_{dip} = \frac{m}{4\pi r^3}[2\cos\theta\hat r + \sin\theta\hat\theta] </math>נשווה מקדמים ונקבל:<math display="block">\frac{m}{4\pi} = \frac{Ma^3}{3} \Rightarrow m = M\cdot\underbrace{\frac{4}{3}\pi a^3}_{V_{ball}} </math>
\end{cases}</math>נבחר פתרון כללי <math>(l=0, n=1)</math>:<math display="block">\phi = (c_1r+\frac{c_2}{r^2})\cos\theta</math><math display="block">\phi_{m_1} =Ar\cos\theta \quad , \quad \phi_{m_2} =\frac{C}{r^2}\cos\theta</math>נציב בתנאי השפה:<math display="block">Aa\cos\theta = \frac{C}{a^2}\cos\theta \Rightarrow a^3A = C</math>מתנאי השפה האחרון:<math display="block">\hat r \cdot [\mu_0\cdot(-\nabla\phi_{m_2}) - \mu_0(-\nabla\phi_{m_1})] = \mu_0M\cos\theta </math><math display="block">-\frac{\partial \phi_{m_2}}{\partial r} + \frac{\partial \phi_{m_1}}{\partial r} = M\cos\theta \Rightarrow -[\frac{-2C}{a^3}\cos\theta]+A\cos\theta=M\cos\theta \Rightarrow \frac{2C}{a^3}+A=M  </math>נקבל את המקדמים:<math display="block">A=\frac{M}{3} \quad, \quad C = a^3\frac{M}{3}</math>נציב את המקדם חזרה בפוטנציאל הראשון:<math display="block">\phi_{m_1} =\frac{M}{3}r\cos\theta \quad \Rightarrow \vec H_1 = -\nabla\phi_{m_1} = -\frac{M}{3}\hat z</math>נמצא את השדה המגנטי:<math display="block">\Rightarrow \vec B_1 = \mu_0\cdot(\vec H_1 +\vec M) = \mu_0\cdot(-\frac{M}{3}\hat z+M\hat z)=\frac{2}{3}\mu_0M\hat z</math>כעת נציב את המקדם בפוטנציאל השני:<math display="block">\phi_{m_2} =\frac{M}{3}\frac{a^3}{r^2}\cos\theta \quad \Rightarrow \vec H_2 = -\nabla\phi_{m_2} = \frac{Ma^3}{3r^3}[2\cos\theta\hat r + \sin\theta\hat\theta] \quad , \vec B_2 = \mu_0\vec H_2 </math>תזכורת - שדה מגנטי של דיפול:<math display="block">\vec H_{dip} = \frac{m}{4\pi r^3}[2\cos\theta\hat r + \sin\theta\hat\theta] </math>נשווה מקדמים ונקבל:<math display="block">\frac{m}{4\pi} = \frac{Ma^3}{3} \Rightarrow m = M\cdot\underbrace{\frac{4}{3}\pi a^3}_{V_{ball}} </math>


=== יחסי חוקה - סוספטביליות מגנטית, פרמאביליות ===
=== יחסי חוקה - סוספטביליות מגנטית, פרמאביליות ===
אם הקשר לינארי:<math display="block">\vec M = \chi_m\vec H \Rightarrow \vec B = \mu_0(\vec H + \vec M) = \mu_0(1+\chi_m)\vec H </math>כאשר <math>\chi_m </math> הסוספטביליות המגנטית.
אם הקשר לינארי:<math display="block">\vec M = \chi_m\vec H \Rightarrow \vec B = \mu_0(\vec H + \vec M) =  
\overbrace{\mu_0
\underbrace{(1+\chi_m)}_{\equiv \mu_r}}^{\equiv \mu}
\vec H </math>כאשר <math>\chi_m </math> הסוספטביליות המגנטית.


==== משוואות מקסוול בחומר לינארי ====
==== משוואות מקסוול בחומר לינארי ====
Line 147: Line 156:
|אופי התגובה
|אופי התגובה
|}
|}
</div>

Revision as of 13:33, 20 April 2022

שדות מגנטיים בחומר

מנגנוני מגנטיזציה

  • Spin Magnetization
  • Orbital Magnetization

דיאמגנטים - Orbital Magnetization

מנגנון התגובה ל-H הוא דרך שינוי השטף. מתוך עקרון לנץ, דיפול התגובה יהיה הפוך לשדה החיצוני.

פאראמגנטים, פרומגנטים - Spin Magnetization

וקטור מגנטיזציה - M

נסתכל על אלמנט מגנטיזציה קטן:dm=m=Mdvdmdv=Mישנם שני מודלים לתיאור המקורות השקולים המייצגים את המגנטיזציה:

  1. מודל הזרם האמפרי
  2. מודל המטען המגנטי

1.מודל הזרם האמפרי

כאשר באזור מסוים משתנה המגנטיזציה, תהיה צפיפות זרם שקולה המייצגת שינוי זה.

נרצה לשכנע שמתקיים: Ja=×M. נתחיל מלהסתכל שוב על אלמנט מגנטיזציה קטן:dm=M(dlda)=(Mdl)daמתקיים I=Mdl ולכן:dm=Idaקיבלנו את התוצאה שקיבלנו דרך מגנטוסטטיקה עבור מומנט הדיפול של לולאת זרם בשטח da.

מה סך הזרם שעובר דרך הלולאה שהגדרנו?I=Mdlמצד אחד, ישנו הקשר בין הזרם לצפיפות הזרם:I=Jadaמצד שני, לפי משפט סטוקס נוכל לומר:Mdl=×Mdaמאחר שאין תלות בלולאה בה נבחר, נקבל את השוויון:Ja=×Mוהוכחנו.

זרמי מגנטיזציה משטחיים

נמצא תנאי שפה במעבר בין תווכים בהם H שונה:×H=J+Dtn^×(H2H1)=kובין תווכים בהם M שונה:×M=Jan^×(M2M1)=ka

משוואות מקסוול בחומר

נוכל לרשום את משוואות מקסוול בנוכחות מגנטיזציה:{×E=(μ0Ha)tD=ρf×Ha=Dt+Jf+×MJa(μ0Ha)=0ותנאי השפה:{n^×(E2E1)=0n^(D2D1)=ηfn^×(Ha,2Ha,1)=Kf+n^×(M2M1)Kan^(μ0Ha,2μ0Ha,1)=0

2. מודל המטענים המגנטיים

נבנה את המודל באמצעות השוואה בין הפולרזיציה לבין המגנטיזציה:Pμ0Mצפיפות המטען הנפחית:ρp= Pρm=(μ0M)צפיפות הזרם:Jp=ptJm=t(μ0M)צפיפות המטען המשטחית:ηp=n^(P2P1)ηm=n^(μ0M2μ0M1)

חוק שימור המטען המגנטי

קיבלנו את הביטוי לצפיפות המטען המשטחית:ηm=n^(μ0M2μ0M1)נגזור אותו בזמן:ηmt=n^(μ0M2tμ0M1t)=n^(Jm2Jm1)וקיבלנו את חוק שימור המטען המגנטי:ηmt=n^(Jm2Jm1)

משוואות מקסוול במודל המטען (אנלוגיה עם מודל הפולריזציה החשמלית)

נרשום את משוואות מקסוול:{(ϵ0E)=ρf+(P)×H=(ϵ0E)t+Jf+Ptn^(ϵ0E2ϵ0E1)=ηf+(n^[P2P1]){(μ0H)=ρmf=0+ρm=ρm×E=t(μ0H)t(μ0M)JmJmf=0n^(μ0H2μ0H1)=ηmf=0+ηm=n^(μ0M2μ0M1)

סיכום המודלים - משוואות מקסוול בחומר

מודל הזרם האמפרי:{×E=(μ0Ha)tD=ρf×(HaM)=Dt+Jf(μ0Ha)=0n^×(E2E1)=0n^(D2D1)=ηfn^×([Ha,2M2][Ha,1M1])=Kfn^(μ0Ha,2μ0Ha,1)=0מודל המטען המגנטי:{×E=(μ0Ha)tD=ρf×H=Dt+Jf(μ0[H+M])=0n^×(E2E1)=0n^(D2D1)=ηfn^×(H2H1)=Kfn^(μ0[H2+M2]μ0[H1+M1])=0נשים לב לכך שאם נגדיר H+M=Ha נקבל בדיוק את אותן משוואות!

משוואות מקסוול בחומר - צפיפות השטף המגנטי

{×E=BtD=ρf×H=Dt+JfB=0n^×(E2E1)=0n^(D2D1)=ηfn^×(H2H1)=Kfn^(B2B1)=0נגדיר B=μ0Ha צפיפות השטף המגנטי.

תזכורת: D=ϵ0E+P

דוגמה 1

גליל קטן בעל מגנטיזציה אחידה M=Mz^.

מודל המטען:ρm=(μ0M)=0צפיפות המטען המשטחית על חלקו העליון של הגליל:ηm,top=z^(0μ0Mz^)=μ0Mצפיפות המטען המשטחית בתחתית הגליל:ηm,bottom=z^(μ0Mz^0)=μ0Mרחוק מאוד מהגליל נראה דיפול בעל מגנטיזציה:m=MV=Mπa2hz^אם נסתכל על הגליל כדיפול נקבל:μ0m=μ0Mπa2hz^קיבלנו את אותו הביטוי! כעת אפשר להציב בביטוי לשדה דיפולי.

מודל הזרם האמפרי:Ja=×M=0Ka=r^×(0Mz^)=Mφ^Ha=H+MB=μ0Ha=μ0(H+M)connection between models

דוגמה 2

כדור בעל מגנטיזציה אחידה. מהו B בכל המרחב?

נשתמש במודל המטען:ηm=n^(MoutMin)μ0=r^(0Mz^μ0)=Mr^z^=Mcosθμ0צפיפות המטען:ρm=(μ0M)=0נפתור באמצעות פוטנציאל סקלרי:×H=Jf=0+Dt=0 static+Ja=0 Not using this field=0H=ϕmנציב ונקבל ממקסוול:(μ0H)=ρm=0(μ0(ϕm))=0קיבלנו את משוואת לפלס:2ϕm=0נפתור את משוואת לפלס עם מקורות משטחיים בלבד:{ϕm(r>>a)0ϕm(r0)<n^×(H2H1)=Kf=0n^(μ0H2μ0H1)=ηm=μ0Mcosθנבחר פתרון כללי (l=0,n=1):ϕ=(c1r+c2r2)cosθϕm1=Arcosθ,ϕm2=Cr2cosθנציב בתנאי השפה:Aacosθ=Ca2cosθa3A=Cמתנאי השפה האחרון:r^[μ0(ϕm2)μ0(ϕm1)]=μ0Mcosθϕm2r+ϕm1r=Mcosθ[2Ca3cosθ]+Acosθ=Mcosθ2Ca3+A=Mנקבל את המקדמים:A=M3,C=a3M3נציב את המקדם חזרה בפוטנציאל הראשון:ϕm1=M3rcosθH1=ϕm1=M3z^נמצא את השדה המגנטי:B1=μ0(H1+M)=μ0(M3z^+Mz^)=23μ0Mz^כעת נציב את המקדם בפוטנציאל השני:ϕm2=M3a3r2cosθH2=ϕm2=Ma33r3[2cosθr^+sinθθ^],B2=μ0H2תזכורת - שדה מגנטי של דיפול:Hdip=m4πr3[2cosθr^+sinθθ^]נשווה מקדמים ונקבל:m4π=Ma33m=M43πa3Vball

יחסי חוקה - סוספטביליות מגנטית, פרמאביליות

אם הקשר לינארי:M=χmHB=μ0(H+M)=μ0(1+χm)μrμHכאשר χm הסוספטביליות המגנטית.

משוואות מקסוול בחומר לינארי

נוכל לעדכן את משוואות מקסוול עבור חומרים לינאריים:{×E=(μH)t(ϵE)=ρf×H=(ϵE)t+Jf(μH)=0n^×(E2E1)=0n^(ϵ2E2ϵ1E1)=ηfn^×(H2H1)=Kfn^(μ2H2+μ1H1)=0

חומרים לא מגנטיים (תגובה חלשה)
פאראמגנטים דיאמגנטים סוג החומר
0<χm<<1 |χm|<<1,χm<0 χm
חומרים מגנטיים (תגובה חזקה)
פרומגנטים פרימגנטים סוג החומר
תגובה חזקה מאוד,

בד"כ לא לינארית

תגובה חזקה מאוד אופי התגובה