Chaste Release::3.1
Nash2004ContractionModel Class Reference

#include <Nash2004ContractionModel.hpp>

Inheritance diagram for Nash2004ContractionModel:
Collaboration diagram for Nash2004ContractionModel:

List of all members.

Public Member Functions

 Nash2004ContractionModel ()
void EvaluateYDerivatives (double time, const std::vector< double > &rY, std::vector< double > &rDY)
void SetInputParameters (ContractionModelInputParameters &rInputParameters)
void SetStretchAndStretchRate (double stretch, double stretchRate)
double GetActiveTension ()
double GetNextActiveTension ()
bool IsStretchDependent ()
bool IsStretchRateDependent ()

Private Attributes

double mScaledVoltage

Static Private Attributes

static const double kTa = 47.9
static const double e0ByT0 = 1.0/100

Detailed Description

Nash2004 contraction model - Nash & Panfilov "Electromechanical mode of excitable tissue to study reentrant cardiac arrthrymias", Progress in Biophysics and Molecular Biology, 2004.

A simple, stretch- and stretch-rate-independent contraction model, just dependent on the voltage. If v is the non-dimensionalised voltage, Ta is given by dTa/dt = eps(v) (kTa*v - Ta); where eps(v) = e0 if v < 0.05, = 10*e0 v >= 0.05.

We use the non-dimensionalisation: V in [-85,40] ---> v in [0,1], ie v=(V+85)/125

Not sure what the appropriate value of e0 is: as the paper uses non-dimensionalised time the e0 above corresponds to "eps0/t0", where eps0 is value used in the paper (1.0) and t0 is the characteristic time. The paper suggests using t0=25.9 (?), which gives Ta growing too quickly at the beginning (as rapidly as the voltage). t0 = 100 coded at the moment..

Definition at line 60 of file Nash2004ContractionModel.hpp.


Constructor & Destructor Documentation

Nash2004ContractionModel::Nash2004ContractionModel ( ) [inline]

Member Function Documentation

void Nash2004ContractionModel::EvaluateYDerivatives ( double  time,
const std::vector< double > &  rY,
std::vector< double > &  rDY 
) [inline, virtual]

Calculate the derivative of the Ta (the state variable)

Parameters:
timetime
rY1D vector containing Ta
rDY1D vector in which dTa/dt is set

Implements AbstractOdeSystem.

Definition at line 87 of file Nash2004ContractionModel.hpp.

References e0ByT0, kTa, and mScaledVoltage.

double Nash2004ContractionModel::GetActiveTension ( ) [inline, virtual]

Get the current active tension

Implements AbstractContractionModel.

Definition at line 115 of file Nash2004ContractionModel.hpp.

References AbstractParameterisedSystem< std::vector< double > >::rGetStateVariables().

double Nash2004ContractionModel::GetNextActiveTension ( ) [inline, virtual]

Get the active tension corresponding to the temporary stored state variables produced by calling RunDoNotUpdate (and before calling UpdateStateVariables())

Implements AbstractContractionModel.

Definition at line 124 of file Nash2004ContractionModel.hpp.

References AbstractOdeBasedContractionModel::mTemporaryStateVariables.

bool Nash2004ContractionModel::IsStretchDependent ( ) [inline, virtual]

This model is stretch-independent

Implements AbstractContractionModel.

Definition at line 133 of file Nash2004ContractionModel.hpp.

bool Nash2004ContractionModel::IsStretchRateDependent ( ) [inline, virtual]

This model is stretch-rate-independent

Implements AbstractContractionModel.

Definition at line 141 of file Nash2004ContractionModel.hpp.

void Nash2004ContractionModel::SetInputParameters ( ContractionModelInputParameters rInputParameters) [inline, virtual]

Set the input parameters. Only the voltage is used

Parameters:
rInputParametersreference to the input parameters

Implements AbstractContractionModel.

Definition at line 97 of file Nash2004ContractionModel.hpp.

References DOUBLE_UNSET, mScaledVoltage, and ContractionModelInputParameters_::voltage.

void Nash2004ContractionModel::SetStretchAndStretchRate ( double  stretch,
double  stretchRate 
) [inline, virtual]

Neither stretch nor stretch rate are used so this method does nothing

Parameters:
stretchstretch
stretchRatestretch rate

Implements AbstractContractionModel.

Definition at line 108 of file Nash2004ContractionModel.hpp.


Member Data Documentation

const double Nash2004ContractionModel::e0ByT0 = 1.0/100 [static, private]

Other parameter. See above and reference.

Definition at line 66 of file Nash2004ContractionModel.hpp.

Referenced by EvaluateYDerivatives().

const double Nash2004ContractionModel::kTa = 47.9 [static, private]

Stiffness parameter. See reference. kPa

Definition at line 63 of file Nash2004ContractionModel.hpp.

Referenced by EvaluateYDerivatives().

Non-dimensionalised voltage. See above.

Definition at line 69 of file Nash2004ContractionModel.hpp.

Referenced by EvaluateYDerivatives(), Nash2004ContractionModel(), and SetInputParameters().


The documentation for this class was generated from the following files: