Beeler-Reuter Action Potential Propagation in a Ring
Description
- In this example, we simulate action potential propagation in a ring. In a small section of the ring, the sodium channels are initialized to be refractory so that propagation is unidirectional.
Start Continuity
- Launch the Continuity 6.4 Client
Go to File→Library→Search, select 2D_ring (ID: 1226), right-click, and select Load.
- A dialog box will pop-up asking if you want to save your current problem (your current problem should be blank).
Click to select Save (without reset) and proceed
Click OK
The shortcut: Execute pre-written script to solve EP problem
Go to File→Scripts→ScriptManager
Select EPscript
Click Execute Selected
- Several minutes should elapse while the Electrophysiology problem executes. You can track the progress of the simulation by opening the Python Shell.
- The solution should render and play automatically.
Execute steps manually
Click lines radio button
Click Render to display mesh lines
Click surfaces radio button
Click Render to display mesh surface
Click Electrophysiology→Solve→Initialize
- This specifies the initial conditions for the ordinary differential equations at each point in the mesh.
Do a Send with File→Send
Click Electrophysiology→Solve→Integration
Set Duration to 40.0 (ms)
Set Step Size to 0.02 (ms)
Under Output tab, select Display solution every 20 steps
click OK to start solving
- You may track progress of the simulation in the Python shell
When solution is complete, verify OpenMesh is the renderer
Press OK when the warning pops up.
If DejaVu is selected, click OpenMesh radio button
Click OK to change renderer
You will have to re-render all objects if you were originally using DejaVu
You may adjust the resolution of the rendered solution is by changing the number at View→Set Divisions...
The Number of Divisions: should be 6 by default, but you may change it to 10.
- Go to Electrophysiology→Render→Render Solution
Set Min value to -80
Set Max value: to 0 or 10
Click OK to view an animated color map of voltage