[jose.ceballos]

I'm working with pemfcMultiphaseModel-4.0 solver, and I'm trying to build my IV Curve. I downloaded the zip file from https://data.mendeley.com/datasets/c743sh73j8/1 and extracted it. I want to build the IV curve without changes. Just indicating the Ibar0 and V in run/constant/cellProperties. But I have some problems with high values of density current. I executed the solver (just like the way it is, no changes to the model) and I obtain this values:

V	I	My run				
1	5	5.5				
0.95	25	27				
0.9	100	110				
0.85	500	510				
0.8	1000	1000				
0.78	2000	2100				
0.74	3000	3200				
0.7	4000	4400				
0.65	5000	5700				
0.6	6000	6900				
0.58	7000	8300				
0.54	8000	9800				
0.5	9000	11000				
0.48	10000	11300				
0.44	11000	No Convergence (Solving Electrochemestry min, mean, max (Sigma))				
0.4	12000	No Convergence (Solving Electrochemestry min, mean, max (Sigma))				
0.38	13000	No Convergence (Solving Electrochemestry min, mean, max (Sigma))				
0.34	14000	No Convergence (Solving Electrochemestry min, mean, max (Sigma))				
0.3	15000	No Convergence (Solving Electrochemestry min, mean, max (Sigma))				
0.25	16000	No Convergence (Solving Electrochemestry min, mean, max (Sigma))

I'm running the solver in Ubuntu 18.04LTS and OpenFOAM 4.x. The initials conditions are the sames as the zip file. The idea is to change the geometry and data to validate others polarization curves, but if I can't build this model that should be already validated I wont be able to use a different geometry. I already have a different geometry but I have the same problems and calculations stop very quickly.

If anyone have experience whit this code built by Kone I would appreciate any help.

[Morteza123]

Originally Posted by jose.ceballos

I'm working with pemfcMultiphaseModel-4.0 solver, and I'm trying to build my IV Curve. I downloaded the zip file from https://data.mendeley.com/datasets/c743sh73j8/1 and extracted it. I want to build the IV curve without changes. Just indicating the Ibar0 and V in run/constant/cellProperties. But I have some problems with high values of density current. I executed the solver (just like the way it is, no changes to the model) and I obtain this values:

V	I	My run				
1	5	5.5				
0.95	25	27				
0.9	100	110				
0.85	500	510				
0.8	1000	1000				
0.78	2000	2100				
0.74	3000	3200				
0.7	4000	4400				
0.65	5000	5700				
0.6	6000	6900				
0.58	7000	8300				
0.54	8000	9800				
0.5	9000	11000				
0.48	10000	11300				
0.44	11000	No Convergence (Solving Electrochemestry min, mean, max (Sigma))				
0.4	12000	No Convergence (Solving Electrochemestry min, mean, max (Sigma))				
0.38	13000	No Convergence (Solving Electrochemestry min, mean, max (Sigma))				
0.34	14000	No Convergence (Solving Electrochemestry min, mean, max (Sigma))				
0.3	15000	No Convergence (Solving Electrochemestry min, mean, max (Sigma))				
0.25	16000	No Convergence (Solving Electrochemestry min, mean, max (Sigma))

I'm running the solver in Ubuntu 18.04LTS and OpenFOAM 4.x. The initials conditions are the sames as the zip file. The idea is to change the geometry and data to validate others polarization curves, but if I can't build this model that should be already validated I wont be able to use a different geometry. I already have a different geometry but I have the same problems and calculations stop very quickly.

If anyone have experience whit this code built by Kone I would appreciate any help.

Using this tool you should be very careful. I think it's full of bugs. As an example, at some point in order to compute Capillary pressure, the code passes angles in degree as an argument to a cosine function. And many other formulas implemented in the code do not match with the paper. I don't know how it's possible to have this level of bugs in an open software.

[ramakant]

Hi...

I am also working on the same code. We can share our problems?

Thank you