A Finite Element Method for Fluctuating Navier--Stokes Equations

तकनीकी सारांश: उतार-चढ़ाव वाले नेवियर-स्टोक्स समीकरणों के लिए एक परिमित तत्व विधि (Finite Element Method)

समस्या विवरण
तापीय उतार-चढ़ाव (Thermal fluctuations) उन भौतिक घटनाओं के महत्वपूर्ण चालक हैं जो ब्राउनियन गति से लेकर कैविटेशन (cavitation), उबाल (boiling) और हाइड्रोडायनामिक अस्थिरता तक फैली हुई हैं। हालांकि मैक्रोस्कोपिक स्तर पर ये अक्सर नगण्य होते हैं, हाल के अध्ययन संकेत देते हैं कि तापीय शोर (thermal noise) और विक्षोभ (turbulence) के बीच का युग्मन, विशेष रूप से कोलमोगोरोव लंबाई पैमाने (Kolmogorov length scale) के पास, प्रणाली के व्यवहार को महत्वपूर्ण रूप से बदल सकता है। इन प्रभावों को मॉडल करने के लिए, शास्त्रीय नेवियर-स्टोक्स समीकरणों को एक स्टोकेस्टिक स्ट्रेस टेंसर के साथ बढ़ाया जाता है, जिससे उतार-चढ़ाव वाला हाइड्रोडायनामिक्स ढांचा (fluctuating hydrodynamics framework) बनता है। यह ढांचा फ्लक्चुएशन-डिसिपेशन थ्योरम (FDT) पर निर्भर करता है, जो यह निर्धारित करता है कि स्टोकेस्टिक फोर्सिंग का सहप्रसरण (covariance), विस्कस डिसिपेशन ऑपरेटर के समानुपाती होना चाहिए ताकि यह सुनिश्चित हो सके कि प्रणाली सही सांख्यिकीय संतुलन तक पहुँचती है।

कंप्यूटेशनल फ्लूइड डायनेमिक्स में प्राथमिक चुनौती इन समीकरणों को विविक्तिकृत (discretize) करना है जबकि विविक्तीकरण स्तर पर FDT को संरक्षित रखा जा सके। मानक संख्यात्मक दृष्टिकोण, जैसे कि फाइनाइट-डिफरेंस या फाइनाइट-वॉल्यूम विधियाँ, जटिल ज्यामिति (complex geometries) के साथ संघर्ष करती हैं। इसके विपरीत, जबकि फाइनाइट एलीमेंट मेथड (FEM) जटिल डोमेन को संभालने के लिए मजबूत विकल्प प्रदान करता है, मानक गैलरकिन विविक्तीकरण (Galerkin discretizations) उतार-चढ़ाव क्षेत्रों में अवास्तविक, दीर्घ-श्रेणी के सहसंबंध (long-range correlations) उत्पन्न करते हैं। यह आर्टिफैक्ट बेसिस फंक्शन्स के वैश्विक समर्थन (global support) और परिणामी गैर-विकर्ण द्रव्यमान मैट्रिसेस (non-diagonal mass matrices) से उत्पन्न होता है, जो स्थानीय संतुलन सांख्यिकी को विकृत कर देता है।

कार्यप्रणाली
लेखक एक समतापी (isothermal), संपीड़ित (compressible) उतार-चढ़ाव वाले नेवियर-स्टोक्स समीकरणों के लिए एक नवीन, संरचना-संरक्षण (structure-preserving) फाइनाइट एलीमेंट फ्रेमवर्क प्रस्तावित करते हैं। कार्यप्रणाली तीन मुख्य स्तंभों पर आधारित है:

1. वेरिएशनल फॉर्मुलेशन और स्टोकेस्टिक फोर्सिंग: यह विधि स्टोकेस्टिक फोर्सिंग टर्म को सीधे वीक फॉर्मुलेशन (weak formulation) के भीतर परिभाषित करती है। विवेकीकृत स्टोकेस्टिक फोर्स का सहप्रसरण स्पष्ट रूप से विवेकीकृत विस्कस डिसिपेशन ऑपरेटर ($a_h(\cdot, \cdot)$) के समानुपाती बनाया गया है। यह सुनिश्चित करता है कि विवेकीकृत प्रणाली निरंतर (continuum) FDT की नकल करे, जिससे अवास्तविक ऊर्जा इंजेक्शन या कृत्रिम अवमंदन (damping) को रोका जा सके।
2. नोडल क्वाड्रैचर (मास लंपिंग): मानक FEM में निहित अवास्तविक मेश-स्केल सहसंबंधों को समाप्त करने के लिए, लेखक प्रमुख आंतरिक उत्पादों (inner products) के लिए एक नोडल क्वाड्रैचर नियम (ट्रैपेज़ॉइडल नियम) का उपयोग करते हैं। यह तकनीक विकर्ण, लंपेड मास मैट्रिसेस प्रदान करती है। फलस्वरूप, विवेकीकृत साम्यावस्था सहप्रसरण (equilibrium covariance) नोड्स पर विकर्ण होता है, जिससे यह सुनिश्चित होता है कि उतार-चढ़ाव स्थानीय और भौतिक रूप से यथार्थवादी बने रहें, जो प्रभावी रूप से फाइनाइट-वॉल्यूम विधियों के सांख्यिकीय गुणों की नकल करते हुए FEM के वेरिएशनल लाभों को बनाए रखते हैं।
3. समय एकीकरण (Time Integration): समय विविक्तीकरण एक सेमी-इम्प्लिसिट $\theta$-स्कीम का उपयोग करता है। लेखक विशेष रूप से क्रैंक-निकोल्सन स्कीम ($\theta = 0.5$) का विश्लेषण और सत्यापन करते हैं। रैखिकीकृत (linearized) समीकरणों का गणितीय विश्लेषण दर्शाता है कि क्रैंक-निकोल्सन स्कीम अद्वितीय है क्योंकि यह विवेकीकृत सहप्रसरण संतुलन समीकरण में $\Delta t$-निर्भर त्रुटि पद (error term) को शून्य कर देती है। यह इसे बिना किसी प्रतिबंधात्मक स्थिरता स्थितियों के, जो समय चरण (time step) को स्थानिक विविक्तीकरण से जोड़ती हैं, स्थिर फ्लक्चुएशन-डिसिपेशन संतुलन को संतुष्ट करने की अनुमति देता है।

प्रमुख योगदान

• संरचना-संरक्षण विविक्तीकरण: यह शोध पत्र एक ऐसी फाइनाइट एलीमेंट विधि प्रस्तुत करता है जो विवेकीकृत स्तर पर फ्लक्चुएशन-डिसिपेशन संतुलन को कड़ाई से संरक्षित करती है, जिससे स्टोकेस्टिक फोर्सिंग सहप्रसरण को विवेकीकृत विस्कस ऑपरेटर के साथ संरेखित किया जाता है।
• कृत्रिम सहसंबंधों का उन्मूलन: नोडल क्वाड्रैचर का उपयोग करके, यह विधि साम्यावस्था सहप्रसरण मैट्रिक्स में गैर-भौतिक ऑफ-डायगोनल सहसंबंधों को सफलतापूर्वक हटा देती है, जो आमतौर पर निरंतर गैलरकिन FEM दृष्टिकोणों को प्रभावित करते हैं।
• क्रैंक-निकोल्सन का सत्यापन: लेखक यह सैद्धांतिक प्रमाण प्रदान करते हैं कि रैखिक शासन (linearized regime) में विवेकीकृत FDT संतुलन को बनाए रखने के लिए क्रैंक-निकोल्सन स्कीम सबसे इष्टतम विकल्प है, क्योंकि यह सहप्रसरण संतुलन में टाइम-स्टेप-निर्भर विचलन को समाप्त करती है।
• व्यापक सत्यापन: इस फ्रेमवर्क को एक, दो और तीन स्थानिक आयामों में त्रिकोणीय (triangular), चतुर्भुज (quadrilateral) और टेट्राहेड्रल (tetrahedral) मेश का उपयोग करके संख्यात्मक रूप से सत्यापित किया गया है।

परिणाम
संख्यात्मक प्रयोग पुष्टि करते हैं कि प्रस्तावित विधि विभिन्न विविक्तीकरण मापदंडों में साम्यावस्था उतार-चढ़ाव सांख्यिकी को सही ढंग से पुनरुत्पादित करती है:

• वैरिएंस स्केलिंग (Variance Scaling): घनत्व और वेग क्षेत्रों के अनुभवजन्य वैरिएंस सांख्यिकीय यांत्रिकी द्वारा अनुमानित सैद्धांतिक साम्यावस्था मानों की ओर अभिसरित (converge) होते हैं। परिणाम सही स्केलिंग व्यवहार ($\langle \delta (\cdot)^2 \rangle \propto |T^i_h|^{-1}$) प्रदर्शित करते हैं, जहाँ मेश आकार (या तत्व आयतन) को आधा करने से वैरिएंस लगभग दोगुना हो जाता है।
• मजबूती (Robustness): विधि विभिन्न मेश रिज़ॉल्यूशन, टाइम स्टेप और तत्व प्रकारों (2D में त्रिकोणीय बनाम चतुर्भुज) में लगातार प्रदर्शन करती है।
• द्रव्यमान संरक्षण (Mass Conservation): एक निरंतर गैलरकिन विविक्तीकरण का उपयोग करने के बावजूद, यह योजना कुल द्रव्यमान को संरक्षित करती है, जो एक महत्वपूर्ण गुण है जिसे अक्सर अपविंडिंग (upwinding) या अन्य स्थिरीकरण तकनीकों द्वारा प्रभावित किया जाता है जो फ्लक्चुएशन-डिसिपेशन संतुलन को बदल सकती हैं।
• सीमाएँ: लेखक नोट करते हैं कि बहुत कम औसत घनत्व (जैसे, 3D में $\rho \approx 10^{-2}$) के लिए, रैखिकीकृत साम्यावस्था सूत्र कम वैध हो जाते हैं, और यह विधि स्पष्ट रूप से धनात्मक घनत्व ($\rho \geq 0$) के भौतिक प्रतिबंध को लागू नहीं करती है। हालांकि, विधि मध्यम उतार-चढ़ाव आयामों के लिए सकारात्मकता-संरक्षण (positivity-preserving) फाइनाइट-डिफरेंस विधियों के तुलनीय परिणाम सफलतापूर्वक पुनरुत्पादित करती है।

महत्व और दावे
यह शोध पत्र दावा करता है कि यह जटिल ज्यामिति पर संपीड़ित तरल पदार्थों में तापीय उतार-चढ़ाव के अनुकरण के लिए एक मजबूत और सटीक उपकरण प्रदान करता है, एक ऐसा क्षेत्र जहाँ पारंपरिक फाइनाइट-डिफरेंस या फाइनाइट-वॉल्यूम विधियाँ चुनौतियों का सामना करती हैं। यह सुनिश्चित करके कि विवेकीकृत योजना मौलिक फ्लक्चुएशन-डिसिपेशन संतुलन का सम्मान करती है और कृत्रिम सहसंबंधों को समाप्त करती है, यह विधि मानक FEM अनुप्रयोगों में सामान्य आर्टिफैक्ट्स के बिना साम्यावस्था सांख्यिकी को कैप्चर करने का एक विश्वसनीय तरीका प्रदान करती है।

लेखक इस कार्य को पूरी तरह से जुड़े हुए, गैर-रैखिक लैंडौ-लिफ़शिट्ज़ (Landau–Lifshitz) उतार-चढ़ाव वाले नेवियर-स्टोक्स समीकरणों तक फाइनाइट एलीमेंट तकनीकों को विस्तारित करने की दिशा में एक महत्वपूर्ण कदम के रूप में देखते हैं। वे इस बात पर जोर देते हैं कि हालांकि वर्तमान कार्य समतापी, एकल-घटक तरल पदार्थों पर केंद्रित है, लेकिन विवेकीकृत ढांचा बहु-घटक और गैर-समतापी प्रणालियों के साथ-साथ एडेप्टिव मेश रिफाइनमेंट और जटिल सीमाओं वाले अनुप्रयोगों के लिए भी लचीला है। शोध पत्र निष्कर्ष निकालता है कि यह विधि वास्तविक सेटिंग्स में उतार-चढ़ाव वाले हाइड्रोडायनामिक्स के मल्टीस्केल सिमुलेशन के लिए एक शक्तिशाली उपकरण है, बशर्ते कि भविष्य के कार्य सशक्त उतार-चढ़ाव के लिए सकारात्मकता प्रतिबंधों (positivity constraints) को लागू करने के प्रश्न को संबोधित करें।