שיטת זיהוי פני שטח חדשה לחישוב MPS תלת־ממדי

מדוע חשוב לעקוב אחרי פני מים נעים

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

דרך שונה להבחין בין מים לאוויר

בשיטת Moving Particle Semi-Implicit (MPS) ממודלים המים כלא־מספר גדול של חלקיקים זעירים שנעים ומתקשרים זה עם זה. גישות מסורתיות קובעות האם חלקיק נמצא על פני השטח על ידי ספירה פשוטה של כמה שכנים יש לו. בדרך כלל פחות שכנים מרמז שהוא על או קרוב לפני השטח. עם זאת, כלל אצבע זה עלול להיכשל כאשר החלקיקים מרווחים באופן לא אחיד, כפי שקורה לעתים בזרמים תזזיתיים או ליד פינות חדות. החוקרים במקום זאת מתמקדים באופן שבו השכנים מסודרים במרחב, ומשתמשים בכמות הנקראת סטייה יחסית של מיקומים (relative position divergence), שמודדת עד כמה השכנים "מפוזרים" סביב חלקיק נתון. אם השכנים נוטים להימצא בעיקר בצד אחד, ערך הסטייה משתנה באופן שמצביע על סיכוי גבוה שחלקיק זה נמצא על פני השטח.

הוספת כיוון למבחן הפנים

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

מבחנים על צורות מורכבות ומים שקטים

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

לכידת גלים, פגיעות ונענוע נוזלים

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

מה זה אומר עבור סימולציות בעולם האמיתי

למספרים מחוץ לתחום של דינמיקת נוזלים חישובית, המסר המעשי הוא שעבודה זו משפרת את אופן שבו מחשבים מחליטים היכן מסתיימים המים והיכן מתחיל האוויר במודלים מבוססי חלקיקים. על ידי שילוב מדד מעודן לצורת סידור השכנים סביב כל חלקיק עם מבחן כיוון "קונוס ריק", שיטת RPD+NV החדשה מקטינה באופן ניכר סיווגים מוטעים של חלקיקי משטח. התוצאה היא חיזויי לחץ נקיים ויציבים יותר במגוון זרימות — ממים עומדים ועד פגיעות אלימות ונענוע תלת־ממדי מורכב. מכיוון שהשיטה מוסיפה עלות חישובית מועטה בלבד, היא מציעה שדרוג מעשי למהנדסים המדמים מבני חוף, תנועת ספינות, מכלי אחסון ומערכות אחרות שבהן לכידת פני שטח נעים של מים בדיוק חיוני.

ציטוט: Geng, C., Wang, Wh., Heng, My. et al. A new free surface identification method for 3D MPS method. Sci Rep 16, 13829 (2026). https://doi.org/10.1038/s41598-026-44218-9

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