Velocity2D.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: Velocity2D.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 VELOCITY2D_H
#define VELOCITY2D_H
 
#include <memory>
 
#include "defs.h"
#include "cudaTreeInfo.h"
#include "OptimizedVelocity2D.h"
 
namespace VM2D
{
 
    class Airfoil;
    class Boundary;
    class WakeDataBase;
    class World2D;
    
 
    struct VortexesParams
    {
        std::vector<Point2D> convVelo;
 
        std::vector<Point2D> diffVelo;
 
        std::vector<double> I0;
 
        std::vector<double> I1;
 
        std::vector<Point2D> I2;
 
        std::vector<Point2D> I3;
 
        std::vector<double> epsastWake;
    };
 
    class Velocity
    {
    protected:
        const World2D& W;
 
    public:
        VortexesParams wakeVortexesParams;
 
        std::vector<VortexesParams> virtualVortexesParams;
 
        OptimizedVelocity optimizedVelocity;
 
        Velocity(const World2D& W_) :
            W(W_)
        {
            virtualVortexesParams.resize(0);
        };
 
        void CalcConvVelo();
 
        virtual void CalcConvVeloToSetOfPointsFromWake(const WakeDataBase& pointsDb, std::vector<Point2D>& velo, std::vector<double>& domainRadius, bool calcVelo, bool calcRadius) = 0;
        virtual void CalcConvVPVeloToSetOfPointsFromWake(const WakeDataBase& pointsDb, std::vector<Point2D>& velo, std::vector<double>& domainRadius, bool calcVelo, bool calcRadius) {};
 
#if defined (USE_CUDA)
        virtual void GPUCalcConvVeloToSetOfPointsFromWake(std::unique_ptr<BHcu::CudaTreeInfo>& cntrTree, const WakeDataBase& pointsDb, std::vector<Point2D>& velo, std::vector<double>& domainRadius, bool calcVelo, bool calcRadius) = 0;
        virtual void GPUCalcConvVelocityToSetOfPointsFromSheets(std::unique_ptr<BHcu::CudaTreeInfo>& cntrTree, const WakeDataBase& pointsDb, std::vector<Point2D>& velo) const = 0;
#endif
 
        virtual void CalcVeloToWakeVP();
 
        void CalcDiffVeloI1I2ToSetOfPointsFromWake(const WakeDataBase& pointsDb, const std::vector<double>& domainRadius, const WakeDataBase& vorticesDb, std::vector<double>& I1, std::vector<Point2D>& I2);
        void CalcDiffVeloI1I2ToSetOfPointsFromSheets(const WakeDataBase& pointsDb, const std::vector<double>& domainRadius, const Boundary& bnd, std::vector<double>& I1, std::vector<Point2D>& I2);
        
        void CalcDiffVeloI1I2ToWakeFromSheets(const WakeDataBase& pointsDb, const std::vector<double>& domainRadius, const Boundary& bnd, std::vector<double>& I1, std::vector<Point2D>& I2);
        void CalcDiffVeloI1I2ToWakeFromWake(const WakeDataBase& pointsDb, const std::vector<double>& domainRadius, const WakeDataBase& vorticesDb, std::vector<double>& I1, std::vector<Point2D>& I2);
#if defined(USE_CUDA)
        void GPUCalcDiffVeloI1I2ToSetOfPointsFromWake(const WakeDataBase& pointsDb, const std::vector<double>& domainRadius, const WakeDataBase& vorticesDb, std::vector<double>& I1, std::vector<Point2D>& I2, bool useMesh = false);
        void GPUCalcDiffVeloI1I2ToSetOfPointsFromSheets(const WakeDataBase& pointsDb, const std::vector<double>& domainRadius, const Boundary& bnd, std::vector<double>& I1, std::vector<Point2D>& I2, bool useMesh = false);
        void GPUDiffVeloFAST(const std::vector<double>& domainRadius, std::vector<double>& I1, std::vector<Point2D>& I2);
#endif
 
        void CalcDiffVeloI1I2();
        void CalcDiffVeloI0I3();
 
 
        void LimitDiffVelo(std::vector<Point2D>& diffVel);
 
        void CalcDiffVelo();
 
        void ResizeAndZero();
 
        void SaveVisStress();
 
 
        void GetWakeInfluenceToRhs(const Airfoil& afl, std::vector<double>& wakeRhs) const;
#if defined(USE_CUDA)
        void GPUGetWakeInfluenceToRhs(const Airfoil& afl, std::vector<double>& wakeVelo) const;
        void GPUFASTGetWakeInfluenceToRhs(const Airfoil& afl, std::vector<double>& wakeVelo) const;
#endif
        void FillRhs(/*Eigen::VectorXd& rhs, Eigen::VectorXd& rhsSkos,*/ Eigen::VectorXd& rhsReord) const;
 
 
        virtual ~Velocity() { };
    };
 
}//namespace VM2D
 
#endif