গ্যালারকিন পদ্ধতির সাথে সীমানা শর্তগুলি কীভাবে যুক্ত করা যায়?


21

আমি পিডিইগুলি সমাধানের জন্য গ্যালারকিন পদ্ধতি সম্পর্কে ওয়েবে কিছু সংস্থান পড়ছি, তবে আমি কিছু সম্পর্কে পরিষ্কার নই not নীচে আমি যা বুঝতে পেরেছি তার নিজস্ব অ্যাকাউন্ট রয়েছে।

নিম্নলিখিত সীমানা মান সমস্যা (BVP) বিবেচনা করুন:

L[u(x,y)]=0on(x,y)∈Ω,S[u]=0on(x,y)∈∂Ω

যেখানে একটি ২ য় অর্ডার লিনিয়ার ডিফারেনটিভেশন অপারেটর, \ ওমেগা set সাবসেট \ ম্যাথবিবি {আর} ^ 2 হল বিভিপির ডোমেন, \ আংশিক \ ওমেগা ডোমেনের সীমানা, এবং এস একটি প্রথম অর্ডার লিনিয়ার ডিফারেন্সিয়াল অপারেটর। আপনার (x, y) ফর্মের একটি সংক্ষিপ্ত বিবরণ হিসাবে প্রত্যাশা করুন :LΩ⊂R2∂ΩSu(x,y)

u(x,y)≈∑i=1Naigi(x,y)

যেখানে gআমি ফাংশনগুলির একটি সেট যা আমরা আপনার আনুমানিক u ব্যবহার করব তোমার দর্শন লগ করা। বিভিপিতে প্রতিস্থাপন:

Σআমিএকটিআমিএল[ছআমি(এক্স,Y)]=আর(একটি1,।।।,একটিএন,এক্স,Y)

যেহেতু আমাদের অনুমানের পরিমাণটি সঠিক নয়, তাই অবশিষ্টাংশগুলি আর শূন্য নয়। গ্যালারকিন-রিটজ-রালেহ পদ্ধতিতে আমরা \ ল্যাংল আর, জি_আই \ রেঙ্গেল = 0 প্রয়োজনের দ্বারা আনুমানিক ফাংশনগুলির সেটকে সম্মান করে আর কে কমিয়ে আনি । অত: পরআর⟨R,gi⟩=0

⟨R,gi⟩=∑j=1Naj⟨L[gj],gi⟩=0

সুতরাং, a_i সহগের সন্ধান করতে একটিআমিআমাদের অবশ্যই ম্যাট্রিক্স সমীকরণটি সমাধান করতে হবে:

(⟨এল[ছ1],ছ1⟩...⟨এল[ছএন],ছ1⟩.........⟨এল[ছ1],ছএন⟩...⟨এল[ছএন],ছএন⟩)(একটি1...একটিএন)=0

আমার প্রশ্ন: আমি কীভাবে সীমানা শর্তকে এর সাথে অন্তর্ভুক্ত করব?

সম্পাদনা: মূলত প্রশ্নটি বলেছিল যে S[u] ২ য় অর্ডার লিনিয়ার ডিফারেন্সিয়াল অপারেটর। আমি এটি 1 ম অর্ডার লিনিয়ার ডিফারেন্সিয়াল অপারেটরে পরিবর্তন করেছি।


1
বেকো, আপনাকে স্বাগতম scicomp! ক্রস পোস্টিং সম্পর্কিত আমাদের নীতি অন্যান্য স্ট্যাক এক্সচেঞ্জ সাইটের অনুসরণ করে । যদি আপনি একই শ্রেনীর (কমবেশি) বিভিন্ন শ্রোতাদের কাছে উপযুক্ত হন তবে ক্রস পোস্ট করা জায়েজ। কিছুক্ষণ পরেই আপনার প্রশ্নটি অন্য সাইটে স্থানান্তরিত হওয়ার জন্য জিজ্ঞাসা করা জায়েয, যদি আপনি মনে করেন যে প্রাথমিকভাবে যেখানে পোস্ট করা হয়েছে সেখানে আপনার প্রশ্নের সন্তোষজনকভাবে (বা মোটেও) উত্তর দেওয়া হচ্ছে না।
— জেফ অক্সবেরি

তবে ক্রস-পোস্টকে এটি সাধারণত আপত্তিজনক আচরণ বলে মনে করা হয়। আপনি যদি অঞ্চল ৫১ এর বিটা সাইটের তালিকার দিকে লক্ষ্য করেন তবে তাদের মধ্যে অনেকগুলি এক বছর পরেও সর্বজনীন বিটাতে রয়েছে। আমরা এখনও কিছু সময়ের জন্য ঘুরতে যাচ্ছি (কমপক্ষে একটি দীর্ঘ সময়ের স্কেল থেকে উত্তর দেওয়ার জন্য এই সাইটের বেশিরভাগ প্রশ্নের দরকার পড়ে)। তদ্ব্যতীত, mathআপনার প্রশ্নের উত্তর দেওয়ার ব্যবহারকারীরাও ব্যবহারকারী না হলে scicompতারা scicompযদি এ থেকে অনুলিপি করে পেস্ট করেন তবে তার বিপরীতে যথাযথ ক্রেডিট বা এট্রিবিউশন পাবেন না math।
— জেফ অক্সবেরি

1
আপনি কি নিশ্চিত যে দ্বিতীয় আদেশের অপারেটর? সাধারণভাবে, এটি একটি সু-পোজযুক্ত সমস্যা নয়। উদাহরণস্বরূপ, যদি , তবে আপনি এর জন্য জিজ্ঞাসা করছেন যা মারাত্মকভাবে অ-অনন্য সমাধান রয়েছে (যেমন কোনও সীমানা শর্তযুক্ত বৃহত্তর ডোমেনে PDE এর সমাধানও একটি সমাধান)। সাধারণত আমরা একটি (সম্ভবত ননলাইনার) প্রথম অর্ডার অপারেটর হতে বলি । এস = এল এল ইউSS=Lˉ Ω ΩLu=0Ω¯S
— জেদ ব্রাউন

2
এমনকি , আপনি এখনও অ-অনন্য সমাধানের দিকে তাকিয়ে রয়েছেন। যদি বিবেচনা করুন laplace অপারেটর এবং অন্য কোন দ্বিতীয় ক্রম রৈখিক অপারেটর। তারপরে কিছু ধ্রুবক ভেক্টর জন্য such মতো যে কোনও অন্য সমাধান তৈরি করার জন্য একটি সমাধানে যুক্ত করা যেতে পারে। এল এস ইউ ∇ ইউ = কে কেS≠LLSu∇u=kট
— ড্যান

1
@ জিফঅক্সবেরি আপনার মানসিক প্রশান্তির জন্য, এর সদৃশ প্রশ্নটি mathমুছে ফেলা হয়েছে। স্পষ্টতই আপনি প্রশ্নটি এখানে রাখার বিষয়ে ঠিক বলেছেন। আমি খুব সহায়ক প্রতিক্রিয়া পেয়েছি।
— বেকো

উত্তর:


15

গাণিতিক বিমূর্ততা ছাড়াই একটি দ্রুত এবং সাধারণ উত্তর। সীমানা শর্ত আরোপ করার জন্য বেশ কয়েকটি বিকল্প রয়েছে, যেমন

কঠোরভাবে গ্যালারকিন পদ্ধতিতে কথা বলতে আপনার বিসি কে সন্তুষ্ট করে এমন ফাংশনগুলির একটি সেট বেছে নেওয়া প্রয়োজন (উদাহরণস্বরূপ ভিত্তিক এবং / অথবা সমাধানের জন্য দায়ী বুদ্ধি এবং একটি আংশিক যোগফল) একজাতীয় শর্তগুলি সন্তুষ্ট করে এমন ভিত্তিতে ফাংশনগুলির উপর নির্ভর করে)uh=u0+uNu0uN

  • পেনাল্টি পদ্ধতি / ল্যাঞ্জরেঞ্জ বহুগুণ হয় যেখানে কেউ মূলত একটি পেনাল্টি শব্দ যুক্ত করে যা সীমানা শর্তকে অন্তর্ভুক্ত করে, যেমন A + যেখানে পৃথক সীমানা অবস্থার জন্য দায়ী ম্যাট্রিক্স এবং পদগুলির জন্য দায়ী। সীমাতে শর্তগুলি strongly়ভাবে চাপিয়ে দেওয়া হয় এবং অন্যথায় এটি দুর্বলভাবে চাপানো হয়। Au পছন্দ সিস্টেমের কন্ডিশনারকে প্রভাবিত করে। বি বি পি τ → ∞ τ ττ∗B=b+τ∗bpBbpτ→∞τ

  • তাউ পদ্ধতি যেখানে সীমানা শর্তের বিচ্ছিন্ন সংস্করণগুলি সহ স্পষ্টভাবে প্রয়োগ করা হয় সেখানে বেশ কয়েকটি সমীকরণ বিনিময় হয় (গ্যালারকিন সিস্টেমে সারিগুলির পরিবর্তন)। দ্রষ্টব্য: অতিরিক্ত বাউন্ডারি শর্তাদির সাহায্যে সিস্টেমগুলি অতিরিক্ত নির্ধারণ করাও একটি বিকল্প।

  • বিবেচনার আগে (রিটজ মেথড) গ্যালস ডাইভারজেন তত্ত্বের মাধ্যমে গ্যালার্কিন গঠনের পুনর্লিখনের জন্য ভলিউম ইন্টিগ্রালগুলিকে সীমানা ইন্টিগ্রালগুলিতে রূপান্তর করতে এবং তারপরে বিচ্যুতির আগে সরাসরি সূচনায় (সঠিক বা আনুমানিক) সীমানা শর্তগুলি অন্তর্ভুক্ত করে।

  • পরিশেষে, নোডাল / মডেল বিস্তারের মধ্যে সংযোগ কাজে লাগিয়ে নোডাল গ্যালারকিন পদ্ধতি নেওয়া সম্ভব যেখানে সিস্টেমের সমাধানটি একটি মডেল ভিত্তির পরিবর্তে ল্যাঞ্জরেন্জ ভিত্তিক সহগ হয়।


1
আমি হয় এটা, তাই না? τλ
— শুহালো


1
এটি কি "গ্যালারকিন পদ্ধতিতে আপনার বিসি সন্তুষ্ট করে এমন বেসিক ফাংশনগুলির একটি সেট বেছে নেওয়া দরকার ?"
— নট করুন

@ কেএনএল: আমিও তাই মনে করি, অন্য বাক্যটিও বোঝায় না। আমি একটি সম্পাদনা করব।
— ডেভিডিঘ

7

এক সম্ভাবনা সিস্টেম ম্যাট্রিক্স জড় করা হয় এবং ডান দিকে ভেক্টর খ , অজানা হিসেবে স্বাধীনতার নির্ধারিত মাত্রার, স্বাধীনতা অন্য কোন ডিগ্রী মত। এর পরে, একটি এবং খ সারি এবং নির্ধারিত dofs সঙ্গে যুক্ত কলাম zeroing, এবং সংশ্লিষ্ট তির্যক এন্ট্রিতে একটি এক নির্বাণ, এবং উপযুক্তভাবে পরিবর্তন RHS ভেক্টর দ্বারা পরিবর্তিত হয় খ ।AbAbb

আপনি যখন সারি শূন্য করেন, তখন একটিকে তির্যকের মধ্যে রাখুন এবং আরএইচএস পরিবর্তন করুন যাতে আপনি নির্ধারিত মান প্রয়োগ করেন, সিস্টেমটি আর প্রতিসাম্যহীন নয়। এজন্যই আপনি কলামগুলি শূন্য করেন এবং নির্ধারিত মানের জন্য অ্যাকাউন্টে rhs ভেক্টর পরিবর্তন করেন modb

অন্য সম্ভাবনাটি হ'ল নির্ধারিত ডফের তির্যকটিতে খুব বেশি সংখ্যক (সাধারণত 1e10) যুক্ত করা এবং তারপরে p * ˉ u এ আরএইচএস এন্ট্রি সেট করুন , যেখানে that ইউ সেই ডফের নির্ধারিত মান।pu¯u¯


6

সীমাবদ্ধ অবস্থার সাথে সীমাবদ্ধ অবস্থার সাথে কাজ করার সাধারণ সমস্যাটি বেশ জটিল হয়ে উঠতে পারে। কিন্তু যদি:

  • যেমন শুধুমাত্র আরোপ করে এস ( U ) = 0 ফর্মে তোলে তোমার দর্শন লগ করা যে কিছু সমান চ ( এক্স , Y ) উপর δ Ω ।S(u)S(u)=0uf(x,y)δΩ

  • তুমি তোমার উপাদানের ঠকান করতে পারেন যাতে করে বিভিন্ন উপাদানের গণ্ডি উপর সম্পূর্ণরূপে হয়δΩ

এটা আসলে খুব সহজ। আপনার সমীকরণ:

সাথে ( ⟨ ) প্রতিস্থাপন করা দরকার

(⟨L[g1],g1⟩…⟨L[gN],g1⟩………⟨L[g1],gN⟩…⟨L[gN],gN⟩)(a1…aN)=0
(⟨L[g1],g1⟩…⟨L[gN],g1⟩………⟨L[g1],gN⟩…⟨L[gN],gN⟩)(a1…aN)=b

যেখানে ডান পাশের ভেক্টর সীমানা পরিস্থিতি উপস্থাপন করে।b

নির্ধারণ , যে মান নির্ধারণ আপনার ভিত্তি উপাদান সেট U উপর δ Ω যাই হোক না কেন মান তারা সীমানা শর্ত সন্তুষ্ট থাকতে হবে। ইন ⟨ এল [ ছ ঞ ] , ছ আমি ⟩ , আপনি তাদের কাছ থেকে বাদ উচিত ছ ঞ কিন্তু ছ আমি (উপাদান একটি এই ফাংশন যে মিলা ইতিমধ্যে নির্ধারিত হয়েছে যাতে তারা ম্যাট্রিক্স মধ্যে অন্তর্ভুক্ত করা উচিত নয় সমীকরণ)। এর পরে, সেট আপ ⟨ আর , জি আমি ⟩ = এন ΣbuδΩ⟨L[gj],gi⟩gjgiaএকটি ম্যাট্রিক্স সমীকরণ, এবং উপাদানের মানখঅভ্যন্তরীণ পণ্য হিসেবে অধিকার পপ আউট উচিতএলউপাদান সঙ্গে আপনার আপনার অভ্যন্তর ভিত্তিতে অপারেটিং আপনার সীমানা ভিত্তিতে।

⟨R,gi⟩=∑j=1Naj⟨L[gj],gi⟩=0
bL

উত্তরের জন্য ধন্যবাদ ড্যান। "অনুধাবন করার জন্য ..." (যা অপরিহার্য অংশ, আমি মনে করি) শুরু হয় এমন অনুচ্ছেদটি আমি বুঝতে পারি না । আপনি কি এটি আরও কিছুটা স্পষ্ট করে তুলতে পারেন? খ
— বেকো

অন্যদিকে, আমি যে সমস্যাটি সমাধান করার চেষ্টা করছি তা কেবল আপনার দ্বিতীয় অবস্থার জন্যই সন্তুষ্ট: সীমানাটি একটি আয়তক্ষেত্র। প্রথম শর্ত হিসাবে, সীমানা শর্তগুলি সীমানায় ফাংশনের মান নির্দিষ্ট করে না। সীমানা শর্ত (ভালো কিছু ফাংশনের দ্বিতীয় ক্রম ডেরাইভেটিভস রৈখিক সমন্বয় মান উল্লেখ , যেখানে একটি , খ ধ্রুবক আছে। এছাড়াও সীমানা শর্তগুলি একজাতীয়। একটি∂2তোমার দর্শন লগ করা/∂এক্স∂Y+ +খ∂2তোমার দর্শন লগ করা/∂এক্স2=0)একটি,খ
— বেকো

@ বেকো: আপনার প্রশ্নে এবং এস কী রয়েছে সে সম্পর্কে আপনি আরও স্পষ্ট হতে চাইতে পারেন । বিভিন্ন ধরণের সীমানা পরিস্থিতি বিভিন্ন উপায়ে পরিচালনা করা যায়। এলএস
— ড্যান

1
আমি প্রশ্নটি আরও স্পষ্ট করার জন্য এটি সম্পাদনা করেছি, আমি মনে করি। আমি যে সঠিক সমস্যাটি সমাধান করতে চাইছি তা পোস্ট করতে চাই না কারণ আমি প্রশ্নটি যতটা পারি সাধারণভাবে রাখতে চাই। আমি মনে করি আমি সেই পদ্ধতিটি আরও ভালভাবে বুঝতে পারি।
— বেকো

@becko: আমরা এই স্থানান্তর করতে চান পারে চ্যাট যেমন দীর্ঘ ধরনের হচ্ছে।
— ড্যান

2

এখানে ভিত্তি পুনঃসংযোগ হিসাবে পরিচিত একটি পদ্ধতি যা বর্তমান থ্রেডে উল্লেখ করা হয়নি। আমি জেপি বয়েড বইটি থেকে উদ্ধৃত করছি, "চেবিশেভ এবং ফুরিয়ার স্পেকট্রাল পদ্ধতি", ২ য় সংস্করণ, অধ্যায় .5.৫ .:

যদি সমস্যাটি এর মধ্যে একজাতীয় সীমানা শর্ত থাকে, তবে এটি সর্বদা সমজাতীয় সীমানা শর্তগুলির সাথে সমতুল্য সমস্যা দ্বারা প্রতিস্থাপিত হতে পারে, যতক্ষণ না সীমানা শর্তটি লিনিয়ার হয়। প্রথম পদক্ষেপটি হ'ল সাধারণ ফাংশন বি ( এক্স ) চয়ন করা যা অসাধারণ সীমানা শর্তকে সন্তুষ্ট করে। এক তারপর নতুন পরিবর্তনশীল সংজ্ঞায়িত পারে বনাম ( এক্স ) এবং নতুন অত্যাচার ফাংশন ছ ( এক্স ) মাধ্যমে তোমার দর্শন লগ করা ( এক্স ) ≡ বনাম ( এক্স ) + + বি ( এক্স

এলতোমার দর্শন লগ করা=চ
বি(এক্স)বনাম(এক্স)ছ(এক্স) যাতে পরিবর্তিত সমস্যাটি এল ভি = জি যেখানে ভি ( এক্স ) একজাতীয় সীমানা শর্ত পূরণ করে। ...
তোমার দর্শন লগ করা(এক্স)≡বনাম(এক্স)+ +বি(এক্স)ছ(এক্স)≡চ(এক্স)-এলবি(এক্স)
এলবনাম=ছ
বনাম(এক্স)

শিফট ফাংশন সীমাবদ্ধতা ব্যতীত স্বেচ্ছাসেবী যে এটি অবশ্যই সমস্ত অসাধারণ সীমানা শর্ত পূরণ করতে পারে। তবে সর্বাধিক সহজ পছন্দটি সেরা: সর্বাধিক অর্ডারের বহুপদী প্রবৃত্তি যা কাজ করে।বি(এক্স)

এর পরে আমার নিজস্ব ব্যাখ্যা আসে:

  • "অসাধারণ সীমানা শর্ত" এর অর্থ এমন একটি শর্ত যা একটি ধ্রুবক থাকে, যেমন

    ∂এক্সতোমার দর্শন লগ করা(এক্স,Y)|এক্স=এক্স0=1।

    উপরের প্রোগ্রাম অনুসারে, সুবিধাজনক ফাংশন করে আপনি এটিটি ∂ x u ( x , y ) এ নেমেছেন | x = x 0 = 0।বি(এক্স)

    ∂এক্সতোমার দর্শন লগ করা(এক্স,Y)|এক্স=এক্স0=0।
  • আপনি একবার এইভাবে সমস্ত সীমানা শর্তগুলি একজাতীয় করে তুললে, আপনি আপনার ভিত্তির প্রসারণের দিকে ফিরে যেতে পারেন (যা আমি ধরে নিই যে পণ্য ভিত্তিতে করা হয়): সংশ্লিষ্ট বিসি অপারেটর প্রয়োগ করে, একটি ∂ x u ( x , y ) = ∑ i j a i j প্রাপ্ত করে

    তোমার দর্শন লগ করা(এক্স,Y) = Σআমিঞএকটিআমিঞφআমি(এক্স)φঞ(Y)
    এবংউপরের উদাহরণ অনুসারে এটি x = x 0 এর জন্য শূন্য হওয়া উচিত।
    ∂এক্সতোমার দর্শন লগ করা(এক্স,Y) = Σআমিঞএকটিআমিঞφআমি'(এক্স)φঞ(Y)
    এক্স=এক্স0
  • φআমি(এক্স)φআমি'(এক্স)|এক্স=এক্স0=0আমি

  • =1φআমি'(এক্স)|এক্স=এক্স0=1

    ∂এক্সতোমার দর্শন লগ করা(এক্স,Y)|এক্স=এক্স0 = Σআমিঞএকটিআমিঞφঞ(Y)
    1Yএকটিআমিঞ

এই সম্পূর্ণ পদ্ধতির সুন্দর জিনিসটি এটি তুলনামূলকভাবে বিমূর্ত স্তরে কাজ করছে। প্রয়োজনীয় উপাদানগুলি হ'ল বিসি অপারেটরের কেবল লিনিয়ারিটি এবং পণ্য ভিত্তিক ক্রিয়াগুলির ক্ষেত্রে একটি আনস্যাটজ। যেমন, এটি আনুমানিক পদ্ধতিতেও প্রযোজ্য।

আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.