VM2D::BoundaryLinLayerAver Class Reference

Класс, определяющий способ удовлетворения граничного условия на обтекаемом профиле More...

#include <Boundary2DLinLayerAver.h>

Inheritance diagram for VM2D::BoundaryLinLayerAver:

Collaboration diagram for VM2D::BoundaryLinLayerAver:

Public Member Functions
	BoundaryLinLayerAver (const World2D &W_, size_t numberInPassport_)
	Конструктор More...
virtual	~BoundaryLinLayerAver ()
	Деструктор More...
virtual void	FillMatrixSelf (Eigen::MatrixXd &matr, Eigen::VectorXd &lastLine, Eigen::VectorXd &lactCol) override
	Генерация блока матрицы More...
virtual void	FillIQSelf (std::pair< Eigen::MatrixXd, Eigen::MatrixXd > &IQ) override
	Генерация блока матрицы, состоящей из интегралов от (r-xi)/|r-xi|^2, влияние профиля самого на себя More...
virtual void	FillMatrixFromOther (const Boundary &otherBoundary, Eigen::MatrixXd &matr) override
	Генерация блока матрицы влияния от другого профиля того же типа More...
virtual void	FillIQFromOther (const Boundary &otherBoundary, std::pair< Eigen::MatrixXd, Eigen::MatrixXd > &IQ) override
	Генерация блока матрицы, состоящей из интегралов от (r-xi)/|r-xi|^2, влияние одного профиля на другой More...
virtual void	SolutionToFreeVortexSheetAndVirtualVortex (const Eigen::VectorXd &sol) override
	Пересчет решения на интенсивность вихревого слоя и на рождаемые вихри на конкретном профиле More...
virtual void	CalcConvVelocityToSetOfPointsFromSheets (const WakeDataBase &pointsDb, std::vector< Point2D > &velo) const override
	Вычисление конвективных скоростей в наборе точек, вызываемых наличием слоев вихрей и источников на профиле More...
virtual void	CalcConvVelocityAtVirtualVortexes (std::vector< Point2D > &velo) const override
	Вычисление конвективной скорости в точках расположения виртуальных вихрей More...
virtual void	ComputeAttachedSheetsIntensity () override
	Вычисление интенсивностей присоединенного вихревого слоя и присоединенного слоя источников More...
virtual void	GetInfluenceFromVorticesToRectPanel (size_t panel, const Vortex2D *ptr, ptrdiff_t count, std::vector< double > &wakeRhs) const override
	Вычисление влияния части подряд идущих вихрей из вихревого следа на прямолинейную панель для правой части More...
virtual void	GetInfluenceFromSourcesToRectPanel (size_t panel, const Vortex2D *ptr, ptrdiff_t count, std::vector< double > &wakeRhs) const override
	Вычисление влияния части подряд источников из области течения на прямолинейную панель для правой части More...
virtual void	GetInfluenceFromVortexSheetAtRectPanelToVortex (size_t panel, const Vortex2D &vtx, Point2D &vel) const override
	Вычисление влияния вихревых слоев (свободный + присоединенный) конкретной прямолинейной панели на вихрь в области течения ///. More...
virtual void	GetInfluenceFromSourceSheetAtRectPanelToVortex (size_t panel, const Vortex2D &vtx, Point2D &vel) const override
	Вычисление влияния слоя источников конкретной прямолинейной панели на вихрь в области течения More...
virtual void	GetInfluenceFromVInfToRectPanel (std::vector< double > &vInfRhs) const override
	Вычисление влияния набегающего потока на прямолинейную панель для правой части More...
size_t	GetUnknownsSize () const
	Возврат размерности вектора решения More...

Public Attributes
const Airfoil &	afl
int	minVortexPerPanel
	Минимальное число вихрей, рождаемых на панели профиля и формирующих виртуальный вихревой след More...
std::vector< std::pair< int, int > >	vortexBeginEnd
	Номера первого и последнего вихрей, рождаемых на каждой панели профиля (формируется после решения СЛАУ) More...
VirtualWake	virtualWake
	Виртуальный вихревой след конкретного профиля More...
size_t	sheetDim
	Размерность параметров каждого из слоев на каждой из панелей More...
Sheet	sheets
	Слои на профиле More...
Sheet	oldSheets
	Слои на профиле с предыдущего шага More...

Protected Attributes
const World2D &	W
	Константная ссылка на решаемую задачу More...
const size_t	numberInPassport
	Номер профиля в паспорте More...

Detailed Description

Класс, определяющий способ удовлетворения граничного условия на обтекаемом профиле

Способ удовлетворения граничного условия:

• вихревой слой представляется кусочно-линейном виде с разрывами в каждом узле профиля;
• условие ортогональности невязки граничного условия константным и линейным проекционным функциям.

Author

Марчевский Илья Константинович

Сокол Ксения Сергеевна

Рятина Евгения Павловна

Колганова Александра Олеговна

\Version 1.12
\date 14 января 2024 г.

Definition at line 65 of file Boundary2DLinLayerAver.h.

Constructor & Destructor Documentation

VM2D::BoundaryLinLayerAver::BoundaryLinLayerAver ( const World2D &  W_, size_t  numberInPassport_  )

inline

Конструктор

Parameters

[in]	W_	константная ссылка на решаемую задачу
[in]	numberInPassport_	номер профиля в паспорте задачи

Definition at line 74 of file Boundary2DLinLayerAver.h.

                                                                          :
            Boundary(W_, numberInPassport_, 2)
        { };

virtual VM2D::BoundaryLinLayerAver::~BoundaryLinLayerAver ( )

inlinevirtual

Деструктор

Definition at line 79 of file Boundary2DLinLayerAver.h.

{};

Here is the call graph for this function:

Member Function Documentation

void BoundaryLinLayerAver::CalcConvVelocityAtVirtualVortexes ( std::vector< Point2D > &  velo ) const

overridevirtual

Вычисление конвективной скорости в точках расположения виртуальных вихрей

Вычисление производится в точках расположения виртуальных вихрей. Скорости находятся как скорости соответствующих точек профиля плюс половина интенсивности присоединенного вихревого слоя, умноженная на касательную к профилю на данной панели.

Parameters

[out]

velo

ссылка на заполняемый список скоростей

Warning

Массив velo заполняется путем присвоения, а не прибавления значений

Implements VM2D::Boundary.

Definition at line 355 of file Boundary2DLinLayerAver.cpp.

{
    std::vector<Point2D>& Vel = velo;

    //Скорости виртуальных вихрей

#pragma omp parallel for default(none) shared(Vel)      
    for (int i = 0; i < (int)Vel.size(); ++i)
        Vel[i] = afl.getV(virtualWake.aflPan[i].second) \
        + virtualWake.vecHalfGamma[i] \
        - W.getPassport().physicalProperties.V0();
}//CalcConvVelocityAtVirtualVortexes(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void BoundaryLinLayerAver::CalcConvVelocityToSetOfPointsFromSheets ( const WakeDataBase &  pointsDb, std::vector< Point2D > &  velo  ) const

overridevirtual

Вычисление конвективных скоростей в наборе точек, вызываемых наличием слоев вихрей и источников на профиле

Вычисляет конвективные скорости в наборе точек, которые вызваны влиянием слоев вихрей и источников на профиле

Parameters

[in]	pointsDb	константная ссылка на базу данных вихрей, в точках которых вычисляются скорости
[out]	velo	ссылка на вектор скоростей, которые приобретают точки из-за влияния слоев вихрей и источников на профиле

Warning

velo — накапливается!

Todo:

Тут надо разобраться, как должно быть...

Todo:

сделать if(move || deform)

Todo:

почему не sheets.freeVortexSheet(j, 0)?

Implements VM2D::Boundary.

Definition at line 268 of file Boundary2DLinLayerAver.cpp.

{
    std::vector<Point2D> selfVelo(pointsDb.vtx.size());

#pragma warning (push)
#pragma warning (disable: 4101)
    //Локальные переменные для цикла
    Point2D velI;
    Point2D tempVel;
#pragma warning (pop)

#pragma omp parallel for default(none) shared(selfVelo, pointsDb, std::cout) private(velI, tempVel)
    for (int i = 0; i < pointsDb.vtx.size(); ++i)
    {
        velI.toZero();

        const Point2D& posI = pointsDb.vtx[i].r();

        for (size_t j = 0; j < afl.getNumberOfPanels(); ++j)
        {
            Point2D dj = afl.getR(j + 1) - afl.getR(j);
            Point2D tauj = dj.unit();

            Point2D s = posI - afl.getR(j);
            Point2D p = posI - afl.getR(j + 1);

            double a = VMlib::Alpha(p, s);

            double lambda = VMlib::Lambda(p, s);

            Point2D u1 = 0.5 / dj.length() * VMlib::Omega(p + s, tauj, tauj) ;

            Point2D skos0 = -a * tauj.kcross() + lambda * tauj;
            Point2D skos1 = -a * u1.kcross() + lambda * u1 - tauj;

            velI += sheets.freeVortexSheet(j, 0) * skos0.kcross() + sheets.freeVortexSheet(j, 1) * skos1.kcross();
            velI += sheets.attachedVortexSheet(j, 0) * skos0.kcross() + sheets.attachedVortexSheet(j, 1) * skos1.kcross();
            velI += sheets.attachedSourceSheet(j, 0) * skos0 + sheets.attachedSourceSheet(j, 1) * skos1;
        }//for j

        velI *= IDPI;
        selfVelo[i] = velI;
    }//for i

    for (size_t i = 0; i < velo.size(); ++i)
        velo[i] += selfVelo[i];
}//CalcConvVelocityToSetOfPointsFromSheets(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void BoundaryLinLayerAver::ComputeAttachedSheetsIntensity ( )

overridevirtual

Вычисление интенсивностей присоединенного вихревого слоя и присоединенного слоя источников

Implements VM2D::Boundary.

Definition at line 370 of file Boundary2DLinLayerAver.cpp.

{
    for (size_t i = 0; i < afl.getNumberOfPanels(); ++i)
    {
        oldSheets.attachedVortexSheet(i, 0) = sheets.attachedVortexSheet(i, 0);
        oldSheets.attachedVortexSheet(i, 1) = sheets.attachedVortexSheet(i, 1);
        oldSheets.attachedSourceSheet(i, 0) = sheets.attachedSourceSheet(i, 0);
        oldSheets.attachedSourceSheet(i, 1) = sheets.attachedSourceSheet(i, 1);
    }

    const Airfoil* oldAfl = (W.getCurrentStep() == 0) ? &afl : &W.getOldAirfoil(numberInPassport);

    for (size_t i = 0; i < afl.getNumberOfPanels(); ++i)
    {
        sheets.attachedVortexSheet(i, 0) = 0.5 * (afl.getV(i + 1) + afl.getV(i)) & afl.tau[i];
        sheets.attachedVortexSheet(i, 1) = (afl.len[i] - oldAfl->len[i]) / W.getPassport().timeDiscretizationProperties.dt;         
            
        sheets.attachedSourceSheet(i, 0) = 0.5 * (afl.getV(i + 1) + afl.getV(i)) & afl.nrm[i];
        sheets.attachedSourceSheet(i, 1) = -afl.len[i] * (afl.phiAfl - oldAfl->phiAfl) / W.getPassport().timeDiscretizationProperties.dt;
                
                 
    }
}//ComputeAttachedSheetsIntensity()

Here is the call graph for this function:

Here is the caller graph for this function:

void BoundaryLinLayerAver::FillIQFromOther ( const Boundary &  otherBoundary, std::pair< Eigen::MatrixXd, Eigen::MatrixXd > &  IQ  )

overridevirtual

Генерация блока матрицы, состоящей из интегралов от (r-xi)/|r-xi|^2, влияние одного профиля на другой

Генерирует блок матрицы влияния от другого профиля того же типа

Parameters

[in]	otherBoundary	константная ссылка на граничное условие на втором профиле
[out]	IQ	ссылка на пару матриц, выражающих взаимные влияния (касательные и нормальные) панелей профиля

Todo:

Пока считается, что граничные условия одинаковые

Implements VM2D::Boundary.

Definition at line 261 of file Boundary2DLinLayerAver.cpp.

{
    afl.calcIQ(2, otherBoundary.afl, IQ);
}//FillIQFromOther(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void BoundaryLinLayerAver::FillIQSelf ( std::pair< Eigen::MatrixXd, Eigen::MatrixXd > &  IQ )

overridevirtual

Генерация блока матрицы, состоящей из интегралов от (r-xi)/|r-xi|^2, влияние профиля самого на себя

Parameters

[out]

IQ

ссылка на генерируемую матрицу

Implements VM2D::Boundary.

Definition at line 235 of file Boundary2DLinLayerAver.cpp.

{
    afl.calcIQ(2, afl, IQ);
}//FillIQSelf(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void BoundaryLinLayerAver::FillMatrixFromOther ( const Boundary &  otherBoundary, Eigen::MatrixXd &  matr  )

overridevirtual

Генерация блока матрицы влияния от другого профиля того же типа

Генерирует блок матрицы влияния от другого профиля того же типа

Parameters

[in]	otherBoundary	константная ссылка на граничное условие на втором профиле
[out]	matr	ссылка на генерируемый блок матрицы

Todo:

Пока считается, что граничные условия одинаковые

Implements VM2D::Boundary.

Definition at line 242 of file Boundary2DLinLayerAver.cpp.

{
    size_t np = afl.getNumberOfPanels();
    size_t npOther = otherBoundary.afl.getNumberOfPanels();

    std::vector<double> res(4, 0.0);

    for (size_t i = 0; i < np; ++i)
    for (size_t j = 0; j < npOther; ++j)
    {
        res = afl.getA(2, i, otherBoundary.afl, j);

        matr(i, j) = res[0];
        matr(i, npOther + j) = res[1];
        matr(np + i, j) = res[2];
        matr(np + i, npOther + j) = res[3];
    }
}//FillMatrixFromOther(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void BoundaryLinLayerAver::FillMatrixSelf ( Eigen::MatrixXd &  matr, Eigen::VectorXd &  lastLine, Eigen::VectorXd &  lactCol  )

overridevirtual

Генерация блока матрицы

Генерирует следующие компоненты матрицы:

• диагональный блок матрицы — влияние данного профиля на самого себя;
• нижнюю строку для матрицы для данного профиля;
• правый столбец матрицы для данного профиля.

Parameters

[out]	matr	ссылка на генерируемую матрицу
[out]	lastLine	ссылка на нижнюю строку
[out]	lactCol	ссылка на правый столбец

Implements VM2D::Boundary.

Definition at line 209 of file Boundary2DLinLayerAver.cpp.

{
    size_t np = afl.getNumberOfPanels();

    std::vector<double> res(4, 0.0);

    for (size_t i = 0; i < np; ++i)
    {
        lastCol(i) = 1.0;
        lastLine(i) = afl.len[i];
        lastCol(np + i) = 0.0;
        lastLine(np + i) = 0.0;
    }

    for (size_t i = 0; i < np; ++i)
    for (size_t j = 0; j < np; ++j)
    {
        res = afl.getA(2, i, afl, j);
        matr(i, j) = res[0];
        matr(i, np + j) = res[1];
        matr(np + i, j) = res[2];
        matr(np + i, np + j) = res[3];
    }
}//FillMatrixSelf(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void BoundaryLinLayerAver::GetInfluenceFromSourceSheetAtRectPanelToVortex ( size_t  panel, const Vortex2D &  vtx, Point2D &  vel  ) const

overridevirtual

Вычисление влияния слоя источников конкретной прямолинейной панели на вихрь в области течения

Parameters

[in]	panel	номер панели профиля, от которой считается влияние
[in]	vtx	ссылка на вихрь
[out]	vel	ссылка на вектор полученной скорости

Implements VM2D::Boundary.

Definition at line 459 of file Boundary2DLinLayerAver.cpp.

{
    vel.toZero();

    const Point2D& posI = ptr.r();

    Point2D dj = afl.getR(panel + 1) - afl.getR(panel);
    Point2D tauj = dj.unit();

    Point2D s = posI - afl.getR(panel);
    Point2D p = posI - afl.getR(panel + 1);

    Point2D u1 = 0.5 / dj.length()* VMlib::Omega(p + s, tauj, tauj);

    double a = VMlib::Alpha(p, s);

    double lambda;
    if ((s.length2() > 1e-16) && (p.length2() > 1e-16))
        lambda = VMlib::Lambda(p, s);
    else
        lambda = 0.0;

    vel += sheets.attachedSourceSheet(panel, 0) * (-a * u1.kcross() + lambda * u1 - tauj);
    vel *= IDPI;
}// GetInfluenceFromSourceSheetAtRectPanelToVortex(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void BoundaryLinLayerAver::GetInfluenceFromSourcesToRectPanel ( size_t  panel, const Vortex2D *  ptr, ptrdiff_t  count, std::vector< double > &  wakeRhs  ) const

overridevirtual

Вычисление влияния части подряд источников из области течения на прямолинейную панель для правой части

Вычисляет влияния части подряд идущих источников из области течения на прямолинейную панель для правой части

Parameters

[in]	panel	номер панели профиля, на которую считается влияние
[in]	ptr	указатель на начало диапазона источников
[in]	count	длина диапазона источников
[out]	wakeRhs	ссылка на вектор полученных влияние для правой части СЛАУ

Implements VM2D::Boundary.

Definition at line 427 of file Boundary2DLinLayerAver.cpp.

{
    double& velI = wakeRhs[0];
    double& velILin = wakeRhs[1];

    const Point2D& posI0 = afl.getR(panel);
    const Point2D& posI1 = afl.getR(panel + 1);
    Point2D di = posI1 - posI0;
    const Point2D& taui = afl.tau[panel];

    for (size_t it = 0; it != count; ++it)
    {
        const Vortex2D& vt = ptr[it];
        const Point2D& posJ = vt.r();
        const double& gamJ = vt.g();

        Point2D s = posJ - posI0;
        Point2D p = posJ - posI1;

        Point2D u1 = 0.5 / di.length() * VMlib::Omega(p + s, taui, taui);

        double alpha = VMlib::Alpha(p, s);
        double lambda = VMlib::Lambda(p, s);

        velI -= gamJ * lambda;

        velILin -= gamJ * (alpha * (u1 &  taui.kcross()) + lambda * (u1 & taui) - 1.0);
    }

}//GetInfluenceFromSourcesToRectPanel(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void BoundaryLinLayerAver::GetInfluenceFromVInfToRectPanel ( std::vector< double > &  vInfRhs ) const

overridevirtual

Вычисление влияния набегающего потока на прямолинейную панель для правой части

Вычисляет влияния набегающего потока на прямолинейную панель для правой части

Parameters

[out]

vInfRhs

ссылка на вектор полученных влияние для правой части СЛАУ

Implements VM2D::Boundary.

Definition at line 512 of file Boundary2DLinLayerAver.cpp.

{
    size_t np = afl.getNumberOfPanels();

    vInfRhs.resize(2 * np);

#pragma omp parallel for default(none) shared(vInfRhs, np)
    for (int i = 0; i < np; ++i)
    {
        vInfRhs[i] = afl.tau[i] & W.getPassport().physicalProperties.V0();
        vInfRhs[np + i] = 0;
    }
}

void BoundaryLinLayerAver::GetInfluenceFromVortexSheetAtRectPanelToVortex ( size_t  panel, const Vortex2D &  vtx, Point2D &  vel  ) const

overridevirtual

Вычисление влияния вихревых слоев (свободный + присоединенный) конкретной прямолинейной панели на вихрь в области течения ///.

Parameters

[in]	panel	номер панели профиля, от которой считается влияние
[in]	vtx	ссылка на вихрь
[out]	vel	ссылка на вектор полученной скорости

Implements VM2D::Boundary.

Definition at line 486 of file Boundary2DLinLayerAver.cpp.

{
    vel.toZero();

    const Point2D& posI = ptr.r();

    Point2D dj = afl.getR(panel + 1) - afl.getR(panel);
    Point2D tauj = dj.unit();

    Point2D s = posI - afl.getR(panel);
    Point2D p = posI - afl.getR(panel + 1);

    Point2D u1 = 0.5 / dj.length()* VMlib::Omega(p + s, tauj, tauj);

    double lambda;
    if ((s.length2() > 1e-16) && (p.length2() > 1e-16))
        lambda = VMlib::Lambda(p, s);
    else
        lambda = 0.0;

    vel += sheets.freeVortexSheet(panel, 0) * (lambda * u1.kcross() - tauj.kcross());
    vel += sheets.attachedVortexSheet(panel, 0) * (lambda * u1.kcross() - tauj.kcross());
    vel *= IDPI;
}// GetInfluenceFromVortexSheetAtRectPanelToVortex(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void BoundaryLinLayerAver::GetInfluenceFromVorticesToRectPanel ( size_t  panel, const Vortex2D *  ptr, ptrdiff_t  count, std::vector< double > &  wakeRhs  ) const

overridevirtual

Вычисление влияния части подряд идущих вихрей из вихревого следа на прямолинейную панель для правой части

Вычисляет влияния части подряд идущих вихрей из вихревого следа на прямолинейную панель для правой части

Parameters

[in]	panel	номер панели профиля, на которую считается влияние
[in]	ptr	указатель на начало диапазона вихрей
[in]	count	длина диапазона вихрей
[out]	wakeRhs	ссылка на вектор полученных влияние для правой части СЛАУ

Implements VM2D::Boundary.

Definition at line 395 of file Boundary2DLinLayerAver.cpp.

{
    double& velI = wakeRhs[0];
    double& velILin = wakeRhs[1];

    const Point2D& posI0 = afl.getR(panel);
    const Point2D& posI1 = afl.getR(panel + 1);
    Point2D di = posI1 - posI0;
    const Point2D& taui = afl.tau[panel];

    for (size_t it = 0; it != count; ++it)
    {
        const Vortex2D& vt = ptr[it];
        const Point2D& posJ = vt.r();
        const double& gamJ = vt.g();

        Point2D s = posJ - posI0;
        Point2D p = posJ - posI1;

        Point2D u1 = 0.5 / di.length() * VMlib::Omega(p + s, taui, taui);

        double alpha = VMlib::Alpha(p, s);
        double lambda = VMlib::Lambda(p, s);

        velI -= gamJ * alpha;

        velILin -= gamJ * (alpha * (u1 & taui) + lambda * (u1 & (-taui.kcross())));
    }
}//GetInfluenceFromVorticesToRectPanel(...)

Here is the call graph for this function:

Here is the caller graph for this function:

size_t Boundary::GetUnknownsSize ( ) const

inherited

Возврат размерности вектора решения

(без учета регуляризирующей переменной)

Returns

размерность вектора решения

Definition at line 71 of file Boundary2D.cpp.

{
    return sheetDim * afl.getNumberOfPanels();
}

void BoundaryLinLayerAver::SolutionToFreeVortexSheetAndVirtualVortex ( const Eigen::VectorXd &  sol )

overridevirtual

Пересчет решения на интенсивность вихревого слоя и на рождаемые вихри на конкретном профиле

1) Приводит решение к интенсивности вихревого слоя и записывает его в sheets.freeVortexSheet:

• если неизвестное — интенсивность вихря, то он "размазывается" по панели;
• если неизвестное — интенсивность слоя, то она передается непосредственно.

2) Приводит интенсивность вихревого слоя к рождаемым вихрям, а также вычисляет их положения

Parameters

[in]

sol

вектор решения СЛАУ

Implements VM2D::Boundary.

Definition at line 57 of file Boundary2DLinLayerAver.cpp.

{
    Vortex2D virtVort;
    Point2D midNorm;

    double delta = W.getPassport().wakeDiscretizationProperties.delta;

    int nVortPerPan = W.getPassport().wakeDiscretizationProperties.minVortexPerPanel;

    //Очистка и резервирование памяти
    virtualWake.vecHalfGamma.clear();
    virtualWake.vecHalfGamma.reserve(afl.getNumberOfPanels() * nVortPerPan);

    //Очистка и резервирование памяти
    virtualWake.aflPan.clear();
    virtualWake.aflPan.reserve(afl.getNumberOfPanels() * nVortPerPan);

    //Очистка и резервирование памяти
    virtualWake.vtx.clear();
    virtualWake.vtx.reserve(afl.getNumberOfPanels() * nVortPerPan);

    //Очистка и резервирование памяти
    vortexBeginEnd.clear();
    vortexBeginEnd.reserve(afl.getNumberOfPanels());

    double maxG = W.getPassport().wakeDiscretizationProperties.maxGamma;

    std::pair<int, int> pair;

    for (size_t i = 0; i < afl.getNumberOfPanels(); ++i)
    {
        midNorm = afl.nrm[i] * delta;

        //число участков, на которые разбивается панель для сброса вихрей
        //определяется через наибольшее значение решения на профиле, т.е. в крайней точке
         

        pair.first = (int)virtualWake.vtx.size();

        double a = sol(i) - 0.5 * sol(afl.getNumberOfPanels() + i);
        double b = sol(i) + 0.5 * sol(afl.getNumberOfPanels() + i);

        if (fabs(a - b) < 1e-10)
        {
            size_t NEWnVortPerPan = (size_t)std::max((int)std::ceil(fabs(sol(i)*afl.len[i]) / maxG), nVortPerPan);
            Point2D dr = 1.0 / NEWnVortPerPan * (afl.getR(i + 1) - afl.getR(i));

            for (size_t j = 0; j < NEWnVortPerPan; ++j)
            {
                virtVort.r() = afl.getR(i) + dr * (j * 1.0 + 0.5) + midNorm;
                virtVort.g() = sol(i) * afl.len[i] / NEWnVortPerPan;
                virtualWake.vtx.push_back(virtVort);

                virtualWake.vecHalfGamma.push_back(0.5 * sol(i)  * afl.tau[i]);
                virtualWake.aflPan.push_back({ numberInPassport, i });
            }
        }
        else if (a * b >= 0.0)
        {
            size_t NEWnVortPerPan = (size_t)std::max((int)std::ceil(fabs(0.5*(a+b)*afl.len[i]) / maxG), nVortPerPan);
            std::vector<double> ds(NEWnVortPerPan+1, 0.0);
            
            for (size_t k = 1; k <= NEWnVortPerPan; ++k)
            {
                if ((a > 0) || (b > 0))
                    ds[k] = afl.len[i] / (a - b) * (a - sqrt((b*b * k + a * a*(NEWnVortPerPan - k)) / NEWnVortPerPan));
                else
                    ds[k] = afl.len[i] / (a - b) * (a + sqrt((b*b * k + a * a*(NEWnVortPerPan - k)) / NEWnVortPerPan));
            }

            for (size_t j = 0; j < NEWnVortPerPan; ++j)
            {
                virtVort.r() = afl.getR(i) + 0.5*(ds[j]+ds[j+1])*afl.tau[i] + midNorm;
                virtVort.g() = 0.5 * (a+b) * afl.len[i] / NEWnVortPerPan;
                virtualWake.vtx.push_back(virtVort);

                virtualWake.vecHalfGamma.push_back(0.5 * (a + 0.5*(ds[j] + ds[j + 1]) * (b-a)/ afl.len[i])  * afl.tau[i]);
                virtualWake.aflPan.push_back({ numberInPassport, i });
            }
        }
        else
        {
            double sast = -a * afl.len[i] / (b - a);

            //from 0 to sast
            {
                size_t NEWnVortPerPan = (size_t)std::max((int)std::ceil(fabs(0.5*(a + 0.0)*sast) / maxG), (int)std::ceil(nVortPerPan * sast / afl.len[i]));
                std::vector<double> ds(NEWnVortPerPan + 1, 0.0);

                for (size_t k = 1; k <= NEWnVortPerPan; ++k)
                {
                    if (a > 0)
                        ds[k] = sast / (a - 0.0) * (a - sqrt((a * a*(NEWnVortPerPan - k)) / NEWnVortPerPan));
                    else
                        ds[k] = sast / (a - 0.0) * (a + sqrt((a * a*(NEWnVortPerPan - k)) / NEWnVortPerPan));
                }

                for (size_t j = 0; j < NEWnVortPerPan; ++j)
                {
                    virtVort.r() = afl.getR(i) + 0.5*(ds[j] + ds[j + 1])*afl.tau[i] + midNorm;
                    virtVort.g() = 0.5 * (a + 0.0) * sast / NEWnVortPerPan;
                    virtualWake.vtx.push_back(virtVort);

                    virtualWake.vecHalfGamma.push_back(0.5 * (a + 0.5*(ds[j] + ds[j + 1]) * (0.0 - a) / sast)  * afl.tau[i]);
                    virtualWake.aflPan.push_back({ numberInPassport, i });
                }
            }

            double sastast = afl.len[i] - sast;

            //from sats to len
            {
                size_t NEWnVortPerPan = (size_t)std::max((int)std::ceil(fabs(0.5*(0.0 + b)*sastast) / maxG), (int)std::ceil(nVortPerPan * sastast / afl.len[i]));
                std::vector<double> ds(NEWnVortPerPan + 1, 0.0);

                for (size_t k = 1; k <= NEWnVortPerPan; ++k)
                {
                    if (b > 0)
                        ds[k] = sastast / (0.0 - b) * (0.0 - sqrt((b*b * k) / NEWnVortPerPan));
                    else
                        ds[k] = sastast / (0.0 - b) * (0.0 + sqrt((b*b * k) / NEWnVortPerPan));
                }

                for (size_t j = 0; j < NEWnVortPerPan; ++j)
                {
                    virtVort.r() = afl.getR(i) + sast * afl.tau[i] + 0.5*(ds[j] + ds[j + 1])*afl.tau[i] + midNorm;
                    virtVort.g() = 0.5 * (0.0 + b) * sastast / NEWnVortPerPan;
                    virtualWake.vtx.push_back(virtVort);

                    virtualWake.vecHalfGamma.push_back(0.5 * (0.0 + 0.5*(ds[j] + ds[j + 1]) * (b - 0.0) / sastast)  * afl.tau[i]);
                    virtualWake.aflPan.push_back({ numberInPassport, i });
                }

            }
        }

        pair.second = (int)virtualWake.vtx.size();
        vortexBeginEnd.push_back(pair);
    }


    for (size_t j = 0; j < afl.getNumberOfPanels(); ++j)
    {
        sheets.freeVortexSheet(j, 0) = sol(j);
        sheets.freeVortexSheet(j, 1) = sol(afl.getNumberOfPanels() + j);
    }

}//SolutionToFreeVortexSheetAndVirtualVortex(...)

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

const Airfoil& VM2D::Boundary::afl

inherited

Константная ссылка на профиль
инициализируется автоматом в конструкторе

Definition at line 76 of file Boundary2D.h.

int VM2D::Boundary::minVortexPerPanel

inherited

Минимальное число вихрей, рождаемых на панели профиля и формирующих виртуальный вихревой след

Definition at line 79 of file Boundary2D.h.

const size_t VM2D::Boundary::numberInPassport

protectedinherited

Номер профиля в паспорте

Definition at line 70 of file Boundary2D.h.

Sheet VM2D::Boundary::oldSheets

inherited

Слои на профиле с предыдущего шага

Definition at line 98 of file Boundary2D.h.

size_t VM2D::Boundary::sheetDim

inherited

Размерность параметров каждого из слоев на каждой из панелей

Указывает, сколькими числами задается интенсивность каждого из слоев на каждой панели:

• 1 — одно число — задается только среднее значение;
• 2 — два числа — задается среднее значение и "наклон".

Definition at line 92 of file Boundary2D.h.

Sheet VM2D::Boundary::sheets

inherited

Слои на профиле

Definition at line 95 of file Boundary2D.h.

VirtualWake VM2D::Boundary::virtualWake

inherited

Виртуальный вихревой след конкретного профиля

Definition at line 85 of file Boundary2D.h.

std::vector<std::pair<int, int> > VM2D::Boundary::vortexBeginEnd

inherited

Номера первого и последнего вихрей, рождаемых на каждой панели профиля (формируется после решения СЛАУ)

Definition at line 82 of file Boundary2D.h.

const World2D& VM2D::Boundary::W

protectedinherited

Константная ссылка на решаемую задачу

Definition at line 67 of file Boundary2D.h.

---

The documentation for this class was generated from the following files:

• VM2D/Boundary2D/Boundary2DLinLayerAver.h
• VM2D/Boundary2D/Boundary2DLinLayerAver.cpp

• VM2D
• BoundaryLinLayerAver
• Generated by 1.8.11