[Dan96]

Dear community,


I want to create a 3D quasi-steady fuel cell model. The model input are two educt fluxes (cathodic and anodic). I want to calculate, which current density distribution and which overpotential distribution result from the given input.



Mass transport is already described by Darcy (convective flux) and Maxwell-Stefan (diffusive flux).


Now I want to add a reaction layer in which the educts convert to products. The reaction is adds a source term in species conservation equation.


What I don't understand is how to calculate the source term (which depends on the current density distribution according to Faraday's law) ...


In literature I found that the kinetics in a fuel cell are mainly described by the butler-volmer-equation. In a 1D-model the current density can be calculated from a given voltage (should be given by the load) or overpotential . In a 3D-Model you need the overpotential distribution. The load should be the mean of the spatial resolved cell voltages ...


My question is, how do you calculate the source term in a 3D spatial resolved fuel cell model?


Best regards


Daniel