חשיבות התנאים הסגורים פנימה נאוייר-טסטוק פתרונות נומרניים

התפקיד האכזרי של תנאי Boundary ב Navier-Stokes Numerical Solutions

משוואות Navier-Stokes מהוות אבן הפינה של דינמיקת נוזל חישובית (CFD), השולטת בהתנהגות של זרמי נוזל מולקווס. בעוד המשוואות עצמן ללכוד את הפיזיקה של שימור התנופה, הן אינן שלמות ללא תיאור מדויק של מה שקורה בגבולות של פתרונות סימולציה חישובית. Boundary תנאים לא רק פרטים חישוביים; הן לאכוף את האינטראקציה פיזית בין נוזל וסביבתה.

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

יסודות מתמטיים של תנאים

לפני צלילה לסוגים ספציפיים של גבולות, חשוב להבין כיצד התנאים הגבול מתאימים למסגרת המתמטית של משוואות Navier-Stokes. משוואות השלטון בצורה לא יעילה הם:

(ב) ויקרא י"ד:

(ב) כאשר הוא ⁇ :0 (uentiFLT:1 ⁇ הוא צפיפות, p הוא לחץ, ⁇ הוא סטיות kinematic מולקולות, וFLT:2freaFLT 3: 3 מייצג כוחות גוף.משוואות אלה דורשות תנאים גבול על כל משטחים מקיפים את התחום, כמו גם תנאים ראשוניים לבעיות טרנספורמטיביות.

ג'נריק סיווג: Dirichlet, Neumann ורובין

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

(FLT:0 (חומר) תנאי גבול: ההרחבה: 1:1 אלה ישירות לציין את הערך של המשתנה תלוי (הרחבה או הלחץ) בגבול.לדוגמה, הגדרת המהירות וקטור לאפס על הקיר המתוצב (ללא מצב) הוא מצב דיריץ' על מהירות.
(FLT:0) תנאי גבול טבעיים: ⁇ 1 (התוצאה הסופית) של המשתנה התלויה (טבעי) עבור משוואות Navier-Stokes, דוגמא נפוצה היא "המצב החופשי של רצף" או "זרימה" על גבול יציאה, שבו נגזרת רגילה של מהירות נקבעת לאפס ( ⁇ u /n=n=n) ו הלחץ הוא לעתים קרובות מוגדר למתח קבוע לשימוש בקירות זעירים.
(FLT:0) רובין (מעורבים) תנאי גבול: ההרחבה 1 (המונחים של Dirichlet ו- Neumann), בדרך כלל כשילוב ליניארי של המשתנה וגזרתו הרגילה.תנאים מופיעים בבעיות העברת חום או כאשר דוגמנות בזרימות מיקרו-scale (למשל, מצב חלקי: u +α ⁇ u /n= g).

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

על מצב הלחץ

קצבה מיוחדת מתעוררת ללחץ.במצבה המעמיק של נאוייר-Stokes, הלחץ אינו משתנה תרמודינמי אלא מכפיל Lagrange המאוחסן שימור המוני.מצבו אינו נגזר ישירות מהפיזיקה אלא מעקביות עם תנאי הגבול המהיר.בדרך כלל, בגבולות זרימה, חד-פעמי אחד (Dirichlet) ולחצים נוספים (Neumann) או לחץ קבוע על ידי משוואה.

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

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

No-Slip and Slip Walls

ה-FLT:0 slipמצב FLT:1 (u= v=0 על הקיר) הוא ברירת המחדל עבור רוב זרמי ה-cous בקנה מידה מאקרוסקופי.זה נובע מן ההנחה כי חלקיקים נוזליים ליד משטח מוצק יש אפס מהירות שיזוף יחסית עקב דבקות מולקולוס.

עם זאת, במשטרים מסוימים – כמו מיקרו-השפעה, זורם מעל פני השטח המרומזים, או זרמים עם שר רשלני – ה-FLT:0slip ModeFLT:1 הופך מתאים.

⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇

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

המונחים: Inlet andOut Conditions

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

(FLT:0)Uniform או מהירות פרופיל אינלט: ראט' 1 (עבור זרימת צינורות laminar, פרופיל פרבולי מפותח במלואו נקבע. עבור אווירודינמיקה חיצונית, זרם חופשי אחיד עם אינטנסיביות סוערת קטנה משמש.
(FLT:0) סיקור אינפורמטיבי: 1FLT לבעיות שבהן הזרם מונע על ידי הבדל לחץ, ניתן לציין לחץ קיפאון במבוא.זה נפוץ בזרימים דחוסים או כאשר המהירות של החדירה אינה ידועה.
(FLT:0Outflow (Neumann): FLT:1 ב שקעים, הגדרת אפס ⁇ רגיל עבור כל המשתנים (כולל מהירות וכמויות זעזוע) הוא סטנדרטי.זה מאפשר זרימה ליציאה ללא השפעה על הזרם, בתנאי שהיציאה ממוקמת די רחוק מאזורי עניין.
(FLT:0) זרימה חיצונית:FLT:1 עבור זרמים לא יציבים עם מערבים היוצאים מהתחום, מצב גבול משותף ( ⁇ / ⁇ t + U ⁇ / ⁇ n = 0) יכול להפחית את ההשתקפות ולשמור דיוק.

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

תנאי תקופתיים וסימטים

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

⁇ (x + L) = ⁇ (x) + ⁇ pphLT:0dledleduaFLT:1 עבור לחץ, ו ⁇ = אותו הדבר עבור מהירות.

יישום תנאים תקופתיים דורש התאמה זהירה של צמתים ברשת וטיפול בקפיצות לחץ (למשל, לחץ ⁇ ) (FLT:0Symmetry תנאים של LT:1 ), לכפות אפס מהירות נורמלית אפס ⁇ נורמלית של קשקשים אחרים ( ⁇ / ⁇ n= 0) הם שימושיים עבור תוכנית או זרמי אלקיסטמי אבל לא צריך לשמש באופן חד-סימטרי.

תנאי פתוח ל- Unbounded Flows

(למשל, זרימה מעל כנף), יש להזיז את התחום.גבולות מלאכותיים – שטח רחוק, זרם חופשי או קרינה – נועדו לתת הפרעות לעבור ללא השתקפות.הגישה הפופולרית ביותר היא ה-FLT:0characteristic ConditionFLT:1 המבוססת על Riemann invariants (עבור זרימה בלתי דחוסה) או LT2 דורש לעתים קרובות ספוגה קרינת חום (Ricial) או ספוגה).

השפעה על יציבות נומרית וכלכלה

הבחירה של תנאי גבול משפיעה ישירות על המיזוג והמיזוג של מערכת הדיסקרטית.עבור משוואות Navier-Stokes המדכאות, הפתרון המזוגי או המבודד (למשל, SIMPLE, PISO) תלוי בתנאי גבול כדי לייצר מאטריקס לחץ לא סמוי.מצב נוימן טהור על לחץ בכל מקום (כמו בתחום סגור) מוביל למערכת טיפול קבועה; יש צורך בלחצים קבועים.

שיקולים

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

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

אמינות ויציבות גריידי

תנאים קשים מציגים שגיאות באמצעות אי-התאמה והתערבות.לדוגמה, כאשר קיר ללא-Slip משתמש ברשת קוהרזה, מהירות ⁇ - ומכאן הלחץ הקיר-הקיר - ניתן להישמר בחומרה.

(FLT:0) טיפול בקיר: FLT:1 שימוש בפונקציות הקיר (למשל, חוק-העץ) אשר כופים מצב חלקלק בהתבסס על המרחק לקיר, ובכך מרגיע את הדרישה של רשת הקיר.
(ב) ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
שיטות תאים:0 (Ghost Celleur: FLT:1) ,ההחלים את הפתרון לתאי רוח באמצעות ערכי גבול מוגדרים מבטיחים דיוק סדר גבוה מבלי לצמצם את הסדר הפנימי.

שיעור ההתכנסות של פותרים רציונאליים גם מתפתל בתנאי גבול.לדוגמה, יחידות פינה (שם קיר ללא-Slip עומד בשקע) יוצר נקודה ייחודית מתמטית המפרקת את ההתכנסות של שיטות מרובותgrid.

טכניקות מתקדמות ומחקרים נוכחיים

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

שיטות סודיות

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

קוד פתוח ופתרון מסחרי

תוכנת CFD מודרנית כגון FLT:0 (OpenFOAMIRLT:1 ו- (FLT:2ANSYS FluentFLT 3: 3) מציעים ספריות נרחבות של תנאי גבול.עם זאת, משתמשים חייבים להבין את ההבדלים ביישום (למשל, זרימה קבועה לעומת תיקון קבוע של אפשרויות דחיסה) ואת התכניות המספריות שמאחוריהם, לדוגמה, ב-AMCREC, בתנאי לחץ על ידי תיקון מהיר (Fure) אשר עשוי להבטיח את ה-reme) באמצעות תצורה של הגבלת הגבולות ה-repericows-Fure, כאשר הם עלולים, כאשר הם עלולים, כאשר הם עלולים להוביל את ה-Fure, כאשר הם עלולים להחלפה, כאשר הם עלולים, כאשר הם עלולים ל-remecremecremed) בתנאי גבול , כאשר הם עלולים ל-Fure, כאשר הם עלולים להתרחש באמצעות , כאשר הם עלולים להתרחש באמצעות , כאשר הם עלולים ל-Fure, כאשר הם עלולים להתרחש, כאשר הם עלולים ל-Fure, בתנאי גבול מוגבל, כאשר הם עלולים להתרחש, כאשר הם עלולים להבטיח את ה-remed, כאשר הם עלולים להתרחש ב-remecremecreme

נתונים-Driven ו- Machine Learning Boundary Conditions

מחקר מתפתח משתמש ברשתות עצביות כדי להפר תנאים של מדידות חלקיות או להחליף פונקציות קיר אמפיריות.לדוגמה, FLT:0physics-מודעות רשתות עצביות (PINNs)OVAFLT:1 יכול לפתור את משוואות הסימולציות של נא-Stoke עם נתונים מלוחים, ביעילות ללמוד את תנאי הגבול המתאימים מן הנתונים, בעוד שעדיין בשלבים המוקדמים, גישות כאלה לפתרונות אמיתיים לזרימה של פתרונות מסחר אלקטרונים (Freplelimate) ו-of-Frepled) נמצאים ב-Freorderd (במסגרת משוואות בקרה דיגיטלית, לעומת שימוש ב-Freplelimates:2, כלומר, כלומר, כלומר, משוואות נוספות, כלומר, כלומר, משוואות של שימוש בגרסאות של שימוש ב-Frepleplelimates:

מלכודות נפוצות ועיסוקים טובים

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

גודלו של LT:0 [Domain Size:FLT:1 Place inlet and outlet Limits far than the Zone of Interest (למשל, לפחות 10 אורכו של ספקטרום ו-20 אורך במורד הזרם עבור גופים bluff).
בדיקה אחרונה ב-17 במאי 2010. ^ FLT:0.2017, App:0.comsistency Check: FLT:1hil, ודא כי סך זרימת המסה נכנסת לשווי של זרימה המונית הכוללת, המשאירה את הזרמים הבלתי ניתנים לדיכוי יציב; התאימו את רמת הלחץ או השתמשו בתיקון לחץ ההתייחסות.
[ה] עצמאות:0 [Grid Independence:] דחיית [ה] לקירות ולעוסקים שבהם ⁇ s הם גבוהים. בדוק כי יישום תנאי הגבול אינו יוצר עששת (לא פיזיקליות) בפתרון.
(FLT:0) רזולוציה של סימולציה חד-ממדית (ד) או סימולציה מספרית ישירה (DNS), תנאי גבול ב Inlets חייבים להיות סוערים.
(ב) טיפול בזרימה:0) בזרימה: 1:1 בכלים, אם זרימה עשויה להפוך לכיוון הפוך, להשתמש בלחיצת לחץ או מצב קולקטיבי כדי להימנע מהבדלים.
(ב) ,0) ,Validation: FLT:1 לעומת תוצאות עם נתונים ניסיוניים או פתרונות אנליטיים (למשל, FLT:2Poiseuille FlowFLT 3 עבור זרימה פנימית).

תגית: Lid-Driven Cavity Flow

אחת הבעיות האיקוניות ביותר של נווארה-Stokes היא המשטח המונע על ידי הלימוזינה. בבעיה זו, חלל מרובע יש שלושה קירות נייחים (ללא-slip, u=v=0) ומכסה נעה למעלה (u= U, v= 0) האינטראקציה ליקולורית יוצרת פינות גבול משותף בפינות העליונות שבהן הקיר והמעבר לפתרונות המשתנים לעיתים קרובות.

מסקנה

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

לקריאה נוספת על טכניקות מתקדמות של תנאי גבול, ייעוץ:0Computational Fluid Dynamics: Principles and ApplicationsFLT:1 על ידי J. Blazek או הטקסט הקלאסי של טקסט FLT:2Computational Methods for Fluidamide DynamicsFLT 3: by Ferziger and Peric יישום פרטים ניתן למצוא בתיעוד של OpenFOAM ו- AnSYSent, אשר מציעים דוגמאות נרחבות של שימוש בעולם.