An all-topology two-fluid model for two-phase flows derived through Hamilton's Stationary Action Principle

1. समस्या विवरण (Problem Statement)

यह शोध पत्र संपीड्य दो-चरणीय प्रवाह (compressible two-phase flows) के लिए एक सुदृढ़ गणितीय मॉडल विकसित करने की दीर्घकालिक चुनौती को संबोधित करता है, जो सभी प्रकार के प्रवाह टोपोलॉजी (पृथक से लेकर विसरित शासन तक) में मान्य हो और शॉक (shocks) तथा टोपोलॉजी परिवर्तनों (जैसे कि फेज इनवर्जन, एटोमाइजेशन) को संभालने में सक्षम हो।

मौजूदा मॉडलों की मुख्य सीमाएँ निम्नलिखित हैं:

• एकल-दबाव मॉडल (Single-pressure models): ये सामान्यतः गैर-हाइपरबोलिक होते हैं, जिससे 'इल-पोज़्डनेस' (ill-posedness) और संख्यात्मक अस्थिरता उत्पन्न होती है।
• बेयर-नुनज़ियाटो (BN) मॉडल: हालांकि ये लोकप्रिय और हाइपरबोलिक हैं, लेकिन इन्हें इंटरफेशियल मात्राओं (इंटरफेशियल वेग $u_I$ और दबाव $p_I$) के लिए क्लोजर धारणाओं (closure assumptions) की आवश्यकता होती है।
  • मानक क्लोजर (जैसे $u_I$ = एक चरण का वेग) जटिल टोपोलॉजी के लिए विफल हो जाते हैं।
  • "ऑल-टोपोलॉजी" क्लोजर (जैसे द्रव्यमान-भारित औसत वेग) अक्सर ऐसे मॉडल की ओर ले जाते हैं जिनमें एक स्पष्ट एंट्रॉपी संरक्षण नियम का अभाव होता है या जिन्हें ऐसे इंटरफेशियल प्रेशर एक्सप्रेशन की आवश्यकता होती है जो चरणों के तापमान पर निर्भर करते हैं (जैसे कि पेपर में समीकरण 1.1), जिसमें मोमेंटम एक्सचेंज के लिए भौतिक औचित्य की कमी होती है।
• ऊष्मागतिक रूप से संगत प्रणालियाँ (Thermodynamically Compatible Systems - TCS): ये मॉडल अच्छा गणितीय ढांचा प्रदान करते हैं लेकिन पूर्ण ऊष्मागतिकी तक विस्तार करने पर मोमेंटम समीकरणों में कृत्रिम हीट फ्लक्स या गैर-भौतिक पद (terms) पेश कर सकते हैं।

मुख्य समस्या एक ऐसे मॉडल के अभाव की है जिसमें साथ में निम्नलिखित गुण हों:

1. गणितीय सु-संरचना (Mathematical well-posedness): हाइपरबोलिसिटी, सिमेट्राइज़ेबिलिटी (symmetrizability), और शॉक के लिए स्पष्ट रूप से परिभाषित जंप स्थितियाँ।
2. भौतिक निरंतरता (Physical consistency): प्रथम सिद्धांतों से व्युत्पन्न इंटरफेशियल क्लोजर, जो बिना किसी मनमाने अनुमान के पदों की स्पष्ट भौतिक व्याख्या सुनिश्चित करता है (कोई कृत्रिम ऊष्मा विनिमय नहीं)।

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

लेखक हैमिल्टन के स्टेशनरी एक्शन सिद्धांत (Hamilton's Stationary Action Principle) पर आधारित एक नवीन व्युत्पत्ति ढांचा प्रस्तावित करते हैं, जिसे दो परिवारों के प्रक्षेप पथों (trajectories) (प्रत्येक चरण के लिए एक) तक विस्तारित किया गया है।

• वैरियेशनल फ्रेमवर्क (Variational Framework): एकल-प्रक्षेप पथ दृष्टिकोण (जो मिश्रण वेग के लिए एक आभासी प्रक्षेप पथ मानता है) के विपरीत, यह विधि दोनों चरणों के प्रक्षेप पथों ($\Phi_1$ और $\Phi_2$) को स्वतंत्र रूप से मानती है।
• लैग्रेंजियन फॉर्मूलेशन (Lagrangian Formulation): एक जेनेरिक लैग्रेंजियन डेंसिटी बनाई गई है जिसमें गतिज ऊर्जा, आंतरिक ऊर्जा और बाधाएं (constraints) शामिल हैं।
• बाधा प्रबंधन (Constraint Handling):
  • वॉल्यूम फ्रैक्शन (Volume fraction) बाधा को द्रव्यमान-भारित औसत वेग ($u = Y_1 u_1 + Y_2 u_2$) को इंटरफेशियल वेग के रूप में लागू करने के लिए संशोधित किया गया है।
  • चूंकि इस बाधा को व्यक्तिगत चरण प्रक्षेप पथों के साथ सटीक रूप से एकीकृत नहीं किया जा सकता है, इसलिए एक लैग्रेंज मल्टीप्लायर ($\zeta$) पेश किया गया है। इससे $\zeta$ के लिए एक नया विकास समीकरण और एक नया चर $v = \nabla \zeta$ प्राप्त होता है।
• व्युत्पत्ति (Derivation): इन विशिष्ट बाधाओं के साथ एक्शन इंटीग्रल पर वैरिएशनल प्रिंसिपल लागू करके, लेखक आly-लेग्रेज समीकरण (Euler-Lagrange equations) व्युत्पन्न करते हैं। यह प्रक्रिया स्वाभाविक रूप से गवर्निंग समीकरणों और इंटरफेशियल मात्राओं के क्लोजर को बिना किसी मनमाने पोस्टुलेशन के प्रदान करती है।

3. मुख्य योगदान (Key Contributions)

A. नया ऑल-टोपोलॉजी मॉडल

यह पेपर एक पूर्णतः क्लोज्ड टू-प्रेशर, टू-फ्लुइड मॉडल व्युत्पन्न करता है जिसमें निम्नलिखित विशेषताएं हैं:

• इंटरफेशियल वेग: द्रव्यमान-भारित औसत वेग ($u_I = u$) के रूप में परिभाषित है, जो सभी प्रवाह शासनों में वैधता सुनिश्चित करता है।
• इंटरफेशियल प्रेशर ($p_I$): इसे $p_I = Y_2 p_1 + Y_1 p_2$ (द्रव्यमान-भारित औसत) के रूप में व्युत्पन्न किया गया है। यह अभिव्यक्ति भौतिक रूप से सुदृढ़ है (बाहरी दबाव बलों का प्रतिनिधित्व करती है) और यह चरणों के तापमान पर निर्भर नहीं है।
• इंटरफेशियल वर्क ($(pu)_I$): वैरियेशनल प्रिंसिपल द्वारा पेश किया गया एक नया परिमाण: $(pu)_I = Y_1 p_2 u_1 + Y_2 p_1 u_2$। यह पद $p_I u_I$ से भिन्न है और ऊर्जा समीकरणों को सही ढंग से क्लोज करने के लिए महत्वपूर्ण है।
• लिफ्ट फोर्सेस (Lift Forces): बहु-आयामी सेटिंग्स में, यह व्युत्पत्ति स्वाभाविक रूप से लैग्रेंज मल्टीप्लायर ($\nabla \zeta$) द्वारा संचालित लिफ्ट जैसे बल उत्पन्न करती है। ये बल दबाव गैर-संतुलन (pressure non-equilibrium) से जुड़े हैं।

B. गणितीय गुण (1D विश्लेषण)

लेखक एक-आयामी सेटिंग में मॉडल का कठोर विश्लेषण करते हैं (विश्लेषण के लिए लिफ्ट बलों को नजरअंदाज करते हुए):

• हाइपरबोलिसिटी (Hyperbolicity): सिस्टम को मानक गैर-अनुनाद स्थिति ($u \neq u_k \pm c_k$) के तहत हाइपरबोलिक सिद्ध किया गया है।
• सिमेट्राइज़ेबिलिटी (Symmetrizability): सिस्टम सिमेट्राइज़ेबल है, जो स्मूथ इनिशियल डेटा के लिए स्थानीय सु-संरचना (local well-posedness) सुनिश्चित करता है।
• एंट्रॉपी संरक्षण (Entropy Conservation): मॉडल कुल एंट्रॉपी के लिए एक पूरक संरक्षण नियम स्वीकार करता है। यह एक महत्वपूर्ण परिणाम है, क्योंकि यह द्वितीय नियम के अनुकूल भौतिक रूप से स्वीकार्य कमजोर समाधानों (शॉक्स) के चयन की अनुमति देता है।
• जंप स्थितियाँ (Jump Conditions): इंटरफेशियल वेव की लीनियर डीजनरेसी के कारण, विच्छिन्नता (discontinuities) के पार गैर-संरक्षण उत्पाद (non-conservative products) अद्वितीय रूप से परिभाषित हैं, जो रैंकिन-ह्यूगोनियट जंप स्थितियों का एक अनूठा सेट प्रदान करते हैं।

C. भौतिक निरंतरता

• मॉडल TCS मॉडलों में पाए जाने वाले "कृत्रिम ऊष्मा विनिमय" से बचता है।
• इंटरफेशियल प्रेशर और वर्क पद शुद्ध रूप से मैकेनिक्स और थर्मोडायनामिक्स से व्युत्पन्न हैं, जो तापमान-निर्भर प्रेशर क्लोजर से बचते हैं जो पिछले ऑल-टोपोलॉजी BN मॉडलों की समस्या थी।
• फेसिक एंट्रॉपी अपने स्वयं के वेगों ($D_k s_k = 0$) पर ट्रांसपोर्ट होती है, जो प्रतिवर्ती प्रवाह (reversible flow) धारणाओं के अनुरूप है।

4. परिणाम (Results)

• व्युत्पत्ति की सफलता: वैरियेशनल प्रिंसिपल ने सिस्टम को सफलतापूर्वक क्लोज किया, जिससे भौतिक रूप से सुदृढ़ अभिव्यक्तियाँ ($p_I$) प्राप्त हुईं और आवश्यक $(pu)_I$ पद पेश हुआ।
• शॉक उपचार (Shock Treatment): मॉडल अन्य BN वेरिएंट में पाई जाने वाली शॉक जंप स्थितियों की अस्पष्टता को हल करता है। इंटरफेशियल वेव की लीनियर डीजनरेसी यह सुनिश्चित करती है कि गैर-संरक्षण पद अद्वितीय समाधानों की ओर न ले जाएं।
• बहु-आयामी प्रभाव: मॉडल सापेक्ष वेग और दबाव अंतर क्षमता के ग्रेडिएंट ($v$) के अनुपात में लिफ्ट बलों की भविष्यवाणी करता है। हालांकि $v$ की भौतिक व्याख्या के लिए आगे के अध्ययन की आवश्यकता है, मॉडल 1D मामले में सही ढंग से घटित होता है जहाँ ये बल शून्य हो जाते हैं।
• सोर्स टर्म्स (Source Terms): मॉडल मानक रिलैक्सेशन सोर्स टर्म्स (मैकेनिकल, किनेमैटिक, थर्मल) के साथ संगत है जो सिस्टम को इक्विलिब्रियम की ओर ले जाते हैं, जिससे एसिम्प्टोटिक लिमिट में ज्ञात रिड्यूस्ड-ऑर्डर मॉडल (जैसे कि कपिला मॉडल) प्राप्त होते हैं।

5. महत्व (Significance)

यह कार्य दो-चरणीय प्रवाह के गणितीय मॉडलिंग में एक महत्वपूर्ण प्रगति है:

1. दृष्टिकोणों का एकीकरण: यह गणितीय दृष्टिकोण (सु-संरचना और एंट्रॉपी संरक्षण की आवश्यकता) को भौतिक दृष्टिकोण (पदों की स्पष्ट यांत्रिक व्याख्या की आवश्यकता) के साथ मेल खाता है।
2. कठोर आधार के साथ पहला ऑल-टोपोलॉजी मॉडल: लेखकों की जानकारी के अनुसार, यह पहला ऑल-टोपोलॉजी मल्टी-फ्लुइड मॉडल है जिसमें पूर्ण भौतिक निरंतरता और शॉक-प्रकार के समाधानों का उचित वर्णन करने के लिए आवश्यक कठोर गणितीय संरचना दोनों मौजूद हैं।
3. न्यूमेरिक्स के लिए आधार: अद्वितीय जंप स्थितियों का व्युत्पत्ति और एंट्रॉपी संरक्षण कानून का अस्तित्व, शॉक और टोपोलॉजी परिवर्तन वाले जटिल दो-चरणीय प्रवाह के लिए मजबूत संख्यात्मक स्कीम्स (जैसे कि एप्रोक्सिमेट रीमैनन सॉल्वर) विकसित करने के लिए एक ठोस सैद्धांतिक आधार प्रदान करता है।
4. नई भौतिक अंतर्दृष्टि: एक अलग चर के रूप में "इंटरफेशियल वर्क" का परिचय और दबाव गैर-संतुलन द्वारा संचालित लिफ्ट बलों की पहचान, मल्टी-फेज डायनेमिक्स को समझने के नए रास्ते खोलती है।

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