Second-Order Area/Volume-Preserving PFEMs for Surface Diffusion via Simpson--Boole Geometric Identities

तकनीकी सारांश: सिम्पसन-बुल (Simpson–Boole) ज्यामितीय पहचानों के माध्यम से सतह विसरण (Surface Diffusion) के लिए द्वितीय-क्रम वाले क्षेत्र/आयतन-संरक्षण करने वाले PFEMs

समस्या विवरण
सतह विसरण (Surface diffusion) एक चतुर्थ-क्रम (fourth-order) वाला ज्यामितीय विकास प्रवाह है जहाँ अभिलंब वेग (normal velocity) वक्रता के लाप्लास-बेल्ट्रामी ऑपरेटर द्वारा निर्धारित होता है। दो आयामों में बंद वक्रों और तीन आयामों में बंद सतहों के लिए, यह प्रवाह स्वाभाविक रूप से घिरे हुए क्षेत्रफल या आयतन को संरक्षित करता है जबकि परिधि या सतह क्षेत्र को कम करता है। यद्यपि विभिन्न संख्यात्मक दृष्टिकोण मौजूद हैं (जैसे, लेवल-सेट, फेज-फील्ड), पैरामीट्रिक परिमित तत्व विधियाँ (PFEM), विशेष रूप से बैरेट, गारके और नूर्बर्ग (BGN) द्वारा विकसित ढांचा, प्रेरित स्पर्शरेखीय गति (induced tangential motion) के माध्यम से अच्छी मेश गुणवत्ता बनाए रखने की अपनी क्षमता के कारण पसंद किए जाते हैं। हालाँकि, मौजूदा द्वितीय-क्रम वाले BGN-प्रकार के स्कीमें अक्सर असतत स्तर (discrete level) पर घिरे हुए क्षेत्रफल या आयतन को सटीक रूप से संरक्षित करने में विफल रहती हैं, या उन्हें इन ज्यामितीय बाधाओं को लागू करने के लिए सहायक लैग्रेंज मल्टीप्लायर्स (Lagrange multipliers) की आवश्यकता होती है, जो कार्यान्वयन को जटिल बनाता है और स्थिरता को प्रभावित कर सकता है।

कार्यप्रणाली
लेखक एक नए वर्ग के द्वितीय-क्रम-इन-टाइम (second-order-in-time) PFEMs का प्रस्ताव करते हैं जो लैग्रेंज मल्टीप्लायर्स को पेश किए बिना सटीक क्षेत्रफल (2D) और आयतन (3D) संरक्षण प्राप्त करते हैं। मुख्य कार्यप्रणाली एक द्विघात टेम्पोरल इंटरपोलेशन पथ (quadratic temporal interpolation path) के साथ सटीक ज्यामितीय भिन्नता पहचानों से सीधे संरक्षण पदों को व्युत्पन्न करने पर आधारित है।

1. द्विघात टेम्पोरल इंटरपोलेशन (Quadratic Temporal Interpolation): यह विधि तीन क्रमिक समय स्तरों ($m-1, m, m+1$) को एक द्विघात बहुपद $X_h(\alpha)$ का उपयोग करके जोड़ती है। यह पथ घिरे हुए ज्यामितीय मात्रा के परिवर्तन की दर के लिए सटीक सूत्रों के व्युत्पन्न की अनुमति देता है।
2. ज्यामितीय पहचान और क्वाड्रैचर (Geometric Identities and Quadrature):
   • 2D (वक्र): द्विघात पथ के साथ घिरे हुए क्षेत्रफल का परिवर्तन एक ऐसा इंटीग्रैंड (integrand) उत्पन्न करता है जो टेम्पोरल पैरामीटर में अधिकतम तीसरे घात का बहुपद है। फलस्वरूप, सिम्पसन का नियम (Simpson's rule) इस भिन्नता का सटीक मूल्यांकन करता है। यह तीन समय स्तरों पर वेग अनुमानों के एक विशिष्ट रैखिक संयोजन को प्रदान करता है जो शून्य शुद्ध क्षेत्रफल परिवर्तन की गारंटी देता है।
   • 3D (सतहें): सतह के पैच के ओरिएंटेड एरिया वेक्टर के माध्यम से सतह के आयतन का भिन्नता शामिल है। इंटीग्रैंड अधिकतम पांच डिग्री का बहुपद बन जाता है। इसलिए, इस भिन्नता का सटीक मूल्यांकन करने के लिए बुल के नियम (Boole's rule) (एक उच्च-क्रम न्यूटन-कोट फार्मूला) की आवश्यकता होती है। इसके लिए दो अतिरिक्त "हाफ-टाइम" (half-time) ज्यामितीय स्तरों ($m \pm 1/2$) की आवश्यकता होती है।
3. प्रेडिक्शन रणनीति (Prediction Strategy): मौजूदा कुशल कार्यान्वयन के साथ अनुकूलता बनाए रखने के लिए, स्कीमें एक BGN-प्रेडिक्टेड ऑक्सिलरी ज्योमेट्री का उपयोग करती हैं। एक प्रथम-क्रम BGN चरण नए-समय वक्र/सतह (और 3D में हाफ-टाइम सतहों) की भविष्यवाणी करता है। रूढ़िवादी अभिलंबों (conservative normal vectors) का निर्माण इन अनुमानित ज्यामितिओं का उपयोग करके किया जाता है, जिससे नॉर्मल वेक्टर की गैर-रैखिकता (nonlinearity) का संरक्षण प्रमाण से अलगाव हो जाता है।
4. असतत स्कीमें (Discrete Schemes):
   • BGN/AC (2D): एक स्कीम जो सिम्पसन के गुणांकों (1/6, 4/6, 1/6) द्वारा भारित BDF2-प्रकार, क्रैंक-निकोलसन-लीपफ्रॉग (Crank–Nicolson-leapfrog), और बैकवर्ड एंडपॉइंट सन्निकटन को जोड़ती है।
   • BGN/VC (3D): एक स्कीम जिसमें पांच पद शामिल हैं जिन्हें बुल के गुणांकों (7/90, 32/90, 12/90, 32/90, 7/90) द्वारा भारित किया गया है, जिसमें मुख्य समय स्तर और अनुमानित हाफ-टाइम स्तर शामिल हैं।
5. सॉल्वर (Solver): परिणामी गैर-रैखिक प्रणालियों को न्यूटन-प्रकार की पुनरावृत्ति विधि का उपयोग करके हल किया जाता है। लीनियरलाइजेशन (linearization) अज्ञात ज्यामिति पर अभिलंब वेक्टर्स की निर्भरता को ध्यान में रखता है, जिसमें 3D में ओरिएंटेड एरिया वेक्टर्स के लीनियरलाइजेशन पर विशेष ध्यान दिया गया है।

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

• ज्यामितीय संरक्षण का विस्तार: यह शोध रैखिक इंटरपोलेशन के बजाय द्विघात इंटरपोलेशन का उपयोग करके बाओ और झाओ के प्रथम-क्रम संरचना-संरक्षण विचार को द्वितीय-क्रम सेटिंग में विस्तारित करता है।
• मल्टीप्लायर्स के बिना सटीक संरक्षण: लेखक सिम्पसन और बुल नियमों के आधार पर सटीक असतत भिन्नता सूत्रों को व्युत्पन्न करते हैं। यह उन्हें स्कीमें में आंतरिक रूप से सटीक संरक्षण गुणों को बनाने की अनुमति देता है, जिससे सहायक लैग्रेंज मल्टीप्लायर्स की आवश्यकता समाप्त हो जाती है।
• BGN ढांचे के साथ अनुकूलता: प्रस्तावित स्कीमें BGN-प्रेडिक्टेड ज्योमेट्री पर असेंबल करने के लिए डिज़ाइन की गई हैं। यह सुनिश्चित करता है कि वे BGN-प्रकार की विधियों की अनुकूल मेश गुणवत्ता और स्पर्शरेखीय गति के गुणों को बनाए रखें, जिससे प्रत्यक्ष एक्सट्रपलेशन से जुड़ी मेश गिरावट से बचा जा सके।
• एकीकृत ढांचा: यह कार्य विभिन्न मौजूदा द्वितीय-क्रम दृष्टिकोणों को एक एकल संरक्षण निर्माण के घटकों के रूप में एकीकृत करता है, जो मुख्य रूप से नॉर्मल वेक्टर्स की परिभाषा में भिन्न होते हैं।

संख्यात्मक परिणाम
लेखक सैद्धांतिक गुणों को मान्य करने के लिए व्यापक संख्यात्मक प्रयोग प्रस्तुत करते हैं:

• अभिसरण (Convergence): 2D (दीर्घवृत्त) और 3D (दीर्घवृत्त) परीक्षण दोनों ही समय में द्वितीय-क्रम अभिसरण प्रदर्शित करते हैं, जो सैद्धांतिक सटीकता के अनुरूप है।
• संरक्षण (Conservation): सभी परीक्षणों में, सापेक्ष क्षेत्रफल हानि (2D) और आयतन हानि (3D) गैर-रैखिक सॉल्वर टॉलरेंस (लगभग $10^{-14}$) के स्तर पर बनी रहती है, जो सटीक संरक्षण की पुष्टि करती है। इसके विपरीत, मानक द्वितीय-क्रम BGN स्कीमें (CNLF, BDF2, PC) $10^{-5}$ के क्रम में विचलन (drifts) दिखाती हैं।
• मेश गुणवत्ता: मेश अनुपात (एलिमेंट साइज की एकरूपता) अच्छी तरह से नियंत्रित रहता है और इक्विडिस्ट्रिब्यूशन (equidistribution) के करीब पहुँच जाता है, जो मानक BGN विधियों के समान है। यह चुनौतीपूर्ण ज्यामिति जैसे कि फ्लावर कर्व्स, पीनट शेप्स और पिंच-ऑफ के करीब पहुँचते सिगार-जैसे सतहों के लिए भी बना रहता है।
• मजबूती (Robustness): स्कीमें जटिल विकास को सफलतापूर्वक संभालती हैं, जिसमें सिगार-नुमा सतह की नेकिंग अस्थिरता (necking instability) शामिल है, जो पिंच-ऑफ तक स्थिरता और मेश गुणवत्ता बनाए रखती है।

महत्व और दावे
यह पेपर दावा करता है कि प्रस्तावित निर्माण सतह विसरण के लिए उच्च-क्रम सटीकता और सटीक ज्यामितीय संरक्षण के लिए एक कठोर मार्ग प्रदान करता है, बिना लैग्रेंज मल्टीप्लायर्स के कम्प्यूटेशनल ओवरहेड या जटिलता के। द्विघात इंटरपोलेशन पथ के साथ सटीक ज्यामितीय पहचानों (सिम्पसन और बुल) का लाभ उठाकर, यह विधि "दोनों दुनियाओं का सर्वश्रेष्ठ" परिदृश्य प्राप्त करती है: यह प्रेडिक्टर-आधारित BGN स्कीम्स के व्यावहारिक मेश व्यवहार को बनाए रखती है और सख्त संरक्षण नियमों को लागू करती है। लेखक इस कार्य को संरक्षण तंत्र पर एक केंद्रित प्रगति के रूप में प्रस्तुत करते हैं, यह नोट करते हुए कि इसे ऊर्जा स्थिरता के लिए परिधि/सतह क्षेत्र न्यूनीकरण कानूनों को सख्ती से लागू करने हेतु लैग्रेंज-मल्टीप्लायर उपचारों के साथ जोड़ा जा सकता है। यह कार्य उच्च-क्रम विस्तार और अधिक जटिल ज्यामितीय प्रवाह के लिए एक आधारभूत कदम के रूप में प्रस्तुत किया गया है।