الکترود ضربان ساز

View Categories

الکترود ضربان ساز

7 min read

PDF

الکترود ضربان ساز
معرفی
این مثال استفاده از COMSOL Multiphysics را برای مدل‌سازی مشکلات توزیع جریان یونی در الکترولیت‌ها، در این مورد در بافت انسان، نشان می‌دهد. این مشکل در یک الکترود ضربان ساز نشان داده شده است، اما می توان آن را در سلول های الکتروشیمیایی مانند سلول های سوختی، باتری ها، حفاظت در برابر خوردگی یا هر فرآیند دیگری که در آن رسانش یونی در غیاب گرادیان غلظت انجام می شود، اعمال کرد.
دستگاه مدل سازی شده یک الکترود ضربان ساز است که در داخل قلب قرار می گیرد و به قلب بیمار کمک می کند تا ریتم طبیعی خود را حفظ کند. این دستگاه به عنوان یک الکترود شناخته می شود، اما در واقع از دو الکترود تشکیل شده است: یک کاتد و یک آند.
شکل 1 یک طرح شماتیک از دو جفت الکترود قرار گرفته در داخل قلب را نشان می دهد. الکترودها با جریان از واحد مولد پالس که در بیمار نیز کاشته می شود، تامین می شود.
شکل 1: ترسیم شماتیک قلب با دو جفت الکترود ضربان ساز.
این مثال به توزیع جریان و پتانسیل در اطراف یک جفت الکترود می پردازد.
تعریف مدل
دامنه مدل شامل خون و بافت اطراف جفت الکترود است. الکترودهای واقعی و تکیه گاه الکترود، مرزهای حوزه مدل شده هستند. شکل 2 الکترود را در سایه تیره‌تر نشان می‌دهد، در حالی که حوزه مدل‌سازی اطراف با سایه روشن‌تر نشان داده شده است.
شکل 2: مدل سازی دامنه و مرزها.
الکترود کار شامل یک نیمکره است که در نوک ساختار استوانه ای نگهدارنده قرار می گیرد. الکترود شمارنده در “کمر” این ساختار قرار می گیرد. تمام سطوح دیگر سازه نگهدارنده عایق بندی شده اند. مرزهای بیرونی به اندازه کافی دور از الکترود قرار می گیرند تا تأثیر کمی بر توزیع جریان و پتانسیل داشته باشند.
در COMSOL Multiphysics، از رابط Electric Currents برای تجزیه و تحلیل الکترود استفاده کنید. این رابط فیزیک برای مدلسازی مواد رسانا که در آن جریان به دلیل میدان الکتریکی اعمال شده جریان دارد مفید است.
معادلات دامنه
جریان در دامنه توسط معادله پیوستگی کنترل می شود که از معادلات ماکسول به دست می آید:
جایی که σ رسانایی بافت انسان است. این معادله از روابط زیر بین پتانسیل الکتریکی و میدان ها استفاده می کند.
شرایط مرزی
شرایط مرزی پتانسیل زمین بر روی کمر نازک‌تر الکترود اعمال می‌شود. نوک الکترود دارای پتانسیل ثابت 1  ولت است. تمام مرزهای دیگر عایق الکتریکی هستند.
نتایج و بحث
این شبیه سازی توزیع پتانسیل را بر روی سطح الکترود و توزیع جریان در قلب انسان را ساده می کند. شکل 3 را ببینید .
شکل 3: نمودار پتانسیل الکترواستاتیک توزیع شده روی سطح الکترود را نشان می دهد. چگالی جریان کل به صورت خطوط جریان نشان داده می شود.
همانطور که انتظار می رود، چگالی جریان در نیمکره کوچک، که باعث تحریک قلب می شود، بالاترین میزان است. چگالی جریان روی الکترود کار تقریباً یکنواخت است. الکترود شمارنده بزرگتر است و همچنین تغییرات بزرگتری در چگالی جریان روی سطح آن وجود دارد. به طور عمده، جریان با افزایش فاصله از الکترود کار کمتر می شود. مدل نشان می دهد که بازوهای لنگر دستگاه تأثیر کمی بر توزیع چگالی جریان دارند.
مسیر کتابخانه برنامه: COMSOL_Multiphysics/Electromagnetics/pacemaker_electrode
دستورالعمل مدلسازی
از منوی File ، New را انتخاب کنید .
جدید
در پنجره جدید ، روی  Model  Wizard کلیک کنید .
مدل جادوگر
1
در پنجره Model  Wizard ، روی  3D کلیک کنید .
2
در درخت Select  Physics ، AC/DC>Electric  Fields  و  Currents>Electric  Currents  (ec) را انتخاب کنید .
3
روی افزودن کلیک کنید .
4
 روی مطالعه کلیک کنید .
5
در درخت انتخاب  مطالعه ، General  Studies>Stationary را انتخاب کنید .
6
 روی Done کلیک کنید .
هندسه 1
دنباله هندسی را از فایل pacemaker_electrode_geom_sequence.mph وارد کنید .
1
در نوار ابزار Geometry ، روی Insert  Sequence کلیک کنید و Insert  Sequence را انتخاب کنید .
2
به پوشه Application Libraries مدل بروید و روی فایل pacemaker_electrode_geom_sequence.mph دوبار کلیک کنید .
3
در نوار ابزار Geometry ، روی  ساختن  همه کلیک کنید .
سپس حجم اطراف الکترود را مشخص کنید. شبیه سازی فقط در این حجم انجام می شود، جایی که مرزهای الکترود بر نتیجه تأثیر می گذارد.
سیلندر 1 (cyl1)
1
در نوار ابزار Geometry ، روی  Cylinder کلیک کنید .
2
در پنجره تنظیمات سیلندر ، بخش اندازه  و  شکل را پیدا کنید .
3
در قسمت نوشتار Radius ، 10[mm] را تایپ کنید .
4
در قسمت متن ارتفاع ، 40[mm] را تایپ کنید .
5
قسمت Position را پیدا کنید . در قسمت متن z ، -20[mm] را تایپ کنید .
تفاوت 1 (dif1)
1
در نوار ابزار Geometry ، روی  Booleans  and  Partitions کلیک کنید و Difference را انتخاب کنید .
2
فقط شی cyl1 را انتخاب کنید.
3
در پنجره تنظیمات برای تفاوت ، بخش تفاوت را پیدا کنید .
4
از لیست اشیاء  برای  تفریق ، الکترود را انتخاب کنید .
فرم اتحادیه (فین)
در نوار ابزار Geometry ، روی  ساختن  همه کلیک کنید .
تعاریف
سپس، یک انتخاب مربوط به الکترود زمین شده را برای استفاده بعدی تعریف کنید.
هندسه 1
شمارنده الکترود
1
در نوار ابزار Geometry ، روی  Selections کلیک کنید و Explicit  Selection را انتخاب کنید .
2
در پنجره تنظیمات برای انتخاب صریح  ، بخش Entities to Select را پیدا کنید .
3
از لیست سطح نهاد هندسی  ، Boundary را انتخاب کنید .
4
تیک گروه  بر اساس  مماس پیوسته را  انتخاب کنید .
5
 روی دکمه Wireframe  Rendering در نوار ابزار Graphics کلیک کنید .
6
در باله شی ، فقط مرزهای 29، 30، 58 و 63 را انتخاب کنید.
7
در قسمت نوشتار Label ، Counter Electrode را تایپ کنید .
8
 روی دکمه Wireframe  Rendering در نوار ابزار Graphics کلیک کنید تا حالت رندر پیش فرض بازیابی شود.
مواد
یک راه راحت برای فهمیدن اینکه کدام پارامتر مواد را باید مشخص کنید، اضافه کردن یک ماده است. سپس COMSOL Multiphysics هرگونه پارامتر از دست رفته را برای رابط های فیزیکی که به مدل اضافه کرده اید نشان می دهد.
بافت قلب
1
در پنجره Model  Builder ، در قسمت Component   (comp1) روی Materials راست کلیک کرده و Blank  Material را انتخاب کنید .
به طور پیش‌فرض، اولین متریالی که اضافه می‌کنید برای همه دامنه‌ها اعمال می‌شود، بنابراین نیازی به تغییر دامنه هندسی نیست.
الکترود به قلب انسان وارد می شود، بنابراین باید رسانایی بافت قلب را مشخص کنید.
2
در پنجره تنظیمات برای Material ، قسمت Material  Contents را پیدا کنید .
3
در جدول تنظیمات زیر را وارد کنید:
 
ویژگی
متغیر
ارزش
واحد
گروه اموال
رسانایی الکتریکی
sigma_iso ; sigmaii = sigma_iso، sigmaij = 0
0.4 [S/m]
S/m
پایه ای
گذر نسبی
epsilonr_iso ; epsilonrii = epsilonr_iso، epsilonrij = 0
1
1
پایه ای
4
در قسمت نوشتار برچسب ، بافت قلب را تایپ کنید .
جریان های الکتریکی (EC)
تنها تنظیمات فیزیک که باید مشخص شود پتانسیل های الکترود هستند.
زمین 1
1
در پنجره Model  Builder ، در قسمت Component   (comp1) روی Electric  Currents  (ec) کلیک راست کرده و Ground را انتخاب کنید .
2
در پنجره تنظیمات برای زمین ، قسمت انتخاب مرز  را پیدا کنید .
3
از لیست Selection ، Counter  Electrode را انتخاب کنید .
پتانسیل الکتریکی 1
1
در نوار ابزار Physics ، روی  Boundaries کلیک کنید و Electric  Potential را انتخاب کنید .
2
در پنجره تنظیمات برای پتانسیل الکتریکی  ، قسمت انتخاب مرز را پیدا کنید .
3
از لیست Selection ، Spherical  Electrode را انتخاب کنید .
4
قسمت Electric  Potential را پیدا کنید . در قسمت متنی 0 ، 1 را تایپ کنید .
مش 1
از مش پیش فرض کنترل شده با فیزیک استفاده کنید.
1
در پنجره Model  Builder ، در قسمت Component   (comp1) روی Mesh  1 کلیک راست کرده و Build  All را انتخاب کنید .
مطالعه 1
از تنظیمات پیش‌فرض برای حل‌کننده ثابت استفاده کنید، که به حل‌کننده گرادیان مزدوج با چندشبکه جبری به عنوان پیش‌شرطی‌کننده می‌دهد.
1
در نوار ابزار صفحه اصلی ،  روی محاسبه کلیک کنید .
نتایج
پتانسیل الکتریکی (EC)
نمودار پیش فرض برش های پتانسیل الکتریکی را نشان می دهد. برای بازتولید نمودار نشان داده شده در شکل 3 ، با پنهان کردن مرزهای بیرونی شروع کنید.
جریان های الکتریکی (EC)
1
در پنجره Model  Builder ، در قسمت Component   (comp1) روی Electric  Currents  (ec) کلیک کنید .
2
 روی دکمه کلیک  و  پنهان کردن در نوار ابزار Graphics کلیک کنید .
3
در نوار ابزار پنجره Graphics ، در کنار  Select  Domains کلیک کنید ، سپس Select  Boundaries را انتخاب کنید .
4
فقط مرزهای 1-4، 45 و 74 را انتخاب کنید.
5
 روی دکمه کلیک  و  پنهان کردن در نوار ابزار Graphics کلیک کنید .
نتایج
گروه سه بعدی پلات 3
1
در نوار ابزار صفحه اصلی ، روی  Add  Plot  Group کلیک کنید و 3D  Plot  Group را انتخاب کنید .
2
در پنجره Settings for 3D  Plot  Group ، قسمت Plot  Settings را پیدا کنید .
3
کادر بررسی لبه های مجموعه داده Plot را  پاک کنید .
سطح 1
1
روی 3D  Plot  Group  کلیک راست کرده و Surface را انتخاب کنید .
2
در نوار ابزار 3D Plot Group 3 ، روی  Plot کلیک کنید .
گروه سه بعدی پلات 3
نمودار سطح پتانسیل را با خطوط جریانی که چگالی جریان کل را تجسم می کنند ترکیب کنید.
ساده 1
1
در پنجره Model  Builder ، روی 3D  Plot  Group  3 کلیک راست کرده و Streamline را انتخاب کنید .
2
در پنجره تنظیمات برای Streamline ، بخش Selection را پیدا کنید .
3
از لیست انتخاب ، همه  مرزها را انتخاب کنید .
4
بخش Streamline  Positioning را پیدا کنید . از لیست موقعیت یابی ، نقطه شروع  کنترل شده را انتخاب کنید .
5
قسمت Coloring  and  Style را پیدا کنید . زیربخش Line  style را پیدا کنید . از لیست نوع ، لوله را انتخاب کنید .
6
در قسمت متن بیان شعاع لوله ،  min(ec.normJ/0.1[mA/mm^2],1)*0.2[mm] را تایپ کنید . ec.normJ متغیر هنجار چگالی جریان است. این عبارت بیان می کند که لوله ها 0.2[mm] عرض در نقاط دارای 0.1 [mA/mm^2] دارند. تابع “min” افزایش شعاع لوله را در جایی که چگالی جریان به 0.1 [mA/mm^2] نزدیک می‌شود، اشباع می‌کند.
برای به دست آوردن خطوط جریان ضخیم مناسب باید ضریب مقیاس را تنظیم کنید.
7
چک باکس Radius  scale  factor را انتخاب کنید .
8
در نوار ابزار 3D Plot Group 3 ، روی  Plot کلیک کنید .
بیان رنگ 1
1
روی Streamline  کلیک راست کرده و Color  Expression را انتخاب کنید .
از عبارت پیش فرض استفاده کنید. از آنجایی که این عبارت همان عبارت سطحی است، می توانید افسانه رنگ را غیرفعال کنید:
2
در پنجره تنظیمات برای بیان رنگ  ، قسمت Coloring and Style را پیدا کنید .
3
تیک Color  legend را پاک کنید .
پس از انجام عملیات چرخش و بزرگنمایی مناسب، باید چیزی شبیه به نمودار در شکل 3 مشاهده کنید .