On the Numerical Treatment of an Abstract Nonlinear System of Coupled Hyperbolic Equations Associated with the Timoshenko Model

तकनीकी सारांश: टिमोशेंको मॉडल से जुड़े एक अमूर्त अरैखिक युग्मित हाइपरबोलिक समीकरणों के निरूपण का संख्यात्मक उपचार

समस्या का निरूपण
यह कार्य एक वास्तविक हिल्बर्ट स्पेस $H$ में परिभाषित एक अमूर्त अरैखिक युग्मित हाइपरबोलिक समीकरणों की कॉची समस्या (Cauchy problem) को संबोधित करता है। यह प्रणाली एक गैर-रैखिक गतिशील टिमोशेंको बीम मॉडल के सामान्यीकरण के रूप में कार्य करती है, जो कतरनी विरूपण (shear deformation) और घूर्णी बंकन (rotational bending) को ध्यान में रखती है, जिससे यह मोटे बीम और उच्च-आवृत्ति उत्तेजनाओं के विश्लेषण के लिए उपयुक्त हो जाती है। अमूर्त प्रणाली इस प्रकार दी गई है:
$$\begin{aligned} \frac{d^2u}{dt^2} + \left(\alpha + \beta \|A^{1/2}u\|^2\right) Au + a_1 Bv &= f_1(t), \\ \frac{d^2v}{dt^2} + \gamma Av + \delta Cv + a_2 Bu &= f_2(t), \end{aligned}$$
जहाँ प्रारंभिक स्थितियाँ $u(0)=\phi_0, u'(0)=\phi_1, v(0)=\psi_0, v'(0)=\psi_1$ हैं। यहाँ, $A$ एक स्व-संयोजक (self-adjoint), धनात्मक-निश्चित ऑपरेटर है; $C$ एक सममित, परिबद्ध गैर-ऋणात्मक ऑपरेटर है; और $B$ एक बंद रैखिक ऑपरेटर है जो $A$ के सापेक्ष विशिष्ट सबऑर्डिनेशन स्थितियों को पूरा करता है। पद $\|A^{1/2}u\|^2$ किरचॉफ अरैखिकता (Kirchhoff nonlinearity) के एक अमूर्त अनुरूप का प्रतिनिधित्व करता है।

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

कार्यप्रणाली
यह शोध पत्र एक एकीकृत संख्यात्मक दृष्टिकोण प्रस्तावित करता जिसमें टेम्पोरल (समय) और स्पेशियल (स्थानिक) विवेकीकरण (discretization) शामिल है:

1. टेम्पोरल विवेकीकरण (Temporal Discretization):
   एक सममित तीन-परत अर्ध-विवेकी समय-चरण (semi-discrete time-stepping) योजना विकसित की गई है। अरैखिक पद का मूल्यांकन टेम्पोरल मिडपॉइंट ($t_k$) पर किया जाता है। यह विशिष्ट उपचार यह सुनिश्चित करता है कि मूल अरैखिक समस्या को प्रत्येक समय चरण पर एक रैखिक समस्या में बदला जा सके। फलस्वरूप, अज्ञात $u_k$ और $v_k$ की गणना समानांतर (parallel) रूप में की जा सकती है। यह योजना समय-चरण आकार $\tau$ के संबंध में द्वितीय-क्रम सटीक (second-order accurate) है।

2. स्पेशियल विवेकीकरण (Spatial Discretization):
   स्थानिक एक-आयामी मामले (विशिष्ट टिमोशेंको बीम प्रणाली) के लिए, एक लेजेंड्रे-गेलरकिन स्पेक्ट्रल विधि (Legendre–Galerkin spectral method) का उपयोग किया गया है। परीक्षण फलन (trial functions) शिफ्टेड लेजेंड्रे बहुपदों के अंतर से निर्मित किए गए हैं। यह चयन सुनिश्चित करता है कि परिणामी रैखिक प्रणाली विरल (sparse) हो। विशेष रूप से, गुणांक मैट्रिक्स एक सममित धनात्मक-निश्चित ट्राइएडीगल (tridiagonal) मैट्रिक्स है जिसमें एक "गैप" है (मुख्य विकर्ण और दूसरे उप- और सुपर-विकर्णों पर ही गैर-शून्य प्रविष्टियाँ हैं)।

• डिकपलिंग (Decoupling): इस विशिष्ट मैट्रिक्स संरचना के कारण, प्रणाली को दो स्वतंत्र ट्राइएडीगल उप-प्रणालियों (एक विषम-सूचकांक वाले अज्ञात के लिए और एक सम-सूचकांक वाले अज्ञात के लिए) में विभाजित किया जा सकता है।
• दक्षता (Efficiency): इन उप-प्रणालियों को थॉमस एल्गोरिदम या स्क्वायर-रूट-मुक्त चोलेस्की अपघटन (Cholesky decomposition) का उपयोग करके हल किया जा सकता है, जिससे प्रति समय चरण $O(N)$ कम्प्यूटेशनल जटिलता प्राप्त होती है, जहाँ $N$ आधार फलनों की संख्या है। यह साइन-आधारित सूत्रीकरण के विपरीत है, जिसमें घने युग्मित मैट्रिक्स के कारण $O(N^2)$ संचालन की आवश्यकता होती है।

प्रमुख योगदान और सैद्धांतिक परिणाम

• अभिसरण विश्लेषण (Convergence Analysis): लेखक प्रस्तावित योजना के अभिसरण को सिद्ध करते हैं। सटीक समाधान पर स्मूथनेस धारणाओं (विशेष रूप से, तृतीय-क्रम अवकलनीयता और डेरिवेटिव्स पर लिप्सचिट्ज़ स्थितियाँ) के तहत, अनुमानित समाधानों $z_{j,k}$ के लिए त्रुटि अनुमान इस प्रकार हैं:
  $$\max_{1 \le k \le m} \|A^{1/2}z_{j,k}\| \le C \tau^2, \quad \max_{1 \le k \le m} \left\| \frac{\Delta z_{j,k-1}}{\tau} \right\| \le C \tau^2$$
  यह समय में द्वितीय-क्रम सटीकता स्थापित करता है।
• एकीकृत परिबद्धता (Uniform Boundedness): प्रमाण का एक महत्वपूर्ण हिस्सा समाधान वेक्टर्स और उनके अंतर एनालॉग्स (जैसे, $A^{1/2}u_k$ और $A u_k$) की एक स्थानीय टेम्पोरल अंतराल पर यूनिफॉर्म बाउंडेडनेस स्थापित करना है। इसे ऊर्जा विधियों और डिस्क्रीट ग्रोनवाल-प्रकार के असमानताओं (Grönwall-type inequalities) का उपयोग करके प्राप्त किया गया है, जो इस गैर-triviality को दूर करता है जहाँ मानक ऊर्जा विधियाँ इस अरैखिक प्रणाली के लिए टेलिस्कोपिंग श्रृंखला रद्दीकरण (telescoping series cancellation) प्रदान नहीं करती हैं।
• स्पेक्ट्रल विधि के लिए त्रुटि अनुमान: शोध पत्र $L^2$ और यूनिफॉर्म नॉर्म्स दोनों में लेजेंड्रे-गेलरकिन सन्निकटन के लिए त्रुटि अनुमान प्रदान करता है, जो दिखाता है कि $C^{p+1}$ समाधानों के लिए त्रुटि $O(N^{-(p-1/2)})$ के रूप में घटती है।
• अमूर्त सेटिंग में विस्तार: 1D टिमोशेंको बीम के लिए व्युत्पन्न परिणामों को वापस अमूर्त हिल्बर्ट स्पेस सेटिंग में विस्तारित किया गया है, जो व्यापक श्रेणी के युग्मित हाइपरबोलिक प्रणालियों के लिए इस दृष्टिकोण को मान्य करता है।

संख्यात्मक परिणाम
प्रस्तावित एल्गोरिदम को पायथन (Python) में लागू किया गया और चार बेंचमार्क समस्याओं पर परीक्षण किया गया:

1. परीक्षण 1 और 2: ज्ञात सटीक विश्लेषणात्मक समाधानों वाली समस्याएं जिनमें दोलनी व्यवहार (oscillatory behavior) शामिल है। संख्यात्मक परिणामों ने प्रदर्शित किया कि आधार फलनों ($N$) की संख्या बढ़ाने से $L^2$ त्रुटि काफी कम हो गई, जिससे क्रमशः $10^{-6}$ और $10^{-4}$ के क्रम की सटीकता प्राप्त हुई। योजना ने समाधानों की गुणात्मक संरचना को सफलतापूर्वक पुनरुत्पादित किया।
2. परीक्षण 3: तेजी से बढ़ती आयाम (amplitude) वाली एक समस्या। परिणामों ने पुष्टि की कि सटीकता बनाए रखने के लिए टेम्पोरल ग्रिड और आधार फलनों की संख्या दोनों को परिष्कृत करना आवश्यक था।
3. परीक्षण 4: वेव पैकेट प्रारंभिक डेटा के साथ एक होमोजेनियस समस्या जिसमें कोई सटीक विश्लेषणात्मक समाधान नहीं है। एक एडेप्टिव रिफाइनमेंट रणनीति (अभिसरण तक टाइम स्टेप्स को दोगुना करना और आधार फलनों को बढ़ाना) का उपयोग किया गया। परिणाम स्थिर तरंग प्रसार (wave propagation) और परावर्तन दिखाए, जो भौतिक अपेक्षाओं के अनुरूप है, और त्रुटियां निर्धारित सहनशीलता ($10^{-4}$) से नीचे स्थिर हो गईं।

महत्व और दावे
लेखक दावा करते हैं कि यह कार्य अरैखिक टिमोशेंको प्रणाली के अमूर्त अनुरूप के लिए संख्यात्मक एल्गोरिदम और अनुमानित समाधानों के प्रथम डिज़ाइन प्रदान करके साहित्य के अंतराल को भरता है। इस कार्य का महत्व निम्नलिखित में निहित है:

• एकीकृत दृष्टिकोण: अमूर्त सूत्रीकरण इस पद्धति को संबंधित समस्याओं की एक श्रेणी में लागू करने की अनुमति देता है, जिसमें विभिन्न सीमा मान समस्याएं (Dirichlet, Neumann, Robin) और लिप्सचिट्ज़ डोमेन पर बहु-आयामी मामले शामिल हैं।
• कम्प्यूटेशनल दक्षता: सममित तीन-परत योजना और विशिष्ट लेजेंड्रे-गेलरकिन आधार चयन का संयोजन समानांतर गणना को सक्षम बनाता है और स्थानिक विवेकीकरण लागत को रैखिक जटिलता $O(N)$ तक कम करता है, जो इस विशिष्ट युग्मित प्रणाली के लिए अन्य स्पेक्ट्रल विधियों की तुलना में एक विशिष्ट लाभ है।
• सत्यापन: द्वितीय-क्रम अभिसरण और स्थिरता के सैद्धांतिक निष्कर्षों की पुष्टि बेंचमार्क समस्याओं पर संख्यात्मक प्रयोगों के माध्यम से की गई है, जो संयुक्त योजना की मजबूती और उच्च व्यावहारिक सटीकता को प्रदर्शित करती है।

यह शोध पत्र निष्कर्ष निकालता है कि प्रस्तावित विधि अरैखिक गतिशील टिमोशेंको बीम प्रणालियों और उनके अमूर्त सामान्यीकरणों के संख्यात्मक उपचार के लिए एक व्यवहार्य और कुशल उपकरण है।