Self-similar structure of non-isothermal variable-density mushy Stefan problems and an improved low-Mach enthalpy method

तकनीकी सारांश: गैर-समतापीय परिवर्तनशील-घनत्व वाले मषी (mushy) स्टेफ़न समस्याओं की स्व-समान संरचना और एक उन्नत लो-मैच एन्थैल्पी विधि

समस्या विवरण
एन्थैल्थी विधि (enthalpy method) निश्चित ग्रिडों पर चरण-परिवर्तन (पिघलना और जमना) की समस्याओं के अनुकरण के लिए व्यापक रूप से अपनाई जाने वाली एक पद्धति है। यह आमतौर पर एक परिमित तापमान अंतराल, $\Delta T = T_{liq} - T_{sol}$ के भीतर चरण संक्रमण को मॉडल करती है, जो एक "मषी" (mushy) क्षेत्र बनाती है। हालांकि इस विधि को सामान्यतः यह मानकर उपयोग किया जाता है कि $\Delta T \to 0$ होने पर यह शास्त्रीय शार्प-इंटरफेस स्टेफ़न समस्या की ओर अभिसरित (converge) होती है, लेकिन यह धारणा समान ठोस और तरल घनत्व ($\rho_S = \rho_L$) की सरलीकृत स्थिति के तहत ही मुख्य रूप से मान्य रही है।

महत्वपूर्ण आयतन परिवर्तन (जैसे, उबालना/boiling, संघनन/condensation, या विशिष्ट धातु जमने) वाले व्यावहारिक अनुप्रयोगों में, घनत्व अनुपात $\rho_L/\rho_S$ एक से काफी भिन्न हो सकता है। ऐसे मामलों में, घनत्व परिवर्तन के कारण चरण परिवर्तन से तरल गति उत्पन्न होती है, जो तापीय समस्या को द्रव्यमान और संवेग संरक्षण के साथ जोड़ती है। मौजूदा साहित्य में उच्च-घनत्व-अनुपात वाली समस्याओं के लिए परिमित $\Delta T$ के साथ एन्थैल्थी विधियों की सटीकता और अभिसरण गुणों के व्यवस्थित अन्वेषण का अभाव है। इसके अलावा, शास्त्रीय स्टेफ़न समस्या (अनंत सूक्ष्म $\Delta T$) के विरुद्ध सत्यापन पर्याप्त नहीं हो सकता है, क्योंकि मषी स्टेफ़न समस्या (परिमित $\Delta T$) मौलिक रूप से शास्त्रीय स्टेफ़न समस्या से भिन्न है।

कार्यप्रणाली
लेखक सैद्धांतिक विश्लेषण और संख्यात्मक विकास के संयोजन के माध्यम से इन कमियों को दूर करते हैं:

1. मषी स्टेफ़न समस्या का सैद्धांतिक निरूपण:
   लेखक एक गैर-समतापीय, परिवर्तनशील-घनत्व वाले मषी स्टेफ़न समस्या का प्रतिपादन करते हैं। वे प्रदर्शित करते हैं कि गवर्निंग समीकरणों में एक स्व-समान संरचना (self-similar structure) है, जो आंशिक अंतर समीकरणों (PDEs) को दो अज्ञात समानता मापदंडों ($\lambda_s$ और $\lambda_l$) से संबंधित एक सीमा-मूल्य समस्या (boundary-value problem) में बदल देती है, जो सॉलिडस और लिक्विडस फ्रंट की स्थितियों को परिभाषित करते हैं।

• समान-घनत्व मामला: एक बंद-रूप (closed-form) विश्लेतिक समानता समाधान व्युत्पन्न किया गया है।
• असमान-घनत्व मामला: यह समस्या तापमान और वेग क्षेत्रों को नियंत्रित करने वाली एक युग्मित गैररेखीय प्रणाली में बदल जाती है। इस प्रणाली को एक बाउंडरी-वैल्यू सॉल्वर (bvp5c का उपयोग करके MATLAB) और एक निरंतरता रणनीति (continuation strategy) के साथ हल किया जाता है ताकि अत्यधिक घनत्व अनुपातों ( $\rho_L/\rho_S = 2700$ तक) और छोटे $\Delta T$ को संभाला जा सके।

2. एन्थैल्थी परिभाषाओं का विश्लेषण:
   यह अध्ययन मषी क्षेत्र में विशिष्ट एन्थैल्थी की दो परिभाषाओं का परीक्षण करता है। "मूल" परिभाषा घनत्व भिन्न होने पर तरल अंश ($\phi$) और तापमान ($T$) के बीच एक गैररेखीय संबंध की ओर ले जाती है, जबकि "उन्नत" परिभाषा एक रैखिक $\phi$-$T$ संबंध प्रदान करती है। लेखक दिखाते हैं कि मूल प्रतिरूपण उच्च घनत्व अनुपातों पर अस्वाभाविक अस्थिरता प्रदर्शित करता है, जबकि उन्नत प्रतिरूपण स्थिर रहता है।

3. उन्नत लो-मैच एन्थैल्थी विधि:
   सैद्धांतिक अंतर्दृष्टि के आधार पर, लेखक एक उन्नत लो-मैच एन्थैल्थी विधि प्रस्तावित करते हैं। मुख्य विशेषताएं शामिल हैं:

• एक संशोधित विशिष्ट एन्थैल्थी परिभाषा जो रैखिक तरल-अंश-तापमान संबंध सुनिश्चित करती है।
• एक काइनेमैटिक वेग-विचलन बाधा (kinematic velocity-divergence constraint) जो द्रव्यमान संरक्षण से सीधे प्राप्त की गई है ताकि घनत्व-परिवर्तन-प्रेरित प्रवाह को नियंत्रित किया जा सके।
• इसे IBAMR लाइब्रेरी के भीतर एक स्टैगर्ड कार्टेशियन ग्रिड, गैररेखीय एन्थैल्थी समीकरण के लिए न्यूटन इटरेशन, और ठोस क्षेत्र में कठोर कारमन-कोज़ेनी (Carman–Kozeny) ड्रैग फोर्स को संभालने के लिए एक प्रेडिक्शन मेथड के साथ लागू किया गया है।

मुख्य परिणाम

• स्व-समानता और अभिसरण: विश्लेषण पुष्टि करता है कि एन्थैल्थी विधि मषी स्टेफ़न समस्या को हल करती है, न कि शास्त्रीय स्टेफ़न समस्या को। मषी क्षेत्र की चौड़ाई $\sqrt{t}$ के समानुपाती रूप से बढ़ती है, जबकि शास्त्रीय स्टेफ़न फ्रंट अनंत सूक्ष्म बना रहता है।
• घनत्व अनुपात संवेदनशीलता:
  • समान और निम्न-घनत्व अनुपात ($\rho_L/\rho_S \approx 1$ से $2$) के लिए, मषी समाधान$\Delta T$ के विस्तृत रेंज में शार्प-इंटरफेस समाधान के करीब होता है। फलस्वरूप, समान-घनत्व बेंचमार्क घनत्व-परिवर्तन-प्रेरित प्रवाह के लिए विधियों को सत्यापित करने हेतु पर्याप्त कठोर नहीं हैं।
  • बड़े घनत्व अनुपात ($\rho_L/\rho_S \gg 1$) के लिए, चरण-फ्रंट का स्थान $\Delta T$ के प्रति अत्यधिक संवेदनशील हो जाता है। मषी समाधान केवल तभी शास्त्रीय शार्प-इंटरफेस सीमा का अनुमान लगाता है जब $\Delta T$ पर्याप्त छोटा हो।
• विधि का प्रदर्शन:
  • उन्नत लो-मैच एन्थैल्थी विधि $\Delta T$ के विस्तृत रेंज और उच्च घनत्व अनुपातों (सत्यापन मामलों में $\rho_L/\rho_S \approx 540$ तक) के लिए व्युत्पन्न मषी स्टेफ़न समानता समाधानों के साथ उत्कृष्ट तालमेल प्रदर्शित करती है।
  • इसके विपरीत, मूल प्रतिरूपण बड़े घनत्व अनुपातों के लिए इंटरफेशियल अस्थिरता और तरल क्षेत्र में अस्वाभाविक भंवर संरचनाएं (vortical structures) उत्पन्न करता है।
  • अभिसरण दरें: उन्नत विधि चरण-फ्रंट स्थान और तापमान क्षेत्रों दोनों के लिए प्रथम और द्वितीय क्रम के बीच स्थानिक-कालिक अभिसरण दर प्राप्त करती है, यहाँ तक कि $O(3)$ के घनत्व अनुपात पर भी।

महत्व और दावे
यह शोध पत्र परिमित $\Delta T$ पर संख्यात्मक विधियों के सत्यापन के लिए एक कठोर सैद्धांतिक बेंचमार्क (गैर-समतापीय परिवर्तनशील-घनत्व मषी स्टेफ़न समस्या का समानता समाधान) प्रदान करके एन्थैल्थी विधि साहित्य में एक लंबे समय से चली आ रही कमी को भरने का दावा करता है।

लेखक तर्क देते हैं कि:

1. बेंचमार्क वैधता: बड़े-घनत्व-अनुपात वाले मामले आधुनिक एन्थैल्थी विधियों के लिए समान-घनत्व मामलों की तुलना में काफी अधिक कठोर सत्यापन बेंचमार्क प्रदान करते हैं, क्योंकि बाद वाले चरण-परिवर्तन अंतराल के प्रति समाधान की संवेदनशीलता को पकड़ने में विफल रहते हैं।
2. पद्धतिगत सुधार: प्रस्तावित उन्नत लो-मैच एन्थैल्थी विधि (विशेष रूप से संशोधित एन्थैल्थी परिभाषा और रैखिक $\phi$-$T$ संबंध) पर्याप्त घनत्व अंतर के साथ चरण परिवर्तन के अनुकरण के लिए स्थिरता और सटीकता के लिए आवश्यक है।
3. विस्तारित प्रयोज्यता: ये सैद्धांतिक और संख्यात्मक विकास एन्थैल्थी विधियों की प्रयोज्यता को पारंपरिक धातु जमने (जहाँ घनत्व अनुपात लगभग 1 होता है) से आगे बढ़ाकर उबालने और संघनन जैसे बहुत बड़े घनत्व अनुपात वाले प्रक्रियाओं तक विस्तारित करते हैं।

यह कार्य नए प्रायोगिक अनुप्रयोगों का प्रस्ताव नहीं करता है, बल्कि व्यापक श्रेणी की परिवर्तनशील-घनत्व वाली चरण-परिवर्तन समस्याओं के लिए एन्थैल्थी-आधारित सिमुलेशन को विस्तारित करने हेतु एक मजबूत संख्यात्मक ढांचा और सैद्धांतिक आधार स्थापित करता है।