Symmetric structure-preserving discretization of N-phase incompressible fluid mixtures with arbitrary density ratios

तकनीकी सारांश: N-फेज असम्पीड्य तरल मिश्रणों का सममित संरचना-संरक्षण विविक्तीकरण (Symmetric Structure-Preserving Discretization)

समस्या विवरण
डिफ्यूज-इंटरफेस मॉडल, विशेष रूप से नेवियर-स्टोक्स-कैन-हिलियर्ड (NSCH) परिवार, जटिल तरल पदार्थों में इंटरफेशियल गतिकी को दर्शाने के लिए व्यापक रूप से उपयोग किए जाते हैं, जहाँ इंटरफेस को सुचारू संक्रमण परतों (smooth transition layers) के रूप में प्रस्तुत किया जाता है। जबकि द्विचर (दो-फेज) प्रवाह के लिए मजबूत संख्यात्मक विधियाँ मौजूद हैं, मनमाने घनत्व अनुपात वाले $N$-फेज असम्पीड्य मिश्रणों तक इसका विस्तार करना महत्वपूर्ण चुनौतियाँ पेश करता है। मौजूदा संरचना-संरक्षण विधियाँ अक्सर द्विचर सूत्रीकरणों पर निर्भर करती हैं या एक संदर्भ चरण (reference phase) को अलग करती हैं, जिससे भौतिक प्रणाली की अंतर्निहित समरूपता (symmetry) टूट जाती है।

$N$-फेज सेटिंग में मुख्य कठिनाई विविक्तीकरण के स्तर पर कई प्रतिस्पर्धी आवश्यकताओं को एक साथ पूरा करने में निहित है:

1. समरूपता (Symmetry): संख्यात्मक योजना को सभी चरणों के साथ समान व्यवहार करना चाहिए बिना सैचुरेशन बाधा (saturation constraint) के माध्यम से चरों को समाप्त किए या किसी विशेषाधिकार प्राप्त संदर्भ चरण को पेश किए।
2. संरक्षण (Conservation): विधि को चरण-वार आयतन और द्रव्यमान, साथ ही कुल आयतन और कुल द्रव्यमान को संरक्षित करना चाहिए।
3. सैचुरेशन बाधा (Saturation Constraint): प्रत्येक बिंदु और समय चरण पर आयतन अंशों का योग इकाई के बराबर होना चाहिए ($\sum \phi_\alpha = 1$), ताकि गैर-भौतिक रिक्त स्थान या ओवरलैप को रोका जा सके।
4. ऊर्जा अपव्यय (Energy Dissipation): विविक्त समाधान को निरंतर ऊष्मागतिकी (continuum thermodynamics) के अनुरूप एक विविक्त ऊर्जा-अपव्यय नियम का पालन करना चाहिए।
5. मनमाने घनत्व अनुपात (Arbitrary Density Ratios): विधि को तब भी स्थिर और सटीक रहना चाहिए जब चरणों के घनत्व में महत्वपूर्ण अंतर हो, जो घनत्व और द्रव्यमान संतुलन समीकरणों के बीच मजबूत युग्मन (coupling) को जन्म देता है।

कार्यप्रणाली
लेखक $N$-फेज असम्पीड्य NSCH मिश्रण मॉडल के लिए एक पूर्णतः-विविक्त परिमित तत्व विधि (fully-discrete finite element method) प्रस्तावित करते हैं। यह दृष्टिकोण एक एकीकृत मिश्रण-सिद्धांत ढांचे पर आधारित है जहाँ प्रणाली एक एकल मिश्रण मोमेंटम समीकरण और $N$ चरण द्रव्यमान संतुलन समीकरणों द्वारा नियंत्रित होती है।

• तुल्य रूपांतरण (Equivalent Reformulation): संरचना संरक्षण को सुगम बनाने के लिए, लेखक निरंतर समीकरणों के एक तुल्य स्ट्रॉन्ग फॉर्मुलेशन (strong formulation) को व्युत्पन्न करते हैं। यह रूपांतरण विशिष्ट परीक्षण फलनों (test functions/weights) का उपयोग करता है ताकि ऊर्जा विकास को प्राप्त किया जा सके। महत्वपूर्ण रूप से, वे द्रव्यमान संतुलन भारों से गैर-रेखीय गतिज पदों (nonlinear kinetic terms) को हटाने के लिए मोमेंटम समीकरण को संशोधित करते हैं, जिससे यह सुनिश्चित होता है कि सभी परीक्षण फलन परिमित तत्व विविक्तीकरण के लिए उपयुक्त मानक सोबोलेव स्थानों (Sobolev spaces) से संबंधित हैं।
• सममित उपचार (Symmetric Treatment): योजना को सीधे चरण चरों ($\phi_\alpha$) के पूर्ण सेट और एक एकल मिश्रण वेग ($v$) के संदर्भ में तैयार किया गया है, बिना किसी चरण को हटाए या संदर्भ चरण पेश किए। यह चरणों की क्रमपरिवर्तन समरूपता (permutation symmetry) को बनाए रखता है।
• विविक्तीकरण रणनीति (Discretization Strategy):
  • समय विविक्तीकरण (Time Discretization): एक पूर्णतः-विविक्त योजना को टाइम-स्टेपिंग दृष्टिकोण का उपयोग करके बनाया गया है। उत्तलता (convexity) और स्थिरता सुनिश्चित करने के लिए ग्रेडिएंट-मुक्त मुक्त ऊर्जा (gradient-free part of the free energy) के लिए टाइम-एवरेजेस का उपयोग किया जाता है। भौतिक अंतराल $[0, 1]$ से अस्थायी रूप से बाहर जाने की स्थिति को संभालने के लिए घनत्व और मोबिलिटी मैट्रिक्स के लिए रेगुलराइजेशन पेश किए गए हैं।
  • स्थानिक विविक्तीकरण (Spatial Discretization): यह विधि $H^1$-कन्फॉर्मिंग स्पेस का उपयोग करते हुए एक परिमित तत्व सन्निकटन (finite element approximation) नियोजित करती है। वेग स्थान को दबाव स्थान के साथ इन्फ-सुप (inf-sup) स्थिर चुना गया है।
  • प्रमुख संरचनात्मक विशेषताएं: योजना कन्वेक्शन टर्म के लिए एक स्क्यू-सिमेट्रिक (skew-symmetric) रूप और द्रव्यमान एवं मोमेंटम समीकरणों के लिए विशिष्ट विविक्त भारों का उपयोग करती है। यह सटीक विविक्त संरक्षण नियमों और एक विविक्त ऊर्जा-अपव्यय असमानता के व्युत्पत्ति की अनुमति देता है।

मुख्य योगदान

1. सममित $N$-फेज विविक्तीकरण: यह पत्र $N$-फेज असम्पीड्य NSCH मॉडलों के लिए पहली पूर्णतः-विविक्त विधि प्रस्तुत करता है जो सभी चरणों के प्रति वास्तव में सममित है, और संदर्भ-चरण सूत्रीकरणों द्वारा पेश की गई विषमता से बचती है।
2. सटीक संरचनात्मक संरक्षण: लेखक सिद्ध करते हैं कि प्रत्येक विविक्त समाधान निम्नलिखित का पालन करता है:
   • सटीक चरण आयतन संरक्षण।
   • सटीक चरण द्रव्यमान संरक्षण।
   • सटीक कुल आयतन और कुल द्रव्यमान संरक्षण।
   • एक विविक्त ऊर्जा-अपव्यय नियम।
   • यदि प्रारंभिक डेटा के लिए सत्य हो, तो बिंदुवार आयतन-सैचुरेशन बाधा ($\sum \phi_\alpha = 1$) का संरक्षण।
3. मनमाने घनत्व अनुपात: विधि मनमाने घनत्व कंट्रास्ट वाले मिश्रणों के लिए मजबूत है, जो वह क्षेत्र है जहाँ मौजूदा स्कीमें घनत्व और द्रव्यमान संतुलन के बीच युग्मन के कारण संघर्ष करती हैं।
4. मिश्रण-जागरूक निरंतरता (Mixture-Aware Consistency): यह सूत्रीकरण "मिश्रण-जागरूक" गुणों का सम्मान करता है, यह सुनिश्चित करता है कि जब चरणों को मिला दिया जाता है या कोई चरण अनुपस्थित होता है, तो प्रणाली सही ढंग से कम (reduce) हो जाती है।

संख्यात्मक परिणाम
लेखक अपने सैद्धांतिक गुणों को कई संख्यात्मक प्रयोगों के माध्यम से मान्य करते हैं:

• अभिसरण परीक्षण ($N=3$): तीन चरणों और मनमाने घनत्वों वाला एक परीक्षण यह प्रदर्शित करता है कि योजना चुने गए परिमित तत्व स्थानों के लिए अपेक्षित अभिसरण क्रम प्राप्त करती है।
• चरण पृथक्करण ($N=3$): टर्नरी (ternary) चरण पृथक्करण का एक सिमुलेशन परस्पर जुड़े ढांचों और ट्रिपल जंक्शनों के निर्माण को दर्शाता है। परिणाम पुष्टि करते हैं कि कुल ऊर्जा गैर-बढ़ने वाली है और द्रव्यमान/आयतन संरक्षण त्रुटियां मशीन परिशुद्धता (machine precision) के स्तर पर हैं।
• उठता हुआ बुलबुला ($N=2$): विधि का परीक्षण एक उठते हुए बुलबुले के लिए मानक हysing et al. बेंचमार्क के विरुद्ध किया गया है। बुलबुले के द्रव्यमान केंद्र और उठने के वेग के लिए परिणाम शार्प-इंटरफेस कोड और अन्य डिफ्यूज-इंटरफेस सूत्रीकरणों से संदर्भ डेटा के साथ अच्छी तरह मेल खाते हैं। यह योजना सटीक रूप से आयतन संरक्षण बनाए रखती है।
• तीन-द्रव वाला उठते हुए बुलबुले का परीक्षण ($N=3$): दो स्तरित तरल परतों के माध्यम से गैस के बुलबुले के उठने के एक जटिल परिदृश्य का अनुकरण किया गया है। विधि बुलबुले के भौतिक व्यवहार को सही ढंग से पकड़ती है—चाहे वह इंटरफेस पर फंसा रहे या ऊपरी परत में प्रवेश करे, जो कि कैपिलरी रेडियस के सापेक्ष बुलबुले की त्रिज्या पर निर्भर करता है।

महत्व और दावे
यह पत्र जटिल इंटरफेशियल गतिकी और मनमाने घनत्व कंट्रास्ट वाले असम्पीड्य $N$-फेज मिश्रण प्रवाह के लिए एक मजबूत कम्प्यूटेशनल फ्रेमवर्क प्रदान करने का दावा करता है। इसकी प्राथमिक महत्ता विषमता को त्याग बिना, निरंतर समीकरणों के परिभाषित ऊष्मागतिक और संरक्षण संरचनाओं को पूर्णतः-विविक्त स्तर पर बनाए रखने की क्षमता में निहित है।

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

लेखक निष्कर्ष निकालते हैं कि हालांकि वर्तमान कार्य सममित संरचना-संरक्षण विविक्तीकरण के लिए आधार स्थापित करता है, भविष्य के कार्यों में इन विचारों को उच्च-क्रम के टेम्पोरल डिस्क्रीटाइजेशन, अधिक सामान्य मिश्रण-जागरूक क्लोजर और संपीड्य $N$-फेज प्रवाह तक विस्तारित किया जा सकता है।