Chaste Release::3.1
ImplicitCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM > Class Template Reference

#include <ImplicitCardiacMechanicsSolver.hpp>

Inheritance diagram for ImplicitCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM >:
Collaboration diagram for ImplicitCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM >:

List of all members.

Public Member Functions

 ImplicitCardiacMechanicsSolver (ContractionModelName contractionModelName, QuadraticMesh< DIM > &rQuadMesh, ElectroMechanicsProblemDefinition< DIM > &rProblemDefinition, std::string outputDirectory)
virtual ~ImplicitCardiacMechanicsSolver ()
void Solve (double time, double nextTime, double odeTimestep)

Private Member Functions

bool IsImplicitSolver ()
void GetActiveTensionAndTensionDerivs (double currentFibreStretch, unsigned currentQuadPointGlobalIndex, bool assembleJacobian, double &rActiveTension, double &rDerivActiveTensionWrtLambda, double &rDerivActiveTensionWrtDLambdaDt)
void InitialiseContractionModels (ContractionModelName contractionModelName)

Friends

class TestImplicitCardiacMechanicsSolver
class TestExplicitCardiacMechanicsSolver

Detailed Description

template<class ELASTICITY_SOLVER, unsigned DIM>
class ImplicitCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM >

Implicit Cardiac Mechanics Solver

The first template parameter should be either IncompressibleNonlinearElasticitySolver or CompressibleNonlinearElasticityAssembler; this will be the class that this class ultimately inherits from.

Solves cardiac mechanics implicitly (together with the contraction models for determining the active tension), taking in the intracellular calcium concentration. See CardiacElectroMechanicsProblem documentation for more detail.

Definition at line 59 of file ImplicitCardiacMechanicsSolver.hpp.


Constructor & Destructor Documentation

template<class ELASTICITY_SOLVER , unsigned DIM>
ImplicitCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM >::ImplicitCardiacMechanicsSolver ( ContractionModelName  contractionModelName,
QuadraticMesh< DIM > &  rQuadMesh,
ElectroMechanicsProblemDefinition< DIM > &  rProblemDefinition,
std::string  outputDirectory 
)

Constructor

Parameters:
contractionModelNameThe contraction model.
rQuadMeshA reference to the mesh.
rProblemDefinitionObject defining body force and boundary conditions
outputDirectoryThe output directory, relative to TEST_OUTPUT

Definition at line 43 of file ImplicitCardiacMechanicsSolver.cpp.

template<class ELASTICITY_SOLVER , unsigned DIM>
ImplicitCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM >::~ImplicitCardiacMechanicsSolver ( ) [virtual]

Destructor

Definition at line 106 of file ImplicitCardiacMechanicsSolver.cpp.


Member Function Documentation

template<class ELASTICITY_SOLVER , unsigned DIM>
void ImplicitCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM >::GetActiveTensionAndTensionDerivs ( double  currentFibreStretch,
unsigned  currentQuadPointGlobalIndex,
bool  assembleJacobian,
double rActiveTension,
double rDerivActiveTensionWrtLambda,
double rDerivActiveTensionWrtDLambdaDt 
) [private, virtual]

A method called by AbstractCardiacMechanicsSolver::AssembleOnElement(), providing the active tension (and other info) at a particular quadrature point. This version uses C to determine the current stretch and stretch rate, and integrates the contraction model ODEs to determine the active tension, and derivatives of the active tension with respect to stretch and stretch rate.

Parameters:
currentFibreStretchThe stretch in the fibre direction
currentQuadPointGlobalIndexQuadrature point integrand currently being evaluated at in AssembleOnElement.
assembleJacobianA bool stating whether to assemble the Jacobian matrix.
rActiveTensionThe returned active tension
rDerivActiveTensionWrtLambdaThe returned dT_dLam, derivative of active tension wrt stretch
rDerivActiveTensionWrtDLambdaDtThe returned dT_dLamDot, derivative of active tension wrt stretch rate

Implements AbstractCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM >.

Definition at line 147 of file ImplicitCardiacMechanicsSolver.cpp.

References DataAtQuadraturePoint_::ContractionModel, EXCEPTION, AbstractContractionModel::GetNextActiveTension(), AbstractContractionModel::RunDoNotUpdate(), AbstractContractionModel::SetStretchAndStretchRate(), DataAtQuadraturePoint_::Stretch, and DataAtQuadraturePoint_::StretchLastTimeStep.

template<class ELASTICITY_SOLVER , unsigned DIM>
void ImplicitCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM >::InitialiseContractionModels ( ContractionModelName  contractionModelName) [private, virtual]

Initialise contraction models for each quadrature point

Parameters:
contractionModelNameThe name of the contraction model (from the enumeration ContractionModel defined in AbstractContactionModel)

Implements AbstractCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM >.

Definition at line 58 of file ImplicitCardiacMechanicsSolver.cpp.

References EXCEPTION.

template<class ELASTICITY_SOLVER, unsigned DIM>
bool ImplicitCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM >::IsImplicitSolver ( ) [inline, private, virtual]

This solver is an implicit solver (overloaded pure method)

Implements AbstractCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM >.

Definition at line 65 of file ImplicitCardiacMechanicsSolver.hpp.

template<class ELASTICITY_SOLVER , unsigned DIM>
void ImplicitCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM >::Solve ( double  time,
double  nextTime,
double  odeTimestep 
) [virtual]

Solve for the deformation using quasi-static nonlinear elasticity. (not dynamic nonlinear elasticity, despite the times taken in - just ONE deformation is solved for. The cell models are integrated implicitly over the time range using the ODE timestep provided, as part of the solve, and updated at the end once the solution has been found, as is lambda.

Parameters:
timethe current time
nextTimethe next time
odeTimestepthe ODE timestep

Implements AbstractCardiacMechanicsSolver< ELASTICITY_SOLVER, DIM >.

Definition at line 112 of file ImplicitCardiacMechanicsSolver.cpp.

References AbstractContractionModel::UpdateStateVariables().


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