The purpose of this paper is to propose a new algorithm for obtaining approximate solutions to the Burgers' equation (BE). Integration in time by a quadratic B-spline collocation method is shown. To the best of our knowledge, B-splines have not previously been used to integrate partial differential equations in both time and space. First, the BE is integrated using quadratic B-spline functions in time, and then the time-integrated BE is further solved in space via the cubic B-spline collocation method. The resulting recursive algebraic equation is used to obtain both shock wave and front propagation solutions of the BE, demonstrating the effectiveness of the space--time B-spline collocation method.
- पेपर ID: 2510.09408
- शीर्षक: A time-space B-spline integrator for the Burgers' equation
- लेखक: इद्रिस डैग, सर्कन उगुरलुओगलु, निहत अदर (एस्किशेहिर ओस्मानगाजी विश्वविद्यालय, तुर्की)
- वर्गीकरण: math.NA cs.NA
- प्रकाशन समय: 25 अक्टूबर 10 (arXiv प्रीप्रिंट)
- पेपर लिंक: https://arxiv.org/abs/2510.09408
यह पेपर बर्गर्स समीकरण के अनुमानित समाधान के लिए एक नया एल्गोरिदम प्रस्तावित करता है। यह विधि समय समाकलन के लिए द्विघात B-स्पलाइन कोलोकेशन विधि का उपयोग करती है। लेखकों के ज्ञान के अनुसार, यह पहली बार है कि B-स्पलाइन का उपयोग समय और स्थान दोनों आयामों में आंशिक अवकल समीकरणों को एकीकृत करने के लिए किया गया है। सबसे पहले समय दिशा में बर्गर्स समीकरण को द्विघात B-स्पलाइन फलनों का उपयोग करके एकीकृत किया जाता है, फिर स्थान दिशा में घन B-स्पलाइन कोलोकेशन विधि के माध्यम से समय-एकीकृत समीकरण को हल किया जाता है। प्राप्त पुनरावर्ती बीजगणितीय समीकरणों का उपयोग बर्गर्स समीकरण के आघात तरंग और अग्रभाग प्रसार समाधान प्राप्त करने के लिए किया जाता है, जो समय-स्थान B-स्पलाइन कोलोकेशन विधि की प्रभावशीलता को प्रदर्शित करता है।
- बर्गर्स समीकरण का महत्व: बर्गर्स समीकरण एक अर्ध-रैखिक परवलयिक आंशिक अवकल समीकरण है, जो इंजीनियरिंग विज्ञान में अशांति, गैस गतिविज्ञान, ऊष्मा चालन, लोच, यातायात प्रवाह और ध्वनि तरंगों जैसी घटनाओं को मॉडल करने में व्यापक रूप से लागू होता है। यह संवहन और प्रसार के पारस्परिक क्रिया का वर्णन करने के लिए गणितीय अभिव्यक्ति प्रदान करता है।
- संख्यात्मक समाधान की आवश्यकता: अरैखिक आंशिक अवकल समीकरणों के पास आमतौर पर विश्लेषणात्मक समाधान नहीं होते हैं, विशेष रूप से जब श्यानता स्थिरांक υ छोटा हो। श्रृंखला समाधान धीरे-धीरे अभिसरित होते हैं और व्यावहारिक उपयोगिता कम होती है। छोटे υ मान आघात तरंग के गठन की ओर ले जाते हैं, जिससे विश्लेषणात्मक और संख्यात्मक दोनों विधियों के लिए उचित समाधान प्रदान करना कठिन हो जाता है।
- मौजूदा विधियों की सीमाएं: पारंपरिक समय इंटीग्रेटर में क्रैंक-निकोलसन प्रारूप, रंज-कुट्टा विधि, ज्यामितीय इंटीग्रेटर, घातीय इंटीग्रेटर आदि शामिल हैं, लेकिन ये विधियां बर्गर्स समीकरण की आघात तरंग और अग्रभाग प्रसार समस्याओं को संभालते समय सटीकता और स्थिरता के संदर्भ में चुनौतियों का सामना करती हैं।
इस पेपर की मूल प्रेरणा एक नई संख्यात्मक विधि विकसित करना है जो:
- समय और स्थान आयामों में B-स्पलाइन का उपयोग करके एकीकृत कर सके
- छोटी श्यानता मान स्थितियों में आघात तरंग समस्याओं को प्रभावी ढंग से संभाल सके
- उच्च सटीकता के संख्यात्मक समाधान प्रदान कर सके
- अग्रणी विधि: पहली बार बर्गर्स समीकरण को हल करने के लिए समय और स्थान इंटीग्रेटर के रूप में B-स्पलाइन का एक साथ उपयोग करने का प्रस्ताव
- दोहरी B-स्पलाइन रणनीति: समय दिशा में द्विघात B-स्पलाइन और स्थान दिशा में घन B-स्पलाइन के नवीन संयोजन
- आघात तरंग और अग्रभाग प्रसार सिमुलेशन: बर्गर्स समीकरण के आघात तरंग समाधान और यात्रा तरंग अग्रभाग समाधान का सफल सिमुलेशन
- पुनरावर्ती एल्गोरिदम ढांचा: B-स्पलाइन कोलोकेशन के आधार पर पुनरावर्ती बीजगणितीय समीकरण समाधान प्रणाली की स्थापना
एक-आयामी समय-निर्भर बर्गर्स समीकरण को हल करना:
ut+uux−νuxx=0,x∈[a,b],t∈(0,T]
प्रारंभिक शर्त के साथ: u(x,0)=f(x)
सीमा शर्तें: u(a,t)=u0,ux(a,t)=u0′,u(b,t)=um,ux(b,t)=um′
जहां ν श्यानता स्थिरांक है।
- जाल विभाजन: समय क्षेत्र [c,d] को जाल बिंदुओं में समान रूप से विभाजित करें c=t0<t1<⋯<tk=d, जहां tj=t0+jΔt
- द्विघात B-स्पलाइन आधार फलन: Bj(t)∈C1[a,b] को परिभाषित करें, j=−1,…,n:
(t_{j+2}-t)^2 - 3(t_{j+1}-t)^2 + 3(t_j-t)^2, & t_{j-1} \leq t \leq t_j \\
(t_{j+2}-t)^2 - 3(t_{j+1}-t)^2, & t_j \leq t \leq t_{j+1} \\
(t_{j+2}-t)^2, & t_{j+1} \leq t \leq t_{j+2} \\
0, & \text{अन्यथा}
\end{cases}$$
3. **अनुमानित समाधान प्रतिनिधित्व**:
$$U(x,t) = \sum_{j=-1}^n \delta(x,t_j)B_j(t)$$
#### स्थान विवेकीकरण (घन B-स्पलाइन)
1. **स्थान जाल**: स्थान क्षेत्र $[a,b]$ को समान रूप से विभाजित करें, $x_i = x_0 + ih$, $i = 0,\ldots,m$
2. **घन B-स्पलाइन आधार फलन**: $B_i(x) \in C^2[a,b]$, $i = -1,\ldots,m$
3. **पैरामीटर विस्तार**:
$$\delta(x,t^j) = \sum_{i=-1}^{m+1} \sigma_i^j B_i(x)$$
### तकनीकी नवाचार बिंदु
1. **स्तरीय एकीकरण रणनीति**: पहले समय दिशा में द्विघात B-स्पलाइन का उपयोग करके एकीकृत करें, फिर स्थान दिशा में घन B-स्पलाइन का उपयोग करके हल करें
2. **कोलोकेशन बिंदु विधि**: जाल बिंदुओं पर अवकल समीकरण को संतुष्ट करने के लिए बाध्य करें
3. **रैखिकीकरण प्रक्रिया**: पिछले समय चरण के मान का चयन करके अरैखिक पदों को रैखिक करें
4. **पुनरावर्ती परिशोधन**: प्रत्येक समय चरण की सटीकता में सुधार के लिए पुनरावर्ती प्रक्रिया अपनाएं
## प्रयोगात्मक सेटअप
### परीक्षण समस्याएं
#### समस्या 1: आघात तरंग समाधान
विश्लेषणात्मक समाधान:
$$u(x,t) = \frac{x/t}{1 + \sqrt{t/t_0}\exp(x^2/(4\nu t))}, \quad t \geq 1, \quad 0 \leq x \leq 1$$
जहां $t_0 = \exp(1/(8\nu))$
#### समस्या 2: यात्रा तरंग अग्रभाग समाधान
विश्लेषणात्मक समाधान:
$$u(x,t) = \frac{\alpha + \mu + (\mu - \alpha)\exp(\eta)}{1 + \exp(\eta)}$$
जहां $\eta = \frac{\alpha(x - \mu t - \gamma)}{\nu}$, $\alpha = 0.4$, $\mu = 0.6$, $\gamma = 0.125$
### मूल्यांकन संकेतक
- **$L_\infty$ मानदंड**: $|U - U_N|_\infty = \max|U_i^j - (U_N)_i^j|$
- **$L_2$ मानदंड**: $|U - U_N|_2 = \sqrt{\sum_i |U_i^j - (U_N)_i^j|^2}$
- **शिखर मान**: आघात तरंग का अधिकतम मान
### कार्यान्वयन विवरण
- समय-स्थान चरण: $h = \Delta t = 0.01, 0.005, 0.001$
- श्यानता पैरामीटर: $\nu = 0.01, 0.005$
- समय श्रेणी: $t = 3.25$ (आघात तरंग समस्या), $t = 1.2$ (अग्रभाग तरंग समस्या)
## प्रयोगात्मक परिणाम
### मुख्य परिणाम
#### आघात तरंग समस्या परिणाम
| $\nu=0.01$ | $h=\Delta t$ | $L_\infty$ | $L_2$ | शिखर मान |
|-------------|--------------|------------|-------|----------|
| | 0.01 | 0.0498 | 0.0211| 0.2233 |
| | 0.005 | 0.0508 | 0.0207| 0.2225 |
| | 0.001 | 0.0516 | 0.0206| 0.2203 |
| $\nu=0.005$ | $L_\infty$ | $L_2$ | शिखर मान |
|--------------|------------|-------|----------|
| | 0.0582 | 0.0166| 0.2503 |
| | 0.0571 | 0.0159| 0.2480 |
| | 0.0560 | 0.0155| 0.2467 |
#### अग्रभाग तरंग समस्या परिणाम
| $\nu=0.01$ | $h=\Delta t$ | $L_\infty$ | $L_2$ |
|-------------|--------------|------------|-------|
| | 0.01 | 0.4875 | 0.1138|
| | 0.005 | 0.2964 | 0.0638|
| | 0.001 | 0.1514 | 0.0354|
| $\nu=0.005$ | $L_\infty$ | $L_2$ |
|--------------|------------|-------|
| | 0.8841 | 0.2532|
| | 0.6906 | 0.1329|
| | 0.2427 | 0.0367|
### प्रयोगात्मक निष्कर्ष
1. **सटीकता और जाल संबंध**: छोटे समय-स्थान चरणों का उपयोग करके अधिक सुचारु समाधान और बेहतर सटीकता प्राप्त की जा सकती है
2. **श्यानता प्रभाव**: छोटी श्यानता स्थिरांक तीव्र तरंग रूप और उच्च त्रुटि की ओर ले जाती है
3. **त्रुटि वितरण**: अधिकतम त्रुटि आमतौर पर आघात तरंग शिखर या तरंग अग्रभाग क्षेत्र के पास होती है
4. **स्थिरता**: विधि विभिन्न पैरामीटर सेटिंग्स के तहत अच्छी संख्यात्मक स्थिरता बनाए रखती है
## संबंधित कार्य
इस क्षेत्र की मुख्य अनुसंधान दिशाएं शामिल हैं:
1. **पारंपरिक संख्यात्मक विधियां**: परिमित अंतर, परिमित तत्व, वर्णक्रम विधि आदि
2. **B-स्पलाइन विधियां**: मुख्य रूप से स्थान विवेकीकरण के लिए उपयोग की जाती हैं, समय एकीकरण के लिए शायद ही कभी
3. **कोलोकेशन विधियां**: गैलेरकिन विधि, कोलोकेशन विधि, जालरहित विधि आदि
4. **विशेष समाधानकर्ता**: बर्गर्स समीकरण के लिए विशेष संख्यात्मक विधियां
इस पेपर का नवाचार समय और स्थान एकीकरण के लिए B-स्पलाइन को पहली बार एक साथ लागू करना है, जो इस क्षेत्र में एक अंतराल को भरता है।
## निष्कर्ष और चर्चा
### मुख्य निष्कर्ष
1. प्रस्तावित समय-स्थान B-स्पलाइन कोलोकेशन विधि बर्गर्स समीकरण को प्रभावी ढंग से हल कर सकती है
2. विधि आघात तरंग और अग्रभाग प्रसार समस्याओं पर अच्छी सटीकता और स्थिरता प्रदर्शित करती है
3. समय-स्थान चरणों को समायोजित करके समाधान की सटीकता और सुचारुता को नियंत्रित किया जा सकता है
### सीमाएं
1. **कम्प्यूटेशनल जटिलता**: प्रत्येक समय चरण पर रैखिक समीकरण प्रणाली को हल करने की आवश्यकता है
2. **अरैखिकता प्रक्रिया**: रैखिकीकरण और पुनरावर्ती प्रक्रिया की आवश्यकता है
3. **पैरामीटर संवेदनशीलता**: श्यानता पैरामीटर और जाल पैरामीटर के प्रति अपेक्षाकृत संवेदनशील
4. **प्रयोज्यता श्रेणी**: मुख्य रूप से एक-आयामी बर्गर्स समीकरण के लिए
### भविष्य की दिशाएं
लेखकों ने भविष्य के कार्य पर ध्यान केंद्रित करने का प्रस्ताव दिया है:
- सटीकता बढ़ाने के लिए उच्च-क्रम B-स्पलाइन का उपयोग
- बहु-आयामी बर्गर्स समीकरण तक विस्तार
- एल्गोरिदम दक्षता का अनुकूलन
## गहन मूल्यांकन
### लाभ
1. **विधि नवाचार**: समय-स्थान B-स्पलाइन एकीकरण विधि पहली बार प्रस्तावित की गई है, मौलिकता है
2. **सैद्धांतिक पूर्णता**: गणितीय व्युत्पत्ति से संख्यात्मक कार्यान्वयन तक काफी पूर्ण है
3. **प्रयोगात्मक पर्याप्तता**: आघात तरंग और अग्रभाग प्रसार दोनों विशिष्ट समस्याओं का परीक्षण किया गया है
4. **सटीकता विश्लेषण**: विस्तृत त्रुटि विश्लेषण और अभिसरण अध्ययन प्रदान किए गए हैं
### कमियां
1. **अपर्याप्त सैद्धांतिक विश्लेषण**: अभिसरण और स्थिरता का कठोर सैद्धांतिक प्रमाण अभाव है
2. **सीमित तुलनात्मक प्रयोग**: अन्य उन्नत विधियों के साथ पर्याप्त तुलना नहीं की गई है
3. **कम्प्यूटेशनल दक्षता विश्लेषण अभाव**: कम्प्यूटेशन समय और मेमोरी उपयोग विश्लेषण प्रदान नहीं किया गया है
4. **पैरामीटर चयन मार्गदर्शन अपर्याप्त**: B-स्पलाइन क्रम चयन के लिए सैद्धांतिक मार्गदर्शन अभाव है
### प्रभाव
1. **शैक्षणिक योगदान**: PDE संख्यात्मक समाधान के लिए नई सोच प्रदान करता है
2. **व्यावहारिक मूल्य**: विधि सरल और कार्यान्वयन में आसान है, इंजीनियरिंग अनुप्रयोग के लिए उपयुक्त है
3. **विस्तारशीलता**: ढांचा अन्य अरैखिक PDE तक विस्तारित किया जा सकता है
### प्रयोज्य परिदृश्य
यह विधि विशेष रूप से उपयुक्त है:
- बर्गर्स समीकरण के उच्च-सटीकता समाधान की आवश्यकता वाले वैज्ञानिक कम्प्यूटिंग के लिए
- आघात तरंग और अग्रभाग प्रसार घटनाओं के संख्यात्मक सिमुलेशन के लिए
- कम्प्यूटेशनल सटीकता के लिए उच्च आवश्यकता वाले लेकिन दक्षता के लिए मध्यम आवश्यकता वाले अनुप्रयोगों के लिए
## संदर्भ
पेपर ने 16 संबंधित संदर्भों का हवाला दिया है, जिसमें बर्गर्स समीकरण की संख्यात्मक विधियां, B-स्पलाइन सिद्धांत, कोलोकेशन विधि आदि संबंधित क्षेत्रों के महत्वपूर्ण कार्य शामिल हैं, जो इस अनुसंधान के लिए एक मजबूत सैद्धांतिक आधार प्रदान करते हैं।