سیستم های معادلات دیفرانسیل-جبری (Differential-Algebraic Equation Systems)، در ریاضیات (Mathematics)
انواع سیستم های دینامیکی (Dynamical Systems را در آموزش زیر شرح دادیم :
سیستم های معادلات دیفرانسیل-جبری (Differential-Algebraic Equation Systems) :
بسیاری از سیستم های فیزیکی و مهندسی، علاوه بر معادلات دیفرانسیل (که نرخ تغییرات را توصیف می کنند)، شامل قیود جبری (constraints) نیز هستند. این قیود ممکن است ناشی از پایستگی انرژی، اتصالات مکانیکی، یا قوانین کیرشهف در مدارها باشند. چنین سیستم هایی با معادلات دیفرانسیل-جبری (DAE) مدل می شوند. شکل کلی یک DAE به صورت زیر است:
\[ F(\dot{x}, x, t) = 0 \]که در آن
\[ \dot{x} \]مشتق زمانی متغیرهاست و تابع
\[ F \]ممکن است نسبت به
\[ \dot{x} \]تکین (singular) باشد. معمولا DAEها به فرم نیمه صریح (semi-explicit) نوشته می شوند:
\[ \begin{cases} \dot{x} = f(x, y, t) \\ 0 = g(x, y, t) \end{cases} \]در اینجا
\[ x \]متغیرهای دیفرانسیلی و
\[ y \]متغیرهای جبری (که مشتق شان در معادلات ظاهر نمی شود) هستند. معادله دوم یک قید جبری است که باید در تمام زمان ها برقرار باشد. این قیود، فضای状態 را به یک زیرخمینه (submanifold) از فضای اصلی محدود می کنند.
یک مثال کلاسیک، معادلات حرکت یک آونگ ساده با مختصات دکارتی است. اگر موقعیت آونگ را با
\[ (x,y) \]و طول ثابت
\[ L \]در نظر بگیریم، قید
\[ x^2 + y^2 = L^2 \]همراه با معادلات حرکت (که از قوانین نیوتن با نیروی کشش میله می آیند) یک DAE مرتبه ۳ (شاخص ۳) ایجاد می کند.
مفهوم شاخص (index) در DAEها اهمیت زیادی دارد. شاخص تعداد مشتقات گیری لازم برای تبدیل DAE به یک ODE است. هرچه شاخص بالاتر باشد، حل عددی دشوارتر است. DAEهای با شاخص ۱ نسبتا ساده هستند، اما شاخص های بالاتر (مانند شاخص ۳ در مکانیک) نیاز به روش های خاص کاهش شاخص (مانند روش های تصحیح) دارند.
کاربردهای DAE در شبیه سازی مدارهای الکتریکی (با استفاده از نرم افزارهایی مانند SPICE)، دینامیک چندجسمی (multibody dynamics)، سیستم های قدرت، و فرآیندهای شیمیایی (تعادل فازی و واکنش ها) فراوان است. در این مسائل، قیود فیزیکی مانند پایستگی بار، قانون ولتاژ کیرشهف، یا اتصال صلب بین اجسام، مستقیما معادلات جبری ایجاد می کنند.
تحلیل پایداری DAEها نیازمند توجه ویژه به قیود است. نقاط تعادل باید هم معادلات دیفرانسیل و هم معادلات جبری را ارضا کنند. پایداری خطی با بررسی یک ماتریس مشخصه تعمیم یافته (pencil) انجام می شود. همچنین مفهوم پایداری شاخص (index stability) مطرح است.
روش های عددی برای DAEها شامل روش های BDF (تفاضل عقب رو) و روش های ضمنی (implicit) رونگه-کوتا هستند. برای DAEهای با شاخص بالا، معمولا با مشتق گیری از قیود، شاخص را کاهش داده و سپس حل می کنند، اما باید مراقب بود که قیود اولیه نقض نشوند.
یکی از چالش های مهم در DAEها، تعیین شرایط اولیه سازگار (consistent initial conditions) است. شرایط اولیه باید همه قیود (و مشتقات آن ها) را ارضا کنند. نرم افزارهای شبیه سازی معمولا الگوریتم هایی برای یافتن چنین شرایطی دارند.
