לפלסיאן

במתמטיקה ופיזיקה, אופרטור לפלס או לפלסיאן, המסומל באמצעות
<math>\Delta\,</math>
או
<math>\nabla^2</math>
ונקרא על שם פייר-סימון לפלס, הוא אופרטור דיפרנציאלי, ובפרט אופרטור אליפטי, בעל שימושים רבים, הפועל על פונקציות סקלריות. בפיזיקה, הוא משמש למשל במודלים מתמטיים של התפשטות גלים ושל הולכת חום, וכן במשוואת הלמהולץ. הלפלסיאן הוא בעל חשיבות מרכזית באלקטרוסטטיקה ובמכניקת הזורמים, ומשמש במשוואת לפלס ומשוואת פואסון. במכניקת הקוונטים, הוא מייצג את רכיב האנרגיה הקינטית במשוואת שרדינגר. במתמטיקה, פונקציה אשר הלפלסיאן שלה מתאפס נקראת פונקציה הרמונית. הלפלסיאן הוא מרכיב ליבה בתורת הודג' ובתוצאותיה של קוהומולוגיית דה רהם.

==הגדרה==
אופרטור לפלס הוא אופרטור דיפרנציאלי מסדר שני במרחב אוקלידי n-ממדי, המוגדר כדיברגנץ (<math>\nabla\cdot</math>) של הגרדיאנט (<math>\vec{\nabla} f</math>). אם <math>\ f </math> היא פונקציה ממשית הגזירה פעמיים, אז הלפלסיאן של <math>\ f </math> מוגדר על ידי

:<math>\Delta f=\nabla^2 f=\nabla \cdot \vec {\nabla} f </math> (1)

באופן שקול, הלפלסיאן של f הוא סכום כל הנגזרות החלקיות השניות הבלתי מעורבות בקואורדינטות הקרטזיות <math>\ x_i</math>:

:<math>\ \Delta f=\sum_{i=1}^n \frac {\partial^2 f}{\partial x^2_i}</math> (2)

כאופרטור דיפרנציאלי מסדר שני, אופרטור לפלס ממפה פונקציות-C<sup>k</sup> לפונקציות- 
<math>\ C^\left( k-2 \right) </math>
עבור
<math>\ k\ge 2 </math>. 
הביטויים לעיל מגדירים אופרטור
<math>\ \triangle : C^k \left( \mathbb{R}^n \right) \rightarrow C^{k-2} \left( \mathbb{R}^n \right) </math>,
או באופן כללי יותר אופרטור
<math>\ \triangle : C^k \left( \Omega \right) \rightarrow C^{k-2} \left( \Omega \right) </math>
לכל קבוצה פתוחה <math>\ \Omega </math>.

הלפלסיאן של פונקציה הוא גם העקבה של מטריצת הסיאן של הפונקציה, הגדרה נוחה לשימוש באלגברה ליניארית ובסטטיסטיקה:
:<math>\ \Delta f=\mathrm{tr}(H(f))</math>

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

בפרט, אם u היא הצפיפות במצב שיווי משקל של גודל כלשהו (כגון ריכוז כימי), אזי השטף של u דרך שפתו של התחום V הוא אפס:

:<math>\iint_{\partial V} \nabla u \cdot \mathbf{n}\, dS=0</math>

כאשר n הוא וקטור יחידה הניצב לשפתו של V. ממשפט הדיברגנץ נקבל:

:<math>\iiint_V \mathrm{div} \nabla u\, dV=\iint_{\partial V} \nabla u\cdot\mathbf{n}\, dS=0</math>

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

:<math>\mathrm{div} \nabla u= \nabla \cdot \nabla u = \Delta u=0</math>

כאן, צד שמאל של המשוואה הוא אופרטור לפלס.

==ביטויים במערכות צירים שונות==

===שני ממדים===
אופרטור לפלס בשני ממדים בקואורדינטות קרטזיות הוא:

:<math>\Delta f=\frac{\partial^2f}{\partial x^2} + \frac{\partial^2 f}{\partial y^2}</math>

ובקואורדינטות קוטביות:

:<math> \Delta f
={1 \over r} {\partial \over \partial r}
 \left( r {\partial f \over \partial r} \right)
+ {1 \over r^2} {\partial^2 f \over \partial \theta^2}
</math>

===שלושה ממדים===
בשלושה ממדים, בקואורדינטות קרטזיות:

:<math>
\Delta f=\frac{\partial^2 f}{\partial x^2} + \frac{\partial^2 f}{\partial y^2} + \frac{\partial^2 f}{\partial z^2}.
</math>

בקואורדינטות גליליות:

:<math> \Delta f
={1 \over \rho} {\partial \over \partial \rho}
 \left( \rho {\partial f \over \partial \rho} \right)
+ {1 \over \rho^2} {\partial^2 f \over \partial \theta^2}
+ {\partial^2 f \over \partial z^2 }. 
</math>

ובקואורדינטות כדוריות:

:<math> \Delta f
={1 \over r^2} {\partial \over \partial r}
 \left( r^2 {\partial f \over \partial r} \right)
+ {1 \over r^2 \sin \theta} {\partial \over \partial \theta}
 \left( \sin \theta {\partial f \over \partial \theta} \right)
+ {1 \over r^2 \sin^2 \theta} {\partial^2 f \over \partial \phi^2}.
</math>

כאשר 
<math> \theta \ </math>
היא הזווית הקוטבית (כלומר, הזווית מטה מציר z ו-
<math>\phi</math>
היא הזווית האזימוטית (מציר x).

את הביטוי 
<math>{1 \over r^2} {\partial \over \partial r} \left( r^2 {\partial f \over \partial r} \right)</math>
ניתן להחליף בביטוי השקול 
<math>{1 \over r} {\partial^2 \over \partial r^2} \left( r f \right)</math>.

===N ממדים===
בקואורדינטות כדוריות ב-<math>N</math> ממדים, עם הפרמטריזציה 
<math>x=r\theta \in {\mathbb R}^N</math>
כאשר
<math>r \in [0,+\infty)</math>
וכן
<math> \theta \in S^{N-1}</math>,

:<math> \Delta f
=\frac{\partial^2 f}{\partial r^2}
+ \frac{N-1}{r} \frac{\partial f}{\partial r}
+ \frac{1}{r^2} \Delta_{S^{N-1}} f
</math>

כאשר
<math>\Delta_{S^{N-1}}</math>
הוא אופרטור לפלס-בלטרמי בכדור <math>N-1</math>-ממדי, או לפלסיאן כדורי.

את הביטוי 
<math>{\partial^2 f \over \partial r^2}
+ \frac{N-1}{r} \frac{\partial f}{\partial r}</math>
ניתן להחליף בביטוי השקול 
<math>\frac{1}{r^{N-1}} \frac{\partial}{\partial r} \Bigl(r^{N-1} \frac{\partial f}{\partial r} \Bigr)</math>.
כתוצאה מכך, ניתן לחשב את הלפלסיאן הכדורי של פונקציה המוגדרת על 
<math>S^{N-1}\subset{\mathbb R}^N </math>
בתור הלפלסיאן הרגיל של הפונקציה, אשר הורחב ל-
<math>{\mathbb R}^N \setminus\{0\}</math>
כך שהוא קבוע לאורך הקרניים.

==זהויות==

אם f ו-g הן פונקציות, אז הלפלסיאן של מכפלתן יהיה

:<math>\Delta(fg)=(\Delta f)g+2((\nabla f)\cdot(\nabla g))+f(\Delta g)</math>

מקרה מיוחד הוא זה שבו f היא פונקציה רדיאלית <math>\ f(r)</math> ו-g היא הרמוניה ספרית, <math>\ Y_{lm}(\theta,\phi)</math>. מקרה זה מופיע במודלים פיזיקליים רבים. הגרדיאנט של <math>\ f(r)</math> הוא וקטור רדיאלי והגרדיאנט של פונקציה זוויתית משיק לווקטור זה, לכן

:<math>2(\nabla f(r))\cdot(\nabla Y_{lm}(\theta,\phi))=0</math>

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

:<math>\Delta Y_{\ell m}(\theta,\phi)=-\frac{\ell(\ell+1)}{r^2} Y_{\ell m}(\theta,\phi)</math>

מכאן,

:<math>\Delta( f(r)Y_{\ell m}(\theta,\phi) )=\left(\frac{d^2f(r)}{dr^2} + \frac{2}{r} \frac{df(r)}{dr} - \frac{\ell(\ell+1)}{r^2} f(r)\right)Y_{\ell m}(\theta,\phi)</math>

==הכללות==
ניתן להכליל את הלפלסיאן למרחבים לא אוקלידים, בהן הוא יכול להיות אופרטור אליפטי, היפרבולי או אולטרה-היפרבולי.

במרחב מינקובסקי הופך הלפלסיאן לאופרטור ד'אלמבר (הנקרא גם ד'אלמברטיאן):

<math>\square
=
\frac {1}{c^2}{\partial^2 \over \partial t^2 }
-
{\partial^2 \over \partial x^2 }
-
{\partial^2 \over \partial y^2 }
-
{\partial^2 \over \partial z^2 }
</math>

הד'אלמברטיאן משמש למשל על מנת לרשום את משוואת הגלים בארבעה ממדים ואת משוואת קליין-גורדון. נציין כי הסימנים לפני הנגזרות המרחביות הם שליליים, בעוד במרחב אוקלידי הם יהיו חיוביים. המקדם המכיל את c נדרש אם המרחב והזמן נמדדים ביחידות שונות; מקדם נוסף ידרש אם, לדוגמה, כיוון x היה נמדד במטרים בעוד כיוון y היה נמדד בסנטימטרים. הפיזיקאים עובדים בדרך כלל במערכת של יחידות טבעיות כגון יחידות פלאנק, בהן c=1, על מנת לפשט את המשוואה.

===אופרטור לפלס-בלטרמי===
הכללה נוספת של הלפלסיאן היא אופרטור אליפטי הנקרא אופרטור לפלס-בלטרמי, המוגדר על יריעה רימנית. בדומה, ניתן להכליל את אופרטור ד'אלמבר לאופרטור היפרבולי על יריעה פסאודו-רימנית. את אופרטור לפלס-בלטרמי ניתן להמשיך ולהכליל לשדה טנזורי.

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

==הלפלסיאן הווקטורי==
בעוד הלפלסיאן הסקלרי שהוגדר לעיל פועל על שדה סקלרי ומחזיר סקלר, הלפלסיאן הווקטורי פועל של שדה וקטורי ומחזיר וקטור.

===הגדרה===
הלפלסיאן הווקטורי של שדה וקטורי <math> \mathbf{A} </math> מוגדר על ידי

:<math> \nabla^2 \mathbf{A}=\nabla(\nabla \cdot \mathbf{A}) - \nabla \times (\nabla \times \mathbf{A}) </math>

בקואורדינטות קרטזיות הוא יוצג כך:

:<math> \nabla^2 \mathbf{A}=(\nabla^2 A_x, \nabla^2 A_y, \nabla^2 A_z) </math>

כאשר <math>\,A_x</math>, <math>\,A_y</math> ו-<math>\,A_z</math> הם רכיבי הווקטור <math>\mathbf{A}</math>. (לביטויים במערכות צירים אחרות ראו דל במערכות צירים שונות)

===הכללה לטנזורים===
הלפלסיאן של שדה טנזורי <math>T</math> (כאשר סקלר ווקטור הם מקרים פרטיים של טנזור) מוגדר כדיברגנץ של הגרדיאנט של הטנזור:

:<math>\nabla^2 T=\nabla \cdot (\nabla T)</math>

במקרה המיוחד בו <math>T</math> הוא סקלר (טנזור מדרגה 0), נקבל את הלפלסיאן הסקלרי הרגיל. 
אם <math>T</math> הוא וקטור (טנזור מדרגה 1), הגרדיאנט שלו הוא נגזרת קו-וריאנטית אשר תוצאה הפעלתה הוא טנזור מדרגה 2, והדיברגנץ של תוצאה זו הוא שוב וקטור. המשוואה ללפלסיאן הווקטורי לעיל עשויה לשמש על מנת להימנע מחשבון טנזורי וניתן להראות כי הוא שקולה לדיברגנץ של הגרדיאנט של הווקטור.

===שימושים בפיזיקה===
דוגמה לשימוש בלפלסיאן הווקטורי ניתן למצוא במשוואות נאוויה-סטוקס לנוזל ניוטוני בלתי דחיס:

:<math>\rho \left(\frac{\partial \mathbf{v}}{\partial t}+ ( \mathbf{v} \cdot \nabla ) \mathbf{v}\right)=\rho \mathbf{f}-\nabla p +\mu\left(\nabla ^2 \mathbf{v}\right)</math>

כאשר הביטוי הכולל את הלפלסיאן הווקטורי של שדה המהירות <math>\mu\left(\nabla ^2 \mathbf{v}\right)</math> מייצג את מאמץ הצמיגות בנוזל.

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

:<math>\nabla^2 \mathbf{E} - \mu_0 \epsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2}=\nabla^2 \mathbf{E} - \frac{1}{c^2} \frac{\partial^2 \mathbf{E}}{\partial t^2}=0</math>

את משוואה זו ניתן לכתוב גם בצורה

:<math>\Box^2\, \mathbf{E}=0</math>
כאשר
:<math>\Box^2=\frac{1}{c^2} \frac{\partial^2}{\partial t^2}-\nabla^2</math>
הוא הד'אלמברטיאן.

==קישורים חיצוניים==
* 

==הערות שוליים==

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