Exact expressions for the unresolved stress in a finite-volume based large-eddy simulation

तकनीकी सारांश: एक परिमित-आयतन आधारित लार्ज-एडी सिमुलेशन (LES) में अनसुलझे तनाव के सटीक व्यंजक

समस्या विवरण
लार्ज-एडी सिमुलेशन (LES) अशांत प्रवाह (turbulent flows) को मॉडल करने के लिए बड़े पैमानों को हल करने और अनसुलझे उप-ग्रिड पैमानों (sub-grid scales) के प्रभावों को मॉडल करने का लक्ष्य रखता है। पारंपरिक रूप से, LES में क्लोजर समस्या (closure problem) स्थानिक फ़िल्टरिंग और गैर-रेखीय संवहन पदों (nonlinear advection terms) के बीच गैर-विनिमेयता (non-commutation) से उत्पन्न होने वाले कम्यूटेटर त्रुटि (commutator error) पर केंद्रित होती है, जिसके परिणामस्वरूप शास्त्रीय अवशिष्ट स्ट्रेस टेंसर (residual stress tensor, RST), $\tau_{ij} = \overline{u_i u_j} - \bar{u}_i \bar{u}_j$ प्राप्त होता है। हालाँकि, व्यावहारिक सिमुलेशन में, LES को संख्यात्मक विविक्तीकरण योजनाओं (numerical discretization schemes) का उपयोग करके हल किया जाता है, जो सबसे अधिकतः 'फिनाइट वॉल्यूम मेथड' (FVM) है।

एक महत्वपूर्ण मुद्दा यह उठता है कि LES में फ़िल्टर की चौड़ाई $\Delta$ अक्सर ग्रिड स्पेसिंग $h$ के समान क्रम की होती है (अर्थात $\Delta \approx h$)। ऐसे शासन (regimes) में, विविक्तीकरण त्रुटियाँ (discretization errors)—जो संख्यात्मक फ्लक्स (numerical fluxes) और विविक्त डाइवर्जेंस ऑपरेटरों (discrete divergence operators) से उत्पन्न होती हैं—शास्त्रीय उप-फ़िल्टर स्ट्रेस के समान ही क्रम की होती हैं। शास्त्रीय दृष्टिकोण आमतौर पर निरंतर समीकरणों के आधार पर क्लोजर मॉडल को परिभाषित करते हैं और फिर उसे विविक्त (discretize) करते हैं, या विविक्तीकरण त्रुटियों को अलग अभिसरण (convergence) समस्याओं के रूप में देखते हैं। यह पृथक्करण लक्ष्य डेटा (वास्तविक अवशिष्ट) और मॉडल के बीच एक बेमेल स्थिति पैदा करता है, विशेष रूप से जब विविक्तीकरण स्वयं महत्वपूर्ण गैर-स्थानीय (non-local) या गैर-सममित (non-symmetric) प्रभाव उत्पन्न करता है जिसे मानक सममित, स्थानीय क्लोजर मॉडल (जैसे स्मैगोरिंस्की) कैप्चर नहीं कर सकते।

कार्यप्रणाली
लेखक एक एकीकृत गणितीय ढांचा प्रस्तावित करते हैं जिसे LES-FVM कहा जाता है, जो LES फ़िल्टरिंग प्रक्रिया को सीधे FVM विविक्तीकरण के साथ एकीकृत करता है। मुख्य कार्यप्रणाली इस प्रकार है:

1. एकीकृत औपचारिकता (Unified Formalism): निरंतर समीकरणों पर LES फ़िल्टरिंग लागू करने और फिर उसे विविक्त करने के बजाय, लेखक LES-फ़िल्टर्ड समीकरणों पर FVM फ़िल्टर (नियंत्रित आयतनों पर एक वॉल्यूम एवरेज) लागू करते हैं। यह एक "डबल-फ़िल्टर्ड" प्रणाली बनाता है जहाँ फ़िल्टर की चौड़ाई $\Delta$ और ग्रिड स्पेसिंग $h$ को स्वतंत्र मापदंडों के रूप में माना जाता है।
2. सटीक अवशिष्टों का व्युत्पन्न (Derivation of Exact Residuals): एक "डिस्क्रीट फ़िल्टर-स्वैप" गुण (जहाँ फाइनाइट डिफरेंस ऑपरेटर एक फ़िल्टर्ड सटीक डेरिवेटिव के समकक्ष है) का उपयोग करते हुए, लेखक विविक्त समीकरणों में अवशिष्ट फ्लक्स के सटीक व्यंजक प्राप्त करते हैं।
3. अवसादों का अपघटन (Decomposition of Residuals): कुल अवशिष्ट फ्लक्स $\tau_{\Delta,h}$ को तीन विशिष्ट घटकों में विभाजित किया गया है:
   • क्लासिक टर्म (Classic Term): डबल फ़िल्टर ($\Delta$ और $h$) से जुड़ा मानक कम्यूटेटर एरर।
   • फ्लक्स टर्म (Flux Term): सटीक फ्लक्स और संख्यात्मक फ्लक्स ($r$ बनाम $r_h$) के बीच अंतर से उत्पन्न त्रुटियाँ।
   • डाइवर्जेंस टर्म (Divergence Term): सटीक डाइवर्जेंस और विविक्त डाइवर्जेंस के बीच अंतर से उत्पन्न त्रुटियाँ।
4. 3D असंपीडित प्रवाह (Incompressible Flow) का विस्तार: इस ढांचे को 3D असंपीडित नेवियर-स्टोक्स समीकरणों तक विस्तारित किया गया है। इसके लिए वेग क्षेत्र को विविक्त सेटिंग में डाइवर्जेंस-मुक्त बनाए रखने हेतु प्रोजेक्शन ऑपरेटरों के माध्यम से दबाव पद (pressure term) को संभालना आवश्यक है। लेखक एक 'प्रोजेक्टेड RST' व्युत्पन्न करते हैं जो दबाव प्रोजेक्शन और विविक्तीकरण के बीच की अंतःक्रिया को ध्यान में रखता है।

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

• सटीक विविक्तीकरण-सूचित RST (Exact Discretization-Informed RST): यह शोध पत्र एक FVM-आधारित LES के लिए अनसुलझे स्ट्रेस टेंसर के लिए एक सटीक व्यंजक व्युत्पन्न करता है। शास्त्रीय RST के विपरीत, यह नया व्यंजक स्पष्ट रूप से संख्यात्मक फ्लक्स और विविक्त डाइवर्जेंस ऑपरेटरों के योगदान को शामिल करता है।
• गैर-सममिति और गैर-स्थानीयता (Non-Symmetry and Non-Locality): असंपीडित प्रवाह के लिए, लेखक प्रदर्शित करते हैं कि प्रस्तावित LES-FVM RST स्वाभाविक रूप से गैर-सममित और गैर-स्थानीय (प्रेशर प्रोजेक्शन के कारण) है। यह शास्त्रीय LES RST के विपरीत है, जिन्हें आमतौर पर सममित और स्थानीय माना जाता है।
• शून्य ए-पोस्टीरियोरियस त्रुटि (Zero A-Posteriori Error): "DNS-सहायता प्राप्त" प्रयोगों में (जहाँ डायरेक्ट न्यूमेरिकल सिमुलेशन डेटा से सटीक अवशिष्ट की गणना की जाती है और उसे LES समीकरण में इंजेक्ट किया जाता है), प्रस्तावित RST व्यंजक शून्य ए-पोस्टीरियोरियस त्रुटि देता है। इसके विपरीत, शास्त्रीय RST व्यंजक त्रुटियों को समय के साथ संचित करते हैं।
• बेहतर क्लोजर मॉडलिंग: जब नए RST का उपयोग स्मैगोरिंस्की मॉडल के गुणांक को फिट करने के लिए किया जाता है, तो प्राप्त गुणांक शास्त्रीय RST की तुलना में अधिक होता है। यह ऊर्जा स्पेक्ट्रा और सापेक्ष त्रुटि के संदर्भ में बेहतर प्रदर्शन की ओर ले जाता है, विशेष रूप से छोटे फ़िल्टर-टू-ग्रिड अनुपात ($\Delta/h$) के लिए।

परिणाम
लेखक अपने सिद्धांत को दो परीक्षण मामलों का उपयोग करके मान्य करते हैं: 1D विस्कस बर्गर्स समीकरण और 3D डीकेइंग आइसोट्रोपिक टर्बुलेंस।

• 1D बर्गर्स समीकरण:

  • ए-प्रायोरी विश्लेषण (A-Priori Analysis): अध्ययन दिखाता है कि सामान्य फ़िल्टर अनुपातों ($\Delta \le 4h$) के लिए, विविक्तीकरण त्रुटियाँ (फ्लक्स और डाइवर्जेंस पद) कुल अवशिष्ट के एक महत्वपूर्ण हिस्से (70-70% तक) का निर्माण करती हैं। इन पदों की उपेक्षा करने से वास्तविक अवशिष्ट का खराब प्रतिनिधित्व होता है।
  • ए-पोस्टीरियोरियस विश्लेषण (A-Posteriori Analysis): क्लोजर टर्म के रूप में सटीक विविक्तीकरण-सूचित अवशिष्ट का उपयोग करने से, फ़िल्टर्ड DNS समाधान का "परफेक्ट" पुनर्निर्माण प्राप्त होता है।
  • स्मैगोरिंस्की मॉडल: नए RST के प्रति स्मैगोरिंस्की गुणांक को फिट करने से शास्त्रीय RST की तुलना में उच्च गुणांक मान प्राप्त होता है। इस सूचित गुणांक का उपयोग करने वाले सिमुलेशन कम सापेक्ष त्रुटियां और बेहतर ऊर्जा स्पेक्ट्रा दिखाते हैं, जो कम $\Delta/h$ अनुपात पर फ़िल्टर्ड DNS के साथ बेहतर मेल खाते हैं।

• 3D नेवियर-स्टोक्स समीकरण:

  • व्युत्पत्ति पुष्टि करती है कि 3D RST गैर-सममित और गैर-स्थानीय है।
  • 3D डीकेइंग टर्बुलेंस पर DNS-सहायता प्राप्त प्रयोग पुष्टि करते हैं कि प्रस्तावित ढांचा अवशिष्ट स्ट्रेस को सही ढंग से कैप्चर करता है, जो पूर्ण असंपीडित प्रवाह के लिए 1D सिद्धांत के विस्तार को मान्य करता है।

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

प्राथमिक महत्व इस तथ्य में निहित है कि विविक्तीकरण-प्रेरित त्रुटियाँ केवल संख्यात्मक शोर नहीं हैं, बल्कि अवशिष्ट स्ट्रेस का अभिन्न अंग हैं जिन्हें क्लोजर मॉडल में शामिल किया जाना चाहिए। लेखक तर्क देते हैं कि शास्त्रीय RST को डेटा-संचालित या ट्यून किए गए क्लोजर मॉडल के लिए लक्ष्य के रूप में उपयोग करना एक मौलिक विसंगति पेश करता है, जिससे उप-इष्टतम मॉडल गुणांक और संचित त्रुटियां होती हैं। उनका प्रस्तावित "विविक्तीकरण-सूचित" RST सही लक्ष्य डेटा प्रदान करता है, जिससे अधिक सटीक और स्थिर क्लोजर मॉडल विकसित करने में सक्षम होता है, विशेष रूप से इम्प्लिसिट LES या मोटे ग्रिड वाले परिदृश्यों के लिए जहाँ $\Delta \approx h$ होता है।

लेखक जटिल औद्योगिक प्रवाह के लिए इन निष्कर्षों के तत्काल अनुप्रयोग के संबंध में विनम्र रहते हैं, यह नोट करते हुए कि हालांकि नया RST सरल स्मैगोरिंकी मॉडल के प्रदर्शन में सुधार करता है, इस ढांचे की पूर्ण क्षमता संभवतः अधिक अभिव्यंजक, गैर-स्थानीय और गैर-सममित क्लोजर मॉडल (जैसे मशीन लर्निंग-आधारित दृष्टिकोण) के साथ साकार होगी जो व्युत्पन्न अवशिष्ट की संरचना का पूरी तरह से लाभ उठा सकते हैं।