معرفی
درست مانند لنزهای نوری که نور را متمرکز می کنند، لنزهای الکتریکی و مغناطیسی نیز می توانند پرتوهای ذرات باردار را متمرکز کنند. سیستمهای عدسیهای چهار قطبی مغناطیسی در فوکوس کردن پرتوهای یونی و ذرهای در شتابدهندههای مراکز هستهای و فیزیک ذرات مانند CERN، SLAC و ISIS کاربرد مشترکی پیدا میکنند. این مدل COMSOL Multiphysics مسیر یون های B 5 + را نشان می دهد که از سه عدسی چهار قطبی مغناطیسی متوالی عبور می کنند. مدل در مقطعی از هندسه تنظیم شده است.
تعریف مدل
چهارقطبی از مجموعه ای از چهار آهنربای دائمی تشکیل شده است، همانطور که در شکل 1 در زیر مشاهده می شود، جایی که آهنرباها با هم کار می کنند تا تقریب خوبی از یک میدان چهار قطبی ارائه دهند. برای تقویت میدان و حفظ آن در داخل سیستم، آهنرباها در یک استوانه آهنی قرار می گیرند.

شکل 1: نمای مقطعی یکی از چهارقطبی های مغناطیسی مورد استفاده در عدسی.
یون ها از طریق سیستمی متشکل از سه مجموعه چهار قطبی متوالی ارسال می شوند. طول وسط دو برابر بقیه است و 90 درجه حول محور مرکزی می چرخد. این بدان معنی است که قطبیت آهنرباهای آن معکوس شده است. شکل 2 نمای کاملی از لنز چهارقطبی مغناطیسی را نشان می دهد.

شکل 2: برش لنز چهار قطبی. قطب دوم (Q2) در مقایسه با Q1 و Q3 معکوس شده است. پس از عبور از لنز، یون ها به اندازه 1 متر رانش می شوند.
یک شتاب دهنده سیستم را با یون هایی تغذیه می کند که با سرعت 0.01c در امتداد محور مرکزی حرکت می کنند . برای مطالعه اثر تمرکز چهارقطبی، تعدادی از یونها را که از فاصله 3 سانتیمتری از محور مرکزی شروع میشوند و به طور مساوی در امتداد محیط یک دایره در صفحه عرضی توزیع میشوند، دنبال کنید. فرض بر این است که همه آنها دارای سرعت عرضی اولیه صفر هستند. هر چهار قطبی پرتو یونی را در امتداد یکی از محورهای عرضی متمرکز می کند و آن را در امتداد محور دیگر از بین می برد. اثر خالص پس از عبور از سیستم سه چهارقطبی و طول رانش در تمام جهات متمرکز است. همانطور که یون ها از سیستم خارج می شوند، همه آنها در شعاع 1 سانتی متر در صفحه عرضی قرار می گیرند.
این مدل در یک مقطع 2 بعدی از هر یک از دو چهار قطبی Q1 و Q3 یکسان تنظیم شده است. با غفلت از میدان های حاشیه ای، میدان مغناطیسی عرضی در یک نقطه معین در صفحه عرضی در Q2 به طور خودکار همان اندازه مربوط به Q1 و Q3 را دارد، اما در جهت مخالف است. بنابراین کافی است که فیلدها را در یکی از چهارقطبی ها مدل کنیم.
معادلات دامنه
میدان مغناطیسی با استفاده از معادله Magnetostatics، برای حل مولفه z پتانسیل مغناطیسی A (Wb/m) توصیف میشود:

در اینجا μ 0 = 4 π · 10 – 7 H/m نشان دهنده نفوذپذیری خلاء است، Br چگالی شار پسماند (T)، σ هدایت (S/m) و v سرعت محیط (m/s) است. ). در این مثال، رسانه در حال حرکت نیست. سمت راست معادله دارای جریان تحمیلی است که بر حسب چگالی جریان خارجی J z e (A/m 2 ) مشخص شده است. هیچ جریانی تحمیل نمی شود. دامنه آهن از فرمول کمی متفاوت از همان معادله استفاده می کند:

که در آن μ r = 4000 نفوذپذیری نسبی است. پتانسیل مغناطیسی در همه جا تعریف شده است به طوری که B = ∇ × A .
شرایط مرزی
میدان مغناطیسی تقریباً موازی با مرز بیرونی استوانه آهنی است. برای اجرای این امر، از شرط مرزی عایق مغناطیسی استفاده کنید ، با بیان اینکه A z = 0 .
نتایج و بحث
چگالی میدان مغناطیسی و خطوط جریان در مقطع Q1 یا Q3 در شکل 3 در زیر نشان داده شده است.

شکل 3: چگالی میدان مغناطیسی و خطوط جریان در مرکز یکی از آهنرباهای چهارقطبی.
هر یونی که از مجموعه عبور می کند، نیروی ماکسول برابر با F = q v × B را تجربه می کند ، که در آن v (m/s) سرعت یون است. سپس فرض کنید که مولفه z سرعت ثابت و بسیار بزرگتر از مولفه های x و y (عرضی) است. بنابراین فقط سهم نیروی مولفه z سرعت را در نظر بگیرید. برای یافتن موقعیت عرضی به عنوان تابعی از زمان، باید قانون دوم نیوتن را برای هر یون حل کنید، q v × B = m a ، که در آن m جرم یون (kg) و a است.نشان دهنده شتاب آن (m/s 2 ) است. اگر چگالی شار مغناطیسی محاسبه شده در Q1 برابر با B ‘ باشد ، و طول چهارقطبی i L i (m) باشد ، چگالی شاری که یون تجربه می کند به دست می آید.

که در آن t (s) زمان پرواز است. این وابستگی چگالی شار مغناطیسی به زمان پرواز به الگوریتم ردیابی ذرات به عنوان یک عبارت منطقی تغذیه می شود. شکل 4 زیر نحوه حرکت یون ها در صفحه عرضی را نشان می دهد.

شکل 4: وقتی یون ها وارد Q1 می شوند، به طور مساوی در اطراف دایره بزرگتر، 3 سانتی متر از محور z توزیع می شوند. Q1 در امتداد محور x تمرکز می کند و در امتداد محور y تمرکز نمی کند. نیروی وارد بر هر یون تقریباً متناسب با فاصله آن از محور z است، بنابراین با ورود یون ها به Q2، یون هایی که در محور x قرار دارند به سرعت به دور خود می چرخند و به سمت مرکز حرکت می کنند. Q3 حرکت را تثبیت می کند و همه یون ها را در مسیر درست قرار می دهد. در نهایت یونها به سمت کمری که کمی بیش از 1 متر فراتر از Q3 قرار دارد، رها میشوند.
مسیر کتابخانه برنامه: COMSOL_Multiphysics/Electromagnetics/چهار قطبی
دستورالعمل مدلسازی
از منوی File ، New را انتخاب کنید .
جدید
در پنجره جدید ، روی
Model Wizard کلیک کنید .

مدل جادوگر
1 | در پنجره Model Wizard روی ![]() |
2 | در درخت Select Physics ، AC/DC>Electromagnetic Fields>Magnetic Fields (mf) را انتخاب کنید . |
3 | روی افزودن کلیک کنید . |
4 | ![]() |
5 | در درخت انتخاب مطالعه ، General Studies>Stationary را انتخاب کنید . |
6 | ![]() |
تعاریف جهانی
پارامترهای 1
1 | در پنجره Model Builder ، در قسمت Global Definitions روی Parameters 1 کلیک کنید . |
2 | در پنجره تنظیمات برای پارامترها ، بخش پارامترها را پیدا کنید . |
3 | در جدول تنظیمات زیر را وارد کنید: |
نام | اصطلاح | ارزش | شرح |
م | 11 | 11 | عدد جرمی یون |
ز | 5 | 5 | شماره شارژ یونی |
L1 | 1[m] | 1 متر | طول چهارقطبی اول |
L2 | 2[m] | 2 متر | طول چهارقطبی دوم |
L3 | 1[m] | 1 متر | طول چهارقطبی سوم |
vz | 0.01*c_const | 2.9979E6 متر بر ثانیه | سرعت یون |
متر | M*mp_const | 1.8399E – 26 کیلوگرم | جرم یون |
q | Z*e_const | 8.0109E-19 C | شارژ یونی |
برادر | 8 [mT] | 0.008 T | چگالی شار چهار قطبی باقی می ماند |
در اینجا، c_const ، mp_const و e_const به ترتیب برای سرعت نور، جرم پروتون و بار ابتدایی ثابت های از پیش تعریف شده هستند.
هندسه 1
مستطیل 1 (r1)
1 | در نوار ابزار Geometry ، روی ![]() |
2 | در پنجره تنظیمات برای Rectangle ، بخش Size and Shape را پیدا کنید . |
3 | در قسمت متن Width ، 0.177 را تایپ کنید . |
4 | در قسمت متن ارتفاع ، 0.07 را تایپ کنید . |
5 | قسمت Position را پیدا کنید . در قسمت متن y ، -0.035 را تایپ کنید . |
6 | ![]() |
چرخش 1 (rot1)
1 | در نوار ابزار Geometry ، روی ![]() |
2 | فقط شی r1 را انتخاب کنید. |
3 | در پنجره تنظیمات برای چرخش ، قسمت چرخش را پیدا کنید . |
4 | در قسمت Angle text عدد 45 را تایپ کنید . |
5 | ![]() |
6 | ![]() |
دایره 1 (c1)
1 | در نوار ابزار Geometry ، روی ![]() |
2 | در پنجره تنظیمات دایره ، بخش اندازه و شکل را پیدا کنید . |
3 | در قسمت متن Radius ، 0.2 را تایپ کنید . |
4 | قسمت Position را پیدا کنید . در قسمت متن x ، 0.2 را تایپ کنید . |
5 | در قسمت متن y ، 0.2 را تایپ کنید . |
6 | ![]() |
7 | ![]() |
تقاطع 1 (int1)
1 | در نوار ابزار Geometry ، روی ![]() |
2 | در پنجره Graphics کلیک کنید و سپس Ctrl+A را فشار دهید تا هر دو شی انتخاب شوند. |
3 | در پنجره تنظیمات برای تقاطع ، روی ![]() |
چرخش 2 (rot2)
1 | در نوار ابزار Geometry ، روی ![]() |
2 | فقط شی int1 را انتخاب کنید. |
3 | در پنجره تنظیمات برای چرخش ، قسمت چرخش را پیدا کنید . |
4 | در قسمت نوشتار Angle ، 90، 180، 270 را تایپ کنید . |
5 | قسمت Input را پیدا کنید . چک باکس Keep input objects را انتخاب کنید . |
6 | ![]() |
7 | ![]() ![]() |
حالا دو دایره در مرکز مبدا درست کنید.
دایره 2 (c2)
1 | در نوار ابزار Geometry ، روی ![]() |
2 | در پنجره تنظیمات دایره ، بخش اندازه و شکل را پیدا کنید . |
3 | در قسمت متن Radius ، 0.2 را تایپ کنید . |
4 | ![]() |
5 | ![]() |
دایره 3 (c3)
1 | در نوار ابزار Geometry ، روی ![]() |
2 | در پنجره تنظیمات دایره ، بخش اندازه و شکل را پیدا کنید . |
3 | در قسمت متن Radius 0.12 را تایپ کنید . |
4 | ![]() |
با کم کردن آهنرباها از اتحاد دو دایره، یک جسم ترکیبی ایجاد کنید. این عملیات پیچیده را می توان با استفاده از عملیات Compose geometry و تعیین فرمول در یک مرحله انجام داد.
ترکیب 1 (co1)
1 | در نوار ابزار Geometry ، روی ![]() |
2 | در پنجره Graphics کلیک کنید و سپس Ctrl+A را فشار دهید تا همه اشیا انتخاب شوند. |
3 | در پنجره تنظیمات برای نوشتن ، بخش نوشتن را پیدا کنید . |
4 | در قسمت متن فرمول تنظیم ، c2+c3-(int1+rot2(1)+rot2(2)+rot2(3)) را تایپ کنید . |
5 | ![]() |
یک دایره دیگر در مرکز مبدا ایجاد کنید.
دایره 4 (c4)
1 | در نوار ابزار Geometry ، روی ![]() |
2 | در پنجره تنظیمات دایره ، بخش اندازه و شکل را پیدا کنید . |
3 | در قسمت متن Radius ، 0.2 را تایپ کنید . |
4 | ![]() |
5 | ![]() ![]() |
مواد را اضافه کنید
1 | در نوار ابزار Home ، روی ![]() |
2 | به پنجره Add Material بروید . |
3 | در درخت، Built-in>Iron را انتخاب کنید . |
4 | روی Add to Component در نوار ابزار پنجره کلیک کنید . |
5 | در نوار ابزار Home ، روی ![]() |
مواد
آهن (حصیر 1)
فقط دامنه 2 را انتخاب کنید.
آهنرباها
1 | در پنجره Model Builder ، روی Materials کلیک راست کرده و Blank Material را انتخاب کنید . |
2 | در پنجره تنظیمات مواد ، Magnets را در قسمت نوشتار Label تایپ کنید . |
3 | قسمت انتخاب موجودیت هندسی را پیدا کنید . روی Paste Selection کلیک کنید . ![]() |
4 | در کادر محاوره ای Paste Selection ، 1، 3، 5-6 را در فیلد متن Selection تایپ کنید. |
5 | روی OK کلیک کنید . |
6 | در پنجره تنظیمات برای Material ، قسمت Material Properties را پیدا کنید . |
7 | در درخت خواص مواد ، Electromagnetic Models>Remanent Flux Density را انتخاب کنید . |
8 | ![]() |
9 | قسمت محتوای مواد را پیدا کنید . در جدول تنظیمات زیر را وارد کنید: |
ویژگی | متغیر | ارزش | واحد | گروه اموال |
نفوذپذیری پس زدگی | murec_iso ; murecii = murec_iso، murecij = 0 | 1.05 | 1 | چگالی شار باقی می ماند |
آنها هنجار چگالی جریان باقی می مانند | شماره هنجار | برادر | تی | چگالی شار باقی می ماند |
میدان های مغناطیسی (MF)
از شرایط پیشفرض، عایق مغناطیسی ، در تمام مرزهای بیرونی استفاده کنید.
قانون آمپر 2
1 | در پنجره Model Builder ، در قسمت Component 1 (comp1) روی Fields Magnetic (mf) راست کلیک کرده و قانون آمپر را انتخاب کنید . |
2 | فقط دامنه 6 را انتخاب کنید. |
3 | در پنجره تنظیمات قانون آمپر ، بخش Constitutive Relation B-H را پیدا کنید . |
4 | از لیست مدل مغناطیسی ، چگالی شار پسماند را انتخاب کنید . |
5 | از لیست μ rec ، From material را انتخاب کنید . |
6 | بردار e را به صورت مشخص کنید |
-1/sqrt(2) | ایکس |
-1/sqrt(2) | y |
0 | z |
7 | بخش Constitutive Relation Jc-E را پیدا کنید . از لیست σ ، User defined را انتخاب کنید . قسمت Constitutive Relation D-E را پیدا کنید . از لیست ε r ، User defined را انتخاب کنید . |
قانون آمپر 3
1 | در نوار ابزار فیزیک ، روی ![]() |
2 | فقط دامنه 3 را انتخاب کنید. |
3 | در پنجره تنظیمات قانون آمپر ، بخش Constitutive Relation B-H را پیدا کنید . |
4 | از لیست مدل مغناطیسی ، چگالی شار پسماند را انتخاب کنید . |
5 | بردار e را به صورت مشخص کنید |
-1/sqrt(2) | ایکس |
1/sqrt (2) | y |
0 | z |
6 | بخش Constitutive Relation Jc-E را پیدا کنید . از لیست σ ، User defined را انتخاب کنید . قسمت Constitutive Relation D-E را پیدا کنید . از لیست ε r ، User defined را انتخاب کنید . |
قانون آمپر 4
1 | در نوار ابزار فیزیک ، روی ![]() |
2 | فقط دامنه 1 را انتخاب کنید. |
3 | در پنجره تنظیمات قانون آمپر ، بخش Constitutive Relation B-H را پیدا کنید . |
4 | از لیست مدل مغناطیسی ، چگالی شار پسماند را انتخاب کنید . |
5 | بردار e را به صورت مشخص کنید |
1/sqrt (2) | ایکس |
1/sqrt (2) | y |
0 | z |
6 | بخش Constitutive Relation Jc-E را پیدا کنید . از لیست σ ، User defined را انتخاب کنید . قسمت Constitutive Relation D-E را پیدا کنید . از لیست ε r ، User defined را انتخاب کنید . |
قانون آمپر 5
1 | در نوار ابزار فیزیک ، روی ![]() |
2 | فقط دامنه 5 را انتخاب کنید. |
3 | در پنجره تنظیمات قانون آمپر ، بخش Constitutive Relation B-H را پیدا کنید . |
4 | از لیست مدل مغناطیسی ، چگالی شار پسماند را انتخاب کنید . |
5 | بردار e را به صورت مشخص کنید |
1/sqrt (2) | ایکس |
-1/sqrt(2) | y |
0 | z |
6 | بخش Constitutive Relation Jc-E را پیدا کنید . از لیست σ ، User defined را انتخاب کنید . قسمت Constitutive Relation D-E را پیدا کنید . از لیست ε r ، User defined را انتخاب کنید . |
قانون آمپر 6
1 | در نوار ابزار فیزیک ، روی ![]() |
2 | فقط دامنه 4 را انتخاب کنید. |
3 | در پنجره تنظیمات قانون آمپر ، بخش Constitutive Relation B-H را پیدا کنید . |
4 | از لیست μ r ، User defined را انتخاب کنید . بخش Constitutive Relation Jc-E را پیدا کنید . از لیست σ ، User defined را انتخاب کنید . قسمت Constitutive Relation D-E را پیدا کنید . از لیست ε r ، User defined را انتخاب کنید . |
مش 1
1 | در پنجره Model Builder ، در قسمت Component 1 (comp1) روی Mesh 1 کلیک کنید . |
2 | در پنجره تنظیمات مش ، قسمت Physics-Controlled Mesh را پیدا کنید . |
3 | از فهرست اندازه عنصر ، Extra fine را انتخاب کنید . |
4 | ![]() ![]() |
مطالعه 1
در نوار ابزار صفحه اصلی ،
روی محاسبه کلیک کنید .

نتایج
هنجار چگالی شار مغناطیسی (mf)
نمودار پیش فرض هنجار چگالی شار مغناطیسی را نشان می دهد. برای مشاهده میدان مغناطیسی دستورالعمل های زیر را دنبال کنید.
سطح 1
1 | در پنجره Model Builder ، گره Magnetic Flux Density Norm (mf) را گسترش دهید ، سپس روی Surface 1 کلیک کنید . |
2 | در پنجره تنظیمات برای Surface ، روی Replace Expression در گوشه سمت راست بالای بخش Expression کلیک کنید . از منو، Component 1 (comp1)>Magnetic Fields>Magnetic>mf.normH – Norma میدان مغناطیسی – A/m را انتخاب کنید . |
3 | در نوار ابزار هنجار چگالی شار مغناطیسی (mf) ، روی ![]() |
4 | ![]() |
برای مشاهده نحوه حرکت یون ها در سیستم چهارقطبی، موارد زیر را انجام دهید:
مسیرهای یون ها
1 | در نوار ابزار صفحه اصلی ، روی ![]() |
2 | در پنجره تنظیمات گروه طرح دوبعدی ، Ions Trajectories را در قسمت نوشتار Label تایپ کنید . |
سطح 1
1 | روی Ions Trajectories کلیک راست کرده و Surface را انتخاب کنید . |
2 | در پنجره تنظیمات برای Surface ، روی Replace Expression در گوشه سمت راست بالای بخش Expression کلیک کنید . از منو، Component 1 (comp1)>Magnetic Fields>Magnetic>mf.normH – Norma میدان مغناطیسی – A/m را انتخاب کنید . |
3 | قسمت Coloring and Style را پیدا کنید . ![]() |
4 | در کادر محاوره ای Color Table ، Thermal>GrayBody را در درخت انتخاب کنید. |
5 | روی OK کلیک کنید . |
6 | ![]() |
7 | در کادر محاورهای Show More Options ، در درخت، کادر را برای گره Results>All Plot Types انتخاب کنید . |
8 | روی OK کلیک کنید . |
مسیرهای یون ها
در پنجره Model Builder ، روی Ions Trajectories کلیک کنید .
ردیابی ذرات با جرم 1
1 | در نوار ابزار Ions Trajectories ، روی ![]() |
عبارات نیروی لورنتس را که بر ذرات وارد می شود وارد کنید.
موقعیت ذره در امتداد محور z را می توان با استفاده از متغیر زمان ذره partt محاسبه کرد .
2 | در پنجره تنظیمات برای ردیابی ذرات با جرم ، بخش نیروی کل را پیدا کنید . |
3 | در قسمت متن Fx ، -q*vz*mf.By*(1-2*(partt>L1/vz)+2*(partt>(L1+L2)/vz)-(partt>(L1+L2) را تایپ کنید +L3)/vz)) . |
4 | در قسمت متن Fy ، q*vz*mf.Bx*(1-2*(partt>L1/vz)+2*(partt>(L1+L2)/vz)-(partt>(L1+L2+) را تایپ کنید L3)/vz)) . |
5 | برای گسترش بخش Mass and Velocity کلیک کنید . در قسمت Mass text، m را تایپ کنید . |
6 | قسمت Particle Positioning را پیدا کنید . در قسمت متن x ، 0.03*cos(range(0,0.05*pi,2*pi)) را تایپ کنید . |
7 | در قسمت متن y ، 0.03*sin(range(0,0.05*pi,2*pi)) را تایپ کنید . |
بیان رنگ 1
1 | روی Particle Tracing with Mass 1 کلیک راست کرده و Color Expression را انتخاب کنید . |
2 | در پنجره تنظیمات برای Color Expression ، بخش Expression را پیدا کنید . |
3 | در قسمت متن Expression ، 1+(partt>L1/vz)+(partt>(L1+L2)/vz)+(partt>(L1+L2+L3)/vz) را تایپ کنید . |
4 | قسمت Coloring and Style را پیدا کنید . ![]() |
5 | در کادر محاوره ای Color Table ، Rainbow>Cyclic را در درخت انتخاب کنید. |
6 | روی OK کلیک کنید . |
ردیابی ذرات با جرم 1
1 | در پنجره Model Builder ، روی Particle Tracing with Mass 1 کلیک کنید . |
2 | در پنجره تنظیمات برای ردیابی ذرات با جرم ، برای گسترش بخش انتشار کلیک کنید . |
3 | برای گسترش بخش Coloring and Style کلیک کنید . برای گسترش بخش کیفیت کلیک کنید . زیربخش تنظیمات حل کننده ODE را پیدا کنید . در قسمت متنی Relative tolerance ، 1e-6 را تایپ کنید . |
4 | برای گسترش بخش Advanced کلیک کنید . زیربخش خاتمه را پیدا کنید . |
5 | کادر بررسی حداکثر تعداد مراحل را انتخاب کنید . در قسمت متن مرتبط، 1e5 را تایپ کنید . |
6 | زیربخش میدان جریان آنی را پیدا کنید . |
7 | کادر بررسی زمان پایان را انتخاب کنید . در قسمت متن مرتبط، 5/3e6 را تایپ کنید . |
8 | در نوار ابزار Ions Trajectories ، روی ![]() |
کانتور 1
1 | در پنجره Model Builder ، روی Ions Trajectories کلیک راست کرده و Contour را انتخاب کنید . |
2 | در پنجره تنظیمات برای Contour ، بخش Expression را پیدا کنید . |
3 | در قسمت متن Expression ، sqrt(x^2+y^2) را تایپ کنید . |
4 | قسمت Levels را پیدا کنید . از لیست روش ورود ، سطوح را انتخاب کنید . |
5 | در قسمت متن Levels 0.01 0.03 را تایپ کنید . |
6 | قسمت Coloring and Style را پیدا کنید . از لیست Coloring ، Uniform را انتخاب کنید . |
7 | از لیست رنگ ، سفارشی را انتخاب کنید . |
8 | در ویندوز، روی نوار رنگی زیر کلیک کنید، یا – اگر از دسکتاپ چند پلتفرمی استفاده می کنید – روی دکمه Color کلیک کنید. |
9 | تیک Color legend را پاک کنید . |
10 | در نوار ابزار Ions Trajectories ، روی ![]() |
روی مرکز هندسه مدل بزرگنمایی کنید.