I am trying to simulate a system in which two electrodes are separated by oxides and an electrolyte. I would like to plot potential profile when I apply voltage and ground to the electrodes. *A picture of the system is attached below.* I know the potential drop in electrolyte will be minimal...therefore, my goal is really to figure out what proportion of the voltage is dropped in each of the two oxides.
I am going with 2D to keep it simpler. I was thinking I would have to use the electrostatics interface for potential drop in the oxides and the primary current distribution interface (electrochemistry module) for potential drop in electrolyte. I know I cannot use electrostatics for potential drop in electrolyte since that potential drop is dependent on the conductivity of electrolyte NOT its permittivity. I think the primary current distribution physics only account for ohmic drop in the electrolyte and NOT the electrode-electrolyte reaction kinetics...so it is the simplest case. I made a model and added the two sets of physics and set up the boundary conditions, but I am having trouble getting it to solve (some issue with a phis entry missing in the LU matrix). So, I am wondering if I am using the right approach.
Could someone please help? I can upload the model if you think my approach is right and the problem is in my model. I really appreciate your time! Thanks!
EDIT: Of course there will also be a potential drop in the double layers between the electrolyte and oxides on both sides. Does "secondary current distribution" model this effect of double layers?
I am going with 2D to keep it simpler. I was thinking I would have to use the electrostatics interface for potential drop in the oxides and the primary current distribution interface (electrochemistry module) for potential drop in electrolyte. I know I cannot use electrostatics for potential drop in electrolyte since that potential drop is dependent on the conductivity of electrolyte NOT its permittivity. I think the primary current distribution physics only account for ohmic drop in the electrolyte and NOT the electrode-electrolyte reaction kinetics...so it is the simplest case. I made a model and added the two sets of physics and set up the boundary conditions, but I am having trouble getting it to solve (some issue with a phis entry missing in the LU matrix). So, I am wondering if I am using the right approach.
Could someone please help? I can upload the model if you think my approach is right and the problem is in my model. I really appreciate your time! Thanks!
EDIT: Of course there will also be a potential drop in the double layers between the electrolyte and oxides on both sides. Does "secondary current distribution" model this effect of double layers?