CC BY
29
УДК 532.526
ПРИМЕНЕНИЕ СХЕМЫ ИНТЕГРИРОВАНИЯ С РАЗНЫМИ ШАГАМИ
В МЕТОДЕ РАСЩЕПЛЕНИЯ ЗАВИХРЕННОСТИ ПРИ МОДЕЛИРОВАНИИ АЭРОДИНАМИЧЕСКОГО СЛЕДА ЗА ПЛОСКОЙ ПЛАСТИНОЙ
© 2008 В.В. Никонов, В.Г. Шахов Самарский государственный аэрокосмический университет
Рассматривается применение метода расщепления завихрённости к моделированию аэродинамического следа за продольно обтекаемой плоской пластиной. Применяется схема интегрирования с разными шагами по времени с учетом разности скоростей протекания процессов диффузии и конвекции. Показано, что схема метода позволяет получать результаты с хорошей точностью в широком диапазоне изменения чисел Рейнольдса.
Обтекание, прямое численное моделирование, метод расщепления завихрённости, аэродинамический след, плоская пластина, диффузия, конвекция, интегрирование с разными шагами по времени, число Рейнольдса
После расчёта поля скоростей течения новые координаты частиц получаются аналогично методу «вихрь в ячейке» численным интегрированием системы обыкновенных дифференциальных уравнений методом Эйлера. ке» путём расщепления завихрённости на её уТ г
■’ пт Новое местоположение частиц не обязатель-
1. Математическая формулировка метода расщепления завихрённости
Метод расщепления завихрённости (МРЗ) был получен из метода «вихрь в ячей-
составляющие [1]. В МРЗ вместо циркуляции вихрей в ячейках используются ве личины
х1 +к/2 У; +к!2
Дц- (х , У; ) = \ |
xi -к/2 у;-к/2 хI +к/2 У; +к/2
Д« (Х, У;) = 1 1
Эи (х, у) Э У
Э у( х, у) Э х
ах йу, ах йу,
xi-к/2 у;-к/2
где к - размер ячейки однородной сетки.
Скорость течения рассчитывается интегрированием вдоль координатных линий, и при этом двумерная задача сводится к нескольким одномерным:
и(і, і + 0,5) = и (і, і - 0,5) + -
К(і, і)
у(і + °А і) = у(і - 0,5 і)е +
К ’
А ух (і, і )
К
Скорость в центрах ячеек сетки рассчитывается следующим образом:
но совпадёт с координатами расчётной сетки, и поэтому процедура перераспределения используется отдельно для величин Дух и Аиу:
Аиу (X , У; ) = Д иу (X , У1 )Л(X - X )Л(У; - У1 ),
ДУх (X , У; ) = Дvx (X, У1 )Л(X - X )Л(У; - У1). (1)
В качестве интерполяционной функции Л использовалась формула «облако в ячейке» [2].
Диффузия в свободном потоке для схемы МРЗ рассчитывается с использованием метода донор-акцептор (Д-А) аналогично [3], но отдельно для ДиУ и Д^:
ДиУ (* + Д *, X , У; ) = ДиУ (*, X, У; ) +
+ X (ДиУ(* , X, У,) 6*(У) -
]'- £ ]'+пё
-ДиУ (*, X , У; ) 61- (У) ) ,
Д vx (* + Д X , У; ) = Д vx (*, X, У; ) +
и(і, І) = и (і, і - 0,5) + 0,5Д у (і, і) = у (і - 0,5, і) + 0,5Д
(Аі)
XV Кх
<Л і) К
+ X (Дух (І, х9 , Уі ) (х)'
і-п4 £+ПЛ
-Дух 0, х9 , Уі ) С*г(х) ) ,
В схеме метода МРЗ частицы движутся с потоком и переносят величины Дис и Диу.
где щ - радиус «диффузионной молекулы», а коэффициенты 0*ч определяются по формуле
G*Pq (-)
- erf
erf
V4V^
4
v41
Аun (^ , nj ) = u(S, , +0)hs Х
un\i> j /
Х
erf
nj + h /2
4
■v*t
erf
n
-h/2
V4^
А vs (si, «j ) = 0
где кц - размер ячейки локальной сетки в направлении координаты б, (б;, nj) - координаты центра рассматриваемой ячейки в локальной системе координат (СК). Величины Диу и Дух в глобальной СК запишутся как
Д иУ (X , У, ) =Дип (Эг, П; ) Уп ,
А vx (x , у, ) =Аun (si, nj ) xn .
(2)
Здесь z - Xp- или _ур-координаты центра ячейки p-го вихря, erf - интеграл вероятности (функция ошибок).
В схеме МРЗ с поверхности тела диффундируют величины Диу и Дух. Для удобства сначала определим Дип и Dvs в системе координат, связанной с панелью:
Здесь (хі, у]) - координаты центра рассматриваемой ячейки в глобальной СК, (х„, Уп} - единичный вектор, нормальный к рассматриваемой панели, записанный в этой же системе. Величина скорости жидкости у поверхности тела определяется следующим образом:
и(я, +0) = ііш и(я, п) ,
П®+0
В общем случае координаты (хі, у]) в (2) не обязательно совпадают с ячейками глобальной сетки. Поэтому для значений Диу и Дух применяется процедура перераспределения (1).
Для учёта уравнения неразрывности в схему метода необходимо включить корректировку поля скорости. Тогда для расположенных над пластиной ячеек вертикальная компонента скорости определится как
- 05к' 157 )%4+\Щ,;-1} (3>
Для аппроксимации пространственных производных в правой части (3) использовалась центральная конечно-разностная схема.
В схеме метода процессы конвекции и диффузии рассматриваются отдельно на каждом шаге по времени. В работе [4] было показано, что для достижения заданной точности шаг по времени для метода Д-А определяется следующим соотношением:
д * = клк2 / V ,
где ка - константа, зависящая только от радиуса па «диффузионной молекулы», и для па = 1 эта константа находится в диапазоне
0.2 < ка < 0.21. Принято, что ка = 0.21, так как в данном случае [5] ошибки методов Д-А и моделирования процесса конвекции будут иметь разные знаки и поэтому будут компенсировать друг друга. Шаг по времени для процесса конвекции определялся с помощью неравенства к
А *с < ка~ ,
и¥
соответствующего критерию Курант-Фрид-рих-Леви [6] с величиной коэффициента кс = 1.5. При этом он не обязательно совпадёт с оптимальным шагом для расчёта процесса диффузии. Поэтому предлагается применять метод интегрирования с раздельными шагами по времени для процессов диффузии и конвекции.
2. Прямое численное моделирование аэродинамического следа за плоской пластиной
Рассматривается случай продольного обтекания плоской пластины конечной длины вязкой несжимаемой жидкостью. В результате прямого численного моделирования получены профили скорости в аэродинамическом следе за пластиной в диа-пазаоне чисел Рейнольдса от 10 до 10 6, которые сравниваются с аналитическим решением Голдстейна [7] для сечений, расположенных на расстояниях: х1 = 0; 0.0135; 0.108;
0.256 от задней кромки пластины. Полученные результаты представлены на рис. 1, 2.
Вертикальная безразмерная координата Л1 определяется следующим образом:
Лі = У
и¥ V /і
где 11 - длина пластины, равная 1.0, так как она принималась за характерный размер.
Представленные результаты позволяют сделать вывод о том, что метод МРЗ применим не только для моделирования ламинарного пограничного слоя на плоской пластине [1], но и для моделирования аэродинамического следа за пластиной.
Лі
Рис. 1 Распределение продольной компоненты скорости и в следе за плоской пластиной в сравнении с
аналитическим решением (Кв = 10, И = 0.025).
Численное решение: О - х1 » 0.0, - х1 » 0.0135, - х1 » 0.108, X - х1 » 0.256;
_______- Голдстейн [7]
Рис. 2. Распределение продольной компоненты скорости и в следе за плоской пластиной в сравнении с
аналитическим решением (Кв = 106, И = 1.563-10'4).
Численное решение: О - х1 » 0.0, - х1 » 0.0135, - х1 » 0.108, X - х1 » 0.256;
______- Голдстейн [7]
Библиографический список
1. Никонов, В.В. Схема расчета скорости для метода «вихрь в ячейке» применительно к моделированию двумерного ламинарного пограничного слоя [Текст] / В.В. Никонов, В Г. Шахов // Известия СНЦ РАН, Самара. - 2005. - Т. 7, № 2. - С. 392-398.
2. Григорьев, Ю.Н. Численные методы «частицы-в-ячейках» [Текст] / Ю.Н. Григорьев, В.А. Вшивков. - Новосибирск: Наука. - Сибирская издательская фирма РАН, 2000. - 184 с.
3. Taranov, A. Development of the Computational Vortex Method for Calculation of Two-Dimensional Ship Sections with Flow Separation [Текст] / A. Taranov, N. Kornev,
A.Leder // Schiffbauforschung. - 2000. - Vol. 39, N 2. - P. 95-105.
4. Nikonov, V. The Ratio between Spatial and Time Resolutions for the Diffusion Substep in 2D Computational Vortex Methods [Tекст] / V. Nikonov, N. Kornev, A. Leder // Schiffbauforschung. - 2002. - Vol. 41, N 3/4. - P. 5-12.
5. Никонов, В.В. Модификация схемы «донор-акцептор» для расчета диффузии завихренности и ее применение в методе «вихрь в ячейке» [Текст] / В.В. Никонов,
В.Г. Шахов // Вестник СГАУ, Самара. -2003. - № 1 (3). - С. 38-46.
6. Ferziger, J. Computational methods for fluid dynamics [Tекст] / J. Ferziger, M. Peric, 3 rev. ed. - Springer-Verlag, 2002. - 423 p.
7. Шлихтинг, Г. Теория пограничного слоя [Текст] / Г. Шлихтинг; пер. с нем. Г. А. Вольперта; под. общ. ред. Л.Г. Лойцян-ского. - М.: Наука, 1974. - 712 с.
References
1. Nikonov, V.V. Calculation velocity scheme for “Vortex-in-Cell” method applying to simulation of two-dimensional laminar boundary layer / V.V. Nikonov, V.G. Shakhov // Samara: Proceedings of Samara Scientific Center of Russian Academy of Sciences. - 2005. - Vol. 7, N 2. - P. 392-398. - [in Russian].
2. Grigoryev, Yu.N. Numerical methods “Particles-in-cells” / Yu.N. Grigoryev, V.A. Vshivkov. - Novosibirsk: “Nauka” (Science). -Siberian publishing firm of RAS. - 2000. - 184 p. - [in Russian].
3. Taranov, A. Development of the Computational Vortex Method for Calculation of Two-Dimensional Ship Sections with Flow Separation / A. Taranov, N. Kornev, A. Leder //
Schiffbauforschung. - 2000. - Vol. 39, N 2. -P.95-105.
4. Nikonov, V. The Ratio between Spatial and Time Resolutions for the Diffusion Substep in 2D Computational Vortex Methods / V. Nikonov, N. Kornev, A. Leder // Schiffbauforschung. - 2002, vol. 41, N 3/4. - pp. 5-12.
5. Nikonov, V.V. “Donor-Acceptof’ scheme modification for vorticity diffusion calculation and it application in “Vortex-in-Cell” method / V.V. Nikonov, V.G. Shakhov // Vest-nik (Bulletin) SSAU. - Samara, 2003. - N 1 (3).
- P. 38-46. - [in Russian].
6. Ferziger, J. Computational methods for fluid dynamics / J. Ferziger, M. Peric. - 3 rev. ed., Springer-Verlag. - 2002. - 423 p.
7. Schlichting, H. Grenzschicht-Theorie / H. Schlichting - Springer, Berlin, 1997. - 851 p.
MULTI-TIME-STEP SCHEME IN THE VORTICITY SPLITTING METHOD
APPLIED TO THE SIMULATION OF THE WAKE BEHIND THE FLAT PLATE
© 2008 V.V. Nikonov, V.G. Shakhov Samara State Aerospace University
The Vorticity Splitting Method is applied to the simulation of the wake behind the longitudinal streamlined flat plate. Moreover the Multi-Time-Step Scheme is used in the numerical algorithm because the differences in the rates of convection and diffusion processes. The method is allowed results obtaining with fine accuracy in the wide range of Reynolds numbers.
Streamline, direct numerical simulation, vorticity splitting method, wake, flat plate, diffusion, convection, multi-time-step integration, Reynolds number
Сведения об авторах
Никонов Валерий Владимирович, инженер НТП «Авиатехнокон» СГАУ, кандидат технических наук. Е-mail: v [email protected]. Область научных интересов: вихревые методы, прямое численное моделирование несжимаемых и сжимаемых течений, пограничный слой.
Шахов Валентин Гаврилович, заведующий кафедрой аэрогидродинамики СГАУ, профессор, кандидат технических наук. Е-mail: [email protected]. Область научных интересов: теория пограничного слоя, турбулентность, численные методы, аэродинамика летательных аппаратов.
Nikonov Valery Vladimirovich, engineer NTP “Aviatechnokon” SSAU, PhD of technical sciences. Е-mail: v [email protected]. Scientific interests: vortex methods, direct numerical simulation of incompressible and compressible flows, boundary layer.
Shakhov Valentin Gavrilovich, Head of Aerohydrodynamic department SSAU, professor, PhD of technical sciences. Е-mail: [email protected]. Scientific interests: boundary layer theory, turbulence, numerical methods, flying vehicles aerodynamics.
