An Object Finite Element Library

Equa_Therm.h

Go to the documentation of this file.

/*==============================================================================

                                    O  F  E  L  I

                           Object  Finite  Element  Library

  ==============================================================================

   Copyright (C) 1998 - 2008 Rachid Touzani

   This program is free software; you can redistribute it and/or modify it under
   the terms of the GNU General Public License as published by the Free
   Software Foundation; Version 2 of the License.

   This program 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 General Public License for more
   details.

   You should have received a copy of the GNU General Public License
   along with this program; if not, write to the :

   Free Software Foundation
   Inc., 59 Temple Place - Suite 330
   Boston, MA  02111-1307, USA

  ==============================================================================

    Definition of abstract class 'Equa_Therm' to solve heat transfer problems

  ==============================================================================*/


#ifndef __EQUA_THERM_H
#define __EQUA_THERM_H

#include "Equation.h"
#include "therm/PhaseChange.h"

namespace OFELI {

#ifndef DOXYGEN_SHOULD_SKIP_THIS
enum PDE_TERMS {
  CAPACITY            =   1,     
  LUMPED_CAPACITY     =   2,     
  DIFFUSION           =   4,     
  CONVECTION          =   8,     
  HEAT_SOURCE         =  12,     
  HEAT_FLUX           =  20      
};
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

template<class T_, size_t NEN_, size_t NEE_, size_t NSN_, size_t NSE_> class Equa_Therm;

template<class T_, size_t NEN_, size_t NEE_, size_t NSN_, size_t NSE_>
class Equa_Therm : virtual public Equation<T_,NEN_,NEE_,NSN_,NSE_>
{

 public:

   using Equation<T_,NEN_,NEE_,NSN_,NSE_>::setMaterialProperty;
   using Equation<T_,NEN_,NEE_,NSN_,NSE_>::_theElement;
   using Equation<T_,NEN_,NEE_,NSN_,NSE_>::_theSide;
   using Equation<T_,NEN_,NEE_,NSN_,NSE_>::_terms;
   using Equation<T_,NEN_,NEE_,NSN_,NSE_>::_analysis;
   using Equation<T_,NEN_,NEE_,NSN_,NSE_>::_time_scheme;

    Equa_Therm()
    {
       _terms = DIFFUSION;
       _analysis = STEADY_STATE;
               _time_scheme = STATIONARY;
    }

    virtual ~Equa_Therm()
    {
#ifdef _OFELI_DEBUG
      std::clog << "An instance of class Equa_Therm is destructed in\n";
      std::clog << "File: " << __FILE__ << ", Line: " << __LINE__ << endl;
#endif
    }

    virtual void LCapacityToLHS(double coef=1) { _coef = coef; }

    virtual void LCapacityToRHS(double coef=1) { _coef = coef; }

    virtual void CapacityToLHS(double coef=1) { _coef = coef; }

    virtual void CapacityToRHS(double coef=1) { _coef = coef; }

    void setLumpedCapacity() 
    {
       LCapacityToLHS(1./AbsEqua<T_>::_time_step);
       LCapacityToRHS(1./AbsEqua<T_>::_time_step);
    }

    void setCapacity() 
    {
       CapacityToLHS(1./AbsEqua<T_>::_time_step);
       CapacityToRHS(1./AbsEqua<T_>::_time_step);
    }

    virtual void Diffusion(double coef=1.) { _coef = coef; }

    virtual void DiffusionToRHS(double coef=1.) { _coef = coef; }

    void setDiffusion()
    {
       if (AbsEqua<T_>::_analysis==STEADY_STATE || AbsEqua<T_>::_time_scheme==BACKWARD_EULER)
          Diffusion(1.);
       else if (AbsEqua<T_>::_time_scheme==CRANK_NICOLSON)
          Diffusion(0.5);
       if (AbsEqua<T_>::_time_scheme==CRANK_NICOLSON)
          DiffusionToRHS(0.5);
       else if (AbsEqua<T_>::_time_scheme==FORWARD_EULER)
          DiffusionToRHS();
       else;
    }

    virtual void Convection(double coef=1.) { _coef = coef; }

    virtual void ConvectionToRHS(double coef=1.) { _coef = coef; }

    void setConvection()
    {
       if (AbsEqua<T_>::_time_scheme==STEADY_STATE || 
           AbsEqua<T_>::_time_scheme==BACKWARD_EULER)
          Convection(1.);
       else if (AbsEqua<T_>::_time_scheme==CRANK_NICOLSON)
          Convection(0.5);
       if (AbsEqua<T_>::_time_scheme==CRANK_NICOLSON)
          ConvectionToRHS(0.5);
       else if (AbsEqua<T_>::_time_scheme==FORWARD_EULER)
          ConvectionToRHS();
       else;
    }

#ifndef DOXYGEN_SHOULD_SKIP_THIS
    void Build() { }
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

 protected:

#ifndef DOXYGEN_SHOULD_SKIP_THIS
   PhaseChange *_phase;
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

    void RhoCp(const double &rhocp) { _rhocp = rhocp; }

    void Conduc(const double &diff) { _diff = diff; }

    void RhoCp(const string &exp) { _rhocp = setMaterialProperty(exp.c_str(),"RhoCp"); }

    void Conduc(const string &exp) { _diff = setMaterialProperty(exp.c_str(),"Conductivity"); }

    void setMaterial()
    {
       _rhocp = _theElement->getMaterial()->Density() * _theElement->getMaterial()->SpecificHeat();
       _diff = _theElement->getMaterial()->ThermalConductivity();
    }

// Coupling fields
#ifndef DOXYGEN_SHOULD_SKIP_THIS
   double       _diff, _rhocp, _coef;
   NodeVect<T_> *_velocity;
#endif /* DOXYGEN_SHOULD_SKIP_THIS */
};

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani