A family of second order, linear, unconditionally stable methods for the Cahn-Hilliard-Navier-Stokes equations

तकनीकी सारांश: Cahn-Hilliard-Navier-Stokes समीकरणों के लिए द्वितीय-क्रम, रैखिक, बिना शर्त स्थिर विधियों का एक परिवार

समस्या विवरण
जटिल बहु-चरणीय द्रव गतिकी (multiphase fluid dynamics) का सटीक संख्यात्मक अनुकरण, विशेष रूप से जो Cahn-Hilliard-Navier-Stokes (CHNS) समीकरणों द्वारा मॉडल किए जाते हैं, महत्वपूर्ण चुनौतियाँ प्रस्तुत करता है। इसका कारण चतुर्थ-क्रम (fourth-order) Cahn-Hilliard ऑपरेटर की गणितीय कठोरता (mathematical stiffness) और असम्पीड्य (incompressible) Navier-Stokes समीकरणों के साथ इसकी गैररेखीय युग्मन (nonlinear coupling) है। मानक स्पष्ट समय-चरण योजनाएं (explicit time-stepping schemes) स्थिरता बनाए रखने के लिए गंभीर समय-चरण प्रतिबंधों (जैसे, $\Delta t \sim O(\Delta x^4)$) का सामना करती हैं। जबकि पूर्णतः अंतर्निहित (fully implicit) या कॉन्वेक्स-स्प्लिटिंग (convex-splitting) विधियाँ बिना शर्त ऊर्जा स्थिरता प्रदान करती हैं, उन्हें अक्सर प्रत्येक समय चरण में बड़े, युग्मित गैररेखीय प्रणालियों को हल करने की आवश्यकता होती है, जिससे कम्प्यूटेशनल लागत बढ़ जाती है। एक ऐसे एल्गोरिदम की आवश्यकता है जो द्वितीय-क्रम की अस्थायी सटीकता (temporal accuracy), बिना शर्त स्थिरता और रैखिक समाधान क्षमता को जोड़ सके, बिना प्रणाली के भौतिक संरक्षण गुणों से समझौता किए।

कार्यप्रणाली (Methodology)
लेखक मेल-मैचड-डेंसिटी (matched-density) दो-चरणीय प्रवाह के लिए द्वितीय-क्रम, रैकी, बिना शर्त स्थिर 'इम्प्लिसिट-एक्सप्लिसिट' (IMEX) विधियों का एक परिवार प्रस्तावित करते हैं। अर्ध-विविक्त योजना (semi-discrete scheme) तीन प्राथमिक घटकों को एकीकृत करती है:

1. गैररेखीय पदों का एक्सट्रपलेशन (Extrapolation of Nonlinear Terms): गैररेखीय संवहनी (convective) और युग्मन पदों को एक्सट्रपोलेटेड मात्राओं (denote by $H_{n+\theta}$) का उपयोग करके स्पष्ट रूप से (explicitly) उपचारित किया जाता है, जिससे गैररेखीय पुनरावृत्तियों (nonlinear iterations) की आवश्यकता समाप्त हो जाती है।
2. सहायक चर सूत्रीकरण (Auxiliary Variable Formulation): गैररेखीय मुक्त-ऊर्जा पद को एक सहायक चर $q$ (जो डबल-वेल पोटेंशियल $F(\phi)$ से व्युत्पन्न है) के माध्यम से संभाला जाता है, जिससे समस्या को रैखिकता बनाए रखने के लिए पुनर्गठित किया जाता है।
3. अस्थायी-वक्रता नियमितीकरण (Temporal-Curvature Regularization): एक स्थिरीकरण तंत्र जिसे पैरामीटर $\epsilon \geq 0$ द्वारा नियंत्रित किया जाता है, पेश किया गया है। इसमें एक विशिष्ट इंटरपोलेशन ऑपरेटर $J^\epsilon_{n+\theta}$ शामिल है जो समाधान के विविक्त अस्थायी वक्रता (discrete temporal curvature) को सम्मिलित करता है, जो Navier-Stokes नियमितीकरण पर आधारित एक अवधारणा है।

परिणामस्वरूप प्राप्त एल्गोरिदम (Algorithm 2.1) प्रत्येक समय चरण में केवल स्थिर गुणांकों वाले रैखिक प्रणालियों को हल करने की आवश्यकता रखता है। यह योजना $\theta \in (1/2, 1]$ और $\epsilon \geq 0$ द्वारा पैरामीटराइज्ड है।

प्रमुख योगदान और सैद्धांतिक परिणाम

• बिना शर्त स्थिरता (Unconditional Stability): शोध पत्र एक कठोर विविक्त ऊर्जा अनुमान (discrete energy estimate) स्थापित करता है, जो यह सिद्ध करता है कि प्रस्तावित योजना $\theta \in (1/2, 1]$ और $\epsilon \geq 0$ के लिए बिना शर्त दीर्घकालिक रूप से स्थिर है। प्रमाण सममित धनात्मक निश्चित मैट्रिक्स (symmetric positive definite matrices) का उपयोग करके विविक्त नॉर्म्स को परिभाषित करता है और यह प्रदर्शित करता है कि योजना समय-चरण प्रतिबंधों के बिना एक विविक्त ऊर्जा-क्षय नियम (energy-dissipation law) का पालन करती है।
• रैखिक समाधान क्षमता (Linear Solvability): पूर्णतः अंतर्निहित या कॉन्वेक्स-स्प्लिटिंग दृष्टिकोणों के विपरीत, जो गैररेखीय सॉल्वर की आवश्यकता में सकते हैं, यह विधि प्रत्येक समय चरण में रैखिक प्रणालियों को उत्पन्न करती है, जिससे कम्प्यूटेशनल जटिलता काफी कम हो जाती है।
• द्वितीय-क्रम सटीकता (Second-Order Accuracy): सैद्धांतिक विश्लेषण और संख्यात्मक प्रयोग पुष्टि करते हैं कि विधि द्वितीय-क्रम की अस्थायी सटीकता बनाए रखती है।

संख्यात्मक परिणाम
लेखक बेंचमार्क गणनाओं के माध्यम से विधि को मान्य करते हैं:

• अभिसरण विश्लेषण (Convergence Analysis): 'मेथड ऑफ मैन्युफैक्चरड सॉल्यूशंस' (MMS) का उपयोग करते हुए, योजना विभिन्न परिमित तत्व विन्यास ($P2-P1-P2-P2$ और $P3-P2-P3-P3$) में वेग, दबाव और फेज फील्ड चरों के लिए लगभग द्वितीय-क्रम का अस्थायी अभिसरण प्रदर्शित करती है।
• मजबूती और नियमितीकरण (Robustness and Regularization): दो-चरणीय लिड-ड्रिवन कैविटी समस्या में, अध्ययन $\theta$ और $\epsilon$ के बीच परस्पर क्रिया का परीक्षण करता है। यह देखा गया है कि क्रैंक-निकोल्सन सीमा (जैसे, $\theta=0.51$) के करीब $\theta$ मानों के लिए $\epsilon=0$ के साथ, कमजोर रूप से डैम्प्ड (damped) अस्थायी दोलनों के कारण योजना अस्थिर हो सकती है। धनात्मक वक्रता नियमितीकरण ($\epsilon > 0$) का परिचय इन दोलनों को प्रभावी ढंग से कम करता है, जिससे बड़े समय चरणों के लिए भी सीमाबद्धता और मजबूती बहाल होती है।
• भौतिक बेंचमार्क (Physical Benchmarks):
  • स्पिनोडल डिकंपोजिशन (Spinodal Decomposition): विधि दीर्घकालिक सिमुलेशन ($T=100$) पर ऊर्जा क्षय और द्रव्यमान संरक्षण को सटीक रूप से कैप्चर करती है, जो प्रभावी चरण पृथक्करण और मोटा होने (coarsening) को दर्शाती है।
  • ड्रॉपलेट रिलैक्सेशन (Droplet Relaxation): एल्गोरिदम सतह-तनाव-प्रेरित इंटरफेशियल गति को सफलतापूर्वक हल करता है, जो द्रव्यमान को संरक्षित करते हुए एक वर्गाकार ड्रॉपलेट को गोलाकार संतुलन में परिवर्तित करता है।
  • लिड-ड्रिवन कैविटी और रेले-टेलर अस्थिरता (Lid-Driven Cavity & Rayleigh-Taylor Instability): यह योजना बिना किसी दृश्य स्प्यूरियस दोलन के, मजबूत शियर (shear), इंटरफेशियल विरूपण और घनत्व-प्रेरित अस्थिरताओं (कम श्यानता वाले शासन में केल्विन-हेल्महोल्ट्ज़ रोल-अप सहित) को संभालती है।

महत्व
यह शोध पत्र दावा करता है कि प्रस्तावित विधियों का परिवार CHNS प्रणालियों के लिए एक व्यावहारिक और कुशल विकल्प प्रदान करता है। एक्सट्रपलेशन और सहायक चरों के माध्यम से गैररेखीयताओं को अलग करके और वक्रता नियमितीकरण के माध्यम से अस्थायी विवक्तीकरण को स्थिर करके, यह विधि द्वितीय-क्रम की सटीकता, बिना शर्त स्थिरता और रैखिक कम्प्यूटेशनल लागत के बीच संतुलन प्राप्त करती है। लेखक इस बात पर जोर देते हैं कि यह योजना मॉडल की मुख्य भौतिक संरचनाओं, जिसमें ऊर्जा क्षय और द्रव्यमान संरक्षण शामिल है, को सुरक्षित रखती है, जो इसे विभिन्न शासन (regimes) में जटिल इंटरफेशियल प्रवाहों के अनुकरण के लिए उपयुक्त बनाती है।