Hello everyone,
I would like to ask you for help. I would like to validate a model that I am working on by calculating the current that is flowing through a simple model. For this purpose I use the Mickey Mouse TDCS model supplied in SCIRun 4.7. I would like to confirm that the flowing current that I measure equals the values given in the settings of the model (the RHS of the linear set of equations). My procedure:
1. Calculate the current density within the model (nothing out of the ordinary here: -gradPhi*sigma)
2. Map the normal component of the current density onto a plane that is a cross-section of the entire model. To simplify the problem I just picked a vertical plane parallel to both electrodes.
3. Multiply the normal component of the current density from the plane by the areas of the elements of the plane to get the flowing current per plane element
4. Sum all such "unit" currents to get the total current flowing between the two electrodes.
As you can imagine, the values that I get by doing this are off from the values set to the input of the model and it is not just a matter of scaling. When I multiply the solution matrix by the TDCS matrix, the RHS that I get equals the one that I set into the model, so here everything is fine. I tried to study the maths behind the TDCS model from:
[1] P. Hua, Eung Je Woo, J. G. Webster, and W. J. Tompkins, “Finite element modeling of electrode-skin contact impedance in electrical impedance tomography,” IEEE Trans. Biomed. Eng., vol. 40, no. 4, pp. 335–343, 1993.
Current shunting along the electrodes was mentioned in the paper. Could it be that some current set into the electrodes simply does not get through into the model? Or am I making some sort of a mistake that has been oblivious to me so far?