Airfoil2D.h

Go to the documentation of this file.

/*--------------------------------*- VM2D -*-----------------*---------------*\
| ##  ## ##   ##  ####  #####   |                            | Version 1.12   |
| ##  ## ### ### ##  ## ##  ##  |  VM2D: Vortex Method       | 2024/01/14     |
| ##  ## ## # ##    ##  ##  ##  |  for 2D Flow Simulation    *----------------*
|  ####  ##   ##   ##   ##  ##  |  Open Source Code                           |
|   ##   ##   ## ###### #####   |  https://www.github.com/vortexmethods/VM2D  |
|                                                                             |
| Copyright (C) 2017-2024 I. Marchevsky, K. Sokol, E. Ryatina, A. Kolganova   |
*-----------------------------------------------------------------------------*
| File name: Airfoil2D.h                                                      |
| Info: Source code of VM2D                                                   |
|                                                                             |
| This file is part of VM2D.                                                  |
| VM2D 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, either version 3 of the License, or           |
| (at your option) any later version.                                         |
|                                                                             |
| VM2D 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 VM2D.  If not, see <http://www.gnu.org/licenses/>.               |
\*---------------------------------------------------------------------------*/


#ifndef AIRFOIL_H
#define AIRFOIL_H

#include "defs.h"

namespace VM2D
{

    class World2D;
    class WakeDataBase;

    class Airfoil
    {
    protected:
        std::vector<Point2D> r_;

        std::vector<Point2D> v_;

    public:
        const World2D& W;

        const size_t numberInPassport;

        Point2D rcm;

        double phiAfl;

        double area;

        double m;

        double J;

        std::vector<double> meanEpsOverPanel;

        const Point2D& getR(size_t q) const
        {
            return (q < r_.size()) ? r_[q] : r_[0];
        };

        const Point2D& getV(size_t q) const
        {
            return (q < v_.size()) ? v_[q] : v_[0];
        };

        void setV(const Point2D& vel)
        {
            v_.clear();
            v_.resize(r_.size(), vel);
        }

        void setV(const std::vector<Point2D>& vel)
        {
            v_.clear();
            v_.insert(v_.end(), vel.begin(), vel.end());                
        }

        size_t getNumberOfPanels() const { return r_.size(); };
        
        bool inverse;

        IFCUDA(mutable double* devRPtr);

        IFCUDA(mutable double* devRhsPtr);

        IFCUDA(mutable double* devRhsLinPtr);

        IFCUDA(mutable std::vector<double> tmpRhs);

        IFCUDA(mutable double* devFreeVortexSheetPtr);

        IFCUDA(mutable double* devAttachedVortexSheetPtr);

        IFCUDA(mutable double* devAttachedSourceSheetPtr);

        IFCUDA(mutable double* devMeanEpsOverPanelPtr);

        IFCUDA(mutable double* devViscousStressesPtr);
               
        IFCUDA(mutable std::vector<double> tmpViscousStresses);

        std::vector<Point2D> nrm;

        std::vector<Point2D> tau;

        std::vector<double> len;

        std::vector<double> viscousStress;

        Point2D lowLeft; 
        Point2D upRight; 

        std::vector<double> gammaThrough;

        Airfoil(const World2D& W_, const size_t numberInPassport_);


        virtual ~Airfoil() { };


        bool isAfter(size_t i, size_t j) const;

        virtual void Rotate(double alpha) = 0;

        virtual void Scale(const Point2D&) = 0;

        virtual void Move(const Point2D& dr) = 0;

        virtual void GetGabarits(double gap = 0.02) = 0;

        bool isInsideGabarits(const Point2D& r) const;

        bool isOutsideGabarits(const Point2D& r) const;

        virtual bool IsPointInAirfoil(const Point2D& point) const = 0;

        virtual void ReadFromFile(const std::string& dir) = 0;

        virtual std::vector<double> getA(size_t p, size_t i, const Airfoil& airfoil, size_t j) const = 0;

        virtual void calcIQ(size_t p, const Airfoil& otherAirfoil, std::pair<Eigen::MatrixXd, Eigen::MatrixXd>& matrPair) const = 0;


        void calcMeanEpsOverPanel();


        virtual void GetDiffVelocityI0I3ToSetOfPointsAndViscousStresses(const WakeDataBase& pointsDb, std::vector<double>& domainRadius, std::vector<double>& I0, std::vector<Point2D>& I3) = 0;
#if defined(USE_CUDA)
        virtual void GPUGetDiffVelocityI0I3ToSetOfPointsAndViscousStresses(const WakeDataBase& pointsDb, std::vector<double>& domainRadius, std::vector<double>& I0, std::vector<Point2D>& I3) = 0;
#endif
        virtual void GetDiffVelocityI0I3ToWakeAndViscousStresses(const WakeDataBase& pointsDb, std::vector<double>& domainRadius, std::vector<double>& I0, std::vector<Point2D>& I3) = 0;

        virtual void GetInfluenceFromVorticesToPanel(size_t panel, const Vortex2D* ptr, ptrdiff_t count, std::vector<double>& panelRhs) const = 0;
                
    
        virtual void GetInfluenceFromSourcesToPanel(size_t panel, const Vortex2D* ptr, ptrdiff_t count, std::vector<double>& panelRhs) const = 0;


        virtual void GetInfluenceFromSourceSheetToVortex(size_t panel, const Vortex2D& vtx, Point2D& vel) const = 0;

        virtual void GetInfluenceFromVortexSheetToVortex(size_t panel, const Vortex2D& vtx, Point2D& vel) const = 0;

        virtual void GetInfluenceFromVInfToPanel(std::vector<double>& vInfRhs) const = 0;


};

}//namespace VM2D

#endif