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

#include <Airfoil2DRect.h>

Inheritance diagram for VM2D::AirfoilRect:

Collaboration diagram for VM2D::AirfoilRect:

[legend]

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

	AirfoilRect (const Airfoil &afl)

virtual	~AirfoilRect ()
	Деструктор More...

void	CalcNrmTauLen ()
	Вычисление нормалей, касательных и длин панелей по текущему положению вершин More...

virtual void	GetGabarits (double gap=0.02) override
	Вычисляет габаритный прямоугольник профиля More...

virtual void	ReadFromFile (const std::string &dir) override
	Считывание профиля из файла More...

virtual void	Rotate (double alpha) override
	Поворот профиля More...

virtual void	Scale (const Point2D &) override
	Масштабирование профиля More...

virtual void	Move (const Point2D &dr) override
	Перемещение профиля More...

virtual std::vector< double >	getA (size_t p, size_t i, const Airfoil &airfoil, size_t j) const override
	Вычисление коэффициентов матрицы A для расчета влияния панели на панель More...

virtual void	calcIQ (size_t p, const Airfoil &otherAirfoil, std::pair< Eigen::MatrixXd, Eigen::MatrixXd > &matrPair) const override
	Вычисление коэффициентов матрицы, состоящей из интегралов от (r-xi)/\|r-xi\|^2. More...

virtual bool	IsPointInAirfoil (const Point2D &point) const override
	Определяет, находится ли точка с радиус-вектором \( \vec r \) внутри профиля More...

virtual void	GetDiffVelocityI0I3ToSetOfPointsAndViscousStresses (const WakeDataBase &pointsDb, std::vector< double > &domainRadius, std::vector< double > &I0, std::vector< Point2D > &I3) override
	Вычисление числителей и знаменателей диффузионных скоростей в заданном наборе точек, обусловленных геометрией профиля, и вычисление вязкого трения More...

virtual void	GetDiffVelocityI0I3ToWakeAndViscousStresses (const WakeDataBase &pointsDb, std::vector< double > &domainRadius, std::vector< double > &I0, std::vector< Point2D > &I3) override

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

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

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

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

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

const Point2D &	getR (size_t q) const
	Возврат константной ссылки на вершину профиля More...

const Point2D &	getV (size_t q) const
	Возврат константной ссылки на скорость вершины профиля More...

void	setV (const Point2D &vel)
	Установка постоянной скорости всех вершин профиля More...

void	setV (const std::vector< Point2D > &vel)
	Установка скоростей всех вершин профиля More...

size_t	getNumberOfPanels () const
	Возврат количества панелей на профиле More...

	IFCUDA (mutable double *devRPtr)
	Указатель на девайсе, где хранятся вершины профиля More...

	IFCUDA (mutable double *devRhsPtr)
	Указатель на девайсе, где хранится правая часть (константная) матрицы More...

	IFCUDA (mutable double *devRhsLinPtr)
	Указатель на девайсе, где хранится правая часть (линейная) матрицы More...

	IFCUDA (mutable std::vector< double > tmpRhs)
	Указатель на хосте, где хранится временная часть матрицы, полученная с девайса More...

	IFCUDA (mutable double *devFreeVortexSheetPtr)
	Указатель на девайсе, где хранятся интенсивности свободного вихревого слоя на панелях More...

	IFCUDA (mutable double *devAttachedVortexSheetPtr)
	Указатель на девайсе, где хранятся интенсивности присоединенного вихревого слоя на панелях More...

	IFCUDA (mutable double *devAttachedSourceSheetPtr)
	Указатель на девайсе, где хранятся интенсивности присоединенного слоя источников на панелях More...

	IFCUDA (mutable double *devMeanEpsOverPanelPtr)
	Указатель на девайсе, где хранятся средние eps на панелях More...

	IFCUDA (mutable double *devViscousStressesPtr)
	Указатель на девайсе, где хранится вектор (по панелям) для силы вязкого трения More...

	IFCUDA (mutable std::vector< double > tmpViscousStresses)
	Указатель на хосте, где хранится временная часть вектора (по панелям) для силы вязкого трения More...

bool	isAfter (size_t i, size_t j) const
	Проверка, идет ли вершина i следом за вершиной j. More...

bool	isInsideGabarits (const Point2D &r) const
	Определяет, находится ли точка с радиус-вектором \( \vec r \) внутри габаритного прямоугольника профиля More...

bool	isOutsideGabarits (const Point2D &r) const
	Определяет, находится ли точка с радиус-вектором \( \vec r \) вне габаритного прямоугольника профиля More...

void	calcMeanEpsOverPanel ()
	Вычисление средних значений eps на панелях More...

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

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

Point2D	rcm
	Положение центра масс профиля More...

double	phiAfl
	Поворот профиля More...

double	area
	Площадь профиля More...

double	m
	Масса профиля More...

double	J
	Полярный момент инерции профиля относительно центра масс More...

std::vector< double >	meanEpsOverPanel
	Средние значения Eps на панелях More...

bool	inverse
	Признак разворота нормалей (для расчета внутренних течений) More...

std::vector< Point2D >	nrm
	Нормали к панелям профиля More...

std::vector< Point2D >	tau
	Касательные к панелям профиля More...

std::vector< double >	len
	Длины панелей профиля More...

std::vector< double >	viscousStress
	Касательные напряжения на панелях профиля More...

Point2D	lowLeft
	Левый нижний угол габаритного прямоугольника профиля More...

Point2D	upRight
	Правый верхний угол габаритного прямоугольника профиля More...

std::vector< double >	gammaThrough
	Суммарные циркуляции вихрей, пересекших панели профиля на прошлом шаге More...

Protected Attributes
std::vector< Point2D >	r_
	Координаты начал панелей More...

std::vector< Point2D >	v_
	Скорости начал панелей More...

Detailed Description

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

Тип профиля:

профиль с прямолинейными панелями.

Author

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

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

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

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

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

Definition at line 64 of file Airfoil2DRect.h.

Constructor & Destructor Documentation

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

inline

Конструктор

Definition at line 70 of file Airfoil2DRect.h.

71 :Airfoil(W_, numberInPassport_)

72 { };

VM2D::Airfoil::Airfoil

Airfoil(const World2D &W_, const size_t numberInPassport_)

Definition: Airfoil2D.cpp:55

VM2D::AirfoilRect::AirfoilRect ( const Airfoil & afl )

inline

Definition at line 74 of file Airfoil2DRect.h.

74 : Airfoil(afl) {};

VM2D::Airfoil::Airfoil

Airfoil(const World2D &W_, const size_t numberInPassport_)

Definition: Airfoil2D.cpp:55

virtual VM2D::AirfoilRect::~AirfoilRect ( )

inlinevirtual

Деструктор

Definition at line 77 of file Airfoil2DRect.h.

77 { };

Here is the call graph for this function:

Member Function Documentation

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

overridevirtual

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

Parameters

[in]	p	размерность матрицы-результата
[in]	otherAirfoil	константная ссылка на профиль, от которого рассчитывается влияние
[out]	matrPair	ссылка на пару матриц, выражающих взаимные влияния (касательные и нормальные) панелей профиля return соответствующий блок матрицы СЛАУ, вытянутый в линию

Implements VM2D::Airfoil.

Definition at line 708 of file Airfoil2DRect.cpp.

 {
     bool self = (&otherAirfoil == this);
 
     size_t aflNSelf = numberInPassport;
     size_t aflNOther = otherAirfoil.numberInPassport;
 
     size_t npI;
     size_t npJ;
 
     numvector<double, 3> alpha, lambda;
 
     //auxillary vectors
     Point2D p1, s1, p2, s2, di, dj, i00, i01, i10, i11;
     numvector<Point2D, 3> v00, v11;
     numvector<Point2D, 2> v01, v10;
 
 #pragma omp parallel for \
     default(none) \
     shared(otherAirfoil, self, aflNSelf, aflNOther, matrPair, p, IDPI, IQPI) \
     private(npI, npJ, alpha, lambda, p1, s1, p2, s2, di, dj, i00, i01, i10, i11, v00, v11, v01, v10) schedule(dynamic, DYN_SCHEDULE)
     for(int i = 0; i < (int)getNumberOfPanels(); ++i)
         for (int j = 0; j < (int)otherAirfoil.getNumberOfPanels(); ++j)
         {
             npI = getNumberOfPanels();
             npJ = otherAirfoil.getNumberOfPanels();
 
             if ((i == j) && self)
             {
 
                 matrPair.first(i, j) = 0.0;
                 matrPair.second(i, j) = 0.0;
 
                 if (p == 2)
                 {
                     matrPair.first(i, npJ + j) = 0.0;
                     matrPair.first(npI + i, j) = 0.0;
 
                     matrPair.second(i, npJ + j) = -IQPI;
                     matrPair.second(npI + i, j) = IQPI;
 
                     matrPair.first(npI + i, npJ + j) = 0.0;
                     matrPair.second(npI + i, npJ + j) = 0.0;
 
                 }
                 if (p > 2)
                     throw (-42);
             }//if i==j
             else
             {
                 const Point2D& taui = tau[i];
                 const Point2D& tauj = otherAirfoil.tau[j];
 
                 p1 = getR(i + 1) - otherAirfoil.getR(j + 1);
                 s1 = getR(i + 1) - otherAirfoil.getR(j);
                 p2 = getR(i) - otherAirfoil.getR(j + 1);
                 s2 = getR(i) - otherAirfoil.getR(j);
                 di = getR(i + 1) - getR(i);
                 dj = otherAirfoil.getR(j + 1) - otherAirfoil.getR(j);
 
                 alpha = { \
                     (self && isAfter(j, i)) ? 0.0 : VMlib::Alpha(s2, s1), \
 					VMlib::Alpha(s2, p1), \
                     (self && isAfter(i, j)) ? 0.0 : VMlib::Alpha(p1, p2) \
                 };
 
                 lambda = { \
                     (self && isAfter(j, i)) ? 0.0 : VMlib::Lambda(s2, s1), \
 					VMlib::Lambda(s2, p1), \
                     (self && isAfter(i, j)) ? 0.0 : VMlib::Lambda(p1, p2) \
                 };
 
                 v00 = {
                     VMlib::Omega(s1, taui, tauj),
                     -VMlib::Omega(di, taui, tauj),
                     VMlib::Omega(p2, taui, tauj)
                 };
 
                 i00 = IDPI / len[i] * (-(alpha[0] * v00[0] + alpha[1] * v00[1] + alpha[2] * v00[2]).kcross() \
                     + (lambda[0] * v00[0] + lambda[1] * v00[1] + lambda[2] * v00[2]));
 
                 matrPair.first(i, j) = i00 & nrm[i];
                 matrPair.second(i, j) = i00 & taui;
 
 
                 if (p == 2)
                 {
                     v01 = {
                         0.5 / (dj.length()) * (((p1 + s1) & tauj) * VMlib::Omega(s1, taui, tauj) - s1.length2() * taui),
                         -0.5 * di.length() / dj.length() * VMlib::Omega(s1 + p2, tauj, tauj)
                     };
 
                     i01 = IDPI / len[i] * (-((alpha[0] + alpha[2]) * v01[0] + (alpha[1] + alpha[2]) * v01[1]).kcross() \
                         + ((lambda[0] + lambda[2]) * v01[0] + (lambda[1] + lambda[2]) * v01[1]) - 0.5 * di.length() * tauj);
 
                     matrPair.first(i, npJ + j) = i01 & nrm[i];
                     matrPair.second(i, npJ + j) = i01 & taui;
 
                     
                     v10 = {
                         -0.5 / di.length() * (((s1 + s2) & taui) * VMlib::Omega(s1, taui, tauj) - s1.length2() * tauj),
                         0.5 * dj.length() / di.length() * VMlib::Omega(s1 + p2, taui, taui)
                     };
 
                     i10 = IDPI / len[i] * (-((alpha[0] + alpha[2]) * v10[0] + alpha[2] * v10[1]).kcross() \
                         + ((lambda[0] + lambda[2]) * v10[0] + lambda[2] * v10[1]) + 0.5 * dj.length() * taui);
 
                     matrPair.first(npI + i, j) = i10 & nrm[i];
                     matrPair.second(npI + i, j) = i10 & taui;
 
 
                     v11 = {
                         1.0 / (12.0 * di.length() * dj.length()) * (2.0 * (s1 & VMlib::Omega(s1 - 3.0 * p2, taui, tauj)) * VMlib::Omega(s1, taui, tauj) - s1.length2() * (s1 - 3.0 * p2)) - 0.25 * VMlib::Omega(s1, taui, tauj),
                         -di.length() / (12.0 * dj.length()) * VMlib::Omega(di, tauj, tauj),
                         -dj.length() / (12.0 * di.length()) * VMlib::Omega(dj, taui, taui)
                     };
 
                     i11 = IDPI / len[i] * (-((alpha[0] + alpha[2]) * v11[0] + (alpha[1] + alpha[2]) * v11[1] + alpha[2] * v11[2]).kcross()\
                         + (lambda[0] + lambda[2]) * v11[0] + (lambda[1] + lambda[2]) * v11[1] + lambda[2] * v11[2] \
                         + 1.0 / 12.0 * (dj.length() * taui + di.length() * tauj - 2.0 * VMlib::Omega(s1, taui, tauj)));
 
                     matrPair.first(npI + i, npJ + j)  = i11 & nrm[i];
                     matrPair.second(npI + i, npJ + j) = i11 & taui;
                 }
 
 
                 if (p > 2)
                     throw (-42);
             }//else(i == j)
 
         }//for(...)
 }//getIQ(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void Airfoil::calcMeanEpsOverPanel ( )

inherited

Вычисление средних значений eps на панелях

Definition at line 82 of file Airfoil2D.cpp.

 {
     meanEpsOverPanel.clear();
     meanEpsOverPanel.resize(getNumberOfPanels());
 
 
     double midEps;
 
     const Boundary& bnd = W.getBoundary(numberInPassport);
     VortexesParams virtVortParams = W.getVelocity().virtualVortexesParams[numberInPassport];
 
     for (size_t i = 0; i < getNumberOfPanels(); ++i)
     {
         /*
         midEps = 0.0;
 
         for (int j = bnd.vortexBeginEnd[i].first; j < bnd.vortexBeginEnd[i].second; ++j)
             midEps += virtVortParams.epsastWake[j];
 
         midEps /= (bnd.vortexBeginEnd[i].second - bnd.vortexBeginEnd[i].first);
 
         meanEpsOverPanel[i] = midEps;
         */
 
         /*
         midEps = 0.0;
 
         for (int j = bnd.vortexBeginEnd[i].first; j < bnd.vortexBeginEnd[i].second; ++j)
             midEps += std::max(virtVortParams.epsastWake[j], 0.5 * len[i] / (bnd.vortexBeginEnd[i].second - bnd.vortexBeginEnd[i].first));
 
         midEps /= (bnd.vortexBeginEnd[i].second - bnd.vortexBeginEnd[i].first);
         meanEpsOverPanel[i] = midEps;
         //*/
 
         midEps = 0.0;
         for (int j = bnd.vortexBeginEnd[i].first; j < bnd.vortexBeginEnd[i].second; ++j)
             midEps += virtVortParams.epsastWake[j];
 
         midEps /= (bnd.vortexBeginEnd[i].second - bnd.vortexBeginEnd[i].first);
         meanEpsOverPanel[i] = midEps;
     }//for i
 }

Here is the call graph for this function:





Here is the caller graph for this function:

void AirfoilRect::CalcNrmTauLen ( )

Вычисление нормалей, касательных и длин панелей по текущему положению вершин

Definition at line 616 of file Airfoil2DRect.cpp.

 {
     if (nrm.size() != r_.size())
     {
         nrm.resize(r_.size());
         tau.resize(r_.size());
         len.resize(r_.size());
     }
 
     Point2D rpan;
 
     for (size_t i = 0; i < r_.size(); ++i)
     {
         rpan = (getR(i + 1) - getR(i));
         len[i] = rpan.length();
         tau[i] = rpan.unit();
 
         nrm[i] = { tau[i][1], -tau[i][0] };
     }
 
     //закольцовываем
     nrm.push_back(nrm[0]);
     tau.push_back(tau[0]);
     len.push_back(len[0]);
 }//CalcNrmTauLen()

Here is the call graph for this function:

Here is the caller graph for this function:

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

overridevirtual

Вычисление коэффициентов матрицы A для расчета влияния панели на панель

Parameters

[in]	p	размерность матрицы - результата
[in]	i	номер панели, на которую оказывается влияние
[in]	airfoil	константная ссылка на профиль, от которого рассчитывается влияние
[in]	j	номер влияющей панели return соответствующий блок матрицы СЛАУ, вытянутый в линию

Implements VM2D::Airfoil.

Definition at line 665 of file Airfoil2DRect.cpp.

 {
     std::vector<double> res(p*p, 0.0);
 
     if ( (i == j) && (&airfoil == this))
     {
         res[0] = 0.5 * (tau[i] ^ nrm[i]);
         if (p == 1)
             return res;
 
         res[3] = (1.0 / 12.0) * res[0];
         if (p == 2)
             return res;
 
         if (p > 2)
             throw (-42);
 
     }//if i==j
         
     const auto& miq = W.getIQ(numberInPassport, airfoil.numberInPassport);
 
     res[0] = miq.first(i, j);
 
     if (p == 1)
         return res;
 
     res[1] = miq.first(i, airfoil.getNumberOfPanels() + j);
 
     res[2] = miq.first(getNumberOfPanels() + i, j);
 
     res[3] = miq.first(getNumberOfPanels() + i, airfoil.getNumberOfPanels() + j);
 
     if (p == 2)
         return res;
 
     if (p > 2)
         throw (-42);
     
     //Хотя сюда никогда и не попадем
     return res;
 }//getA(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void AirfoilRect::GetDiffVelocityI0I3ToSetOfPointsAndViscousStresses	(	const WakeDataBase &	pointsDb,
		std::vector< double > &	domainRadius,
		std::vector< double > &	I0,
		std::vector< Point2D > &	I3
	)

overridevirtual

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

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

Parameters

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

Warning: Векторы I0, I3 — накапливаются!

Implements VM2D::Airfoil.

Definition at line 303 of file Airfoil2DRect.cpp.

 {
     std::vector<double> selfI0;
     std::vector<Point2D> selfI3;
 
     std::vector<double> currViscousStress;
     currViscousStress.resize(r_.size(), 0.0);
 
 
     selfI0.resize(pointsDb.vtx.size(), 0.0);
     selfI3.resize(pointsDb.vtx.size(), { 0.0, 0.0 });
 
 #pragma warning (push)
 #pragma warning (disable: 4101)
     //Локальные переменные для цикла
     Point2D q, xi, xi_m, v0;
     double lxi, lxi_m, lenj_m;
 
     Point2D mn;
     int new_n;
     Point2D h;
 
     Point2D vec;
     double s, d;
 
     double vs;
     double iDDomRad, domRad, expon;
 #pragma warning (pop)
 
 #pragma omp parallel for \
     default(none) \
     shared(domainRadius, pointsDb, selfI0, selfI3, currViscousStress, PI) \
     private(xi, xi_m, lxi, lxi_m, lenj_m, v0, q, new_n, mn, h, d, s, vec, vs, expon, domRad, iDDomRad) schedule(dynamic, DYN_SCHEDULE)
     for (int i = 0; i < pointsDb.vtx.size(); ++i)
     {
         const Vortex2D& vtxI = pointsDb.vtx[i];
         
         domRad = std::max(domainRadius[i], W.getPassport().wakeDiscretizationProperties.getMinEpsAst());
         iDDomRad = 1.0 / domRad;
 
         for (size_t j = 0; j < r_.size(); ++j)
         {
             vs = 0.0;
             q = vtxI.r() - 0.5 * (getR(j) + getR(j + 1));
             vec = tau[j];
 
             s = q & vec;
             d = fabs(q & nrm[j]);
 
             if ((d < 50.0 * len[j]) && (fabs(s) < 50.0 * len[j]))
             {
                 v0 = vec * len[j];
 
                 if ((d > 5.0 * len[j]) || (fabs(s) > 5.0 * len[j]))
                 {
                     xi = q * iDDomRad;
                     lxi = xi.length();
 
                     expon = exp(-lxi) * len[j];
                     mn = nrm[j] * expon;
 
                     if (selfI0[i] != -PI * domRad)
                     {
                         selfI0[i] += (xi & mn) * (lxi + 1.0) / (lxi * lxi);
                         selfI3[i] += mn;
                     }
 
                     vs = vtxI.g() * expon / (PI * sqr(meanEpsOverPanel[j]));
                 }
                 else if ((d >= 0.1 * len[j]) || (fabs(s) > 0.5 * len[j]))
                 {
                     vs = 0.0;
                     double den = (fabs(s) < 0.5 * len[j]) ? d : (fabs(s) + d - 0.5 * len[j]);
 
                     new_n = std::max(1, static_cast<int>(ceil(5.0 * len[j] / den)));
                     //new_n = 100;
                     h = v0 * (1.0 / new_n);
 
                     for (int m = 0; m < new_n; ++m)
                     {
                         xi_m = (vtxI.r() - (getR(j) + h * (m + 0.5))) * iDDomRad;
                         lxi_m = xi_m.length();
 
                         lenj_m = len[j] / new_n;
                         expon = exp(-lxi_m) * lenj_m;
 
                         mn = nrm[j] * expon;
                         if (selfI0[i] != -PI * domRad)
                         {
                             selfI0[i] += (xi_m & mn) * (lxi_m + 1.0) / (lxi_m * lxi_m);
                             selfI3[i] += mn;
                         }
                         vs += expon;
                     }//for m
                     vs *= vtxI.g() / (PI * sqr(meanEpsOverPanel[j]));
                 }
                 else
                 {
 
                     selfI0[i] = -PI * domRad;
                     double mnog = 2.0 * domRad * (1.0 - exp(-len[j] * iDDomRad / 2.0) * cosh(fabs(s) * iDDomRad));
                     selfI3[i] = nrm[j] * mnog;
                     vs = mnog * vtxI.g() / (PI * sqr(meanEpsOverPanel[j]));
                 }
             }//if d<50 len 
 
 #pragma omp atomic
             currViscousStress[j] += vs;
         }//for j
     }//for i
 
 
     if (&pointsDb == &(W.getWake()))
         for (size_t i = 0; i < viscousStress.size(); ++i)
         {
             viscousStress[i] += currViscousStress[i];
         }
 
     for (size_t i = 0; i < I0.size(); ++i)
     {
         domRad = std::max(domainRadius[i], W.getPassport().wakeDiscretizationProperties.getMinEpsAst());
 
         if (I0[i] != -PI * domRad)
         {
             if (selfI0[i] == -PI * domRad)
             {
                 I0[i] = selfI0[i];
                 I3[i] = selfI3[i];
             }
             else
             {
                 I0[i] += selfI0[i];
                 I3[i] += selfI3[i];
             }
         }
     }
 }; //GetDiffVelocityI0I3ToSetOfPointsAndViscousStresses(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void AirfoilRect::GetDiffVelocityI0I3ToWakeAndViscousStresses	(	const WakeDataBase &	pointsDb,
		std::vector< double > &	domainRadius,
		std::vector< double > &	I0,
		std::vector< Point2D > &	I3
	)

overridevirtual

Implements VM2D::Airfoil.

Definition at line 442 of file Airfoil2DRect.cpp.

 {
     if (&pointsDb == &(W.getWake()))
     {
         viscousStress.clear();
         viscousStress.resize(r_.size(), 0.0);
     }
 
     GetDiffVelocityI0I3ToSetOfPointsAndViscousStresses(pointsDb, domainRadius, I0, I3);
 }//GetDiffVelocityI0I3ToWakeAndViscousStresses(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void AirfoilRect::GetGabarits ( double gap = 0.02 )

overridevirtual

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

Implements VM2D::Airfoil.

Definition at line 643 of file Airfoil2DRect.cpp.

 {
     lowLeft = { 1E+10, 1E+10 };
     upRight = { -1E+10, -1E+10 };
 
     for (size_t i = 0; i < r_.size(); ++i)
     {
         lowLeft[0] = std::min(lowLeft[0], r_[i][0]);
         lowLeft[1] = std::min(lowLeft[1], r_[i][1]);
 
         upRight[0] = std::max(upRight[0], r_[i][0]);
         upRight[1] = std::max(upRight[1], r_[i][1]);
     }
 
     Point2D size = upRight - lowLeft;
     lowLeft -= size*gap;
     upRight += size*gap;
 }//GetGabarits(...)

Here is the caller graph for this function:

void AirfoilRect::GetInfluenceFromSourceSheetToVortex	(	size_t	panel,
		const Vortex2D &	vtx,
		Point2D &	vel
	)		const

overridevirtual

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

Parameters

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

Implements VM2D::Airfoil.

Definition at line 855 of file Airfoil2DRect.cpp.

 {
     W.getBoundary(numberInPassport).GetInfluenceFromSourceSheetAtRectPanelToVortex(panel, vtx, vel);
 }//GetInfluenceFromSourceSheetToVortex(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void AirfoilRect::GetInfluenceFromSourcesToPanel	(	size_t	panel,
		const Vortex2D *	ptr,
		ptrdiff_t	count,
		std::vector< double > &	panelRhs
	)		const

overridevirtual

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

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

Parameters

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

Implements VM2D::Airfoil.

Definition at line 849 of file Airfoil2DRect.cpp.

 {
     W.getBoundary(numberInPassport).GetInfluenceFromSourcesToRectPanel(panel, ptr, count, panelRhs);
 }//GetInfluenceFromSourcesToPanel(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void AirfoilRect::GetInfluenceFromVInfToPanel ( std::vector< double > & vInfRhs ) const

overridevirtual

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

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

Parameters

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

Implements VM2D::Airfoil.

Definition at line 866 of file Airfoil2DRect.cpp.

 {
     W.getBoundary(numberInPassport).GetInfluenceFromVInfToRectPanel(vInfRhs);
 }//GetInfluenceFromVInfToPanel(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void AirfoilRect::GetInfluenceFromVortexSheetToVortex	(	size_t	panel,
		const Vortex2D &	vtx,
		Point2D &	vel
	)		const

overridevirtual

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

Parameters

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

Implements VM2D::Airfoil.

Definition at line 861 of file Airfoil2DRect.cpp.

 {
     W.getBoundary(numberInPassport).GetInfluenceFromVortexSheetAtRectPanelToVortex(panel, vtx, vel);
 }//GetInfluenceFromVortexSheetToVortex(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void AirfoilRect::GetInfluenceFromVorticesToPanel	(	size_t	panel,
		const Vortex2D *	ptr,
		ptrdiff_t	count,
		std::vector< double > &	panelRhs
	)		const

overridevirtual

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

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

Parameters

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

Implements VM2D::Airfoil.

Definition at line 842 of file Airfoil2DRect.cpp.

 {
     W.getBoundary(numberInPassport).GetInfluenceFromVorticesToRectPanel(panel, ptr, count, panelRhs);
 }//GetInfluenceFromVorticesToPanel(...)

Here is the call graph for this function:

Here is the caller graph for this function:

size_t VM2D::Airfoil::getNumberOfPanels ( ) const

inlineinherited

Возврат количества панелей на профиле

Returns: количество панелей на профиле

Definition at line 139 of file Airfoil2D.h.

139 { return r_.size(); };

VM2D::Airfoil::r_

std::vector< Point2D > r_

Координаты начал панелей

Definition: Airfoil2D.h:64

Here is the caller graph for this function:

const Point2D& VM2D::Airfoil::getR ( size_t q ) const

inlineinherited

Возврат константной ссылки на вершину профиля

Организовано "зацикливание" в сторону увеличения индекса, т.е. getR[size()] = getR[0];
Это позволяет удобно обращаться к getR(i) и getR(i+1) как к началу и концу i-й панели

Parameters

[in] q номер вершины профиля return константную ссылку на вершину профиля

Definition at line 101 of file Airfoil2D.h.

         {
             return (q < r_.size()) ? r_[q] : r_[0];
         };

Here is the caller graph for this function:

const Point2D& VM2D::Airfoil::getV ( size_t q ) const

inlineinherited

Возврат константной ссылки на скорость вершины профиля

Организовано "зацикливание" в сторону увеличения индекса, т.е. getV[size()] = getV[0];
Это позволяет удобно обращаться к getV(i) и getV(i+1) как к скоростям начала и конца i-й панели

Parameters

[in] q номер вершины профиля return константную ссылку на скорость вершины профиля

Definition at line 113 of file Airfoil2D.h.

         {
             return (q < v_.size()) ? v_[q] : v_[0];
         };

Here is the caller graph for this function:

VM2D::Airfoil::IFCUDA ( mutable double * devRPtr )

inherited

Указатель на девайсе, где хранятся вершины профиля

VM2D::Airfoil::IFCUDA ( mutable double * devRhsPtr )

inherited

Указатель на девайсе, где хранится правая часть (константная) матрицы

VM2D::Airfoil::IFCUDA ( mutable double * devRhsLinPtr )

inherited

Указатель на девайсе, где хранится правая часть (линейная) матрицы

VM2D::Airfoil::IFCUDA ( mutable std::vector< double > tmpRhs )

inherited

Указатель на хосте, где хранится временная часть матрицы, полученная с девайса

VM2D::Airfoil::IFCUDA ( mutable double * devFreeVortexSheetPtr )

inherited

Указатель на девайсе, где хранятся интенсивности свободного вихревого слоя на панелях

VM2D::Airfoil::IFCUDA ( mutable double * devAttachedVortexSheetPtr )

inherited

Указатель на девайсе, где хранятся интенсивности присоединенного вихревого слоя на панелях

VM2D::Airfoil::IFCUDA ( mutable double * devAttachedSourceSheetPtr )

inherited

Указатель на девайсе, где хранятся интенсивности присоединенного слоя источников на панелях

VM2D::Airfoil::IFCUDA ( mutable double * devMeanEpsOverPanelPtr )

inherited

Указатель на девайсе, где хранятся средние eps на панелях

VM2D::Airfoil::IFCUDA ( mutable double * devViscousStressesPtr )

inherited

Указатель на девайсе, где хранится вектор (по панелям) для силы вязкого трения

VM2D::Airfoil::IFCUDA ( mutable std::vector< double > tmpViscousStresses )

inherited

Указатель на хосте, где хранится временная часть вектора (по панелям) для силы вязкого трения

bool Airfoil::isAfter	(	size_t	i,
		size_t	j
	)		const

inherited

Проверка, идет ли вершина i следом за вершиной j.

Parameters

[in]	i	проверяемая вершина
[in]	j	контрольная вершина

Returns: true, если i-я вершина следует зп j-й в порядке обхода профиля

Definition at line 62 of file Airfoil2D.cpp.

 {
     return ((i == j + 1) || (i == 0 && j == r_.size() - 1));
 }//isAfter(...)

Here is the caller graph for this function:

bool Airfoil::isInsideGabarits ( const Point2D & r ) const

inherited

Определяет, находится ли точка с радиус-вектором \( \vec r \) внутри габаритного прямоугольника профиля

Parameters

[in] r константная ссылка на текущее положение точки

Returns: true, если точка внутри габаритного прямоугольника

Definition at line 69 of file Airfoil2D.cpp.

 {
     return (r[0] <= upRight[0] && (r[0] >= lowLeft[0] && r[1] >= lowLeft[1] && r[1] <= upRight[1]));
 }//isInsideGabarits(...)

Here is the caller graph for this function:

bool Airfoil::isOutsideGabarits ( const Point2D & r ) const

inherited

Определяет, находится ли точка с радиус-вектором \( \vec r \) вне габаритного прямоугольника профиля

Parameters

[in] r константная ссылка на текущее положение точки

Returns: true, если точка вне габаритного прямоугольника

Definition at line 76 of file Airfoil2D.cpp.

 {
     return (r[0] > upRight[0] || (r[0] < lowLeft[0] || r[1] < lowLeft[1] || r[1] > upRight[1]));
 }//isOutsideGabarits(...)

Here is the caller graph for this function:

bool AirfoilRect::IsPointInAirfoil ( const Point2D & point ) const

overridevirtual

Определяет, находится ли точка с радиус-вектором \( \vec r \) внутри профиля

Parameters

[in] point константная ссылка на текущее положение точки

Returns: true, если точка внутри профиля

Implements VM2D::Airfoil.

Definition at line 560 of file Airfoil2DRect.cpp.

 {
     double sumAngle = 0.0;
 
     Point2D v1, v2;
 
     for (size_t i = 0; i < r_.size(); ++i)
     {
         v1 = getR(i) - point;
         v2 = getR(i+1) - point;
         sumAngle += atan2(v1^v2, v1 & v2);
     }
 
     if (fabs(sumAngle) < 0.1)
         return false;
     else 
         return true;
 }//IsPointInAirfoil(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void AirfoilRect::Move ( const Point2D & dr )

overridevirtual

Перемещение профиля

Плоскопараллельно перемещает профиль на вектор \( \overrightarrow{dr} \)

Parameters

[in] dr константная ссылка на вектор перемещения

Implements VM2D::Airfoil.

Definition at line 580 of file Airfoil2DRect.cpp.

 {
     for (size_t i = 0; i < r_.size(); ++i)
         r_[i] += dr;
     rcm += dr;
     CalcNrmTauLen();
     GetGabarits();
 }//Move(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void AirfoilRect::ReadFromFile ( const std::string & dir )

overridevirtual

Считывание профиля из файла

Считывание геометрии профиля из файла, вычисление всех прочих параметров профиля
После загрузки из файла профиль поворачивается на нужный угол и масштабируется на нужный коэффициент

Warning: Сейчас масса, момент инерции и скорости вершин зануляются.

Parameters

[in] dir константная ссылка на строку — имя каталога, где лежит cчитываемый файл

Implements VM2D::Airfoil.

Definition at line 59 of file Airfoil2DRect.cpp.

 {
     const AirfoilParams& param = W.getPassport().airfoilParams[numberInPassport];
     std::string filename = dir + param.fileAirfoil;
     std::ifstream airfoilFile;
     
     if (fileExistTest(filename, W.getInfo(), { "txt", "TXT" }))
     {
         std::stringstream airfoilFile(VMlib::Preprocessor(filename).resultString);
 
         VMlib::StreamParser airfoilParser(W.getInfo(), "airfoil file parser", airfoilFile);
 
         m = 0.0; //TODO
         J = 0.0; //TODO
         rcm = { 0.0, 0.0 }; //TODO
         phiAfl = 0.0;
 
         inverse = param.inverse;
 
 
 
         //*
 
         //Проверяем отсутствие "задвоенных" подряд идущих точек
         //airfoilParser.get("r", r_);
         std::vector<Point2D> rFromFile;
         airfoilParser.get("r", rFromFile);
 
         if (rFromFile.size() > 0)
         {
             if (param.requiredNPanels <= rFromFile.size()) {
 
                 r_.reserve(rFromFile.size());
                 r_.push_back(rFromFile[0]);
                 for (size_t i = 1; i < rFromFile.size(); ++i)
                     if ((rFromFile[i] - rFromFile[i - 1]).length2() > 1e-12)
                         r_.push_back(rFromFile[i]);
 
                 //Если первая совпадает с последней, то убираем последнюю
                 if ((r_.back() - r_.front()).length2() < 1e-12)
                     r_.resize(r_.size() - 1);
             }
             else 
             { 
                 W.getInfo('e') << "Airfoil shape is given explicitely, it can not be automatically split!" << std::endl; 
                 exit(200);
             }
         }
         else
         {
             //std::vector<GeomPoint> geomFromFile(X.size()); //new
             std::vector<GeomPoint> geomFromFile; //VM2D
 
             //for (int q = 0; q < X.size(); ++q)            // new
             //{                                             //
             //  geomFromFile[q].r = { X[q], Y[q] };         //
             //  geomFromFile[q].type = key[q];              //
             //}                                             //
             airfoilParser.get("geometry", geomFromFile);  //VM2D
 
             if ((geomFromFile.front().r - geomFromFile.back().r).length2() < 1e-12)
                 geomFromFile.resize(geomFromFile.size() - 1);
 
 
             //size_t reqN = NumberOfPoints;                 //new
             size_t reqN = param.requiredNPanels;        //VM2D
 
             std::vector<double> L(geomFromFile.size());
             double totalLength = 0.0;
             std::vector<size_t> nc = {};
             std::vector<size_t> ni = {};
 
             for (size_t i = 0; i < geomFromFile.size(); ++i)
             {
                 const Point2D& p1 = geomFromFile[i].r;
                 const Point2D& p2 = geomFromFile[(i + 1) % geomFromFile.size()].r;
 
                 L[i] = (p2 - p1).length();
                 totalLength += L[i];
             }
 
 
             //Ищем первый сплайн
             size_t i = 0;
             std::vector<size_t> splineStart, splineFinish;
 
             while (i < geomFromFile.size())
             {
                 while (i < geomFromFile.size() && !(geomFromFile[i].type == "c"))
                     ++i;
 
                 if (i < geomFromFile.size())
                 {
                     splineStart.push_back(i);
                     ++i;
                     //ищем второй конец
                     while (i < geomFromFile.size() && !(geomFromFile[i].type == "c"))
                         ++i;
 
                     splineFinish.push_back(i);
                 }
             }
 
             if ((splineStart.size() > 0) && (splineStart[0] != 0))
                 splineFinish.back() = splineStart[0];
 
             //std::vector<Point2D> r_; //new
 
             if (reqN == 0)
             {
                 r_.reserve(geomFromFile.size());
                 for (size_t i = 0; i < geomFromFile.size(); ++i)
                     r_.push_back(geomFromFile[i].r);
             }
             else
             {
                 double hRef = totalLength / reqN;
 
                 std::vector<int> nPanels;
 
                 bool cyclic = (splineStart.size() == 0);
 
                 if (cyclic)
                 {
                     splineStart.push_back(0);
                     splineFinish.push_back(geomFromFile.size());
                 }
 
                 for (size_t s = 0; s < splineStart.size(); ++s)
                 {
                     //Сплайн                
                     std::vector<double> T, X, Y;
 
                     size_t splineLegs = splineFinish[s] - splineStart[s];
                     if (splineFinish[s] <= splineStart[s])
                         splineLegs += geomFromFile.size();
 
                     T.reserve(splineLegs + 1);
                     X.reserve(T.size());
                     Y.reserve(T.size());
 
                     double tbuf = 0.0;
 
                     for (size_t i = splineStart[s]; i <= splineStart[s] + splineLegs; ++i)
                     {
                         const Point2D& pt = geomFromFile[i % geomFromFile.size()].r;
 
                         T.push_back(tbuf);
                         X.push_back(pt[0]);
                         Y.push_back(pt[1]);
 
                         if ((i == splineStart[s]) && (splineFinish[s] - splineStart[s] == 1))
                         {
                             double halfL = (i < geomFromFile.size()) ? 0.5 * L[i] : 0.5 * (geomFromFile.front().r - geomFromFile.back().r).length();
                             tbuf += halfL;
                             T.push_back(tbuf);
                             Point2D nextPt = geomFromFile[(i + 1) % geomFromFile.size()].r;
 
                             X.push_back(0.5 * (pt[0] + nextPt[0]));
                             Y.push_back(0.5 * (pt[1] + nextPt[1]));
 
                             tbuf += halfL;
                         }
                         else
                             tbuf += (i < geomFromFile.size()) ? L[i] : (geomFromFile.front().r - geomFromFile.back().r).length();
                     }
 
 
 
                     tk::spline s1, s2;
                     if (cyclic)
                     {
                         s1.set_boundary(tk::spline::cyclic);
                         s2.set_boundary(tk::spline::cyclic);
                     }
                     else
                     {
                         s1.set_boundary(tk::spline::second_deriv, 0.0, tk::spline::second_deriv, 0.0);
                         s2.set_boundary(tk::spline::second_deriv, 0.0, tk::spline::second_deriv, 0.0);
                     }
 
                     s1.set_points(T, X);
                     s2.set_points(T, Y);
 
                     int NumberOfPanelsSpline = (int)ceil(T.back() / hRef);
                     double hSpline = T.back() / NumberOfPanelsSpline;
 
                     for (int i = 0; i < NumberOfPanelsSpline; ++i)
                         r_.push_back({ s1(hSpline * i), s2(hSpline * i) });
 
                     nPanels.push_back(NumberOfPanelsSpline);
 
                 }
             }
         }//else geometry
 
 
 
         //std::ofstream of("airfoil.txt");      
         //for (size_t i = 0; i < r_.size(); ++i)
         //  of << r_[i] << "," << std::endl;
         //of.close();
         
         if (r_.size() == 0)
         {
             W.getInfo('e') << "No points on airfoil contour!" << std::endl;
             exit(200);
         }       
 
         if (inverse)
             std::reverse(r_.begin(), r_.end());
 
 
         v_.resize(0);
         for (size_t q = 0; q < r_.size(); ++q)
             v_.push_back({ 0.0, 0.0 });     
 
         Move(param.basePoint);
         Scale(param.scale);
         
         double rotationAngle = param.angle;
         if (W.getPassport().geographicalAngles)
             rotationAngle = -param.angle - 0.5 * PI;
 
         Rotate(-rotationAngle);
         //в конце Rotate нормали, касательные и длины вычисляются сами
 
         gammaThrough.clear();
         gammaThrough.resize(r_.size(), 0.0);
 
         //Вычисляем площадь
         area = 0.0;
         for (size_t q = 0; q < r_.size(); ++q)
         {
             const Point2D& cntq = 0.5 * (getR(q) + getR(q + 1));
             const Point2D& drq = getR(q + 1) - getR(q);
             area += fabs(cntq[0] * drq[1] - cntq[1] * drq[0]);
         }
         area *= 0.5;
     }
 }//ReadFromFile(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void AirfoilRect::Rotate ( double alpha )

overridevirtual

Поворот профиля

Поворачивает профиль на угол \( \alpha \) вокруг центра масс

Parameters

[in] alpha угол поворота против часовой стрелки в радианах

Implements VM2D::Airfoil.

Definition at line 591 of file Airfoil2DRect.cpp.

 {
     phiAfl += alpha;
     nummatrix<double, 2, 2> rotMatrix = { { cos(alpha), -sin(alpha) }, { sin(alpha), cos(alpha) } };
 
     for (size_t i = 0; i < r_.size(); ++i)
         r_[i] = rcm + (rotMatrix & (r_[i] - rcm));
 
 
     CalcNrmTauLen();
     GetGabarits();
 }//Rotate(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void AirfoilRect::Scale ( const Point2D & )

overridevirtual

Масштабирование профиля

Масштабирует профиль в factor раз относительно центра масс

Parameters

[in] factor масштабный коэффициент

Implements VM2D::Airfoil.

Definition at line 606 of file Airfoil2DRect.cpp.

 {
     for (size_t i = 0; i < r_.size(); ++i)
         r_[i] = rcm + Point2D{ factor[0] * (r_[i] - rcm)[0], factor[1] * (r_[i] - rcm)[1] };
 
     CalcNrmTauLen();
     GetGabarits();
 }//Scale(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void VM2D::Airfoil::setV ( const Point2D & vel )

inlineinherited

Установка постоянной скорости всех вершин профиля

Parameters

[in] vel константная ссылка на величину устанавливаемой скорости

Definition at line 121 of file Airfoil2D.h.

         {
             v_.clear();
             v_.resize(r_.size(), vel);
         }

Here is the caller graph for this function:

void VM2D::Airfoil::setV ( const std::vector< Point2D > & vel )

inlineinherited

Установка скоростей всех вершин профиля

Parameters

[in] vel константная ссылка на вектор величин скоростей вершин профиля

Definition at line 130 of file Airfoil2D.h.

         {
             v_.clear();
             v_.insert(v_.end(), vel.begin(), vel.end());                
         }

Member Data Documentation

double VM2D::Airfoil::area

inherited

Площадь профиля

Definition at line 83 of file Airfoil2D.h.

std::vector<double> VM2D::Airfoil::gammaThrough

inherited

Суммарные циркуляции вихрей, пересекших панели профиля на прошлом шаге

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

Definition at line 196 of file Airfoil2D.h.

bool VM2D::Airfoil::inverse

inherited

Признак разворота нормалей (для расчета внутренних течений)

Definition at line 139 of file Airfoil2D.h.

double VM2D::Airfoil::J

inherited

Полярный момент инерции профиля относительно центра масс

Definition at line 89 of file Airfoil2D.h.

std::vector<double> VM2D::Airfoil::len

inherited

Длины панелей профиля

Definition at line 185 of file Airfoil2D.h.

Point2D VM2D::Airfoil::lowLeft

inherited

Левый нижний угол габаритного прямоугольника профиля

Definition at line 190 of file Airfoil2D.h.

double VM2D::Airfoil::m

inherited

Масса профиля

Definition at line 86 of file Airfoil2D.h.

std::vector<double> VM2D::Airfoil::meanEpsOverPanel

inherited

Средние значения Eps на панелях

Definition at line 92 of file Airfoil2D.h.

std::vector<Point2D> VM2D::Airfoil::nrm

inherited

Нормали к панелям профиля

Нормали задаются внешними, нормированными на единицу

Definition at line 177 of file Airfoil2D.h.

const size_t VM2D::Airfoil::numberInPassport

inherited

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

Definition at line 74 of file Airfoil2D.h.

double VM2D::Airfoil::phiAfl

inherited

Поворот профиля

Definition at line 80 of file Airfoil2D.h.

std::vector<Point2D> VM2D::Airfoil::r_

protectedinherited

Координаты начал панелей

Definition at line 64 of file Airfoil2D.h.

Point2D VM2D::Airfoil::rcm

inherited

Положение центра масс профиля

Definition at line 77 of file Airfoil2D.h.

std::vector<Point2D> VM2D::Airfoil::tau

inherited

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

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

Definition at line 182 of file Airfoil2D.h.

Point2D VM2D::Airfoil::upRight

inherited

Правый верхний угол габаритного прямоугольника профиля

Definition at line 191 of file Airfoil2D.h.

std::vector<Point2D> VM2D::Airfoil::v_

protectedinherited

Скорости начал панелей

Definition at line 67 of file Airfoil2D.h.

std::vector<double> VM2D::Airfoil::viscousStress

inherited

Касательные напряжения на панелях профиля

Definition at line 188 of file Airfoil2D.h.

const World2D& VM2D::Airfoil::W

inherited

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

Definition at line 71 of file Airfoil2D.h.

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

VM2D/Airfoil2D/Airfoil2DRect.h
VM2D/Airfoil2D/Airfoil2DRect.cpp

Public Member Functions

Public Attributes

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation