fsi/examples/example_SmoothAneurysm/flowConditions.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 boundary conditions for the Monolithic Test
 *
 *  @date 2009-04-09
 *  @author Paolo Crosetto <crosetto@iacspc70.epfl.ch>
 *
 *  @contributor Cristiano Malossi <cristiano.malossi@epfl.ch>
 *  @maintainer Paolo Crosetto <crosetto@iacspc70.epfl.ch>
 */

#ifndef __FLOWCONDITIONS_HPP
#define __FLOWCONDITIONS_HPP


#include <Epetra_SerialDenseVector.h>


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

namespace LifeV
{

class FlowConditions
{
public:

    FlowConditions();

    void initParameters      ( FSIOperator&  oper, const int&    outflowFlag);

    void renewParameters     ( FSISolver&  oper, const int&    outflowFlag, const FSIOperator::vector_Type& fluidSolution);

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

    static Real outPressure0         (const Real& t, const Real& x, const Real& y, const Real& z, const ID& i);
    static Real outPressure1         (const Real& t, const Real& x, const Real& y, const Real& z, const ID& i);
    static Real outPressure2         (const Real& t, const Real& x, const Real& y, const Real& z, const ID& i);
    static Real outPressure3         (const Real& t, const Real& x, const Real& y, const Real& z, const ID& i);
    static Real outPressure4         (const Real& t, const Real& x, const Real& y, const Real& z, const ID& i);
    static Real outPressure5         (const Real& t, const Real& x, const Real& y, const Real& z, const ID& i);
    static Real outPressure6         (const Real& t, const Real& x, const Real& y, const Real& z, const ID& i);

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


private:
    Real pi;

    bool bcOnFluid;

    Real M_outflux;
    Real M_influx;
    Real M_outP;

    Real M_area0;
    Real M_inRadius0;
    Real M_outRadius0;
    Real M_inDeltaRadius;
    Real M_outDeltaRadius;

    Real M_beta;
    Real M_rhos;
    static    std::vector<Real> outputVector;
    UInt conditionNumber;
};
}

#endif /* __FLOWCONDITIONS_HPP */