---

sinflow.h

Go to the documentation of this file.

/*
 * Copyright (c) 2001 The AUTHORS 
 *   Romeu Andre' Pieritz, Ph.D.        - romeu_pieritz@hotmail.com
 *   Rafael Mendes, Eng.                – mendes_rafael@yahoo.com
 *   Rodrigo Ferraz de Andrade, Eng.    – rferraz@iname.com 
 * All rights reserved.
 *
 * Permission to use, copy and distribute this software and its
 * documentation for educational and personal use, without fee is hereby granted, 
 * provided that the above copyright notice and the following
 * two paragraphs appear in all copies of this software.
 *
 * IN NO EVENT SHALL THE AUTHORS BE LIABLE TO ANY PARTY FOR
 * DIRECT, INDIRECT, SPECIAL, INCIDENTAL, OR CONSEQUENTIAL DAMAGES ARISING OUT
 * OF THE USE OF THIS SOFTWARE AND ITS DOCUMENTATION, EVEN IF THE AUTHORS
 * HAS BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 *
 * THE AUTHORS SPECIFICALLY DISCLAIMS ANY WARRANTIES,
 * INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY
 * AND FITNESS FOR A PARTICULAR PURPOSE.  THE SOFTWARE PROVIDED HERE UNDER IS
 * ON AN "AS IS" BASIS, AND THE AUTHORS HAVE NO OBLIGATION TO
 * PROVIDE MAINTENANCE, SUPPORT, UPDATES, ENHANCEMENTS, OR MODIFICATIONS.
 *
 * SINMEC Lab. - CFD Sinflow Project - http://www.sinmec.ufsc.br/cfd
 */


// File Define
#ifndef __CSFL_MATH_EQUATION_SINFLOW_H__
#define __CSFL_MATH_EQUATION_SINFLOW_H__


// Include
#include <csfl/math/equation/eqparam/physical.h>
#include <csfl/math/equation/eqparam/interpolparam.h>
#include <csfl/math/equation/equation.h>
#include <csfl/field/value/scalar.h>
#include <csfl/field/value/vector.h>
#include <csfl/field/value/vector_face.h>
#include <csfl/field/value/scalar_face.h>
#include <csfl/sys/math/linearsys9d.h>



// Namespace
namespace csfl {

//==============================================================================
//              Class IEquationSinflow
//              Description: Base Class for Sinflow's Equations
//==============================================================================

class IEquationSinflow : 
        public IEquation


{
        
public:
        IEquationSinflow( 
                const IString &_name,
                IGrid *_grid, 
                IContainerBCond *_contBCond, 
                IContainerParameter *contParam,
        IFieldVector *_vel,
                TEquationType _et = etNULL );

        IEquationSinflow( IEquationSinflow *_e );       
        ~IEquationSinflow();

        IParamInterpol * Interpol() { return interpol; }
        
        IGrid * Grid() { return grid; }

        IMetric * Metric() { return metric; }

        ILinearSystem * LinearSystem() { return lsys; }
 
        inline void  SetLinearSystem( ILinearSystem *_lsys ); 

        virtual void SetFieldVelocity( IFieldVector *_v );  
        
    inline IFieldVector * Velocity() { return vel; }   

    virtual void SetVelocity( Vector * ) {}

    virtual void SetBoundaryConditions();
        
    virtual void SetCoefficients( double = 1.0 ) {}

        virtual void SetCoefficientsSolidBlocks( double = 1.0 ) {}

    virtual void CalculateSourceTerm( double = 1.0 ) {}

    double VelU( IFieldVectorFace *_fvel, const IPoint & _p, TDirection _dir = diNULL );
 
        inline double VelU( IFieldVectorFace *_fvel, int _i, int _j, TDirection _dir = diNULL );

        double VelV( IFieldVectorFace *_fvel, const IPoint &_p, TDirection _dir = diNULL );
 
        inline double VelV( IFieldVectorFace *_fvel, int _i, int _j, TDirection _dir = diNULL );

        inline IFieldVectorFace * FaceVelocity() { return fvel; }

    inline void SetFaceVelocity( IFieldVectorFace *_fv );

        void InitFaceVelocity( IFieldVector *_vel );

        void SetInitialFieldVelocity( IFieldVector *_v );

        void SetFaceVelBCondPhi( IFieldVectorFace *_fvel, IFieldVector *_vel );

        void SetFaceVelBCondOutlet( IFieldVectorFace *_fvel, IFieldVector *_vel );

        void SetFaceVelBCondSymmetryX( IFieldVectorFace *_fvel, IFieldVector *_vel );

        void SetFaceVelBCondSymmetryY( IFieldVectorFace *_fvel, IFieldVector *_vel );

        void SetFaceVelSolidBlock( IFieldVectorFace *_fvel, IFieldVector *_vel  );

        void SetFaceVelBCondWall( IFieldVectorFace *_fvel, IFieldVector *_vel  );
        
        void SetFaceVelBCondOpenWall( IFieldVectorFace *_fvel, IFieldVector *_vel  );


protected:
        IParamInterpol   *interpol;     
        ICoefficient     coef;      
    IFieldVector         *vel, *vel0;              
    IFieldVectorFace *fvel;    
        ILinearSystem9D  *lsys; 
        IGrid               *grid;
        IMetric             *metric;

        bool isOutlet, isInlet, isPhi, isSolidBlock, isWall;
        bool isOpenWall, isSymmetryX, isSymmetryY;
        bool isPressureN, isPressureS, isPressureE, isPressureW;

        double Pn, Ps, Pe, Pw;
         
    virtual void SetupParameters(){}
    virtual void SetupParametersSolidBlocks() {}
  
    virtual double GammaPhi( const IPoint &, TDirection = diNULL )
        { return 0.0; }
                
        virtual void SetBCondEdges() {}
        virtual void SetBCondNorth() {}
        virtual void SetBCondSouth() {}
        virtual void SetBCondEast() {}
        virtual void SetBCondWest() {}
        
        void SetSolidBlockDomainRegionID();
  
};

// =============================================================================
// Inline functions
// =============================================================================

inline void IEquationSinflow::SetFaceVelocity( IFieldVectorFace *_fv )
{
        int i, j;
    IVector v;
    for( i=0; i<vel->Size().width+1; i++ )
        for( j=0; j<vel->Size().height+1; j++ ) {
                        v = _fv->Vector( i, j, diSouth );
            fvel->SetVector( i, j, diSouth, v );
                        v = _fv->Vector( i, j, diWest );
            fvel->SetVector( i, j, diWest, v );
        }   
}    
    
inline void  IEquationSinflow::SetLinearSystem( ILinearSystem * _lsys )
{
        if(!_lsys)
                return;
        
        ILinearSystem9D *l9d = dynamic_cast<ILinearSystem9D * >(_lsys);
        if (l9d){
                if (lsys) 
                        delete lsys; 
                lsys = new ILinearSystem9D( l9d ); 
        }
}

inline double IEquationSinflow::VelU( IFieldVectorFace *_fvel, int _i, int _j, TDirection _dir )
{
    return VelU( _fvel, IPoint( _i, _j ), _dir );
}

inline double IEquationSinflow::VelV( IFieldVectorFace *_fvel, int _i, int _j, TDirection _dir )
{
    return VelV( _fvel, IPoint( _i, _j ), _dir );
}
} // csfl namespace

#endif// __CSFL_MATH_EQUATION_SINFLOW_H__

 

CFD Project | Documents | Downloads | Contact us | Use Terms
© SINMEC/EMC/UFSC, 2001.	All rights reserved.	Last Update: Jan. 18, 2002	Webmaster