An Object Finite Element Library

DC2DT3.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 class DC2DT3
   for 2-D diffusion-convection scalar equation using 3-node triangular element

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


#ifndef __DC2DT3_H
#define __DC2DT3_H


#include "Equa_Therm.h"
#include "Triang3.h"
#include "Line2.h"

namespace OFELI {

class DC2DT3 : public Equa_Therm<double,3,3,2,2>
{

 public :

    DC2DT3();

    DC2DT3(const Element *el);

    DC2DT3(const Side *sd);

    DC2DT3(const Element *el, const Vect<double> &u, double time=0.);

    DC2DT3(const Element *el, const Vect<double> &u, double time, double deltat, int scheme);

    DC2DT3(const Side *sd, const Vect<double> &u, double time=0.);

    DC2DT3(const Side *sd, const Vect<double> &u, double time, double deltat, int scheme);

#ifndef DOXYGEN_SHOULD_SKIP_THIS
    DC2DT3(Mesh &mesh, Vect<double> &b, double &t, double &ts) : Equation<double,3,3,2,2>(mesh,b,t,ts) { }
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

    ~DC2DT3();

#ifndef DOXYGEN_SHOULD_SKIP_THIS
//    void withConvection(int fl) { fl=0; _with_convection=true; }
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

    int runOneStep();

    int run() { return runOneStep(); }

    void build();

    void LCapacityToLHS(double coef=1);

    void LCapacityToRHS(double coef=1);

    void LCapacity(double coef) { LCapacityToLHS(coef); LCapacityToRHS(coef); }

    void CapacityToLHS(double coef=1);

    void CapacityToRHS(double coef=1);

    void Capacity(double coef=1);

    void Diffusion(double coef=1);

    void Diffusion(const DMatrix<double> &diff, double coef=1);

    void DiffusionToRHS(double coef=1);

    void Convection(const Point<double> &v, double coef=1);

    void Convection(const Vect<double> &v, double coef=1);

    void Convection(double coef=1);

    void ConvectionToRHS(const Point<double> &v, double coef=1);

    void ConvectionToRHS(double coef=1);

    void setVelocity(NodeVect<double> &v);

    void BodyRHS(UserData<double> &ud, double coef=1);

    void BodyRHS(const Vect<double> &b, int opt=LOCAL_ARRAY);

    void BoundaryRHS(UserData<double> &ud, double coef=1);

    void BoundaryRHS(double flux);

    void BoundaryRHS(const Vect<double> &b, int opt=LOCAL_ARRAY);

    void Periodic(double coef=1.e20);

#ifndef DOXYGEN_SHOULD_SKIP_THIS
    double Energy(const Vect<double> &u, UserData<double> &ud);
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

#ifndef DOXYGEN_SHOULD_SKIP_THIS
    void EnerGrad(const Vect<double> &u, UserData<double> &ud, Vect<double> &g);
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

    Point<double> &Flux() const;

    Point<double> &Grad(const LocalVect<double,3> &u) const;

    Point<double> & Grad(const Vect<double> &u) const;

    void setInput(int opt, Vect<double> &u);

 protected :
   void ElementSet(const Element *el);
   void SideSet(const Side *sd);

 private :
   Triang3               *_tr;
   Line2                 *_ln;
   mutable Point<double> _f;
   Vect<double>          _vx, _vy;
};

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani