Hamiltonian Interface Dynamics for Reduced-Order Optimization of Incompressible Mixing

तकनीकी सारांश: असम्पीड्य मिश्रण के कम-क्रम अनुकूलन (Reduced-Order Optimization) के लिए हैमिल्टनियन इंटरफ़ेस गतिकी

समस्या विवरण
यह शोध पत्र दो-आयामी असम्पीड्य प्रवाहों (incompressible flows) में मिश्रण के इष्टतम नियंत्रण (optimal control) को संबोधित करता है। शासक मॉडल एक निष्क्रिय स्केलर $\theta$ के लिए परिवहन समीकरण $\partial_t \theta + v \cdot \nabla \theta = 0$ है, जिसे एक सीमित समतलीय डोमेन $\Omega$ में एक अपसरण-मुक्त (divergence-free) वेग क्षेत्र $v$ द्वारा विस्थापित किया जाता है। उद्देश्य एक स्टिरिंग फील्ड $v$ (जो एक स्ट्रीम फंक्शन $\psi$ से व्युत्पन्न है) को डिजाइन करना है जो मिश्रण दक्षता को अधिकतम करे।

पारंपरिक रूप से, मिश्रण को ऋणात्मक सोबोलेव नॉर्म्स (negative Sobolev norms), विशेष रूप से होमोजेनियस $\dot{H}^{-1}(\Omega)$ नॉर्म का उपयोग करके मापा जाता है। हालांकि यह प्रभावी है, लेकिन इस नॉर्म के सापेक्ष सीधे अनुकूलन करने के लिए प्रत्येक पुनरावृत्ति (iteration) पर पूर्ण परिवहन आंशिक अंतर समीकरण (PDE) और उसके एडजॉइंट (adjoint) को हल करने की आवश्यकता होती है, जो गणनात्मक रूप से अत्यधिक महंगा है। लेखक एक ज्यामितीय विकल्प प्रस्तावित करते हैं: सीधे स्केलर क्षेत्र को अनुकूलित करने के बजाय, वे एक प्रवाहित होने वाले सामग्री इंटरफ़ेस (material interface) की लंबाई को अनुकूलित करते हैं। यह दृष्टिकोण इस तथ्य का लाभ उठाता है कि 2D असम्पीड्य प्रवाहों में, वेग क्षेत्र एक स्ट्रीम-फंक्शन निरूपण ($v = \nabla^\perp \psi$) स्वीकार करता है, जिससे विशेषता गतिकी (characteristic dynamics) हैमिल्टनियन हो जाती है।

कार्यप्रणाली
लेखक एक कम-क्रम (reduced-order) ढांचा विकसित करते हैं जो यूलेरियन PDE अनुकूलन को एक लैग्रेंजियन हैमिल्टनियन नियंत्रण समस्या से बदल देता है।

1. हैमिल्टनियन इंटरफ़ेस गतिकी: विधि एक सामग्री इंटरफ़ेस $\gamma_0$ को ट्रैक करती है जो $s$ द्वारा पैरामीटराइज्ड है। इंटरफ़ेस की स्थिति $X(t, s)$ और स्पर्श रेखा सदिश (tangent vector) $a(t, s) = \partial_s X(t, s)$ स्ट्रीम फंक्शन $\psi$ द्वारा संचालित युग्मित साधारण अंतर समीकरणों (ODEs) के अनुसार विकसित होते हैं:
   $$\partial_t X = J \nabla_x \psi(t, X), \quad \partial_t a = J \nabla^2_x \psi(t, X) a$$
   जहाँ $J$ एक सिम्प्लेक्टिक मैट्रिक्स है। उद्देश्य कार्यात्मक (objective functional) इस इंटरफ़ेस की अंतिम लंबाई है, $L(\psi; T) = \int_0^1 |a(T, s)| ds$।

2. इष्टतम नियंत्रण सूत्रीकरण (Optimal Control Formulation): स्ट्रीम फंक्शन को चिकने आधार फलनों (smooth basis functions) $\{h_k\}_{k=1}^N$ द्वारा स्पैन किए गए एक परिमित-आयामी उपसमष्टि (finite-dimensional subspace) तक सीमित किया गया है। नियंत्रण चर समय-निर्भर गुणांक $u(t) \in \mathbb{R}^N$ हैं। यह एक नियमितीकृत लागत कार्यात्मक (regularized cost functional) को न्यूनतम करता है:
   $$J_{\lambda, \epsilon}(\psi) = -\int_0^1 \ell_\epsilon(a(T, s)) ds + \frac{\lambda}{2} \int_0^T \|\psi(t, \cdot)\|_{L^2(\Omega)}^2 dt$$
   जहाँ $\ell_\epsilon$ यूक्लिडियन नॉर्म का एक चिकना सन्निकटन (smooth approximation) है ताकि अवकलनीयता सुनिश्चित हो सके। इस बाधाबद्ध समस्या के लिए एक मिनिमाइज़र की उपस्थिति सिद्ध की गई है।

3. एडजॉइंट-आधारित ग्रेडिएंट गणना: लेखक इंटरफ़ेस स्थिति और स्पर्श रेखा सदिश के लिए निरंतर एडजॉइंट सिस्टम $(p, q)$ को व्युत्पन्न करते हैं। यह मार्कर प्रक्षेप पथों (marker trajectories) के साथ ले जाए जाने वाले एडजॉइंट चरों के आधार पर नियंत्रण गुणांकों $u$ के सापेक्ष रिड्यूस्ड ग्रेडिएंट की गणना करने की अनुमति देता है। ग्रेडिएंट को एडजॉइंट चरों के संदर्भ में स्पष्ट रूप से व्यक्त किया गया है।

4. संख्यात्मक विवेकीकरण (Numerical Discretization): फॉरवर्ड और एडजॉइंट ODEs को इम्प्लिसिट मिडपॉइंट नियम (implicit midpoint rule) का उपयोग करके विवेकीकृत किया गया है, जो एक सिम्प्लेक्टिक और स्व-एडजॉइंट इंटीग्रेटर है। यह विकल्प सुनिश्चित करता है कि असतत एडजॉइंट (discrete adjoint), पूर्णतः विवेकीकृत उद्देश्य फलन के व्युत्पन्न के बीजगणितीय रूप से सुसंगत है। अनुकूलन को पोलाक-रिबिएर कंजुगेट-ग्रेडिएंट विधि का उपयोग करके किया जाता है जिसमें एक आर्मिजो लाइन सर्च और $L^2$-बॉल बाधा पर प्रोजेक्शन शामिल है।

प्रमुख योगदान

• ज्यामितीय कम-क्रम मॉडल: यह शोध पत्र स्थापित करता है कि 2D असम्पीड्य मिश्रण को इंटरफ़ेस स्थिति और स्पर्श रेखा को नियंत्रित करने वाले ODEs के एक बंद सिस्टम के माध्यम से मॉडल किया जा सकता है, जो पूर्ण यूलेरियन परिवहन PDE को प्रतिस्थापित करता है।
• कम-क्रम मॉडल पर एडजॉइंट-आधारित अनुकूलन: एक परिमित-आयामी स्ट्रीम-फंक्शन एंसेटल (ansatz) के भीतर इंटरफ़ेस लंबाई को अनुकूलित करने के लिए एक कठोर ढांचा प्रस्तुत किया गया है, जिसमें अस्तित्व के प्रमाण और निरंतर एडजॉइंट सिस्टम का व्युत्पन्न शामिल है।
• सिम्प्लेक्टिक डिस्क्रीट एडजॉइंट: लेखक "डिसक्रेटाइज-देन-ऑप्टिमाइज़" दृष्टिकोण का उपयोग करते हुए इम्प्लिसिट मिडपॉइंट नियम का कार्यान्वयन करते हैं। यह गारंटी देता है कि डिस्क्रीट ग्रेडिएंट, डिस्क्रीट कॉस्ट फलन के संबंध में बीजगणितीय रूप से सटीक है, जो नॉनलीन सॉल्वर की सहनशीलता (tolerance) तक सुसंगत है।
• सत्यापन और तुलना: कम-क्रम नियंत्रणों का पूर्ण परिवहन PDE पर मूल्यांकन किया गया है। विधि का व्यवस्थित रूप से एक अत्याधुनिक यूलेरियन सोबोलेव-नॉर्म ऑप्टिमाइज़र के विरुद्ध मूल्यांकन किया गया है।

परिणाम
संख्यात्मक प्रयोग दो बेंचमार्क पर किए गए: एक स्थिर सेलुलर फ्लो और एक डोसवेल फ्रंटोजेनेसिस मॉडल।

• स्थिर प्रवाह (Stationary Flows): दोनों बेंचमार्क ने इंटरफ़ेस लंबाई का बहुपद (polynomial) विकास और $\dot{H}^{-1}$ नॉर्म का बहुपद क्षय प्रदर्शित किया, जो स्वायत्त (autonomous) 2D प्रवाहों के लिए सैद्धांतिक अपेक्षाओं के अनुरूप है।
• अनुकूलित समय-निर्भर प्रवाह: एडजॉइंट-आधारित ऑप्टिमाइज़र ने ऐसे समय-निर्भर नियंत्रणों को सफलतापूर्वक डिजाइन किया जो इंटरफ़ेस के निकट-चरघातांकीय (near-exponential) खिंचाव और $\dot{H}^{-1}$ नॉर्म के चरघातांकीय क्षय को प्रेरित करते हैं।
• गणना दक्षता: जब इसकी तुलना यूलेरियन सोबोलेव-नॉर्म ऑप्टिमाइज़र (जो प्रत्येक चरण में पूर्ण 2D PDE और एडजॉइंट PDE को हल करता है) से की गई, तो इंटरफ़ेस-लंबाई ऑप्टिमाइज़र काफी तेज़ था। सेलुलर फ्लो बेंचमार्क में, इंटरफ़ेस विधि ने गणना समय को 13–21 गुना कम कर दिया जबकि तुलनीय या बेहतर $\dot{H}^{-1}$ क्षय दर प्राप्त की। डोसवेल बेंचमार्क में, गति वृद्धि लगभग 2.7–2.9 गुना थी।
• सरोगेट की सीमाएं: एक प्रयोग जिसने नियंत्रण आधार को $N=2$ से $N=4$ मोड तक बढ़ाया, एक सीमा को प्रकट करता है। जबकि $N=4$ ऑप्टिमाइज़र ने बहुत अधिक लंबी इंटरफ़ेस लंबाई (लगभग 16 गुना लंबी) प्राप्त की, परिणामी $\dot{H}^{-1}$ क्षय दर केवल मध्यम रूप से सुधरी। यह इंगित करता है कि हालांकि इंटरफ़ेस लंबाई मिश्रण के लिए एक प्रभावी प्रॉक्सी है, लेकिन यह जटिल ज्यामितीय शासन (regimes) में पूरी तरह से वफादार नहीं है जहाँ स्थानीय खिंचाव वैश्विक मिश्रण दक्षता में परिवर्तित नहीं हो सकता है।

महत्व और दावे
लेखक दावा करते हैं कि यह पहला ढांचा है जो एक हैमिल्टनियन मार्कर/टैंगेंट कम-क्रम मॉडल, इंटरफ़ेस लंबाई के एडजॉइंट-आधारित अनुकूलन, एक सिम्प्लेक्टिक डिस्क्रीट एडजॉइंट और एक व्यवस्थित मैच-कंडीशन वैलिडेशन को जोड़ता है।

यह शोध पद्धति को $\dot{H}^{-1}$ (मिक्स-नॉर्म) अनुकूलन के सार्वभौमिक सरोगेट के रूप में नहीं, बल्कि एक गणनात्मक रूप से कुशल साथी के रूप में प्रस्तुत करता है। यह ग्लोविंस्की (Glowinski) के "टू-ग्रिड फिलॉसफी" के अनुरूप है, जहाँ एक कम-क्रम मॉडल पर अनुकूलन—जो आवश्यक ज्यामितीय तंत्रों को पकड़ता है—ऐसे नियंत्रण प्रदान कर सकता है जो मूल PDE पर बहुत कम लागत पर अच्छा प्रदर्शन करते हैं। यह कार्य रेखांकित करता है कि इंटरफ़ेस लंबाई मिश्रण के लिए एक शक्तिशाली ज्यामितीय प्रॉक्सी है, लेकिन चेतावनी देता है कि इसकी भविष्य कहने वाली शक्ति का आकलन पूर्ण परिवहन समीकरण पर ए पोस्टेरियो (a posteriori) किया जाना चाहिए, विशेष रूप से जटिल ज्यामितीय संरचनाओं वाले शासन में।