2025-11-22T13:22:16.398648

A time-space B-spline integrator for the Burgers' equation

Dag, Uğurluoğlu, Adar
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.
academic

বার্গার্স সমীকরণের জন্য একটি সময়-স্থান B-স্প্লাইন ইন্টিগ্রেটর

মৌলিক তথ্য

  • পেপার আইডি: 2510.09408
  • শিরোনাম: A time-space B-spline integrator for the Burgers' equation
  • লেখক: ইদ্রিস ড্যাগ, সেরকান উগুরলুওগলু, নিহাত আদার (এস্কিশেহির অসমানগাজি বিশ্ববিদ্যালয়, তুরস্ক)
  • শ্রেণীবিভাগ: math.NA cs.NA
  • প্রকাশনার সময়: ২০২৫ সালের ১০ অক্টোবর (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_t + uu_x - \nu u_{xx} = 0, \quad x \in [a,b], \quad t \in (0,T]

প্রাথমিক শর্ত সহ: u(x,0)=f(x)u(x,0) = f(x)

সীমানা শর্ত: u(a,t)=u0,ux(a,t)=u0′,u(b,t)=um,ux(b,t)=um′u(a,t) = u_0, u_x(a,t) = u'_0, u(b,t) = u_m, u_x(b,t) = u'_m

যেখানে ν\nu সান্দ্রতা ধ্রুবক।

মডেল আর্কিটেকচার

সময় বিচ্ছেদ (দ্বিঘাত B-স্প্লাইন)

१. গ্রিড বিভাজন: সময় ডোমেইন [c,d][c,d] কে গ্রিড পয়েন্টে সমানভাবে বিভক্ত করা c=t0<t1<⋯<tk=dc = t_0 < t_1 < \cdots < t_k = d, যেখানে tj=t0+jΔtt_j = t_0 + j\Delta t

२. দ্বিঘাত B-স্প্লাইন ভিত্তি ফাংশন: Bj(t)∈C1[a,b]B_j(t) \in C^1[a,b], j=−1,…,nj = -1,\ldots,n সংজ্ঞায়িত করা: Bj(t)=1Δt{(tj+2−t)2−3(tj+1−t)2+3(tj−t)2,tj−1≤t≤tj(tj+2−t)2−3(tj+1−t)2,tj≤t≤tj+1(tj+2−t)2,tj+1≤t≤tj+20,অন্যথায়B_j(t) = \frac{1}{\Delta t} \begin{cases} (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}

३. আনুমানিক সমাধান প্রতিনিধিত্ব: U(x,t)=∑j=−1nδ(x,tj)Bj(t)U(x,t) = \sum_{j=-1}^n \delta(x,t_j)B_j(t)

স্থান বিচ্ছেদ (ত্রিঘাত B-স্প্লাইন)

१. স্থান গ্রিড: স্থান ডোমেইন [a,b][a,b] কে সমানভাবে বিভক্ত করা, xi=x0+ihx_i = x_0 + ih, i=0,…,mi = 0,\ldots,m

२. ত্রিঘাত B-স্প্লাইন ভিত্তি ফাংশন: Bi(x)∈C2[a,b]B_i(x) \in C^2[a,b], i=−1,…,mi = -1,\ldots,m

३. প্যারামিটার সম্প্রসারণ: δ(x,tj)=∑i=−1m+1σijBi(x)\delta(x,t^j) = \sum_{i=-1}^{m+1} \sigma_i^j B_i(x)

প্রযুক্তিগত উদ্ভাবন পয়েন্ট

१. স্তরযুক্ত সমন্বয় কৌশল: প্রথমে সময় দিকে দ্বিঘাত B-স্প্লাইন দিয়ে সমন্বয় করা, তারপর স্থান দিকে ত্রিঘাত B-স্প্লাইন দিয়ে সমাধান করা २. সমন্বয় পয়েন্ট পদ্ধতি: গ্রিড পয়েন্টে অন্তরীয় সমীকরণ সন্তুষ্ট করতে বাধ্য করা ३. রৈখিকীকরণ প্রক্রিয়াকরণ: অরৈখিক পদগুলি রৈখিক করতে পূর্ববর্তী সময় পদক্ষেপের মান নির্বাচন করা ४. পুনরাবৃত্তিমূলক পরিমার্জন: প্রতিটি সময় পদক্ষেপের নির্ভুলতা উন্নত করতে পুনরাবৃত্তিমূলক প্রক্রিয়া গ্রহণ করা

পরীক্ষামূলক সেটআপ

পরীক্ষা সমস্যা

সমস্যা ১: শক ওয়েভ সমাধান

বিশ্লেষণাত্মক সমাধান: u(x,t)=x/t1+t/t0exp⁡(x2/(4νt)),t≥1,0≤x≤1u(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 যেখানে t0=exp⁡(1/(8ν))t_0 = \exp(1/(8\nu))

সমস্যা २: ভ্রমণকারী তরঙ্গ ফ্রন্ট সমাধান

বিশ্লেষণাত্মক সমাধান: u(x,t)=α+μ+(μ−α)exp⁡(η)1+exp⁡(η)u(x,t) = \frac{\alpha + \mu + (\mu - \alpha)\exp(\eta)}{1 + \exp(\eta)} যেখানে η=α(x−μt−γ)ν\eta = \frac{\alpha(x - \mu t - \gamma)}{\nu}, α=0.4\alpha = 0.4, μ=0.6\mu = 0.6, γ=0.125\gamma = 0.125

মূল্যায়ন সূচক

  • L∞L_\infty নর্ম: ∣U−UN∣∞=max⁡∣Uij−(UN)ij∣|U - U_N|_\infty = \max|U_i^j - (U_N)_i^j|
  • L2L_2 নর্ম: ∣U−UN∣2=∑i∣Uij−(UN)ij∣2|U - U_N|_2 = \sqrt{\sum_i |U_i^j - (U_N)_i^j|^2}
  • শিখর মান: শক ওয়েভের সর্বোচ্চ মান

বাস্তবায়ন বিবরণ

  • সময়-স্থান পদক্ষেপ: h=Δt=0.01,0.005,0.001h = \Delta t = 0.01, 0.005, 0.001
  • সান্দ্রতা প্যারামিটার: ν=0.01,0.005\nu = 0.01, 0.005
  • সময় পরিসীমা: t=3.25t = 3.25 (শক ওয়েভ সমস্যা), t=1.2t = 1.2 (ফ্রন্ট তরঙ্গ সমস্যা)

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

প্রধান ফলাফল

শক ওয়েভ সমস্যার ফলাফল

ν=0.01\nu=0.01h=Δth=\Delta tL∞L_\inftyL2L_2শিখর মান
0.010.04980.02110.2233
0.0050.05080.02070.2225
0.0010.05160.02060.2203
ν=0.005\nu=0.005L∞L_\inftyL2L_2শিখর মান
0.05820.01660.2503
0.05710.01590.2480
0.05600.01550.2467

ফ্রন্ট তরঙ্গ সমস্যার ফলাফল

ν=0.01\nu=0.01h=Δth=\Delta tL∞L_\inftyL2L_2
0.010.48750.1138
0.0050.29640.0638
0.0010.15140.0354
ν=0.005\nu=0.005L∞L_\inftyL2L_2
0.88410.2532
0.69060.1329
0.24270.0367

পরীক্ষামূলক আবিষ্কার

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

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

এই ক্ষেত্রের প্রধান গবেষণা দিকগুলির মধ্যে রয়েছে: १. ঐতিহ্যবাহী সংখ্যাসূচক পদ্ধতি: সীমিত পার্থক্য, সীমিত উপাদান, বর্ণালী পদ্ধতি ইত্যাদি २. B-স্প্লাইন পদ্ধতি: প্রধানত স্থান বিচ্ছেদের জন্য ব্যবহৃত, সময় সমন্বয়ের জন্য খুব কম ব্যবহৃত ३. সমন্বয় পদ্ধতি: গ্যালারকিন পদ্ধতি, সমন্বয় পদ্ধতি, জালবিহীন পদ্ধতি ইত্যাদি ४. বিশেষায়িত সমাধানকারী: বার্গার্স সমীকরণের জন্য বিশেষ সংখ্যাসূচক পদ্ধতি

এই পেপারের উদ্ভাবন সময় এবং স্থান সমন্বয়ে B-স্প্লাইন একযোগে প্রয়োগ করা, যা এই ক্ষেত্রের একটি ফাঁক পূরণ করে।

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

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

१. প্রস্তাবিত সময়-স্থান B-স্প্লাইন সমন্বয় পদ্ধতি বার্গার্স সমীকরণ কার্যকরভাবে সমাধান করতে পারে २. পদ্ধতি শক ওয়েভ এবং ফ্রন্ট প্রসারণ সমস্যায় ভাল নির্ভুলতা এবং স্থিতিশীলতা প্রদর্শন করে ३. সময়-স্থান পদক্ষেপ সামঞ্জস্য করে সমাধানের নির্ভুলতা এবং মসৃণতা নিয়ন্ত্রণ করা যায়

সীমাবদ্ধতা

१. গণনাগত জটিলতা: প্রতিটি সময় পদক্ষেপে রৈখিক সমীকরণ ব্যবস্থা সমাধান করতে হবে २. অরৈখিক প্রক্রিয়াকরণ: রৈখিকীকরণ এবং পুনরাবৃত্তিমূলক প্রক্রিয়া প্রয়োজন ३. প্যারামিটার সংবেদনশীলতা: সান্দ্রতা প্যারামিটার এবং গ্রিড প্যারামিটারের প্রতি তুলনামূলকভাবে সংবেদনশীল ४. প্রয়োগের পরিসীমা: প্রধানত একমাত্রিক বার্গার্স সমীকরণের জন্য

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

লেখকরা ভবিষ্যত কাজ নিম্নলিখিত বিষয়গুলিতে ফোকাস করবে বলে প্রস্তাব করেছেন:

  • উচ্চতর অর্ডার B-স্প্লাইন ব্যবহার করে নির্ভুলতা উন্নত করা
  • বহুমাত্রিক বার্গার্স সমীকরণে সম্প্রসারণ করা
  • অ্যালগরিদম দক্ষতা অপ্টিমাইজ করা

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

সুবিধা

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

অপূর্ণতা

१. তাত্ত্বিক বিশ্লেষণ অপর্যাপ্ত: সংযোগ এবং স্থিতিশীলতার কঠোর তাত্ত্বিক প্রমাণ অনুপস্থিত २. তুলনামূলক পরীক্ষা সীমিত: অন্যান্য উন্নত পদ্ধতির সাথে পর্যাপ্ত তুলনা নেই ३. গণনাগত দক্ষতা বিশ্লেষণ অনুপস্থিত: গণনা সময় এবং মেমরি ব্যবহার বিশ্লেষণ প্রদান করা হয়নি ४. প্যারামিটার নির্বাচন নির্দেশনা অপর্যাপ্ত: B-স্প্লাইন অর্ডার নির্বাচনের তাত্ত্বিক নির্দেশনা অনুপস্থিত

প্রভাব

१. একাডেমিক অবদান: PDE সংখ্যাসূচক সমাধানের জন্য নতুন চিন্তাভাবনা প্রদান করা २. ব্যবহারিক মূল্য: পদ্ধতি সহজ এবং বাস্তবায়ন সহজ, প্রকৌশল প্রয়োগের জন্য উপযুক্ত ३. সম্প্রসারণযোগ্যতা: ফ্রেমওয়ার্ক অন্যান্য অরৈখিক PDE-তে সাধারণীকরণ করা যায়

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

এই পদ্ধতি বিশেষভাবে উপযুক্ত:

  • বার্গার্স সমীকরণ উচ্চ নির্ভুলতার সাথে সমাধান করার প্রয়োজন এমন বৈজ্ঞানিক গণনা
  • শক ওয়েভ এবং ফ্রন্ট প্রসারণ ঘটনার সংখ্যাসূচক সিমুলেশন
  • গণনা নির্ভুলতার উচ্চ প্রয়োজন কিন্তু দক্ষতার মধ্যম প্রয়োজন এমন প্রয়োগ

সংদর্ভ

পেপারটি ১৬টি সম্পর্কিত সংদর্ভ উদ্ধৃত করেছে, যা বার্গার্স সমীকরণের সংখ্যাসূচক পদ্ধতি, B-স্প্লাইন তত্ত্ব, সমন্বয় পদ্ধতি এবং অন্যান্য সম্পর্কিত ক্ষেত্রের গুরুত্বপূর্ণ কাজ অন্তর্ভুক্ত করে, যা এই গবেষণার জন্য একটি দৃঢ় তাত্ত্বিক ভিত্তি প্রদান করে।