|March 28, 2013, 19:14||
Current density visualisation (PEM fuel cell add-on module)
Join Date: Aug 2011
Posts: 13Rep Power: 5
I am using a PEMFC add-on module to simulate a simple pemfc with a serpentine flow field. I believe I have achieved a converged solution and want to visualise the local current density distribution.
What really confuses me is that, every single published paper refers to different location (surface) when they show the predicted local current density distribution. Where (which surface) exactly should I be looking at?
I have investigated all these surfaces (interfaces);
1. mid-plane of the membrane
2. anode CL/membrane interface
3. cathode CL/membrane interface
4. anode CL/anode GDL interface
5. cathode CL/cathode GDL interface
6. anode GDL/flow field interface (no current collector)
7. cathode GDL/flow field interface (no current collector)
All I can see so far is that the current (I am using User Defined Memory, UDM3: Current Flux Density Magnitude [A/m^2, or current per area] as the variable of interest) at the anode side is a lot less than that at the cathode side, so regardless of which cathode-side interfaces, we should be interested in only the current density at the cathode side right? Please correct me if I am wrong. And then, which interface should I display the contour of local current density distribution exactly? (Interfaces 3 VS 5 VS 7?)
My another question is regarding the variable of interest; the CURRENT DENSITY as referred to as the performance parameter of the PEMFC in the open literature. Which UDS, UDM do we mean by this performance parameter CURRENT DENSITY?
I am using UDM3 as the unit is A/m^2 which matches the unit of current density in the published papers. What about UDM13: Transfer Current [A/m^3, or current per volume]? Can someone point out the difference between these two and where and when shall we look at each of these parameters?
Additionally, I have found a bizarre contour as a result from using UDM3 VS UDM13; I picked the interface between cathode CL/membrane and plotted these contours;
Please note that the oxygen inlet is located at the bottom-left and exits at the top right of the cell.
pic1: contour of UDM3: Current Flux Density Magnitude
pic2: contour of UDM13: Transfer Current
pic3: contour of mass fraction of oxygen
pic4: contour of mass fraction of water VAPOUR
pic5: contour of volume fraction of LIQUID water (UDS2: liquid saturation, s)
Many papers suggest that their current density (here I do not know which variable to use yet) follows the distribution of oxygen and opposes the distribution of liquid water. This means at higher current density will occur near the inlet area where oxygen is present and less current density will occur at the exit (area with high water saturation which is the indication of water flooding).
From the discussion above, using UDM13 gives a sensible result and agrees with the contour of oxygen, water vapour, and liquid water. The only frustrating issue is its unit of current per volume, not current per area.
On the other hand, UDM3 contour seems to be regardless of species distribution. What happens here?
One more question about the liquid/vapour water. From my understanding, the contour of water vapour in pic4 is not the indication of water flooding in the cell, as the name implies, it only shows the region (in red) where the mass fraction of water is very high due to the fact that oxygen gas is consumed in the upstream region. Is this right? Can someone assure me that UDS2: saturation is equivalent to the volume fraction of liquid water, s in multiphase model calculation?
Another question which is more generic. Tutorial given from ANSYS suggest a way to check for a good converged solution by checking mass flux balance through the REPORT > REPORT FLUX tool. When I calculate the net mass flux, the net flux is in the same order of magnitude of the inlet flux (ANSYS suggest net flux should be 0.1% of the inlet flux to ensure good convergence). However, the point monitors throughout the domain show that the variables have converged though. What is happening here? What is wrong here? Is this because the simulation involves species consumption?
When the programs iterates, it returns the value of averaged current density for both anode and cathode (why both? and what is the physical meaning of this averaged current density?), how does it calculate this value? I guess it must be something about surface integral over a certain surface, which interface is this? How exactly do we obtain this averaged current density?
Last question is that if the multiphase option is activated (under PEMFC model sub-tick box, not a separate multiphase model), if my reactant gases are both fully humidified. What is the correct boundary condition at the inlet for water saturation under UDS tab? Should this be "specified value" or "specified flux" and what value?
Please feel free to ask if you need some clarification on my questions. Thank you very much.
Last edited by pchoopanya; March 28, 2013 at 19:34.
|March 29, 2013, 08:43||
Join Date: Jan 2013
Posts: 175Rep Power: 3
I'm working on PEM fuel cell as well. you asked a couple of questions. Regarding exact place for displaying contours, I suppose that most of the papered used membrane-catalyst interface for contour. You also talked about current density, when you start iteration, current density is reported in console window. you can also get the same value for current density from report-surface- integral type and The surface you need to select is wall terminal whether in anode or cathode not important.
|January 20, 2014, 05:14||
Join Date: Jan 2014
Posts: 7Rep Power: 2
To plot the IV curve, can I just use the current density (unit : A/cm^2) that show in every iteration so that no need to divide the value from surface integral function by the area of the cell. and which one that i can use? cathode or anode? because the value is not the same.
and the final question, the massage "reversed flow" at outlet is shown and I have tried many ways to fix this but it still occurs. What does it mean? something flow back from outlet? how to fix it? (this thing has occurred since high voltage (1.0 0.9 V) and I haven't tried any low voltage yet) Is it mean that the solution has not converged yet?
Any comments would be highly appreciated
Thank you very much
Last edited by P.Charoen-amornkitt; January 20, 2014 at 08:48. Reason: wrong spelling
|Thread||Thread Starter||Forum||Replies||Last Post|
|FvMatrix coefficients||shrina||OpenFOAM Running, Solving & CFD||10||October 3, 2013 14:38|
|2D Fuel Cell Modeling||Nehal||FLUENT||2||August 7, 2012 07:47|
|Write UDF to add a sink term for PEM fuel cell design||kitrax||Fluent UDF and Scheme Programming||0||March 15, 2012 05:56|
|how to add OpenFoam in module environment||smpark||OpenFOAM||4||May 12, 2011 14:00|
|electrochemical reacation in the fuel cell?||Kim||FLUENT||5||June 23, 2004 16:44|