Stimulus not propagating

Question

Hi,

I am trying to stimulate an external mesh and get activation map. My assumption was that consider the rheobase curve, given a stimulus of strong enough strength and sufficient duration, electrical activity will propagate to neighboring cardiac tissue. I modified the code from 02_EP_tissue/00_simple and 02_EP_tissue/01_basic_usage by changing the stimulus and mesh parameters.

However, I was only able to stimulate the example block mesh. I cannot stimulate an external mesh under same condition, and the activation remains localized like in the picture, not propagating outwards. I varied the stimulus strength, location, duration and geometry, but I am getting strange result combinations. Can anyone help me to suggest how to check the numerics? Or are there physiology parts I should look deeper into? 

Thank you!

Dongcheng


 

Answer 1

Your mesh is very coarse. Average edge length: 1.1mm. (meshtool query edges -msh ellipsoid)

You should aim for at max 400µm (see this openCARP example for background).

Comments

I have resampled the mesh and the propagation problem has been solved, thank you so much!

Answer 2

If you provide a minimal executable example that reproduces your behavior, we can look into it.

Comments

Thanks in advance for helping. Here are the modified code and the ellipsoid mesh.

cangdongcheng/cardiac_stimulation: Minimal working example