2025-11-18T01:22:13.423549

Existence and numerical approximation of a one-dimensional Boussinesq system with variable coefficients on a finite interval

Grajales, Pizo
In this paper, we investigate the well-posedness of a nonlinear dispersive model with variable coefficients that describes the evolution of surface waves propagating through a one-dimensional shallow water channel of finite length with irregular bottom topography. To complement the theoretical analysis, we utilize the numerical solver developed by the authors in \cite{PizoMunoz} to approximate solutions of the model on a finite spatial interval, considering various parameter values and forms of the variable coefficients in the Boussinesq system under study. Additionally, we present preliminary numerical experiments addressing an inverse problem: the reconstruction of the initial wave elevation and fluid velocity from measurements taken at a final time. This is achieved by formulating an optimization problem in which the initial conditions are estimated as minimizers of a functional that quantifies the discrepancy between the observed final state and the numerical solution evolved from a trial initial state.
academic

एक-आयामी Boussinesq प्रणाली का अस्तित्व और संख्यात्मक सन्निकटन परिमित अंतराल पर परिवर्तनशील गुणांकों के साथ

मूल जानकारी

  • पेपर ID: 2510.10829
  • शीर्षक: Existence and numerical approximation of a one-dimensional Boussinesq system with variable coefficients on a finite interval
  • लेखक: Deissy Marcela Pizo, Juan Carlos Muñoz Grajales (Universidad del Valle, कोलंबिया)
  • वर्गीकरण: math.NA, cs.NA, math.AP
  • प्रकाशन समय: 14 अक्टूबर, 2025 (arXiv प्रीप्रिंट)
  • पेपर लिंक: https://arxiv.org/abs/2510.10829

सारांश

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

अनुसंधान पृष्ठभूमि और प्रेरणा

  1. समाधान की जाने वाली समस्याएं:
    • परिवर्तनशील गुणांकों के साथ एक-आयामी Boussinesq प्रणाली की गणितीय सुस्थिरता (अस्तित्व, अद्वितीयता, सतत निर्भरता) का अध्ययन
    • इस परिवर्तनशील गुणांक प्रणाली को हल करने के लिए प्रभावी संख्यात्मक विधियों का विकास
    • अंतिम समय के अवलोकन डेटा के आधार पर प्रारंभिक शर्तों के पुनर्निर्माण की व्युत्क्रम समस्या की खोज
  2. समस्या की महत्ता:
    • Boussinesq प्रणाली उथले जल तरंग प्रसार का वर्णन करने वाली एक महत्वपूर्ण गणितीय मॉडल है, जिसका समुद्री इंजीनियरिंग, तटीय सुरक्षा आदि क्षेत्रों में महत्वपूर्ण अनुप्रयोग है
    • वास्तविक जल निकायों में अक्सर असमान तलीय स्थलाकृति होती है, इसलिए परिवर्तनशील गुणांक मॉडल वास्तविक स्थिति के अधिक करीब है
    • व्युत्क्रम समस्याएं तरंग पूर्वानुमान, समुद्री निगरानी आदि में महत्वपूर्ण व्यावहारिक मूल्य रखती हैं
  3. मौजूदा विधियों की सीमाएं:
    • अधिकांश अनुसंधान स्थिर गुणांक Boussinesq प्रणाली पर केंद्रित है
    • परिवर्तनशील गुणांक स्थिति के सैद्धांतिक विश्लेषण में अपेक्षाकृत कमी है
    • इस प्रकार की फैलाववादी प्रणालियों के लिए व्युत्क्रम समस्या अनुसंधान अभी भी रिक्त है
  4. अनुसंधान प्रेरणा:
    • परिवर्तनशील गुणांक Boussinesq प्रणाली के सैद्धांतिक विश्लेषण में अंतराल को भरना
    • वास्तविक समुद्री इंजीनियरिंग अनुप्रयोगों के लिए अधिक सटीक गणितीय मॉडल और संख्यात्मक उपकरण प्रदान करना
    • फैलाववादी तरंग गतिशील प्रणालियों की व्युत्क्रम समस्या अनुसंधान में नई दिशा खोलना

मूल योगदान

  1. परिवर्तनशील गुणांक एक-आयामी Boussinesq प्रणाली के लिए सुस्थिरता सिद्धांत की स्थापना: समाधान के अस्तित्व, अद्वितीयता और प्रारंभिक डेटा पर सतत निर्भरता को सिद्ध किया
  2. Green फलन तकनीक के आधार पर सैद्धांतिक विश्लेषण विधि का विकास: समस्या को अरैखिक संचालक की निश्चित बिंदु समस्या में परिवर्तित किया
  3. ऊर्जा संरक्षण कानून प्रदान किया: प्रणाली के ऊर्जा कार्यात्मक को समय पर संरक्षित सिद्ध किया
  4. कुशल परिमित तत्व संख्यात्मक समाधानकर्ता का कार्यान्वयन: FEniCS लाइब्रेरी के साथ स्थिर संख्यात्मक योजना विकसित की
  5. इस प्रकार की फैलाववादी प्रणाली की व्युत्क्रम समस्या का पहली बार अध्ययन: परिवर्तनशील अनुकूलन के आधार पर प्रारंभिक शर्तों के पुनर्निर्माण की विधि प्रस्तावित की

विधि विवरण

कार्य परिभाषा

निम्नलिखित परिवर्तनशील गुणांक Boussinesq प्रणाली का अध्ययन करें:

N_t + [(1 + \alpha c^2(\xi)N)V]_\xi - \frac{\beta}{6}N_{\xi\xi t} = 0 \\ V_t + (c(\xi)N)_\xi + \frac{1}{2}\alpha(c^2V^2)_\xi - \frac{\beta}{6}V_{\xi\xi t} = 0 \end{cases}$$ जहां $(\xi,t) \in [0,L] \times [0,T]$, प्रारंभिक शर्तों के साथ: $$N(\xi,0) = N_0(\xi), \quad V(\xi,0) = V_0(\xi)$$ और Dirichlet सीमा शर्तें: $$N(0,t) = N(L,t) = 0, \quad V(0,t) = V(L,t) = 0$$ ### सैद्धांतिक विश्लेषण विधि **1. Green फलन तकनीक** प्रणाली को निम्नलिखित रूप में पुनः लिखें: $$\left(I - \frac{\beta}{6}\partial_\xi^2\right)N_t = -\partial_\xi[(1 + \alpha c^2(\xi)N)V]$$ $$\left(I - \frac{\beta}{6}\partial_\xi^2\right)V_t = -\partial_\xi\left[c(\xi)N + \frac{1}{2}\alpha c^2V^2\right]$$ संचालक $P = I - \frac{\beta}{6}\partial_\xi^2$ के Green फलन का उपयोग करें: $$G(\xi,s) = \frac{1}{2\sqrt{\beta/6}} \frac{\cosh\left(\frac{L-|\xi-s|}{\sqrt{\beta/6}}\right) - \cosh\left(\frac{L-(\xi+s)}{\sqrt{\beta/6}}\right)}{\sinh\left(\frac{L}{\sqrt{\beta/6}}\right)}$$ **2. निश्चित बिंदु प्रमेय** समस्या को समाकल समीकरण समूह में परिवर्तित करें, अरैखिक संचालक $A$ को परिभाषित करें, Banach निश्चित बिंदु प्रमेय के माध्यम से समाधान की अस्तित्व और अद्वितीयता को सिद्ध करें। ### संख्यात्मक विधि **परिमित तत्व विवेकीकरण** Galerkin विधि का उपयोग करें, कमजोर रूप: $$\partial_t\left(\langle N,v_1\rangle_{L^2} + \frac{\beta}{6}\langle N_\xi, v_{1\xi}\rangle_{L^2}\right) = \langle[(1+\alpha c^2N)V], v_{1\xi}\rangle_{L^2}$$ **समय विवेकीकरण** समय एकीकरण के लिए $\theta$-विधि ($\theta = 1/2$) का उपयोग करें, संख्यात्मक स्थिरता सुनिश्चित करें। ### व्युत्क्रम समस्या विधि लक्ष्य कार्यात्मक को परिभाषित करें: $$J(N,V) := \frac{1}{2}\int_0^L \left(|N(\xi,T;N_0,V_0) - N_T(\xi)|^2 + |V(\xi,T;N_0,V_0) - V_T(\xi)|^2\right)d\xi$$ इस कार्यात्मक को न्यूनतम करने के लिए L-BFGS-B एल्गोरिदम का उपयोग करें, Dolfin-Adjoint लाइब्रेरी के साथ स्वचालित अवकलन को लागू करें। ## प्रयोग सेटअप ### संख्यात्मक पैरामीटर - **स्थान डोमेन**: $[-20, 40]$ - **जाल बिंदु संख्या**: 3000 स्थान बिंदु - **स्थान संकल्प**: $\Delta\xi \approx 0.02$ - **समय चरण**: $\Delta t = 8/3000 \approx 2.7 \times 10^{-3}$ - **मॉडल पैरामीटर**: $\alpha = \beta = 0.1$ ### प्रारंभिक शर्तें स्थानीयकृत गाऊसी समोच्च का उपयोग करें: $$N(\xi,0) = V(\xi,0) = e^{-(\xi-18)^2}$$ ### परिवर्तनशील गुणांक फलन **स्थिति 1**: चिकनी दोलन-गाऊसी विक्षोभ $$c(\xi) = 1 + 0.3\sin\left(\frac{\pi}{5}\xi\right) + 0.6e^{-(\xi-8)^2}$$ **स्थिति 2**: खंडशः स्थिर फलन (स्तरीकृत माध्यम का अनुकरण) $$c(\xi) = \begin{cases} 0.8, & \xi < 20 \\ 1.5, & 20 < \xi \leq 22 \\ 2.0, & 22 < \xi \leq 27 \\ \vdots \end{cases}$$ ## प्रयोग परिणाम ### सैद्धांतिक परिणाम **प्रमेय 3.1 (स्थानीय सुस्थिरता)**: $T_0 > 0$ मौजूद है, जैसे कि प्रणाली $L^2_{T_0} \times L^2_{T_0}$ में अद्वितीय समाधान $(N,V)$ रखती है, और अनुमान को संतुष्ट करती है: $$\|(N,V)\|_{L^2_{T_0} \times L^2_{T_0}} \leq \|(N_0,V_0)\|_{H^1_0 \times H^1_0} T_0^{1/2} \exp\left(\frac{L^{1/2}}{(\beta/6)^{1/2}}DT_0\right)$$ **प्रमेय 3.2 (ऊर्जा संरक्षण)**: $$\|V\|_{L^2}^2 + \alpha\|c|N|^{1/2}V\|_{L^2}^2 + \||c|^{1/2}N\|_{L^2}^2 = \text{स्थिरांक}$$ ### संख्यात्मक प्रयोग परिणाम 1. **अग्रगामी समस्या समाधान**: - परिवर्तनशील गुणांक माध्यम में तरंग प्रसार का सफलतापूर्वक अनुकरण किया - तलीय स्थलाकृति परिवर्तन के कारण तरंग गति परिवर्तन, आयाम मॉड्यूलेशन और फैलाववादी प्रभाव देखे गए - संख्यात्मक विधि अच्छी स्थिरता और सटीकता प्रदर्शित करती है 2. **व्युत्क्रम समस्या पुनर्निर्माण**: - L-BFGS-B एल्गोरिदम 4 पुनरावृत्तियों में प्रारंभिक शर्तों के पुनर्निर्माण में महत्वपूर्ण सुधार किया - पुनर्निर्मित तरंग ऊंचाई और द्रव वेग क्रमशः लक्ष्य समोच्च के करीब आते हैं - परिवर्तनशील सहअनुकूलन रणनीति की प्रभावशीलता को सिद्ध किया ### संख्यात्मक सत्यापन - ऊर्जा संरक्षण कानून संख्यात्मक गणना में अच्छी तरह से संरक्षित है - विभिन्न जाल संकल्पों के तहत परिणामों की सामंजस्यता संख्यात्मक विधि के अभिसरण को सत्यापित करती है ## संबंधित कार्य 1. **Boussinesq प्रणाली सिद्धांत**: - Bona & Chen (1998, 2002, 2004): स्थिर गुणांक स्थिति के लिए मूल सिद्धांत की स्थापना - यह पेपर परिवर्तनशील गुणांक स्थिति तक विस्तारित करता है, सैद्धांतिक अंतराल को भरता है 2. **संख्यात्मक विधियां**: - फैलाववादी तरंग समीकरणों में परिमित तत्व विधि का अनुप्रयोग - यह पेपर FEniCS लाइब्रेरी के साथ कुशल समाधानकर्ता को लागू करता है 3. **व्युत्क्रम समस्याएं**: - तरंग समीकरण व्युत्क्रम समस्या की पारंपरिक विधियां - यह पेपर पहली बार Boussinesq प्रणाली व्युत्क्रम समस्या पर परिवर्तनशील अनुकूलन लागू करता है ## निष्कर्ष और चर्चा ### मुख्य निष्कर्ष 1. परिवर्तनशील गुणांक एक-आयामी Boussinesq प्रणाली के गणितीय सिद्धांत की नींव सफलतापूर्वक स्थापित की 2. प्रभावी संख्यात्मक समाधान विधि विकसित की 3. इस प्रकार की फैलाववादी प्रणाली व्युत्क्रम समस्या अनुसंधान में नई दिशा खोली 4. वास्तविक समुद्री इंजीनियरिंग अनुप्रयोग के लिए सैद्धांतिक और कम्प्यूटेशनल उपकरण प्रदान किए ### सीमाएं 1. **सैद्धांतिक पहलू**: - केवल स्थानीय सुस्थिरता सिद्ध की गई है, वैश्विक अस्तित्व अनुसंधान की प्रतीक्षा में है - व्युत्क्रम समस्या की अद्वितीयता और स्थिरता सिद्धांत विश्लेषण अनुपस्थित है 2. **संख्यात्मक पहलू**: - व्युत्क्रम समस्या प्रयोग केवल आदर्श शोर-मुक्त स्थिति पर विचार करते हैं - अवलोकन शोर संवेदनशीलता विश्लेषण अनुपस्थित है 3. **व्यावहारिक पहलू**: - एक-आयामी स्थिति तक सीमित, वास्तविक अनुप्रयोग के लिए दो-आयामी तक विस्तार की आवश्यकता है - परिवर्तनशील गुणांक फलन का चयन अपेक्षाकृत सरल है ### भविष्य की दिशाएं 1. **सैद्धांतिक विकास**: - वैश्विक सुस्थिरता सिद्धांत - व्युत्क्रम समस्या की अद्वितीयता और स्थिरता विश्लेषण - नियमितकरण सिद्धांत 2. **संख्यात्मक सुधार**: - शोर डेटा के लिए मजबूत पुनर्निर्माण विधि - उच्च-आयामी स्थितियों का विस्तार - मशीन लर्निंग तकनीकों का एकीकरण 3. **व्यावहारिक अनुप्रयोग**: - वास्तविक समुद्री डेटा के साथ सत्यापन - समुद्री अवलोकन प्रणालियों के साथ एकीकरण ## गहन मूल्यांकन ### लाभ 1. **सैद्धांतिक योगदान महत्वपूर्ण**: परिवर्तनशील गुणांक Boussinesq प्रणाली की सुस्थिरता का पहली बार व्यवस्थित अध्ययन 2. **विधि नवाचार**: परिवर्तनशील गुणांक समस्या को संभालने के लिए Green फलन तकनीक का चतुर उपयोग 3. **संख्यात्मक कार्यान्वयन पूर्ण**: संपूर्ण संख्यात्मक समाधान ढांचा प्रदान किया 4. **समस्या नवीन**: फैलाववादी प्रणाली व्युत्क्रम समस्या अनुसंधान में नया क्षेत्र खोला 5. **लेखन स्पष्ट**: गणितीय व्युत्पत्ति कठोर, अभिव्यक्ति स्पष्ट ### कमियां 1. **सैद्धांतिक गहराई सीमित**: केवल स्थानीय सिद्धांत स्थापित किया, वैश्विक विश्लेषण अनुपस्थित 2. **प्रयोग सत्यापन अपर्याप्त**: व्युत्क्रम समस्या प्रयोग अपेक्षाकृत सरल, जटिल परिदृश्य सत्यापन अनुपस्थित 3. **व्यावहारिकता सत्यापन प्रतीक्षा में**: अधिक वास्तविक डेटा सत्यापन की आवश्यकता 4. **कम्प्यूटेशनल जटिलता विश्लेषण अनुपस्थित**: एल्गोरिदम की कम्प्यूटेशनल जटिलता का विश्लेषण नहीं किया ### प्रभाव 1. **शैक्षणिक मूल्य**: परिवर्तनशील गुणांक फैलाववादी प्रणाली सिद्धांत विकास के लिए आधार स्थापित करता है 2. **अनुप्रयोग संभावनाएं**: समुद्री इंजीनियरिंग, तरंग पूर्वानुमान आदि क्षेत्रों में संभावित अनुप्रयोग मूल्य 3. **पद्धति योगदान**: Green फलन विधि और परिवर्तनशील अनुकूलन का संयोजन संबंधित समस्याओं के लिए नई सोच प्रदान करता है ### लागू परिदृश्य 1. उथले जल तरंग प्रसार मॉडलिंग 2. तटीय इंजीनियरिंग में तरंग विश्लेषण 3. समुद्री पर्यावरण निगरानी में व्युत्क्रम समस्याएं 4. संख्यात्मक समुद्र विज्ञान में डेटा सहअनुकूलन ## संदर्भ पेपर 20 महत्वपूर्ण संदर्भों का हवाला देता है, जिसमें Boussinesq प्रणाली सिद्धांत, संख्यात्मक विधियां, कार्यात्मक विश्लेषण आदि संबंधित क्षेत्रों के शास्त्रीय कार्य शामिल हैं, जो अनुसंधान के लिए एक मजबूत सैद्धांतिक आधार प्रदान करते हैं।