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

#include <Boundary2DLinLayerAver.h>

Inheritance diagram for VM2D::BoundaryLinLayerAver:

Collaboration diagram for VM2D::BoundaryLinLayerAver:

[legend]

Public Member Functions
	BoundaryLinLayerAver (const World2D &W_, size_t numberInPassport_)
	Конструктор

virtual	~BoundaryLinLayerAver ()
	Деструктор

virtual void	FillMatrixSelf (Eigen::MatrixXd &matr, Eigen::VectorXd &lastLine, Eigen::VectorXd &lactCol) override
	Генерация блока матрицы

virtual void	FillIQSelf (std::pair< Eigen::MatrixXd, Eigen::MatrixXd > &IQ) override
	Генерация блока матрицы, состоящей из интегралов от (r-xi)/\|r-xi\|^2, влияние профиля самого на себя

virtual void	FillMatrixFromOther (const Boundary &otherBoundary, Eigen::MatrixXd &matr) override
	Генерация блока матрицы влияния от другого профиля того же типа

virtual void	FillIQFromOther (const Boundary &otherBoundary, std::pair< Eigen::MatrixXd, Eigen::MatrixXd > &IQ) override
	Генерация блока матрицы, состоящей из интегралов от (r-xi)/\|r-xi\|^2, влияние одного профиля на другой

virtual void	SolutionToFreeVortexSheetAndVirtualVortex (const Eigen::VectorXd &sol) override
	Пересчет решения на интенсивность вихревого слоя и на рождаемые вихри на конкретном профиле

virtual void	CalcConvVelocityToSetOfPointsFromSheets (const WakeDataBase &pointsDb, std::vector< Point2D > &velo) const override
	Вычисление конвективных скоростей в наборе точек, вызываемых наличием слоев вихрей и источников на профиле

virtual void	ComputeAttachedSheetsIntensity () override
	Вычисление интенсивностей присоединенного вихревого слоя и присоединенного слоя источников

virtual void	GetInfluenceFromVorticesToRectPanel (size_t panel, const Vortex2D *ptr, ptrdiff_t count, std::vector< double > &wakeRhs) const override
	Вычисление влияния части подряд идущих вихрей из вихревого следа на прямолинейную панель для правой части

virtual void	GetInfluenceFromSourcesToRectPanel (size_t panel, const Vortex2D *ptr, ptrdiff_t count, std::vector< double > &wakeRhs) const override
	Вычисление влияния части подряд источников из области течения на прямолинейную панель для правой части

virtual void	GetInfluenceFromVortexSheetAtRectPanelToVortex (size_t panel, const Vortex2D &vtx, Point2D &vel) const override
	Вычисление влияния вихревых слоев (свободный + присоединенный) конкретной прямолинейной панели на вихрь в области течения ///.

virtual void	GetInfluenceFromSourceSheetAtRectPanelToVortex (size_t panel, const Vortex2D &vtx, Point2D &vel) const override
	Вычисление влияния слоя источников конкретной прямолинейной панели на вихрь в области течения

virtual void	GetInfluenceFromVInfToRectPanel (std::vector< double > &vInfRhs) const override
	Вычисление влияния набегающего потока на прямолинейную панель для правой части

size_t	GetUnknownsSize () const
	Возврат размерности вектора решения

Public Attributes
const Airfoil &	afl

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

std::vector< std::pair< int, int > >	vortexBeginEnd
	Номера первого и последнего вихрей, рождаемых на каждой панели профиля (формируется после решения СЛАУ)

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

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

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

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

Protected Attributes
const World2D &	W
	Константная ссылка на решаемую задачу

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

Detailed Description

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

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

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

Author: Марчевский Илья Константинович; Сокол Ксения Сергеевна; Рятина Евгения Павловна; Колганова Александра Олеговна

\Version 1.14

Date: 6 марта 2026 г.

Definition at line 65 of file Boundary2DLinLayerAver.h.

Constructor & Destructor Documentation

◆ BoundaryLinLayerAver()

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)
        { };

◆ ~BoundaryLinLayerAver()

virtual VM2D::BoundaryLinLayerAver::~BoundaryLinLayerAver ( )

inlinevirtual

Деструктор

Definition at line 79 of file Boundary2DLinLayerAver.h.

79{};

Member Function Documentation

◆ CalcConvVelocityToSetOfPointsFromSheets()

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

overridevirtual

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

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

Parameters

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

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

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

Todo:

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

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

Implements VM2D::Boundary.

Definition at line 269 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, IDPI) 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:

◆ ComputeAttachedSheetsIntensity()

void BoundaryLinLayerAver::ComputeAttachedSheetsIntensity ( )

overridevirtual

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

Implements VM2D::Boundary.

Definition at line 373 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 AirfoilGeometry* 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:

◆ FillIQFromOther()

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 262 of file Boundary2DLinLayerAver.cpp.

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

Here is the call graph for this function:

◆ FillIQSelf()

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 236 of file Boundary2DLinLayerAver.cpp.

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

Here is the call graph for this function:

◆ FillMatrixFromOther()

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

overridevirtual

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

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

Parameters

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

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

Implements VM2D::Boundary.

Definition at line 243 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:

◆ FillMatrixSelf()

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 208 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:

◆ GetInfluenceFromSourceSheetAtRectPanelToVortex()

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 462 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:

◆ GetInfluenceFromSourcesToRectPanel()

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 430 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:

◆ GetInfluenceFromVInfToRectPanel()

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

overridevirtual

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

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

Parameters

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

Implements VM2D::Boundary.

Definition at line 515 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.getV0();
        vInfRhs[np + i] = 0;
    }
}

Here is the call graph for this function:

◆ GetInfluenceFromVortexSheetAtRectPanelToVortex()

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 489 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:

◆ GetInfluenceFromVorticesToRectPanel()

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 398 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:

◆ GetUnknownsSize()

size_t Boundary::GetUnknownsSize ( ) const

inherited

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

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

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

Definition at line 73 of file Boundary2D.cpp.

{
    return sheetDim * afl.getNumberOfPanels();
}

Here is the call graph for this function:

Here is the caller graph for this function:

◆ SolutionToFreeVortexSheetAndVirtualVortex()

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

overridevirtual

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

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

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

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

Parameters

[in] sol вектор решения СЛАУ

Implements VM2D::Boundary.

Definition at line 56 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:

Member Data Documentation

◆ afl

const Airfoil& VM2D::Boundary::afl

inherited

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

Definition at line 77 of file Boundary2D.h.

◆ minVortexPerPanel

int VM2D::Boundary::minVortexPerPanel

inherited

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

Definition at line 80 of file Boundary2D.h.

◆ numberInPassport

const size_t VM2D::Boundary::numberInPassport

protectedinherited

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

Definition at line 71 of file Boundary2D.h.

◆ oldSheets

Sheet VM2D::Boundary::oldSheets

inherited

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

Definition at line 99 of file Boundary2D.h.

◆ sheetDim

size_t VM2D::Boundary::sheetDim

inherited

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

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

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

Definition at line 93 of file Boundary2D.h.

◆ sheets

Sheet VM2D::Boundary::sheets

inherited

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

Definition at line 96 of file Boundary2D.h.

◆ virtualWake

VirtualWake VM2D::Boundary::virtualWake

inherited

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

Definition at line 86 of file Boundary2D.h.

◆ vortexBeginEnd

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

inherited

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

Definition at line 83 of file Boundary2D.h.

◆ W

const World2D& VM2D::Boundary::W

protectedinherited

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

Definition at line 68 of file Boundary2D.h.

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

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

Public Member Functions

Public Attributes

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ BoundaryLinLayerAver()

◆ ~BoundaryLinLayerAver()

Member Function Documentation

◆ CalcConvVelocityToSetOfPointsFromSheets()

◆ ComputeAttachedSheetsIntensity()

◆ FillIQFromOther()

◆ FillIQSelf()

◆ FillMatrixFromOther()

◆ FillMatrixSelf()

◆ GetInfluenceFromSourceSheetAtRectPanelToVortex()

◆ GetInfluenceFromSourcesToRectPanel()

◆ GetInfluenceFromVInfToRectPanel()

◆ GetInfluenceFromVortexSheetAtRectPanelToVortex()

◆ GetInfluenceFromVorticesToRectPanel()

◆ GetUnknownsSize()

◆ SolutionToFreeVortexSheetAndVirtualVortex()

Member Data Documentation

◆ afl

◆ minVortexPerPanel

◆ numberInPassport

◆ oldSheets

◆ sheetDim

◆ sheets

◆ virtualWake

◆ vortexBeginEnd

◆ W