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

החידה: משהו חוצה, ושום דבר אינו עושה את הדרך¶

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

קשרו חבל באורך 20 מ’ לקיר, מתחו אותו, והצליפו בקצה. משהו רץ לאורך החבל ומכה בקיר שבריר שנייה אחר כך — אם תניחו את היד השנייה על הקיר תרגישו אותו מגיע. אף חלק של החבל לא עשה את הדרך הזו.

לנבא לפני שמחשבים¶

התחייבו לתשובה לכל אחת לפני שתריצו משהו.

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

לחקור את המודל¶

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

פולס רץ 2.6 מ’ לאורך מיתר במהירות 20 מ’/ש’. נקודה אחת של המיתר, 2 מ’ מן הקצה, מסומנת בכתום, עם שובל דהוי של כל מקום שבו הייתה. השובל הוא קו אנכי. הפאנל התחתון רושם את גובה הנקודה כנגד הזמן: היא עולה 10 מ"מ כשחזית הפולס מגיעה, יורדת כשגבו עוזב, ומסיימת במרחק של פחות מ- $6\times10^{-5}$ מ"מ מן המקום שבו התחילה. מה שחצה את המיתר, דבר מן המיתר לא הלך איתו.

מיתר שנמשך למשולש בגובה 10 מ"מ ושוחרר ממנוחה. הוא מתפצל לשני משולשים בחצי הגובה, הרצים זה מזה. המילויים השקופים שמתחת הם שני החצאים של ד’אלמבר, $y_0(x - vt)/2$ בכחול ו- $y_0(x + vt)/2$ בכתום: בפריים הראשון הם מונחים זה על זה והמיתר הוא סכומם; כשהם מחליקים זה מזה המיתר עוקב. הפותר תואם את הנוסחה עד $2\times 10^{-14}$ מן הגובה לאורך כל הריצה — הוא נמצא בצעד הזמן המיוחד שנגזר למטה, שבו פריטה מועברת במדויק.

אותה פריטה מחושבת פעמיים, עם צעד זמן הקטן ב-1% מגבול היציבות (כחול) וגדול ממנו ב-1% (כתום). במשך מאה צעדים שני המיתרים נראים זהים. הפאנל התחתון מראה מדוע זה מטעה: על ציר לוגריתמי הוא עוקב אחר הגל הקצר ביותר שהרשת יכולה להחזיק, ובריצה הכתומה הוא מטפס בקו ישר משגיאת העיגול, ומוכפל פי 1.33 בכל צעד. בצעד 137 הוא מגיע לגודל הפולס עצמו, ותוך תריסר צעדים פאנל המיתר מוצף. שום דבר פיזיקלי לא קרה.

מפרט המודל

המערכת: שדה ההזזה הרוחבית $y(x,t)$ של מיתר אחיד בעל מתיחות $T$ וצפיפות מסה אורכית $\mu$ , באורך $L$ . נצפים: ההזזה והמהירות הרוחבית של כל נקודה, וזמני ההגעה של פולסים לגלאים קבועים.
הדינמיקה: $\mu\,y_{tt} = T\,y_{xx}$ , מוערכת במדויק בפתרון ד’אלמבר ומקודמת נומרית בסכמת leapfrog, שהסטנסיל שלה הוא חוק ניוטון לשרשרת מסות $\mu\,\Delta x$ על קפיצים $T/\Delta x$ .
תנאי השפה: קצוות קבועים, מחוץ להישג ידם של הפולסים; מה שקורה כשפולס מגיע לשם הוא נושא מודול 10.
האנסמבל: דטרמיניסטי; רעש גלאי, היכן שהוא מתווסף, הוא גאוסי, בלתי תלוי ובעל זרע קבוע.
מה מוזנח: קשיחות כפיפה, ריסון, שקיעה בכובד, תנועה אורכית ואי-ליניאריות — המיתר גמיש לחלוטין, חסר הפסדים, חסר משקל בשקיעתו וקטן בשיפועיו.
תקף כאשר: השיפועים קטנים, $|y_x| \ll 1$ ; נומרית, מספר קורנט $S = v\,\Delta t/\Delta x \le 1$ וכל מאפיין של פולס משתרע על נקודות רשת רבות.
אופני כשל: שיפועים מסדר גודל אחד, שבהם נשבר המודל ולא הקוד; $S > 1$ , שבו הסכמה הורסת את עצמה מקנה המידה של הרשת ומעלה; נפיצת רשת במאפיינים שלא נפתרו מספיק, שנקראת כפיזיקה.

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

צפו בשרשרת של 5, 20 ו-100 חרוזים הופכת למיתר, ושימו לב כמה לאט הפינות של פריטה משולשת מגיעות לשם.
מדדו זמן של פולס בין שני שערים אופטיים בזמן שהמתיחות משתנה, ואז ערכו מרוץ בין פולסים בגבהים וברוחבים שונים על אותו מיתר.
עקבו אחר גרגר האבק, והשוו כמה מהר הוא נע לכמה מהר הפולס נע.
קבעו צעד זמן ארוך ב-1% מדי וספרו כמה צעדים המיתר שורד.

לגזור את התוצאה¶

דרך (א): השרשרת של מודול 07, צעד אחד הלאה¶

הניחו חרוזים בעלי מסה $m$ על חוט, במרחק $a$ זה מזה, ומתחו את החוט למתיחות $T$ . הזיזו את חרוז $n$ הצידה ב- $y_n$ . קטע החוט שמימינו נוטה בזווית שהסינוס שלה הוא $(y_{n+1} - y_n)/\sqrt{a^2 + (y_{n+1}-y_n)^2}$ , והוא מושך את החרוז בכוח רוחבי של $T$ כפול הסינוס הזה. עבור נטיות קטנות השורש הוא פשוט $a$ , ושני הקטעים יחד נותנים

m\,\ddot y_n = \frac{T}{a}\left(y_{n+1} - y_n\right) - \frac{T}{a}\left(y_n - y_{n-1}\right) = \frac{T}{a}\left(y_{n+1} - 2y_n + y_{n-1}\right).

(1)

זו בדיוק השרשרת של מודול 07, עם קפיצים בקשיחות $k_s = T/a$ — המילון $T = k_s a$ שבו מודול 07 הסתיים, עכשיו עם סיבה מאחוריו.

עכשיו תנו לחרוזים להצטופף. התייחסו להזזות כדגימות של פונקציה חלקה, $y_n(t) = y(na, t)$ , ופתחו את השכנים סביב $x = na$ :

y_{n\pm1} = y \pm a\,y_x + \frac{a^2}{2}\,y_{xx} \pm \frac{a^3}{6}\,y_{xxx} + \frac{a^4}{24}\,y_{xxxx} + \cdots

(2)

בסכום $y_{n+1} + y_{n-1}$ כל איבר אי-זוגי מתבטל, ולכן ההפרש השני הוא

y_{n+1} - 2y_n + y_{n-1} = a^2\,y_{xx} + \frac{a^4}{12}\,y_{xxxx} + \cdots

(3)

השאירו את האיבר המוביל, חלקו ב- $a$ , וכתבו $\mu = m/a$ למסה ליחידת אורך:

\mu\,\frac{\partial^2 y}{\partial t^2} = T\,\frac{\partial^2 y}{\partial x^2}.

(4)

מודול 07 כבר ידע את התשובה לגבי המהירות. יחס הנפיצה שלו $\omega(k) = 2\sqrt{k_s/m}\,|\sin(ka/2)|$ הוא קו ישר ב- $ka$ קטן, בשיפוע $a\sqrt{k_s/m} = \sqrt{(T/a)\,a^2/m} = \sqrt{T/\mu}$ . האיבר $a^4 y_{xxxx}/12$ שהושלך הוא בדיוק מה שמכופף את הקו הזה ב- $k$ גדול יותר: פיתוח הסינוס לסדר הבא נותן $\omega^2 \approx (T/\mu)\,k^2\left(1 - (ka)^2/12\right)$ , אותו שנים-עשר.

דרך (ב): חוק ניוטון לקטע מיתר¶

שכחו מחרוזים. קחו קטע קצר של מיתר רציף, בין $x$ ל- $x + \mathrm{d}x$ , ומנו את הכוחות הפועלים עליו. בכל קצה שאר המיתר מושך בכיוון שלו עצמו במתיחות $T$ . אם המיתר נוטה בזווית $\theta(x)$ בקצה השמאלי וב- $\theta(x + \mathrm{d}x)$ בימני, הרכיבים האנכיים של שתי המשיכות הם $-T\sin\theta(x)$ ו- $+T\sin\theta(x + \mathrm{d}x)$ , והם מתבטלים רק אם המיתר ישר.

עבור שיפועים קטנים, $\sin\theta \approx \tan\theta = y_x$ . הכוח האנכי השקול הוא אז $T\left[y_x(x + \mathrm{d}x) - y_x(x)\right] = T\,y_{xx}\,\mathrm{d}x$ , מסת הקטע היא $\mu\,\mathrm{d}x$ , והחוק השני של ניוטון קורא

\mu\,\mathrm{d}x\;y_{tt} = T\,y_{xx}\,\mathrm{d}x \qquad\Longrightarrow\qquad \mu\,y_{tt} = T\,y_{xx}.

(5)

שיפועים קטנים, בשימוש שלוש פעמים

שתי הדרכים זקוקות ל- $|y_x| \ll 1$ , ודרך (ב) משתמשת בזה שלוש פעמים. היא מחליפה את $\sin\theta$ בשיפוע. היא לוקחת את מסת הקטע כ- $\mu\,\mathrm{d}x$ , כשאורכו האמיתי הוא $\sqrt{1 + y_x^2}\,\mathrm{d}x$ . והיא מתייחסת ל- $T$ כזהה בשני הקצוות: הרכיבים האופקיים של המשיכות, $T\cos\theta$ , חייבים להתאזן מפני שהקטע אינו זז הצידה, ו- $\cos\theta$ נבדל מ-1 רק בסדר שני בשיפוע. דרך (א) משתמשת באותו קירוב פעם אחת, מוסתר במשיכה של כל חרוז.

כמה קטן הוא קטן? למיתר גיטרה שנפרט 3 מ"מ מן הקו באמצע של 650 מ"מ יש שיפועים של כ-0.01, וכל אחד מן התיקונים האלה הוא חלק ב-10⁴. לחבל קפיצה שמסובב בקשת רחבה יש שיפועים מסדר גודל אחד, והמודל הזה פשוט שגוי לגביו. זו טענה על המודל, לא על הפותר — שום כמות של חישוב אינה מצילה מודל שמשמש מחוץ לתחום שבו נגזר.

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

\boxed{\;v = \sqrt{\frac{T}{\mu}}\;}

(6)

והמשוואה מקבלת את הצורה שכל חלק מאוחר יותר של הקורס יזהה:

\frac{\partial^2 y}{\partial t^2} = v^2\,\frac{\partial^2 y}{\partial x^2}.

(7)

פתרון ד’אלמבר: כל פתרון הוא שתי צורות נעות¶

נסו $y = f(x - vt)$ עבור כל פונקציה $f$ גזירה פעמיים. לפי כלל השרשרת, כל נגזרת בזמן מורידה גורם $-v$ וכל נגזרת במרחב גורם 1:

y_{tt} = v^2 f''(x - vt), \qquad y_{xx} = f''(x - vt),

(8)

ולכן $y_{tt} = v^2 y_{xx}$ , תהא $f$ אשר תהא. הגרף של $f(x - vt)$ הוא הגרף של $f$ מוזז ימינה ב- $vt$ : צורה הנעה לעבר $+x$ במהירות $v$ , ללא שינוי. באותו טיעון $g(x + vt)$ היא צורה הנעה לעבר $-x$ . הדבר המופלא הוא שאין שום דבר אחר.

משפט ד’אלמבר

כל פתרון של $y_{tt} = v^2 y_{xx}$ על הישר כולו הוא

y(x,t) = f(x - vt) + g(x + vt)

(9)

עבור פונקציות $f$ ו- $g$ כלשהן: צורה אחת הנעה ימינה, צורה אחת הנעה שמאלה.

ההוכחה היא החלפת משתנים. עם $u = x - vt$ ו- $w = x + vt$ , כלל השרשרת נותן $\partial_x = \partial_u + \partial_w$ ו- $\partial_t = -v\,\partial_u + v\,\partial_w$ , ולכן

y_{tt} - v^2 y_{xx} = v^2\left(\partial_u - \partial_w\right)^2 y - v^2\left(\partial_u + \partial_w\right)^2 y = -4v^2\,\frac{\partial^2 y}{\partial u\,\partial w}.

(10)

משוואת הגלים אומרת אפוא $\partial^2 y/\partial u\,\partial w = 0$ . בצעו אינטגרציה על $w$ : $\partial y/\partial u$ תלויה ב- $u$ בלבד. בצעו אינטגרציה על $u$ : $y$ היא פונקציה של $u$ ועוד פונקציה של $w$ .

קראו מה המשפט משאיר בחוץ. המהירות $v$ נמצאת במשוואה; הצורות $f$ ו- $g$ אינן. התווך בוחר את המהירות, והיד בוחרת רק את הצורה. פולס גבוה יותר הוא אותו פתרון מוכפל בקבוע; פולס חד יותר הוא $f$ אחרת; שניהם נעים במהירות $\sqrt{T/\mu}$ . זה ניבוי 1 — אף אחד אינו מגיע ראשון — וניבוי 3: הכפלת המתיחות מכפילה את $v$ ב- $\sqrt 2$ , ולכן זמן החצייה יורד ל- $1/\sqrt 2$ , כ-71%, ולא לחצי.

בעיית ההתחלה: פריטות והקשות¶

מיתר אמיתי מתחיל עם צורה כלשהי $y_0(x)$ ומהירות כלשהי $w_0(x)$ . אילו $f$ ו- $g$ הן בוחרות? ב- $t = 0$ ,

f(x) + g(x) = y_0(x), \qquad -v\,f'(x) + v\,g'(x) = w_0(x).

(11)

בצעו אינטגרציה למשוואה השנייה, $g - f = \frac{1}{v}\int_0^x w_0(s)\,\mathrm{d}s$ עד כדי קבוע שמתבטל, ופתרו את הזוג:

\boxed{\; y(x,t) = \frac{y_0(x - vt) + y_0(x + vt)}{2} + \frac{1}{2v}\int_{x - vt}^{x + vt} w_0(s)\,\mathrm{d}s \;}

(12)

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

להקשה יש מהירות ואין צורה: פטיש שמכה במיתר פסנתר, מלבן של מהירות $w$ על פני רוחב $2\ell$ . עכשיו רק האינטגרל שורד. בנקודה נתונה הוא סופר כמה מיתר מוכה נמצא בטווח $vt$ ממנה, ולכן המיתר מתרומם לרמה שגדלה בקצב $w$ עד $t = \ell/v$ ואז נשארת בגובה $w\ell/v$ , כשקצותיה רצים החוצה במהירות $v$ . המיתר נשאר מוזז לצמיתות — שום דבר אינו מושך אותו בחזרה, כי הזזת קטע ישר שלם הצידה אינה מותחת דבר.

גלים סינוסיים, ו- $\omega = vk$ ¶

הגל המתקדם של הקורס הוא

y(x,t) = \Real\!\left[A\,e^{\ii(kx - \omega t)}\right] = |A|\cos(kx - \omega t + \arg A).

(13)

הוציאו את $k$ מחוץ לסוגריים: $kx - \omega t = k\,(x - (\omega/k)\,t)$ , ולכן זה $f(x - vt)$ עם $v = \omega/k$ . הוא פותר את משוואת הגלים בדיוק כאשר

\omega = v\,k,

(14)

יחס הנפיצה הישר שמודול 07 מצא בתחתית הסרט של השרשרת. כל אורך גל נע באותה מהירות: זה מה שפירושו שהמיתר האידאלי חסר נפיצה, וזו הסיבה שפולס — סופרפוזיציה של אורכי גל רבים — שומר על צורתו. ביחידות יומיומיות, $v = \omega/k = f\lambda$ .

מה נע, ומה לא¶

עמדו בנקודה אחת, $x_0$ , וצפו בפולס הנע ימינה חולף. ההזזה שם היא $y(x_0, t) = f(x_0 - vt)$ : ככל ש- $t$ רץ, הנקודה קוראת את צורת הפולס לאחור בזמן, עולה כשהחזית מגיעה ויורדת כשהגב עוזב. כשהפולס עבר, $f$ חזרה לאפס, וכך גם הנקודה. זה ניבוי 2: הגרגר מסיים בדיוק במקום שבו התחיל.

מה שנע היה התבנית — עובדת היותו מוזז, שמועברת מכל קטע של מיתר לבא אחריו — וכפי שמודול 09 מראה, האנרגיה שבאה איתה. הקטעים עצמם נעים רק לרוחב, ובמהירות שאין לה קשר ל- $v$ :

y_t = -v\,f'(x - vt) = -v\,y_x .

(15)

המהירות הרוחבית היא $v$ כפול השיפוע, ומכיוון שהשיפוע קטן, המיתר נע לאט בהרבה מן הגל. בפולס של המעבדה, המהירות הגבוהה ביותר של הגרגר היא 1.14 מ’/ש’ בעוד הפולס נע ב-20 מ’/ש’.

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

איך המחשב פותר זאת: השרשרת, נבנית מחדש¶

כדי לבצע אינטגרציה נומרית של משוואת הגלים, דגמו את המיתר בנקודות $x_i = i\,\Delta x$ ובזמנים $t_n = n\,\Delta t$ , והחליפו כל נגזרת שנייה בהפרש המרכזי שלה:

\frac{y_i^{n+1} - 2y_i^n + y_i^{n-1}}{\Delta t^2} = v^2\,\frac{y_{i+1}^n - 2y_i^n + y_{i-1}^n}{\Delta x^2}.

(16)

כשפותרים אותה עבור העתיד, זוהי סכמת leapfrog,

y_i^{n+1} = 2y_i^n - y_i^{n-1} + S^2\left(y_{i+1}^n - 2y_i^n + y_{i-1}^n\right), \qquad S = \frac{v\,\Delta t}{\Delta x},

(17)

כאשר $S$ הוא מספר קורנט. עכשיו הביטו באגף ימין של משוואת ההפרשים לפני החלוקה: הוא $\mu\,\Delta x\;\ddot y_i = (T/\Delta x)\left(y_{i+1} - 2y_i + y_{i-1}\right)$ . זה חוק ניוטון למסה $\mu\,\Delta x$ על קפיצים בקשיחות $T/\Delta x$ — השרשרת של דרך (א), כשהחרוזים במרחק צעד רשת זה מזה. המחשב פותר את משוואת הרצף בכך שהוא בונה מחדש, בשקט, את השרשרת שממנה היא באה. הדיסקרטיזציה מריצה את גבול הרצף לאחור.

תנאי אחד הופך את הסכמה לאמינה. בצעד זמן אחד העדכון בנקודה $i$ מתייעץ רק בנקודות $i-1$ , $i$ ו- $i+1$ , כך שמידע על הרשת נע לכל היותר מרווח אחד בכל צעד. הפתרון האמיתי ב- $x_i$ תלוי בכל מה שנמצא בטווח $v\,\Delta t$ ממנו. אם $v\,\Delta t > \Delta x$ , הפתרון תלוי בנקודות שהסכמה מעולם לא הביטה בהן, ושום סידור של מספרים לא יוכל לקלוע. תנאי קורנט–פרידריכס–לוי הוא

\boxed{\; S = \frac{v\,\Delta t}{\Delta x} \le 1 \;}

(18)

והפרתו אינה אי-דיוק קטן. הסעיף הבא מראה מה היא כן.

לאמת חישובית¶

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

1. הפותר מול ד’אלמבר. פריטה גאוסית על מיתר באורך 4 מ’ מקודמת במספר קורנט קבוע $S = 0.5$ ומושווית, בכל נקודה, לפתרון המדויק. חיצוי מרווח הרשת — ואיתו הצעד — מקטין את שגיאת $L^2$ בגורמים 3.996, 4.002, 4.001 ו-4.000, והמעריך המותאם הוא -2.000: הסכמה מסדר שני, כפי שההפרשים המרכזיים מבטיחים. מבחן ההתכנסות הזה רץ לפני שכל אנימציה בעמוד הזה מרונדרת.

2. הסטנסיל הוא השרשרת. ל-simulate_string נמסר מיתר בן 201 נקודות, ול-coupled.simulate_coupled — האינטגרטור של מודול 07, שאינו יודע דבר על רשתות — נמסרת שרשרת המסות $\mu\,\Delta x$ על קפיצים $T/\Delta x$ . בהתנעה זהה, הם מסכימים לאורך 2000 צעדים עד כ- $1\times 10^{-14}$ מן ההזזה הגדולה ביותר. הטענה החותמת של הגזירה אינה אנלוגיה.

3. הצעד הקסום. ב- $S = 1$ בדיוק, הסטנסיל מתפשט ל- $y_i^{n+1} = y_{i+1}^n + y_{i-1}^n - y_i^{n-1}$ , שמזיז כל צורה נעה במרווח רשת אחד בדיוק בכל צעד.

מדויק ב-S = 1, ורק שם

פריטה שמקודמת 400 צעדים ב- $S = 1$ תואמת את ד’אלמבר עד $5.9\times 10^{-15}$ מגובהה — שגיאת עיגול. אותה פריטה ב- $S = 0.99$ , על אותה רשת, שגויה ב- $1.0\times 10^{-5}$ : עשרה סדרי גודל עבור שינוי של אחוז אחד בצעד.

היא מדויקת רק בהינתן ערכים מדויקים בשני הצעדים הראשונים. פריטה מספקת אותם, מפני שלמיתר במנוחה יש צעד שני מדויק. פולס שהותנע בתנועה אינו מספק: הצעד השני שלו שגוי בסדר שלישי, ההעברה נושאת את הפגם הזה בנאמנות, והשגיאה מתייצבת בסדר שני — עדיין קטנה כפי 7.5 מזו שב- $S = 0.5$ על אותה רשת, אבל כבר לא מדויקת. הטריק אמיתי וצר.

4. מעבר לגבול. הפריטה מורצת ב- $S = 1.01$ .

ההתפוצצות, וקצבה

ההזזה עוברת פי עשרה מגובה הפריטה בצעד 140 וממשיכה להתרבות. הגידול אינו אקראי: הוא שוכן בגל הקצר ביותר שהרשת יכולה להחזיק, שכנים בפאזה הפוכה בדיוק, והוא מתחיל משגיאת העיגול כבר בצעד הראשון. הגידול הנמדד בכל צעד תואם את ניבוי פון נוימן שנגזר בחלק המתקדם עד פחות מ-0.5% ב- $S = 1.01$ , 1.05 ו-1.1. ב- $S = 0.99$ , אלפיים צעדים משאירים את גובה הפולס במקומו.

המיתר יציב בכל $S$ . מה שנכשל הוא סכמה שנדרשה להעביר מידע רחוק יותר ממרווח רשת אחד בכל צעד — וזו הסיבה ש-simulate_string מסרב ל- $S > 1$ אלא אם נאמר לו שהצפייה בכישלון היא העניין.

5. שערים אופטיים, ומה הם מודדים. שני גלאים רושמים את המיתר במקום שבו הם עומדים, והגעת הפולס לכל אחד נמדדת לפי מרכז הכובד של הרישום שלו. כשצעד הזמן קבוע עבור המיתר המהיר ביותר, המתיחות נסרקת לאורך עשור: המעריך המותאם של $v$ כנגד $T$ הוא 0.5000000, וכנגד $\mu$ הוא -0.5000000, וכל מהירות נמצאת במרחק של פחות מ- $3\times 10^{-9}$ מ- $\sqrt{T/\mu}$ . תשעה פולסים בשלושה גבהים ושלושה רוחבים על מיתר אחד מגיעים באותה מהירות עד 10^-7 — המפריך של “ניעור חזק יותר מאיץ אותו”, שהורץ תשע פעמים.

מדידת זמן לפי שיא כל רישום קוראת במקום זאת מהירות נמוכה מדי, בשיעור היחסי $(1 - S^2)(\Delta x/\sigma)^2/4$ עבור גאוסיאן ברוחב $\sigma$ — נמדד בטווח של 0.4% מן הנוסחה הזו ב- $S = 0.25$ , 0.5 ו-0.75, ובדיוק $v$ ב- $S = 1$ . זו נפיצת רשת: על הרשת, אורכי גל קצרים נעים מעט לאט יותר מארוכים, ושיא הפולס נגרר לאחור. מרכז הכובד חסין, מפני שהוא מגיב רק לשיפוע יחס הנפיצה ב- $k = 0$ , שם הרשת מדויקת. חלק המדידה של המעבדה הופך את אותה הבחנה לסרגל שגיאה.

6. אנרגיה. לאורך שש-עשרה חציות בין שני קצוות קבועים, האנרגיה הכוללת של הפותר מתנדנדת ב- $4.7\times 10^{-4}$ מעצמה ב- $S = 0.5$ ואינה נסחפת; חצו את הצעד והתנודה יורדת פי ארבעה. מהי האנרגיה הזו, והיכן על המיתר היא שוכנת, זה מודול 09.

להעביר את הרעיון¶

משוואה אחת, תווכים רבים

הגזירה שלמעלה מעולם לא נזקקה לכך שהמיתר יהיה מיתר. היא נזקקה לגודל שאפשר להזיז, להתמד שמתנגד להזזה, ולהשפעה מחזירה ביחס לעקמומיות. בכל מקום ששלושת אלה נפגשים, $y_{tt} = v^2 y_{xx}$ נובעת, כש- $v^2$ הוא הקשיחות המחזירה מחולקת בהתמד:

קול בגז: הפרעות לחץ, $v = \sqrt{\gamma p/\rho}$ — 343 מ’/ש’ באוויר בטמפרטורת החדר.
גלים סייסמיים: גלי גזירה בסלע עם $v = \sqrt{G/\rho}$ , וגלי דחיסה, מהירים יותר, עם $v = \sqrt{(K + \tfrac43 G)/\rho}$ ; הפער בין זמני ההגעה שלהם הוא הדרך שבה סייסמוגרף מודד את מרחקו מרעידת אדמה.
מים רדודים: צונאמי באוקיינוס בעומק 4000 מ’ נע ב- $\sqrt{g h}$ , כ-200 מ’/ש’ — מהירות של מטוס נוסעים.
אור: משוואות מקסוול נותנות את אותה משוואה עבור השדה החשמלי, עם $v = 1/\sqrt{\mu_0\varepsilon_0} = c$ . מודול 14 גוזר זאת.

גלים עומדים. הציבו קירות בשני הקצוות והחצאים הנעים חוזרים. 11-standing-waves בונה את אופני המיתר משני הכיוונים של ד’אלמבר, ומשחזר את $\sin(p\pi x/L)$ של מודול 07 כגבול $N \to \infty$ של צורות האופנים של השרשרת.
נפיצה. המיתר האידאלי חסר נפיצה מעצם בנייתו: $\omega = vk$ הוא קו ישר. 13-dispersion שואל מה קורה כשאינו כזה, שם פולס מתפשט בדרכו — וכבר ראיתם זאת קורה, על הרשת, בבדיקת השערים האופטיים שלמעלה.
בחזרה למודול 07. קצה הסרט של השרשרת היה תקרה שאין למיתר. משוואת המיתר היא מה שנשאר מן השרשרת כשכל אורך גל שמעניין אותנו ארוך בהרבה מן המרווח, ודרך (א) אומרת בדיוק כמה ארוך.
גלים אורכיים. השרשרת של מודול 07 צוירה לרוחב אבל מודלה לאורך: מסות המחליקות לאורך הקו שלהן. כשקוראים אותה כך, דרך (א) מדויקת בלי שום הנחת שיפוע קטן, ונותנת קול במוט מוצק, כש- $k_s a$ ממלא את תפקיד הקשיחות של המוט ולא של מתיחות.

בדקו את הבנתכם¶

שני פולסים נשלחים לאורך אותו חבל מתוח, זה אחר זה. השני נוצר בהצלפה חזקה פי שניים בקצה, ולכן הוא גבוה פי שניים, ומהירה פי שניים, ולכן הוא צר בחצי. איך המהירות שלו לאורך החבל משתווה למהירות הפולס הראשון?
A. זהה. המתיחות של החבל והמסה ליחידת אורך שלו קובעות את המהירות; ההצלפה קובעת רק את הצורה. B. מהירה יותר. הפולס הגבוה והחד נושא יותר אנרגיה ולכן עוקף את הראשון. C. איטית יותר. פולס חד יותר מותח את החבל יותר, ואת המתיחה הנוספת צריך למשוך בחזרה. D. מהירה יותר, פי שניים, שהרי היד נעה מהר פי שניים.
בדקו את התשובה
התשובה הנכונה: A — זהה. המתיחות של החבל והמסה ליחידת אורך שלו קובעות את המהירות; ההצלפה קובעת רק את הצורה.
- A. נכון. משוואת הגלים mu y_tt = T y_xx מכילה את T ואת mu ולא שום דבר אחר, ולכן כל פתרון שלה נע במהירות v = sqrt(T/mu) תהא צורתו אשר תהא. הגובה אינו יכול להשפיע מפני שהמשוואה ליניארית: פולס גבוה יותר הוא אותה תנועה מוכפלת בקבוע. הרוחב אינו יכול להשפיע מפני שפתרון ד’אלמבר f(x - vt) מזיז כל f שהיא, חדה או רחבה. במדידה על המיתר המדומה, תשעה פולסים בשלושה גבהים ובשלושה רוחבים מגיעים באותה מהירות עד כדי חלק אחד לעשרה מיליון.
- B. הוא אכן נושא יותר אנרגיה, ויותר ממנה מגיע רק מפני שיותר ממנה נשלח, לא מפני שהיא נעה מהר יותר. אנרגיה ומהירות הן כאן שאלות נפרדות: המהירות נקבעת על ידי התווך, ומודול 09 מראה שהאנרגיה נישאת אז באותה מהירות בדיוק. הצלפה חזקה יותר משנה את מה שמגיע, לעולם לא את מתי.
- C. המתיחה הנוספת אמיתית, אבל היא מסדר שני בשיפוע, והמודל שנותן v = sqrt(T/mu) כבר השליך אותה כשהניח שהשיפוע קטן. בתוך המודל הזה המתיחות זהה בכל מקום, וכך גם המהירות. פולס תלול מספיק כדי שהמתיחה תשפיע הוא פולס שהמיתר האידאלי כבר אינו מתאר.
- D. מהירות היד קובעת כמה מהר כל חלק של החבל נע מעלה ומטה — המהירות הרוחבית — ולא כמה מהר ההפרעה נעה לאורכו. השתיים בלתי תלויות: פולס מדומה על מיתר של 20 m/s הזיז את חלקיו ב-1.1 m/s. האיצו את היד והחלקים ינועו מהר יותר לרוחב בעוד הפולס נע באותם 20 m/s בדיוק.
גרגר אבק קטן מונח על מיתר אופקי ארוך. פולס יחיד בצורת גבשושית נע לאורך המיתר, עובר את הגרגר, וממשיך הלאה. היכן הגרגר אחרי שהפולס עבר?
A. בדיוק במקום שבו התחיל. הוא עלה וירד כשהפולס עבר, ושום דבר לא נשא אותו הלאה. B. מעט הלאה, בכיוון שבו הפולס נע. C. מורם מעט, מפני שהפולס הרים אותו ושום דבר לא דחף אותו בחזרה למטה. D. נישא עם הפולס עד קצה המיתר.
בדקו את התשובה
התשובה הנכונה: A — בדיוק במקום שבו התחיל. הוא עלה וירד כשהפולס עבר, ושום דבר לא נשא אותו הלאה.
- A. נכון. במקומו של הגרגר ההזזה היא y = f(x0 - vt): היא עוקבת אחר צורת הפולס בזמן, עולה ויורדת, וחוזרת ל-f = 0 אחרי שהפולס עבר. הפולס הוא תבנית שמועברת מכל חלק של המיתר לבא אחריו, והחלקים נשארים בבית. בסימולציה הגרגר עלה 10 mm וסיים במרחק של פחות ממיליונית מכך מן האפס.
- B. זו התמונה שצריך לוותר עליה: שגל הוא חומר שנע. הגרגר מחובר לחלק אחד של המיתר, והחלק הזה נמשך למעלה על ידי שכנו מצד אחד ובחזרה למטה על ידי שכנו מצד שני. אף אחת מן המשיכות אינה מכוונת לאורך המיתר, ולמיתר האידאלי אין שום תנועה לאורכו. מה שנע הוא הצורה — כמו שגל ביציעים נע סביב האצטדיון בזמן שכל אוהד נשאר במושבו.
- C. משהו כן דחף אותו בחזרה למטה: החצי האחורי של הפולס. כשהשפה האחורית עוברת, השיפוע של המיתר מתהפך ושכניו מושכים אותו אל הקו. מיתר שנשאר מוזז לצמיתות אכן קיים — הקשה משאירה רמה — אבל לשם כך צריך מהירות פתיחה, ופולס שפשוט חולף אינו מעביר כזו במצטבר.
- D. אילו הפולס היה נושא את המיתר איתו, פולס שנשלח לאורך חבל היה צובר חבל ליד הקיר הרחוק, והחבל היה מתקצר אחרי כל הצלפה. מה שמגיע לקיר הוא הפרעה והאנרגיה שלה. החבל נשאר במקומו.
מיתר ארוך ומתוח נמשך למשולש בגובה 12 mm באמצעו, מוחזק בלי תנועה, ומשוחרר. אחרי שלתנועה היה זמן להיפרד, איך המיתר נראה?
A. שני משולשים, כל אחד בגובה 6 mm, הנעים זה מזה במהירות הגל. B. משולש אחד בגובה 12 mm, הנע בכיוון אחד במהירות הגל. C. משולש בגובה 12 mm שנשאר באמצע ומתנדנד מעלה ומטה במקומו. D. שני משולשים בגובה 12 mm, הנעים זה מזה.
בדקו את התשובה
התשובה הנכונה: A — שני משולשים, כל אחד בגובה 6 mm, הנעים זה מזה במהירות הגל.
- A. נכון. פתרון ד’אלמבר לצורה ששוחררה ממנוחה הוא y = [y0(x - vt) + y0(x + vt)]/2: חצי מן הצורה רץ ימינה, חצי רץ שמאלה. ב-t = 0 שני החצאים מונחים זה על זה ומסתכמים למקור. הוא חייב להתפצל, מפני שלמיתר ששוחרר ממנוחה אין מהירות שתעדיף כיוון כלשהו, והוא חייב להתפצל שווה בשווה מאותה סיבה.
- B. משולש יחיד שנע ימינה זקוק ל-y_t = -v y_x כבר מן הרגע הראשון: חזית המשולש כבר עולה, גבו כבר יורד. למיתר שהוחזק בלי תנועה ושוחרר יש y_t = 0 בכל מקום, ולכן הוא אינו יכול להיות צורה נעה יחידה. מהירות הפתיחה בוחרת את הכיוון; בלעדיה, כל כיוון מקבל חצי.
- C. זה מה שמסה יחידה על קפיץ הייתה עושה, וזו התמונה הנכונה למיתר כולו רק כששני הקצוות קרובים מספיק כדי להחזיר את החצאים ישר בחזרה — גלים עומדים, מודול 11. על מיתר ארוך שני החצאים פשוט עוזבים.
- D. הכיוון נכון והגובה לא. שני משולשים בגובה מלא היו מכילים פי שניים מן ההזזה שהמיתר התחיל איתה, והמשוואה ליניארית: ב-t = 0 שני החלקים הנעים חייבים להסתכם בדיוק לצורה המקורית, ולכן כל אחד הוא חצי.
חבל באורך 12 מ’ בעל מסה כוללת של 0.60 kg מוחזק במתיחות של 80 N. כמה זמן לוקח לפולס לעבור מקצה אחד לשני, בשניות?
בדקו את התשובה
התשובה הנכונה: 0.3 (± 0.01)
המסה ליחידת אורך היא mu = 0.60 kg / 12 m = 0.050 kg/m, ולכן v = sqrt(T/mu) = sqrt(80 / 0.050) = sqrt(1600) = 40 m/s, והחצייה לוקחת 12 m / 40 m/s = 0.30 s. המסה הכוללת נכנסת רק דרך mu: הכפילו את אורך החבל באותו mu והמהירות לא תשתנה בעוד הדרך תיקח פי שניים. בדקו את היחידות תוך כדי — N/(kg/m) הוא (kg m/s^2)(m/kg) = m^2/s^2, שהשורש שלו מהירות.
מיתר עם מתיחות 50 N וצפיפות מסה אורכית 5.0 g/m מדומה בסכמת leapfrog על רשת במרווח 2.0 mm. מהו צעד הזמן הגדול ביותר, במיקרו-שניות, שעבורו הסכמה יציבה?
בדקו את התשובה
התשובה הנכונה: 20.0 (± 0.5)
קודם מהירות הגל, ביחידות SI: mu = 5.0 g/m = 0.0050 kg/m, ולכן v = sqrt(50 / 0.0050) = 100 m/s. תנאי קורנט S = v dt / dx <= 1 נותן אז dt <= dx / v = 0.0020 m / 100 m/s = 2.0e-5 s = 20 מיקרו-שניות. פיזיקלית: בצעד אחד, שום דבר אינו יכול לנוע יותר ממרווח רשת אחד, מפני שכל נקודה מתעדכנת רק משכנותיה הקרובות. הריצו ב-1% מעל זה והגל הקצר ביותר שהרשת מחזיקה גדל פי 1.33 בכל צעד, משגיאת העיגול, עד שהוא מציף הכול.
שתי תמונות של מיתר שצולמו בהפרש 0.10 s מראות את אותה גבשושית, בלי שינוי בצורתה, שהשיא שלה זז מ-x = 3.0 m ל-x = 2.4 m. איזה תיאור מתאים?
A. y = g(x + vt) עם v = 6.0 m/s — נע שמאלה. B. y = f(x - vt) עם v = 6.0 m/s — נע ימינה. C. y = f(x - vt) עם v = -6.0 m/s. D. y = g(x + vt) עם v = 0.6 m/s.
בדקו את התשובה
התשובה הנכונה: A — y = g(x + vt) עם v = 6.0 m/s — נע שמאלה.
- A. נכון. השיא זז 0.6 m לכיוון x- ב-0.10 s, ולכן התבנית נעה ב-6.0 m/s בכיוון x-. פונקציה של x + vt שומרת על ערכה כש-x קטן בזמן ש-t גדל, וזו בדיוק התנועה הזו: השיא יושב היכן ש-x + vt קבוע. כש-v חיובי, הסימן בתוך הארגומנט נושא את הכיוון ו-v נושא רק את המהירות.
- B. פונקציה של x - vt שומרת על ערכה היכן ש-x - vt קבוע, ולכן השיא שלה זז ל-x גדול יותר ככל ש-t גדל. השיא הזה זז ל-x קטן יותר. הסימן בתוך הארגומנט הוא הכיוון; הפכו אותו וכל בעיית החזרה במודול 10 תצא הפוכה.
- C. אלגברית זו אותה פונקציה כמו g(x + vt) עם v = 6.0 m/s, ולכן זה לא כל כך שגוי כמו מלכודת חשבונאית. הקורס שומר את v = sqrt(T/mu) חיובי, כתכונה של התווך, וכותב את הכיוון בארגומנט. “מהירות” שלילית פירושה גודל שלפעמים הוא מהירות ולפעמים לא.
- D. הכיוון נכון והמהירות שגויה בגורם של פרק הזמן. מהירות היא מרחק חלקי זמן: 0.6 m ב-0.10 s הם 6.0 m/s. מרחק של 0.6 m היה לוקח שנייה שלמה ב-0.6 m/s.
אל משוואת הגלים d^2y/dt^2 = v^2 d^2y/dx^2 אפשר להגיע משרשרת המסות של מודול 07 ומחוק ניוטון לקטע קטן של מיתר. תארו את שתי הדרכים במילים, אמרו היכן הן נפגשות, וציינו את הקירוב ששתיהן זקוקות לו.
בדקו את התשובה
דרך השרשרת מתחילה מ-N חרוזים בעלי מסה m, במרווח a, על חוט במתיחות T. כל חרוז נמשך אל שכניו, ועבור שיפועים קטנים המשיכה הרוחבית של כל קטע חוט היא T כפול הפרש הגבהים חלקי a: קפיץ בקשיחות T/a. זה נותן את משוואת מודול 07, m y_n’’ = (T/a)(y_{n+1} - 2 y_n + y_{n-1}). עכשיו תנו לחרוזים להצטופף: כתבו y_n = y(na, t) ופתחו את השכנים לטור טיילור סביב x = na. הצירוף y_{n+1} - 2 y_n + y_{n-1} הופך ל-a^2 כפול הנגזרת השנייה ב-x, ולכן y_tt = (T a/m) y_xx, ועם mu = m/a זה y_tt = (T/mu) y_xx. יחס הנפיצה של השרשרת אומר את אותו דבר באופן בלתי תלוי: ב-ka קטן השיפוע שלו הוא a sqrt(k_s/m), שהוא sqrt(T/mu).

דרך המיתר לעולם אינה מזכירה חרוזים. קחו קטע מיתר בין x ל-x + dx. המתיחות מושכת לאורך המיתר בכל קצה, ואם המיתר נוטה באופן שונה בשני הקצוות הרכיבים האנכיים של המשיכות האלה אינם מתבטלים. עבור שיפועים קטנים הסינוס של זווית הנטייה הוא השיפוע עצמו, ולכן הכוח האנכי השקול הוא T כפול השינוי בשיפוע לאורך הקטע, T y_xx dx. חוק ניוטון למסה שלו mu dx נותן mu y_tt = T y_xx. הן נפגשות באותה משוואה עם אותה מהירות, v = sqrt(T/mu). ההסכמה הזו היא העיקר: משוואת הגלים אינה עוסקת בחרוזים או במיתרים, אלא בכל תווך שחלקיו נמשכים בחזרה אל הקו על ידי שכניהם ביחס למרחקם מן הקו. שתי הדרכים זקוקות לשיפועים קטנים, |dy/dx| << 1. בדרך השרשרת זה מוסתר במשיכה של החרוז, שם sin(theta) הוחלף בהפרש הגבהים חלקי a; בדרך המיתר זה מחליף את sin(theta) בשיפוע, מאפשר לאורך הקטע להיות dx, ושומר על המתיחות זהה בשני הקצוות. פולס תלול שובר את המודל, לא את המתמטיקה.

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

להסביר במילים שלכם¶

ענו בכמה משפטים על כל אחת. שום מספר לא יציל אתכם.

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

מתקדם: יציבות, הצעד הקסום, דמויות, ומהי משוואה דיפרנציאלית חלקית ראשונה¶

הטיעון של פון נוימן. מאין בדיוק מגיע $S \le 1$ ? נסו אופן פורייה יחיד על הרשת, $y_i^n = g^n e^{\ii k i \Delta x}$ , ושאלו כיצד המשרעת שלו $g$ משתנה בכל צעד. הצבה בסטנסיל leapfrog נותנת

g^2 - 2\left[1 - 2S^2\sin^2\!\left(\frac{k\Delta x}{2}\right)\right] g + 1 = 0 .

(19)

מכפלת שני השורשים היא 1. כל עוד הסוגריים נמצאים בין -1 ל-1 הם זוג צמודים מרוכבים על מעגל היחידה, $|g| = 1$ , והאופן אינו גדל ואינו דועך. ברגע שהסוגריים יורדים מתחת ל--1 הם ממשיים, לאחד מהם $|g| > 1$ , והאופן גדל מעריכית. הסוגריים נשארים בתחום לכל $k$ אם ורק אם $S \le 1$ . האופן הגרוע ביותר הוא הקצר ביותר שהרשת יכולה להחזיק, $k\Delta x = \pi$ , שבו עבור $S > 1$

|g| = \left(2S^2 - 1\right) + \sqrt{\left(2S^2 - 1\right)^2 - 1} ,

(20)

1.3266 ב- $S = 1.01$ . זה המספר שההתפוצצות בסעיף האימות תאמה — ומכיוון ששגיאת עיגול מכילה כל $k$ בכמויות קטנות, תמיד יש לה זרע לגדול ממנו.

מדוע S = 1 מדויק. עם $|g| = 1$ , כתבו $g = e^{-\ii\omega\Delta t}$ ואותה משוואה הופכת ליחס הנפיצה של הרשת,

\sin\!\left(\frac{\omega\Delta t}{2}\right) = S\,\sin\!\left(\frac{k\Delta x}{2}\right).

(21)

ב- $S = 1$ שני הסינוסים שווים, ולכן $\omega\Delta t = k\Delta x$ ו- $\omega = vk$ בדיוק — כל אורך גל שהרשת יכולה להחזיק נע במהירות $v$ בדיוק, ואין נפיצת רשת כלל. עבור $S < 1$ , פיתוח שני הסינוסים נותן $\omega \approx vk\left[1 - (1 - S^2)(k\Delta x)^2/24\right]$ : גלים קצרים מפגרים, וזו ההטיה של השערים האופטיים שנמדדה למעלה. הצעד המדויק והצעד הלא יציב רחוקים $1\%$ זה מזה.

חצי הישר, בדמויות. פתרון ד’אלמבר חי על ישר אינסופי. הציבו קצה קבוע ב- $x = 0$ , שבו $y(0, t) = 0$ לכל זמן, והוא עדיין עובד, בזכות תעלול. הרחיבו את צורת הפתיחה אל $x$ שלילי כפונקציה אי-זוגית, $y_0(-x) = -y_0(x)$ , ופתרו על הישר כולו. ב- $x = 0$ שני החצאים שווים ומנוגדים תמיד, ולכן ההזזה שם אפס לנצח — תנאי השפה, מתקיים מכוח סימטריה ולא נכפה. פולס שמתקרב לקיר פוגש את הדמות ההפוכה שלו שבאה מהכיוון השני, ואחרי שהם עוברים זה דרך זה הדמות היא שיוצאת: קצה קבוע מחזיר פולס הפוך. הרחיבו באופן זוגי במקום זאת והשיפוע ב- $x = 0$ מתאפס — קצה חופשי — שמחזיר את הפולס זקוף. מודול 10 גוזר את שתי התוצאות כהלכה, כשני הגבולות של נוסחה אחת.

מה הופך משוואה דיפרנציאלית חלקית ל"ראשונה". משוואת הגלים היא המשוואה ההיפרבולית הפשוטה ביותר, ואופייה נמצא במשתנים של ד’אלמבר: $x - vt$ ו- $x + vt$ קבועים לאורך קווים במישור $(x, t)$ שנקראים אופייניים, ומידע נע רק לאורכם. הפתרון ב- $(x, t)$ תלוי בנתוני הפתיחה רק בקטע $[x - vt, x + vt]$ — תחום התלות שלו — ובשום דבר מחוצה לו. תנאי CFL הוא ההד הבדיד של העובדה הזו: תחום התלות של הרשת עצמה, מרווח אחד בכל צעד, חייב להכיל את האמיתי. משוואת החום, $y_t = D\,y_{xx}$ , נראית כמעט זהה ומתנהגת אחרת לגמרי: אין לה אופייניים, הפרעה מורגשת בכל מקום מיד, וצורות נמרחות במקום לנוע. נגזרת אחת בזמן במקום שתיים היא כל ההבדל.