Asymptotic preserving scheme for the shallow water equations with non-flat bottom topography and Manning friction term

तकनीकी सारांश: गैर-समतल तल स्थलाकृति (Non-Flat Bottom Topography) और मैनिंग घर्षण (Manning Friction) के साथ उथले जल समीकरणों (Shallow Water Equations) के लिए एसिम्प्टोटिक प्रिजर्विंग स्कीम (Asymptotic Preserving Scheme)

1. समस्या विवरण

यह शोध पत्र उथले जल समीकरणों (SWEs) के संख्यात्मक सिमुलेशन को संबोधित करता है, जो दो जटिल भौतिक घटनाओं के अधीन हैं: गैर-समतल तल स्थलाकृति (non-flat bottom topography) और मैनिंग घर्षण (Manning friction)।

यह प्रणाली एक पैरामीटर $\epsilon$ द्वारा नियंत्रित बहु-पैमाने (multiscale) की प्रकृति से युक्त है। जब $\epsilon$ बड़ा ($O(1)$) होता है, तो प्रणाली संवहन-प्रधान (convection-dominated/hyperbolic) शासन में होती है, जहाँ शॉक (shocks) और विच्छेदन (discontinuities) प्रचलित होते हैं। जैसे-जैसे $\epsilon \to 0$ होता है, प्रणाली एक घर्षण-प्रधान (friction-dominated/diffusive) शासन में परिवर्तित हो जाती है, जो एक गैर-रैखिक, अपभ्रंश (degenerate) परवलयिक समीकरण द्वारा शासित होती है जिसे $p$-लैपलेसियन ($p$-Laplacian) कहा जाता है।

मानक संख्यात्मक विधियाँ एक दुविधा का सामना करती हैं:

• स्पष्ट विधियाँ (Explicit methods): ये सरल हैं लेकिन $\epsilon \to 0$ होने पर गंभीर स्थिरता बाधाओं (stiffness) से ग्रस्त होती हैं, जिसके लिए अत्यंत छोटे समय चरणों (time steps) की आवश्यकता होती है।
• निहित विधियाँ (Implicit methods): ये स्थिर हैं लेकिन गणनात्मक रूप से महंगी हैं।
• पूर्ववर्ती AP (Asymptically Preserving) विधियाँ: (विशेष रूप से लेखकों के पिछले कार्य में) स्थिरता सुनिश्चित करने के लिए "पेनलाइजेशन तकनीक" (भारित विसरित पदों को जोड़ना) का उपयोग किया गया था, लेकिन इससे गणनात्मक लागत बढ़ गई और मध्यवर्ती शासन में अभिसरण (convergence) धीमा हो गया।

2. कार्यप्रणाली (Methodology)

लेखक एक नए वर्ग के उच्च-क्रम सेमी-इम्प्लिसिट इम-एक्स रनगे-कुट्टा (SI-IMEX-RK) स्कीम्स का प्रस्ताव करते हैं। मुख्य नवाचार पेनलाइजेशन पद को हटाना है जबकि AP गुण को बनाए रखना है।

कार्यप्रणाली के प्रमुख घटक शामिल हैं:

• समय विविक्तकरण (Time Discretization - SI-IMEX-RK): योगात्मक पेनलाइजेशन के बजाय, लेखक एक सेमी-इम्प्लिसिट ढांचे का उपयोग करते हैं। वे कठिन (stiff) पदों (घर्षण और दबाव प्रवणता) को निहित (implicitly) रूप से और गैर-कठिन (non-stiff) पदों (संवहन) को स्पष्ट (explicitly) रूप से उपचारित करते हैं। उच्च-क्रम सटीकता सुनिश्चित करने के लिए, वे स्टिफली एक्यूरेट (SA) DIRK (डायगोनली इमप्लिसिट रनगे-कुट्टा) स्कीम्स का उपयोग करते हैं।
• स्थानिक विविक्तकरण (Spatial Discretization):
  • प्रथम-क्रम के व्युत्पन्न (First-order derivatives): शॉक को सटीक रूप से पकड़ने और दोलनों (oscillations) को रोकने के लिए उच्च-क्रम WENO (Weighted Essentially Non-Oscillatory) पुनर्निर्माण का उपयोग करके विविक्त किया गया है।
  • द्वितीय-क्रम के व्युत्पन्न (Second-order derivatives): विसरित सीमा (diffusive limit) को संभालने के लिए एक सेंट्रल डिफरेंस स्कीम का उपयोग किया गया है।
• वेल-बैलेंस्ड गुण (Well-Balanced Property): "लेक एट रेस्ट" (lake at rest) स्थिर अवस्था ($h+B = \text{const}, m=0$) को बनाए रखने के लिए, लेखक एक विशिष्ट फ्लक्स-स्प्लिटिंग तकनीक और एक सुसंगत तल स्थलाकृति प्रवणता (bottom topography gradient) उपचार लागू करते हैं।
• मोमेंटम नॉर्म का निहित उपचार (Implicit Treatment of Momentum Norm): एक महत्वपूर्ण तकनीकी निष्कर्ष यह है कि मैनिंग घर्षण पद में संवेग (momentum) के $L_2$ नॉर्म ($|m|$) को AP गुण बनाए रखने के लिए निहित (implicitly) रूप से उपचारित किया जाना चाहिए।

3. मुख्य योगदान

1. पेनलाइजेशन का उन्मूलन: एक ऐसी स्कीम विकसित की गई है जो कृत्रिम भारित विसरित पदों की आवश्यकता के बिना AP गुण प्राप्त करती है, जिससे गणनात्मक ओवरहेड कम हो जाता है।
2. उच्च-क्रम सटीकता: यह सिद्ध किया गया है कि यह स्कीम न केवल AP (सीमा समीकरण के अनुरूप) है, बल्कि एसिम्प्टोटिकली एक्यूरेट (AA) भी है, जिसका अर्थ है कि यह $\epsilon \to 0$ होने पर भी अपनी उच्च-क्रम कालिक सटीकता (temporal accuracy) बनाए रखती है।
3. मध्यवर्ती शासन में दक्षता: यह प्रदर्शित किया गया है कि SI-IMEX दृष्टिकोण पूर्ववर्ती योगात्मक IMEX विधियों की तुलना में काफी अधिक कुशल है, विशेष रूप से संवहन और विसरण के बीच संक्रमण के दौरान।
4. सैद्धांतिक कठोरता: चेपमैन-एनस्कोग (Chapman-Enskog) विस्तार और गणितीय आगमन (mathematical induction) का उपयोग करके AP और AA गुणों के लिए औपचारिक प्रमाण प्रदान किए गए।

4. परिणाम और सत्यापन

लेखकों ने कई संख्यात्मक प्रयोगों के माध्यम से अपनी स्कीम को सत्यापित किया:

• 1D सटीकता परीक्षण: विभिन्न $\epsilon$ स्केल्स (1 से $10^{-6}$ तक) में रैखिक और गैर-रैखिक मैनिंग घर्षण दोनों के लिए तृतीय-क्रम सटीकता की पुष्टि की।
• AP गुण सत्यापन: दिखाया गया कि जैसे-जैसे $\epsilon$ घटता है, संख्यात्मक समाधान $p$-लैपलेसियन सीमा समीकरण के समाधान की ओर अभिसरित होता है।
• योगात्मक IMEX के साथ तुलना: प्रयोगों ने सिद्ध किया कि SI-IMEX स्कीम को अभिसरण तक पहुँचने के लिए काफी कम पुनरावृत्तियों (iterations) की आवश्यकता होती है, विशेष रूप से मध्यवर्ती और विसरित शासन में।
• 2D परीक्षण: जटिल 2D सुचारू और विच्छिन्न प्रारंभिक स्थितियों को संभालने की स्कीम की क्षमता, साथ ही तल स्थलाकृति के चारों ओर तरंग प्रसार (तरंग प्रसार) को पकड़ने की इसकी क्षमता (वेल-बैलेंस्ड गुण) को प्रदर्शित किया।
• निहित $|m|$ की आवश्यकता: परीक्षणों ने दिखाया कि $|m|$ को स्पष्ट (explicitly) रूप से उपचारित करने (SI-S2 स्कीम) से स्कीम विसरित सीमा में विफल हो जाती है, जिससे पुष्टि होती है कि स्थिरता और AP संरक्षण के लिए निहित उपचार अनिवार्य है।

5. महत्व

यह कार्य भू-भौतिक तरल गतिकी (geophysical fluid dynamics) के लिए एक मजबूत, उच्च-क्रम संख्यात्मक ढांचा प्रदान करता है। पेनलाइजेशन की "ब्रूट फोर्स" पद्धति के बिना इस कठोरता (stiffness) की समस्या को हल करके, यह स्कीम महासागरीय और नदी प्रवाह के कुशल, बड़े पैमाने के सिमुलेशन की अनुमति देती है जहाँ तल स्थलाकृति और घर्षण प्रभावी होते हैं। यह हाइपरबोलिक और परवलयिक शासनों के बीच के अंतर को पाटता है, जिससे यह बहु-स्तरीय पर्यावरणीय मॉडलिंग के लिए अत्यधिक उपयुक्त बन जाता है।