---

momentum_u.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_MOMENTUM_U_H__
#define __CSFL_MATH_EQUATION_MOMENTUM_U_H__


// Include
#include <csfl/math/equation/sinflow.h>


// Namespace
namespace csfl {

//==============================================================================
//              Class IEquationMomentumU
//              Description: Class for Sinflow's Momentum Equation
//==============================================================================

class IEquationMomentumU : 
        public IEquationSinflow

{
        
public:
        IEquationMomentumU( 
                const IString & _name,
                IGrid *_grid, 
                IContainerBCond *_contBCond, 
                IContainerParameter *_contParam,
                IFieldScalar *_pre = NULL,             //Initial Pressure Field   P
        IFieldVector *_vel = NULL,             //Initial Velocity Field   u & v
                TEquationType _et = etMomU );

        IEquationMomentumU( IEquationMomentumU *_e );   
        ~IEquationMomentumU();
        
        IFieldScalar * Pressure() { return pre; } 
        
        void SetPressure( IFieldScalar *_pre );

    void SetVelocity( Vector *_v );

    void SetFieldVelocity( IFieldVector *_v );  

        void SetBoundaryConditions();

        void SetCoefficients( double _timestep );

        inline IParamPhysical * Mi()   { return mi; }

        inline IParamPhysical * Rho()  { return rho; }

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


protected:
        IParamPhysical  *mi, *rho;      
        IFieldScalar    *pre;       

        void SetupParameters();
        void SetupParametersSolidBlocks();

    void CalculateSourceTerm( double, int, int );

    double Phi( const IPoint & );
        double GammaPhi( const IPoint &, TDirection = diNULL );
                                
        void SetBCondEdges();
        void SetBCondNorth();
        void SetBCondSouth();
        void SetBCondEast();
        void SetBCondWest();

};
} // csfl namespace

#endif// __CSFL_MATH_EQUATION_MOMENTUM_U_H__

 

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