Экспериментально-теоретическое исследование обтекания сферы с учетом вдува газа с ее поверхности
Представлены новые экспериментально-теоретические данные по исследованию обтекания твердой сферы в условиях истечения потока массы с ее поверхности. Предложен новый способ определения коэффициента аэродинамического сопротивления твердой сферы при истечении потока воздуха с ее поверхности, включающий измерение аэродинамической силы, действующей на обдуваемую равномерным газовым потоком полую сферу с пористой оболочкой, при подаче во внутреннюю полость сферы сжатого газа под давлением.
Experimental and theoretical studies of the flow around a sphere with account for gas injection from its surface.pdf Введение Обтекание сферической частицы дисперсной фазы (твердой или жидкой) газовым потоком реализуется в различных технических системах и технологических процессах. Исследованию влияния истечения потока газа с поверхности твердой сферы посвящены в основном расчетно-теоретические работы [1-9]. В частности, в работах J.K. Dukowicz [2, 3] проведено численное моделирование обтекания твердой сферической частицы с равномерным истечением потока массы с ее поверхности при малых числах Рейнольдса (Re < 1) и получены теоретические зависимости коэффициента аэродинамического сопротивления от Re обтекающего потока и параметра вдува, характеризующего скорость оттекающего от поверхности частицы потока газа. В работе T.R. Jayawickrama и соавт. [4] проведено теоретическое исследование влияния истечения потока массы с поверхности частицы на коэффициент аэродинамического сопротивления при ее движении в однородном потоке в изотермических условиях в диапазоне чисел Рейнольдса Re = 0.2 ^ 14. Такая же задача решалась в диапазоне чисел Рейнольдса Re =1 ^ 200 численным методом конечных разностей для уравнения Навье-Стокса в работах [5, 6]. В работе M. Watanabe, J. Yahagi [7] представлены результаты численного моделирования динамики движения частицы горящего топлива 58 Архипов В.А., Басалаев С.А., Костюшин К.В. и др. Экспериментально-теоретическое исследование с неоднородным вдувом газового потока с поверхности в диапазоне чисел Рейнольдса Re = 1 + 200. Исследовано изменение двух составляющих коэффициента аэродинамического сопротивления (коэффициенты давления и трения) при разном направлении истечения газового потока с поверхности частицы. В большинстве рассмотренных работ показано снижение коэффициента сопротивления при увеличении скорости вдуваемого потока газа с поверхности частицы. Экспериментальному исследованию обтекания сфероида в условиях вдува газа с его поверхности посвящено ограниченное число публикаций [8, 9]. Анализ известных публикаций по данной тематике показал, что в большинстве работ исследования проводились, как правило, при ускорении частицы, что не позволяет выделить влияние только вдува газа с поверхности частицы на коэффициент аэродинамического сопротивления. При этом данные разных авторов имеют большой разброс (иногда противоречат друг другу). Поэтому для практических расчетов известные зависимости по влиянию эффекта вдува на коэффициент сопротивления требуют дальнейшего тщательного экспериментально-теоретического подтверждения и уточнения. Цель данной работы - комплексное (экспериментальное и теоретическое) исследование обтекания твердой сферы в условиях оттока массы от ее поверхности и получение зависимости для силы аэродинамического сопротивления. Влияние силы, действующей на твердую сферу при обтекании в автомодельном и переходном режимах, экспериментально исследовано в работах [10, 11]. Для повышения чувствительности метода измерения коэффициента сопротивления в переходном режиме разработана экспериментальная установка, включающая рычаг с закрепленной на нем твердой сферой [12]. Экспериментальная установка и методика исследования Схема экспериментальной установки для исследования влияния истечения потока газа с поверхности одиночной твердой сферы на коэффициент аэродинамического сопротивления приведена на рис. 1 [12]. Твердая полая перфорированная сфера 1 закреплялась на длинном плече 2 рычага. Рычаг устанавливался на неподвижной точке опоры 3 с возможностью его вращения вокруг горизонтальной оси 4. На коротком плече 5 рычага крепился стержень 6, контактирующий с приемной площадкой 7 датчика силы 8, в качестве которого использовались высокоточные электронные аналитические весы марки AND GX-200 с погрешностью ± 0.002 гс (1 гс = 9.80665 мН). В ходе эксперимента предварительно на коротком плече 5 рычага размещался дополнительный груз 9 для уравновешивания моментов силы тяжести, приложенных к плечам рычага. При открытии запорного вентиля 10 во внутреннюю полость твердой сферы 1 подавался через тонкостенную трубку 11 сжатый воздух из батареи баллонов 12 через расходомер 13, давление перед которым регистрировалось манометром 14. После открытия запорного вентиля 15 сжатый воздух из батареи баллонов 12 через расходомер 16 и гибкий шланг 17 поступал внутрь патрубка 18, расположенного под полой твердой сферой 1. Давление перед расходомером регистрировалось манометром 19. В качестве расходомеров 13, 16 использовались ротаметры марки РМ-0.63Г с погрешностью ± 2.5%, которые были отградуированы для объемного расхода. 59 Механика / Mechanics Рис. 1. Схема экспериментальной установки для исследования обтекания твердой сферы в условиях вдува газа с ее поверхности: 1 - полая сфера; 2 - длинное плечо рычага; 3 - точка опоры; 4 - ось; 5 - короткое плечо рычага; 6 - стержень; 7 - приемная площадка; 8 - датчик силы; 9 - груз; 10 - запорный вентиль; 11 - трубка; 12 - батарея баллонов; 13 - расходомер; 14 - манометр; 15 - запорный вентиль; 16 - расходомер; 17 - гибкий шланг; 18 - патрубок; 19 - манометр Fig. 1. Diagram of an experimental setup for studying the flow around a solid sphere under conditions of gas injection from its surface: (1) hollow sphere; (2) long lever arm; (3) fulcrum; (4) axis; (5) short lever arm; (6) pin; (7) receiving platform; (8) force sensor; (9) block; (10) shutoff valve; (11) tube; (12) bank of air flasks; (13) flow meter; (14) manometer; (15) shut-off valve; (16) flow meter; (17) flexible hose; (18) branch pipe; and (19) manometer Закрепление твердой полой сферы на длинном плече рычага и стержня, контактирующего с приемной площадкой датчика силы, на коротком плече рычага обеспечивает увеличение измеренной датчиком силы F в n раз по сравнению с аэродинамической силой Fa, действующей на обдуваемую равномерным газовым потоком сферу: F = Fa-n, где n = l\\/l2 > 1 - передаточное отношение рычага, равное отношению длин плеч. Предварительное уравновешивание моментов силы тяжести, приложенных к плечам рычага, путем помещения дополнительного груза на коротком плече рычага обеспечивает исключение влияния массы твердой сферы, тонкостенной трубки для подачи сжатого газа во внутреннюю полость сферы, стержня, контактирующего с приемной площадкой датчика силы и другими элементами рычага, на результаты измерений. При этом датчик силы регистрирует только усиленную в n раз аэродинамическую силу, действующую на твердую полую сферу. Предложенный способ определения влияния истечения потока газа с поверхности твердой перфорированной сферы на коэффициент аэродинамического сопротивления позволяет повысить точность измерения аэродинамической силы при малых числах Рейнольдса, соответствующих малым скоростям обдува сферы. Для получения достоверных данных о влиянии истечения потока массы через перфорированную поверхность твердой сферы в экспериментах измерены основные параметры, описывающие рассматриваемый процесс: диаметр твердой сферы D; скорость обдувающего потока и; скорость вдуваемого с поверхности твердой сферы потока us. 60 Архипов В.А., Басалаев С.А., Костюшин К.В. и др. Экспериментально-теоретическое исследование В экспериментах использовались перфорированные твердые сферы, изготовленные на 3.0-принтере. Диаметр используемых в экспериментах сфер измерялся микрометром с погрешностью 0.01 мм. Эксперименты проведены для сфер диаметрами D = 10 и D = 20 мм; количество перфораций диаметром 1 мм на поверхностях сфер n = 11 и n = 44 соответственно. Длина тонкостенной трубки составляла 65 см. Расстояние от центра сферы до точки опоры рычага l1 = 44 см, а расстояние от стержня, контактирующего с приемной площадкой датчика силы, до точки опоры рычага l2 = 7 см. Скорость обдувающего и истекающего потока воздуха с поверхности твердой сферы измерялась с помощью цифрового анемометра АП-1 с погрешностью ± 0.1 м/c. В экспериментах значения скорости истекающего с поверхности потока воздуха варьировали в диапазоне us = 0.5 ч 1.7 м/с, а значения обдувающего потока воздуха составляли и = 0.7ч1.4 м/с. Экспериментальное значение коэффициента аэродинамического сопротивления от числа Рейнольдса обдувающего потока и плотности потока вдуваемого газа определялось из соотношения F (Re, qs ) F0 (Re) , Со (Re, qs ) = CD0 (Re) где в качестве CD0 (Re) использовалась стандартная кривая сопротивления [13] 24 Re Re < 1, Re = (1 ^ 700), CD0 (Re) = " 24 4 - + ^=, Re 0.44 Re > (700 ^ 3 -105), и число Рейнольдса Re определялось из соотношения Re =S!°, Ц где .F(Re, qS) - сила, измеренная датчиком при заданных значениях Re, qs; Fo(Re) -сила, измеренная датчиком при заданном значение Re при отсутствии вдува с поверхности сферы; D - диаметр полой сферы; р, ц - плотность и коэффициент динамической вязкости газа, обдувающего сферу; и - скорость обдувающего сферу газового потока. При численном моделировании плотность потока газа, равномерно вдуваемого через пористую оболочку, определялась по формуле qs =Ps Q_ nD2 (1) где ps - плотность газа, вдуваемого с поверхности сферы; Qs - объемный расход газа. В экспериментах величина qs рассчитывалась из выражения (2) qs =ps а S ’ где S - суммарная площадь перфораций. 61 Механика / Mechanics Математическая постановка задачи и метод решения Для математического описания задачи обтекания твердой сферы при наличии газоприхода с ее поверхности использовалась система уравнений Навье-Стокса, осредненная по Фавру [14, 15] в декартовой системе координат хг, x2, х3. Уравнение неразрывности имеет вид: (3) (4) Эр д г- -| - +- р и = О, dt дх,1 іА где р - плотность; t - время; Uj - скорость; Xj - координата. Уравнения импульса имеют вид: 3(рм,.) д г_ dt где i - 1, 2, 3; p - давление; 5у - символ Кронекера. Уравнение энергии запишется в виде: д(рЁ) 5 г- " - 1 ■+L puJE+ujP- иіху \\=0: dt (5) где Е - полная внутренняя энергия. Система уравнений (4), (5) замыкается уравнением состояния идеального газа р = рRT, где Т - температура; R - газовая постоянная. Тензор вязких напряжений имеет следующий вид: ( дйі дх, V J дй] дх. 2 ди, -р-- 3 дх,г 8„. -ри, и где р - коэффициент динамической вязкости. Для замыкания системы уравнений вязкого газа использовалась модель турбулентности SST к-щ [15]. Расчеты проведены для двух модельных вариантов подачи газового потока с поверхности твердой сферы (рис. 2). Для первого варианта (см. рис. 2, а), соответствующего физической постановке рассматриваемой задачи, равномерный вдув газового потока задавался по всей поверхности пористой сферы; плотность потока газа определялась формулой (1). Для второго варианта (см. рис. 2, b) газ подавался через трубку в полость твердой сферы с перфорированной поверхностью, что соответствовало условиям проведения экспериментов; плотность потока газа при этом определялась формулой (2). На рис. 2 представлены виды расчетной области и обозначения границ для постановки граничных условий. В качестве граничных условий на границе 1 задавались вектор скорости обтекающего твердую сферу потока и = uinn (где n - вектор внутренней нормали к границе 1), давление и температура окружающей среды: P = 1 атм, T = 296 К, кинетическая энергия турбулентности и диссипация кинетической энергии турбулентности: к = к0, ю=ю0 . Для границ 2 и 3 использовались мягкие граничные условия: d / dn = 0 . На поверхности вдува (граница 4) задавались газоприход G = Gvd 62 Архипов В.А., Басалаев С.А., Костюшин К.В. и др. Экспериментально-теоретическое исследование и полная энтальпия H = H0. Направление вдува - по нормали к поверхности. Для конфигурации со вдувом через отверстия на всех непроницаемых поверхностях твердой сферы и трубки (граница 5) для скорости задавались условия непротекания и прилипания: и = [0 0 0]г , для температуры - условия теплоизолированности: grad(T) = 0, для давления: grad(P) = 0. Для постановки граничных условий для характеристик турбулентности на стенке используется метод пристеночных функций. В качестве начальных условий во всей области течения задавались параметры набегающего потока. Газ - воздух с параметрами: к = 1.4, R = 287 Дж/(кгК), динамическая вязкость: 1.8105 Пас. Рис. 2. Расчетная область и обозначения граничных условий: а - вариант с равномерным по пористой поверхности твердой сферы вдувом; b - вариант с вдувом через перфорированные отверстия на поверхности твердой сферы Fig. 2. Computational domain and notations of boundary conditions: (a) a case with uniform injection on the porous surface of a solid sphere and (b) a case with injection through perforated holes on the solid sphere surface Численный расчет проводился методом установления. Для повышения порядка точности по пространству использовался метод кусочно-линейной реконструкции решения, с ограничителями Venkatakrishnan [16]. Для определения параметров на гранях расчетных ячеек использовался метод Roe с энтропийной коррекцией [17]. Результаты экспериментального и численного исследования Результаты экспериментальных исследований получены в диапазоне чисел Рейнольдса обдувающего потока Re = 133 ^ 1811, что соответствует переходному и автомодельному режимам. Качественная картина истечения газового потока из перфорированных отверстий на поверхности твердой сферы при отсутствии и в присутствии обдувающего потока (и = 0.2 м/с) показана на рис. 3. Экспериментальные значения коэффициента аэродинамического сопротивления твердой сферы CD при истечении газового потока с ее поверхности в зависимости от значения плотности потока газа qs приведены в табл. 1-3. Эксперименты проведены для двух значений qs: qs = 0.60 кг/(м2 с) и qs = 1.09 кг/(м2 с). В табл. 1 -3 приведены измеренные значения скорости обдувающего потока газа (воздуха) и, рассчитанные по формуле (3) значения плотности потока вдува-63 Механика / Mechanics емого газа qs, значения числа Рейнольдса Re, значения коэффициента аэродинамического сопротивления от числа Рейнольдса Ст (Re), рассчитанные по стандартной зависимости (2) и измеренные значения коэффициент аэродинамического сопротивления твердой сферы CD (Re, qs), определенные по формуле (1). Из приведенных в табл. 1 и 2 результатов следует, что при вдуве газа с поверхности сферы в исследованном диапазоне чисел Re = 133^667 коэффициент аэродинамического сопротивления уменьшается. По мере увеличения плотности вдуваемого потока газа при постоянном значении числа Рейнольдса обтекающего потока наблюдается наибольшее изменение (уменьшение) коэффициента аэродинамического сопротивления сферы. а b Рис. 3. Фотографии процесса истечения воздушного потока из перфорированных отверстий на поверхности твердой сферы (D = 20 мм) при отсутствии (а) и в присутствии (b) обдувающего потока (и = 0.2 м/с) Fig. 3. Pictures of the air flow injection from perforated holes on the solid sphere surface (D = 20 mm): (a) in the absence and (b) in the presence of a blowing flow (u = 0.2 m/s) Результаты измерений коэффициента аэродинамического сопротивления сферы D = 10 мм при qs = 0.60 кг/(м2-с) Таблица 1 Параметр Значения и, м/с 0.20 0.30 0.40 0.50 0.60 0.80 1.00 Re 133 200 267 333 400 533 667 Cd0 0.96 0.80 0.71 0.65 0.60 0.54 0.49 Cd 0.87 0.72 0.71 0.63 0.60 0.52 0.48 Результаты измерений коэффициента аэродинамического сопротивления сферы D = 10 мм при qs = 1.09 кг/(м2-с) Таблица 2 Параметр Значения и, м/с 0.20 0.30 0.40 0.50 0.60 0.80 1.00 Re 133 200 267 333 400 533 667 Cd0 0.96 0.80 0.71 0.65 0.60 0.54 0.49 Cd 0.98 0.65 0.62 0.56 0.60 0.51 0.49 64 Архипов В.А., Басалаев С.А., Костюшин К.В. и др. Экспериментально-теоретическое исследование Таблица 3 Результаты измерений коэффициента аэродинамического сопротивления сферы D = 20 мм при qs = 0.60 кг/(м2-с) Параметр Значения u, м/с 0.20 0.30 0.40 0.50 0.60 0.80 1.00 Re 267 400 533 667 800 1066 1333 Cd0 0.71 0.60 0.54 0.49 0.44 0.44 0.44 Cd 0.61 0.56 0.53 0.48 0.44 0.43 0.46 Из табл. 3 следует, что по мере увеличения числа Re обтекающего твердую сферу потока газа коэффициент аэродинамического сопротивления при истечении газа с ее поверхности в промежуточном режиме уменьшается до некоторого переходного значения, соответствующего началу автомодельного режима. При наступлении автомодельного режима (при Re > 700) коэффициент аэродинамического сопротивления соответствует стандартному значению Cm = const = 0.44. Получены результаты численного моделирования обтекания твердой сферы диаметром 10 и 20 мм для случая равномерного вдува воздуха с ее поверхности и для случая вдува потока воздуха через перфорированные отверстия на поверхности твердой сферы. В расчетах варьировали скорость обтекающего твердую сферу потока и газоприход с поверхности сферы в соответствии со значениями параметров, используемых в экспериментах. Расчеты проведены в трехмерной постановке на неструктурированной расчетной сетке с числом ячеек ~ 950 тыс. На рис. 4 показана локальная картина течения для двух рассмотренных расчетных вариантов: с вдувом газа через перфорированные отверстия на поверхности твердой сферы (см. рис. 4, а) и с равномерным вдувом газа (см. рис. 4, b). При вдуве газового потока в полость сферы реализуется струйная картина течения, которая приводит к сильному возмущению течения в окрестности сферы. При равномерном вдуве газа наблюдается лишь небольшое отклонение линий тока набегающего потока. 3.00 2.70 2.40 2.10 a b 1.80 1.50 1.20 0.90 0.60 0.30 0.00 Рис. 4. Поле модуля вектора скорости \\Щ в окрестности твердой сферы (D = 20 мм) при Re = 1 120, м/с: а - расчет с вдувом воздуха через перфорированные отверстия на поверхности сферы; b - расчет с равномерным вдувом воздуха Fig. 4. The field of the velocity vector modulus in the vicinity of a solid sphere (D = 20 mm) at the Reynolds number Re = 1120, m/s: (a) calculation with air injection through perforated holes on the sphere surface and (b) calculation with uniform gas injection 65 Механика / Mechanics Локальная картина поля скоростей в окрестности твердой сферы диаметром 10 мм при равномерном вдуве потока воздуха с ее поверхности при значении числа Рейнольдса Re = 506 показана на рис. 5. По мере увеличения газоприхода с поверхности твердой сферы увеличивается область торможения набегающего потока и существенно расширяется зона рециркуляционного следа. На рис. 5 изображены изолинии скорости, в каждой точке которых величина модуля скорости сохраняет постоянное значение. Рис. 5. Поле модуля вектора скорости |u| в окрестности твердой сферы (D = 10 мм) при Re = 506 (и = 0.8 м/с), м/с: а - us = 0.5 м/с; b - us = 1.7 м/с Fig. 5. The field of the velocity vector modulus |u | in the vicinity of a solid sphere (D = 10 mm) at the Reynolds number Re = 506 (u = 0.8 m/s), m/s: us = (а) 0.5 and (b) 1.7 m/s 0.81 0.77 0.73 0.69 0.64 0.60 0.56 0.51 0.47 0.43 0.38 0.34 0.30 0.25 0.21 0.17 0.13 0.08 0.04 0.00 Распределение давления по поверхности твердой сферы диаметром 10 мм при различных числах Рейнольдса при фиксированной скорости вдуваемого потока воздуха приведено на рис. 6 (ф - зенитный угол в сферической системе координат). По мере увеличения скорости обтекающего твердую сферу потока уровень давления на поверхности сферы существенно увеличивается, но область повышенного давления остается неизменной. На рис. 7 показано распределение давления по поверхности сферы диаметром 10 мм при значении числа Рейнольдса Re = 905. Видно, что при наличии вдува с поверхности твердой сферы, при увеличении газоприхода, максимальное значение давления изменяется незначительно. Стоит отметить, что по сравнению с результатами расчетов, полученными при отсутствии вдуваемого потока воздуха (us = 0), площадь области повышенного давления увеличивается на 10-15%. Эти эффекты наблюдаются для всех рассмотренных расчетных вариантов. На рис. 8 приведены зависимости коэффициента сопротивления твердой сферы диаметром 20 мм от скорости вдуваемого потока воздуха us, рассчитанные для случая равномерного истечения газа с поверхности сферы (см. рис. 8, а) и при вдуве через перфорированные отверстия (см. рис. 8, b) при разных числах Рейнольдса. Отличия в значениях CD на рис. 8, а, б при скорости вдува us = 0 объяс-66 Архипов В.А., Басалаев С.А., Костюшин К.В. и др. Экспериментально-теоретическое исследование няется различием структуры поверхности моделируемой твердой сферы. Результаты расчета показывают, что при равномерном вдуве по мере увеличения скорости вдува воздуха с поверхности твердой сферы наблюдается постепенное снижение коэффициента сопротивления, а при значении us > 1.3 м/с происходит монотонное увеличение CD относительно условий без вдува. Этот эффект связан, по-видимому, с изменением соотношения составляющих коэффициента полного сопротивления (коэффициентов сопротивления давления и трения) при изменении скорости вдуваемого потока [1]. Рис. 6. Расчетное распределение давления по поверхности твердой сферы (D = 10 мм) при скорости вдуваемого потока воздуха us = 0.9 м/с: 1 - Re = 905; 2 - Re = 766; 3 - Re = 633; 4 - Re = 559; 5 - Re = 506; 6 - Re = 453 Fig. 6. Calculated pressure distribution on the solid sphere surface (D = 10 mm) at the blowing air flow velocity us = 0.9 m/s: Re = (1) 905; (2) 766; (3) 633; (4) 559; (5) 506; and (6) 453 Рис. 7. Расчетное распределение давления по поверхности твердой сферы (D = 10 мм) при Re = 905: 1 - us = 0; 2 - us = 1.7 м/с Fig. 7. Calculated pressure distribution on the solid sphere surface (D = 10 mm) at the Reynolds number Re=905: us = (1) 0 and (2) 1.7 m/s 67 1.2 b Рис. 8. Расчетные зависимости коэффициента сопротивления твердой сферы (D = 20 мм) от скорости вдува при равномерном вдуве с ее поверхности (а) и вдувом через перфорированные отверстия (b): 1 - Re = 1 266; 2 - Re = 1 811 Fig. 8. Calculated drag coefficient of a solid sphere (D = 20 mm) as a function of blowing velocity at (a) uniform injection from the sphere surface and (b) injection through perforated holes: Re = (1) 1266 and (2) 1811 Механика / Mechanics a При вдуве воздуха через перфорированную поверхность сферы наблюдается монотонное уменьшение CD с увеличением скорости вдува. Это объясняется отличием интеграла сил давления по поверхности твердой сферы в области торможения потока при вдуве воздуха и его отсутствии. Отличие вызвано влиянием струйного течения при вдуве воздуха через перфорированные отверстия. Рис. 9. Зависимость коэффициента сопротивления твердой сферы (D = 20 мм) от скорости вдува при Re = 1 120 (точки - экспериментальные данные, сплошная линия - численный расчет) Fig. 9. Drag coefficient as a function of blowing velocity for a solid sphere (D = 20 mm) at Re = 1120 (the dotted line indicates experimental data; the solid line, the numerical calculation) На рис. 9 приведены расчетная и экспериментальная зависимости коэффициента сопротивления CD твердой сферы диаметром 20 мм от скорости вдува через ее перфорированную поверхность. Из приведенных на рис. 9 данных следует их количественное и качественное соответствие. 68 Архипов В.А., Басалаев С.А., Костюшин К.В. и др. Экспериментально-теоретическое исследование Заключение - Предложены новая экспериментальная установка и методика исследования влияния истечения воздуха через перфорированную поверхность полой твердой сферы на коэффициент аэродинамического сопротивления в обдувающем потоке. - В исследованном диапазоне чисел Рейнольдса обдувающего потока Re = 133 ^ 1 811 экспериментально показано, что при истечении воздуха через перфорированную поверхность твердой сферы коэффициент аэродинамического сопротивления CD уменьшается. - Численным решением системы уравнений Навье-Стокса получена картина обтекания твердой сферы при равномерном вдуве газа с ее поверхности и при вдуве газа через перфорированную поверхность сферы. Показано, что при равномерном вдуве наблюдается небольшое отклонение лини тока набегающего потока; при вдуве через перфорированную поверхность реализуется струйная картина обтекания, которая приводит к сильному возмущению течения в окрестности сферы. - Проведены численные расчеты коэффициента сопротивления сферы при равномерном вдуве газа с ее поверхности и при вдуве через перфорированные отверстия. Результаты расчетов показали, что при равномерном вдуве по мере увеличения скорости вдува наблюдается монотонное снижение CD, а при значении скорости вдува более 1.3 м/с - увеличение CD по сравнению с его значением при в отсутствие вдува. Этот эффект связан, по-видимому, с изменением соотношения составляющих коэффициента полного сопротивления (коэффициентов сопротивления давления и трения) при изменении скорости вдуваемого потока. - Сравнение расчетных и экспериментальных зависимостей CD(us) при вдуве газа через перфорированную поверхность показало их количественное и качественное соответствие.
Скачать электронную версию публикации
Загружен, раз: 53
Ключевые слова
твердая сфера, обтекание сферы, вдув газа, коэффициент сопротивления, экспериментальное исследование, численное моделированиеАвторы
ФИО | Организация | Дополнительно | |
Архипов Владимир Афанасьевич | Томский государственный университет | доктор физико-математических наук, профессор, заведующий отделом газовой динамики и физики взрыва Научно-исследовательского института прикладной математики и механики | leva@niipmm.tsu.ru |
Басалаев Сергей Александрович | Томский государственный университет | кандидат физико-математических наук, инженер-исследователь Научно-исследовательского института прикладной математики и механики | tarm@niipmm.tsu.ru |
Костюшин Кирилл Владимирович | Томский государственный университет | младший научный сотрудник Научно-исследовательского института прикладной математики и механики | kostushink@niipmm.tsu.ru |
Перфильева Ксения Григорьевна | Томский государственный университет | кандидат физико-математических наук, инженер-исследователь Научно-исследовательского института прикладной математики и механики | k.g.perfiljeva@yandex.ru |
Усанина Анна Сергеевна | Томский государственный университет | кандидат физико-математических наук, инженер-исследователь Научно-исследовательского института прикладной математики и механики | usaninaanna@mail.ru |
Ссылки

Экспериментально-теоретическое исследование обтекания сферы с учетом вдува газа с ее поверхности | Вестник Томского государственного университета. Математика и механика. 2023. № 81. DOI: 10.17223/19988621/81/6
Скачать полнотекстовую версию
Загружен, раз: 147