Distributed optimal control problems governed by poroelasticity equations

तकनीकी सारांश: पोरोइलास्टिसिटी समीकरणों द्वारा नियंत्रित वितरित इष्टतम नियंत्रण समस्याएँ

समस्या विवरण
यह शोध पत्र बायोट के कंसोलिडेशन मॉडल (Biot's consolidation model) में पोरोइलास्टिसिटी (poroelasticity) के अंतर्गत नियंत्रित वितरित इष्टतम नियंत्रण समस्याओं के संख्यात्मक सन्निकटन (numerical approximation) को संबोधित करता है। यह भौतिक प्रणाली अंतरालीय द्रव प्रवाह (interstitial fluid flow) और एक छिद्रयुक्त माध्यम (porous medium) के विरूपण के बीच की अंतःक्रिया का वर्णन करती है, जो भू-भौतिकी, जैविक ऊतक यांत्रिकी और कार्बन पृथक्करण जैसे अनुप्रयोगों के लिए एक महत्वपूर्ण मॉडल है। विशिष्ट नियंत्रण समस्या का उद्देश्य एक लागत फलन (cost functional) को कम करना है जो वांछित लक्ष्य मानों ($u_C, p_C$) से ठोस विस्थापन ($u$) और द्रव दबाव ($p$) के विचलन को मापता है, जबकि नियंत्रण चर के परिमाण को दंडित (penalize) करता है। यहाँ नियंत्रण को डोमेन के भीतर वितरित द्रव स्रोत/सिंक पद ($m$) के रूप में पहचाना गया है। यह अध्ययन अर्ध-स्थैतिक (quasi-static) कंसोलिडेशन पर केंद्रित है जहाँ त्वरण प्रभावों की उपेक्षा की जाती है, और यह विशेष रूप से "पोरोइलास्टिक लॉकिंग" (poroelastic locking) से जुड़ी गणितीय और संख्यात्मक चुनौतियों को संबोधित करता है, जो तब होती है जब विशिष्ट भंडारण गुणांक ($s_0$) शून्य के करीब होता है, जिससे द्रव और लोचदार माध्यम लगभग असंपीड्य (incompressible) हो जाते हैं।

कार्यप्रणाली
लेखक एक नवीन द्वि-क्षेत्र सममित सूत्रीकरण (two-field symmetric formulation) प्रस्तावित करते हैं और उसका विश्लेषण करते हैं जहाँ प्राथमिक अज्ञात चर ठोस विस्थापन और द्रव दबाव हैं। कार्यप्रणाली निम्नलिखित चरणों के माध्यम से आगे बढ़ती है:

1. सतत सूत्रीकरण और सु-स्थापना (Well-Posedness):

   • लेखक मोमेंटम संतुलन के समय व्युत्पन्न (time derivative) को लेकर गवर्निंग समीकरणों को पुनर्गठित करते हैं, जिससे सिस्टम विभेदक समीकरणों के एक सेट में परिवर्तित हो जाता है।
   • वे उपयुक्त फलन स्थान (function spaces) ($V$ विस्थापन के लिए, $Q$ दबाव के लिए) को परिभाषित करते हैं और सैडल-पॉइंट सिद्धांत (saddle-point theory) और लेडीज़ेन्स्काया-बाबुष्का-ब्रेज़ी (LBB) स्थिति का उपयोग करके कमजोर समाधान (weak solution) की सु-स्थापना स्थापित करते हैं।
   • वे प्रत्यक्ष विधि (direct method in the calculus of variations) का उपयोग करके इष्टतम नियंत्रण की अस्तित्व और विशिष्टता को सिद्ध करते हैं, जो लागत फलन की उत्तलता (convexity) और अवस्था ऑपरेटर (state operator) की रैखिकता का लाभ उठाता है।
   • प्रथम-क्रम आवश्यक और पर्याप्त इष्टतमता स्थितियों को एक एडजॉइंट स्टेट दृष्टिकोण के माध्यम से व्युत्पन्न किया जाता है, जो नियंत्रण के लिए एक परिवर्तनशील असमानता (variational inequality) की ओर ले जाता है।

2. विविक्त योजना (Discrete Scheme):

   • स्थानिक विवक्तीकरण (Spatial Discretization): डोमेन को परिमित तत्वों ($V_h \times Q_h$) का उपयोग करके विविक्त किया जाता है जो इन्फ-सुप (inf-sup) स्थिति को संतुष्ट करते हैं (जैसे, MINI-elements या Taylor-Hood elements)।
   • कालिक विवक्तीकरण (Temporal Discretization): एक बैकवर्ड यूलर स्कीम (विशेष रूप से समय में dG(0) विधि) का उपयोग किया जाता है, जिसके परिणामस्वरूप एक पूर्णतः विविक्त प्रणाली प्राप्त होती है।
   • नियंत्रण विवक्तीकरण: यह शोध पत्र नियंत्रण चर के लिए एक विवरणात्मक विवक्तीकरण दृष्टिकोण (variational discretization approach) (हिंज़े द्वारा प्रस्तावित) का उपयोग करता है। इसका अर्थ है कि नियंत्रण स्थान को स्पष्ट रूप से विविक्त नहीं किया जाता है; इसके बजाय, विविक्त नियंत्रण को स्वीकार्य नियंत्रण सेट पर विविक्त एडजॉइंट अवस्था को प्रोजेक्ट करके प्राप्त किया जाता है। यह उन सन्निकटन त्रुटियों से बचता है जो आमतौर पर नियंत्रण स्थान को विविक्त करने से उत्पन्न होती हैं।

3. त्रुटि विश्लेषण:

   • लेखक पूर्णतः विविक्त योजना के लिए ए प्रायर त्रुटि अनुमान (a priori error estimates) व्युत्पन्न करते हैं।
   • विश्लेषण कुल त्रुटि को विवक्तीकरण त्रुटियों (निश्चित नियंत्रण के लिए) और निरंतर एवं विविक्त इष्टतम नियंत्रणों के बीच के अंतर के कारण होने वाली विक्षोभ त्रुटियों (perturbation errors) में विभाजित करता है।
   • फॉरवर्ड स्टेट और बैकवर्ड एडजॉइंट दोनों प्रणालियों के लिए त्रुटियों को सीमित करने हेतु स्थिरता अनुमान स्थापित किए गए हैं।

प्रमुख योगदान
यह शोध पत्र साहित्य में तीन प्राथमिक योगदान देता है:

• नवीन सममित सूत्रीकरण: यह एक द्वि-क्षेत्र सममित सूत्रीकरण (विस्थापन और दबाव) पेश करता है जो विशेष रूप से बायोट मॉडल के भीतर इष्टतम नियंत्रण समस्याओं के लिए तैयार किया गया है, जो पिछले चार-क्षेत्र या मिश्रित सूत्रीकरणों से भिन्न है।
• अस्तित्व और विशिष्टता: यह नए सूत्रीकरण की सु-स्थापना को कड़ाई से सिद्ध करता है और इष्टतम नियंत्रण के अस्तित्व और विशिष्टता को स्थापित करता है, साथ ही प्रथम-क्रम इष्टतमता स्थितियों का व्युत्पन्न भी करता है।
• ए प्रायर त्रुटि अनुमान: यह पूर्णतः विविक्त योजना के लिए एक पूर्ण ए प्रायर त्रुटि विश्लेषण प्रदान करता है। विशेष रूप से, यह विस्थापन और दबाव के लिए $L^\infty(0, T; H^1(\Omega))$ और $L^2(0, T; L^2(\Omega))$ नॉर्म्स में, और नियंत्रण चर के लिए $L^2(0, T; L^2(\Omega))$ नॉर्म में त्रुटि अनुमान व्युत्पन्न करता है, जो मानक नियमितता धारणाओं के तहत इष्टतम अभिसरण दर (convergence rates) प्राप्त करता है।

परिणाम

• सैद्धांतिक: शोध पत्र सिद्ध करता है कि संख्यात्मक योजना $O(h^s + \Delta t)$ के क्रम के साथ अभिसरित होती है, जहाँ $h$ स्थानिक मेश आकार है, $\Delta t$ समय अंतराल है, और $s$ समाधान की नियमितता से संबंधित है। विशेष रूप से, यदि समाधान पर्याप्त नियमित ($s=1$) है, तो अभिसरण स्थान और समय दोनों में रैखिक ($O(h + \Delta t)$) होता है।
• संख्यात्मक: एक यूनिट स्क्वायर डोमेन पर एक संख्यात्मक प्रयोग सैद्धांतिक निष्कर्षों की पुष्टि करता है। ओपन-सोर्स लाइब्रेरी FEniCS का उपयोग करके, लेखकों ने एक फॉरवर्ड-बैकवर्ड एल्गोरिदम का उपयोग करके इष्टतमता प्रणाली को हल किया। परिणाम दर्शाते हैं कि प्रस्तावित योजना निम्नतम-क्रम के MINI और Taylor-Hood तत्वों का उपयोग करते समय विस्थापन, दबाव और नियंत्रण चरों के लिए अनुमानित इष्टतम अभिसरण दर प्राप्त करती है।

महत्व
लेखक इस कार्य को पोरोइलास्टिक प्रणालियों द्वारा नियंत्रित इष्टतम नियंत्रण समस्याओं के संख्यात्मक विश्लेषण के लिए एक आधारभूत कदम के रूप में प्रस्तुत करते हैं। वे उल्लेख करते हैं कि जबकि विरूपित छिद्रयुक्त माध्यमों में द्रव प्रवाह के लिए इष्टतम नियंत्रण का अध्ययन किया गया है (जैसे, बोसियू आदि द्वारा लुप्त भंडारण गुणांक की धारणा के तहत), सामान्य बायोट मॉडल के लिए ऐसी समस्याओं के संख्यात्मक सन्निकटन की पहले जांच नहीं की गई थी। एक नए सममित सूत्रीकरण की सु-स्थापना को स्थापित करके और एक पूर्णतः विविक्त योजना के लिए कठोर त्रुटि अनुमान प्रदान करके, यह शोध पत्र साहित्य के एक अंतराल को भरता है, जो जैविक ऊतक मॉडलिंग और भू-भौतिकीय सिमुलेशन जैसे अनुप्रयोगों में पोरोइलास्टिक प्रक्रियाओं को नियंत्रित करने के लिए एक मजबूत सैद्धांतिक ढांचा प्रदान करता है।