פרק 3א - מבוא לקווזיסטטיקה - גלים

From EM Fields - TAU
Jump to navigation Jump to search

מבוא

משוואות מקסוול שקיבלנו בהרצאות הקודמות הינן:

(1) ×E=μ0Ht(2) ×H=ϵ0Et+J(3) (ϵ0E)=ρ(4) (μ0H)=0 נשים לב כי זו מערכת משוואות בה השדות הן פונקציה של 4 משתנים - וקטור המיקום r והזמן t. המשוואות מצומדות זו לזו, כלומר - הפעלת שדה מגנטי משתנה בזמן יוצרת שדה חשמלי משתנה בזמן (משוואה 1), אך משינוי בשדה החשמלי, ישתנה גם השדה המגנטי (משוואה 2). דהיינו לא ניתן למצוא את השדה החשמלי, בלי לדעת מהו השדה המגנטי, ולהיפך.

כפי שהוזכר בהרצאות קודמות, הצימוד דרך המשוואות הוא זה שמאפשר פתרונות גליים. פתרונות אלו, המשתנים בזמן, הם אלו שמאפשרים את האפליקציות הטכנולוגיות ותופעות הטבע כמו אור השמש שמגיע אלינו (פיזור), ולכן אנו לא יכולים להתעלם מהשפעתן של הנגזרות הזמניות.

אחת הדרכים "להסיר" צימוד זה, היא להניח שהשדות סטטיים, ואז הנגזרות הזמניות מתאפסות. הנחה זו מובילה בד"כ לבעיה פשוטה הרבה יותר, אבל כמובן שסטטיקה מושלמת היא מקרה תאורטי, ובנוסף אינה מאפשרת טיפול באופי הגלי של השדה האלקרטומגנטי. עם זאת, היא מהווה נדבך בסיסי בטיפול במשטר ביניים שימושי מאוד - הפתרון הקוואזי-סטטי.

משוואות מקסוול - משטר קווזי סטטי

במשטר קווזיסטטי אנו מניחים שהשדות משתנים בזמן, אך לאט מאוד. מה הכוונה ב"משתנים לאט מאוד"? נראה בהמשך הפרק.

משוואת הגלים - תווך חסר מקורות

בתווך חסר מקורות:

J=0,ρ=0כעת, נפעיל רוטור על שני האגפים של משוואה (1) (חוק פארדיי):

אגף שמאל:

×(×E)=(E)ρ=02E=2Eכאשר במעבר הראשון השתמשנו בזהות וקטורית:

×(×a)=(a)2aאגף ימין:

×(μ0Ht)=μ0×Htמכיוון שהנחנו שכל השדות שאנו עובדים איתם הינם גזירים ורציפים, נחליף את הסדר בין הרוטור המרחבי שפועל על השדה המגנטי, לנגזרת בזמן:

=μ0t(×H)=μ0t(ϵ0Et+JJ=0)Eq. 2 - Ampere's law=μ0ϵ02Et2בסך הכל נקבל:

2E=μ0ϵ02Et2

נגדיר את מהירות הגל להיות c1μ0ϵ0, נעביר אגפים ונקבל את משוואת הגלים:(5) 2E1c22Et2=0זוהי מערכת של 3 משוואות גלים סקלריות לכל רכיב (x^,y^,z^).

משוואת הגלים - תווך חסר מקורות - פתרונות

כעת נרצה לפתור את משוואת הגלים שקיבלנו, ולהבין כיצד נראים השדות החשמליים והמגנטיים שמקיימים אותה.

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

השדה החשמלי

על מנת להבין אינטואיטיבית את פיתרון משוואת הגלים, נניח את הפיתרון הבא:

E=Ex(z,t)x^כלומר, שדה חשמלי שתלוי בזמן ובקורדינטה Z, וכיוונו בקורדינטה X.

נציב את התוצאה בנוסחה (5):

2E1c22Et2=2Ex2=0+2Ey2=0+2Ez21c22Et2=2Ez21c22Et2

כעת קיבלנו משוואה פשוטה יותר, שהפיתרון שלה הוא:

תרשים 1

(6) E=Ex(z,t)x^=(f1[tzc]+f2[t+zc])x^ כאשר f1,f2 פונקציות שרירותיות.


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

השדה המגנטי

כעת, כדי למצוא את השדה המגנטי, נציב את השדה החשמלי שמצאנו לתוך משוואה (1) (חוק פאראדיי):

מצד שמאל:

×E=|x^y^z^xyzExEyEz|=|x^y^z^00zEx00|=Exzy^כאשר איפסנו את הנגזרות בכיוון x,y ואת השדות בכיוונים y,z משום שהשדה שלנו תלוי ב z וכיוונו הוא ב x.

נציב את תוצאת משוואה (6) לתוך תוצאת הרוטור

Exz=f1(tzc)z+f2(t+zc)z= chain rulef1(tzc)(tzc)f'1(tzc)z1c+f2(t+zc)(t+zc)f'2(t+zc)z1c=f'1c+f'2c=1c(f'1f'2)=μ0ϵ0(f'1f2)כעת נציב בחוק פאראדיי:

μ0ϵ0(f'1f2)y^=μ0Htנחלץ את השדה המגנטי:

H=ϵ0μ0(f1(tzc)f2(t+zc))t y^H=ϵ0μ0(f1(tzc)(tzc)f2(t+zc)(t+zc))t y^נשים לב כי:

t=(tzc)=(t+zc)לכן הנגזרות הזמניות של f1,f2 יתבטלו עם הנגזרת הזמנית t. אחרי אינטגרציה בזמן נקבל (נאפס את קבוע האינטגרציה c):

H=ϵ0μ0(f1f2)y^ואכן, כצפוי קיבלנו שהשדה המגנטי ניצב לשדה החשמלי.

תכונות הפתרונות

נרצה לשאול: מהי "מהירות" הגל? כלומר, באיזו מהירות עלינו לנוע על מנת להישאר עם השיא (ראו תרשים 2)?


תרשים 2

נביט בשרטוט בצד שמאל, של f1(tzc):

ברגע t=0 הארגומנט הוא: const=zc

ולפיכך, בזמן כלשהו t , נקבל: tzc=const.

לכן:

z=c(tconst)לכן נבין, כי הגל נע במהירות c, שהגדרנו אותה כמהירות האור:

c1ϵ0μ03108[ms]והשדות הינם:

E=f1(tzc)x^H=ϵ0μ0f1(tzc)y^1ηf1(tzc)y^כאשר η הוא קבוע שנקרא "אימפדנס הואקום", והוא שווה ל:

η377[Ω]עד כה, דנו בפיתרונות כללים למשוואת הגלים (f1,f2 פונקציות ארביטרריות), כעת נרצה למצוא פיתרונות הרמוניים למשוואת הגלים.

משוואת הגלים - פתרונות הרמוניים

הפתרונות ההרמוניים הם מהפתרונות החשובים והשימושיים ביותר למשוואת הגלים, ננחש:

תרשים 3 - שלשה ימנית
תרשים 4

f1(t)=A1cos(ωt)ולכן נקבל את השדות:

E=f1(tzc)x^=A1cos(ω(tzc))x^H=1ηf1(tzc)y^=1ηA1cos(ω(tzc))y^כאשר לעתים קרובות נהוג להגדיר את מספר הגל, כ:

kωc[1m]

כפי שניתן לראות בשרטוטים מצד שמאל, השדה החשמלי (הגל הסגול) מאונך לשדה החשמלי (הגל הצהוב), ושניהם מאונכים לכיוון התקדמות הגל (z), לכן נקבל שלשה ציקלית ימנית: E,H,z^ (תרשימים 3 ו- 4).


נשאל, מדוע הפיתרון ההרמוני שהצגנו נקרא "גל מישורי"?

מכיוון שכל הנקודות שעליהן השדה החשמלי קבוע, נמצאות על מישור שניצב לציר z.

תכונות

זמן מחזור

אם נעמוד בנקודה מסויימת z=z0, ונמדוד את השדה החשמלי כתלות בזמן, קיים זמן כלשהו, שנקרא זמן המחזור ומוגדר:

תרשים 5

T=2πω שעבורו מתקיים:

Ex(z,t0)=Ex(z,t0+T)

כלומר זהו הזמן בקצר ביותר שלאחריו השדה חוזר בדיוק לערכו המקורי (תרשים 5) כאשר השדה מחזורי בזמן.

אורך גל

אם "נקפיא" את הזמן ונמדוד את השדה החשמלי בכיוון x, נקבל את אורך הגל:

תרשים 6

λ=2πk=2πcω=cf

זהו המרחק הקרוב ביותר בני שתי נקודות בהן ערך השדה זהה, כאשר השדה מחזורי במרחב (תרשים 6).

מהירות הפאזה

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

משוואת הגלים - פתרונות הרמוניים - ייצוג פאזורי

לעתים מסובך אלגברית להשתמש בייצוג של סינוסים וקוסינוסים לפונקציות מחזוריות (במיוחד עבור גזירה ואינטגרציה), לכן כדי לפתור בעיה זו נציג את הייצוג הפאזורי של הפתרונות ההרמוניים: E={E~(r,ω)ejωt}=12{E~ejωt+E~*ejωt} H={H~(r,ω)ejωt}=12{H~ejωt+H~*ejωt}

משוואת הגלים - משוואת הלמהולץ

נציב את הייצוג הפאזורי למשוואת הגלים:

2E1c22Et2=122[E~ejωt+E~*ejωt]=2E~ejωt+2E~*ejωt12c22t2[E~ejωt+E~*ejωt]=ω2(E~ejωtE~*ejωt)=0נכפול את הביטוי ב 2, ונכנס אגפים, לפי חזקת האקספוננט:

ejωt(2E~+ω2c2E~)+ejωt(2E~*ω2c2E~*)=0מאחר והאקספוננטים לעולם לא יתאפסו, כל אחד מהאיברים בסוגריים מתאפס זהותית:

2E~+k2E~=0 ; 2E~*k2E~*=0כאשר K הוא מספר הגל (ראינו אותו כבר).

נקבל לבסוף את משוואת הלמהולץ - משוואת הגלים בתחום התדר:2E~+k2E~=0

משוואת הגלים - משוואת הלמהולץ - גל מישורי כללי

נכתוב את השדה החשמלי כ:

E~=E~0ej(kr)נשים לב שאם נבחר את וקטור הגל להיות: k=kz^, נקבל את הפיתרון שראינו מקודם.

נציב את השדה החשמלי החדש למשוואת הלמהולץ, ונקבל:

2E~+ω2c2E~=2(E~0ejkr)+ω2c2E~=(jk)(jk)E~0ejkr+ω2c2E~=(kk)E~0ejkr+ω2c2E~=(kk+ω2c2)E=0

על מנת לקיים את המשוואה, נוכל מחד לאפס את השדה החשמלי, אבל אז נקבל פתרונות לא מעניינים, מנגד ניתן לאפס את הביטוי בסוגריים:

kx2+ky2+kz2=(ωc)2 קיבלנו משוואה שמזכירה משוואה של כדור, ולכן כל הווקטורים האפשריים נמצאים על שפה של כדור, שנקרא Ewald sphere (תרשים 7).

תרשים 7

נציב בחוק גאוס:

(ϵ0E)=(ϵ0E~0ejkr)=jkϵ0E~0ejkr=0ולכן:

kE~0=0כלומר, ווקטור K ניצב לשדה החשמלי.

ניתן לקבל באופו אופן:

kH~0=0ניתן לראות שקיבלנו שוב את השלשה הימנית.

משוואת הגלים - גל מישורי כללי - מישורים שווי פאזה

נביט בתרשים 8, ונחזור חזרה לתחום הזמן:

תרשים 8

E={E~0ejkrejωt}=E0cos(ωtkr)כמו שראינו מקודם, על מנת למצוא את הנקודות שבהם ווקטור השדה החשמלי קבוע, נשווה את הארגומנט של הקוסינוס למספר קבוע. נקבל:

kr=ωtconstהפעם, עבור המקרה הכללי, קיבלנו כי המישור ניצב לוקטור k.

משוואות מקסוול עבור הפאזורים של השדות

אם נרשום את השדות:

E={E~ejωt} ; H={H~ejωt}נציב את השדות החדשים במשוואות מקסוול, ונקבל את משואוות מקסוול הפאזוריות: ×E~=jωμ0H~×H~=jωϵ0E~+J~ϵ0E~=ρ~μH~=0

קוואזי סטטיקה

תרשים 9

על מנת לפתור את הבעיה, שדנו בה בתחילת הפרק, נצטרך להגדיר ולהבדיל בין שדות שמתשנים "מהר" לשדות שמשתנים "לאט". בתרשים 9a ניתן לראות שבמערכת שבה מתקיים L>λ השדות משתנים בצורה משמעותית לאורך המערכת. במקרה זה, נאלץ לספק פיתרון מלא למשוואות מקסוול. בתרשים 9b אנו רואים שתחת התנאי L<<λ השדה משתנה "לאט" ביחס לגודל המערכת.

מה יהיה הזמן, במקרה זה, שלוקח לגל להתפשט במערכת ולחזור לנקודת התחלה? tpropagating=2Lc<<2λc=2c2πcω=22πω=2T

דוגמא

נתון מעגל חשמלי שגודלו 1cm.

התדר האופייני של הסיגנל במערכת הוא f=1KHz.

מה אורך הגל האופייני?

λ=cf=3108103=300Km>>1cmהקירוב שלנו מתקיים, ולכן נוכל לפתור מערכת זו בקירוב קווזי סטטי.