Axial Symmetric Navier Stokes Equations and the Beltrami /anti Beltrami spectrum in view of Physics Informed Neural Networks

यहाँ पिएत्रो फ्रे (Pietro Fré) के शोध पत्र "Axial Symmetric Navier Stokes equations and the Beltrami/anti Beltrami spectrum in view of Physics Informed Neural Networks" का विस्तृत तकनीकी सारांश दिया गया है।

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

यह शोध पत्र विस्कस (viscous) और अघुलनशील (incompressible) तरल पदार्थों के लिए नेवियर-स्टोक्स (NS) समीकरणों के सटीक या उच्च-परिशुद्धता वाले सन्निकट समाधान खोजने की चुनौती को संबोधित करता है। विशेष रूप से, यह एक बेलनाकार ट्यूब (cylindrical tube) के भीतर अक्षीय सममित प्रवाह (axially symmetric flows) पर केंद्रित है, जिसमें अनुदैर्ध्य दिशा (longitudinal direction) में आवधिक सीमा स्थितियाँ (periodic boundary conditions) लागू हैं।

• संदर्भ: पारंपरिक कम्प्यूटेशनल फ्लूइड डायनेमिक्स (CFD) विभेदक ऑपरेटरों (differential operators) के विवेकीकरण (discretization) पर निर्भर करता है, जो अक्सर अंतर्निहित विश्लेत्मक संरचनाओं और समरूपताओं (symmetries) को अस्पष्ट कर देता है।
• लक्ष्य: लेखक एक नवीन दृष्टिकोण का प्रस्ताव करते हैं जिसका उपयोग फिजिक्स इन्फॉर्म्ड न्यूरल नेटवर्क (PINNs) के माध्यम से किया जाता है। हालाँकि, यहाँ लक्ष्य नेटवर्क को ग्रिड पर सीधे समाधान का अनुमान लगाने के लिए प्रशिक्षित करना नहीं है, बल्कि NS समीकरणों को एक विशिष्ट कार्यात्मक आधार (functional basis) के गुणांकों (coefficients) से जुड़ी एक बीजगणितीय समस्या में बदलना है। यह छिपे हुए विश्लेक नियमों और समरूपताओं की खोज करने की अनुमति देता है।
• विशिष्ट फोकस: यह अध्ययन टर्बुलेंस (turbulence) की शुरुआत में बेल्ट्रमी प्रवाह (जहाँ वेग भंवर/vorticity के समानांतर होता है, $\nabla \times \mathbf{u} = \lambda \mathbf{u}$) और एंटी-बेल्ट्रमी प्रवाह की भूमिका की जांच करता है, जो 3-टोरी ($T^3$) पर किए गए पिछले कार्यों का विस्तार है।

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

A. ज्यामितीय पुनर्गठन (Geometric Reformulation)

यह शोध पत्र डिफरेंशियल ज्योमेट्री (एक्सटीरियर कैलकुलस) का उपयोग करके नेवियर-स्टोक्स समीकरणों को पुनर्गठित करता है:

• वेग क्षेत्र $\mathbf{u}$ को एक 1-फॉर्म $\Omega$ के रूप में माना जाता है।
• NS समीकरण को इस प्रकार लिखा जाता है:
  $$-d(p + \frac{1}{2}\|\mathbf{u}\|^2) = \partial_t \Omega + i_{\mathbf{u}} d\Omega - \nu \Delta_g \Omega - \mathbf{f}$$
• स्थिर प्रवाह (steady flows) के लिए, जहाँ बरनौली फलन (Bernoulli function) स्थिर रहता है, समीकरण गैर-रेखीय पद ($i_{\mathbf{u}} d\Omega$) और विस्कस पद ($\nu \Delta_g \Omega$) के बीच संतुलन में बदल जाता है।

B. कार्यात्मक आधार का निर्माण (Construction of the Functional Basis)

मुख्य गणितीय योगदान एक कॉम्पैक्ट ट्यूब पर अक्षीय सममित वर्ग-समावेशी (square-integrable) वेक्टर क्षेत्रों ($L^2_{axial}$) के लिए एक पूर्ण ऑर्थोगोनल आधार का निर्माण है।

• टोपोलॉजी: डोमेन को एक सिलेंडर के रूप में मॉडल किया गया है जिसके आधार आपस में जुड़े हुए हैं (टोपोलॉजिकल रूप से एक टोरस का हिस्सा), जो $z$-दिशा में आवधिकता लागू करता है और $r=1$ पर रेडियल वेग को शून्य करता है।
• आइजनफंक्शंस (Eigenfunctions): आधार का निर्माण हारमोनिक 1-फॉर्म्स (लैपलेस-बेल्ट्रमी ऑपरेटर के आइजनस्टेट्स) से किया गया है।
• विघटन (Decomposition): प्रत्येक ऊर्जा शेल (energy shell) के लिए, जिसे स्पेक्ट्रल इंडिसेस $(n, k)$ द्वारा परिभाषित किया गया है (जहाँ $n$ बेसेल ज़ीरोज़ से संबंधित है और $k$ अनुदैर्ध्य वेवनंबर से), 6-आयामी आइजनस्पेस तीन भागों में विभाजित होता है:
  1. दो बेल्ट्रमी मोड (बाएँ हाथ की घूर्णन, आइजनवैल्यू $+\varpi$)।
  2. दो एंटी-बेल्ट्रमी मोड (दाएँ हाथ की घूर्णन, आइजनवैल्यू $-\varpi$)।
  3. दो क्लोज्ड (अघूर्णी/Irrotational) मोड (शून्य आइजनवैल्यू, ग्रेडिएंट फ्लो)।
• आधार फलन (Basis Functions): रेडियल निर्भरता को बेसेल फलनों ($J_0, J_1$) का उपयोग करके व्यक्त किया जाता है, जिन्हें सीमा स्थितियों को पूरा करने के लिए उनके ज़ीरोज़ ($j_{1,n}$) द्वारा स्केल किया जाता है। अनुदैर्ध्य निर्भरता के लिए ट्रिगोनोमेट्रिक फलनों ($\sin, \cos$) का उपयोग किया जाता है।

C. "डायमंड प्रोडक्ट" (The "Diamond Product")

NS समीकरण में गैर-रेखीय संवहन पद ($\mathbf{u} \cdot \nabla \mathbf{u}$) को डायमंड प्रोडक्ट ($\diamond$) नामक एक बीजगणितीय क्रिया में मैप किया जाता है।

• यह उत्पाद मैनिफोल्ड पर सामान्यीकृत प्रति-सममित (antisymmetric) वेक्टर उत्पाद का प्रतिनिधित्व करता है।
• दो आधार तत्वों का उत्पाद उसी आधार में विस्तारित किया जाता है:
  $$U_{s_1} \diamond U_{s_2} = \sum_s C(s_1, s_2 | s) U_s$$
• गुणांक $C(s_1, s_2 | s)$ विश्लेषणात्मक रूप से प्राप्त किए जाते हैं और बेसेल फलनों के उत्पादों के इंटीग्रल्स पर निर्भर करते हैं।

D. बीजगणितीय पुनरावृत्ति में न्यूनीकरण (Reduction to Algebraic Recurrence)

वेग क्षेत्र $\mathbf{u} = \sum c_s U_s$ को इस आधार में विस्तारित करके, विभेदक NS समीकरण गुणांकों $c_s$ के लिए द्विघाती बीजगणितीय समीकरणों (quadratic algebraic equations) की एक श्रेणी में बदल जाता है:
$$c[s] = \frac{1}{\nu \varpi(s)^2} \sum_{s_1, s_2} \varpi(s_2) c[s_1] c[s_2] C(s_1, s_2 | s)$$
यह प्रणाली एक रिकर्सिव ऑप्टिमाइज़ेशन एल्गोरिदम (एक विशेषीकृत PINN) के माध्यम से हल की जाने वाली प्रस्तावित है, जो इन द्विघाती संबंधों के अवशेष (residual) पर आधारित लॉस फंक्शन को कम करता है।

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

1. नया कार्यात्मक आधार: यह शोध पत्र एक कॉम्पैक्ट ट्यूब पर अक्षीय सममित हार्मोनिक 1-फॉर्म्स (बेल्ट्रमी, एंटी-बेल्ट्रमी और क्लोज्ड) का एक पूर्ण, स्पष्ट आधार प्रदान करता है, जो सीमा स्थितियों को कठोरता से संभालने के लिए बेसेल फलनों के एक विशिष्ट रीस्केलिंग का उपयोग करता है।
2. त्रि-स्तरीय बेल्ट्रमी इंडेक्स (Tripartite Beltrami Index): 3-टोरी पर पिछले कार्यों के विपरीत, जहाँ बेल्ट्रमी इंडेक्स बाइनरी (चिरल) था, लेखक प्रदर्शित करते हैं कि बेलनाकार ज्यामिति में, स्पेक्ट्रल इंडेक्स त्रि-स्तरीय (tripartite) है, जो इनके बीच अंतर करता है:
   • बाएँ-घूर्णन (बेल्ट्रमी)
   • दाएँ-घूर्णन (एंटी-बेल्ट्रमी)
   • अघूर्णी (क्लोज्ड फॉर्म्स)
     यह उस भौतिक तंत्र को स्पष्ट करता है जहाँ बेल्ट्रमी और एंटी-बेल्ट्रमी मोड गैर-रेखीय अंतःक्रिया के माध्यम से अघूर्णी प्रवाह में विलीन हो सकते हैं।
3. NS का बीजगणितीय निरूपण: यह शोध पत्र सफलतापूर्वक फ्लूइड डायनेमिक्स की जटिल PDEs को विस्तार गुणांकों पर द्विघाती पुनरावृत्ति संबंधों (quadratic recurrence relations) के सेट में बदल देता है। यह समस्या को संख्यात्मक PDE समाधान से बीजगणितीय अनुकूलन (algebraic optimization) की ओर ले जाता है।
4. डायमंड प्रोडक्ट अपघटन: लेखक डायमंड प्रोडक्ट के आधार तत्वों के स्ट्रक्चर कांस्टेंट्स (गुणांकों) की गणना स्पष्ट रूप से करते हैं, जो प्रस्तावित न्यूरल नेटवर्क एल्गोरिदम के लिए आवश्यक "इंटरैक्शन नियम" प्रदान करते हैं।
5. जीरो-मोड विश्लेषण: ज़ीरो-मोड्स (जहाँ $k=0$ या $n=0$) का विस्तृत विश्लेषण प्रदान किया गया है, जो यह पहचानता है कि कौन से मोड सीमा स्थितियों के तहत भौतिक रूप से स्वीकार्य हैं (विशेष रूप से, केवल कुछ हेलिकल मोड ही जीवित रहते हैं)।

4. परिणाम (Results)

• विश्लेषणात्मक व्युत्पत्ति: शोध पत्र प्रत्येक स्पेक्ट्रल शेल $(n,k)$ के लिए 6 आधार वेक्टर्स के स्पष्ट रूप को व्युत्पन्न करता है और सामान्य हार्मोनिक गुणांकों को बेल्ट्रमी/एंटी-बेल्ट्रमी/क्लोज्ड घटकों से संबंधित करने वाला ट्रांसफॉर्मेशन मैट्रिक्स प्रदान करता है।
• उत्पादों का अभिसरण (Convergence of Products): परिशिष्ट B में संख्यात्मक परीक्षणों से पता चलता है कि दो आधार फलनों का उत्पाद (जैसे, $G^{[n_1]}_0 \times G^{[n_2]}_1$) अत्यधिक सटीकता के साथ उसी आधार में पुनः विस्तारित किया जा सकता है (विशेष रूप से, $n \approx n_1 + n_2$ तक ट्रंकेट करके)। यह प्रस्तावित बीजगणितीय दृष्टिकोण की व्यवहार्यता की पुष्टि करता है।
• स्ट्रीमलाइन विज़ुअलाइज़ेशन: शोध पत्र एकल मोड और सुपरपोजिशन के लिए स्ट्रीमलाइन्स को विज़ुअलाइज़ करता है, जो दिखाता है कि जबकि एकल मोड टोरोइडल सतह बनाते हैं, सुपरपोजिशन जटिल, विकृत संरचनाएं बनाते हैं, जो अक्षीय सममित बाधा के भीतर भी समाधान स्थान की समृद्धि को दर्शाता है।

5. महत्व और भविष्य का कार्य (Significance and Future Work)

• वैज्ञानिक अंतर्दृष्टि: यह दृष्टिकोण केवल "ब्लैक बॉक्स" संख्यात्मक समाधानों से आगे बढ़कर NS समाधानों की विश्लेत्मक संरचना को उजागर करने का लक्ष्य रखता है। अनुकूलन एल्गोरिथ्म द्वारा पाए गए गुणांकों का विश्लेषण करके, लेखक सटीक समाधानों को नियंत्रित करने वाले नियमों का अनुमान लगाने की आशा करते हैं।
• PINN रणनीति: यह शोध पत्र एक भविष्य के प्रकाशन [10] के लिए सैद्धांतिक आधार तैयार करता है जहाँ एक फिजिक्स इन्फॉर्म्ड न्यूरल नेटवर्क बनाया जाएगा। मानक PINNs के विपरीत, जो नेटवर्क आउटपुट को डिफरेंशिएट करते हैं, यह विधि नेटवर्क का उपयोग एक पूर्व-निर्धारित, सिमेट्री-एडेप्टेड आधार के गुणांकों को अनुकूलित करने के लिए करती है।
• सिमेट्री ब्रेकिंग (Symmetry Breaking): यह ढांचा विस्तार योग्य है। एक बार अक्षीय सममित समाधान मिल जाने के बाद, टूटी हुई सिमेट्री के इर्रड्यूसिबल रिप्रेजेंटेशन (irreducible representations) में गड़बड़ी (perturbations) पेश करके सिमेट्री को चरण-दर-चरण तोड़ा जा सकता है, जो संभावित रूप से पूर्णतः 3D टर्बुलेंट समाधानों की ओर ले जा सकता है।
• भौतिक प्रासंगिकता: प्रवाह स्पेक्ट्रम की त्रि-स्तरीय प्रकृति की पहचान, सीमित प्रवाह (confined flows) में ऊर्जा कैसे चिरल (घूर्णन) और अघूर्णी मोड के बीच स्थानांतरित होती है, इस पर एक नया दृष्टिकोण प्रदान करती है, जो पाइपलाइन, रक्त वाहिकाओं और औद्योगिक रिएक्टरों जैसे अनुप्रयोगों के लिए प्रासंगिक है।

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