A penalized ϕ-FEM scheme for the Poisson Dirichlet problem

तकनीकी सारांश: एक दंडित (Penalized) $\phi$-FEM योजना: पॉइसन डिरिचलेट समस्या के लिए

समस्या विवरण
यह शोध पत्र जटिल या विकसित होती ज्यामिति (evolving geometries) पर होमोजेनियस डिरिचलेट बाउंड्री कंडीशंस के साथ पॉइसन समीकरण के संख्यात्मक समाधान (numerical solution) को संबोधित करता है। ऐसी समस्याओं में मुख्य चुनौती बाउंड्री-फिटेड मेश (boundary-fitted meshes) का निर्माण करना है, जो काफी प्रतिबंधात्मक हो सकता है। लेखक पॉइसन समस्या पर ध्यान केंद्रित करते हैं:
$$\begin{cases} -\Delta u = f & \text{in } \Omega, \\ u = 0 & \text{on } \Gamma, \end{cases}$$
जहाँ डोमेन $\Omega$ और उसकी सीमा $\Gamma$ को एक लेवल-सेट फंक्शन $\phi$ द्वारा अप्रत्यक्ष रूप से (implicitly) वर्णित किया गया है, जहाँ $\Omega = \{\phi < 0\}$ और $\Gamma = \{\phi = 0\}$ है।

कार्यप्रणाली (Methodology)
लेखक $\phi$-FEM (phi-Finite Element Method) के एक दंडित (penalized) संस्करण का प्रस्ताव करते हैं, जो एक अनफिटेड (unfitted) परिमित तत्व विधि है और बाउंड्री-फिटेड मेश की आवश्यकता को समाप्त करती है। मूल $\phi$-FEM फॉर्मूलेशन के विपरीत, जो पूरे डोमेन पर अज्ञात चर को $u = \phi w$ के रूप में पुनर्व्याख्यायित करता है, प्रस्तावित योजना दंडन पद (penalization term) के माध्यम से बाउंड्री कंडीशंस को दुर्बल रूप से (weakly) लागू करती है।

मुख्य कार्यप्रणाली संबंधी विशेषताएं शामिल हैं:

1. स्थानीय लेवल-सेट आवश्यकता: इस विधि को केवल सीमा के निकटवर्ती सेल्स ($\Omega_h^\Gamma$) पर ही लेवल-सेट फंक्शन $\phi$ की आवश्यकता होती है, न कि वैश्विक स्तर पर। यह एक सहायक चर $p$ को बाउंड्री-निकटवर्ती सेल्स पर परिभाषित करके प्राप्त किया जाता है ताकि $u = \phi p$ को दंडन के माध्यम से लागू किया जा सके।
2. वैरियेशनल फॉर्मूलेशन (Variational Formulation): योजना $(u_h, p_h) \in V_h^{(k)} \times Q_h^{(k)}$ की खोज करती है जो एक वैरियेशनल समीकरण को संतुष्ट करती है जिसमें शामिल हैं:
   • बैकग्राउंड मेश $\Omega_h$ पर मानक ग्रेडिएंट पद।
   • एक बाउंड्री इंटीग्रल पद।
   • एक दंडन पद $\frac{\gamma}{h^2} \int_{\Omega_h^\Gamma} (u_h - \frac{1}{h}\phi_h p_h)(v_h - \frac{1}{h}\phi_h q_h)$ जो $u \approx \phi p$ संबंध को लागू करने के लिए है।
   • एक "घोस्ट पेनल्टी" (Ghost Penalty) स्थिरीकरण पद ($G_h^{lhs}$) जिसमें बाउंड्री-निकटवर्ती फैसेट्स और सेल्स पर नॉर्मल डेरिवेटिव्स और लैप्लासियन्स के जंप शामिल हैं, ताकि कोएर्सिविटी (coercivity) सुनिश्चित की जा सके।
3. स्थिरीकरण (Stabilization): योजना स्थिरता सुनिश्चित करने के लिए, विशेष रूप से बाउंड्री द्वारा कटे हुए तत्वों के लिए, पैरामीटर्स $\gamma$ और $\sigma_D$ के साथ एक घोस्ट पेनल्टी तकनीक का उपयोग करती है।

मुख्य योगदान और सैद्धांतिक परिणाम
यह शोध पत्र इस दंडित योजना का कठोर संख्यात्मक विश्लेषण प्रदान करता है, जिसे पहले पेश तो किया गया था लेकिन विश्लेषित नहीं किया गया था। मुख्य सैद्धांतिक योगदान हैं:

• कोएर्सिविटी (Coercivity): लेखक यह सिद्ध करते हैं कि उपयुक्त मेश और लेवल-सेट फंक्शन के तहत संबद्ध द्विपद रूप (bilinear form) कोएर्सिव है, बशर्ते कि स्थिरीकरण पैरामीटर्स $\gamma$ और $\sigma_D$ पर्याप्त बड़े हों।
• ए प्रियोरी एरर एस्टीमेट्स (A Priori Error Estimates):
  • $H^1$ सेमी-नॉर्म: योजना $O(h^k)$ की इष्टतम अभिसरण (convergence) दर प्राप्त करती है, जहाँ $k$ परिमित तत्व स्थान की बहुपद डिग्री है।
  • $L^2$ नॉर्म: विश्लेषण $O(h^{k+1/2})$ की एक अर्ध-इष्टतम (quasi-optimal) अभिसरण दर प्रदान करता है।
  • रेगुलैरिटी (Regularity): विश्लेषण के लिए पूर्व साहित्य में पाए गए मूल $\phi$-FEM विश्लेषण की तुलना में समाधान $u$ पर कम रेगुलैरिटी की आवश्यकता होती है।
• कंडीशनिंग (Conditioning): सैद्धांतिक परिणाम और संख्यात्मक प्रयोग सुझाव देते हैं कि संबद्ध मैट्रिक्स का कंडीशन नंबर $O(h^{-2})$ के रूप में स्केल करता है।

संख्यात्मक परिणाम
सैद्धांतिक निष्कर्षों को Dolfinx लाइब्रेरी का उपयोग करके पायथन (Python) में कार्यान्वित संख्यात्मक प्रयोगों के माध्यम से सत्यापित किया गया है:

1. 2D जटिल ज्यामिति: एक लिवर जैसी ज्यामिति का परीक्षण किया गया। दंडित योजना ("Dual $\phi$-FEM") ने $L^2$ और $H^1$ दोनों नॉर्म्स में इष्टतम अभिसरण दर प्रदर्शित की, जो मानक फिटेड FEM और मूल "Direct $\phi$-FEM" के तुल्य है।
2. 3D ज्यामिति: एक मैन्युफैक्चरड सॉल्यूशन के साथ एक गोलाकार डोमेन का परीक्षण किया गया। यहाँ भी, इष्टतम अभिसरण देखा गया।
3. प्रदर्शन तुलना: जबकि दोनों $\phi$-FEM वेरिएंट मानक फिटेड FEM की तुलना में तेज़ थे (जटिल मेश जनरेशन से बचने के कारण), दंडित योजना ने अतिरिक्त चर $p$ के कारण बढ़े हुए सिस्टम साइज के कारण डायरेक्ट $\phi$-FEM की तुलना में अधिक कम्प्यूटेशनल लागत उठाई। हालांकि, दंडित योजना लेवल-सेट फंक्शन के संबंध में मजबूती (robustness) प्रदान करती है, क्योंकि इसे वैश्विक स्तर पर $\phi$ के ग्रेडिएंट या दूसरे डेरिवेटिव की आवश्यकता नहीं होती है।

महत्व और दावे
यह शोध पत्र दावा करता है कि यह दंडित दृष्टिकोण मिश्रित फॉर्मूलेशन की मजबूती और इमर्स्ड (immersed) परिमित तत्व विधियों के लचीलेपन को सफलतापूर्वक जोड़ता है। इसका प्राथमिक महत्व निम्नलिखित है:

• कम लेवल-सेट निर्भरता: लेवल-सेट फंक्शन (और उसके डेरिवेटिव्स) की आवश्यकता को बाउंड्री-निकटवर्ती सेल्स तक सीमित करके, यह विधि लेवल-सेट फंक्शन (जैसे कि साइन्ड डिस्टेंस फंक्शन्स में) की सिंगुलैरिटी के प्रति अधिक मजबूत है और इसके लिए समाधान अज्ञात के वैश्विक पुनर्निर्माण की आवश्यकता नहीं होती है।
• अनुकूलता (Compatibility): लेवल-सेट आवश्यकता की स्थानीय प्रकृति इस योजना को न्यूमैन बाउंड्री कंडीशंस के अनुकूल बनाती है, जैसा कि संबंधित कार्यों में प्रस्तुत किया गया है।
• वैधीकरण (Validation): यह कार्य इस विशिष्ट दंडित संस्करण के लिए पहला कठोर त्रुटि विश्लेषण प्रदान करता है, जो इसके सैद्धांतिक अभिसरण गुणों और जटिल ज्यामिति के लिए इसकी व्यावहारिक व्यवहार्यता की पुष्टि करता है।

लेखक नोट करते हैं कि जबकि सैद्धांतिक $L^2$ अनुमान अर्ध-इष्टतम है, अंतिम खंड में संख्यात्मक प्रयोग व्यवहार में $O(h^{k+1})$ की इष्टतम अभिसरण दर का सुझाव देते हैं।