Stimulation

Module: devtests.bidomain.stimulation.run

Section author: Gernot Plank <gernot.plank@medunigraz.at>

Example of different types of stimulation either through the application of extracellular voltages (Dirichlet) or through the injection/withdrawal of extracellular currents (Neumann).

Problem Setup

A thin strand of tissue of 1 cm lenght is generated. Electrodes are located at both caps of the strand with an additional auxiliary electrode in the very center of the strand. Electrode locations are shown belows:

images/strand_electrodes.png

Each electrode can be used to inject/withdraw current, prescribe an extracellular potential or ground the extracellular potential.

Experiments

Several experiments are defined:

  • extra_V - Stimulation through application of extracellular voltage. The potential on the left cap is clamped to 2000 mV for a duration of 2 ms.
  • extra_V_bal - Stimulation through application of extracellular voltage. The potential on the left cap is clamped to 1000 mV for a duration of 2 ms and the right hand electrode to -1000mV.
  • extra_V_OL - Stimulation through application of extracellular voltage. The potential on the left cap is clamped to 2000 mV for a duration of xx ms. In contrast to the extra_V case the left hand electrode is allowed to float after the end of the shock.
  • extra_I - Stimulation through application of extracellular current. Extracellular current is injected into the left hand electrode to cause an extracellular potential drop of 2000 mV across the strand.
  • extra_I_bal - Stimulation through application of extracellular current. Extracellular current is injected into the left hand electrode and the same amount is withdrawn from the right hand side electrode. Overall this also causes an extracellular potential drop of 2000 mV across the strand, albeit in a symmetric way.

Optional Arguments

  • grounded - A ground can be used or not. In the latter case, no grounding electrode
    is used, that is, a pure Neumann-problem is solved. A rank-one stabilization is used in this case which essentially enforces the integral over the extracellular potential to be zero at each instant in time, but no specific point is clamped.

Tests

extra_V

Tags: FAST SERIAL

Checks:

  • Compare against stored reference: max_error(phie.igb)
  • Compare against stored reference: max_error(vm.igb)

extra_V_OL

Tags: FAST SERIAL

Checks:

  • Compare against stored reference: max_error(phie.igb)
  • Compare against stored reference: max_error(vm.igb)

extra_I

Tags: FAST SERIAL

Checks:

  • Compare against stored reference: max_error(phie.igb)
  • Compare against stored reference: max_error(vm.igb)

extra_IbalGnd

Tags: FAST SERIAL

Checks:

  • Compare against stored reference: max_error(phie.igb)
  • Compare against stored reference: max_error(vm.igb)

extra_IbalNGnd

Tags: FAST SERIAL

Checks:

  • Compare against stored reference: max_error(phie.igb)
  • Compare against stored reference: max_error(vm.igb)

overstim_protect

Tags: FAST SERIAL

Checks:

  • Compare against stored reference: max_error(vm.igb)