I have a 2D Model of a Brain with a scalp layer, stroke layer, brain layer and skull.
Each region has its own conductivity. I have specified conductivity, for Laplaces equation, as:
model = createpde();
specifyCoefficients(model,'m',0,'d',0,'c',1.1,'a',0,'f',0,'Face',2);
specifyCoefficients(model,'m',0,'d',0,'c',-0.03,'a',0,'f',0,'Face',1);
specifyCoefficients(model,'m',0,'d',0,'c',-1,'a',0,'f',0,'Face',3);
%Stroke
specifyCoefficients(model,'m',0,'d',0,'c',-3,'a',0,'f',0,'Face',4);
generateMesh(model,'Hgrad',1,'Hmax',0.05)
----------------------------------------------------------------------------------------------------------------
I solve this with solvepde()
When I apply a current pattern to this EIT setup (applyBoundaryConditions), I get voltages in the range of 1e12, which is very inaccurate. I also get badly scaled vector warnings. I am not sure what is happening. Do I need to individually assign the conductivity to each node of the electrode? Is my current pattern incorrect? Please feel free to ask questions if this is unclear as I was unsure of how much to elaborate.