Singularities in phase separation models: a spectral element approach for the nonlocal Cahn-Hilliard equation

यहाँ शोध पत्र "Singularities in phase separation models: a spectral element approach for the nonlocal Cahn–Hilliard equation" का विस्तृत तकनीकी सारांश दिया गया है।

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

यह शोध पत्र नॉनलोकल (nonlocal) काह्न–हिलियर्ड (Cahn–Hilliard - CH) समीकरण के संख्यात्मक सन्निकटन (numerical approximation) को संबोधित करता है, जो द्विआधारी मिश्रणों (जैसे मिश्र धातु) में चरण पृथक्करण (phase separation) का वर्णन करने वाला एक स्थूल मॉडल है। शास्त्रीय स्थानीय CH समीकरण के विपरीत, नॉनलोकल संस्करण एक सिंगुलर कर्नेल (singular kernel) के साथ एक कन्वोल्यूशन टर्म (convolution term) के माध्यम से दीर्घ-परासी अंतःक्रियाओं (long-range interactions) को शामिल करता है।

अध्ययन किया गया विशिष्ट तंत्र है:
$$\frac{\partial \rho}{\partial t} = \nabla \cdot \left( m(\rho) \nabla \mu \right), \quad \mu = F'(\rho) - K * \rho$$
जहाँ:

• $\rho$: फेज़ फील्ड वेरिएबल (सांद्रता/concentration)।
• $m(\rho)$: स्थिर गतिशीलता (constant mobility) (इसे $m=1$ पर सेट किया गया है)।
• $F(\rho)$: एक लॉगैरिद्मिक पोटेंशियल (logarithmic potential) (शुद्ध अवस्थाओं $\pm 1$ पर सिंगुलर), जिसे $F(s) = \frac{\theta}{2}[(1+s)\log(1+s) + (1-s)\log(1-s)]$ के रूप में परिभाषित किया गया है। यह सुनिश्चित करता है कि समाधान भौतिक अंतराल $(-1, 1)$ के भीतर रहे।
• $K$: एक न्यूटनियन इंटरेक्शन कर्नेल (Newtonian interaction kernel) (2D में, $K(x) = \frac{1}{2\pi}\log\|x\|$), जो मूल बिंदु (origin) पर सिंगुलर है लेकिन समाकलनीय (integrable) है।
• सीमा शर्तें (Boundary Conditions): एक परिबद्ध डोमेन $\Omega$ पर नो-फ्लक्स ($\partial_n \mu = 0$)।

मुख्य चुनौतियाँ:

1. सिंगुलैरिटी (Singularity): कर्नेल $K$ में मूल बिंदु पर एक सिंगुलैरिटी है, जिससे बिना विशेष उपचार के मानक कोलोकेशन (collocation) या परिमित अंतर (finite difference) विधियों को सीधे लागू करना कठिन हो जाता है।
2. नॉन-लोकैलिटी (Non-locality): कन्वोल्यूशन टर्म $K * \rho$ घने ऑपरेटरों (dense operators) का परिणाम देता है, जिससे मानक विविक्तकरण (discretization) के लिए उच्च कम्प्यूटेशनल लागत ($O(N^2)$) उत्पन्न होती है।
3. सिंगुलर पोटेंशियल (Singular Potential): जैसे-जैसे $\rho \to \pm 1$ होता है, $F'(\rho)$ अनियंत्रित (blow up) हो जाता है, जिसके लिए समाधान को भौतिक डोमेन से बाहर जाने से रोकने हेतु सख्त संख्यात्मक नियंत्रण की आवश्यकता होती है।

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

लेखक सिंगुलर कन्वोल्यूशन को कुशलतापूर्वक संभालने के लिए एक प्सूडोस्पेक्ट्रल मल्टीशेप विधि (pseudospectral multishape method) को कटऑफ डिकंपोजिशन (cutoff decomposition) तकनीक के साथ प्रस्तावित करते हैं।

A. गणितीय ढांचा और अस्तित्व सिद्धांत (Mathematical Framework & Existence Theory)

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

• डोमेन सामान्यीकरण (Domain Generalization): वे केवल चिकने (smooth) डोमेन के बजाय लिप्सचिट्ज़ डोमेन (कोणीय बहुभुज डोमेन सहित) पर कमजोर समाधानों (weak solutions) के अस्तित्व और अद्वितीयता को सिद्ध करते हैं।
• सख्त पृथक्करण (Strict Separation): वे यह सिद्ध करने के लिए "सख्त पृथक्करण गुण" का उपयोग करते हैं कि किसी भी $t > 0$ के लिए समाधान शुद्ध अवस्थाओं ($\pm 1$) से समान रूप से दूर रहता है, यहाँ तक कि 3D में भी।
• दीर्घकालिक व्यवहार (Long-time Behavior): वे Łojasiewicz–Simon inequality का उपयोग करके एक अद्वितीय संतुलन अवस्था की ओर अभिसरण दरों को स्थापित करते हैं, जो $L^2$ नॉर्म में बहुपद क्षय (polynomial decay) को दर्शाता है।

B. सिंगुलर कन्वोल्यूशन का संख्यात्मक सन्निकटन

$(K * \rho)(x)$ की गणना करने के लिए, इंटीग्रल को एक कटऑफ त्रिज्या $\varepsilon$ का उपयोग करके दो भागों में विभाजित किया जाता है:
$$(K * \rho)(x) = \underbrace{\int_{\Omega \setminus B_\varepsilon(x)} K(x-y)\rho(y) dy}_{I_1} + \underbrace{\int_{\Omega \cap B_\varepsilon(x)} K(x-y)\rho(y) dy}_{I_2}$$

1. सिंगुलैरिटी को संभालना ($I_2$):

   • छोटे बॉल $B_\varepsilon(x)$ पर आंतरिक इंटीग्रल का अनुमान एक शून्य-क्रम टेलर विस्तार (zero-order Taylor expansion) (piecewise constant approximation) का उपयोग करके लगाया जाता है:
     $$I_2 \approx \rho(x) \int_{x - (\Omega \cap B_\varepsilon(x))} K(u) du$$
   • लेखक इन ज्यामितीय क्षेत्रों पर न्यूटनियन कर्नेल के इंटीग्रल के लिए क्लोज्ड-फॉर्म एनालिटिकल एक्सप्रेशन (closed-form analytical expressions) व्युत्पन्न करते हैं (ध्रुवीय निर्देशांकों और विशिष्ट ज्यामितीय विखंडन का उपयोग करके)। यह पद एक सुधार कारक (correction factor) और त्रुटि संकेतक के रूप में कार्य करता है।
   • इस सन्निकटन की त्रुटि $O(\varepsilon^2 \log \varepsilon)$ दिखाई गई है।

2. नियमित भाग को संभालना ($I_1$):

   • सिंगुलैरिटी को छोड़कर डोमेन पर बाहरी इंटीग्रल को स्पेक्ट्रल एलिमेंट मेथड (SEM) का उपयोग करके कंप्यूट किया जाता है।
   • वे MultiShape सॉफ्टवेयर पैकेज का उपयोग करते हैं, जो चतुर्भुज तत्वों (quadrilateral elements) को जोड़कर जटिल डोमेन का निर्माण करता है।
   • प्रत्येक कोलोकेशन बिंदु $x_i$ के लिए, डोमेन $\Omega \setminus B_\varepsilon(x_i)$ को चतुर्भुजों के "मल्टीशेप" में विभाजित किया जाता है।
   • कन्वोल्यूशन को $\rho$ के विविक्तकृत (discretized) रूप पर कार्य करने वाले एक मैट्रिक्स $M$ के रूप में प्रीकंप्यूट किया जाता है। यह मैट्रिक्स एक बार कंप्यूट किया जाता है (चूंकि $K$ समय-स्वतंत्र है) और टाइम-स्टेपिंग के दौरान पुन: उपयोग किया जाता है, जिससे लागत काफी कम हो जाती है।

C. समय एकीकरण (Time Integration)

• परिणामी सेमी-डिस्क्रीट सिस्टम को परिवर्तनीय स्टेप साइज के साथ एक स्टिफ सॉल्वर (ode15s MATLAB में) का उपयोग करके हल किया जाता है।
• शुद्ध अवस्थाओं के पास $F'$ द्वारा उत्पन्न स्टिफनेस (stiffness) को संभालने के लिए, लेखक एक नियमित पोटेंशियल ($F_\omega$) का भी परीक्षण करते हैं जो सीमाओं के पास $F$ के टेलर विस्तार से मेल खाता है, जिससे दीर्घकालिक सिमुलेशन में बड़े टाइम स्टेप्स की अनुमति मिलती है।

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

1. सैद्धांतिक विस्तार: लेखकों ने 2D और 3D में लिप्सचिट्ज़ (बहुभुज) डोमेन पर नॉनलोकल CH सिस्टम के लिए सुव्यवस्थितता (well-posedness) और संतुलन की ओर अभिसरण को सिद्ध किया, जो पिछले परिणामों का विस्तार है जो केवल चिकने डोमेन तक सीमित थे।
2. नवीन संख्यात्मक योजना: एक सिंगुलर कर्नेल के लिए स्पेक्ट्रल एलिमेंट दृष्टिकोण विकसित किया जो निम्नलिखित को जोड़ता है:
   • सिंगुलर भाग का विश्लेषणात्मक एकीकरण (टेलर विस्तार के माध्यम से)।
   • नियमित भाग का उच्च-क्रम स्पेक्ट्रल एकीकरण (मल्टीशेप डिकंपोजिशन के माध्यम through)।
     यह नो-फ्लक्स स्थितियों के तहत सिंगुलर कर्नेल और लॉगैरिद्मिक पोटेंशियल के साथ नॉनलोकल CH सिस्टम का पहला उच्च-रिज़ॉल्यूशन कंप्यूटेशन है।
3. त्रुटि नियंत्रण: उन्होंने त्रुटि अनुमान प्रदान किए जो दिखाते हैं कि विधि नियमित भाग के लिए स्पेक्ट्रल अभिसरण और सिंगुलर सुधार के लिए $O(\varepsilon^2 \log \varepsilon)$ प्राप्त करती है, जिसे संख्यात्मक रूप से सत्यापित किया गया है।
4. ओपन सोर्स: पुनरुत्पादकता (reproducibility) को सुगम बनाने के लिए कोड और डेटा (GitHub: DDFT-Modelling/NLCH) जारी किया गया है।

4. संख्यात्मक परिणाम (Numerical Results)

लेखकों ने न्यूटनियन कर्नेल $K_\eta = \eta K$ के साथ एक यूनिट स्क्वायर डोमेन ($\Omega = (0,1)^2$) पर व्यापक सिमुलेशन किए।

• पैरामीटर $\eta$ (स्केलिंग):
  • $\eta > 0$ (प्रतिकर्षी/Repulsive): सिस्टम विसरित (diffusively) व्यवहार करता है, एक होमोजेनियस अवस्था ($\rho \equiv 0$) की ओर अभset होता है।
  • $\eta < 0$ (आकर्षक/Attractive): सिस्टम चरण पृथक्करण (phase separation) प्रदर्शित करता है।
    • मध्यम नकारात्मक $\eta$ के लिए, धीमा चरण पृथक्करण होता है।
    • बड़े नकारात्मक $\eta$ के लिए, तीव्र चरण पृथक्करण होता है जिससे तीक्ष्ण इंटरफेस और स्टेप-फंक्शन जैसे प्रोफाइल बनते हैं।
• प्रारंभिक स्थितियाँ (Initial Conditions):
  • आवधिक तरंगों (periodic waves) और कॉम्पैक्टली सपोर्टेड फंक्शनों का परीक्षण किया गया।
  • देखा गया कि सिस्टम विशिष्ट पैटर्न (जैसे कोनों पर सांद्रता या स्टेप फंक्शन का निर्माण) वाली स्थिर अवस्थाओं की ओर विकसित होता है।
• नियमितीकरण (Regularization): यह प्रदर्शित किया गया कि नियमित पोटेंशियल $F_\omega$ का उपयोग करने से बिना भौतिक परिणाम बदले, दीर्घकालिक स्थिर अवस्थाओं (मशीन परिशुद्धता तक पहुँचने के लिए) के लिए काफी तेज़ गणना संभव है, क्योंकि समाधान स्वाभाविक रूप से सिंगुलैरिटी से दूर रहता है।
• कर्नेल मिश्रण (Kernel Mixture): न्यूटनियन कर्नेल और एक मोलिफायर (mollifier) के रैखिक संयोजन का परीक्षण किया गया। इसने एक ऐसे कर्नेल का निर्माण किया जिसमें सकारात्मक/नकारात्मक प्रविष्टियों का मिश्रण था, जिससे जटिल, पैटर्न-युक्त संरचनाएं उत्पन्न हुईं जो शुद्ध न्यूटनियन कर्नेल के साथ नहीं देखी गईं।

5. महत्व (Significance)

• कम्प्यूटेशनल दक्षता: प्रस्तावित विधि सीमित कम्प्यूटेशनल संसाधनों के साथ उच्च-रिज़ॉल्यूशन समाधान की अनुमति देती है। सिंगुलैरिटी के लिए विश्लेषणात्मक सुधारों और कन्वोल्यूशन मैट्रिक्स को प्रीकंप्यूट करके, यह प्रत्येक टाइम-स्टेप पर घने मैट्रिक्स-वेक्टर गुणन की $O(N^2)$ लागत से बचती है।
• भौतिक सटीकता: लॉगैरिद्मिक पोटेंशियल का उपयोग यह सुनिश्चित करता है कि भौतिक प्रतिबंध $|\rho| < 1$ का सम्मान किया जाए, जो द्विआधारी मिश्र धातुओं और जैविक प्रणालियों (जैसे ट्यूमर विकास मॉडल) के मॉडलिंग के लिए महत्वपूर्ण है।
• बहुमुखी प्रतिभा (Versatility): "मल्टीशेप" दृष्टिकोण लचीला है और इसे विभिन्न डोमेन (जैसे डिस्क, क्यूब) या विभिन्न सिंगुलर कर्नेल के लिए अनुकूलित किया जा सकता है।
• भविष्य के अनुप्रयोग: लेखक काह्न–हिलियर्ड–नेवियर-स्टोक्स (Cahn–Hilliard–Navier–Stokes) समीकरणों या ट्यूमर विकास के लिए हॉकिन्स–दारुड (Hawkins–Daruud) मॉडल जैसे युग्मित प्रणालियों पर इस ढांचे को लागू करने की क्षमता को रेखांकित करते हैं, जहाँ नॉनलोकल अंतःक्रियाएं और सिंगुलर पोटेंशियल महत्वपूर्ण होते हैं।

संक्षेप में, यह शोध पत्र सिंगुलर नॉनलोकल PDEs के कठोर गणितीय विश्लेषण और व्यावहारिक, उच्च-निष्ठा (high-fidelity) संख्यात्मक सिमुलेशन के बीच के अंतर को पाटता है, जो जटिल चरण पृथक्करण घटनाओं का अध्ययन करने के लिए एक मजबूत उपकरण प्रदान करता है।