Airfoil2D.h

Go to the documentation of this file.

/*--------------------------------*- VM2D -*-----------------*---------------*\
| ##  ## ##   ##  ####  #####   |                            | Version 1.14   |
| ##  ## ### ### ##  ## ##  ##  |  VM2D: Vortex Method       | 2026/03/06     |
| ##  ## ## # ##    ##  ##  ##  |  for 2D Flow Simulation    *----------------*
|  ####  ##   ##   ##   ##  ##  |  Open Source Code                           |
|   ##   ##   ## ###### #####   |  https://www.github.com/vortexmethods/VM2D  |
|                                                                             |
| Copyright (C) 2017-2026 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"
 
#ifdef USE_CUDA
#ifdef USE_CUDNN
    #include "cuMlp.cuh"
#endif
#endif
 
namespace VM2D
{
 
    class AirfoilGeometry
    {
    protected:
        std::vector<Point2D> r_;

        std::vector<Point2D> v_;
 
    public:
        bool inverse;

        std::vector<Point2D> nrm;

        std::vector<std::pair<Point2D, Point2D>> psn;

        std::vector<Point2D> tau;

        std::vector<double> len;

        Point2D rcm;

        double phiAfl;

        double area;
 

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

        Point2D& setR(size_t q)
        {
            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(); };
    };
 
 
 
 
    class World2D;
    class WakeDataBase;

    class Airfoil : public AirfoilGeometry
    {
    public:
        const World2D& W;

        const size_t numberInPassport;
 
 

        double m;

        double J;

        Point2D startOfWays;
        std::vector<std::vector<Point2D>> possibleWays;
        
        std::vector<int> wayToVertex; 

        void lightningTest();
 

        std::vector<double> meanEpsOverPanel;
 
        

        IFCUDA(mutable double* devRPtr);

        IFCUDA(mutable double* devPsnPtr);

        IFCUDA(mutable double* devRhsPtr);

        IFCUDA(mutable double* devRhsLinPtr);

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

        IFCUDA(mutable double* devFreeVortexSheetPtr);
        
        IFCUDA(mutable double* devFreeVortexSheetLinPtr);

        IFCUDA(mutable double* devAttachedVortexSheetPtr);
        
        IFCUDA(mutable double* devAttachedVortexSheetLinPtr);

        IFCUDA(mutable double* devAttachedSourceSheetPtr);
        
        IFCUDA(mutable double* devAttachedSourceSheetLinPtr);

        IFCUDA(mutable double* devMeanEpsOverPanelPtr);

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

#ifdef USE_CUDNN
        IFCUDA(CudaNN NN);
 
        IFCUDA(std::vector<int> vorticesToProcess);
        IFCUDA(int vtp);
        IFCUDA(std::vector<float> batch);     //(vtp * 3 * 5);
        IFCUDA(std::vector<int> signy);     //(vtp * 3 * 5);
        IFCUDA(std::vector<float> nnout);     //(vtp * 2 * 5);
        IFCUDA(std::vector<int> batchPanels); // (vtp * 5);
#endif

        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;

        void CalcNrmTauLen();

        virtual void Rotate(double alpha);

        virtual void Scale(const Point2D&);

        virtual void Move(const Point2D& dr);

        virtual void GetGabarits(double gap = 0.02);

        bool isInsideGabarits(const Point2D& r) const;

        bool isOutsideGabarits(const Point2D& r) const;

        virtual bool IsPointInAirfoil(const Point2D& point) const;

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

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

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

        void calcMeanEpsOverPanel();
 

        virtual void GetDiffVelocityI0I3ToSetOfPointsAndViscousStresses(const WakeDataBase& pointsDb, std::vector<double>& domainRadius, std::vector<double>& I0, std::vector<Point2D>& I3);
 
#if defined(USE_CUDA)
        virtual void GPUGetDiffVelocityI0I3ToSetOfPointsAndViscousStresses(std::vector<double>& domainRadius, std::vector<double>& I0, std::vector<Point2D>& I3);
#endif
        
        virtual void GetInfluenceFromVorticesToPanel(size_t panel, const Vortex2D* ptr, ptrdiff_t count, std::vector<double>& panelRhs) const;
                
    
        virtual void GetInfluenceFromSourcesToPanel(size_t panel, const Vortex2D* ptr, ptrdiff_t count, std::vector<double>& panelRhs) const;
 

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

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

        virtual void GetInfluenceFromVInfToPanel(std::vector<double>& vInfRhs) const;
};
 
}//namespace VM2D
 
#endif