An Object Finite Element Library

NSP2DQ41.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 NSP2DQ41
             for 2-D Navier-Stokes equations using penalty formulation and
                             Quadrilateral finite element
                          (4-nodes velocity, 1-node pressure)

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


#ifndef __NSP2DQ41_H
#define __NSP2DQ41_H

#include "Equa_Fluid.h"
#include "Quad4.h"
#include "Line2.h"
#include "UserData.h"

namespace OFELI {

class NSP2DQ41  : virtual public Equa_Fluid<double,4,8,2,4>
{

public :

    NSP2DQ41()
    {
       _quad = NULL;
       _ln = NULL;
#ifdef _OFELI_DEBUG
       std::clog << "An instance of class NSP2DQ41 is constructed.\n";
       std::clog << "File : " << __FILE__ << ", Line : " << __LINE__ << endl;
#endif
    }

    NSP2DQ41(const Element *el);

    NSP2DQ41(const Side *sd);

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

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

    ~NSP2DQ41();

    void Viscosity(double visc) { _visc = visc; }

    void Density(double dens) { _dens = dens; }

    void LMass(double coef=1.);

    void Mass(double coef=1.);

    void Viscous(double coef=1.);

    void RHS_Viscous(double coef=1.);

    void Penal(double coef=1.);

    void LHS1_Convection(double coef=1.);

    void LHS2_Convection(double coef=1.);

    void RHS_Convection(double coef=1.);

    void BodyRHS(UserData<double>  &ud);

    void BoundaryRHS(UserData<double> &ud);

    void Periodic(double coef=1.e20);

    double Pressure(double coef=1.);

    void PresMat();

    LocalVect<double,4> PM() { return _PM; }

 private :

   Quad4               *_quad;
   Line2               *_ln;
   double              _gauss[2], _xl[4], _yl[4];
   double              _cgx, _cgy;
   LocalVect<double,4> _PM;
};

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani