Explicit Discrete Solution for Some Optimization Problems and Estimations with Respect to the Exact Solution

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

यह शोध पत्र एक बहुआयामी परिबद्ध डोमेन $\Omega \subset \mathbb{R}^n$ (विशेष रूप से एक आयताकार डोमेन में विश्लेषित) में पॉइसन समीकरण (Poisson equation) द्वारा नियंत्रित स्थिर-अवस्था ऊष्मा चालन (steady-state heat conduction) समस्याओं को संबोधित करता है। यह अध्ययन दो विशिष्ट प्रणालियों पर केंद्रित है:

• सिस्टम (S): इसमें मिश्रित सीमा स्थितियाँ (mixed boundary conditions) शामिल हैं: $\Gamma_1$ पर डिरिचलेट (निश्चित तापमान $b$), $\Gamma_2$ पर न्यूमैन (निश्चित ऊष्मा प्रवाह $q$), और $\Gamma_3$ पर एडियाबेटिक (शून्य प्रवाह)।
• सिस्टम (S$_\alpha$): यह (S) के समान है, लेकिन $\Gamma_1$ पर डिरिचलेट स्थिति को एक रोबिन सीमा स्थिति (संवहनी ऊष्मा प्रवाह/convective heat flux) द्वारा प्रतिस्थापित किया गया है, जिसमें एक गुणांक $\alpha > 0$ शामिल है।

इन प्रणालियों से जुड़ी तीन इष्टतम नियंत्रण समस्याएँ (optimal control problems) प्रत्येक सिस्टम के लिए ($P_i$ और $P_{i\alpha}$ के लिए $i=1,2,3$) हैं, जहाँ नियंत्रण चर (control variables) निम्नलिखित हैं:

1. $g$: वितरित आंतरिक ऊर्जा स्रोत।
2. $q$: $\Gamma_2$ पर सीमा ऊष्मा प्रवाह।
3. $b$: $\Gamma_1$ पर पर्यावरणीय तापमान।

इसका उद्देश्य एक द्विघातीय लागत फलन (quadratic cost functional) $J$ को न्यूनतम करना है, जो वांछित लक्ष्य $z_d$ से अवस्था (तापमान) के विचलन और नियंत्रण चर के परिमाण (नियमितीकरण/regularization) के बीच संतुलन बनाता है।

यद्यपि आयताकार डोमेन पर इन समस्याओं के लिए स्पष्ट निरंतर समाधान (explicit continuous solutions) पहले से ज्ञात थे, इस शोध पत्र का लक्ष्य परिमित अंतर विधियों (finite difference methods) का उपयोग करके स्पष्ट विविक्त समाधान (explicit discrete solutions) प्राप्त करना और उनके अभिसरण (convergence) और त्रुटि अनुमानों का कठोरता से विश्लेषण करना है।

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

लेखक निरंतर प्रणालियों और अनुकूलन समस्याओं को विविक्त (discretize) करने के लिए परिमित अंतर विधि (Finite Difference Method - FDM) का उपयोग करते हैं।

• विवक्तीकरण (Discretization):

  • डोमेन को एक समान मेश आकार $h = x_0/n$ के साथ विविक्त किया गया है।
  • आंतरिक नोड्स (Interior Nodes): द्वितीय व्युत्पन्न ($\partial^2 u / \partial x^2$) के लिए एक शास्त्रीय केंद्रित द्वितीय-क्रम परिमित अंतर सन्निकटन (centered second-order finite difference approximation) का उपयोग किया गया है।
  • सीमा स्थितियाँ (Boundary Conditions - Standard):
    • न्यूमैन स्थितियों का सन्निकटन एक प्रथम-क्रम पश्च अंतर (first-order backward difference) का उपयोग करके किया गया है।
    • रोबिन स्थितियों का सन्निकटन एक प्रथम-क्रम अग्र अंतर (first-order forward difference) का उपयोग करके किया गया है।
  • इससे ट्राइएडिगोनल (tridiagonal) गुणांक आव्यूह वाली रैखिक प्रणालियाँ ($S_h$ और $S_{h\alpha}$) प्राप्त होती हैं।

• अनुकूलन (Optimization):

  • निरंतर लागत फलनों को विविक्त लागत फलनों ($J^h_i$ और $J^h_{i\alpha}$) के निर्माण के लिए विविक्त किया गया है।
  • चूंकि नियंत्रण चर को स्थिर माना गया है, इसलिए विविक्त अनुकूलन समस्याएँ एकल चर के द्विघातीय फलनों को न्यूनतम करने में बदल जाती हैं, जिससे स्पष्ट सूत्र (explicit formulas) प्राप्त करना संभव होता है (जैसे कि इष्टतम विविक्त नियंत्रण $g^h_{op}, q^h_{op}, b^h_{op}$)।

• अभिसरण विश्लेषण (Convergence Analysis):

  • शोध पत्र विविक्त समाधानों ($u^h, g^h_{op}$, आदि) की सटीक निरंतर समाधानों ($u, g_{op}$, आदि) के साथ तुलना करते हुए स्पष्ट त्रुटि अनुमान प्राप्त करता है।
  • यह दोहरे अभिसरण (double convergence) व्यवहार का विश्लेषण करता है: जैसे-जैसे मेश आकार $h \to 0$ और संवहनी गुणांक $\alpha \to \infty$ होता है।

• बेहतर सन्निकटन (Improved Approximation - Section 7):

  • प्रथम-क्रम सीमा सन्निकटन के कारण होने वाले निम्न अभिसरण क्रम ($O(h)$) को संबोधित करने के लिए, लेखक न्यूमैन और रोबिन सीमा स्थितियों के लिए एक तीन-बिंदु परिमित अंतर योजना (three-point finite difference scheme) प्रस्तावित करते हैं।
  • यह संशोधन एक "घोस्ट पॉइंट" (ghost point) फॉर्मूलेशन पेश करता है, जो प्रभावी रूप से सीमा उपचार की सटीकता को बढ़ाता है।

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

1. स्पष्ट विविक्त समाधान (Explicit Discrete Solutions): यह शोध पत्र सभी तीन प्रकार की अनुकूलन समस्याओं (वितरित स्रोत, फ्लक्स, और तापमान) के लिए विविक्त अवस्था चरों और इष्टतम नियंत्रण चरों के लिए बंद-रूप (closed-form) विश्लेषणात्मक व्यंजक प्रदान करता है। यह एक महत्वपूर्ण योगदान है क्योंकि यह अनुकूलन चरण के लिए पुनरावृत्ति संख्यात्मक सॉल्वर (iterative numerical solvers) की आवश्यकता को समाप्त करता है।
2. कठोर त्रुटि अनुमान (Rigorous Error Estimation): लेखक सिद्ध करते हैं कि मानक सीमा सन्निकटन का उपयोग करने पर विविक्त समाधान $L^2$ नॉर्म में निरंतर समाधानों के साथ प्रथम-क्रम सटीकता ($O(h)$) के साथ अभिसरित होते हैं। वे इन त्रुटि सीमाओं के लिए स्पष्ट स्थिरांक (constants) भी प्रदान करते हैं।
3. दोहरा अभिसरण विश्लेषण (Double Convergence Analysis): यह अध्ययन सीमाओं के क्रमविनिमेयता (commutativity) को स्थापित करता है: रोबिन समस्या का विविक्त समाधान $\alpha \to \infty$ होने पर डिरिचलेट समस्या के विविक्त समाधान की ओर अभिसरित होता है, और दोनों ही $h \to 0$ होने पर निरंतर समाधान की ओर अभिसरित होते हैं।
4. क्रम सुधार (Order Improvement): सीमा स्थितियों के लिए तीन-बिंदु सन्निकटन का उपयोग करके, लेखक यह प्रदर्शित करते हैं कि अवस्था चर के लिए वैश्विक अभिसरण का क्रम $O(h)$ से सुधरकर $O(h^2)$ हो जाता है, जबकि व्युत्पन्न त्रुटि $O(h)$ बनी रहती है।

4. परिणाम (Results)

• अभिसरण दरें (Convergence Rates):
  • मानक FDM: त्रुटि $\|u - u^h\|_H$ का व्यवहार $C_1 h$ जैसा है।
  • सुधरा हुआ FDM (3-point): त्रुटि $\|u - \tilde{u}^h\|_H$ का व्यवहार $D_1 h^2$ जैसा है।
• इष्टतम नियंत्रण (Optimal Controls): इष्टतम नियंत्रणों के स्पष्ट सूत्र दर्शाते हैं कि विविक्त इष्टतम नियंत्रण ($g^h_{op}$, आदि) निरंतर वाले चरों से $h$ के समानुपाती पदों द्वारा भिन्न होते हैं।
• संख्यात्मक सिमुलेशन (Numerical Simulations):
  • यूनिट स्क्वायर डोमेन ($x_0=y_0=1$) पर सिमुलेशन सैद्धांतिक अभिसरण दरों की पुष्टि करते हैं।
  • तालिकाएँ और चित्र दर्शाते हैं कि मेश आकार $h$ को आधा करने से $L^2$ त्रुटि लगभग आधी हो जाती है, जो $O(h)$ व्यवहार की पुष्टि करती है।
  • सिमुलेशन यह सत्यापित करते हैं कि जैसे-जैसे $\alpha$ बढ़ता है, रोबिन सिस्टम ($S_\alpha$) के समाधान डिरिचलेट सिस्टम ($S$) के समाधानों की ओर अभिसरित होते हैं।

5. महत्व (Significance)

• बेंचमार्किंग (Benchmarking): प्राप्त स्पष्ट विविक्त समाधान, अनुकूल नियंत्रण समस्याओं के लिए सामान्य प्रयोजन संख्यात्मक विधियों (जैसे फाइनाइट एलीमेंट मेथड्स) की सटीकता और विश्वसनीयता का परीक्षण करने के लिए कठोर बेंचमार्क के रूप में कार्य करते हैं, जहाँ स्पष्ट समाधान दुर्लभ होते हैं।
• दक्षता (Efficiency): इष्टतम नियंत्रणों को स्पष्ट रूप से (पुनरावृत्ति के बिना) गणना करने की क्षमता इस श्रेणी की समस्याओं के लिए कम्प्यूटेशनल लागत को काफी कम कर देती है।
• सीमा उपचार अंतर्दृष्टि (Boundary Treatment Insight): यह शोध पत्र वैश्विक अभिसरण क्रमों पर सीमा स्थिति विवक्तीकरण के महत्वपूर्ण प्रभाव को रेखांकित करता है। यह प्रदर्शित करता है कि भले ही आंतरिक भाग के लिए द्वितीय-क्रम योजना का उपयोग किया जाए, फिर भी प्रथम-क्रम सीमा योजना वैश्विक सटीकता को $O(h)$ तक सीमित कर देती है, और $O(h^2)$ वैश्विक सटीकता प्राप्त करने के लिए उच्च-क्रम सीमा योजनाओं की आवश्यकता होती है।
• सीमाएँ और भविष्य का कार्य (Limitations and Future Work): वर्तमान विश्लेषण स्पष्ट समाधान प्राप्त करने की सुविधा के लिए आयताकार डोमेन तक ही सीमित है। लेखक भविष्य के कार्य के रूप में इस पद्धति को सामान्य डोमेन और समन्वय प्रणालियों (ध्रुवीय, गोलाकार) तक विस्तारित करने का सुझाव देते हैं।

संक्षेप में, यह शोध पत्र सटीक विविक्त सूत्र प्रदान करके और उनके अभिसरण गुणों को सिद्ध करके, निरंतर अनुकूल नियंत्रण सिद्धांत और विविक्त संख्यात्मक कार्यान्वयन के बीच के अंतर को पाटता है, साथ ही बेहतर सीमा सन्निकटन के माध्यम से संख्यात्मक सटीकता बढ़ाने के लिए एक व्यावहारिक तकनीक भी प्रस्तुत करता है।