כמהלכי מים עיליים בזרם, נהר, או אזור גאות ושפל הוא יוצר הדרגתיים ראש המניע את המים לתוך ומחוץ למשקעי 1. במערכות נהרות החלק ממשקעי הנחל שבו חילופי דברים מתרחש ידוע כ2,3 אזור hyporheic. אזור זה חשוב כי חומרים מזינים רבים ומזהמים מאוחסנים, הופקדו, או הפכו בתוך אזור hyporheic 4-9. כמות הזמן נותב מבלה במשקעים נקראת זמן מגורים. בשתי פעמים המגורים ואת מיקומם של מסלולי הזרימה משפיעים על תהליכי השינוי. יש צורך בהבנה משופרת של התהליכים המשפיעים על זרימה דרך המשקעים לחזות תחבורה מומסת בנהרות ולטפל בבעיות סביבתיות גדולות כתוצאה מהתפשטות של חומרים כגון חומרים מזינים (למשל, היפוקסיה החוף 10,11). למרות החשיבות של חילופי hyporheic, הוא לעתים קרובות אינם מתואר בקורסים לתואר ראשון בהידרולוגיה,מכניקת זורמים, הידראוליקה, וכו 'מחנכים המבקשים להוסיף חילופי hyporheic לקורסים שלהם יכולים למצוא את זה שימושי לי חזותי ניסיוני ומספרי שמראה בבירור את התהליך הזה.
פיתולי נחל ערוץ, רמות תהום שמסביב, וטופוגרפיה הנחל (כלומר, בארים, bedforms, ותלים הביוגניים) כל להשפיע חילופי hyporheic בדרגות שונות 12-17. מחקר זה התמקד בbedforms, כגון דיונות ואדוות, אשר בדרך כלל תכונות עיקריות המשפיעות על זרימת גאומורפיות hyporheic 14,15. יצרנו ניסוי סימולציה ומעבדה מספרי כדי להמחיש זרימה דרך סדרה קבועה של bedforms. סימולציה זה מבוססת על גוף של מחקרים קודמים הנוגע נתיבי זרימת hyporheic מאפייני מערכת לצפייה בקלות 15,18-21. כמו מחקר זה מהווה את הרקע המדעי לסימולציה, סיכום קצר של ההיבטים המרכזיים של התאוריה הבאה. הטופוגרפיה Bedform, T (x),ניתן על ידי:
משוואת 1:

כאשר H הוא פעמיים את המשרעת של bedform, k הוא wavenumber, וx הוא מקביל ממד האורך אל פני השטח הנחל הממוצע. דוגמא של הטופוגרפיה bedform זה מוצגת באיור 1.

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

| שם פרמטר | יחידות | הַגדָרָה | מִמְשָׁק | Mousedrop |
| למבדה (λ) | ס"מ | אורך גל של bedform (ראה איור 1) |  |  |
| BedformHeight (H) | ס"מ | פעמיים משרעת bedform (ראה איור 1) |  |  |
| BedDepth (D) | ס"מ | עומק של המשקעים (ראה איור 1) |  |  |
| HydrCond (K) | סנטימטר / s | מוליכות הידראוליות |  |  |
| נקבוביות (θ) | | נַקבּוּבִיוּת |  |  |
| ChannelVelocity (U) | סנטימטר / s | מהירות ממוצעת במים העיליים או ערוץ |  |  |
| עומק (ד) | ס"מ | עומק מים (ראה איור 1) |  |  |
| מדרון (S) | | שיפוע של bedforms ומשטח מים |  |
| NumParticles | | מספר החלקיקים שוחרר למערכת. |  | |
| טיימקס (TIME1, TIME2 ..) | דקות | זמן שבו כל שינוי צבע מתרחש | |  |
| לחצני סימולציה | | הַגדָרָה | מִמְשָׁק | Mousedrop |
| להכין | | הגדר של את הסימולציה באמצעות פרמטרים המוצגת |  |  |
| ללכת / להפסיק | | מתחיל ומפסיק את הסימולציה |  |  |
| שָׁלָב | | לחיצה צעד צעד גורמת פעם אחת כדי לעבור. זה מאפשר למשתמשים להאט את הקוד ולראות בדיוק מה קורה ב100 שניות. |  | |
| שבילים ברורים | | מנקה את כולו נתיבי חלקיקים כחולים מהמסך |  |  |
| להתקדם לפעם הבאה | | זה גורם לתכנית לפעול עד לשינוי צבע בפעם הבא (Timex) | |  |
| עכבר-טיפה | | כפתור זה חייב להיות לחץ לפני חלקיקים עשויים להיות ממוקמים בתת הקרקע על ידי לחיצה מקומות בתת הקרקע ב. | |  |
| נתיבים-הופעה? | | אם נתיבים-הופעה? הוא "על" חלקיקי המים להשאיר עקבות של הקרנה הכחולה שבו הם היו (ראה איור 1). |  |  |
| מחדש לרדת? | | אם מחדש ירידה? הוא "על" החלקיקים מוחלפים באופן משוקלל שטף לכל חלקיק, שיוצא למערכת, והעלילה מצטברת לא עובדת. כאשר Parti CLE יוצא אזור hyporheic מספר החלקיקים במערכת יורד אם מחדש ירידה? הוא "מעל" (ראה איור 1). |  | |
בקרת הטבלה 1. פרמטרי Hyporheic וסימולציה. כל פרמטר, לחצן, ומחוון שיכול להיות מותאם על ידי המשתמש מקבלים בטבלה זו יחד עם הגדרה.
בסימולציה זו, שני תהליכים לגרום למהירות נוזל במיטת החול. הראשון הוא בשל האינטראקציות של זרימת הזרם עם bedforms. ראש המהירות בממשק מים / משקעים הנגרמים על ידי bedforms גם כ סינוסי, ועבר על ידי גל רבעון מbedform עצמו 22. המשרעת של פונקצית ראש המהירות בממשק משטח-מתחת לפני הקרקע כבר מקורבות ממדידות 16:
GE = "תמיד"> משוואה 2:

בי U היא מהירות מים העילית הממוצעת, גרם הוא קבוע הכבידה, וד הוא העומק של המים (מוצגים באיור 1). פונקצית ראש המהירות לאחר מכן ניתנה על ידי:
משוואה 3:

פונקצית ראש זה לאחר מכן ניתן להשתמש כדי לחשב את הרכיב מבוסס bedform של פונקציות מהירות מתחת לפני הקרקע על ידי פתרון משוואת Laplace עם עומק מיטת חול קבוע 20. המרכיב השני של מהירות porewater נקבע על ידי השיפוע של המערכת, S, אשר תואמת את שיפוע ראש הכבידה שתשואות לזרום בכיוון מורד הזרם הפרופורציונלי לs / ftp_upload / 53,285 / 53285eq_S_inline.jpg "/> הפונקציות הסופיות למהירות porewater הן.:
משוואה 4:

משוואה 5:

איפה u הוא רכיב המהירות האורך, נ 'הוא רכיב המהירות האנכי, K הוא המוליכות הידראוליות הממוצעת של המשקעים, הוא נקבוביות הממוצעת של המשקעים, y הוא לתאם האנכי, ו- D הוא העומק של המשקעים.
סימולציות מעקב חלקיקים נוצרו, שמשתמשות בשפת דוגמנות NetLogo ופלטפורמת סימולציה 23. שני יישומים (Mousedrop.nlogo וInterface.nlogo) להשתמש משוואות אלה למודל HYPזרימת orheic עם אותו ליבת הסימולציה. ההבדל העיקרי הוא במקומות הראשוניים של החלקיקים נותב. Mousedrop מאפשר למשתמש למקם נותב מדומה בכל מקום בתוך תת הקרקע. מהירות מתחת לפני הקרקע משוואות 4 ו -5 משמשות כדי להזיז את נותב כדי לדמות ניסויי הזרקת צבע. בממשק, נותב תמיד להציב לאורך גבול משטח / מתחת לפני הקרקע באופן משוקלל שטף. זה מחקה את המשלוח של חומר מומס והושעה מהמים העיליים לporewater, שהוא חיוני להבנת חילופי hyporheic. נותב אז נע בתוך תת הקרקע עד ששוב מגיע למים הנחל. איתור נתיבי הצבע במתעל ומדמה את הנתיבים באמצעות NetLogo מניב המייעל של flowfield, כל עוד תנאי הזרימה ומורפולוגיה bedform יישארו יציבים בתקופה של התבוננות. Interface.nlogo יוצר חלוקת זמן מגורים מצטברת, אשר מציגה את יחס של מספרחלקיקים נותב שנותרו במשקעים למספר הראשוני של חלקיקים נותב הוצבו בזמן 0 כפונקציה של זמן.
כפי שנאמר בסקר הספרות האחרון 24, יש עדיין ויכוח ניכר בתוך קהילת המחקר החינוכית על היתרונות היחסי של ידיים על ניסויי מעבדה לעומת מעבדות מדומה ומודלים ממוחשבים. מצד אחד, חלק מרגיש ש" התנסות היא בלב של למידה "25, וזהירות שעלות חיסכון טיעונים ניתן תדלוק החלפת ידיים על מעבדה פעילויות על ידי סימולציות מבוססות מחשב, לרעת הבנה של תלמידי 26. מצד השני, חלק מהחוקרים בחינוך מדעי / הנדסה טוענים כי סימולציות הן יעילות לפחות כמו מסורתי ידיים על מעבדות 27, או לדון את היתרונות של מחשב סימולציה בטיפוח "למידת גילוי" ממוקד-תלמיד 28. בעוד קונסנסוס לא מחדשכאב, חוקרים רבים הגיעו למסקנה כי, באופן אידיאלי, הדמיות מחשב צריכים להשלים, ולא להחליף, ידיים על ניסויי מעבדה 29,30. היו גם יוזמות בתוך מדע וחינוך הנדסה לבו זמנית כמה ניסויים פיזיים ועולם האמיתי חישה עם סימולציות מחשב של התופעות; ראה, למשל, "דוגמנות דו מוקדי" 31.
סטודנטים יכולים להרוויח יותר ידע רעיוני והבנה טובה יותר של תהליך המחקר המדעי על ידי אינטראקציה עם שתי מערכת פיזית, וסימולציה מבוססת מחשב של מערכת ש. הליך זה כרוך שיש תלמידים לבצע ניסוי תחבורה מומסת המדגים זרימת מטבע hyporheic הכבידה ומושרה bedform, ולהתאים את תוכנית ההתקנה שלהם ניסיונית ותוצאות עם הדמיית מחשב של אותה התופעה. השוואה זו מאפשרת תוצאות חשובות ללמוד-תלמיד, ודיון עמוק יותר של tהוא שיטה מדעית, ויחסי גומלין בין מודל / תיאורית הבנייה ותיקוף אמפירי באמצעות איסוף נתונים. לאחר ביצוע השוואה זו, תלמידים יכולים גם לנצל את היתרונות של סימולציה ממוחשבת כדי לחקור במהירות מספר רב של תרחישים חלופיים על ידי שינוי פרמטרים מודל.