A Temperature-Coupled Cahn-Hilliard-Stokes-Heat Model for Thermally Driven Phase Separation

यहाँ शोध पत्र "A Temperature-Coupled Cahn–Hilliard–Stokes–Heat Model for Thermally Driven Phase Separation" का विस्तृत तकनीकी सारांश दिया गया है।

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

यह शोध पत्र श्यान (viscous), असंपीड्य (incompressible) द्विआधारी मिश्रणों में तापीय रूप से प्रेरित चरण पृथक्करण (thermally driven phase separation) के गणितीय मॉडलिंग और संख्यात्मक सिमुलेशन को संबोधित करता है, जो विशेष रूप से इलेक्ट्रिक वाहन बैटरी हाउसिंग के भीतर जल संघनन (water condensation) से प्रेरित है।

• भौतिक संदर्भ: बैटरी एनक्लोजर में प्रवेश करने वाली नमी ठंडी आंतरिक सतहों पर संघनित हो सकती है, जिससे क्षरण (corrosion) और इन्सुलेशन में गिरावट आती है। विस्तृत ज्यामितीय और परिचालन श्रेणियों में प्रयोगात्मक परीक्षण महंगा है, जिसके लिए मजबूत संख्यात्मक उपकरणों की आवश्यकता है।
• अंतराल (The Gap): मौजूदा मॉडल अक्सर समतापी (isothermal) सेटिंग्स या सामान्य उबलने/संघनन कॉन्फ़िगरेशन पर ध्यान केंद्रित करते हैं। एक ऐसे ढांचे का अभाव है जो सीमित ज्यामिति (जैसे बैटरी केसिंग) के लिए प्रासंगिक यथार्थवादी, मिश्रित सीमा स्थितियों (mixed boundary conditions) के तहत परिवर्तनशील-घनत्व वाले स्टोक्स प्रवाह (variable-density Stokes flow) और ऊष्मा परिवहन के साथ चरण-क्षेत्र गतिशीलता (phase-field dynamics) को युग्मित (couple) कर सके।
• मुख्य चुनौती: इस प्रणाली में तीन अलग-अलग भौतिक तंत्रों के बीच गहरा युग्मन शामिल है:
  1. चरण पृथक्करण (Phase Separation): Cahn–Hilliard समीकरण (CHE) द्वारा नियंत्रित।
  2. द्रव प्रवाह (Fluid Flow): परिवर्तनशील घनत्व और श्यानता के साथ स्टोक्स समीकरणों द्वारा नियंत्रित।
  3. ऊष्मा परिवहन (Heat Transport): संवहन (convection) के साथ ऊष्मा समीकरण द्वारा नियंत्रित।
     यह युग्मन द्विदिश (bidirectional) है: तापमान मुक्त ऊर्जा परिदृश्य (free energy landscape) को प्रभावित करता है (जो चरण परिवर्तन को संचालित करता है), जबकि चरण क्षेत्र (सांद्रता) द्रव के घनत्व और श्यानता को प्रभावित करता है (जो प्रवाह को संचालित करता है)।

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

A. गणितीय मॉडलिंग

लेखक एक डिफ्यूज-इंटरफेस मॉडल (diffuse-interface model) प्रस्तावित करते हैं जिसमें एक परिबद्ध डोमेन $\Omega \subset \mathbb{R}^d$ ($d=2,3$) में तीन युग्मित आंशिक अंतर समीकरण (PDEs) शामिल हैं:

1. संवहनी Cahn–Hilliard समीकरण (Convective Cahn–Hilliard Equation):
   यह ऑर्डर पैरामीटर $c$ (जो आयतन अंश/सांद्रता का प्रतिनिधित्व करता है) के विकास का वर्णन करता है।
   $$\frac{\partial c}{\partial t} + (u \cdot \nabla)c = \frac{1}{Pe} \Delta \mu$$
   जहाँ $\mu$ एक तापमान-निर्भर लैंडौ मुक्त ऊर्जा (temperature-dependent Landau free energy) से प्राप्त रासायनिक क्षमता है।

2. स्टोक्स समीकरण (Stokes Equations - Quasi-static):
   यह परिवर्तनशील गुणांकों के साथ क्रीपिंग फ्लो ($Re \ll 1$) के लिए वेग $u$ और दबाव $p$ का वर्णन करता है।
   $$-\nabla p + \nabla \cdot (2\eta(c) E(u)) + G \rho(c) \hat{g} = 0, \quad \nabla \cdot u = 0$$
   यहाँ, घनत्व $\rho(c)$ और श्यानता $\eta(c)$ $c$ के अफ़ाइन फलन (affine functions) हैं। एक उछाल शब्द (buoyancy term) $G \rho(c) \hat{g}$ गुरुत्वाकर्षण को दर्शाता है।

3. ऊष्मा समीकरण (Heat Equation):
   यह तापमान क्षेत्र $\Theta$ का वर्णन करता है।
   $$\frac{\partial \Theta}{\partial t} + u \cdot \nabla \Theta = \frac{1}{Pe_\Theta} \Delta \Theta$$

मुक्त ऊर्जा में मुख्य नवाचार (Key Innovation in Free Energy):
बल्क फ्री एनर्जी डेंसिटी $f(c, \Theta)$ को तापमान-निर्भर बनाने के लिए संशोधित किया गया है:
$$f(c, \Theta) = 2c^4 - (1 - \beta(\Theta))c^2$$
जहाँ $\beta(\Theta) = \Theta / \Theta_S$ स्थानीय तापमान और संतृप्ति तापमान का अनुपात है।

• यदि $\Theta > \Theta_S$ ($\beta > 1$): क्षमता एक सिंगल-वेल (single-well) है, जो एक सजातीय चरण (homogeneous phase) के पक्ष में है।
• यदि $\Theta < \Theta_S$ ($\beta < 1$): क्षमता एक डबल-वेल (double-well) बन जाती है, जो चरण पृथक्करण को प्रेरित करती है।

B. सीमा स्थितियाँ (Boundary Conditions)

मॉडल यथार्थवादी बैटरी ज्यामिति को प्रतिबिंबित करने के लिए मिश्रित सीमा स्थितियों (mixed boundary conditions) का उपयोग करता है:

• वेग (Velocity): कुछ सीमाओं ($\Gamma_D$) पर नो-स्लिप (no-slip) और अन्य पर ट्रैक्शन-फ्री (traction-free)।
• तापमान (Temperature): कुछ सीमाओं पर डिरिचलेट (Dirichlet - स्थिर तापमान) और अन्य पर न्यूमैन (Neumann - इंसुलेटेड)।
• चरण क्षेत्र (Phase Field): $c$ और $\mu$ दोनों के लिए होमोजेनियस न्यूमैन (no-flux)।

C. संख्यात्मक योजना (Numerical Scheme)

एक पूर्णतः विविक्त परिमित तत्व योजना (fully discrete finite element scheme) प्रस्तावित है जो एक अनुक्रमिक विभाजन रणनीति (sequential splitting strategy) का उपयोग करती है:

1. समय विवक्तीकरण (Time Discretization): प्रथम-क्रम सेमी-इम्प्लिसिट (semi-implicit)।
   • ऊष्मा समीकरण: बैकवर्ड युलर (Backward Euler)।
   • स्टोक्स: अपडेटेड गुणांकों के साथ वेग/दबाव के लिए इमplicit समाधान।
   • Cahn–Hilliard: कॉन्वेक्स-स्प्लिटिंग (convex-splitting) (Eyre-type) विवक्तीकरण। डिफ्यूसिव स्टेप के लिए अनकंडीशनल एनर्जी स्थिरता सुनिश्चित करने के लिए मुक्त ऊर्जा को उत्तल (convex - इमप्लिसिट रूप से उपचारित) और अवतल (concave - एक्सप्लिसिट रूप से उपचारित) भागों में विभाजित किया गया है।
2. स्थानिक विवक्तीकरण (Spatial Discretization):
   • वेग/दबाव: टेलर-हुड तत्व ($P_2/P_1$)।
   • स्केलर ($c, \mu, \Theta$): लीनियर तत्व ($P_1$)।
3. सॉल्वर (Solver): Cahn–Hilliard स्टेप के लिए नॉन-लीनियर बीजगणितीय प्रणालियों को डैम्प्ड न्यूटन इटरेशन (damped Newton iteration) के माध्यम से हल किया जाता है।

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

1. ऊष्मागतिक रूप से प्रेरित क्षमता (Thermodynamically Motivated Potential): एक लैंडौ-टाइप क्षमता का परिचय, जहाँ चरण व्यवहार (सजातीय बनाम पृथक) एड-हॉक स्रोत पदों के बजाय सीधे स्थानीय तापमान अनुपात $\beta(\Theta)$ द्वारा नियंत्रित होता है।
2. कठोर अस्तित्व सिद्धांत (Rigorous Existence Theory):
   • 2D और 3D में युग्मित प्रणाली के लिए कमजोर समाधानों (weak solutions) के स्थानीय-इन-टाइम अस्तित्व को सिद्ध किया।
   • ऑर्डर पैरामीटर पर वैश्विक $L^\infty$ सीमाओं की कमी को संबोधित करने के लिए कार्यात्मक ढांचे $H^2 \to L^2 \to H^{-2}$ का उपयोग किया।
   • एक कम्पोजिट ऑपरेटर $T(u) = S(C(u, H(u)))$ पर शौडर फिक्स्ड-पॉइंट तर्क (Schauder fixed-point argument) का उपयोग किया, जिसमें हीट, Cahn–Hilliard और स्टोक्स ऑपरेटरों की निरंतरता और कॉम्पैक्टनेस का अलग-अलग विश्लेषण किया गया।
   • मिश्रित सीमा स्थितियों को संभाला जो मानक ऊर्जा अनुमानों को जटिल बनाती हैं।
3. संख्यात्मक स्थिरता (Numerical Stability): एक ऐसी योजना विकसित की जो द्रव्यमान का संरक्षण करती है (अभेद्य सीमाओं के तहत) और आइसोथर्मल डिफ्यूसिव लिमिट में Cahn–Hilliard सबस्टेप के लिए अनकंडीशनली एनर्जी-स्टेबल है।
4. बैटरी सुरक्षा के लिए अनुप्रयोग (Application to Battery Safety): वास्तविक बैटरी-प्रासंगिक ज्यामिति में सीमित स्थानों के भीतर संघनन के सिमुलेशन के लिए एक विशिष्ट ढांचा प्रदान किया।

4. परिणाम (Results)

लेखकों ने चार संख्यात्मक प्रयोगों के माध्यम से अपने मॉडल को सत्यापित किया:

1. बेंचमार्क (Lid-Driven Cavity):

   • टू-फेज स्टोक्स फ्लो के लिए मौजूदा साहित्य (Khatavkar et al.) के विरुद्ध तुलना की।
   • परिणाम: योजना ने वेग प्रोफाइल और इंटरफ़ेस विरूपण रुझानों को सटीक रूप से पुनरुत्पादित किया, जिससे हाइड्रोडायनामिक विवक्तीकरण की पुष्टि हुई।

2. समतापी चरण पृथक्करण (Isothermal Phase Separation):

   • तापमान के प्रति क्षमता की प्रतिक्रिया का परीक्षण किया।
   • परिणाम:
     • $\Theta > \Theta_S$ पर: सिस्टम एक सजातीय अवस्था की ओर वापस आया (सिंगल-वेल व्यवहार)।
     • $\Theta < \Theta_S$ पर: सिस्टम ने तीव्र स्पिनोडल डिकंपोजिशन और कोआर्सेनिंग (coarsening) प्रदर्शित की (डबल-वेल व्यवहार)।

3. गुरुत्वाकर्षण-प्रेरित चरण पृथक्करण (Gravity-Driven Phase Separation):

   • परिवर्तनशील घनत्व ($\lambda_\rho \approx 0.001$) और श्यानता के साथ उछाल (buoyancy) को सक्रिय किया।
   • परिणाम: सघन तरल डोमेन बने और गुरुत्वाकर्षण के कारण नीचे की ओर बढ़े, जो नीचे की सीमा पर जमा हुए और एक भारी परत में मिल गए, जो चरण पृथक्करण और उछाल के बीच संबंध को दर्शाता है।

4. तापीय रूप से प्रेरित चरण परिवर्तन (Heated/Cooled Boundaries):

   • एक डोमेन में थर्मल ग्रेडिएंट का सिमुलेशन किया।
   • गर्म सीमा (Heated Boundary): स्थानीय हीटिंग ने चरण पृथक्करण को दबा दिया, जिससे गर्म दीवार के पास के तरल डोमेन एक सजातीय वाष्प चरण में घुल गए।
   • ठंडी सीमा (Cooled Boundary): स्थानीय शीतलन ने ठंडी दीवार के पास तरल बूंदों के न्यूक्लिएशन (nucleation) को ट्रिगर किया, जो बाद में बढ़ीं और मोटे (coarsen) हुईं।
   • महत्व: इसने सफलतापूर्वक बैटरी हाउसिंग के भीतर ठंडी सतहों पर संघनन की भौतिक प्रक्रिया की नकल की।

5. महत्व (Significance)

• सैद्धांतिक प्रगति: यह शोध पत्र तापमान, परिवर्तनशील घनत्व/श्यानता और मिश्रित सीमा स्थितियों को संभालने के माध्यम से डिफ्यूज-इंटरफेस मॉडल के गणितीय विश्लेषण में एक अंतराल को भरता है, जिन्हें अक्सर सरल आइसोथर्मल मॉडलों में छोड़ दिया जाता है।
• इंजीनियरिंग उपयोगिता: यह इलेक्ट्रिक वाहन बैटरी पैक में संघनन पैटर्न (condensation patterns) की भविष्यवाणी करने के लिए एक मजबूत, गणितीय रूप से आधारित उपकरण प्रदान करता है। यह बैटरी के जीवन और सुरक्षा को बढ़ाने के लिए थर्मल मैनेजमेंट सिस्टम को डिजाइन करने हेतु महत्वपूर्ण है।
• कंप्यूटेशनल मजबूती: प्रस्तावित कॉन्वेक्स-स्प्लिटिंग योजना जटिल चरण संक्रमणों के स्थिर सिमुलेशन की अनुमति देती है, जो अक्सर एक्सप्लिसिट विधियों या सिंगुलर पोटेंशियल की कठोरता (stiffness) से जुड़ी समय-चरण प्रतिबंधों से मुक्त होती है।

संक्षेप में, यह कार्य तापीय रूप से प्रेरित चरण पृथक्करण को मॉडल करने के लिए थर्मोडायनामिक्स, फ्लूइड मैकेनिक्स और हीट ट्रांसफर को एकीकृत करने वाला एक व्यापक ढांचा प्रस्तुत करता है, जो कठोर गणितीय प्रमाण और वास्तविक दुनिया के इंजीनियरिंग अनुप्रयोगों के लिए व्यावहारिक संख्यात्मक सिमुलेशन दोनों प्रदान करता है।