Quantitative Homogenization of a Cahn--Hilliard System with Source Term in Periodically Perforated Domains

तकनीकी सारांश: आवधिक छिद्रित डोमेन में स्रोत पद (Source Term) के साथ एक काह्न-हिलियर्ड (Cahn–Hilliard) प्रणाली का मात्रात्मक समरूपीकरण (Quantitative Homogenization)

समस्या विवरण
यह शोध पत्र एक आवधिक छिद्रित डोमेन में एक गैर-संरक्षण (nonconservative), एकदिष्ट (monotone) स्रोत पद वाली काह्न-हिलियर्ड प्रणाली के गुणात्मक और मात्रात्मक समरूपीकरण की जांच करता है। सूक्ष्म मॉडल एक सरंध्र माध्यम (porous medium) में चरण पृथक्करण (phase separation) का वर्णन करता है जहाँ तरल एक जुड़े हुए पोर स्पेस $\Omega^\varepsilon_p \subset \Omega \subset \mathbb{R}^d$ ($d \in \{2,3\}$) में स्थित है, जिसे एक सीमित लिप्सचिट्ज़ डोमेन (Lipschitz domain) से ठोस समावेशन (solid inclusions) की एक आवधिक व्यवस्था को हटाकर प्राप्त किया गया है। यह प्रणाली ऑर्डर पैरामीटर $\phi^\varepsilon$ और रासायनिक क्षमता (chemical potential) $\mu^\varepsilon$ को निम्नलिखित के माध्यम से जोड़ती है:
$$\begin{cases} \partial_t \phi^\varepsilon - \Delta \mu^\varepsilon + G(\phi^\varepsilon) = 0, \\ \mu^\varepsilon = -\Delta \phi^\varepsilon + F'(\phi^\varepsilon), \end{cases}$$
जो बाहरी सीमा $\partial \Omega$ और पोर दीवारों $\Gamma^\varepsilon_s$ दोनों पर सजातीय न्यूमैन सीमा स्थितियों (homogeneous Neumann boundary conditions) के अधीन है। विभव (Potential) $F$ बहुपद वृद्धि (polynomial growth) और कोएर्सिविटी (coercivity) (जैसे कि शास्त्रीय डबल-वेल पोटेंशियल) का पालन करता है, जबकि स्रोत पद $G$ वैश्विक रूप से लिप्सचिट्ज़, एकदिष्ट है और $G(0)=0$ को संतुष्ट करता है। मुख्य चुनौती यह है कि जैसे-जैसे $\varepsilon \to 0$ होता है, प्रभावी स्थूल (macroscopic) समीकरणों को व्युत्पन्न करना और सूक्ष्म तथा समरूपीकृत समाधानों के बीच त्रुटि के लिए तीक्ष्ण अभिसरण दर (sharp convergence rates) स्थापित करना है।

कार्यप्रणाली (Methodology)
विश्लेषण आवधिक छिद्रित डोमेन के लिए अनुकूलित आवधिक अनफोल्डिंग विधि (periodic unfolding method) का उपयोग करता है, जिसमें दोलनशील ज्यामिति (oscillating geometry) को संभालने के लिए अनफोल्डिंग ऑपरेटर $T^\varepsilon_*$ का उपयोग किया जाता है। कार्यप्रणाली दो मुख्य चरणों में आगे बढ़ती है:

1. गुणात्मक समरूपीकरण (Qualitative Homogenization):

   • ऊर्जा अनुमान (Energy Estimates): सूक्ष्म समाधानों के लिए समान ए प्रियोरी अनुमान (uniform a priori estimates) प्राप्त किए जाते हैं। एक प्रमुख संरचनात्मक अवलोकन यह है कि एकदिष्ट स्रोत पद $G$ (जहाँ $G(0)=0$) ऊर्जा पहचान (energy identity) में एक विसरक (dissipative) पद में योगदान देता है। विशेष रूप से, रासायनिक क्षमता समीकरण को $G(\phi^\varepsilon)$ के साथ परीक्षण करने से एक धनात्मक ग्रेडिएंट पद प्राप्त होता है, जो प्रतिस्पर्धा करने के बजाय अपव्यय (dissipation) को सुदृढ़ करता है। यह बिना किसी एक्सपोनेंशियल टाइम फैक्टर्स वाले ग्रोनवालल तर्क (Gronwall argument) की आवश्यकता के, $\varepsilon$ से स्वतंत्र समान बाउंड्स की अनुमति देता है।
   • टू-स्केल अभिसरण (Two-Scale Convergence): अनफोल्डिंग ऑपरेटर और कॉम्पैक्टनेस तर्क (Aubin–Lions–Simon lemma) का उपयोग करते हुए, लेखक $(\phi^\varepsilon, \mu^\varepsilon)$ के एक प्रभावी युग्म $(\phi, \mu)$ की ओर अभिसरण को सिद्ध करते हैं। सीमांत प्रणाली एक प्रभावी डिफ्यूजन टेंसर $B_{\text{hom}}$ के साथ मैक्रोस्कोपिक डोमेन $\Omega$ पर एक काह्न-हीलिएड समीकरण है।
   • सेल समस्याएँ (Cell Problems): प्रभावी टेंसर $B_{\text{hom}}$ को संदर्भ पोर सेल $Y_p$ पर स्केलर न्यूमैन सेल समस्याओं द्वारा अभिलक्षित किया जाता है। सुधारक (correctors) $\chi_i$, $H^1_{\text{per},0}(Y_p)$ में स्थित हैं।

2. मात्रात्मक विश्लेषण (संशोधक अनुमान) (Quantitative Analysis - Corrector Estimates):

   • स्केल-स्प्लिटिंग ऑपरेटर (Scale-Splitting Operator): मात्रात्मक दर प्राप्त करने के लिए, प्रथम-क्रम सुधारक सन्निकटन $\Phi^\varepsilon$ और $M^\varepsilon$ का निर्माण $Q_1$ स्केल-स्प्लिटिंग ऑपरेटर $Q_\varepsilon$ (Griso [29, 30] से) का उपयोग करके किया जाता है। यह ऑपरेटर एक निश्चित मैक्रोस्कोपिक डोमेन $\Omega$ पर कार्य करता है और तत्पश्चात इसे $\Omega^\varepsilon_p$ पर प्रतिबंधित किया जाता है। महत्वपूर्ण रूप से, यह दृष्टिकोण यह सुनिश्चित करता है कि सेल सुधारक $H^1_{\text{per}}(Y_p)$ में रहें, जिससे शास्त्रीय औपचारिक विस्तार (formal expansions) के लिए आवश्यक प्रतिबंधात्मक $W^{1,\infty}$ धारणाओं से बचा जा सके।
   • वैरिएशनल रेसिड्यूल विधि (Variational Residual Method): त्रुटि विश्लेषण एक नेगेटिव-नॉर्म फ्रेमवर्क में किया जाता है। ऑर्डर पैरामीटर त्रुटि $e^\varepsilon_\phi = \phi^\varepsilon - \Phi^\varepsilon$ को न्यूमैन इनवर्स लैप्लासियन $N_\varepsilon e^\varepsilon_\phi$ के विरुद्ध परीक्षण किया जाता है, जबकि रासायनिक क्षमता त्रुटि को $e^\varepsilon_\phi$ के विरुद्ध परीक्षण किया जाता है। यह अवशेषों (residuals) की उच्च नियमितता (regularity) की आवश्यकता को समाप्त करता है।
   • संगति अनुमान (Consistency Estimates): प्रमाण вариेशनल रेसिड्यूल अनुमानों (लेम्मा 4.6) पर निर्भर करता है जो $\partial \Omega$ के पास अपूर्ण कोशिकाओं (incomplete cells) के बाउंड्री लेयर को ध्यान में रखते हैं।

प्रमुख परिणाम

• गुणात्मक अभिसरण (प्रमेय 3.7): सूक्ष्म समाधान (उप-अनुक्रमों तक) समरूपीकृत प्रणाली की ओर अभिसरित होते हैं:
  $$\begin{cases} \partial_t \phi - \text{div}(B_{\text{hom}}\nabla \mu) + G(\phi) = 0, \\ \mu = -\text{div}(B_{\text{hom}}\nabla \phi) + F'(\phi), \end{cases}$$
  जहाँ प्रभावी न्यूमैन सीमा स्थितियाँ लागू हैं। टेंसर $B_{\text{hom}}$ सममित (symmetric) और यूनिफॉर्मली एलिप्टिक है।

• मात्रात्मक दरें (प्रमेय 4.7): एक अतिरिक्त $H^2$-नियमितता धारणा और वेल-प्रिपेयर्ड प्रारंभिक डेटा के तहत, यह शोध पत्र $\varepsilon^{1/2}$ की क्रम की त्रुटि अनुमान स्थापित करता है:
  $$\|\phi^\varepsilon - \phi\|_{L^2(0,T; L^2(\Omega^\varepsilon_p))} + \|\nabla \phi^\varepsilon - (\nabla \phi + \sum Q_\varepsilon(\partial_{x_i}\phi)\nabla_y \chi_i(\tfrac{x}{\varepsilon}))\|_{L^2((0,T)\times\Omega^\varepsilon_p)} + \|\mu^\varepsilon - \mu\|_{L^2(0,T; H^1(\Omega^\varepsilon_p)')} \leq C\varepsilon^{1/2}.$$

  • अन-करेक्टेड ऑर्डर पैरामीटर त्रुटि को $L^2(0,T; L^2(\Omega^\varepsilon_p))$ में नियंत्रित किया जाता है।
  • करेक्टेड ग्रेडिएंट त्रुटि को $L^2((0,T)\times\Omega^\varepsilon_p)$ में नियंत्रित किया जाता है।
  • रासायनिक क्षमता त्रुटि को ड्यूल नॉर्म $L^2(0,T; H^1(\Omega^\varepsilon_p)')$ में नियंत्रित किया जाता है।

• बाउंड्री लेयर प्रभाव (Boundary Layer Effects): $\partial \Omega$ के पास अपूर्ण कोशिकाओं के बाउंड्री लेयर के कारण सीमित डोमेन पर $\varepsilon^{1/2}$ की दर प्राप्त होती है। फ्लैट टॉरस $\mathbb{T}^d$ पर, जहाँ ऐसा कोई बाउंड्री लेयर नहीं होता, दर सुधरकर $O(\varepsilon)$ हो जाती है (टिप्पणी 4.10)।

महत्व और योगदान
यह शोध पत्र छिद्रित मीडिया में फेज-फील्ड समीकरणों के होमोजेनाइजेशन सिद्धांत में दो प्राथमिक योगदान देता है:

1. गुणात्मक विस्तार (Qualitative Extension): यह एक सामान्य एकदिष्ट गैर-संरक्षण स्रोत पद के साथ काह्न-हिलियर्ड प्रणाली को कठोरता से व्युत्पन्न करता है, यह प्रदर्शित करते हुए कि स्रोत पद की विसरक संरचना सामान्य गैर-संरक्षण बल की तुलना में ऊर्जा अनुमानों को सरल बनाती है।
2. मात्रात्मक सुधार (Quantitative Improvement): यह छिद्रित डोमेन में काह्न-हिलियर्ड प्रणालियों के लिए ज्ञात अभिसरण दर को $\varepsilon^{1/4}$ (शमकर-कलियाडासिस [49] द्वारा स्थापित) से बढ़ाकर छिद्रित डोमेन में द्वितीय-क्रम एलिप्टिक समस्याओं के लिए स्वाभाविक दर $\varepsilon^{1/2}$ तक सुधारता है।
   • मेथोडोलॉजिकल प्रगति: $W^{1,\infty}$ में सेल सुधारकों और $C^1(0,T; W^{k,\infty})$ ($k \geq 4$) में मैक्रोस्कोपिक समाधानों की आवश्यकता वाले पिछले कार्यों के विपरीत, यह दृष्टिकोण $Q_\varepsilon$ स्केल-स्प्लिटिंग ऑपरेटर का उपयोग करता है। यह मैक्रोस्कोपिक समाधान की आवश्यक नियमितता को $H^2$ तक कम करता है और सेल सुधारकों को $H^1_{\text{per}}$ में रखता है।
   • वैरिएशनल संगति: संगति अनुमान पूरी तरह से वैरिएशनल रूप में व्युत्पन्न किए गए हैं, जो सेल समस्या से विरासत में मिले कमजोर नो-फ्लक्स (no-flux) की स्थिति का उपयोग करते हैं। यह उन बाउंड्री-टर्म योगदानों से बचता है जो पहले $\varepsilon^{1/4}$ की दर गिरावट का कारण बनते थे।

लेखक नोट करते हैं कि वर्तमान वैरिएशनल ढांचे के भीतर $\varepsilon^{1/2}$ की दर बाउंड्री लेयर के कारण शार्प है, और सीमित डोमेन पर $O(\varepsilon)$ की इष्टतम दर प्राप्त करने के लिए संभवतः एक द्वैत तर्क (Duality argument) (जैसा कि केनिग-लिन-शेन [40] में है) या एक ऑपरेटर-अनुमान दृष्टिकोण (बिरमन-सुस्लिना [7–9]) की आवश्यकता होगी, जो भविष्य के कार्यों के लिए छोड़ दिए गए हैं। यह शोध पत्र यह दावा नहीं करता है कि यह इस क्रम पर रासायनिक क्षमता $\mu^\varepsilon$ के लिए ग्रेडिएंट दर प्राप्त करता है, क्योंकि अवशेषों की वर्तमान नियमितता आवश्यक परीक्षण को न्यायसंगत नहीं बनाती है।