আমি যখন বিভিন্ন সীমানা শর্ত প্রয়োগ করি তখন আমি অ্যাডভেকশন-বিস্তারের সমীকরণের ভিন্ন আচরণ বুঝতে পারি না। আমার অনুপ্রেরণা হল প্রসারণ এবং advection অধীনে একটি বাস্তব শারীরিক পরিমাণ (কণা ঘনত্ব) এর অনুকরণ। কণার ঘনত্বটি প্রান্তগুলি থেকে প্রবাহিত না হওয়া অবধি অভ্যন্তরে সংরক্ষণ করা উচিত। এই যুক্তি অনুসারে, যদি আমি নিউম্যান সীমানা শর্তগুলি প্রয়োগ করি তবে সিস্টেমের শেষ যেমন (বাম এবং ডানদিকে) তারপর সিস্টেমটি"বন্ধ"হওয়া উচিতঅর্থাৎসীমান্তেরফ্লাক্সশূন্য হলে কোনও কণা পালাতে পারে না।
নীচের সমস্ত সিমিউলেশন জন্য, আমি advection-আশ্লেষ সমীকরণ খামখেয়াল-Nicolson discretization আবেদন করেছেন এবং সমস্ত সিমুলেশন আছে সীমানা শর্ত। যাইহোক, ম্যাট্রিক্স (সীমাবস্থা সারি) আমি অনুমতি দেয় প্রথম ও শেষ সারিরঅভ্যন্তর মান স্বাধীনভাবে পরিবর্তন করা। এটি শেষের পয়েন্টগুলি সম্পূর্ণ অন্তর্নিহিত হতে দেয় allows
আমি নীচে 4 টি আলাদা কনফিগারেশন নিয়ে আলোচনা করব, তার মধ্যে কেবলমাত্র আমি প্রত্যাশা করেছি। শেষে আমি আমার বাস্তবায়ন নিয়ে আলোচনা করব।
বিচ্ছিন্নতা কেবল সীমাবদ্ধ
এখানে অভিগমন শর্তাবলী গতিবেগ শূন্যে সেট করে বন্ধ করা আছে।
কেবলমাত্র points বোল্ডসাইম্বল = 0.5 (ক্র্যাঙ্ক-নিসকোলসন) এর সাথে বিস্তৃতি
পরিমাণটি সংরক্ষণ করা হয় না যেমন ডাল অঞ্চল হ্রাস করে দেখা যায়।
শুধুমাত্র আশ্লেষ সঙ্গে 0.5 (খামখেয়াল-Niscolson) এ অভ্যন্তর পয়েন্ট =, এবং = 1 (পূর্ণ অন্তর্নিহিত) গণ্ডি এβ
সীমানায় সম্পূর্ণ অন্তর্নিহিত সমীকরণ ব্যবহার করে আমি যা প্রত্যাশা করি তা অর্জন করি: কোনও কণা এড়ায় না । আপনি এটিকে অঞ্চলটি কণা বিচ্ছুরণ হিসাবে সংরক্ষণ করা দেখতে পাচ্ছেন। সীমানা পয়েন্টগুলিতে পছন্দটি কেন পরিস্থিতিটির পদার্থবিদ্যাকে প্রভাবিত করবে? এটি কি বাগ বা প্রত্যাশিত?
বিবর্তন এবং উত্সাহ
যখন অ্যাডভেকশন শব্দটি অন্তর্ভুক্ত করা হয় তখন সীমানায় মান সমাধানটিকে প্রভাবিত করে না বলে মনে হয়। যাইহোক, সমস্ত ক্ষেত্রে যখন সীমানাগুলি "উন্মুক্ত" বলে মনে হয় অর্থাৎ কণা সীমানা থেকে বাঁচতে পারে। কেন এই ক্ষেত্রে?
সমস্ত পয়েন্টে mb = 0.5 (ক্র্যাঙ্ক-নিকসোলসন) এর সাথে অ্যাডভিশন এবং ডিফিউশন
Advection এবং আশ্লেষ = 0.5 (খামখেয়াল-Niscolson) অভ্যন্তর বিন্দুতে, এবং = 1 (পূর্ণ অন্তর্নিহিত) গণ্ডি এβ
অ্যাডভেকশন-ডিফিউশন সমীকরণের বাস্তবায়ন
অ্যাডভেকশন-ডিফিউশন সমীকরণ দিয়ে শুরু করা,
ক্র্যাঙ্ক-নিকোলসন ব্যবহার করে লেখা দেয়,
লক্ষ্য করুন = 0.5 খামখেয়াল-Nicolson জন্য, = 1 সম্পূর্ণরূপে অন্তর্নিহিত জন্য, এবং, সম্পূর্ণরূপে স্পষ্ট জন্য = 0।β β
স্বরলিপিটি সহজ করার জন্য আসুন প্রতিস্থাপন করা যাক,
এবং সময়ের প্রাপ্ত জ্ঞাত মান move সরান ,
ফ্যাক্টরিং পদ দেয়,
যা আমরা ম্যাট্রিক্স আকারে as হিসাবে লিখতে পারি ,
নিউম্যান সীমানা শর্ত প্রয়োগ করা
এনবি আবারও ডাইরিভিশনের মাধ্যমে কাজ করছে বলে আমি মনে করি যে আমি ত্রুটিটি চিহ্নিত করেছি। সীমানা শর্তের সীমাবদ্ধ পার্থক্যটি লেখার সময় আমি একটি সম্পূর্ণ অন্তর্নিহিত স্কিম ( = 1) ধরেছিলাম। আপনি যদি এখানে ক্র্যাঙ্ক-নিকসোলসন স্কিম ধরে নেন তবে জটিলতা খুব দুর্দান্ত হয়ে উঠবে এবং আমি ডোমেনের বাইরে থাকা নোডগুলি সরিয়ে ফেলার ফলস্বরূপ সমীকরণগুলি সমাধান করতে পারিনি। তবে এটি প্রদর্শিত হবে, দুটি অজানা সঙ্গে দুটি সমীকরণ আছে, কিন্তু আমি এটি পরিচালনা করতে পারি না। এটি সম্ভবত উপরের প্রথম এবং দ্বিতীয় প্লটের মধ্যে পার্থক্য ব্যাখ্যা করে। আমি মনে করি আমরা উপসংহারে পৌঁছাতে পারি যে সীমানা পয়েন্টগুলিতে কেবলমাত্র = 0.5 সহ প্লটগুলি বৈধ।
বাম-প্রান্তে ফ্লাক্স ধরে নেওয়া পরিচিত (পুরোপুরি অন্তর্নিহিত ফর্ম অনুমান করে),
এটিকে কেন্দ্রিক-পার্থক্য হিসাবে লিখতে দেয়,
অতএব,
নোট করুন যে এটি একটি নোড প্রবর্তন করে যা সমস্যার ডোমেনের বাইরে। এই নোড দ্বিতীয় সমীকরণ ব্যবহার করে নির্মূল করা যেতে পারে। আমরা নোড হিসাবে লিখতে পারি ,
সীমানা অবস্থা থেকে পাওয়া of এর মান প্রতিস্থাপন = 1 সারির জন্য নিম্নলিখিত ফলাফল দেয় ,
চূড়ান্ত সারিতে ( = ) ফলনের জন্য একই পদ্ধতি সম্পাদন করা ,
শেষ অবধি সীমানা সারি অন্তর্ভুক্ত করে (সেটিং = 1) দেয়,
সুতরাং নিউমান সীমানা শর্তের সাথে আমরা ম্যাট্রিক্স সমীকরণ লিখতে পারি, ,
কোথায়,
আমার বর্তমান বোঝাপড়া
আমি মনে করি প্রথম এবং দ্বিতীয় প্লটের মধ্যে পার্থক্যটি উপরে বর্ণিত ত্রুটিটি চিহ্নিত করে ব্যাখ্যা করা হয়েছে।
শারীরিক পরিমাণ সংরক্ষণ সম্পর্কিত। আমি বিশ্বাস করি যে কারণটি এখানে উল্লেখ করা হয়েছে যে, আমি যে ফর্মটি লিখেছি তাতে অ্যাডভেশন সমীকরণ বিপরীত দিকে প্রসারণ করতে দেয় না তাই তরঙ্গ কেবল শূন্য-প্রবাহের সীমানা শর্তের সাথেও পার হয়ে যায় । সংরক্ষণ সম্পর্কিত আমার প্রাথমিক স্বীকৃতি কেবল তখনই প্রয়োগ হয় যখন অ্যাডভেকশন শব্দটি শূন্য হয় (এটি ক্ষেত্রের 2 নম্বর প্লটে সমাধান হয়)।এমনকি নিউমান শূন্য-প্রবাহের সীমানা শর্তগুলির সাথে এখনও সিস্টেমটি ছেড়ে যেতে পারে, কারণ এটি এই ক্ষেত্রে সঠিক সীমানা শর্ত রবিন সীমানা শর্ত যা মোট প্রবাহ নির্দিষ্ট করা । অধিকন্তু নিউনম্যান শর্তটি নির্দিষ্ট করে যে ভর বিস্তারের মাধ্যমে ডোমেনটি ছাড়তে পারে না , এটি অ্যাডভিশন সম্পর্কে কিছুই বলে না। সংক্ষেপে আমরা যা শুনেছি তা হ'ল প্রসারণের বদ্ধ সীমানা শর্ত এবং মেনে চলা সীমানা শর্ত। আরও তথ্যের জন্য উত্তরটি এখানে দেখুন, অ্যাডভেকশন-ডিফিউশন সমীকরণে গ্রেডিয়েন্ট জিরো বাউন্ডারি কনডিটন বাস্তবায়ন।
তুমি কি রাজি?