Capacity mismatch

[venkichav]

I am active user of COMSOL. recently i started exploring PyBamm. When i tried to build the same model in PyBamm even though i am defining my cell nominal capacity as 25Ah. I am getting just 0.68Ah as a discharge capacity even at 0.1C. Below are my parameter list
I feel there could be some issue with the parameter values.

parameter_values=pybamm.ParameterValues('Chen2020')
parameter_values['Ambient temperature [K]']= 298.15
parameter_values['Boltzmann constant [J.K-1]']= 1.380649e-23
parameter_values['Bulk solvent concentration [mol.m-3]']= 2636.0 #need to check
parameter_values['Cation transference number']= 0.2594 # need to check
parameter_values['Cell cooling surface area [m2]']= 0.0115
parameter_values['Cell thermal expansion coefficient [m.K-1]']= 1.1e-06
parameter_values['Cell volume [m3]']= 1.142e-04 # need to check
parameter_values['Contact resistance [Ohm]']= 0 # need to check
parameter_values['Current function [A]']= 25.13
parameter_values['EC diffusivity [m2.s-1]']= 2e-18 # need to check
parameter_values['EC initial concentration in electrolyte [mol.m-3]']= 4541.0 # need to check
parameter_values['Electrode height [m]']= 0.069
parameter_values['Electrode width [m]']= 3.908
#parameter_values['Electrolyte conductivity [S.m-1]']= <function electrolyte_conductivity_Nyman2008 at 0x7c78b30e8f40>
#parameter_values['Electrolyte diffusivity [m2.s-1]']= <function electrolyte_diffusivity_Nyman2008 at 0x7c78b30e8ea0>
parameter_values['Electron charge [C]']= 1.602176634e-19
parameter_values['Faraday constant [C.mol-1]']= 96485.33212
parameter_values['Ideal gas constant [J.K-1.mol-1]']= 8.314462618
parameter_values['Initial SEI thickness [m]']= 5e-09
parameter_values['Initial concentration in electrolyte [mol.m-3]']= 1200.0
parameter_values['Initial concentration in negative electrode [mol.m-3]']=29866.0 # need to check
parameter_values['Initial concentration in positive electrode [mol.m-3]']=17038.0 # need to check
parameter_values['Initial temperature [K]']= 298.15
parameter_values['Lithium interstitial reference concentration [mol.m-3]']= 15.0 # need to check
parameter_values['Lower voltage cut-off [V]']= 2.5
parameter_values['Maximum concentration in negative electrode [mol.m-3]']= 33133.0 # need to check
parameter_values['Maximum concentration in positive electrode [mol.m-3]']= 63104.0 # need to check
parameter_values['Negative current collector conductivity [S.m-1]']= 58411000.0
parameter_values['Negative current collector density [kg.m-3]']= 8960.0
parameter_values['Negative current collector specific heat capacity [J.kg-1.K-1]']= 385.0
parameter_values['Negative current collector thermal conductivity [W.m-1.K-1]']= 401.0
parameter_values['Negative current collector thickness [m]']= 8e-06
parameter_values['Negative electrode Bruggeman coefficient (electrode)']= 0
parameter_values['Negative electrode Bruggeman coefficient (electrolyte)']= 1.5#2.3 #From eqn Tort=(Por)**(1-Brugexpo)
#parameter_values['Negative electrode OCP [V]']= <function graphite_LGM50_ocp_Chen2020 at 0x7c78b30e8c20>
parameter_values['Negative electrode OCP entropic change [V.K-1]']= 0.0
parameter_values['Negative electrode active material volume fraction']= 0.669 #Including electrolyte
parameter_values['Negative electrode charge transfer coefficient']= 0.5
parameter_values['Negative electrode conductivity [S.m-1]']= 53.47
parameter_values['Negative electrode density [kg.m-3]']= 1580
parameter_values['Negative electrode double-layer capacity [F.m-2]']= 0.2 #need to check
#parameter_values['Negative electrode exchange-current density [A.m-2]']= <function graphite_LGM50_electrolyte_exchange_current_density_Chen2020 at 0x7c78b30e8cc0>
parameter_values['Negative electrode porosity']= 0.303
parameter_values['Negative electrode reaction-driven LAM factor [m3.mol-1]']= 0.0 #need to check
parameter_values['Negative electrode specific heat capacity [J.kg-1.K-1]']= 700.0
parameter_values['Negative electrode thermal conductivity [W.m-1.K-1]']= 1.5#1.7
parameter_values['Negative electrode thickness [m]']= 201e-06
parameter_values['Negative particle diffusivity [m2.s-1]']= 3.3e-14 #need to check
parameter_values['Negative particle radius [m]']= 1.35e-05
parameter_values['Nominal cell capacity [A.h]']= 25.13
parameter_values['Number of cells connected in series to make a battery']= 1.0
parameter_values['Number of electrodes connected in parallel to make a cell']= 1.0
parameter_values['Open-circuit voltage at 0% SOC [V]']= 2.5
parameter_values['Open-circuit voltage at 100% SOC [V]']= 4.2
parameter_values['Positive current collector conductivity [S.m-1]']= 36914000.0
parameter_values['Positive current collector density [kg.m-3]']= 2700.0
parameter_values['Positive current collector specific heat capacity [J.kg-1.K-1]']= 897.0
parameter_values['Positive current collector thermal conductivity [W.m-1.K-1]']= 237.0
parameter_values['Positive current collector thickness [m]']= 1.2e-05
parameter_values['Positive electrode Bruggeman coefficient (electrode)']= 0
parameter_values['Positive electrode Bruggeman coefficient (electrolyte)']= 1.70
#parameter_values['Positive electrode OCP [V]']= <function nmc_LGM50_ocp_Chen2020 at 0x7c78b30e8d60>
parameter_values['Positive electrode OCP entropic change [V.K-1]']= 0.0
parameter_values['Positive electrode active material volume fraction']= 0.792
parameter_values['Positive electrode charge transfer coefficient']= 0.5
parameter_values['Positive electrode conductivity [S.m-1]']= 0.18
parameter_values['Positive electrode density [kg.m-3]']= 3580 # need to check its press density or true density
parameter_values['Positive electrode double-layer capacity [F.m-2]']= 0.2
#parameter_values['Positive electrode exchange-current density [A.m-2]']= <function nmc_LGM50_electrolyte_exchange_current_density_Chen2020 at 0x7c78b30e8e00>
parameter_values['Positive electrode porosity']= 0.206
parameter_values['Positive electrode reaction-driven LAM factor [m3.mol-1]']= 0.0
parameter_values['Positive electrode specific heat capacity [J.kg-1.K-1]']= 700.0
parameter_values['Positive electrode thermal conductivity [W.m-1.K-1]']= 2.1
parameter_values['Positive electrode thickness [m]']= 141e-06
parameter_values['Positive particle diffusivity [m2.s-1]']= 4e-15
parameter_values['Positive particle radius [m]']= 9.83e-06
parameter_values['Ratio of lithium moles to SEI moles']= 2.0
parameter_values['Reference temperature [K]']= 298.15
parameter_values['SEI electron conductivity [S.m-1]']= 8.95e-14
parameter_values['SEI growth activation energy [J.mol-1]']= 0.0
parameter_values['SEI kinetic rate constant [m.s-1]']= 1e-12
parameter_values['SEI lithium interstitial diffusivity [m2.s-1]']= 1e-20
parameter_values['SEI open-circuit potential [V]']= 0.4
parameter_values['SEI partial molar volume [m3.mol-1]']= 9.585e-05
parameter_values['SEI reaction exchange current density [A.m-2]']= 1.5e-07
parameter_values['SEI resistivity [Ohm.m]']= 200000.0
parameter_values['SEI solvent diffusivity [m2.s-1]']= 2.5e-22
parameter_values['Separator Bruggeman coefficient (electrolyte)']= 1.5#1.92
parameter_values['Separator density [kg.m-3]']= 960
parameter_values['Separator porosity']= 0.44
parameter_values['Separator specific heat capacity [J.kg-1.K-1]']= 700.0
parameter_values['Separator thermal conductivity [W.m-1.K-1]']= 0.16
parameter_values['Separator thickness [m]']= 1.2e-05
parameter_values['Thermodynamic factor']= 1.0
parameter_values['Total heat transfer coefficient [W.m-2.K-1]']= 10.0
parameter_values['Upper voltage cut-off [V]']= 4.2

Below is my code
model=pybamm.lithium_ion.DFN()
#experiment=pybamm.Experiment(

['Discharge at 1C for 1 hours']

#)
experiment=pybamm.Experiment(
[
'Discharge at 0.1C for 10 hours'
]
)
sim=pybamm.Simulation(model,experiment=experiment)

[github-actions[bot]]

Hi! We have now moved our discussions to Discourse. Please post your question there.