Error analysis of a fully discrete structure-preserving finite element scheme for a diffuse-interface model of tumour growth

तकनीकी सारांश: ट्यूमर वृद्धि के एक डिफ्यूज-इंटरफेस मॉडल के लिए पूर्णतः विविक्त संरचना-संरक्षण परिमित तत्व योजना का त्रुटि विश्लेषण

समस्या विवरण
यह कार्य मिश्रण के निरंतर सिद्धांत (continuum theory of mixtures) पर आधारित ट्यूमर वृद्धि के एक ऊष्मप्रवैगिकी रूप से सुसंगत, चार-प्रजाति मॉडल के संख्यात्मक सन्निकटन को संबोधित करता है। मॉडल ट्यूमर सेल वॉल्यूम फ्रैक्शन ($u$) के लिए एक काह्न-हिलियर्ड (Cahn–Hilliard) प्रकार के समीकरण को पोषक तत्व सांद्रता ($n$) के लिए एक गैररेखीय प्रतिक्रिया-विसरण समीकरण के साथ जोड़ता है। यह प्रणाली हेल्महोल्ट्ज़ मुक्त ऊर्जा कार्यात्मक (Helmholtz free energy functional) द्वारा नियंत्रित होती है और निरंतर स्तर पर एक विसरणात्मक ऊर्जा नियम (dissipative energy law) को स्वीकार करती है। जबकि पिछले अध्ययनों ने विभिन्न संख्यात्मक योजनाओं (गैररेखीय कॉनवेक्स स्प्लिटिंग और समय-अर्ध-विविक्त SAV दृष्टिकोणों सहित) का प्रस्ताव दिया है, इस विशिष्ट युग्मित प्रणाली के लिए तीन आयामों में एक पूर्णतः विविक्त, रैखिक और संरचना-संरक्षण योजना के साथ कठोर त्रुटि विश्लेषण का साहित्य में अभाव था।

कार्यप्रणाली
लेखक एक रैखिक, पूर्णतः विविक्त परिमित तत्व विधि प्रस्तावित करते हैं जो निरंतर मॉडल के संरचनात्मक गुणों को संरक्षित करती है। कार्यप्रणाली निम्नलिखित को संयोजित करती है:

1. स्केलर ऑक्जिलरी वेरिएबल (SAV) सूत्रीकरण: गैररेखीय मुक्त-ऊर्जा घनत्व $g(u)$ को संभालने के लिए, प्रणाली को एक ऑक्जिलरी वेरिएबल $r(t) = \sqrt{E_1[u] + B}$ का उपयोग करके पुनर्गठित किया जाता है, जहाँ $E_1$ मुक्त ऊर्जा का गैर-द्विघात भाग है। यह ऊर्जा स्थिरता बनाए रखते हुए एक रैखिक योजना के निर्माण की अनुमति देता है।
2. स्थानिक विवक्तीकरण (Spatial Discretization):
   • काह्न-हिलियर्ड भाग के लिए एक मिश्रित परिमित तत्व विधि (mixed finite element method) का उपयोग किया जाता है ताकि चौथे क्रम के समीकरण को युग्मित द्वितीय-क्रम समीकरणों के निकाय में कम किया जा सके।
   • प्रतिक्रिया-विसरण समीकरण के लिए मानक अनुरूप परिमित तत्वों (standard conforming finite elements) का उपयोग किया जाता है।
3. कालिक विवक्तीकरण (Temporal Discretization): एक प्रथम-क्रम अर्ध-निहित (semi-implicit) टाइम-स्टेपिंग योजना का उपयोग किया जाता है।

परिणामी बीजगणितीय प्रणाली प्रत्येक समय चरण (time step) पर रैखिक होती है। SAV पद के कार्यान्वयन को शेरमैन-मोरिसन सूत्र (Sherman–Morrison formula) का उपयोग करके स्पार्स परिमित तत्व ब्लॉक मैट्रिक्स में एक रैंक-वन सुधार के माध्यम से कुशलतापूर्वक संभाला जाता है, जिससे घने मैट्रिक्स संयोजन या गैररेखीय सॉल्वर की आवश्यकता नहीं होती है।

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

• संरचनात्मक गुण: योजना अनियंत्रित रूप से ऊर्जा-स्थिर (एक संशोधित ऊर्जा कार्यात्मक के लिए असतत ऊर्जा अपव्यय नियम को संतुष्ट करते हुए) और द्रव्यमान-संरक्षण (mass-preserving) है।
• सुव्यवस्थितता (Well-posedness): प्रत्येक समय चरण में रैखिक बीजगणितीय प्रणाली के समाधान का अस्तित्व और विशिष्टता स्थापित की गई है।
• त्रुटि अनुमान: लेखक $L^2$, $H^1$, और $L^\infty$ मानदंडों में कठोर त्रुटि अनुमान प्राप्त करते हैं। विश्लेषण प्रदर्शित करता है कि:
  • समय में प्रथम-क्रम सटीकता ($O(\tau)$) है।
  • स्थान में इष्टतम सटीकता ($L^2$ और $H^1$ के लिए $O(h^{q+1})$, जहाँ $q$ बहुपद डिग्री है) है।
• विश्लेषण में तकनीकी नवाचार:
  • पूर्ववर्ती समय-अर्ध-विविक्त विश्लेषणों (जैसे 2D में) के विपरीत, यह कार्य 3D उत्तल डोमेन (convex domain) में पूर्णतः विविक्त सेटिंग को संबोधित करता है।
  • लेखक पूर्णतः विविक्त समाधानों के लिए $H^2$ नियमितता अनुमानों की कमी (जो केवल $H^1$ हैं) को दूर करने के लिए असतत सोबोलेव एम्बेडिंग और असतत गैलियर्डो-निएरेनबर्ग असमानताओं का उपयोग करके सूक्ष्म अनुमान प्राप्त करते हैं।
  • एक प्रमुख चरण $L^\infty$ में संख्यात्मक समाधानों की परिबद्धता को सिद्ध करना और परिमित तत्व समाधान तथा इसके रစ် प्रोजेक्शन (Ritz projection) के बीच सुपरकन्वर्जेंस अनुमान स्थापित करना है।

परिणाम
सैद्धांतिक परिणामों को FEniCS पैकेज का उपयोग करके संख्यात्मक सिमुलेशन के माध्यम से मान्य किया गया है। प्रयोगों में शामिल हैं:

1. एकल ट्यूमर वृद्धि: स्थिर वृद्धि और द्रव्यमान संरक्षण का प्रदर्शन।
2. एकाधिक ट्यूमरों का एकत्रीकरण: दो ट्यूमरों के विलय का अनुकरण।
3. अभिसरण दरें: सैद्धांतिक प्रथम-क्रम कालिक अभिसरण और इष्टतम स्थानिक अभिसरण दरों को सत्यापित करना, यहाँ तक कि कम नियमितता वाले प्रारंभिक डेटा के लिए भी।
4. 3D कीमोटैक्टिक वृद्धि: पोषक तत्वों के स्रोतों की ओर कीमोटैक्टिक व्यवहार प्रदर्शित करने वाले एक तीन-आयामी ट्यूमर का अनुकरण। विशेष रूप से, सिमुलेशन तब भी स्थिर रहता है जब पर्याप्त कोएर्सिविटी (coercivity) स्थिति ($\lambda > \chi_0^2 \delta$) जिसे सैद्धांतिक विश्लेषण में उपयोग किया गया था, सख्ती से संतुष्ट नहीं होती है, जो सिद्ध सीमाओं से परे इसकी मजबूती का सुझाव देती है।

महत्व
यह शोध पत्र दावा करता है कि यह इस विशिष्ट ट्यूमर वृद्धि मॉडल के लिए पहले पूर्णतः विविक्त, रैखिक और संरचना-संरक्षण संख्यात्मक योजना को प्रदान करके साहित्य के अंतराल को भरता है, जो कि 3D में कठोर त्रुटि विश्लेषण के साथ पूर्ण है। यह विधि महत्वपूर्ण है क्योंकि यह द्रव्यमान संरक्षण और ऊर्जा स्थिरता की गारंटी देते हुए गैररेखीय सॉल्वरों की कम्प्यूटेशनल लागत से बचती है। इसका विश्लेषण उच्च आयामों में पूर्णतः विविक्त सेटिंग्स में संरचना-संरक्षण योजनाओं की समझ को विस्तृत करता है, जहाँ गैररेखीय पदों को नियंत्रित करना समय-अर्ध-विविक्त या निम्न-आयामी मामलों की तुलना में काफी अधिक सूक्ष्म होता है। संख्यात्मक प्रयोग पुष्टि करते हैं कि विधि मजबूत है और एकत्रीकरण एवं कीमोटैक्सिस जैसी विशिष्ट घटनाओं को पकड़ने में सक्षम है।