---

atraction.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_EQPARAM_ATRACTION_H__
#define __CSFL_MATH_EQUATION_EQPARAM_ATRACTION_H__


// Include
#include <csfl/geom/metric/metric.h>


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


// Namespace
namespace csfl {


//==============================================================================
//      enum TLineAtracted
//==============================================================================
typedef
        enum _TLineAtracted {
                laKsi,
                laEta
        }
        TLineAtracted;

//==============================================================================
//              Class IAtraction
//              Description: Base class for mesh Atraction
//==============================================================================
class IAtraction :      public IParameter
{
public:
        IAtraction( const IString & = IString::empty,  
                TParameterType = ptAtract );
        IAtraction( IAtraction * );
        ~IAtraction();

public:
        virtual double EvaluateX( const int &, const int & ){ return 0.0; }
        virtual double EvaluateY( const int &, const int & ){ return 0.0; }

};

//==============================================================================
//              Class ILineAtraction
//              Description: class for line Atraction
//==============================================================================

class ILineAtraction :  public IAtraction
{
public:
        ILineAtraction( const int &, const TLineAtracted &, 
                const double &, const double & = 1.0, const IString & = IString::empty);
        ILineAtraction( ILineAtraction * );
        ~ILineAtraction();

public:
        virtual inline double EvaluateX( const int &, const int &);
        virtual inline double EvaluateY( const int &, const int &);

private:
        int line;
        TLineAtracted linetype;
        double a, c;

};

//==============================================================================
//              Class INodeAtraction
//              Description:  class for node Atraction
//==============================================================================

class INodeAtraction :  public IAtraction
{
public:
        INodeAtraction( const IPoint &, const double &, const double & = 1.0, 
                const IString & = IString::empty);
        INodeAtraction( INodeAtraction * );
        ~INodeAtraction();

public:
        virtual inline double EvaluateX( const int &, const int & );
        virtual inline double EvaluateY( const int &, const int & );

private:
        double b, d;
        IPoint patract;

};
// =============================================================================
// Inline functions
// =============================================================================

inline double ILineAtraction::EvaluateX( const int & ksi, const int & eta )
{
        if( linetype == laEta )
                return 0.0;

        else{
                if( ksi == line ) return 0.0;
                double expo = exp( c * fabs( double( ksi - line )));
                double A = a * ((ksi > line )? 1.0: -1.0 );
                return A / expo;
                //return(a * ((ksi > line )? 1.0: -1.0 ) / //(double( ksi - line )/fabs( double(ksi - line ))) / 
                //      exp( c * fabs( double( ksi - line ))));
        }
}
        
inline double ILineAtraction::EvaluateY( const int & ksi, const int & eta )
{
        if( linetype == laKsi )
                return 0.0;
        else{
                if( eta == line ) return 0.0;
                double expo = exp( c * fabs( double( eta - line )));
                double A = a * ((eta > line )? 1.0: -1.0 );
                return A / expo;
                //return( a * ((eta > line )? 1.0 : -1.0 ) / //( double( eta - line )/fabs(double(eta - line) )) / 
                //      exp( c * fabs( double( eta - line ))));
        }
}

inline double INodeAtraction::EvaluateX( const int & ksi, const int & eta )
{
        if( ksi == patract.x ) 
                return 0.0 ;//b / exp( d * ( eta - patract.y ));

        return( b * ((ksi > patract.x)? 1.0 : -1.0 ) / //(( double(ksi - patract.x ))/fabs( double( ksi - patract.x ))) / 
                        exp( d * sqrt( pow(double( ksi - patract.x ), 2.0 ) + 
                        pow(double( eta - patract.y ), 2.0 ))));
}

inline double INodeAtraction::EvaluateY( const int & ksi, const int & eta )
{
        if( eta == patract.y ) 
                return 0.0; //b / exp( d * ( ksi - patract.x ));

        return( b * ((eta > patract.y)? 1.0 : -1.0 ) / //( double( eta - patract.y )/fabs( double(eta - patract.y) )) / 
                        exp( d * sqrt( pow(double( ksi - patract.x ), 2.0 ) + 
                        pow(double( eta - patract.y ), 2.0 ))));
}
} // csfl namespace

#endif// __CSFL_MATH_EQUATION_EQPARAM_ATRACTION_H__

 

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