داوود دومیری گنجیاستاد مهندسی مکانیکEnglishارتباط با استاد
روش‌ها/HPM
از اینجا آغاز کنید · یادداشت آموزشی 10 دقیقه‌ای

روش اغتشاش هموتوپی

حل یک مسئلهٔ دشوار را از دنباله‌ای از معادلات ساده‌تر بسازید.

u(0.5) = 2/3مدل این یادداشتu(t) = 1 / (1 + t)10.51tuحل دقیق مرجع · u′ + u² = 0, u(0) = 1
حل دقیق مرجع که در سراسر یادداشت‌های آموزشی به کار می‌رود.

ایدهٔ روش، به زبان ساده.

پارامتر تعبیهٔ p را معرفی کنید تا معادله‌ای ساده‌تر را به‌طور پیوسته به معادلهٔ هدف پیوند دهد. حل مجهول را بر حسب توان‌های p بسط دهید، مسائل حاصل را به‌ترتیب حل کنید و تنها پس از بررسی سری حاصل، p=1 قرار دهید. پارامتر مصنوعی ابزاری برای ساختن حل است، نه اثبات همگرایی.

u′+pu2=0u(t;p)=∑n=0∞pnun(t)\begin{gathered}u^{\prime}+pu^2=0\\u(t;p)=\sum_{n=0}^{\infty}p^nu_n(t)\end{gathered}

خاستگاه و ارتباط با پژوهش‌های گنجی

این روش را Ji-Huan He توسعه داده است. گنجی و همکارانش HPM را در علوم حرارتی، مکانیک سیالات و دینامیک غیرخطی به کار گرفته و بررسی کرده‌اند.

ابزاری پرکاربرد در پژوهش‌های منتشرشدهٔ گنجی، از جمله جریان بورگرز، انتقال حرارت غیرخطی و مسائل نوسانی.

پیش‌نیازها: معادلات دیفرانسیل مرتبهٔ اول · انتگرال‌گیری · سری‌های توانی · تطبیق ضرایب.

روش را گام‌به‌گام دنبال کنید.

  1. ساختن هموتوپیبرای افت غیرخطی، u′+u²=0 را با u′+p u²=0 جایگزین کنید. در p=0، مسئلهٔ مقدار اولیه حل ثابت 1 دارد؛ در p=1، معادلهٔ هدف بازیابی می‌شود.
  2. بسط و تطبیقرابطهٔ u=Σpⁿuₙ را جای‌گذاری و ضرایب توان‌های یکسان p را برابر کنید. برای حفظ شرط اولیهٔ اصلی، u₀(0)=1 و برای n≥1، شرط uₙ(0)=0 را قرار دهید.
  3. حل سلسله‌مراتب معادلاتهر معادلهٔ ضریب را با استفاده از جمله‌های محاسبه‌شدهٔ قبلی انتگرال بگیرید: u₀=1، u₁=−t، u₂=t² و u₃=−t³.
  4. بررسی دامنهٔ اعتباردر p=1، این یک سری هندسی است. دامنهٔ همگرایی آن را بررسی کنید، با حل دقیق مقایسه کنید و باقی‌ماندهٔ هر تقریب متناهی را بسنجید.
مثال حل‌شده

همگرایی سری، مرزی دارد

u′ + u² = 0, u(0) = 1حل دقیق مرجع: u(t) = 1/(1+t)، t > −1
  1. معادلات ضرایب عبارت‌اند از u₀′=0 و برای n≥1، رابطهٔ uₙ′=−Σ[j=0…n−1]uⱼuₙ₋₁₋ⱼ.
  2. از آن‌ها uₙ(t)=(−t)ⁿ به دست می‌آید. بنابراین S_N(t)=1−t+t²−⋯+(−t)ᴺ.
  3. اتحاد مجموع متناهی هندسی نتیجه می‌دهد S_N−u_exact=−(−t)ᴺ⁺¹/(1+t).
  4. سری برای |t|<1 همگراست. در t=1، مجموع‌های جزئی بین 1 و 0 نوسان می‌کنند، درحالی‌که مقدار دقیق 1/2 است. برای t>1 جمله‌ها رشد می‌کنند.

بررسی نتیجه

در t=0.5، مقادیر S₄=0.6875 و u_exact=2/3 به دست می‌آیند و خطای مطلق برابر 1/48≈0.0208333 است.

ممکن است حل دقیق و همواری بیرون از بازهٔ همگرایی تقریب وجود داشته باشد. در خارج از این بازه، افزودن جمله‌های بیشتر می‌تواند پاسخ را بدتر کند.

کاوش در دفتر تعاملی

مرزهای کاربرد را بشناسید.

  • انتخاب حدس اولیه و عملگر خطی اهمیت دارد.
  • این ساختار به پارامتر فیزیکیِ کوچک نیاز ندارد؛ اما این به معنای همگرایی برای همهٔ مقادیر پارامترها یا زمان‌ها نیست.
  • مقایسه‌ها را با استفاده از حل‌های دقیق یا محاسبات عددیِ مستقلاً همگراشده اعتبارسنجی کنید.

زمینه‌های کاربرد

جریان از نوع بورگرزانتقال حرارتیارتعاشات غیرخطی

منابع اصلی را بخوانید.

این متن توضیحی آموزشی و اختصاصی است و ادعا نمی‌کند که همهٔ روش‌های ارائه‌شده را استاد گنجی ابداع کرده است. مدل حل‌شده عمداً در یادداشت‌ها مشترک انتخاب شده تا بتوان نتایج و محدودیت‌ها را مقایسه کرد.

بازگشت به همهٔ روش‌ها