حل مسائل الکترومغناطیس در کامسول با یک مثال عملی 

شبیه‌سازی پدیده‌های الکترومغناطیسی یکی از ارکان اساسی در طراحی، تحلیل و بهینه‌سازی ادوات نوین مهندسی برق، الکترونیک، مخابرات و فیزیک کاربردی است. نرم‌افزار کامسول مالی‌فیزیکس (COMSOL Multiphysics) به عنوان یکی از قدرتمندترین بسترهای شبیه‌سازی مبتنی بر روش اجزای محدود (FEM)، امکان مدلسازی دقیق برهم‌کنش میدان‌های الکتریکی و مغناطیسی با محیط‌های مختلف را فراهم می‌سازد. مسائلی نظیر تلفات ژول، تشعشع آنتن‌ها، ادوات مایکروویو، ماشین‌های الکتریکی و پدیده‌های فرکانس بالا با بهره‌گیری از ماژول‌های تخصصی این نرم‌افزار با دقت بسیار بالا قابل حل هستند.

آنچه کامسول را از سایر نرم‌افزارهای تحلیلی متمایز می‌کند، قابلیت اتصال هم‌زمان پدیده‌های فیزیکی (Multiphysics Coupling) مانند تاثیر حرارت بر رفتار الکترومغناطیسی مواد است. در این مقاله از سایت دیتک، با رویکردی گام‌به‌گام و کاربردی، به بررسی ساختار حل مسائل الکترومغناطیس پرداخته و با پیاده‌سازی یک پروژه عملی، شما را با فرآیند کامل مدل‌سازی، تعریف فیزیک، مش‌بندی و تحلیل نتایج آشنا خواهیم کرد.

کسب اطلاعات بیشتر: رفع خطاهای همگرایی (Convergence) در کامسول: راهنمای عیب‌یابی 

۱. مبانی تئوری الکترومغناطیس و معادلات حاکم در نرم‌افزار کامسول

پیش از شروع هرگونه شبیه‌سازی در محیط نرم‌افزاری، تسلط بر مبانی فیزیکی و معادلاتی که نرم‌افزار به صورت عددی حل می‌کند، کاملاً حیاتی است. تمامی شبیه‌سازی‌های الکترومغناطیسی بر پایه معادلات بنیادی ماکسول (Maxwell’s Equations) بنا شده‌اند. این معادلات رفتار میدان الکتریکی ($E$)، میدان جابه‌جایی الکتریکی ($D$)، شدت میدان مغناطیسی ($H$) و چگالی شار مغناطیسی ($B$) را در فضا و زمان به همراه چگالی بار و جریان توصیف می‌کنند.

کامسول با استفاده از روش المان محدود، معادلات دیفرانسیل جزئی حاکم بر این میدان‌ها را بر روی گره‌ها و المان‌های مش‌بندی‌شده تقریب می‌زند. شناخت فرکانس کاری و رژیم حاکم بر مساله (پایا، گذرا، شبه‌استاتیک یا فرکانس بالا) گام نخست در تعیین شکل ساده‌شده معادلات ماکسول است؛ چرا که انتخاب نادرست فرضیات تئوری می‌تواند منجر به نتایج غیرفیزیکی یا مصرف بیش‌ازحد منابع پردازشی سیستم شود.

تفاوت رژیم‌های فرکانسی و انتخاب فرمولاسیون مناسب

در مسائل الکترومغناطیس، مقایسه ابعاد هندسه مورد بررسی با طول موج کاری ($\lambda$) تعیین‌کننده رژیم حاکم است. زمانی که ابعاد فیزیکی سیستم بسیار کوچک‌تر از طول موج باشد، پدیده‌های تاخیری و انتشار موج قابل چشم‌پوشی بوده و سیستم در دسته مسائل شبه‌استاتیک یا فرکانس پایین قرار می‌گیرد. در این حالت، معادلات الکترواستاتیک یا مغناطواستاتیک برای تحلیل رفتارهایی نظیر خازن‌ها، ترانسفورماتورها و موتورهای الکتریکی به کار می‌روند.

در نقطه مقابل، هنگامی که ابعاد ساختار با طول موج میدان قابل مقایسه باشد، اثرات تشعشع، امپدانس موج و تلفات انتشاری ظاهر می‌شوند. در این دسته از مسائل، مانند موج‌برها، آنتن‌ها و خطوط انتقال فرکانس بالا، فرمولاسیون کامل موج تمام‌برداری (Full-Wave) باید مورد استفاده قرار گیرد تا تمامی مدهای انتشاری به درستی در فضای حل عددی محاسبه گردند.

بررسی نقش شرایط ساختاری مواد (Constitutive Relations)

رفتار محیط انتقال از طریق تانسورها یا اسکالرهای گذردهی الکتریکی ($\epsilon$)، تراوایی مغناطیسی ($\mu$) و هدایت الکتریکی ($\sigma$) در معادلات ظاهر می‌شود. کامسول به کاربر اجازه می‌دهد خواص مواد را به صورت خطی، غیرخطی (مانند منحنی اشباع B-H در هسته‌های فرومغناطیس)، ناهمسانگرد یا وابسته به دما و فرکانس تعریف کند.

در نظر گرفتن وابستگی متقابل این پارامترها به سایر متغیرهای فیزیکی، دقت حل عددی را به شکل چشمگیری افزایش می‌دهد. برای نمونه، تغییر رسانایی الکتریکی هادی‌ها ناشی از تلفات حرارتی ژول، می‌تواند توزیع چگالی جریان را تغییر دهد؛ بنابراین تعریف دقیق مدل ماده، گامی غیرقابل تفکیک از فرمولاسیون تئوری مساله به شمار می‌رود.

۲. انتخاب ماژول‌های تخصصی الکترومغناطیس در کامسول

نرم‌افزار کامسول برای پوشش طیف گسترده پدیده‌های الکترومغناطیسی، ابزارهای خود را در چندین ماژول تخصصی دسته‌بندی کرده است. هر ماژول شامل مجموعه‌ای از رابط‌های فیزیکی (Physics Interfaces) بهینه‌سازی‌شده برای حل سریع‌تر و دقیق‌تر معادلات خاص آن حوزه است. انتخاب ماژول اشتباه نه تنها روند مدل‌سازی را پیچیده می‌کند، بلکه ممکن است تنظیمات فیزیکی ضروری مانند شرایط مرزی تشعشعی یا پورت‌های ورودی موج را در دسترس کاربر قرار ندهد.

با درک صحیح از نیازهای پروژه، کاربر می‌تواند میان ماژول‌های فرکانس پایین و فرکانس بالا تفکیک قائل شود. ماژول AC/DC و ماژول RF پرکاربردترین بسترهای الکترومغناطیسی هستند که در کنار ماژول‌های تخصصی‌تری همچون Wave Optics و Plasma برای تحلیل سیستم‌های نوری و محیط‌های یونیزه به کار می‌روند. در ادامه، ویژگی‌ها و کاربردهای اصلی این دو ماژول اساسی را بررسی می‌کنیم.

ماژول AC/DC برای تحلیل‌های الکترواستاتیک و فرکانس پایین

ماژول AC/DC برای تحلیل میدان‌های استاتیک، شبه‌استاتیک و کاربردهای فرکانس پایین مهندسی برق توسعه یافته است. این ماژول رابط‌های کاربری متعددی از جمله Electrostatics (ES)، Electric Currents (EC) و Magnetic Fields (MF) را در اختیار مهندس قرار می‌دهد که برای طراحی ماشین‌های دوار، ترانسفورماتورها، عایق‌های فشارقوی، سنسورهای خازنی و القایی و سیستم‌های محافظت کاتدی کاربرد دارند.

یکی از قابلیت‌های کلیدی این ماژول، توانایی شبیه‌سازی سیم‌پیچ‌های چندرشته‌ای (Multi-Turn Coils) بدون نیاز به مدل‌سازی هندسی تک‌تک رساناها است. همچنین امکان اتصال مدارات خارجی الکتریکی (SPICE Circuit) به محیط المان محدود در این ماژول تعبیه شده است تا رفتار گذرای قطعات در یک مدار واقعی به دقت مورد تحلیل قرار گیرد.

ماژول RF برای ادوات مایکروویو، آنتن‌ها و انتشار امواج

برای مسائلی که با انتشار امواج الکترومغناطیسی، خطوط انتقال فرکانس رادیویی و ساختارهای تشدیدگر سر و کار دارند، ماژول RF انتخاب استاندارد است. این ماژول با بهره‌گیری از رابط فیزیکی Electromagnetic Waves, Frequency Domain (emw)، ابزارهای اختصاصی برای محاسبه ماتریس پراکندگی (S-Parameters)، الگوی تشعشعی دوردست (Far-Field Radiation Pattern) و امپدانس ورودی ادوات فراهم می‌آورد.

علاوه بر این، ماژول RF قابلیت مدل‌سازی ساختارهای باز و نامحدود را از طریق لایه‌های جاذب امواج (PML) میسر می‌سازد. تحلیل تلفات دی‌الکتریک، اثر پوسته‌ای در فرکانس‌های گیگاهرتزی و بهینه‌سازی فیلترهای موج‌بری از دیگر قابلیت‌های بنیادین این محیط نرم‌افزاری محسوب می‌شوند که در صنایع هوافضا و مخابرات نقش کلیدی ایفا می‌کنند.

۳. تعریف شرایط مرزی و لایه‌های انطباق‌پذیر در شبیه‌سازی

یکی از چالش‌های بنیادین در شبیه‌سازی‌های عددی، محدود کردن دامنه حل بی‌نهایت فیزیکی به یک دامنه محاسباتی محدود است. شرایط مرزی (Boundary Conditions) نحوه رفتار متغیرهای میدانی را بر روی مرزهای خارجی و سطوح داخلی مدل دیکته می‌کنند. اعمال نادرست یک شرط مرزی می‌تواند کل فیزیک حاکم بر میدان را منحرف ساخته و خطاهایی فاحش در نتایج محاسبه ظرفیت، اندوکتانس یا پترن تشعشعی ایجاد کند.

در کامسول، شرایط مرزی متنوعی تعبیه شده است که از ساده‌ترین حالت‌ها نظیر هادی کامل الکتریکی (PEC) تا شرایط پیشرفته بازتابش صفر و پورت‌های موج‌بری را پوشش می‌دهند. تسلط بر نحوه استفاده از این مرزها به کاربر اجازه می‌دهد اثرات محیط نامحدود اطراف دستگاه یا تقارن‌های موجود در هندسه را با کمترین هزینه محاسباتی پیاده‌سازی نماید.

هادی‌های ایده‌آل الکتریکی و مغناطیسی (PEC و PMC)

شرط مرزی هادی الکتریکی کامل (Perfect Electric Conductor – PEC) بیانگر مماس بودن میدان الکتریکی بر سطح و صفر بودن مولفه مماس آن ($n\times E=0$) است. این شرط مرزی به عنوان یک تقریب مهندسی بسیار عالی برای دیواره‌های فلزی با ضخامت بالا و هدایت الکتریکی بسیار زیاد به کار می‌رود و نیاز به مش‌بندی داخل حجم فلز را به طور کامل از بین می‌برد.

در نقطه مقابل، شرط مرزی هادی مغناطیسی کامل (Perfect Magnetic Conductor – PMC) مولفه مماس میدان مغناطیسی را بر روی سطح صفر در نظر می‌گیرد ($n\times H=0$). این شرط معمولاً در صفحات تقارن سیستم‌هایی به کار می‌رود که خطوط میدان مغناطیسی بر صفحه مرزی عمود هستند، که این امر به کاهش ۵۰ درصدی یا بیشتر حجم محاسباتی مدل کمک شایانی می‌کند.

پیاده‌سازی لایه‌های کاملاً منطبق (PML) و شرایط تشعشعی

هنگام مدل‌سازی مسائل باز مانند انتشار امواج از یک آنتن، امواج منتشرشونده در صورت برخورد به مرز خارجی دامنه ممکن است به داخل بازتاب شده و در محاسبات تداخل غیرواقعی ایجاد کنند. برای حل این مشکل، لایه‌های کاملاً منطبق (Perfect Matched Layers – PML) به صورت یک پوسته مجازی پیرامون دامنه تعریف می‌شوند که امواج را بدون هیچ‌گونه بازتابشی جذب و میرا می‌کنند.

شرایط مرزی تشعشعی (Scattering Boundary Condition) نیز گزینه‌ای سبک‌تر برای جذب امواج با زاویه تابش عمود بر مرز است. استفاده صحیح از PML نیازمند رعایت ضخامت حداقل نیم طول موج و استفاده از مش‌بندی ساختاریافته (Mapped/Swept) در آن لایه است تا عملکرد جذب موج دچار خطای بازتاب عددی نشود.

۴. مثال عملی: شبیه‌سازی گام‌به‌گام تشعشع و افت ولتاژ در موج‌بر مستطیلی

برای درک عمیق فرآیند مدل‌سازی، یک موج‌بر مستطیلی استاندارد باند X (با ابعاد مقطع ۲۲.۸۶ در ۱۰.۱۶ میلی‌متر) را به عنوان مثال عملی پیاده‌سازی می‌کنیم. هدف از این شبیه‌سازی، تحلیل انتشار مد پایه‌ای TE10TE10​، محاسبه پارامترهای پراکندگی (S11S11​ و S21S21​) در بازه فرکانسی ۸ تا ۱۲ گیگاهرتز و بررسی توزیع میدان الکتریکی در داخل ساختار است.

این سناریو نحوه کار با هندسه سه‌بعدی، تعریف خواص فیزیکی هوا و فلز، تنظیم پورت‌های موج‌بری و اعمال شرایط مرزی دیواره‌های هادی را به وضوح نمایش می‌دهد. با دنبال کردن این مراحل، ساختار استاندارد اجرای پروژه‌های مشابه در کامسول برای شما تثبیت خواهد شد.

ساخت هندسه، متریال‌گذاری و اعمال پورت‌های ورودی و خروجی

در مرحله نخست، در فضای سه‌بعدی نرم‌افزار، یک مکعب مستطیل (Block) با ابعاد مقطع باند X و طول ۵۰ میلی‌متر ایجاد می‌کنیم که نمایانگر حجم هوای داخل موج‌بر است. در بخش مواد (Materials)، متریال Air را از کتابخانه استاندارد به دامنه حجم نسبت می‌دهیم. از آنجا که تلفات دیواره‌های موج‌بر در طول‌های کوتاه ناچیز است، دیواره‌های جانبی را با شرط مرزی PEC مدل‌سازی می‌کنیم.

سپس در رابط فیزیکی Electromagnetic Waves, Frequency Domain، دو شرط مرزی Port بر روی سطوح انتهایی موج‌بر تعریف می‌کنیم. پورت شماره ۱ را به عنوان ورودی تحریک (Wave Excitation On) با مد مستطیلی TE10TE10​ و توان ۱ وات تنظیم کرده و پورت شماره ۲ را به عنوان پورت خروجی بدون تحریک (Wave Excitation Off) برای اندازه‌گیری امواج عبوری پیکربندی می‌نماییم.

اجرای گام‌های شبیه‌سازی در پنجره Model Builder

برای اجرای صحیح این پروژه در محیط نرم‌افزار کامسول، مراحل زیر را به ترتیب در درخت پروژه (Model Builder) پیاده‌سازی کنید:

1. انتخاب مدل: در پنجره Model Wizard، فضای 3D و رابط فیزیکی Electromagnetic Waves, Frequency Domain (emw)را برگزینید.
2. نوع مطالعه: نوع حل را بر روی Frequency Domain قرار داده و دکمه Done را بزنید.
3. تعریف پارامترها: در بخش Global Definitions -> Parameters، متغیرهای هندسی و فرکانسی را تعریف کنید:
   • a = 22.86[mm] (عرض موج‌بر)
   • b = 10.16[mm] (ارتفاع موج‌بر)
   • L = 50[mm] (طول موج‌بر)
4. ساخت هندسه: در شاخه Geometry 1 یک المان Block با ابعاد a، b و L بسازید و دکمه Build All Objects را کلیک کنید.
5. تخصیص متریال: از شاخه Add Material، گزینه Air را به کل دامنه حجم اختصاص دهید.
6. تنظیم پورت‌ها: در فیزیک emw، دو گره Port اضافه کنید:
   • Port 1: نوع مد Rectangular، سطح ورودی انتخاب شود، Wave Excitation فعال باشد.
   • Port 2: نوع مد Rectangular، سطح خروجی انتخاب شود، Wave Excitation غیرفعال باشد.
7. اجرای حل: در شاخه Study 1 -> Step 1: Frequency Domain، بازه فرکانسی را بر روی range(8[GHz], 0.2[GHz], 12[GHz]) تنظیم کرده و گزینه Compute را بزنید.

۵. استراتژی‌های مش‌بندی، تنظیمات حل‌کننده عددی و تحلیل نتایج

دقت نتایج شبیه‌سازی الکترومغناطیسی وابستگی مستقیمی به کیفیت مش‌بندی دارد. در مسائل موجی، اگر اندازه المان‌های مش بزرگ‌تر از کسری از طول موج باشد، پدیده پراکندگی عددی رخ داده و نتایج فاقد اعتبار خواهند بود. قاعده سرانگشتی در امواج الکترومغناطیسی این است که حداکثر اندازه المان (hmaxhmax​) باید حداقل کمتر از $\lambda /5$ یا ترجیحاً $\lambda /8$ در کوتاه‌ترین طول موج (بالاترین فرکانس) باشد.

پس از پایان مش‌بندی و اجرای حل‌کننده، مرحله پس‌پردازش (Post-Processing) آغاز می‌شود. کامسول نمودارهای متعددی را به صورت پیش‌فرض برای نمایش شدت میدان الکتریکی ($|E|$)، تلفات توان و بردارهای فلش چگالی توان ارائه می‌دهد. استخراج داده‌های دقیق ماتریس S و رسم نمودارهای دوبعدی بر حسب فرکانس، عملکرد موج‌بر را به طور کامل اعتبارسنجی می‌کند.

استقلال از مش و کنترل خطای عددی در فرکانس‌های بالا

برای تضمین درستی نتایج، اجرای مطالعه استقلال از مش (Mesh Independence Study) ضروری است. با تعریف اندازه المان‌ها بر اساس فرمول c_const / (freq_max * 8) در بخش مش، نرم‌افزار به صورت خودکار با افزایش فرکانس، اندازه المان‌ها را تطبیق می‌دهد. در نزدیکی لبه‌های تیز یا نواحی با گرادیان بالای میدان، استفاده از مش‌بندی ریزتر (Refinement) مانع از انحراف عددی می‌شود.

همچنین انتخاب حل‌کننده مستقیم (مانند MUMPS یا PARDISO) برای شبکه‌های با ابعاد متوسط، همگرایی پایدار را تضمین می‌کند. در مدل‌های بسیار بزرگ سه‌بعدی که با محدودیت حافظه RAM مواجه هستند، استفاده از حل‌کننده‌های تکراری (Iterative Solvers) نظیر GMRES به همراه پیش‌شرط‌سازهای مناسب مانند Multigrid الزامی خواهد بود تا از سرریز حافظه جلوگیری شود.

تحلیل خروجی‌های شبیه‌سازی و استخراج پارامترهای S

پس از اتمام محاسبات، برای ارزیابی عملکرد موج‌بر، یک گروه پلات یک‌بعدی (1D Plot Group) ایجاد کرده و نمودار پارامترهای emw.S11dBو emw.S21dB را بر حسب فرکانس رسم می‌کنیم. همان‌طور که انتظار می‌رود، برای فرکانس‌های بالاتر از فرکانس قطع مد پایه‌ای (حدود ۶.۵۶ گیگاهرتز)، پارامتر S21S21​ نزدیک به صفر دسی‌بل بوده که نشان‌دهنده عبور کامل موج بدون تضعیف از ساختار موج‌بر است.

همچنین با رسم پلات سه‌بعدی توزیع میدان الکتریکی (Multislice Plot)، توزیع سینوسی میدان الکتریکی در عرض موج‌بر به وضوح مشاهده می‌شود؛ جایی که میدان در مرکز به حداکثر مقدار خود رسیده و در دیواره‌های هادی به صفر میل می‌کند. این تطابق دقیق میان تئوری تحلیلی و شبیه‌سازی عددی، صحت مراحل پیاده‌سازی شده در کامسول را تایید می‌نماید.

جمع‌بندی و چک‌لیست فرآیند شبیه‌سازی الکترومغناطیس

تسلط بر شبیه‌سازی‌های الکترومغناطیسی در نرم‌افزار کامسول نیازمند درک یکپارچه از مفاهیم فیزیکی، تنظیمات شرایط مرزی و اصول کنترل خطای عددی است. با طی کردن اصولی مراحل معرفی‌شده، می‌توانید ساختارهای پیچیده‌تر نظیر آنتن‌های میکرواستریپ، سنسورهای مایکروویو و اتصالات مخابراتی را با اطمینان بالا تحلیل و بهینه‌سازی کنید.

در جدول زیر، چک‌لیست گام‌های اصلی برای شبیه‌سازی مسائل الکترومغناطیس در کامسول به صورت خلاصه آورده شده است: