در پایان نامه حاضر معادلات الکتروشیمیایی حاکم بر پیل های سوختی غشاپلیمری در فضای دو بعدی مدل سازی شد و با استفاده از روش عددی المان مرزی در یک هندسه چند ناحیه ای فرمول بندی شد. این معادلات عموما در قالب کلی معادلات جابجایی-پخشی هستند که حل آنها به طور کلاسیک با استفاده از روش های دینامیک سیالات محاسباتی مانند روش تفاضلات محدود، روش حجم محدود و روش المان محدود صورت می گرفته است؛ اما فرمول بندی و حل معادلات حاکم بر پیل های سوختی برای اولین بار در این پژوهش توسط روش المان مرزی مورد مطالعه قرار گرفت. روش المان مرزی در سال های اخیر به عنوان یکی از ابزارهای بسیار مفید در تحلیل بسیاری از مسائل مهندسی شناخته شده است. عدم نیاز به گسسته سازی میدان محاسباتی و حل مسأله با استفاده از اطلاعات روی مرز مهم ترین ویژگی و نقطه قوت این روش به شمار می آید. با مطرح شدن روش تقابل دوگانه استفاده از این روش برای حل مسائل غیر خطی و همچنین مسائل خطی که درگیر عبارات انتگرالی روی میدان ناشی از عبارت چشمه می باشند اهمیت بیشتری یافته است. در این پژوهش چگونگی اعمال این روش برای حل معادلات حاکم بر پیل سوختی مورد بحث قرار گرفته است.