Dear All,
I am trying to solve an electrical potential distribution problem studied on a 1D axial symmetric model using the PDE Coefficient Form module.
My model is comprised of five coupled problems: three of them use the electrical potential V as an unknown function on distance, x, while the other two problems have the electrical charge (rho) as unknown function. The problem occurs at the V(x) electrical potential distribution along the 0-x direction. The electrical potential should have a negative value in the intracellular domain, it should rise from this value to zero in the membrane domain, and it should remain zero in the extracellular domain. Instead, my electrical potential gets a negative value in the intracellular domain, rises up to a certain value without reaching zero, and in the extracellular domain it continues to rise up until it reaches zero. The electrical charge (rho) needs to be zero along the o-x direction in the intra- and extra-cellular computational domains (where defined as rho1 and rho3).
The mathematical model and the boundary conditions are explicitly described in the .pdf file attached to ths topic. I also attached an .mph file that contains the Comsol model and three .png files that includes the ideal and my V(x) and rho(x) dependencies.
Does anyone encountered a similar problem? Any advice would be much appreciated.
I am trying to solve an electrical potential distribution problem studied on a 1D axial symmetric model using the PDE Coefficient Form module.
My model is comprised of five coupled problems: three of them use the electrical potential V as an unknown function on distance, x, while the other two problems have the electrical charge (rho) as unknown function. The problem occurs at the V(x) electrical potential distribution along the 0-x direction. The electrical potential should have a negative value in the intracellular domain, it should rise from this value to zero in the membrane domain, and it should remain zero in the extracellular domain. Instead, my electrical potential gets a negative value in the intracellular domain, rises up to a certain value without reaching zero, and in the extracellular domain it continues to rise up until it reaches zero. The electrical charge (rho) needs to be zero along the o-x direction in the intra- and extra-cellular computational domains (where defined as rho1 and rho3).
The mathematical model and the boundary conditions are explicitly described in the .pdf file attached to ths topic. I also attached an .mph file that contains the Comsol model and three .png files that includes the ideal and my V(x) and rho(x) dependencies.
Does anyone encountered a similar problem? Any advice would be much appreciated.