lumpedHeart.hpp

Go to the documentation of this file.

//@HEADER
/*
*******************************************************************************

    Copyright (C) 2004, 2005, 2007 EPFL, Politecnico di Milano, INRIA
    Copyright (C) 2010 EPFL, Politecnico di Milano, Emory University

    This file is part of LifeV.

    LifeV is free software; you can redistribute it and/or modify
    it under the terms of the GNU Lesser General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    LifeV is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    Lesser General Public License for more details.

    You should have received a copy of the GNU Lesser General Public License
    along with LifeV.  If not, see <http://www.gnu.org/licenses/>.

*******************************************************************************
*/
//@HEADER

/*!
 *  @file
 *  @brief File containing the lumped heart for the Monolithic Test
 *
 *  @date 2009-06-03
 *  @author Paolo Crosetto <crosetto@iacspc70.epfl.ch>
 *
 *  @contributor Cristiano Malossi <cristiano.malossi@epfl.ch>
 *  @maintainer Paolo Crosetto <crosetto@iacspc70.epfl.ch>
 */

#ifndef __LUMPEDHEART_HPP
#define __LUMPEDHEART_HPP

// LifeV includes
#include <lifev/core/LifeV.hpp>
#include <lifev/core/filter/GetPot.hpp>
#include <lifev/core/fem/BCHandler.hpp>
#include <lifev/fsi/solver/FSISolver.hpp>
#include <lifev/fsi/solver/FSIOperator.hpp>

// Mathcard includes
//#include <life/lifesolver/BCInterface3D.hpp>

namespace LifeV
{

class LumpedHeart
{
public:
    //    typedef BCInterface3D< BCHandler, FSIOperator >                                        bc_type;
    typedef FSIOperator::vector_Type                                                       vector_Type;
    typedef FSIOperator::vectorPtr_Type                                                    vectorPtr_Type;

    LumpedHeart()
        :
        M_time (0.),
        //M_BC(),
        M_ODEscheme(),
        M_dt (0.),
        M_T_max()  ,
        M_E_max()  ,
        M_V_0()    ,
        M_RV_art() ,
        M_RA_V()   ,
        M_LV_art() ,
        M_LA_V()   ,
        M_PV(),
        M_intFlux(),
        M_Vt_ao()
    {
        M_ODEscheme.setup (1);
    }

    void initParameters      ( FSIOperator&  Oper,
                               const std::string&    FileName);

    static Real& outPressure         (const Real& t, const Real& x, const Real& y, const Real& z, const ID& i);

    void renewParameters     ( FSIOperator&  Oper , const int& flag, const Real& time, const Real& flux);

    Real fZero               (const Real& t, const Real& x, const Real& y, const Real& z, const ID& i);

    static Real                               M_pressure;

private:

    //! @name Private Methods
    //@{

    Real                                      M_elastance (const Real& t);

    //! Short description of this method
    /*!
        Add more details about the method.
        NOTE: short description is automatically added before this part.
     */


    //@}

    Real                                      M_time;
    //std::shared_ptr< bc_type >              M_BC;
    TimeAdvanceBDF<Real>                      M_ODEscheme;
    Real                                      M_dt;
    Real                                      M_T_max  ;
    Real                                      M_E_max  ;
    Real                                      M_V_0    ; //reference ventricular volume
    Real                                      M_RV_art ;
    Real                                      M_RA_V   ;
    Real                                      M_LV_art ;
    Real                                      M_LA_V   ;
    Real                                      M_PV     ; //ventricular pressure
    Real                                      M_intFlux;
    Real                                      M_Vt_ao;
    Real                                      M_Tpb  ;
    Real                                      M_Tpw  ;
};
}

#endif /* __LUMPEDHEART_HPP */