পাঠ ৩৫ · ৫৭-এর মধ্যে · মডিউল ৭
Home / Courses / Numerical Methods / PDE শ্রেণীবিভাগ

PDE পরিচিতি — শ্রেণীবিভাগ

Intro to PDEs — classification
৯ মিনিট পড়া মধ্যম · Intermediate Python কোডসহ সম্পূর্ণ বাংলায়

এই পাঠে যা শিখবেন

  • PDE কী এবং ODE/BVP থেকে এটি কীভাবে কাঠামোগতভাবে ভিন্ন
  • ইলিপ্টিক, প্যারাবোলিক ও হাইপারবোলিক PDE-এর সংজ্ঞা এবং প্রতিটির ভৌত ব্যাখ্যা
  • তিনটি ক্যানোনিক্যাল সমীকরণ — ল্যাপ্লাস, হিট, ওয়েভ — এবং তারা বাস্তব-বিশ্বে কী প্রতিনিধিত্ব করে
  • ডিসক্রিমিন্যান্ট-ভিত্তিক শ্রেণীবিভাগ ফর্মুলা এবং কেন এই শ্রেণীবিভাগ সংখ্যাসূচক পদ্ধতি বাছাইয়ে গুরুত্বপূর্ণ

১ · PDE কী — ODE থেকে পার্থক্য

M6-M7-এর প্রতিটি সমীকরণে (ODE ও BVP) অজানা ফাংশন y নির্ভরশীল ছিল শুধু একটি স্বাধীন চলকের উপর — যেমন সময় t বা স্থান x। কিন্তু বাস্তব-বিশ্বের অনেক প্রপঞ্চ (তাপ ছড়ানো, শব্দ তরঙ্গ, তড়িৎচুম্বকীয় ক্ষেত্র) একই সাথে একাধিক চলকের উপর নির্ভরশীল — যেমন তাপমাত্রা u(x, t) নির্ভর করে স্থান x এবং সময় t দুটোর উপরেই।

যখন একটি ডিফারেনশিয়াল সমীকরণে একাধিক স্বাধীন চলকের সাপেক্ষে পার্শিয়াল ডেরিভেটিভ জড়িত থাকে (যেমন $\dfrac{\partial u}{\partial t}$ এবং $\dfrac{\partial^2 u}{\partial x^2}$ একই সমীকরণে), সেটি একটি পার্শিয়াল ডিফারেনশিয়াল ইকুয়েশন (PDE)। এই কোর্সের বাকি অংশে (L36-L37) আমরা দুটি ক্যানোনিক্যাল PDE-এর জন্য ফাইনাইট ডিফারেন্স স্কিম বাস্তবায়ন করব।

ODE
একটি স্বাধীন চলক (যেমন শুধু t) — M6-M7-এ কভার করা হয়েছে।
PDE
দুই বা তার বেশি স্বাধীন চলক (যেমন x ও t) — সমাধান একটি বহুমাত্রিক ফাংশন/গ্রিড।
কেন কঠিন
একটি ODE-এর সমাধান একটি বক্ররেখা; একটি ২D PDE-এর সমাধান একটি সম্পূর্ণ পৃষ্ঠতল (surface) — কম্পিউটেশনাল খরচ দ্রুত বাড়ে।

২ · তিনটি শ্রেণী — ইলিপ্টিক, প্যারাবোলিক, হাইপারবোলিক

দ্বিতীয়-অর্ডার লিনিয়ার PDE-কে সাধারণ আকারে লেখা যায় (দুই চলক x, y ধরে):

$$A\,u_{xx} + B\,u_{xy} + C\,u_{yy} + (\text{নিম্ন-অর্ডার পদ}) = 0$$

এই সমীকরণের আচরণ সম্পূর্ণভাবে নির্ধারিত হয় ডিসক্রিমিন্যান্ট $B^2 - 4AC$-এর চিহ্ন দিয়ে — ঠিক যেমন দ্বিঘাত সমীকরণের ডিসক্রিমিন্যান্ট নির্ধারণ করে মূল বাস্তব না জটিল।

ইলিপ্টিক ($B^2-4AC < 0$)
স্থিতিশীল অবস্থা (steady-state) বর্ণনা করে — কোনো সময়-নির্ভরতা নেই, সমাধান একটি সীমার সব বিন্দুতে একসাথে নির্ভরশীল। উদাহরণ: ল্যাপ্লাসের সমীকরণ $u_{xx}+u_{yy}=0$।
প্যারাবোলিক ($B^2-4AC = 0$)
ডিফিউশন/স্মুথিং প্রপঞ্চ বর্ণনা করে — সময়ের সাথে সমাধান মসৃণ হয়ে যায়, তথ্য অসীম গতিতে ছড়ায়। উদাহরণ: হিট ইকুয়েশন $u_t = \alpha u_{xx}$ (L36)।
হাইপারবোলিক ($B^2-4AC > 0$)
তরঙ্গ প্রচার বর্ণনা করে — তথ্য একটি সসীম গতিতে ভ্রমণ করে, শার্প ফিচার/অসংগতি সংরক্ষিত থাকতে পারে। উদাহরণ: ওয়েভ ইকুয়েশন $u_{tt} = c^2 u_{xx}$ (L37)।
কেন শ্রেণীবিভাগ গুরুত্বপূর্ণ

প্রতিটি শ্রেণীর জন্য উপযুক্ত সংখ্যাসূচক স্কিম আলাদা — একটি হাইপারবোলিক PDE-তে ভুল স্কিম ব্যবহার করলে (যেমন প্যারাবোলিক-উপযোগী স্কিম) সমাধান অস্বাভাবিকভাবে ছড়িয়ে যেতে পারে (numerical diffusion) বা অস্থিতিশীল হতে পারে। এই শ্রেণীবিভাগ জানা মানেই "কোন মেথড কোন সমস্যায়" প্রশ্নের একটি বড় অংশের উত্তর আগে থেকেই জানা — M11/L52-এর সিদ্ধান্ত-কাঠামোর সাথে সরাসরি সম্পর্কিত।

৩ · তিনটি ক্যানোনিক্যাল সমীকরণ ডিসক্রিমিন্যান্ট দিয়ে শ্রেণীবদ্ধ করা

নিচের কোড সেলে তিনটি ক্যানোনিক্যাল সমীকরণের (A, B, C) সহগ বসিয়ে ডিসক্রিমিন্যান্ট গণনা করে দেখানো হয়েছে (হিট ও ওয়েভ ইকুয়েশনে দ্বিতীয় চলকটি স্থান x নয়, সময় t, তবে একই ফর্মুলা প্রযোজ্য — y-এর জায়গায় t বসিয়ে):

Python
def classify(name, A, B, C):
    disc = B**2 - 4*A*C
    if disc < 0:
        kind = "ইলিপ্টিক (elliptic)"
    elif disc == 0:
        kind = "প্যারাবোলিক (parabolic)"
    else:
        kind = "হাইপারবোলিক (hyperbolic)"
    print(f"{name:28s} A={A:5.1f} B={B:5.1f} C={C:5.1f}  ->  B^2-4AC = {disc:7.2f}  =>  {kind}")

# ল্যাপ্লাসের সমীকরণ: u_xx + u_yy = 0
classify("ল্যাপ্লাস (Laplace)", A=1, B=0, C=1)

# হিট ইকুয়েশন: u_t = alpha*u_xx  ->  alpha*u_xx - u_t = 0  (u_tt সহগ C=0)
alpha = 1.0
classify("হিট ইকুয়েশন (heat)", A=alpha, B=0, C=0)

# ওয়েভ ইকুয়েশন: u_tt = c^2*u_xx  ->  c^2*u_xx - u_tt = 0
c = 1.0
classify("ওয়েভ ইকুয়েশন (wave)", A=c**2, B=0, C=-1)

    


    

      কোড সেলটি চালালে দেখা যায় — ল্যাপ্লাসের সমীকরণে B²-4AC = 0 - 4(1)(1) = -4 < 0 ⇒
      ইলিপ্টিক; হিট ইকুয়েশনে B²-4AC = 0 - 4(1)(0) = 0 ⇒ প্যারাবোলিক;
      আর ওয়েভ ইকুয়েশনে B²-4AC = 0 - 4(1)(-1) = 4 > 0 ⇒ হাইপারবোলিক। তিনটি
      সংখ্যাই সরাসরি সূত্র থেকে গণনা করা, এবং তিনটি ভিন্ন চিহ্ন তিনটি ভিন্ন ভৌত আচরণের সাথে হুবহু মেলে।
    


    

      সামনে কী আসছে
      

        L36-এ আমরা প্যারাবোলিক হিট ইকুয়েশনের জন্য একটি এক্সপ্লিসিট ফাইনাইট ডিফারেন্স স্কিম বাস্তবায়ন করব —
        একটি ছোট গ্রিডে তাপমাত্রা প্রোফাইল সত্যিকারভাবে সময়ের সাথে ডিফিউজ হতে দেখব। L37-এ হাইপারবোলিক ওয়েভ
        ইকুয়েশনের জন্য একটি স্কিম বাস্তবায়ন করব — একটি তরঙ্গ পালস সত্যিকারভাবে প্রচার ও প্রতিফলিত হতে দেখব।
        উভয় ক্ষেত্রেই একটি CFL-সদৃশ স্থিতিশীলতা শর্ত (M6/L32-এর স্টিফনেস আলোচনার সাথে সম্পর্কিত)
        মেনে চলতে হবে, নয়তো এক্সপ্লিসিট স্কিম সংখ্যাসূচকভাবে বিস্ফোরিত হবে।
      

    


  


  

    
    
  


  

    
 ভাবনার প্রশ্ন

    
প্রতিটি প্রশ্ন নিজে কিছুক্ষণ ভাবুন — তারপর "→ উত্তর" চাপুন।


    

      
        প্র ০১
        ল্যাপ্লাসের সমীকরণে কোনো সময়-চলক নেই — শুধু স্থানিক চলক (x, y)। তাহলে এটি
        "স্থিতিশীল অবস্থা" (steady-state) কেন বর্ণনা করে বলা হয়?
      
      

        

          ল্যাপ্লাসের সমীকরণ প্রায়ই একটি সময়-নির্ভর প্রপঞ্চের (যেমন হিট ইকুয়েশন) "চূড়ান্ত" বা "সাম্যাবস্থা"
          (equilibrium) অবস্থা বর্ণনা করে — যখন তাপমাত্রা প্রোফাইল আর সময়ের সাথে পরিবর্তিত হয় না
          ($\partial u/\partial t = 0$), তখন হিট ইকুয়েশন $u_t = \alpha u_{xx} + \alpha u_{yy}$ সরল হয়ে যায়
          $0 = \alpha(u_{xx}+u_{yy})$, অর্থাৎ ল্যাপ্লাসের সমীকরণে। তাই ল্যাপ্লাসের সমাধান বলে দেয় একটি সিস্টেম
          অসীম সময় পর কোন অবস্থায় "স্থির" হবে — এই কারণেই এটি স্বাভাবিকভাবে সময়-চলকবিহীন।
        

      

    


    

      
        প্র ০২
        বলা হয়েছে হাইপারবোলিক PDE-তে "তথ্য একটি সসীম গতিতে ভ্রমণ করে," কিন্তু প্যারাবোলিক PDE-তে
        "তথ্য অসীম গতিতে ছড়ায়।" হিট ইকুয়েশনের ক্ষেত্রে এর ভৌত অর্থ কী?
      
      

        

          তাত্ত্বিকভাবে, একটি রডের এক প্রান্তে তাৎক্ষণিকভাবে তাপ প্রয়োগ করলে, হিট ইকুয়েশনের সমাধান অনুযায়ী
          রডের অন্য প্রান্তের তাপমাত্রা সাথে সাথেই (যদিও অতি সামান্য পরিমাণে) পরিবর্তিত হয় — এটি
          পদার্থবিজ্ঞানগতভাবে কিছুটা অস্বাভাবিক (আলোর গতির চেয়েও দ্রুত প্রভাব), কিন্তু হিট ইকুয়েশনের একটি
          পরিচিত গাণিতিক আদর্শীকরণ (idealization)। বিপরীতে, ওয়েভ ইকুয়েশনে (L37) একটি পালস একটি নির্দিষ্ট,
          সসীম গতিতে (c) ভ্রমণ করে — এই গতির বাইরে কোনো বিন্দু তাৎক্ষণিকভাবে প্রভাবিত হয় না,
          যা শব্দ বা আলোক তরঙ্গের বাস্তব আচরণের সাথে অনেক বেশি মেলে।
        

      

    


    

      
        প্র ০৩
        ডিসক্রিমিন্যান্ট ফর্মুলা $B^2-4AC$ ব্যবহার করে হিট ইকুয়েশনকে "প্যারাবোলিক" বলা হচ্ছে, যদিও
        এটি সময়ের সাপেক্ষে মাত্র প্রথম-অর্ডার (u_t, u_tt নয়)। এই ভিন্নতা
        শ্রেণীবিভাগ প্রয়োগে কী প্রভাব ফেলে?
      
      

        

          ডিসক্রিমিন্যান্ট ফর্মুলাটি সাধারণ দ্বিতীয়-অর্ডার আকার $Au_{xx}+Bu_{xy}+Cu_{yy}+\ldots=0$-এর জন্য
          — হিট ইকুয়েশনে u_tt পদ নেই, তাই সেই পদের সহগ C=0 বসানো হয়েছে (নিম্ন-অর্ডার
          পদ u_t ফর্মুলার বাইরে থাকে, কারণ ডিসক্রিমিন্যান্ট শুধু সর্বোচ্চ-অর্ডার পদগুলোর উপর
          ভিত্তি করে গঠিত)। C=0 বসালে ডিসক্রিমিন্যান্ট ঠিক 0 হয়ে যায় — যা প্যারাবোলিক
          শ্রেণীর সংজ্ঞাগত বৈশিষ্ট্য। এটি দেখায় শ্রেণীবিভাগ ফর্মুলাটি নমনীয়ভাবে প্রযোজ্য, এমনকি যখন একটি
          সমীকরণের অর্ডার দুই চলকে সমান নয়।
        

      

    

  


  

    
 অনুশীলন

    
    
      
  1. 
        চিন্তা করুন: একটি সমীকরণ $2u_{xx} + 5u_{xy} + 2u_{yy} = 0$ দেওয়া থাকলে
        (A=2, B=5, C=2), এটি কোন শ্রেণীর PDE হবে বলে আপনার ধারণা — কোড চালানোর আগে গণনা করে
        দেখুন?
        
    
          
          
    
            
    
              $B^2-4AC = 5^2 - 4(2)(2) = 25 - 16 = 9 > 0$ — তাই এটি হাইপারবোলিক। কোড
              সেলে classify("উদাহরণ", A=2, B=5, C=2) যোগ করে চালালে ঠিক এই ফলাফলই আসবে, যা এই
              হাতে-করা গণনাকে নিশ্চিত করে।
            
    
          
    
        
    
      
    2. 
      
  2. 
        পরীক্ষা করুন: কোড সেলে classify ফাংশনে A=1, B=2, C=1
        বসিয়ে একটি নতুন কল যোগ করে Run চাপুন। ফলাফল কোন শ্রেণী দেখায়, এবং এই সীমান্ত-কেস (boundary case)
        কেন আকর্ষণীয়?
        
    
          
          
    
            
    
              A=1, B=2, C=1-এ B²-4AC = 4 - 4 = 0 — ফলাফল প্যারাবোলিক।
              এটি আকর্ষণীয় কারণ এটি দেখায় ডিসক্রিমিন্যান্ট ঠিক শূন্য হওয়া একটি "সীমান্ত" ঘটনা নয় — এটি
              একটি পুরোপুরি বৈধ, সুনির্দিষ্ট তৃতীয় শ্রেণী (প্যারাবোলিক), ইলিপ্টিক ও হাইপারবোলিকের মাঝামাঝি কোনো
              "অস্পষ্ট" অবস্থা নয়।
            
    
          
    
        
    
      
    3. 
    

  


  

    
 আরও পড়ুন · ABCL TECH-এ আপনার পরবর্তী পদক্ষেপ

    
  







  
  
     এই পাঠে ভুল পেয়েছেন? জানান
  





  
    
 আগের পাঠ

    
BVP-এর জন্য ফাইনাইট ডিফারেন্স মেথড