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

একটি সীমিত ব্যবধিতে পরিবর্তনশীল সহগ সহ এক-মাত্রিক বুসিনেস্ক সিস্টেমের অস্তিত্ব এবং সংখ্যাসূচক আনুমানিকতা

মৌলিক তথ্য

  • পেপার আইডি: 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
  • প্রকাশনার সময়: অক্টোবর ১৪, ২০২৫ (arXiv প্রাক-প্রিন্ট)
  • পেপার লিঙ্ক: https://arxiv.org/abs/2510.10829

সারসংক্ষেপ

এই পেপারটি একটি অরৈখিক বিচ্ছুরণ মডেলের সুস্থতার সমস্যা অধ্যয়ন করে যা অনিয়মিত তলদেশের ভূগোল সহ এক-মাত্রিক সীমিত দৈর্ঘ্যের অগভীর জল চ্যানেলে পৃষ্ঠ তরঙ্গের প্রচার বিবর্তন বর্ণনা করে। তাত্ত্বিক বিশ্লেষণ পরিপূরক করার জন্য, লেখকরা তাদের পূর্ববর্তী কাজে উন্নত সংখ্যাসূচক সমাধানকারী (সীমিত উপাদান পদ্ধতি) ব্যবহার করে সীমাবদ্ধ স্থান ব্যবধিতে মডেল সমাধানের আনুমানিক মূল্যায়ন করেন, বুসিনেস্ক সিস্টেমে পরিবর্তনশীল সহগের বিভিন্ন প্যারামিটার মান এবং ফাংশন ফর্ম বিবেচনা করেন। অধিকন্তু, এই পেপারটি চূড়ান্ত সময়ের পরিমাপ ডেটা থেকে প্রাথমিক তরঙ্গ উচ্চতা এবং তরল বেগ পুনর্নির্মাণের জন্য একটি বিপরীত সমস্যা সমাধানের জন্য প্রাথমিক সংখ্যাসূচক পরীক্ষা প্রস্তাব করে। এটি একটি অপ্টিমাইজেশন সমস্যা নির্মাণের মাধ্যমে অর্জিত হয়, যেখানে প্রাথমিক শর্তগুলি একটি উদ্দেশ্য ফাংশনাল ন্যূনতমকরণের সমাধান হিসাবে অনুমান করা হয়, যা পর্যবেক্ষিত চূড়ান্ত অবস্থা এবং পরীক্ষামূলক প্রাথমিক অবস্থা থেকে বিকশিত সংখ্যাসূচক সমাধানের মধ্যে পার্থক্য পরিমাপ করে।

গবেষণার পটভূমি এবং প্রেরণা

  1. সমাধান করার সমস্যা:
    • পরিবর্তনশীল সহগ সহ এক-মাত্রিক বুসিনেস্ক সিস্টেমের গাণিতিক সুস্থতা (অস্তিত্ব, অনন্যতা, ক্রমাগত নির্ভরতা) অধ্যয়ন করা
    • এই পরিবর্তনশীল সহগ সিস্টেম সমাধানের জন্য কার্যকর সংখ্যাসূচক পদ্ধতি বিকাশ করা
    • চূড়ান্ত সময়ের পর্যবেক্ষণ ডেটার উপর ভিত্তি করে প্রাথমিক শর্ত পুনর্নির্মাণের বিপরীত সমস্যা অন্বেষণ করা
  2. সমস্যার গুরুত্ব:
    • বুসিনেস্ক সিস্টেম অগভীর জল তরঙ্গ প্রচার বর্ণনার জন্য একটি গুরুত্বপূর্ণ গাণিতিক মডেল, সামুদ্রিক প্রকৌশল, উপকূলীয় সুরক্ষা এবং অন্যান্য ক্ষেত্রে গুরুত্বপূর্ণ প্রয়োগ রয়েছে
    • বাস্তব জল সংস্থাগুলি প্রায়শই অসমান তলদেশের ভূগোল রয়েছে, তাই পরিবর্তনশীল সহগ মডেল বাস্তব পরিস্থিতির কাছাকাছি
    • বিপরীত সমস্যা তরঙ্গ পূর্বাভাস, সামুদ্রিক পর্যবেক্ষণ এবং অন্যান্য দিকগুলিতে গুরুত্বপূর্ণ ব্যবহারিক মূল্য রয়েছে
  3. বিদ্যমান পদ্ধতির সীমাবদ্ধতা:
    • বেশিরভাগ গবেষণা ধ্রুবক সহগ বুসিনেস্ক সিস্টেমে কেন্দ্রীভূত
    • পরিবর্তনশীল সহগ ক্ষেত্রে তাত্ত্বিক বিশ্লেষণ তুলনামূলকভাবে অভাব
    • এই ধরনের বিচ্ছুরণ সিস্টেমের বিপরীত সমস্যা গবেষণা এখনও শূন্য
  4. গবেষণা প্রেরণা:
    • পরিবর্তনশীল সহগ বুসিনেস্ক সিস্টেমের তাত্ত্বিক বিশ্লেষণের ফাঁক পূরণ করা
    • বাস্তব সামুদ্রিক প্রকৌশল প্রয়োগের জন্য আরও নির্ভুল গাণিতিক মডেল এবং সংখ্যাসূচক সরঞ্জাম প্রদান করা
    • বিচ্ছুরণ তরঙ্গ গতিশীল সিস্টেমের বিপরীত সমস্যা গবেষণার নতুন দিকনির্দেশনা খোলা

মূল অবদান

  1. পরিবর্তনশীল সহগ এক-মাত্রিক বুসিনেস্ক সিস্টেমের সুস্থতার তত্ত্ব প্রতিষ্ঠা করা: সমাধানের অস্তিত্ব, অনন্যতা এবং প্রাথমিক ডেটার ক্রমাগত নির্ভরতা প্রমাণ করা
  2. গ্রীন ফাংশন কৌশলের উপর ভিত্তি করে তাত্ত্বিক বিশ্লেষণ পদ্ধতি বিকাশ করা: সমস্যাটিকে অরৈখিক অপারেটরের নির্দিষ্ট বিন্দু সমস্যায় রূপান্তরিত করা
  3. শক্তি সংরক্ষণ আইন প্রদান করা: সিস্টেমের শক্তি ফাংশনাল সময়ে সংরক্ষিত থাকে তা প্রমাণ করা
  4. দক্ষ সীমিত উপাদান সংখ্যাসূচক সমাধানকারী বাস্তবায়ন করা: FEniCS লাইব্রেরির সাথে স্থিতিশীল সংখ্যাসূচক স্কিম বিকাশ করা
  5. এই ধরনের বিচ্ছুরণ সিস্টেমের বিপরীত সমস্যা প্রথমবার অধ্যয়ন করা: পরিবর্তনশীল অপ্টিমাইজেশনের উপর ভিত্তি করে প্রাথমিক শর্ত পুনর্নির্মাণ পদ্ধতি প্রস্তাব করা

পদ্ধতির বিস্তারিত ব্যাখ্যা

কাজের সংজ্ঞা

নিম্নলিখিত পরিবর্তনশীল সহগ বুসিনেস্ক সিস্টেম অধ্যয়ন করা: {Nt+[(1+αc2(ξ)N)V]ξβ6Nξξt=0Vt+(c(ξ)N)ξ+12α(c2V2)ξβ6Vξξt=0\begin{cases} 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}

যেখানে (ξ,t)[0,L]×[0,T](\xi,t) \in [0,L] \times [0,T], প্রাথমিক শর্ত সহ: N(ξ,0)=N0(ξ),V(ξ,0)=V0(ξ)N(\xi,0) = N_0(\xi), \quad V(\xi,0) = V_0(\xi)

এবং ডিরিশলেট সীমানা শর্ত: N(0,t)=N(L,t)=0,V(0,t)=V(L,t)=0N(0,t) = N(L,t) = 0, \quad V(0,t) = V(L,t) = 0

তাত্ত্বিক বিশ্লেষণ পদ্ধতি

১. গ্রীন ফাংশন কৌশল সিস্টেমটি পুনরায় লিখুন: (Iβ6ξ2)Nt=ξ[(1+αc2(ξ)N)V]\left(I - \frac{\beta}{6}\partial_\xi^2\right)N_t = -\partial_\xi[(1 + \alpha c^2(\xi)N)V](Iβ6ξ2)Vt=ξ[c(ξ)N+12αc2V2]\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β6ξ2P = I - \frac{\beta}{6}\partial_\xi^2 এর গ্রীন ফাংশন ব্যবহার করুন: G(ξ,s)=12β/6cosh(Lξsβ/6)cosh(L(ξ+s)β/6)sinh(Lβ/6)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)}

२. নির্দিষ্ট বিন্দু উপপাদ্য সমস্যাটিকে একীভূত সমীকরণ সিস্টেমে রূপান্তরিত করুন, অরৈখিক অপারেটর AA সংজ্ঞায়িত করুন, বানাচ নির্দিষ্ট বিন্দু উপপাদ্যের মাধ্যমে সমাধানের অস্তিত্ব এবং অনন্যতা প্রমাণ করুন।

সংখ্যাসূচক পদ্ধতি

সীমিত উপাদান বিচ্ছেদ গ্যালারকিন পদ্ধতি ব্যবহার করুন, দুর্বল ফর্ম: t(N,v1L2+β6Nξ,v1ξL2)=[(1+αc2N)V],v1ξL2\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-পদ্ধতি (θ=1/2\theta = 1/2) ব্যবহার করুন সময় একীকরণের জন্য।

বিপরীত সমস্যা পদ্ধতি

উদ্দেশ্য ফাংশনাল সংজ্ঞায়িত করুন: J(N,V):=120L(N(ξ,T;N0,V0)NT(ξ)2+V(ξ,T;N0,V0)VT(ξ)2)dξ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][-20, 40]
  • জাল বিন্দু সংখ্যা: ৩০০০ স্থান গ্রিড পয়েন্ট
  • স্থান রেজোলিউশন: Δξ0.02\Delta\xi \approx 0.02
  • সময় ধাপ: Δt=8/30002.7×103\Delta t = 8/3000 \approx 2.7 \times 10^{-3}
  • মডেল প্যারামিটার: α=β=0.1\alpha = \beta = 0.1

প্রাথমিক শর্ত

স্থানীয়করণ গাউসিয়ান প্রোফাইল ব্যবহার করুন: N(ξ,0)=V(ξ,0)=e(ξ18)2N(\xi,0) = V(\xi,0) = e^{-(\xi-18)^2}

পরিবর্তনশীল সহগ ফাংশন

ক্ষেত্র ১: মসৃণ দোলনশীল-গাউসিয়ান বিঘ্ন c(ξ)=1+0.3sin(π5ξ)+0.6e(ξ8)2c(\xi) = 1 + 0.3\sin\left(\frac{\pi}{5}\xi\right) + 0.6e^{-(\xi-8)^2}

ক্ষেত্র २: অংশবিশেষ ধ্রুবক ফাংশন (স্তরযুক্ত মাধ্যম অনুকরণ) c(ξ)={0.8,ξ<201.5,20<ξ222.0,22<ξ27c(\xi) = \begin{cases} 0.8, & \xi < 20 \\ 1.5, & 20 < \xi \leq 22 \\ 2.0, & 22 < \xi \leq 27 \\ \vdots \end{cases}

পরীক্ষামূলক ফলাফল

তাত্ত্বিক ফলাফল

উপপাদ্য ३.१ (স্থানীয় সুস্থতা): T0>0T_0 > 0 বিদ্যমান, যাতে সিস্টেমটি LT02×LT02L^2_{T_0} \times L^2_{T_0} এ অনন্য সমাধান (N,V)(N,V) রয়েছে, এবং অনুমান সন্তুষ্ট করে: (N,V)LT02×LT02(N0,V0)H01×H01T01/2exp(L1/2(β/6)1/2DT0)\|(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)

উপপাদ্য ३.२ (শক্তি সংরক্ষণ): VL22+αcN1/2VL22+c1/2NL22=ধ্রুবক\|V\|_{L^2}^2 + \alpha\|c|N|^{1/2}V\|_{L^2}^2 + \||c|^{1/2}N\|_{L^2}^2 = \text{ধ্রুবক}

সংখ্যাসূচক পরীক্ষার ফলাফল

१. সামনের দিকের সমস্যা সমাধান:

  • পরিবর্তনশীল সহগ মাধ্যমে তরঙ্গ প্রচার সফলভাবে অনুকরণ করা হয়েছে
  • তলদেশের ভূগোল পরিবর্তনের কারণে তরঙ্গ গতি পরিবর্তন, বিস্তার সমন্বয় এবং বিচ্ছুরণ প্রভাব পর্যবেক্ষণ করা হয়েছে
  • সংখ্যাসূচক পদ্ধতি ভাল স্থিতিশীলতা এবং নির্ভুলতা প্রদর্শন করেছে

२. বিপরীত সমস্যা পুনর্নির্মাণ:

  • L-BFGS-B অ্যালগরিদম ৪ পুনরাবৃত্তির মধ্যে প্রাথমিক শর্তের পুনর্নির্মাণ উল্লেখযোগ্যভাবে উন্নত করেছে
  • পুনর্নির্মাণ করা তরঙ্গ উচ্চতা এবং তরল বেগ ধীরে ধীরে লক্ষ্য প্রোফাইলের কাছাকাছি পৌঁছেছে
  • পরিবর্তনশীল সমন্বয় কৌশলের কার্যকারিতা প্রমাণ করেছে

সংখ্যাসূচক যাচাইকরণ

  • শক্তি সংরক্ষণ আইন সংখ্যাসূচক গণনায় খুব ভালভাবে বজায় রাখা হয়েছে
  • বিভিন্ন জাল রেজোলিউশনের অধীনে ফলাফলের সামঞ্জস্য সংখ্যাসূচক পদ্ধতির সংযোগ যাচাই করেছে

সম্পর্কিত কাজ

१. বুসিনেস্ক সিস্টেম তত্ত্ব:

  • Bona & Chen (१९९८, २००२, २००४): ধ্রুবক সহগ ক্ষেত্রে ভিত্তি তত্ত্ব প্রতিষ্ঠা করেছেন
  • এই পেপার পরিবর্তনশীল সহগ ক্ষেত্রে সম্প্রসারিত করেছে, তাত্ত্বিক ফাঁক পূরণ করেছে

२. সংখ্যাসূচক পদ্ধতি:

  • বিচ্ছুরণ তরঙ্গ সমীকরণে সীমিত উপাদান পদ্ধতির প্রয়োগ
  • এই পেপার FEniCS লাইব্রেরি সংযুক্ত করে দক্ষ সমাধানকারী বাস্তবায়ন করেছে

३. বিপরীত সমস্যা:

  • তরঙ্গ সমীকরণ বিপরীত সমস্যার ঐতিহ্যবাহী পদ্ধতি
  • এই পেপার প্রথমবার বুসিনেস্ক সিস্টেম বিপরীত সমস্যায় পরিবর্তনশীল অপ্টিমাইজেশন প্রয়োগ করেছে

উপসংহার এবং আলোচনা

প্রধান উপসংহার

१. পরিবর্তনশীল সহগ এক-মাত্রিক বুসিনেস্ক সিস্টেমের গাণিতিক তাত্ত্বিক ভিত্তি সফলভাবে প্রতিষ্ঠা করেছে २. কার্যকর সংখ্যাসূচক সমাধান পদ্ধতি বিকাশ করেছে ३. এই ধরনের বিচ্ছুরণ সিস্টেমের বিপরীত সমস্যা গবেষণার নতুন দিকনির্দেশনা খুলেছে ४. বাস্তব সামুদ্রিক প্রকৌশল প্রয়োগের জন্য তাত্ত্বিক এবং গণনামূলক সরঞ্জাম প্রদান করেছে

সীমাবদ্ধতা

१. তাত্ত্বিক দিক:

  • শুধুমাত্র স্থানীয় সুস্থতা প্রমাণ করেছে, বৈশ্বিক অস্তিত্ব গবেষণা অপেক্ষা করছে
  • বিপরীত সমস্যার অনন্যতা এবং স্থিতিশীলতা তাত্ত্বিক বিশ্লেষণ অনুপস্থিত

२. সংখ্যাসূচক দিক:

  • বিপরীত সমস্যা পরীক্ষা শুধুমাত্র আদর্শ শব্দহীন ক্ষেত্রে বিবেচনা করেছে
  • পর্যবেক্ষণ শব্দের প্রতি সংবেদনশীলতা বিশ্লেষণ অনুপস্থিত

३. প্রয়োগ দিক:

  • এক-মাত্রিক ক্ষেত্রে সীমাবদ্ধ, বাস্তব প্রয়োগ দ্বি-মাত্রিক সম্প্রসারণ প্রয়োজন
  • পরিবর্তনশীল সহগ ফাংশনের নির্বাচন তুলনামূলকভাবে সহজ

ভবিষ্যত দিকনির্দেশনা

१. তাত্ত্বিক উন্নয়ন:

  • বৈশ্বিক সুস্থতা তত্ত্ব
  • বিপরীত সমস্যার অনন্যতা এবং স্থিতিশীলতা বিশ্লেষণ
  • নিয়মিতকরণ তত্ত্ব

२. সংখ্যাসূচক উন্নতি:

  • শব্দ ডেটার শক্তিশালী পুনর্নির্মাণ পদ্ধতি
  • উচ্চ-মাত্রিক ক্ষেত্রে সম্প্রসারণ
  • মেশিন লার্নিং প্রযুক্তির সংমিশ্রণ

३. বাস্তব প্রয়োগ:

  • বাস্তব সামুদ্রিক ডেটার যাচাইকরণ
  • সামুদ্রিক পর্যবেক্ষণ সিস্টেমের সাথে সংমিশ্রণ

গভীর মূল্যায়ন

শক্তি

१. তাত্ত্বিক অবদান উল্লেখযোগ্য: পরিবর্তনশীল সহগ বুসিনেস্ক সিস্টেমের সুস্থতা প্রথমবার সিস্টেমেটিকভাবে অধ্যয়ন করেছে २. পদ্ধতি উদ্ভাবনী: পরিবর্তনশীল সহগ সমস্যা পরিচালনার জন্য গ্রীন ফাংশন কৌশল চতুরভাবে প্রয়োগ করেছে ३. সংখ্যাসূচক বাস্তবায়ন সম্পূর্ণ: সম্পূর্ণ সংখ্যাসূচক সমাধান কাঠামো প্রদান করেছে ४. সমস্যা নতুন: বিচ্ছুরণ সিস্টেম বিপরীত সমস্যা গবেষণার নতুন ক্ষেত্র খুলেছে ५. লেখা স্পষ্ট: গাণিতিক অনুমান কঠোর, প্রকাশ স্পষ্ট

অপূর্ণতা

१. তাত্ত্বিক গভীরতা সীমিত: শুধুমাত্র স্থানীয় তত্ত্ব প্রতিষ্ঠা করেছে, বৈশ্বিক বিশ্লেষণ অনুপস্থিত २. পরীক্ষামূলক যাচাইকরণ অপর্যাপ্ত: বিপরীত সমস্যা পরীক্ষা তুলনামূলকভাবে সহজ, জটিল পরিস্থিতি যাচাইকরণ অনুপস্থিত ३. ব্যবহারিকতা যাচাইকরণ অপেক্ষা করছে: আরও বাস্তব ডেটা যাচাইকরণ প্রয়োজন ४. গণনামূলক জটিলতা বিশ্লেষণ অনুপস্থিত: অ্যালগরিদমের গণনামূলক জটিলতা বিশ্লেষণ করেনি

প্রভাব

१. একাডেমিক মূল্য: পরিবর্তনশীল সহগ বিচ্ছুরণ সিস্টেম তত্ত্ব উন্নয়নের ভিত্তি স্থাপন করেছে २. প্রয়োগ সম্ভাবনা: সামুদ্রিক প্রকৌশল, তরঙ্গ পূর্বাভাস এবং অন্যান্য ক্ষেত্রে সম্ভাব্য প্রয়োগ মূল্য রয়েছে ३. পদ্ধতিগত অবদান: গ্রীন ফাংশন পদ্ধতি এবং পরিবর্তনশীল অপ্টিমাইজেশনের সংমিশ্রণ সম্পর্কিত সমস্যার জন্য নতুন চিন্তাভাবনা প্রদান করেছে

প্রযোজ্য পরিস্থিতি

१. অগভীর জল তরঙ্গ প্রচার মডেলিং २. উপকূলীয় প্রকৌশলে তরঙ্গ বিশ্লেষণ ३. সামুদ্রিক পরিবেশ পর্যবেক্ষণে বিপরীত সমস্যা ४. সংখ্যাসূচক সামুদ্রিক বিজ্ঞানে ডেটা সমন্বয়

সংদর্ভ

পেপারটি ২০টি গুরুত্বপূর্ণ সংদর্ভ উদ্ধৃত করেছে, যা বুসিনেস্ক সিস্টেম তত্ত্ব, সংখ্যাসূচক পদ্ধতি, ফাংশনাল বিশ্লেষণ এবং অন্যান্য সম্পর্কিত ক্ষেত্রের ক্লাসিক কাজ অন্তর্ভুক্ত করে, গবেষণার জন্য দৃঢ় তাত্ত্বিক ভিত্তি প্রদান করে।