A consistent-splitting generalized scalar auxiliary variable scheme for the perturbed Boussinesq system

तकनीकी सारांश: विक्षुब्ध बुसिनस्क (Boussinesq) प्रणाली के लिए एक सुसंगत-विभाजन (Consistent-Splitting) सामान्यीकृत स्केलर ऑक्जिलरी वेरिएबल (GSAV) योजना

समस्या विवरण
यह कार्य द्वि-आयामी विक्षुब्त बुसिनस्क प्रणाली के संख्यात्मक सिमुलेशन को संबोधित करता है, जो स्थिर रूप से स्तरित (stably stratified) वातावरण (जैसे वायुमंडलीय और महासागरीय गतिकी) में उछाल-संचालित (buoyancy-driven) प्रवाह का मॉडल तैयार करती है। यह प्रणाली मानक बुसिनस्क समीकरणों से एक स्थिर, रैखिक रूप से स्तरित हाइड्रोस्टेटिक साम्यावस्था को घटाकर व्युत्पन्न की गई है। परिणामी विक्षुब्त प्रणाली (समीकरण 6) में एक विशिष्ट युग्मन (coupling) है जहाँ उछाल (buoyancy) मोमेंटम समीकरण में ऊर्ध्वाधर गति को संचालित करता है, जबकि ऊर्ध्वाधर गति तापमान समीकरण में रैखिक पृष्ठभूमि तापमान प्रवणता (linear background temperature gradient) का अधोवहन (advection) करती है। यह युग्मन आंतरिक गुरुत्वाकर्षण तरंगों को उत्पन्न करता है और विक्षोभों को हाइड्रोस्टेटिक साम्यावस्था की ओर शिथिल करने (relaxation) को प्रेरित करता है। प्राथमिक गणितीय चुनौती असंपीड्यता बाधा (incompressibility constraint) और वेग-दबाव प्रणाली की सैडल-पॉइंट संरचना है, जो गैररेखीय संवहन (nonlinear convection) और रैखिक स्तरीकरण पदों की उपस्थिति में उच्च-क्रम सटीकता और स्थिरता की आवश्यकता से और भी जटिल हो जाती है।

कार्यप्रणाली (Methodology)
लेखक एक द्वितीय-क्रम सुसंगत-विभाजन योजना प्रस्तावित करते हैं जो सामान्यीकृत स्केलर ऑक्जिलरी वेरिएबल (GSAV) दृष्टिकोण पर आधारित है। कार्यप्रणाली कई प्रमुख घटकों को एकीकृत करती है:

1. GSAV ढांचा: एक स्केलर ऑक्जिलरी वेरिएबल $r(t)$ को प्रणाली के निरंतर ऊर्जा नियम (continuous energy law) की नकल करने के लिए पेश किया गया है। यह एक संशोधित विविक्त ऊर्जा (discrete energy) के निर्माण की अनुमति देता है जो बिना शर्त स्थिर (unconditionally stable) रहती है।
2. सुसंगत-विभाजन समय विविक्तकरण (Consistent-Splitting Time Discretization): यह योजना ह्वांग और शेन [17] द्वारा नेवियर-स्टोक्स समीकरणों के लिए विकसित 'टाइम-शिफ्टेड सेकंड-ऑर्डर BDF2' ढांचे का विस्तार करती है।
   • स्पष्ट उपचार (Explicit Treatment): गैररेखीय संवहन ($u \cdot \nabla u$), अधोवहन ($u \cdot \nabla \theta$), और रैखिक युग्मन (उछाल $\theta e_2$ और स्तरीकरण $\alpha v$) को शिफ्ट किए गए समय स्तर $t_{n+k}$ पर द्वितीय-क्रम एक्सट्रपलेशन (extrapolation) का उपयोग करके स्पष्ट रूप से उपचारित किया जाता है।
   • विघटन (Decoupling): यह स्पष्ट उपचार प्रत्येक टाइम स्टेप को स्थिर गुणांकों वाले कुछ छोटे विघटित रैखिक प्रणालियों को हल करने तक सीमित कर देता है, जिससे गैररेखीय सॉल्वर या मोनोलिथिक युग्मित प्रणालियों की आवश्यकता समाप्त हो जाती है।
3. दबाव सुधार (Pressure Correction): एक दबाव पॉइसन समीकरण को एक सुसंगत विभाजन दृष्टिकोण का उपयोग करके हल किया जाता है जो इष्टतम द्वितीय-क्रम सटीकता बनाए रखता है, जिससे मानक प्रोजेक्शन विधियों में पाई जाने वाली सीमा स्थिति विसंगतियों (boundary condition inconsistencies) से बचा जा सके।
4. स्थिरता और त्रुटि विश्लेषण: लेखक एक बिना शर्त कमजोर स्थिरता प्रमेय (unconditional weak stability theorem) सिद्ध करते हैं और वेग, दबाव, और तापमान के लिए इष्टतम द्वितीय-क्रम त्रुटि अनुमान प्राप्त करते हैं। विश्लेषण कठोरता से किनेमैटिक श्यानता ($\nu$) और थर्मल डिफ्यूज़िविटी ($\gamma$) पर त्रुटि स्थिरांकों की निर्भरता को ट्रैक करता है।

प्रमुख योगदान
बुसिनस्क प्रणालियों पर मौजूदा साहित्य की तुलना में यह शोध पत्र तीन विशिष्ट योगदान देता है:

1. विक्षुब्ध प्रणाली विश्लेषण: पिछले कार्यों के विपरीत, जो मानक बुसिनस्क प्रणाली का विश्लेषण करते हैं, यह शोध पत्र विशेष रूप से तापमान समीकरण में रैखिक स्तरीकरण युग्मन पद $\alpha v$ वाले विक्षुब्ध प्रणाली (6) को संबोधित करता है।
2. ढांचा एकीकरण: जबकि वाग्नर एट अल. [36] ने मानक बुसिनस्क प्रणाली पर GSAV दृष्टिकोण लागू किया, यह कार्य ह्वांग और शेन [17] के विशिष्ट टाइम-शिफ्टेड सेकंड-ऑर्डर BDF सुसंगत-विभाजन ढांचे को अपनाता है, लेकिन मानक GSAV अपडेट तंत्र को बनाए रखता है, जिससे सभी युग्मन पदों के सुसंगत एक्सट्रपलेशन के माध्यम से वेग और तापमान दोनों समीकरणों में पूर्ण द्वितीय-क्रम सटीकता सुनिश्चित होती है।
3. स्पष्ट त्रुटि स्थिरांक निर्भरता: एक महत्वपूर्ण योगदान $\nu^{-1}$ और $\gamma^{-1}$ पर त्रुटि स्थिरांक की निर्भरता का स्पष्ट अनुरेखण (tracing) है। विश्लेषण प्रकट करता है कि त्रुटि स्थिरांक इन मापदंडों पर श्यानता और डिफ्यूज़िविटी की ऋणात्मक घातों को समाहित करता है, जिससे इन मापदंडों पर चतुर्गुणी-नेस्टेड घातांकीय (quadruply-nested exponential) निर्भरता उत्पन्न होती है।

परिणाम

• सैद्धांतिक स्थिरता: यह सिद्ध किया गया है कि योजना बिना शर्त रूप से कमजोर रूप से स्थिर (unconditionally weakly stable) है। ऑक्जिलरी वेरिएबल और स्केलिंग फैक्टर के अनुक्रम गैर-ऋणात्मक रहते हैं, और संख्यात्मक समाधान यूनिफॉर्म बाउंड्स का पालन करते हैं, बशर्ते कि स्थिरीकरण स्थिरांक $C$, फोर्सिंग टर्म्स द्वारा निर्धारित थ्रेशोल्ड से अधिक हो।
• अभिसरण (Convergence): उपयुक्त नॉर्म्स में वेग, दबाव और तापमान के लिए इष्टतम द्वितीय-क्रम अभिसरण दरें स्थापित की गई हैं।
• संख्यात्मक प्रयोग:
  • अभिसरण परीक्षण: एक सीमित डोमेन पर एक मैन्युफैक्चरड सॉल्यूशन (manufactured solution) दोनों $L^2$ और $H^1$ नॉर्म्स में सैद्धांतिक द्वितीय-क्रम अभिसरण की पुष्टि करता है।
  • स्तरित प्रवाह सिमुलेशन: प्रारंभिक गॉसियन तापमान विक्षोभ के साथ एक स्तरित प्रवाह का दीर्घकालिक सिमुलेशन ($T=10^4$) अपेक्षित भौतिक गतिकी को पुनरुत्पादित करता है:
    • एक उछाल क्षणिक शासन (buoyancy transient regime) जहाँ संभावित ऊर्जा गतिज ऊर्जा में परिवर्तित होती है।
    • एक आंतरिक-तरंग दोलन शासन जिसका आवर्त काल सैद्धांतिक भविष्यवाणियों ($T \approx \pi\sqrt{2/\alpha}$) के अनुरूप है।
    • एक घातांकीय क्षय शासन जहाँ विक्षोभ हाइड्रोस्टेटिक साम्यावस्था की ओर शिथिल होते हैं, जिसमें क्षय दर लेम्मा 2 में प्राप्त विश्लेषणात्मक सीमाओं के अनुरूप होती है।

महत्व और सीमाएँ
यह शोध पत्र दावा करता है कि प्रस्तावित योजना स्तरित प्रवाहों के सिमुलेशन के लिए एक गणनात्मक रूप से कुशल विधि प्रदान करती है, जो प्रत्येक टाइम स्टेप को विघटित रैखिक समाधानों तक कम कर देती है जबकि द्वितीय-क्रम सटीकता बनाए रखती है। हालाँकि, लेखक शून्य श्यानता या डिफ्यूज़िविटी की सीमा में योजना की मजबूती (robustness) के संबंध में विनम्र हैं।

विश्लेषण दर्शाता है कि त्रुटि स्थिरांक $\nu^{-1}$ और $\gamma^{-1}$ के चतुर्गुणी-नेस्टेड घातांकीय (quadruply-nested exponential) के रूप में बढ़ता है। फलस्वरूप, योजना $\nu, \gamma \to 0$ होने पर मजबूत (robust) नहीं है। संख्यात्मक प्रयोग पुष्टि करते हैं कि कम श्यानता और थर्मल डिफ्यूज़िविटी के लिए, योजना को स्थिर रहने के लिए तेजी से सूक्ष्म स्थानिक और टेम्पोरल रेजोल्यूशन की आवश्यकता होती है, और यह मोटे रेजोल्यूशन पर अस्थिर (blow up) हो सकती है। लेखक निष्कर्ष निकालते हैं कि जबकि यह योजना पूरी तरह से विसरित (dissipative) शासन के लिए प्रभावी है, $\nu$ और $\gamma$ में यूनिफॉर्म एरर बाउंड प्राप्त करने के लिए संवहन पदों के लिए पूरी तरह से इम्प्लिसिट या स्टेबलाइज्ड उपचार की आवश्यकता होगी। यह सीमा स्पष्ट-विभाजन (explicit consistent-splitting) दृष्टिकोण की गणनात्मक दक्षता और कम-डिफ्यूजन शासन में इसकी स्थिरता के बीच के समझौते (trade-off) को रेखांकित करती है।