انتقال گونه ها در لایه های انتشار گاز یک PEM

View Categories

انتقال گونه ها در لایه های انتشار گاز یک PEM

31 min read

PDF

انتقال گونه ها در لایه های انتشار گاز یک PEM
معرفی
این مثال بر حمل و نقل گونه ها در لایه های انتشار گاز (GDLs) یک سلول سوختی غشای تبادل پروتون (PEM) تمرکز دارد. هندسه یک سلول را با دو کانال جریان مجاور با فشارهای مختلف مدل می‌کند، وضعیتی که ممکن است در سلولی با کانال‌های جریان مارپیچ، یا در سلولی با استفاده از طراحی میدان جریان درهم رخ دهد. این مدل از ترازهای جریان، معادلات انتقال جرم (نشر ماکسول-استفان برای واکنش دهنده ها، آب و گاز نیتروژن)، و انتقال تکانه (قانون دارسی برای جریان گاز) برای شبیه سازی رفتار پیل سوختی PEM استفاده می کند.
تعریف مدل
بخش مدل‌سازی‌شده پیل سوختی از سه حوزه تشکیل شده است: یک آند ( Ω a )، یک غشای تبادل پروتون ( Ω m ) و یک کاتد ( Ωc ) همانطور که در شکل 1 نشان داده شده است .
شکل 1: هندسه مدل با برچسب های دامنه و مرز.
هر یک از الکترودها (لایه‌های انتشار گاز) در تماس با یک توزیع‌کننده گاز است که دارای یک کانال ورودی ( ∂Ω a، ورودی ) و یک کلکتور جریان ( ∂Ω a,cc ) است. از همین نماد برای سمت کاتد استفاده می شود.
هیدروژن و هوای مرطوب شده به ترتیب به کانال های ورودی آند و کاتد عرضه می شوند. هیدروژن واکنش نشان می دهد و در لایه فعال آندی مصرف می شود تا پروتون هایی تشکیل دهد که جریان یونی را به کاتد منتقل می کنند.
(1)
در کاتد، اکسیژن همراه با پروتون ها واکنش می دهد تا آب را در لایه فعال بر اساس:
(2)
هر دو گاز خوراک (هیدروژن مرطوب و هوای مرطوب) به عنوان ایده آل در نظر گرفته می شوند و از طریق انتشار و همرفت منتقل می شوند. الکترودها به عنوان محیط متخلخل همگن با خواص مورفولوژیکی یکنواخت مانند تخلخل و نفوذپذیری درمان می شوند. گاز درون هر یک از الکترودها به صورت یک فاز پیوسته وجود دارد، بنابراین قانون دارسی اعمال می شود.
یک مدل آگلومرا واکنش های الکترود را در لایه های کاتالیست فعال توصیف می کند. آگلومراها متشکل از کاتالیزور و ذرات کربن هستند که در الکترولیت پلیمری تعبیه شده اند. معادلات مدل آگلومرا از حل تحلیلی یک مسئله انتشار-واکنش در یک ذره متخلخل کروی سرچشمه می‌گیرد ( مرجع 5 و مرجع 6 ). در لایه کاتالیزور فعال آندی، هیدروژن گونه انتشار دهنده و واکنش دهنده در آگلومراها است، در حالی که اکسیژن گونه انتشار دهنده و واکنش دهنده در آگلومراها در کاتد است. یک مدل آگلومرا از لایه کاتالیست فعال کاتدی یک پیل سوختی PEM توسط بروکا و دیگران ارائه شده است ( مرجع 7 و مرجع 8 ).
موجودی شارژ
رابط پیل سوختی هیدروژنی برای مدل‌سازی توزیع‌های پتانسیل در سه حوزه، با چگالی جریان انتقال بار لایه‌های فعال به عنوان شرایط مرزی داخلی، با استفاده از گره‌های الکترود انتشار گاز نازک H2/O2 استفاده می‌شود. عبارات چگالی جریان محلی برای آند و کاتد، a و c ، با توجه به معادله مشخص می شوند.
در جایی که شاخص e مخفف “a” (آند) یا “c” (کاتد) است، agg شعاع آگلومرا (واحد SI: m) و agg,a و agg,c (واحد SI: A/m 3) است. ) چگالی جریان داده شده توسط مدل آگلومرا هستند. علاوه بر این، act ضخامت لایه فعال (الکترود انتشار گاز نازک) (m) و ε mac تخلخل آن (تخلخل ماکروسکوپی) است. سطح ویژه الکترود نازک انتشار گاز (واحد SI: 1/m) مطابق معادله زیر مشخص می شود.
مدل آگلومره برای آند و کاتد
مدل آگلومرا چگالی جریان را در یک لایه فعال متشکل از آگلومراهای مواد هادی یونی و ذرات رسانای الکتریکی که تا حدی با کاتالیزور پوشانده شده است، توصیف می کند. چگالی جریان محلی را می توان با حل ترکیبی از معادله انتشار و معادله جنبشی الکترود باتلر-ولمر برای یک آگلومرا با پتانسیل الکتریکی و یونی ثابت به صورت تحلیلی بیان کرد. معادلات حاصل برای چگالی جریان در آند و کاتد عبارتند از ( مراجعه 7 )
جایی که، دوباره، شاخص e مخفف “a” (آند) یا “c” (کاتد) است، و
در این معادلات، agg انتشار گاز آگلومرا است (واحد SI: m 2 / s)، e یک عدد “انتقال بار” است (1 برای آند و -2 برای کاتد)، S منطقه ویژه است. کاتالیزور داخل آگلومرا (واحد SI: 1/m) و F ثابت فارادی است (واحد SI: C/mol). علاوه بر این، i، ref غلظت مرجع گونه است (واحد SI: mol/m 3 )، ، agg غلظت های مربوطه در سطح آگلومرا هستند (واحد SI: mol/m 3 )، 0a و 0cچگالی جریان مبادله (واحد SI: A/m 2 )، R ثابت گاز، T دما (واحد SI: K)، و اضافه ولتاژها در آند و کاتد توسط
که در آن Eq (واحد SI: V) نشان دهنده ولتاژ تعادل است .
شما حالت مرجع آندی و کاتدی را برابر با کسر مولی در کانال های ورودی آند و کاتد به ترتیب برابر با 1 atm تنظیم می کنید. غلظت هیدروژن و اکسیژن محلول در سطح آگلومراها با کسر مولی گونه مربوطه در فاز گاز از طریق قانون هنری مرتبط است.
که در آن K ثابت هنری است (واحد SI: Pa·m 3 /mol).
موجودی شارژ، ادامه دارد
اختلاف پتانسیل بین کلکتورهای جریان کاتد و آند با ولتاژ کل سلول مطابقت دارد. با قرار دادن آن روی صفر، پتانسیل موجود در جمع کننده جریان آند را به عنوان سطح مرجع انتخاب کنید. سپس ولتاژ کل سلول به عنوان شرایط مرزی در کلکتور جریان کاتد عمل می کند:
برای سایر مرزها شرایط مرزی عایق الکتریکی دارید.
جریان سیال رسانه متخلخل
برای مدل‌سازی جریان گاز در پشت‌های گاز، این مثال از قانون دارسی در رابط پیل سوختی هیدروژنی استفاده می‌کند. سرعت گاز با معادله پیوستگی مطابق با داده می شود
که ρ چگالی مخلوط فاز گاز (واحد SI: kg/m 3 ) و u نشان دهنده سرعت گاز (واحد SI: m/s) است. قانون دارسی برای محیط های متخلخل بیان می کند که گرادیان فشار، ویسکوزیته سیال و ساختار محیط متخلخل سرعت را تعیین می کند:
در اینجا p نشان‌دهنده نفوذپذیری الکترود است (واحد SI: m2 ) ، η نشان‌دهنده ویسکوزیته گاز (واحد SI: Pa·s)، و p فشار است (واحد SI: Pa). قانون گاز ایده آل چگالی مخلوط فاز گاز را می دهد، ρ :
(3)
در این معادله، R ثابت گاز (واحد SI: J/(mol·K))، T دما (واحد SI: K)، M جرم مولی (واحد SI: kg/mol) و x است. کسر مول
در ورودی ها فشار را به صورت زیر مشخص می کنید:
در مرز نازک الکترود انتشار گاز برای آند و کاتد، سرعت گاز به طور خودکار توسط رابط پیل سوختی هیدروژن، از جریان جرم کل داده شده توسط نرخ واکنش الکتروشیمیایی و ضرایب استوکیومتری معادله 1 و معادله 2، با استفاده از فارادی محاسبه می شود . قانون
حمل و نقل انبوه ماکسول-استفان
این مدل دو گونه را در آند – H 2 و H 2 O – و سه گونه در کاتد – O 2 ، H 2 O و N 2 را در نظر می گیرد . رابط پیل سوختی هیدروژنی از انتشار چند جزئی Maxwell-Stefan استفاده می کند که توسط معادلات کنترل می شود.
برای حل کسرهای جرمی، i . این مدل پیل سوختی PEM خاص فرض می‌کند که انتشار ناشی از دما ناچیز است و عبارت منبع R را صفر می‌کند. ضرایب انتشار باینری Maxwell-Stefan، ij (واحد SI: m2 / s)، به طور خودکار توسط رابط محاسبه می شود.
کسرهای مول خوراک گاز در ورودی ها مشخص می شوند. در مرزهای نازک الکترود انتشار گاز، شار جرمی گونه ها به طور خودکار از روی نرخ واکنش الکتروشیمیایی و ضرایب استوکیومتری با استفاده از قانون فارادی تعیین می شود. کمک های سرعت استفان نیز به طور خودکار توسط رابط محاسبه می شود.
ویژگی‌های حمل و نقل غشایی رابط پیل سوختی هیدروژنی برای مدل‌سازی انتقال آب در فاز یونومر در حوزه غشاء استفاده می‌شود. شار مولکولی آب با استفاده از پارامترهای تخمینی تجربی موجود در مواد Nafion در کتابخانه مواد سلول سوختی و الکترولیزگر، هم به نفوذ شیب پتانسیل شیمیایی و هم به نیروی کشش الکترواسموتیک بستگی دارد.
نتایج و بحث
شکل 2 توزیع جریان را در پیل سوختی PEM نشان می دهد. در گوشه های کلکتورهای فعلی، نوک های جریان قابل توجهی وجود دارد.
شکل 2: چگالی جریان (نقشه سطحی) و میدان برداری جریان (نقشه پیکان) در پیل سوختی که با ولتاژ 0.7 ولت کار می کند. آند در سمت چپ و کاتد در سمت راست است.
برای تجزیه و تحلیل بیشتر رفتار سلول، چگالی جریان در لایه فعال را به عنوان تابعی از ارتفاع سلول ( y ) مانند شکل 3 رسم کنید .
شکل 3: توزیع چگالی جریان در لایه فعال در آند.
چگالی جریان ناهموار با بیشترین چگالی در ناحیه بالایی سلول است. این بدان معنی است که سرعت واکنش کاهش اکسیژن در کاتد توزیع چگالی جریان را تعیین می کند. حداکثر چگالی جریان نزدیک به ورودی هوا ایجاد می شود.
شارهای همرفتی به طور کلی بر حمل و نقل جرم در سلول غالب هستند. برای مطالعه اثرات همرفتی، میدان سرعت را مانند شکل 4 رسم کنید .
شکل 4: میدان سرعت گاز در محفظه آند و کاتد.
بزرگی سرعت جریان در گوشه های کلکتور فعلی به بالاترین مقدار خود می رسد.
شکل 5 کسر وزنی واکنش دهنده (اکسیژن و هیدروژن) را در گازهای کاتد و آند نشان می دهد.
شکل 5: کسر جرمی واکنش دهنده در سمت آند (چپ) و سمت کاتد (راست). واکنش دهنده در آند هیدروژن و در کاتد اکسیژن است.
اگرچه مصرف اکسیژن کم است، اضافه ولتاژ غلظت در آگلومراها سهم قابل توجهی در اضافه ولتاژ غلظت دارد. یک تغییر کوچک در جریان اکسیژن یک تغییر اساسی در قطبش سلول ایجاد می کند.
شکل 6 کسر جرمی آب را در گازهای آند و کاتد نشان می دهد. آب از طریق انتشار و همرفت به غشا در سمت آند منتقل می شود. اگر گاز آند بیش از حد خشک شود، غشاء خشک می‌شود و در نتیجه هدایت یونی کاهش می‌یابد و سلول متعاقباً از کار می‌افتد. از طرف دیگر، در سمت کاتد، قطرات آب می‌توانند منافذ را مسدود کرده و به طور موثر مانع انتقال گاز به لایه فعال شوند.
شکل 6: کسر جرمی آب در آند (چپ، نمودار سطحی) و کاتد (راست، نمودار کانتور).
منابع
1. W. He، JS Yi و TV Nguyen، “مدل جریان دو فازی کاتد پیل‌های سوختی PEM با استفاده از میدان‌های جریان درهم،” AIChe J. ، جلد. 46، صفحات 2053-2063، 2000.
2. C. Marr and X. Li, “Composition and Performance Modeling of Catalyst Layer in a Proton Exchange Membrane Fuel Cell,” J. Power Sources , vol. 77، صفحات 17-27، 1999.
3. P. Futerko و I.-M. Hsing، “مطالعه روش المان محدود دوبعدی مقاومت غشاها در سلولهای سوختی الکترولیت پلیمری،” Electrochimica Acta ، جلد. 45، صفحات 1741-1751، 2000.
4. DM Bernardi و MW Verbrugge، “مدل ریاضی یک الکترود انتشار گاز متصل به یک الکترولیت پلیمری”، AIChe J. ، جلد. 37، صص 1151-1163، 1991.
5. H. Scott Fogler, Elements of Chemical Reaction Engineering , ed. 3, Prentice Hall, 1999.
6. RB Bird، WE Stewart، و EN Lightfoot، پدیده های حمل و نقل ، جان وایلی و پسران، 1960.
7. K. Broka و P. Ekdunge، “مدل سازی کاتد پیل سوختی PEM”، J. Appl. الکتروشیمی. ، جلد 27، صفحات 281-289، 1997.
8. K. Dannenberg، P. Ekdunge، و G. Lindbergh، “مدل ریاضی PEMFC”، J. Appl. Electrochem ., vol. 30، صص 1377–1387، 2000.
مسیر کتابخانه برنامه: Fuel_Cell_and_Electrolyzer_Module/Fuel_Cells/pem_gdl_species_transport_2d
دستورالعمل های مدل سازی
از منوی File ، New را انتخاب کنید .
جدید
در پنجره جدید ، روی  Model  Wizard کلیک کنید .
مدل جادوگر
1
در پنجره Model  Wizard روی  2D کلیک کنید .
2
در درخت Select  Physics ، Electrochemistry>Hydrogen  Fuel  Cells>Proton  Exchange  (fc) را انتخاب کنید .
3
روی افزودن کلیک کنید .
4
 روی مطالعه کلیک کنید .
5
در درخت Select  Study ، Preset  Studies  for  Selected  Physics  Interfaces>Stationary  with  Initialization را انتخاب کنید .
6
 روی Done کلیک کنید .
هندسه 1
هندسه را با استفاده از مستطیل ایجاد کنید.
مستطیل 1 (r1)
1
در نوار ابزار Geometry ، روی  Rectangle کلیک کنید .
2
در پنجره تنظیمات برای Rectangle ، بخش Size  and  Shape را پیدا کنید .
3
در قسمت متن Width ، 2.5e-4 را تایپ کنید .
4
در قسمت متن ارتفاع ، 2e-3 را تایپ کنید .
مستطیل 2 (r2)
1
در نوار ابزار Geometry ، روی  Rectangle کلیک کنید .
2
در پنجره تنظیمات برای Rectangle ، بخش Size  and  Shape را پیدا کنید .
3
در قسمت متن Width ، 1e-4 را تایپ کنید .
4
در قسمت متن ارتفاع ، 2e-3 را تایپ کنید .
5
قسمت Position را پیدا کنید . در قسمت متن x ، 2.5e-4 را تایپ کنید .
مستطیل 3 (r3)
1
در نوار ابزار Geometry ، روی  Rectangle کلیک کنید .
2
در پنجره تنظیمات برای Rectangle ، بخش Size  and  Shape را پیدا کنید .
3
در قسمت متن Width ، 2.5e-4 را تایپ کنید .
4
در قسمت متن ارتفاع ، 2e-3 را تایپ کنید .
5
قسمت Position را پیدا کنید . در قسمت متن x ، 3.5e-4 را تایپ کنید .
مستطیل 4 (r4)
1
در نوار ابزار Geometry ، روی  Rectangle کلیک کنید .
2
در پنجره تنظیمات برای Rectangle ، بخش Size  and  Shape را پیدا کنید .
3
در قسمت متن Width ، 1e-4 را تایپ کنید .
4
در قسمت متن ارتفاع ، 5e-4 را تایپ کنید .
5
قسمت Position را پیدا کنید . در قسمت متن x ، -1e-4 را تایپ کنید .
آرایه 1 (arr1)
1
در نوار ابزار Geometry ، روی  Transforms کلیک کنید و Array را انتخاب کنید .
2
فقط شی r4 را انتخاب کنید.
3
در پنجره تنظیمات آرایه ، بخش Size را پیدا کنید .
4
در قسمت متن  اندازه 2 را تایپ کنید .
5
در قسمت متنی اندازه y ،  2 را تایپ کنید .
6
قسمت Displacement را پیدا کنید . در قسمت متن x ، 7e-4 را تایپ کنید .
7
در قسمت متن y ، 1.5e-3 را تایپ کنید .
اتحادیه 1 (uni1)
1
در نوار ابزار Geometry ، روی  Booleans  and  Partitions کلیک کنید و Union را انتخاب کنید .
2
فقط اشیاء arr1(1،1) ، arr1(1،2) و r1 را انتخاب کنید.
3
در پنجره تنظیمات برای Union ، بخش Union را پیدا کنید .
4
کادر تیک Keep  interior  borders را پاک کنید .
اتحادیه 2 (uni2)
1
در نوار ابزار Geometry ، روی  Booleans  and  Partitions کلیک کنید و Union را انتخاب کنید .
2
فقط اشیاء arr1(2،1) ، arr1(2،2) و r3 را انتخاب کنید.
3
در پنجره تنظیمات برای Union ، بخش Union را پیدا کنید .
4
کادر تیک Keep  interior  borders را پاک کنید .
فیله 1 (fil1)
1
در نوار ابزار Geometry ، روی  Fillet کلیک کنید .
2
در پنجره تنظیمات برای Fillet ، بخش Points را پیدا کنید .
3
زیربخش Vertices  to  fillet را پیدا کنید . برای پاک کردن دکمه  فعالسازی  انتخاب کلیک کنید.
4
 روی دکمه Zoom  Extents در نوار ابزار Graphics کلیک کنید .
5
در شی uni1 ، فقط نقاط 6 و 7 را انتخاب کنید.
6
در شی uni2 ، فقط نقاط 4 و 5 را انتخاب کنید.
7
برای انتخاب دکمه ضامن  فعال کردن  انتخاب کلیک کنید .
8
در شی uni1 ، فقط نقاط 6 و 7 را انتخاب کنید.
9
در شی uni2 ، فقط نقاط 4 و 5 را انتخاب کنید.
10
بخش Radius را پیدا کنید . در قسمت متن Radius ، 1e-5 را تایپ کنید .
11
در نوار ابزار Geometry ، روی  ساختن  همه کلیک کنید .
تعاریف جهانی
به بارگیری مجموعه ای از پارامترهای مدل جهانی از یک فایل متنی ارائه شده با Application Library ادامه دهید.
پارامترهای 1
1
در پنجره Model  Builder ، در قسمت Global  Definitions روی Parameters  1 کلیک کنید .
2
در پنجره تنظیمات برای پارامترها ، بخش پارامترها را پیدا کنید .
3
 روی Load  from  File کلیک کنید .
4
به پوشه Application Libraries مدل بروید و روی فایل pem_gdl_species_transport_2d_parameters.txt دوبار کلیک کنید .
تعاریف
تعدادی انتخاب را روی هندسه تنظیم کنید. اینها بعداً هنگام تنظیم فیزیک استفاده خواهند شد.
آند GDL
1
در نوار ابزار تعاریف ، روی  Explicit کلیک کنید .
2
فقط دامنه 1 را انتخاب کنید.
3
در پنجره تنظیمات برای Explicit ، Anode GDL را در قسمت نوشتار Label تایپ کنید .
غشاء
1
در نوار ابزار تعاریف ، روی  Explicit کلیک کنید .
2
فقط دامنه 2 را انتخاب کنید.
3
در پنجره تنظیمات برای Explicit ، در قسمت نوشتار Label، Membrane را تایپ کنید .
کاتد GDL
1
در نوار ابزار تعاریف ، روی  Explicit کلیک کنید .
2
فقط دامنه 3 را انتخاب کنید.
3
در پنجره تنظیمات برای Explicit ، Cathode GDL را در قسمت نوشتار Label تایپ کنید .
آند GDE
1
در نوار ابزار تعاریف ، روی  Explicit کلیک کنید .
2
در پنجره تنظیمات برای Explicit ، بخش Input  Entities را پیدا کنید .
3
از لیست سطح نهاد هندسی  ، Boundary را انتخاب کنید .
4
فقط مرز 10 را انتخاب کنید.
5
در قسمت Label text، Anode GDE را تایپ کنید .
کاتد GDE
1
در نوار ابزار تعاریف ، روی  Explicit کلیک کنید .
2
در پنجره تنظیمات برای Explicit ، بخش Input  Entities را پیدا کنید .
3
از لیست سطح نهاد هندسی  ، Boundary را انتخاب کنید .
4
فقط مرز 13 را انتخاب کنید.
5
در قسمت نوشتار Label ، Cathode GDE را تایپ کنید .
متغیرهای آند GDL
تعدادی متغیر مخصوص دامنه و مرز را تعریف کنید. اینها هم برای تنظیم فیزیک و هم در پس پردازش استفاده خواهند شد.
1
در نوار ابزار تعاریف ، روی متغیرهای  محلی  کلیک کنید .
2
در پنجره تنظیمات برای متغیرها ، Anode GDL Variables را در قسمت نوشتار Label تایپ کنید .
3
قسمت انتخاب موجودیت هندسی  را پیدا کنید . از لیست سطح نهاد هندسی ، دامنه را انتخاب کنید .
4
از لیست انتخاب ، Anode  GDL را انتخاب کنید .
5
قسمت Variables را پیدا کنید . در جدول تنظیمات زیر را وارد کنید:
 
نام
اصطلاح
واحد
شرح
واکنش نشان دهید
fc.wH2
کسر جرمی، واکنش دهنده
متغیرهای کاتدی GDL
1
در نوار ابزار تعاریف ، روی متغیرهای  محلی  کلیک کنید .
2
در پنجره تنظیمات برای متغیرها ، متغیرهای Cathode GDL را در قسمت متن برچسب تایپ کنید .
3
قسمت انتخاب موجودیت هندسی  را پیدا کنید . از لیست سطح نهاد هندسی ، دامنه را انتخاب کنید .
4
از لیست انتخاب ، Cathode  GDL را انتخاب کنید .
5
قسمت Variables را پیدا کنید . در جدول تنظیمات زیر را وارد کنید:
 
نام
اصطلاح
واحد
شرح
واکنش نشان دهید
fc.wO2
کسر جرمی، واکنش دهنده
متغیرهای آند GDE
1
در نوار ابزار تعاریف ، روی متغیرهای  محلی  کلیک کنید .
2
در پنجره تنظیمات برای متغیرها ، Anode GDE Variables را در قسمت نوشتار Label تایپ کنید .
3
قسمت انتخاب موجودیت هندسی  را پیدا کنید . از لیست سطح نهاد هندسی ، Boundary را انتخاب کنید .
4
از لیست انتخاب ، Anode  GDE را انتخاب کنید .
5
قسمت Variables را پیدا کنید . در جدول تنظیمات زیر را وارد کنید:
 
نام
اصطلاح
واحد
شرح
cH2_agg
fc.p*fc.xH2/KH2
قانون هانری غلظت آگلومرا هیدروژن
eta_a
fc.phis-fc.phil-E_eq_a
V
مازاد پتانسیل آندی
beta_a
cH2_agg-cH2_ref*exp(-2*F_const*eta_a/(R_const*T))
lda_a
sqrt(i0_a*S*R_agg^2/(2*F_const*cH2_ref*D_agg))
زیر بیان چگالی جریان آندی
i_a
K*(1-lda_a*coth(lda_a))*beta_a*(R_agg/3)
چگالی جریان آند
متغیرهای GDE کاتدی
1
در نوار ابزار تعاریف ، روی متغیرهای  محلی  کلیک کنید .
2
در پنجره تنظیمات متغیرها ، Cathode GDE Variables را در قسمت نوشتار Label تایپ کنید .
3
قسمت انتخاب موجودیت هندسی  را پیدا کنید . از لیست سطح نهاد هندسی ، Boundary را انتخاب کنید .
4
از لیست انتخاب ، Cathode  GDE را انتخاب کنید .
5
قسمت Variables را پیدا کنید . در جدول تنظیمات زیر را وارد کنید:
 
نام
اصطلاح
واحد
شرح
cO2_agg
fc.p*fc.xO2/KO2
قانون هانری غلظت آگلومره اکسیژن
eta_c
fc.phis-fc.phil-E_eq_c
V
اضافه ولتاژ کاتدی
lda_c
sqrt(i0_c*S*R_agg^2*exp(-F_const*eta_c/(2*R_const*T))/(4*F_const*cO2_ref*D_agg))
زیر بیان چگالی جریان کاتدی
مدار مجتمع
-2*K*(1-lda_c*coth(lda_c))*cO2_agg*(R_agg/3)
چگالی جریان کاتدی
مواد
این مدل از مواد الکترولیت پلیمری (Nafion 1100، با بخار متعادل) استفاده می کند که در کتابخانه مواد موجود است. مواد را دو بار اضافه کنید و ابتدا آن را به حوزه غشاء و در کنار مرزهای غشاء مجاور گره های فاز گاز اختصاص دهید.
مواد را اضافه کنید
1
در نوار ابزار Home ، روی  Add  Material کلیک کنید تا پنجره Add  Material باز شود .
2
به پنجره Add  Material بروید .
3
در درخت، Fuel  Cell  and  Electrolyzer>Polymer  Electrolytes>Nafion,  EW  1100,  Vapor  Equilibrated,  Protonated را انتخاب کنید .
4
کلیک راست کرده و Add  to  Component   (comp1) را انتخاب کنید .
مواد
Nafion، EW 1100، بخار متعادل، پروتونه (mat1)
1
در پنجره Model  Builder ، در قسمت Component   (comp1)>Materials روی Nafion،  EW  1100،  Vapor  Equilibrated،  Protonated  (mat1) کلیک کنید .
2
در پنجره تنظیمات برای مواد ، قسمت انتخاب موجودیت هندسی  را پیدا کنید .
3
 روی Clear  Selection کلیک کنید .
4
فقط دامنه 2 را انتخاب کنید.
مواد را اضافه کنید
1
به پنجره Add  Material بروید .
2
در درخت، Fuel  Cell  and  Electrolyzer>Polymer  Electrolytes>Nafion,  EW  1100,  Vapor  Equilibrated,  Protonated را انتخاب کنید .
3
کلیک راست کرده و Add  to  Component   (comp1) را انتخاب کنید .
4
در نوار ابزار Home ، روی  Add  Material کلیک کنید تا پنجره Add  Material بسته شود .
مواد
Nafion، EW 1100، بخار متعادل، پروتونه 1 (mat2)
1
در پنجره تنظیمات برای مواد ، قسمت انتخاب موجودیت هندسی  را پیدا کنید .
2
از لیست سطح نهاد هندسی  ، Boundary را انتخاب کنید .
3
فقط مرزهای 10 و 13 را انتخاب کنید.
پیل سوختی هیدروژنی (FC)
مدل توزیع و حمل و نقل فعلی را تنظیم کنید. شامل انتقال جرم با استفاده از انتشار ماکسول-استفان و انتقال تکانه با استفاده از قانون دارسی در هر دو مخلوط گاز آند و کاتد، همراه با کشیدن آب الکترواسموتیک در غشاء. همچنین سطح فشار مرجع را در خصوصیات رابط تعریف کنید. توجه داشته باشید که گونه های گاز پیش فرض هیدروژن و آب در سمت آند و اکسیژن، نیتروژن و آب در سمت کاتد هستند. با افزودن گره های دامنه مربوطه شروع کنید.
1
در پنجره Model  Builder ، در قسمت Component   (comp1) روی پیل سوختی هیدروژنی  (fc) کلیک کنید .
2
در پنجره تنظیمات پیل سوختی هیدروژنی  ، قسمت مخلوط گاز H2 را پیدا کنید .
3
زیربخش مکانیزم های حمل و نقل را  پیدا کنید . کادر بررسی استفاده از قانون دارسی برای انتقال حرکت را انتخاب کنید .
4
قسمت O2  Gas  Mixture را پیدا کنید . کادر بررسی استفاده  از قانون دارسی  برای انتقال حرکت را انتخاب کنید .
5
کلیک کنید تا قسمت Membrane  Transport گسترش یابد . چک باکس Electroosmotic  water  drag را انتخاب کنید .
6
برای گسترش بخش Reference  Pressure  Level کلیک کنید . در قسمت متن ref عدد 0 را تایپ کنید .
غشاء 1
1
در نوار ابزار Physics ، روی  Domains کلیک کنید و Membrane را انتخاب کنید .
2
در پنجره Settings for Membrane ، قسمت Domain  Selection را پیدا کنید .
3
از لیست انتخاب ، غشاء را انتخاب کنید .
لایه انتشار گاز H2 1
1
در نوار ابزار Physics ، روی  Domains کلیک کنید و H2  Gas  Diffusion  Layer را انتخاب کنید .
2
در پنجره تنظیمات برای لایه انتشار گاز H2  ، بخش انتخاب دامنه را پیدا کنید .
3
از لیست انتخاب ، Anode  GDL را انتخاب کنید .
لایه انتشار گاز O2 1
1
در نوار ابزار Physics ، روی  Domains کلیک کنید و O2  Gas  Diffusion  Layer را انتخاب کنید .
2
در پنجره تنظیمات برای لایه انتشار گاز O2  ، بخش انتخاب دامنه را پیدا کنید .
3
از لیست انتخاب ، Cathode  GDL را انتخاب کنید .
الکترود نازک انتشار گاز H2 1
1
در نوار ابزار Physics ، روی  Boundaries کلیک کنید و Thin  H2  Gas  Diffusion  Electrode را انتخاب کنید .
2
در پنجره تنظیمات برای الکترود نازک  انتشار گاز H2  ، بخش انتخاب مرز را پیدا کنید .
3
از لیست انتخاب ، Anode  GDE را انتخاب کنید .
الکترود نازک انتشار گاز O2 1
1
در نوار ابزار Physics ، روی  Boundaries کلیک کنید و Thin  O2  Gas  Diffusion  Electrode را انتخاب کنید .
2
در پنجره تنظیمات برای الکترود نازک  انتشار گاز O2  ، بخش انتخاب مرز را پیدا کنید .
3
از لیست انتخاب ، Cathode  GDE را انتخاب کنید .
در گره فاز الکترولیت ، هدایت الکترولیت تنظیم شده است که از گره مواد گرفته شود . تنظیمات را در گره های فاز گاز H2 و فاز گاز O2 بررسی کنید . توجه داشته باشید که چگالی و ویسکوزیته مخلوط گاز و ضرایب انتشار باینری زمانی که تنظیمات پیش فرض مربوطه استفاده می شود به طور خودکار محاسبه می شوند.
ویژگی‌های کشیدن آب الکترواسموتیک در گره غشایی و در گره‌های فرزند که به‌طور پیش‌فرض اضافه شده‌اند، به‌طور خودکار از گره مواد گرفته می‌شوند .
مقادیر اولیه 1
1
در پنجره Model  Builder ، گره Membrane  1 را گسترش دهید ، سپس روی مقادیر اولیه  1 کلیک کنید .
2
در پنجره تنظیمات برای مقادیر اولیه  ، قسمت مقادیر اولیه را پیدا کنید .
3
در قسمت متن 0 ، T را تایپ کنید .
لایه انتشار گاز H2 1
ویژگی های گره های لایه انتشار گاز H2 و لایه انتشار گاز O2 را تنظیم کنید .
1
در پنجره Model  Builder ، در قسمت Component   (comp1)>Hydrogen  Fuel  Cell  (fc) روی H2  Gas  Diffusion  Layer  1 کلیک کنید .
2
در پنجره تنظیمات برای لایه انتشار گاز H2  ، بخش حمل و نقل شارژ الکترود را پیدا کنید .
3
در قسمت متن σ s ، kappa_s را تایپ کنید .
4
بخش حمل و نقل گاز  را پیدا کنید . در قسمت متن ε g ، eps_mac را تایپ کنید .
5
در قسمت متن κ g ، kappa_p را تایپ کنید .
لایه انتشار گاز O2 1
1
در پنجره Model  Builder ، روی O2  Gas  Diffusion  Layer  1 کلیک کنید .
2
در پنجره تنظیمات برای لایه انتشار گاز O2  ، بخش حمل و نقل شارژ الکترود را پیدا کنید .
3
در قسمت متن σ s ، kappa_s را تایپ کنید .
4
بخش حمل و نقل گاز  را پیدا کنید . در قسمت متن ε g ، eps_mac را تایپ کنید .
5
در قسمت متن κ g ، kappa_p را تایپ کنید .
الکترود نازک انتشار گاز H2 1
ویژگی های گره الکترود انتشار گاز نازک H2 را تنظیم کنید . جزئیات سینتیک الکترود در گره فرزند تنظیم می شود.
1
در پنجره Model  Builder ، روی Thin  H2  Gas  Diffusion  Electrode  1 کلیک کنید .
2
در پنجره تنظیمات برای الکترود نازک  انتشار گاز H2  ، بخش ضخامت الکترود را پیدا کنید .
3
در قسمت متنی gde ، l_act را تایپ کنید .
واکنش الکترود انتشار گاز نازک H2 1
1
در پنجره Model  Builder ، روی Thin  H2  Gas  Diffusion  Electrode  Reaction  1 کلیک کنید .
2
در پنجره تنظیمات برای واکنش الکترود انتشار گاز نازک  H2  ، بخش پتانسیل تعادل را پیدا کنید .
3
از لیست Eq ، User defined را انتخاب کنید . با توجه به بیان سینتیک تعریف شده توسط کاربر، پارامتر پتانسیل تعادل در این مدل استفاده نمی شود. بنابراین مقدار پیش فرض 0 ولت را حفظ کنید.
4
بخش سینتیک الکترود  را پیدا کنید . از لیست loc,expr ، User defined را انتخاب کنید . در قسمت متن مرتبط، i_a را تایپ کنید .
5
قسمت Active  Specific  Surface  Area را پیدا کنید . در قسمت متن v ، Av را تایپ کنید .
الکترود نازک انتشار گاز O2 1
به طور مشابه، ویژگی های گره الکترود انتشار گاز نازک O2 را تنظیم کنید . جزئیات سینتیک الکترود در گره فرزند تنظیم می شود.
1
در پنجره Model  Builder ، در قسمت Component   (comp1)>Hydrogen  Fuel  Cell  (fc) روی Thin  O2  Gas  Diffusion  Electrode  1 کلیک کنید .
2
در پنجره تنظیمات برای الکترود انتشار گاز نازک  O2  ، بخش ضخامت الکترود را پیدا کنید .
3
در قسمت متنی gde ، l_act را تایپ کنید .
واکنش الکترود انتشار گاز نازک O2 1
1
در پنجره Model  Builder ، روی Thin  O2  Gas  Diffusion  Electrode  Reaction  1 کلیک کنید .
2
در پنجره تنظیمات برای واکنش الکترود انتشار گاز نازک  O2  ، بخش پتانسیل تعادل را پیدا کنید .
3
از لیست Eq ، User defined را انتخاب کنید . همچنین برای واکنش کاهش اکسیژن، پارامتر پتانسیل تعادل در این مدل استفاده نمی شود. بنابراین مقدار پیش فرض 0 ولت را حفظ کنید.
4
بخش سینتیک الکترود  را پیدا کنید . از لیست loc,expr ، User defined را انتخاب کنید . در قسمت متن مرتبط، i_c را تایپ کنید .
5
قسمت Active  Specific  Surface  Area را پیدا کنید . در قسمت متن v ، Av را تایپ کنید .
در نهایت، شرایط مرزی خارجی و مقادیر اولیه را تنظیم کنید.
فاز 1 هدایت الکترونیکی
مقدار اولیه پتانسیل الکتریکی در الکترود کاتد را روی پتانسیل سلول تنظیم کنید. (مقادیر اولیه صفر پیش فرض برای هر دو پتانسیل در بقیه هندسه استفاده می شود.)
1
در پنجره Model  Builder ، در قسمت Component   (comp1)>Hydrogen  Fuel  Cell  (fc) روی Electronic  Conducting  Phase  1 کلیک کنید .
مقادیر اولیه، دامنه های O2 1
1
در نوار ابزار Physics ، روی  Attributes کلیک کنید و مقادیر اولیه ،  O2 Domains را انتخاب کنید .
2
فقط دامنه 3 را انتخاب کنید.
3
در پنجره تنظیمات برای مقادیر اولیه ،  دامنه های O2 ، قسمت مقادیر اولیه را پیدا کنید .
4
در قسمت متن φ s ، V_cell را تایپ کنید .
فاز 1 هدایت الکترونیکی
در پنجره Model  Builder ، روی Electronic  Conducting  Phase  1 کلیک کنید .
زمین برق 1
1
در نوار ابزار Physics ، روی  Attributes کلیک کنید و Electric  Ground را انتخاب کنید .
2
فقط مرزهای 8، 23 و 24 را انتخاب کنید.
فاز 1 هدایت الکترونیکی
در پنجره Model  Builder ، روی Electronic  Conducting  Phase  1 کلیک کنید .
پتانسیل الکتریکی 1
1
در نوار ابزار Physics ، روی  Attributes کلیک کنید و Electric  Potential را انتخاب کنید .
2
فقط مرز 17 را انتخاب کنید.
3
در پنجره تنظیمات Electric  Potential ، قسمت Electric  Potential را پیدا کنید .
4
در قسمت متن φ s,bnd ، V_cell را تایپ کنید .
فاز 1 H2 گازی
در پنجره Model  Builder ، در قسمت Component   (comp1)>Hydrogen  Fuel  Cell  (fc) روی H2  Gas  Phase  1 کلیک کنید .
ورودی H2 1
1
در نوار ابزار Physics ، روی  Attributes کلیک کنید و H2  Inlet را انتخاب کنید .
2
فقط مرز 1 را انتخاب کنید.
3
در پنجره تنظیمات برای ورودی H2  ، بخش شرایط مرزی جریان را پیدا کنید .
4
از لیست، فشار را انتخاب کنید .
5
در قسمت متن 0 ، p_a_in را تایپ کنید .
فاز 1 H2 گازی
در پنجره Model  Builder ، روی H2  Gas  Phase  1 کلیک کنید .
ورودی H2 2
1
در نوار ابزار Physics ، روی  Attributes کلیک کنید و H2  Inlet را انتخاب کنید .
2
فقط مرز 4 را انتخاب کنید.
مقادیر اولیه 1
1
در پنجره Model  Builder ، روی مقادیر اولیه  1 کلیک کنید .
2
در پنجره تنظیمات مقادیر اولیه  ، قسمت ترکیب اولیه را پیدا کنید .
3
در قسمت متن 0,H2O ، xH2Oa_in را تایپ کنید .
4
قسمت فشار اولیه  را پیدا کنید . در قسمت متن 0 ، p_ref را تایپ کنید .
فاز 1 گاز O2
در پنجره Model  Builder ، در قسمت Component   (comp1)>Hydrogen  Fuel  Cell  (fc) روی O2  Gas  Phase  1 کلیک کنید .
ورودی O2 1
1
در نوار ابزار Physics ، روی  Attributes کلیک کنید و O2  Inlet را انتخاب کنید .
2
فقط مرز 22 را انتخاب کنید.
3
در پنجره تنظیمات برای ورودی O2  ، بخش شرایط مرزی جریان را پیدا کنید .
4
از لیست، فشار را انتخاب کنید .
5
در قسمت متن 0 ، p_c_in را تایپ کنید .
فاز 1 گاز O2
در پنجره Model  Builder ، روی O2  Gas  Phase  1 کلیک کنید .
ورودی O2 2
1
در نوار ابزار Physics ، روی  Attributes کلیک کنید و O2  Inlet را انتخاب کنید .
2
فقط مرز 21 را انتخاب کنید.
مقادیر اولیه 1
1
در پنجره Model  Builder ، روی مقادیر اولیه  1 کلیک کنید .
2
در پنجره تنظیمات مقادیر اولیه  ، قسمت ترکیب اولیه را پیدا کنید .
3
در قسمت متن 0,H2O ، xH2Oc_in را تایپ کنید .
4
در قسمت متنی 0,N2 ، xN2_in را تایپ کنید .
5
قسمت فشار اولیه  را پیدا کنید . در قسمت متن 0 ، p_ref را تایپ کنید .
تعاریف جهانی
ورودی های مدل پیش فرض
از آنجایی که در همه جای مدل از یک دما استفاده خواهیم کرد، دما را فقط یک بار در گره ورودی های مدل پیش فرض تعریف می کنیم . این گره ممکن است توسط چندین گره فیزیک قابل دسترسی باشد.
1
در پنجره Model  Builder ، در قسمت Global  Definitions، روی Default  Model  Inputs کلیک کنید .
2
در پنجره تنظیمات برای ورودی های مدل پیش فرض  ، بخش Browse Model Inputs را پیدا کنید .
3
در درخت، General>Temperature  (K)  –  minput.T را انتخاب کنید .
4
زیربخش عبارت  برای  انتخاب باقیمانده  را پیدا کنید . در قسمت متن دما ، T را تایپ کنید .
مش 1
شیب تند نزدیک به الکترودها و در نقاط گوشه ای که دنده های صفحات جریان در GDL سرکوب می شوند، انتظار می رود. در این قسمت های هندسه یک مش با وضوح ریزتر ایجاد کنید.
سایز 1
1
در پنجره Model  Builder ، در قسمت Component   (comp1) روی Mesh  1 کلیک راست کرده و Size را انتخاب کنید .
2
در پنجره تنظیمات برای اندازه ، قسمت انتخاب موجودیت هندسی  را پیدا کنید .
3
از لیست سطح نهاد هندسی  ، Boundary را انتخاب کنید .
4
فقط مرزهای 10 و 13 را انتخاب کنید.
5
بخش اندازه عنصر  را پیدا کنید . روی دکمه Custom کلیک کنید .
6
قسمت پارامترهای اندازه عنصر  را پیدا کنید .
7
کادر انتخاب حداکثر  اندازه عنصر را  انتخاب کنید . در قسمت متن مرتبط، 2.5e-5 را تایپ کنید .
سایز ۲
1
در پنجره Model  Builder ، روی Mesh  1 کلیک راست کرده و Size را انتخاب کنید .
2
در پنجره تنظیمات برای اندازه ، قسمت انتخاب موجودیت هندسی  را پیدا کنید .
3
از لیست سطح نهاد هندسی  ، نقطه را انتخاب کنید .
4
فقط نقاط 5، 6، 16 و 17 را انتخاب کنید.
5
بخش اندازه عنصر  را پیدا کنید . روی دکمه Custom کلیک کنید .
6
قسمت پارامترهای اندازه عنصر  را پیدا کنید .
7
کادر انتخاب حداکثر  اندازه عنصر را  انتخاب کنید . در قسمت متن مرتبط، 1e-5 را تایپ کنید .
مثلثی رایگان 1
1
در نوار ابزار Mesh ، روی  Free  Triangular کلیک کنید .
2
روی Mesh  کلیک راست کرده و Build  All را انتخاب کنید .
مش نهایی اکنون باید به صورت زیر باشد:
مطالعه 1
مشکل اکنون برای حل آماده است. در مرحله اول، مقدار اولیه توزیع جریان ثانویه را حل کنید.
مرحله 1: راه اندازی توزیع فعلی
1
در پنجره Model  Builder ، در بخش مطالعه  1 ، روی مرحله  1:  راه‌اندازی توزیع فعلی  کلیک کنید .
2
در پنجره تنظیمات برای راه‌اندازی توزیع فعلی  ، بخش تنظیمات مطالعه را پیدا کنید .
3
از لیست نوع توزیع فعلی ،  ثانویه را انتخاب کنید .
4
در نوار ابزار صفحه اصلی ،  روی محاسبه کلیک کنید .
نتایج
مراحل زیر ارقام موجود در بخش نتایج و بحث را بازتولید می کند .
توزیع چگالی جریان GDL
1
در نوار ابزار صفحه اصلی ، روی  Add  Plot  Group کلیک کنید و 2D  Plot  Group را انتخاب کنید .
2
در پنجره تنظیمات برای گروه طرح دوبعدی  ، GDL Current Density Distribution را در قسمت نوشتاری Label تایپ کنید .
سطح 1
1
روی GDL  Current  Density  Distribution کلیک راست کرده و Surface را انتخاب کنید .
2
در پنجره تنظیمات برای Surface ، روی Replace  Expression در گوشه سمت راست بالای بخش Expression کلیک کنید . از منو، Component   (comp1)>Hydrogen  Fuel  Cell>fc.IsMag  –  قدر چگالی جریان الکترود  – A/m² را انتخاب کنید .
3
برای گسترش بخش Range کلیک کنید . تیک گزینه Manual  color  range را انتخاب کنید .
4
در قسمت Maximum text، 1e4 را تایپ کنید .
سطح پیکان 1
1
در پنجره Model  Builder ، روی GDL  Current  Density  Distribution کلیک راست کرده و Arrow  Surface را انتخاب کنید .
2
در پنجره تنظیمات برای Arrow  Surface ، روی Replace  Expression در گوشه سمت راست بالای بخش Expression کلیک کنید . از منو، Component   (comp1)>Hydrogen  Fuel  Cell>fc.Isx,fc.Isy  -بردار  چگالی جریان الکترود  را انتخاب کنید .
3
در نوار ابزار GDL Current Density Distribution ، روی  Plot کلیک کنید .
4
 روی دکمه Zoom  Extents در نوار ابزار Graphics کلیک کنید .
چگالی جریان واکنش آند
1
در نوار ابزار صفحه اصلی ، روی  Add  Plot  Group کلیک کنید و 1D  Plot  Group را انتخاب کنید .
2
در پنجره تنظیمات برای گروه طرح 1 بعدی  ، چگالی جریان واکنش آند را در قسمت نوشتار برچسب تایپ کنید .
نمودار خطی 1
1
روی Anode  Reaction  Current  Density کلیک راست کرده و Line  Graph را انتخاب کنید .
2
در پنجره تنظیمات برای نمودار خط  ، بخش Selection را پیدا کنید .
3
از لیست انتخاب ، Anode  GDE را انتخاب کنید .
4
روی Replace  Expression در گوشه سمت راست بالای بخش y-Axis  Data کلیک کنید . از منو، Component   (comp1)>Definitions>Variables>i_a  –  چگالی جریان آند  – A/m² را انتخاب کنید .
5
قسمت y-Axis  Data را پیدا کنید . در قسمت متن Expression ، i_a*Av*l_act را تایپ کنید .
6
چک باکس Description را انتخاب کنید . در قسمت متن مرتبط، چگالی جریان واکنش آند را تایپ کنید .
7
روی Replace  Expression در گوشه سمت راست بالای بخش x-Axis  Data کلیک کنید . از منو، Component   (comp1)>Geometry>Coordinate>y  –  y-coordinate را انتخاب کنید .
8
در نوار ابزار Anode Reaction Current Density ، روی  Plot کلیک کنید .
9
 روی دکمه Zoom  Extents در نوار ابزار Graphics کلیک کنید .
میدان سرعت
1
در نوار ابزار صفحه اصلی ، روی  Add  Plot  Group کلیک کنید و 2D  Plot  Group را انتخاب کنید .
2
در پنجره تنظیمات برای گروه طرح دوبعدی  ، فیلد Velocity را در قسمت نوشتار برچسب تایپ کنید .
سطح 1
1
روی Velocity  Field کلیک راست کرده و Surface را انتخاب کنید .
2
در پنجره تنظیمات برای Surface ، روی Replace  Expression در گوشه سمت راست بالای بخش Expression کلیک کنید . از منو، Component   (comp1)>Hydrogen  Fuel  Cell>fc.U  –  Velocity  magnitude  –  m/s را انتخاب کنید .
سطح پیکان 1
1
در پنجره Model  Builder ، روی Velocity  Field کلیک راست کرده و Arrow  Surface را انتخاب کنید .
2
در پنجره تنظیمات برای Arrow  Surface ، روی Replace  Expression در گوشه سمت راست بالای بخش Expression کلیک کنید . از منو، Component   (comp1)>Hydrogen  Fuel  Cell>fc.u,fc.v  –  Velocity را  انتخاب کنید .
3
در نوار ابزار Velocity Field ، روی  Plot کلیک کنید .
4
 روی دکمه Zoom  Extents در نوار ابزار Graphics کلیک کنید .
کسر جرمی واکنش دهنده
1
در نوار ابزار صفحه اصلی ، روی  Add  Plot  Group کلیک کنید و 2D  Plot  Group را انتخاب کنید .
2
در پنجره تنظیمات برای گروه طرح دو بعدی  ، Reactant Mass Fraction را در قسمت نوشتاری Label تایپ کنید .
سطح 1
1
روی Reactant  Mass  Fraction کلیک راست کرده و Surface را انتخاب کنید .
2
در پنجره تنظیمات برای Surface ، بخش Expression را پیدا کنید .
3
در قسمت متن Expression ، wReact را تایپ کنید .
کانتور 1
1
در پنجره Model  Builder ، روی Reactant  Mass  Fraction کلیک راست کرده و Contour را انتخاب کنید .
2
در پنجره تنظیمات برای Contour ، بخش Expression را پیدا کنید .
3
در قسمت متن Expression ، wReact را تایپ کنید .
4
قسمت Coloring  and  Style را پیدا کنید .  روی تغییر  جدول رنگ  کلیک کنید .
5
در کادر محاوره ای Color  Table ، Thermal>ThermalDark را در درخت انتخاب کنید.
6
روی OK کلیک کنید .
7
در پنجره تنظیمات برای Contour ، بخش Coloring  and  Style را پیدا کنید .
8
تیک Color  legend را پاک کنید .
سطح پیکان 1
1
روی Reactant  Mass  Fraction کلیک راست کرده و Arrow  Surface را انتخاب کنید .
2
در پنجره تنظیمات برای Arrow  Surface ، روی Replace  Expression در گوشه سمت راست بالای بخش Expression کلیک کنید . از منو، Component   (comp1)>Hydrogen  Fuel  Cell>fc.u,fc.v  –  Velocity را  انتخاب کنید .
3
در نوار ابزار Reactant Mass Fraction ، روی  Plot کلیک کنید .
4
 روی دکمه Zoom  Extents در نوار ابزار Graphics کلیک کنید .
کسر جرمی آب
1
در نوار ابزار صفحه اصلی ، روی  Add  Plot  Group کلیک کنید و 2D  Plot  Group را انتخاب کنید .
2
در پنجره تنظیمات برای گروه طرح دو بعدی  ، کسر جرمی آب را در قسمت نوشتار برچسب تایپ کنید .
سطح 1
1
روی Water  Mass  Fraction کلیک راست کرده و Surface را انتخاب کنید .
2
در پنجره تنظیمات برای Surface ، بخش Expression را پیدا کنید .
3
در قسمت متن Expression ، fc.wH2O_H2 را تایپ کنید .
کانتور 1
1
در پنجره Model  Builder ، روی Water  Mass  Fraction کلیک راست کرده و Contour را انتخاب کنید .
2
در پنجره تنظیمات برای Contour ، بخش Expression را پیدا کنید .
3
در قسمت متن Expression ، fc.wH2O_O2 را تایپ کنید .
4
در نوار ابزار Water Mass Fraction ، روی  Plot کلیک کنید .
5
 روی دکمه Zoom  Extents در نوار ابزار Graphics کلیک کنید .