Страница на редактировании !!!

.: программа InfSWing

Рассчитывает распределенные и интегральные характеристики бесконечного стреловидного крыла, на заданном числе Маха, угле атаки и стреловидности, включая трансзвуковые режимы. Расчет обтекания осуществляется путем численного интегрирования консервативной формы полного уравнения для потенциала методом неполной факторизации с приближенным учетом неизэнтропичности потока на скачках уплотнения.
Для расчета сетки используется простой алгебраический сеточный генератор.

Интерфейс построен на основе библиотеки Fox ToolKit.

.: инструкция к программе InfSWing

.: панель InfSWing Calculation:

Из меню [File] -> [Open] нужно выбрать исходный файл.

после выбора исходного файла будет доступна вкладка [Grid]

.: вкладка Grid:

NPA число узлов на поверхности профиля
NCO число узлов по нормали к поверхности профиля
NBH число узлов в следовой части
Kgrid (=1,2 или 4) коэффициент измельчения сетки
Vortex sheet Hj+1/Hj в следовой части
Mid-chord Hj+1/Hj в средней части профиля
Leading edge Hj+1/Hj в области передней кромки

Кнопка [Grid Calculation] - расчет сетки.

Кнопка [Grid View ] - визуализация сетки.

после расчета сетки будет доступна вкладка [Calculation]

.: вкладка Calculation:

Alpha угол атаки
Mach число Маха
Swept максимальная толщина
Diagonal dom. диагональное доминирование
Max. residual максимальная невязка
Linear iterations число линейных итераций
GMRES dim. GMRES размерность
Clift коэффициент подъемной силы
Cdragкоэффициент сопротивления (интеграл Cp)

Кнопка [Calculation] - расчет

Кнопка [Cp View ] - визуализация Cp.

Кнопка [M View ] - визуализация числа Маха в расчетной области.

Кнопка [Cp Edit] - окрывает панель [Cp Editor 2D] редактирования Cp для обратной задачи.

Левая кнопка мыши смещение
Левая кнопка мыши + Shift увеличение

Этими кнопками можно упрравлять видом сетки и М

.: панель редактирования Cp:

Cyt коэффициент подъемной силы (исходный Cp)
Cyd коэффициент подъемной силы (целевой Cp)
Nspl число узлов сплайна на интервале редактирования
U1 U2 интервал редактирования на верхней поверхности
L1 L2 интервал редактирования на нижней поверхности
левый клик мыши определяет узел редактирования,
выбранный узел становится красным
движение нажатой левой клавиши мыши на графике смещает соответствующий узел

Кнопка [Cp design save] - сохраняет результат редактирования и возвращает на вкладку [Inverse] основной панели

.: вкладка Inverse:

Leading boundary целевое распределение Cp
Global iteration целевое распределение Cp
Cp deviation максимальное отклонение Cp от целевого
Clif current текущее значение Cy
Clift target целевое значение Cy
Режим решения обратной задачи
Design both surfaces проектирование верхней и нижней поверхностей
Maintain thickness distribution проектирование верхней поверхности и сохранение распределения толщины
Maintain maximum thickness проектирование верхней поверхности и сохранение максимальной толщины

Black line - исходное распределение Cp

Red line - целевое распределение Cp

Каждое нажатие на кнопку [Calculation] добавляет одну глобальную итерацию


Глобальная итерация состоит из трех этапов:

  1. решение обратной задачи с помощью нормальной скорости на поверхности профиля обеспечивающего заданное распределеня давления .
  2. модификация геометрии профиля.
  3. прямой расчет модифицированного профиля.

.: Пример исходного файла:

RAE-2822
MACH ALFA HEE
0.72 1.5 0.00000
BETA ERM CLU
.01 1.E-5 0.04
FNPA FNCO FNBH FKGR
16. 8. 12. 2.
FNLN FMGT FICY FIDG
5. 3. 0. 0.
FNF FNP
0. 33.
XU YU
1.000000 0.000935
0.986047 0.004450
0.953994 0.009902
0.905037 0.019071
0.840999 0.031107
0.764266 0.043588
0.677698 0.055109
0.584520 0.064556
0.488203 0.071295
0.392337 0.074861
0.300493 0.074846
0.216094 0.071192
0.142284 0.064244
0.081813 0.054361
0.036934 0.041518
0.009320 0.023419
0.000000-0.000000
0.009320-0.021822
0.036934-0.039221
0.081813-0.052897
0.142284-0.063740
0.216094-0.071305
0.300493-0.074654
0.392337-0.073152
0.488203-0.066250
0.584520-0.053968
0.677698-0.037751
0.764266-0.021036
0.840999-0.008324
0.905037-0.001757
0.953994 0.000043
0.986047-0.000044
1.000000-0.000935
G01G04G06G07G08G10G13
1.00 1.13 1.09 1.09 1.09 1.13 1.00

.: Описание параметров исходного файла:

строка 1
название профиля
строка 3
MACH число Маха
ALFA угол атаки
HEE угол стреловидности
строка 5
BETA диагональное доминирование
ERM максимальная невязка
CLU (0.04) неизменяемая часть профиля в окрестности передней кромки
строка 7
FNPA число узлов сетки на поверхности профиля
FNCO число узлов сетки по нормали к профилю
FNBH число узлов сетки в следе
FKGR (1,2 или 4) параметр измельчения сетки
строка 9
FNLN число линейных итераций
FMGT размерность GMRES
FICY icy=1 - коррекция Сpd для сохранения Су
FIDG idg=0 - проектирование верхней и нижней поверхностей,
idg=1 - проектирование верхней поверхности и сохранение распределения толщины,
idg=2 - проектирование верхней поверхности и сохранение максимальной толщины
строка 11
FNF определяет способ задания геометрии.

FNF=0 - геометрия профиля задана последовательно: задняя кромка, верхняя поверхность, передняя кромка, нижняя поверхность, задняя кромка.

FNF=1 - геометрия профиля задана от передней кромки к задней для верхней и нижней поверхностей, Х - координаты должны быть равны для обеих поверхностей.
FNSP количество координат задания геометрии профиля
строка 13...(13+NSP-1)
X Х - координаты (1:NSP)
YB Y - координаты верхней поверхности (1:NSP)
YH Y - координаты нижней поверхности (1:NSP)
строка 13+NSP+1
G01 Hj+1/Hj в следовой части (верхняя часть сетки)
G04 Hj+1/Hj у задней кромки (верхняя часть сетки)
G06 Hj+1/Hj в средней части профиля (верхняя часть сетки)
G07 Hj+1/Hj в области передней кромки
G08 Hj+1/Hj в средней части профиля (нижняя часть сетки)
G10 Hj+1/Hj у задней кромки (нижняя часть сетки)
G13 Hj+1/Hj в следовой части (нижняя часть)
Количество скачиваний: