1 Численные методы Автор курса лекций: Породнов Борис Трифонович, д. ф.-м. н., профессор кафедры молекулярной физики УГТУ-УПИ Екатеринбург 2008 решения. - презентация

1 1 Численные методы Автор курса лекций: Породнов Борис Трифонович, д. ф.-м. н., профессор кафедры молекулярной физики УГТУ-УПИ Екатеринбург 2008 решения задач механики сплошных сред 4. Метод подобия и анализа размерностей

2 2 Лекция 17 Аналитические коэффициенты сопротивления. Анализ размерностей физических величин Численные методы…Лекция 17

3 2 3 Цели изучения: Вычисление аналитически аэродинамических сил и коэффициентов сопротивления при движении вязкой несжимаемой среды в круглой цилиндрической трубе, а также при медленном обтекании шара. Метод анализа размерностей физических величин, позволяющий в простейших случаях установить вид функциональной зависимости характеристик движения от его параметров Численные методы…Лекция 17

4 4 Содержание 9.5. Аналитические коэффициенты сопротивления 9.6. Аэродинамические трубы. Бассейны 9.7. Анализ размерностей физических величин Основные и производные величины Определяющие параметры Период колебаний маятника Сила сопротивления шара Выбор определяющих параметров и основных величин Численные методы…Лекция 17

5 Аналитические коэффициенты сопротивления Во всех случаях, когда сила сопротивления может быть вычислена аналитически, предоставляется возможность аналитически вычислить и коэффициенты сопротивления. Рассмотрим движение вязкой несжимаемой среды в круглой цилиндрической трубе. Воспользуемся общей формулой (9.3.1) для вычисления коэффициентов сопротивления c x (Re). (9.5.1) В данном случае плотность среды постоянная, поэтому ρ = ρ. В качестве характерной скорости υ можно выбрать среднюю скорость движения по трубе. За характерную площадь S можно принять площадь «трущейся» боковой поверхности трубы радиусом r 0 и длиной l. Как было показано в предыдущем разделе, вязкая сила, действующая на рассматриваемый участок трубы со стороны среды, равна разности давлений на этом участке, умноженной на площадь поперечного сечения трубы, т.е. (9.5.2) Средняя по сечению трубы скорость согласно (8.3.19) равна Подставляя ΔP как функцию в формулу для F x (9.5.2), получим: (9.5.3)

6 6 Коэффициент гидравлического сопротивления Из последнего соотношения видно, что коэффициент сопротивления является, действительно, функцией только числа Рейнольдса. Тогда, используя формулу для силы сопротивления F x (9.5.2), можно определить разность давлений на участке трубы длиной l в следующем виде:Рейнольдса (9.5.4) Принято называть коэффициентом гидравлического сопротивления трубы не c x, а величину λ. Формула (9.5.4) является универсальной и позволяет вычислять сопротивление трубы ΔP на участке длиной l, если известен коэффициент сопротивления λ, при любом характере движения среды в трубе. Коэффициенты сопротивления участков трубы, имеющих некоторые особенности (расширение, сужение, кран, задвижка, поворот и т.д.), измерены экспериментально и вычисляются как функции числа Рейнольдса в трубе по графикам или эмпирическим формулам.

7 7 Коэффициент сопротивления шара при медленном обтекании вязкой несжимаемой средой Нетрудно вычислить коэффициент сопротивления шара при медленном его обтекании вязкой несжимаемой жидкостью. В соответствии с формулой Стокса сила сопротивления шара определяется как:Стокса (9.5.5) Обычно в качестве характерной выбирают площадь максимального поперечного сечения шара плоскостью, перпендикулярной к набегающему потоку. Тогда из (9.5.5) следует: (9.5.6) Аналогично можно вычислить коэффициент сопротивления цилиндра единичной длины и т.д. Следует ещё раз подчеркнуть, что удобство введения некоторых безразмерных коэффициентов сопротивления заключается в том, что они являются функциями только безразмерных критериев подобия в любом режиме движения среды как в ламинарном, так и турбулентном. Они могут быть вычислены аналитически только для некоторых простейших ламинарных движений. Во всех других случаях для их определения необходимо обращаться к опыту, и тогда их зависимость от безразмерных критериев подобия устанавливаются на основании опытных данных, представляемых в виде таблиц, графиков или эмпирических формул.

8 Аэродинамические трубы. Бассейны Опыты на моделях проводят для летательных аппаратов в аэродинамических трубах, а для судов, винтов и других подобных устройств в бассейнах. Современные аэродинамические трубы - это гигантские сооружения, потребляющие мощности в десятки тысяч КВт. Увеличение размеров современных аэродинамических труб и соответственно их механической прочности диктуется следующими обстоятельствами. Во-первых, если бы было возможно создание аэродинамических труб, позволяющих обдувать натурные объекты, то, точность научных прогнозов значительно увеличилась бы, и моделирование во многих случаях могло быть полным. Во-вторых, трудно создать очень маленькую модель, являющуюся в точности геометрически подобной натуре вплоть до шероховатостей поверхности модели, заклепок и других тонких, но немаловажных деталей. В третьих, даже для частичного приближенного подобия необходима плотность, а следовательно, и давление модельного потока газа увеличить в 8 раз. В этом случае трудно обеспечить механическую прочность трубы. Если бы можно было увеличить модель, то не потребовалось бы сильно увеличивать давление в аэродинамической трубе или вместо воздуха использовать другой дорогостоящий газ. Эти три обстоятельства и заставляют идти по линии увеличения габаритов современных аэродинамических труб, которые в настоящее время наряду с крупными ускорителями являются уникальными сооружениями наиболее крупных, технически развитых стран. Вышеуказанные соображения о динамическом подобии справедливы и для турбулентных потоков, чисто математического, достаточно точного анализа которых не существует в настоящее время. В этом случае их моделирование является единственным способом решения практических задач.

9 Анализ размерностей физических величин Основные и производные величины Рассмотрим ещё один метод исследования движений вязкой жидкости, позволяющий в простейших случаях установить вид функциональной зависимости характеристик движения от параметров, определяющих движение, и способствующего, тем самым. рациональной постановке эксперимента. Этот метод называют методом анализа размерностей физических величин. Между всеми физическими величинами существует качественное (сила, скорость, длина и т.д.) и количественное различие (две различных длины, скорости и т.д.). Подразумевается, что каждая физическая величина может быть измерена. Мерой физической величины определяют некоторое число, которое характеризует её в количественном отношением в сравнении с аналогичной величиной, принятой за единицу измерения. Все физические величины можно разбить на две категории - основные и производные. Величины, для которых единицы измерения вводятся при помощи эталонов или специальных опытов, называют основными. Величины, единицы измерения которых зависят от единиц измерения основных величин, называют производными.

10 10 Размерность основных величин Так, в механике за основные величины принимаются - масса, длина, время. Их единицы измерения в системе CИ - килограмм (кг), метр (м), секунда (с). Все остальные производные единицы могут быть выражены через них. Число основных единиц может быть и больше трёх. Так, в задачах гидродинамики, связанных с теплообменом, используются градус (0К) и калория (кал). В магнитной гидродинамике можно ввести ещё некоторые основные единицы. Для простоты рассмотрим лишь три основные величины: массу, длину и время. Обозначим их размерность как М, L, и T,соответственно. Рассмотрим размерность некоторой производной физической величины Р, которая может быть выражена через размерность основных величин. Очевидно, размерность любой произвольной физической величины Р может быть записана всегда в виде степенного одночлена: (9.7.1) Здесь m, l, t - показатели размерностей – это некоторые целые или дробные вещественные числа. Если увеличить единицу измерения массы в α раз, длины - в β раз и времени - в γ раз, то единица измерения величины Р будет больше первоначальной в α m, β l, γ t раз. Физическую величину нулевой размерности называют безразмерной величиной. Безразмерная величина инвариантна по отношению к изменениям масштаба измерения основных величин.

11 Определяющие параметры Пусть некоторая физическая величина А является функцией ряда определяющих параметров (9.7.2) Определяющие параметры a 1, a 2, … a n - это некоторые размерные физические величины., от которых может зависеть величина А, в том числе, и основные величины. Среди определяющих параметров a 1, a 2, … a n могут быть параметры с независимыми и зависимыми размерностями. Размерность определяющего параметра является независимой в данном наборе, если она не может быть представлена в виде степенного одночлена размерностей остальных определяющих параметров. Например, размерности L, LT -1 и ML 2 T -2 - независимы, а размерности L, LT -1 и LT -2 - зависимы, так как L = (LT -1 ) 2 /(LT -2 ). Очевидно, что число параметров с независимыми размерностями не может быть больше числа основных величин (в рассматриваемом случае оно не может быть больше трёх). Равенство (9.7.2), выражающее некоторый физический закон, должно быть справедливо для любых масштабов измерения основных величин. Следовательно, размерность величины А должна быть равна размерности функции f.

12 12 Пример анализа размерностей Рассмотрим сначала простейший пример. Пусть определяющими параметрами являются основные величины, причем их число равно числу основных величин: Размерность величины А известна: (9.7.3) Тогда функциональная зависимость величины А должна иметь единственный вид: (9.7.4) так как только таким образом можно соблюсти размерность величины А. Здесь K - некоторый постоянный безразмерный коэффициент. Очевидно, что можно поступить аналогично и в том случае, когда в качестве определяющих параметров выступают три параметра, имеющие произвольные, но зависимые размерности. В этом случае уравнение размерностей имеет вид: (9.7.5) Снова, приравнивая показатели трёх основных величин правой и левой части уравнения (9.7.5), можно из трёх уравнений определить неизвестные показатели α, β, γ. Тогда функциональная зависимость величины А от определяющих параметров a 1, a 2 и a 3 будет иметь вид: (9.7.6)

13 13 Из полученных уравнений следует: Тогда функциональная искомая зависимость периода колебаний маятника от определяющих параметров согласно уравнению (9.7.6) имеет вид: (9.7.7) Безразмерная постоянная K в (9.7.7), конечно, не может быть определена этим методом. Аналитическое решение задачи дает K = 2π Период колебаний маятника Определим функциональную зависимость периода колебаний математического маятника от определяющих параметров. Предполагается, что определяющими параметрами в данном случае являются l - длина маятника, m – его масса и g - ускорение силы тяжести. Поэтому функциональную зависимость периода колебаний маятника от трёх определяющих параметров с их независимыми размерностями можно записать в виде: Рис. 9.1 Следуя формуле (9.7.5) имеем Приравнивая показатели степени при одинаковых размерностях M, L, и T имеем соответственно:

14 Сила сопротивления шара Если число определяющих параметров больше трёх, то можно выбрать из них три с независимыми размерностями. Тогда размерности остальных определяющих параметров могут быть выражены через независимые размерности выбранных трёх параметров в виде степенных одночленов. Ясно, что в этом случае система уравнений для показателей размерностей становится неопределённой, а, следовательно, исключена возможность установления однозначной искомой функции вида (9.7.6). Очевидно, что в случае трёх определяющих параметров с независимыми размерностями из них невозможно организовать ни одной безразмерной комбинации. Если число определяющих параметров равно четырём, причем, из них только три имеют независимые размерности, то нетрудно видеть, что из имеющихся четырёх параметров можно составить лишь одну одночленную степенную комбинацию, не имеющую размерности. Если имеется пять определяющее параметров, то таких безразмерных комбинаций можно составить только две, если предыдущая безразмерная комбинация рассматривается как некоторый четвертый определяющий параметр и т.д. В этом случае функциональная зависимость физической величины (9.7.5) может быть записана в следующем виде: (9.7.7) Здесь произвольная функция f от независимых безразмерных комбинаций (δ 1, δ 2, … δ n-3 ) остается неопределенной. Показатели степеней должны выбираться при помощи каких-то дополнительных соображений. Сказанное выше лучше всего иллюстрировать примером.

15 15 Уравнение размерностей для шара Найдём функциональную зависимость силы сопротивления обтекаемого вязкой несжимаемой жидкостью шара от определяющих параметров. Очевидно, определяющими параметрами будут: ρ - плотность жидкости, υ - скорость набегающего потока, r 0 - радиус шара, ν - коэффициент кинематической вязкости. Запишем размерности определяемой величины и параметров (9.7.8) В данной задаче число определяющих параметров равно четырём. Из них первые три могут быть приняты за основные, т.к. они имеют независимые размерности. Найдем безразмерную комбинацию параметров. Очевидно, для этого достаточно составить следующее уравнение размерностей: (9.7.9) Приравнивая в этом уравнении степени одинаковых основных величин, получим уравнения для показателей степени: Решение системы уравнений имеет вид: Следовательно, показатели,, могут быть любыми одинаковыми целыми или дробными вещественными числами. Безразмерная комбинация параметров имеет вид:

16 16 Обтекание шара при малых числах Рейнольдса В последней формуле показатель степени остается неопределённым, поэтому можно заменить степенную формулу некоторой функцией от числа Рейнольдса. Тогда общий вид формулы для силы сопротивления шара будет иметь вид: (9.7.10) Однако, поскольку, приравнивая показатели степеней правой и левой части уравнения размерности (9.7.2), получено лишь три уравнения для четырёх показателей, то задача определения показателей становится неопределенной. Можно их определить лишь в двух предельных случаях движения жидкости. Рассмотрим сначала случай движения при малых числах Рейнольдса. Движение жидкости при малых числах Рейнольдса можно рассматривать как движение жидкости, обладающей очень большой вязкостью. В этом случае, очевидно, плотность жидкости не должна играть существенной роли и её можно исключить из числа определяющих параметров. Тогда определяющими параметрами будут υ, r 0, η, имеющие независимые размерности: Уравнение размерностей имеет вид: (9.7.11) Естественно, что никаких безразмерных комбинаций из этих трёх параметров составить нельзя. Для нахождения показателей степеней имеем уравнения: ;

17 17 Обтекание шара при больших числах Рейнольдса Из этой системы уравнений находим: α = 1, β = 1, γ = 1. Следовательно, функциональная зависимость (9.7.3) имеет вид: (9.7.12) Таким образом, получена формула Стокса для силы сопротивления шара с неопределенным безразмерным численным коэффициентом K. Аналитическое решение, как показано ранее, дает K = 6π.Стокса Рассмотрим другой предельный случай обтекания шара - обтекание его при очень больших числах Рейнольдса (Re ). Движение жидкости при очень больших числах Рейнольдса можно рассматривать как движение маловязкой жидкости. Тогда из определяющих параметров можно исключить коэффициент вязкости η. Размерности остальных определяющих параметров ρ, υ, r 0 независимы и из них нельзя образовать некоторую безразмерную комбинацию. Уравнение размерностей имеет вид (9.7.13)Рейнольдса Для показателей степени имеем уравнения: Решение полученной системы уравнений даёт: α = 1, β = 2, γ = 2. Тогда сила сопротивления зависит от определяющих параметров следующим образом: (9.7.14) Таким образом, сила сопротивления шара при больших числах Рейнольдса (в том числе, и при турбулентном движении) должна быть пропорциональна квадрату скорости, а не первой степени, как в предыдущем случае (Re ). Единственный неизвестный безразмерный коэффициент K может быть определен из одного единственного опыта.

18 18 При промежуточных числах Рейнольдса метод размерностей не может однозначно определить функциональную зависимость силы сопротивления от всех определяющих параметров. Однако, и в этом случае метод размерностей может быть весьма полезен. Действительно, при любых числах Рейнольдса уравнение размерностей имеет вид: (9.7.15) Система уравнений для нахождения показателей степеней содержит три уравнения для нахождения четырех неизвестных: Решение системы неоднозначно. Если положить, то имеем: α = 1, β = 1, γ = 1, = 1. Тогда функциональная зависимость силы сопротивления от определяющих параметров имеет вид: (9.7.16) Если положить β = 2, то α = 1, ε = 0, γ = 2 и аналогичное выражение нужно записать в следующем виде: (9.7.17) К сожалению, теория размерностей не может указать, какая из этих формул является правильной. Но во всяком случае имеем дело в данном случае лишь с двумя функциональными зависимостями, и дело экспериментатора определить неизвестные коэффициенты, функциональную зависимость от безразмерных комбинаций определяющих параметров, а также пределы применимости той или иной формулы. Уравнение размерностей при любых числах Рейнольдса

19 19 Если известна математическая постановка задачи, то выбор определяющих параметров следует из уравнений движения и граничных условий. В этом случае определяющие параметры - это все данные, которые необходимо задавать для вычисления искомых физических величин. Если же уравнения, описывающие данное явление, неизвестны, то остается, полагаясь на физическую интуицию и предварительные экспериментальные данные, испытать несколько гипотез о предполагаемой системе определяющих параметров. Здесь были рассмотрены лишь три основных физических величины: масса, длина и время. Поэтому для нахождения показателей степеней имелось всегда только три уравнения и, если число определяющих параметров было больше трёх, это ставило исследователя в затруднительное положение. Введение некоторых дополнительных основных единиц измерения могло бы помочь в некоторых случаях решать задачу однозначно, хотя вряд ли это будет полезно в методологическом отношении. В заключение отметим, что полезность теории размерности заключается: во- первых, в определении точного вида функциональной зависимости (если число независимых определяющих параметров меньше или равно числу основных величин) и, во-вторых, в определении функциональных альтернативных зависимостей (если число определяющих параметров больше, чем число основных величин). Тем самым, она помогает, не решая уравнений движения, правильно установить вид функциональных зависимостей и целенаправленно поставить эксперимент, значительно сокращая объём экспериментальной работы Выбор определяющих параметров и основных величин.

20 20 Выводы Аналитически вычислены коэффициенты сопротивления при движении вязкой несжимаемой среды в круглой цилиндрической трубе, а также коэффициент сопротивления шара при медленном обтекании. Рассмотрен метод анализа размерностей физических величин. Введены понятия основных и производных величин, размерность основных величин, безразмерной величины, определяющих параметров. Проведен анализ размерностей на примерах колебаний математического маятника, а также зависимости силы сопротивления шара, обтекаемого вязкой несжимаемой средой, от определяющих параметров Численные методы…Лекция 17

21 21 Информационное обеспечение лекции Литература по теме: Ландау Л.Д., Лившиц Е.М.. Гидродинамика. М.: Наука с. Седов Л.И. Механика сплошной среды. М.: Наука Т с.; Т.2, 568с. Фабрикант Н.Я. Аэродинамика. М.: ГИТТЛ с. Лойцянский Л.Г. Механика жидкости и газа. М.: Наука с Численные методы…Лекция 17

22 22 Справочные данные Курс лекций является частью учебно-методического комплекса «Численные методы расчета задач механики сплошных сред. 1. Теория упругости и идеальная среда». Автор: Породнов Борис Трифонович, д. ф. – м. н., профессор кафедры молекулярной физики УГТУ-УПИ. Учебно-методический комплекс подготовлен на кафедре МФ ФТФ ГОУ ВПО УГТУ-УПИ. электронный адрес: