CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > Siemens > STAR-CCM+

Modeling a Fuel Cell

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 22, 2018, 13:07
Default Modeling a Fuel Cell
  #1
New Member
 
Andrea Vaini
Join Date: Mar 2018
Posts: 5
Rep Power: 8
Jonsy is on a distinguished road
Hi!
I always used Star CCM+ for internal combustion engine applications but now, for my post-graduate degree thesis, I switched to fuel cell simulations.
I read many studies about PEMFC modeling but it seems that nobody actually use Star CCM for this kind of applications. So i 'studied' user guide chapters about electrochemistry and surface reaction models but, in one month, i still can not run a decent simulation.

I would like to simulate a discharge with a Cell Voltage of 0.77 V.
The only simulation involving electrochemical surface reactions that seemed to reach convergence is the following:

1. Geometry:
The simplest possible cell, formed by 3 closed parallelepipedons representing: anode (multi-component flowing gas region: H2 + H2O), cathode (multi-component flowing gas region: O2 + H2O) and MEA (Membrane Electrode Assembly solid region)
Mesh.PNG

2. Physics
- Anode and Cathode:
3D, steady, multi-component gas, non reacting, segregated flow and species, gradients, costant density, laminar flow, electrochemistry, electrochemical reactions, electromagnetic, electrodynamic potential.
- MEA:
3D, steady, solid, gradients, electrochemistry, electrochemical reactions, electromagnetics, electrodynamic potential, solid ion.
I'd like to model a Nafion MEA, so i thought to costumize the standard solid material modifying its electrical conductivity according to the proton conductivity of Nafion.

3. Boundary Conditions
- Anode and Cathode:
mass flow inlet, pressure outlet

- Interface Anode/MEA
oxidation surface reaction: 2 H2 -> 4e- + 4H+
Equilibrium Potential: Field function
Specific Exchange Current Density: 1500 A/m^2

- Interface Cathode/MEA
reduction surface reaction: O2 + 4H+ + 4e- -> 2 H2O
Equilibrium Potential: Field function
Specific Exchange Current Density: 1.5 A/m^2

- External Anode Boundary (wall)
Electric Potential: 0 V

- External Cathode Boundary (wall)
Electric Potential: 0.77 V

4. Results
Scalar Velocity.PNG
Scalar Electric Potential.PNG
Scalar Electric Current Density.PNG
Residuals.PNG

I know that this model is far too simple and result values don't mean anything, but i'm trying to understand how it works and what i could improve.

I'll be thankful to anyone could give me some tips and advices!

Last edited by Jonsy; March 22, 2018 at 17:45.
Jonsy is offline   Reply With Quote

Old   April 9, 2018, 08:41
Default
  #2
New Member
 
Join Date: Apr 2018
Posts: 1
Rep Power: 0
floaty is on a distinguished road
Hello Jonsy,

the 13.02 tutorial guide has an SOFC tutorial which should represent a good start.

Regards, Christian
floaty is offline   Reply With Quote

Old   April 9, 2018, 08:49
Default
  #3
New Member
 
Andrea Vaini
Join Date: Mar 2018
Posts: 5
Rep Power: 8
Jonsy is on a distinguished road
Quote:
Originally Posted by floaty View Post
Hello Jonsy,

the 13.02 tutorial guide has an SOFC tutorial which should represent a good start.

Regards, Christian

actually i used 12.02 user guide (according to my software version) but unfortunately it only talks about SOFC but with no details neither tutorials. I’ll check it out the new guide version for sure, so thank you Christian! ����
Jonsy is offline   Reply With Quote

Old   April 9, 2018, 09:46
Default
  #4
New Member
 
Andrea Vaini
Join Date: Mar 2018
Posts: 5
Rep Power: 8
Jonsy is on a distinguished road
Excuse me can you tell me where to find and to download 13.02 user guide?
i checked on internet but finding that is more difficult than i thought...i don’t have credentials to access to steveportal so can anyone share pdf of it (or just the part about sofc tutorial)?
Jonsy is offline   Reply With Quote

Old   April 9, 2018, 13:39
Default
  #5
New Member
 
Andrea Vaini
Join Date: Mar 2018
Posts: 5
Rep Power: 8
Jonsy is on a distinguished road
Ok I found it! Thanks again
Jonsy is offline   Reply With Quote

Old   September 20, 2018, 10:39
Default
  #6
New Member
 
Dominik
Join Date: Aug 2017
Posts: 2
Rep Power: 0
cfdominik is on a distinguished road
Hi everyone!

I am also trying to set up a PEMFC model with star ccm+ 13.04. to analyze the generation of liquid water at the cathode.
In the "whats new in version 13.02." it says you can use a Mixture Multiphase (MMP) Model to simulate liquid water in the fuel cell.

To do so, I chose the following physics models for the air (cathode side) continuum:
- Electrodynamic Potential
- Electromagnetism
- Electrochemical Reactions
- Electrochemistry
- Laminar
- Gradients
- Multiphase Equation of State
- Multiphase Segregated Flow
- Multiphase Interaction
- Eulerian Multiphase
- Steady
- Three Dimensional

I then started setting up my two phases for the Eulerian Multiphase Model:
Gas Phase (containing nitrogen (N2) and oxygen (O2))
- Multi-Component Gas
- Flow
- Laminar
- Non-reacting
- Segregated Species
- Constant Density
- Segregated Fluid Isothermal


Liquid Phase (Water)
- Liquid
- Flow
- Laminar
- Constant Density
- Segregated Fluid Isothermal

I then want to set up the surface reactions at the proton exchange membrane (PEM) but I can not select any of the defined gases or liquids as reactants or products, more precisely there is no "Reactants" or "Products" folder in the simulation tree, only the "Properties".

However, when you create a physics continuum with just the Multi-Component Gas Model as Material it's no problem to select the defined components as reactants or products.

Could you successfully develop your model, Jonsy?
I am very thankful for any kind of advice!


dominik
cfdominik is offline   Reply With Quote

Old   July 28, 2024, 13:44
Default contour lines
  #7
DRC
New Member
 
Rhei
Join Date: Mar 2024
Posts: 2
Rep Power: 0
DRC is on a distinguished road
Hi! Sorry for disturbing everyone here despite having this forum a few years back. I hope anyone from here can explain to me what could have caused why the lines for the contour generated are not connecting through the plane. I was thinking that it should be connected both for in-plane and through-plane. I am generating results that are much lower than the experimental data so the interface problem might have caused it.

Also, how did you generate your polarization curves? As I consider that I might be doing it incorrectly?

Your guidance will be highly appreciated.
Attached Images
File Type: jpg Picture1.jpg (96.7 KB, 5 views)
File Type: png polarization.png (72.6 KB, 5 views)
DRC is offline   Reply With Quote

Reply

Tags
fuel cell, fuel cell model, pem, pem fuel cell, pemfc modeling

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[Other] refineWallLayer Error Yuby OpenFOAM Meshing & Mesh Conversion 2 November 11, 2021 12:04
problem with PEM fuel cell modeling and I-V curve AmiN.D FLUENT 3 July 24, 2017 02:39
Fuel cell in star ccm Xuekun Siemens 2 January 27, 2016 20:01
FvMatrix coefficients shrina OpenFOAM Running, Solving & CFD 10 October 3, 2013 15:38
Warning 097- AB Siemens 6 November 15, 2004 05:41


All times are GMT -4. The time now is 14:15.