סחיפת סטוקס

[ממוזער|250px|שמאל|עצים שנסחפו מן הים לאורך ה[חוף הצפוני של מדינת וושינגטון. סחיפת סטוקס היא אחת התהליכים החשובים ביותר בהסעה של פסולת ימית.]]
בהקשר של תנועה גלית טהורה במכניקת הזורמים, **מהירות סחיפת סטוקס** (באנגלית: **Stokes drift**) היא המהירות הממוצעת של אלמנט זורם כאשר הוא נע עם שדה הזרימה של הגל. לדוגמה, מצוף שצף על פני השטח של גלי מים, יחווה סחיפה נטו בכיוון התקדמות הגל. התופעה הלא-ליניארית הזאת נקראת על שם ג'ורג' גבריאל סטוקס, שגזר ביטויים מתמטיים שמתארים את הסחיפה במחקרו המעמיק מ-1847 על גלי מים.

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

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

## הסבר אינטואיטיבי
[ממוזער|250px|שמאל|סחיפת סטוקס (לחצו [:קובץ:Deep water wave.gif|כאן להפעלת האנימציה.)]]
סחיפת סטוקס היא תופעה בה המהירות האופקית הממוצעת של אלמנט מים שונה מאפס, כך שמתקבלים מסלולים פתוחים ונוצרת "סחיפה" של אלמנטי הזורם בזמן - מעבר מים נטו בכיוון התקדמות הגל. חלקיקים הנמצאים בגל מים בעל מהירות מופע  C_p  ותלילות סופית יפגינו תנועה שאריתית קלה בגלל שהם שוהים מעט יותר זמן בתנועה קדימה מאשר בתנועה אחורנית, אודות להתקדמות הגל עצמו. תופעה זו אינה ייחודית רק לגלי מים, אף כי במקור סטוקס גילה אותה ב-1847 במחקרו המתמטי המעמיק על גלי ים לא ליניארים. בגלי קול קיימת תופעה דומה של סחיפה כאשר מניחים שאמפליטודת הלחץ של הקול אינה זניחה ביחס ללחץ האטמוספירי - או לחלופין שמהירות אלמנטי האוויר אינה זניחה ביחס למהירות הקול; גם אם המהירות הממוצעת של שדה המהירות בנקודה כלשהי היא אפס, המסלול של חלקיק הנע בכל נקודה במהירות הזרימה המקומית לא יהיה בהכרח חסום.

קל יותר להבין זאת כאשר מסתכלים על שני מקרי קיצון - הראשון כאשר אמפליטודת המהירות בנקודה *(x,t)* אפסית והשני כאשר אמפליטודת המהירות שווה למהירות הקול. במקרה הראשון נקבל שהאלמנט מתנודד הרמונית מסביב לנקודת שיווי המשקל שלו, ואילו במקרה השני נקבל שהחלקיק מתקדם במהירות הגל הנושא של שדה הזרימה - כך שמבחינתו הגל נייח - ולכן מהירותו לא צריכה להשתנות כלל, אלא שהיא עומדת על ערך קבוע ששווה למהירות הסחיפה. סחיפת סטוקס מתייחסת לנקודת ביניים בין שני מקרי הקיצון הללו; היחס בין המהירות של האלמנט למהירות התקדמות הגל עצמו אינו זניח, ומנגד עדיין קטן מספיק כדי שאפשר יהיה להזניח אפקטים מסדרים גבוהים יותר.

[שמאל|250px|ממוזער|סחיפת סטוקס במים רדודים, כאשר אורך הגל גדול משמעותית מעומק המים. לחצו [:קובץ:Shallow water wave.gif|כאן להפעלת האנימציה.]]
במילים אחרות, כאשר החלקיק נע בכיוון התקדמות הגל מהירותו משתנה לאט יותר מאשר כאשר הוא נע בניגוד לתנועת הגל, והיא יורדת בקצב שתלוי במהירות היחסית בינו לבין הגל. כיוון שכך, החלקיק שוהה זמן רב יותר בתנועה עם הגל מאשר בתנועה בניגוד לגל. הדבר מאפשר לו לגמוע מרחק רב יותר בהתקדמות עם הגל לעומת המרחק שהוא חוזר אחורנית כאשר הוא נע בניגוד אליו. התוצאה של האסימטריה הזאת היא התקדמות נטו של החלקיק בכיוון האופקי בכל מחזור גלי.

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

### סחיפת סטוקס חד-ממדית
במקרה החד־ממדי, הגורם המכריע במידת הסחיפה שחווה אלמנט זורם הוא כמה זמן הוא שוהה בתנועה במגמה נתונה. הגודל הפיזיקלי שמכתיב כמה זמן האלמנט יהיה בעל מופע נתון הוא המהירות היחסית בינו לגל:  V - u_x = V - u_0 sin\theta  (כאן *V* היא מהירות הגל), שממנו נגזרת התדירות הזוויתית "האפקטיבית" של תנועת האלמנט: \omega_{eff} = k(V - u_0 sin\theta). סחיפת סטוקס שווה לאינטגרל:  \int_0^{T}u_x dt.

את הדיפרנציאל  *dt* נוח יותר להביע במונחי הפאזה  \theta  כך: dt = \frac  = \frac } . מתקיים מכאן: 

 u_x dt = \frac d\theta . כיוון שאנו חותרים לקבלת הביטוי לסחיפת סטוקס מהסדר הנמוך ביותר, אפשר להניח שמהירות האלמנט  u_x  נמוכה מאוד בהשוואה למהירות הגל *V*, ולכן מותר לקרב את הביטוי על ידי קירוב מסדר ראשון (קירוב ליניארי):  \frac  \approx \frac  + \frac  . על כן נקבל איפה:

 \int_0^{T}u_x dt =  \int_0^{2\pi}(\frac  + \frac )d\theta = \frac .

במעבר האחרון נעזרנו בעובדה ש-:  = \frac . כדי לקבל את מהירות הסחיפה יש לחלק בזמן המחזור: T = \frac  , וכך נקבל: u_{drift,1} = \frac  = \frac  k\omega a^2, כאשר במעבר האחרון נעזרנו בקשרי הגל הבסיסיים:  u_0 = a\omega, V = \frac  .

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

### סחיפת סטוקס בגלי מים (דו-ממדית)
במסגרת פיתוח הסדר הנמוך ביותר של סחיפת סטוקס, יש לחשב מה המרחק שגומא אלמנט נתון עקב גרדיאנט המהירות האופקית של הזרימה יחסית למשתנה העומק:
\frac  . כמו בחישוב החד-ממדי, נעריך את האינטגרל  \int_0^{T}u_x dt הפעם תוך התחשבות בשינוי הגודל של 
 u_x  עם השינוי בעומק. ניתן לפרק את האינטגרל לסכום שני אינטגרלים:  \int_0^{T}u_x dt = \int_0^{T/2}u_x dt + \int_{T/2}^{T}u_x dt , אותו ניתן לפרש גם כ-:  \int_0^{T/2}\frac  d\theta, כאשר  \Delta u_x היא ההפרש בין מהירויות האלמנט במופעים  \theta  ו- -\theta . הפרש זה הוא גם: 
 \Delta u_x = \frac \cdot \Delta z = a\omega k sin\theta \cdot 2a sin \theta = 2\omega k a^2 sin^2\theta , כאשר במעברים האחרונים נסמכנו על העובדה שמסלול האלמנט הוא מעגל בעל רדיוס *a* (במסגרת הקירוב הראשוני) ועל התוצאה לגרדיאנט מהירות הזרימה בפני המים במסגרת התאוריה של איירי. תוצאת האינטגרל היא:  \int_0^{\pi} \frac  d\theta = \pi k a^2.

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

 u_{drift,2} = \frac  = \frac  = \frac  k\omega a^2.

מהירות סחיפת סטוקס הכוללת שווה במקרה הדו-ממדי (על שפת המים) לסכום התרומות של כל אחד מהאפקטים:

u_{drift} = u_{drift,1} + u_{drift,2} = \frac  k\omega a^2 + \frac  k\omega a^2 = k\omega a^2 

## ראו גם
* תאוריית הגלים של איירי

קטגוריה:מכניקת הזורמים