Global existence, blow-up behavior and numerical simulations for a class of chemotaxis-driven fish-mussel systems

तकनीकी सारांश: एक वर्ग के कीमोटैक्सिस-संचालित मछली-मसल (Fish-Mussel) तंत्र के लिए वैश्विक अस्तित्व, ब्लो-अप व्यवहार और संख्यात्मक सिमुलेशन

समस्या विवरण
यह कार्य एक कीमोटैक्सिस-संचालित पारिस्थितिकी तंत्र मॉडल की जांच करता है जो एक परिबद्ध डोमेन $\Omega \subset \mathbb{R}^n$ के भीतर पोषक तत्वों ($u$), मछली की आबादी ($v$), और मसल (mussel) आबादी ($w$) की स्थानिक-कालिक गतिशीलता का वर्णन करता है। मॉडल को होमोजेनियस न्यूमैन सीमा स्थितियों (homogeneous Neumann boundary conditions) के अधीन युग्मित आंशिक अवकल समीकरणों (PDEs) के रूप में तैयार किया गया है। यह तंत्र पोषक तत्वों के प्रसार, पोषक तत्वों से समृद्ध क्षेत्रों की ओर मछली के कीमोटैक्सिस (chemotaxis), और पोषक तत्वों के पुनर्चक्रण में फ़िल्टर-फीडिंग भूमिका निभाने वाले मसल्स को शामिल करता है। नियंत्रित समीकरण (1.1) में एक पैरामीटर $\tau \in \{0, 1\}$ शामिल है, जो एक पैराबोलिक-पैराबोलिक-पैराबोलिक प्रणाली ($\tau=1$) और एक पैराबोलिक-इलिप्टिक-पैराबोलिक प्रणाली ($\tau=0$) के बीच अंतर करता है। अध्ययन का उद्देश्य शास्त्रीय समाधानों (classical solutions) के वैश्विक अस्तित्व और परिबद्धता (boundedness) को स्थापित करना, ब्लो-अप (blow-up) की ओर ले जाने वाली स्थितियों का विश्लेषण करना (विशेष रूप से तीन आयामों में), और इन सैद्धांतिक निष्कर्षों को मान्य करने के लिए एक संख्यात्मक योजना विकसित करना है।

कार्यप्रणाली
लेखक विश्लेषणात्मक तकनीकों और संख्यात्मक सिमुलेशन के संयोजन का उपयोग करते हैं:

1. विश्लेषणात्मक ढांचा:

   • स्थानीय अस्तित्व (Local Existence): एलिप्टिक और पैराबोलिक रेगुलैरिटी थ्योरी पर निर्भर करते हुए, सेमग्रुप विधियों और बानाच फिक्स्ड-पॉइंट प्रमेय का उपयोग करके समय-सीमित अद्वितीय स्थानीय शास्त्रीय समाधान का अस्तित्व स्थापित किया गया है।
   • ए प्रियोरी अनुमान (A Priori Estimates): वैश्विक अस्तित्व को सिद्ध करने के लिए, लेखक $L^\infty$ और $L^p$ बाउंड्स प्राप्त करते हैं। $\tau=1$ के मामले के लिए, कीमोटैक्टिक टर्म $-\chi\nabla \cdot (v\nabla u)$ को संभालने के लिए एक विशिष्ट टेस्ट फंक्शन $\phi(x)$ (लेम्मा 2.5 के माध्यम से निर्मित) का उपयोग किया जाता है। यह कीमोटैक्टिक संवेदनशीलता $\chi$ और पोषक तत्व सीमाओं पर विशिष्ट बाधाओं के तहत मछली की आबादी $v$ के लिए $L^p$ अनुमान प्राप्त करने की अनुमति देता है।
   • विस्तार क्षमता मानदंड (Extensibility Criteria): वैश्विक अस्तित्व की सिद्धि इस बात पर निर्भर करती है कि समाधान का $L^\infty$ मान परिबद्ध रहे, जिससे ब्लो-अप को रोका जा सके। यह न्यूमैन हीट सेमग्रुप के स्मूथिंग गुणों के साथ $L^p$ अनुमानों को जोड़कर प्राप्त किया जाता है।
   • ब्लो-अप विश्लेषण (Blow-Up Analysis): तीन-आयामी मामले ($n=3$) के लिए, लेखक ब्लो-अप समय $T_{max}$ के लिए एक निचला स्तर (lower bound) प्राप्त करते हैं। इसमें एक सहायक ऊर्जा कार्यात्मक $\Psi_\tau$ (1.8 में परिभाषित) का निर्माण करना और इसके विकास के लिए एक डिफरेंशियल इनइक्वैलिटी व्युत्पन्न करना शामिल है। इस इनइक्वैलिटी को एकीकृत करके, ब्लो-अप समय के लिए एक स्पष्ट निचला स्तर प्राप्त किया जाता है।

2. संख्यात्मक दृष्टिकोण:

   • विविक्तकरण (Discretization): स्थानिक विविक्तकरण के लिए पीसवाइज लीनियर ($P1$) टेट्राहेड्रल तत्वों और सामयिक एकीकरण के लिए बैकवर्ड यूलर स्कीम का उपयोग करके एक परिमित तत्व विधि (FEM) विकसित की गई है।
     ो
   • कार्यान्वयन: पूर्णतः युग्मित गैर-रेखीय प्रणाली को PETSc SNES न्यूटन लाइन-सर्च सॉल्वर का उपयोग करके हल किया जाता है, जिसमें रेखीयकृत प्रणालियों को LU गुणनखंडन के साथ GM即可 GMRES प्रीकंडीशनर द्वारा संभाला जाता है। सिमुलेशन FEniCSx का उपयोग करके किए जाते हैं और MATLAB में दृश्यित (visualize) किए जाते हैं।
   • मान्यकरण (Validation): संख्यात्मक योजना को मेश कन्वर्जेंस अध्ययनों और MATLAB के ode45 द्वारा हल किए गए स्थानिक रूप से समरूप ODE मॉडल के विरुद्ध दीर्घकालिक गतिशीलता की तुलना करके मान्य किया जाता है।

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

• वैश्विक अस्तित्व और परिबद्धता (प्रमेय 1.1): यह शोध पत्र आयामों $n \geq 3$ के लिए प्रणाली के अद्वितीय शास्त्रीय समाधानों के वैश्विक अस्तित्व को स्थापित करता है। विशेष रूप से, यदि प्रारंभिक डेटा कुछ रेगुलैरिटी शर्तों को पूरा करता है और पैरामीटर निम्नलिखित असमानता को संतुष्ट करते हैं:
  $$\max \left\{ \tau\|u_0\|_{L^\infty(\Omega)}, \frac{M}{\alpha} \right\} < \min \left\{ \frac{\pi}{\chi}\sqrt{\frac{2}{n}}, \frac{\delta_2}{\sigma_2} \right\}$$
  तो समाधान $(u, v, w)$ समय के साथ समान रूप से परिबद्ध है। यह परिणाम $\tau=0$ और $\tau=1$ दोनों के लिए मान्य है।

• ब्लो-अप समय अनुमान (प्रमेय 1.2): उस स्थिति में जहाँ समाधान सीमित समय में ब्लो-अप होते हैं (विशेष रूप से एक 3D उत्तल डोमेन में), लेखक ब्लो-अप समय $T_{max}$ के लिए एक गणना योग्य निचला स्तर प्राप्त करते हैं। यह स्तर ऊर्जा कार्यात्मक $\Psi_\tau(0)$ और सिस्टम पैरामीटर्स को शामिल करने वाले एक इंटीग्रल के रूप में व्यक्त किया गया है। एक स्पष्ट सरलीकृत निचला स्तर भी प्रदान किया गया है:
  $$T_{max} \geq \frac{1}{2E_\tau \Psi_\tau(0)^2}$$

• संख्यात्मक अभिसरण (Numerical Convergence): प्रस्तावित परिमित तत्व योजना को लीनियर तत्वों के उपयोग के अनुरूप, $L^2(\Omega)$ और $L^\infty(\Omega)$ दोनों नॉर्म्स में दूसरे क्रम (second-order) के अभिसरण को प्राप्त करने के लिए दिखाया गया है।

• संख्यात्मक सिमुलेशन:

  • स्थिरता (Stability): मध्यम पैरामीटर्स के साथ सिमुलेशन दिखाते हैं कि प्रणाली सह-अस्तित्व की स्थिर अवस्था (coexistence steady state) की ओर अभिसरित होती है, जो ODE विश्लेषण द्वारा अनुमानित दीर्घकालिक स्थिरता को मान्य करती है।
  • ब्लो-अप व्यवहार: उच्च प्रारंभिक सांद्रता और विशिष्ट पैरामीटर चयन के साथ एक 3D क्यूबिक डोमेन में संख्यात्मक प्रयोग, मछली की आबादी के घनत्व $v$ के तीव्र, अनियंत्रित विकास को दर्शाते हैं, जो परिमित-समय ब्लो-अप की सैद्धांतिक संभावना की पुष्टि करते हैं। सिमुलेशन दिखाते हैं कि जबकि $v$ ब्लो-अप होता है, पोषक तत्व $u$ और मसल $w$ के घनत्व परिबद्ध रहते हैं, जो व्युत्पन्न सैद्धांतिक सीमाओं के अनुरूप है।

महत्व और दावे
यह शोध पत्र कीमोटैक्टिक मछली और फ़िल्टर-फीडिंग मसल्स वाले एक्वाकल्चर पारिस्थितिकी तंत्र की स्थिरता और संभावित अस्थिरता को समझने के लिए एक कठोर गणितीय आधार प्रदान करने का दावा करता है। पिछले ODE-आधारित मॉडलों को PDE ढांचे तक विस्तारित करके, यह कार्य पारिस्थितिकी तंत्र की गतिशीलता में स्थानिक विषमता (spatial heterogeneity) और निर्देशित गति (chemotaxis) की महत्वपूर्ण भूमिका को उजागर करता है।

लेखक अपने कार्य को सैद्धांतिक विश्लेषण और कम्प्यूटेशनल मॉडलिंग के बीच एक सेतु के रूप में प्रस्तुत करते हैं। वे प्रदर्शित करते हैं कि जबकि विशिष्ट पैरामीटर व्यवस्था के तहत वैश्विक अस्तित्व की गारंटी दी जा सकती है (जो यूट्रोफिकेशन-प्रेरित पतन को रोकता है), उच्च आयामों में यदि ये स्थितियाँ उल्लंघन होती हैं, तो प्रणाली परिमित-समय ब्लो-अप के प्रति संवेदनशील होती है। एक मजबूत परिमित तत्व योजना का विकास जटिल गतिशीलता को देखने की अनुमति देता है, जो सैद्धांतिक ब्लो-अप मानदंडों को सत्यापित करने और स्थिर सह-अस्तित्व से जनसंख्या विस्फोट के संक्रमण को सिम्युलेट करने के लिए एक उपकरण प्रदान करता है। यह अध्ययन नए पारिस्थितिक प्रबंधन रणनीतियों का प्रस्ताव नहीं करता है बल्कि बदलते पर्यावरणीय परिस्थितियों के तहत ऐसे तंत्रों की स्थिरता का आकलन करने के लिए आवश्यक विश्लेषणात्मक और संख्यात्मक उपकरण प्रदान करता है।