ГИДРОАЭРОМЕХАНИКА
наука о движении и равновесии жидкостей и газов. При планировании физических экспериментов или при их проведении необходимо создавать теоретические модели, которые либо предсказывают возможные результаты этих экспериментов, либо объясняют уже полученные. Только в тесном взаимодействии теории и эксперимента можно понять то, что происходит в окружающем нас физическом мире. Для создания той или иной количественной или качественной модели физического явления необходим математический фундамент, на основе которого строятся такие модели. Под математическим фундаментом в данном случае подразумеваются те дифференциальные уравнения и те граничные и начальные условия, с помощью которых можно было бы описывать рассматриваемое физическое явление. Гидромеханика и предлагает модели и аппарат для иcследования явлений, происходящих в жидкостях и газах.
О гипотезе сплошности среды.
Гидроаэромеханика изучает движения жидкостей и газов в приближении, когда они могут рассматриваться как сплошные среды, т.е. среды, непрерывным образом заполняющие рассматриваемое пространство течения. Чтобы решать математические проблемы, связанные с расчетом движения различных объектов (самолетов, ракет, кораблей и др.) в воздухе или воде, с изучением волновых процессов в жидкостях и газах, с их течениями по трубам и каналам и т.п., необходим математический аппарат, описывающий эти явления. Этим аппаратом и являются уравнения гидроаэромеханики, которые опираются на гипотезу сплошности среды, т.е. на гипотезу о том, что частицы жидкости или газа непрерывным образом заполняют занимаемую ими часть физического пространства.
Возникает естественный вопрос: при каких предположениях справедлива эта гипотеза? Если для жидкостей (воды, жидких металлов и т.п.) эта гипотеза более или менее очевидна, то для достаточно разреженных газов (например, занимающих космическое пространство, включая атмосферы звезд, планет и Солнца), которые состоят из отдельных атомов или молекул, а также других физических объектов, к которым применим аппарат гидроаэромеханики, она требует своего обоснования. Так, например, при расчете торможения искусственных спутников Земли использование математического аппарата гидроаэромеханики не представляется возможным, в то время как именно этот аппарат используется при расчете торможения космических объектов, входящих в плотные слои атмосфер Земли и планет (например, метеоритов или возвращаемых на Землю космических кораблей и пр.). На этот вопрос легко ответить при выводе уравнений. Однако из этого вывода следует, что гипотеза сплошности среды справедлива, в частности, в том случае, когда характерный размер обтекаемого тела L (например, радиус сферического спутника) много больше длины свободного пробега атомов или молекул газа l, т.е. длины между последовательными их столкновениями. Замкнутая система уравнений гидроаэромеханики. Уравнения гидроаэромеханики в их упрощенном виде представляют собой сложную систему нелинейных дифференциальных уравнений для массовой плотности r (масса жидкости или газа в единице объема), вектора скорости V и давления p, которые, в свою очередь, являются функциями пространственных координат (например, x, y и z в декартовой системе координат) и времени t. Не вдаваясь в математические подробности вывода этих уравнений, можно рассмотреть основные идеи этого вывода, тем более, что эти уравнения представляют собой известные даже из школьных учебников законы сохранения массы, импульса и энергии. Для этого рассматривается некоторый физический объем, непрерывным образом заполненный жидкостью или газом. На рис. 1 изображена движущаяся жидкость (или газ), непрерывным образом заполняющая некоторую часть физического пространства. Выделим из нее некоторый объем U (ограниченный поверхностью S), который в течение всего времени движения состоит из одних и тех же частиц жидкости (этот объем заштрихован).Очевидно, что при своем движении масса жидкости, заключенная в объеме U, остается постоянной (если, конечно, нет каких-либо дополнительных источников этой массы), хотя сам объем может сильно деформироваться, поскольку частицы не скреплены жестко, как в твердом теле. Если выделить из рассматриваемого объема бесконечно малый элемент DU, то очевидно, что в этом элементе масса жидкости или газа будет равна rDU. Тогда закон сохранения массы, заключенной в выделенном объеме U, можно записать в виде![]() ![]() где оператор дивиргенции ![]() ![]() ![]() Закон сохранения массы в интегральной форме справедлив как для непрерывных, так и для разрывных функций r и V. Для непрерывных функций закон сохранения массы можно записать в дифференциальной форме![]() Аналогично можно записать теперь закон сохранения импульса. Импульс единицы объема жидкости, равен rV, в элементарном объеме rDU, а в выделенном объеме U ![]() ![]() ![]() ![]() ![]() ![]() а величины u, v и w, а также ![]() ![]() Это уравнение, называемое уравнением Навье Стокса, выписано в наиболее простой форме для несжимаемой жидкости, где поверхностные силы сводятся к нормальному давлению р, а последний член справа представляет собой «вязкие» силы ( m коэффициент вязкости) в предположении, что r = const. Впервые уравнение движения было выведено в середине 18 в. Л.Эйлером, когда он работал в Петербургской Академии наук. Поскольку эффекты вязкости в жидкости в то время еще не были известны, то Эйлер получил это уравнение при m = 0. В честь его эти уравнения были названы уравнениями Эйлера. Только в 1822 французским инженером Навье в уравнения Эйлера были введены силы, связанные с вязкостью, определяемой коэффициентом m. В общей форме, справедливой и для сжимаемого газа, уравнение получено Стоксом и получило название уравнения Навье Стокса.Для несжимаемой жидкости дифференциальные уравнения неразрывности и импульса (одно скалярное и одно векторное) являются замкнутой системой уравнений для определения вектора скорости V и скалярного давления р (r = const). Если же r № const, то требуется дополнительное уравнение. Это уравнение получается из закона сохранения энергии.Обобщение закона сохранения энергии на случай движения жидкостей и газов получается аналогично обобщению второго закона Ньютона, однако, в силу наличия теплового движения в жидкостях и газах, энергия, приходящаяся на единицу объема, состоит из кинетической энергии rV2/2 и внутренней энергией re, связанной с тепловым движением частиц газа или жидкости. Полная энергия в элементе объема DU равна r(V2/2 + e)DU.Изменение полной энергии в выделенном объеме U равно притоку тепла через поверхность S за счет теплопроводности, а также работе массовых и поверхностных сил, т.е. вместо закона сохранения импульса, получается уравнение![]() ![]() ![]() p = r R T, где R = (ср сv ) газовая постоянная, а ср теплоемкость при постоянном давлении, и законом Фурье![]() ![]() Уравнения гидростатики получаются из уравнений гидроаэромеханики при V = 0. В частности, уравнения сохранения импульса дает Откуда, в частности, следует известный еще из школьных учебников закон Паскаля, согласно которому при отсутствии внешних массовых сил (F = 0) давление всюду является постоянным (p = const ). Равновесие совершенного газа в поле сил тяжести. Пусть есть газ в центральном поле сил тяжести. Уравнения равновесия в сферической системе координат будут в этом случае записываются как:![]() ![]() Можно, например, рассчитать распределения давления в атмосфере Земли до расстояний в 11 км от ее поверхности. Если выбрать декартову систему координат с началом на поверхности Земли и направить ось Oz вертикально вверх, тогда в барометрической формуле вместо координаты r нужно брать координату z = r RЕ, где RЕ радиус Земли. Поскольку этот радиус много больше толщины атмосферы (z << RЕ), то барометрическую формулу для плоской атмосферы можно переписать в виде![]()
![]() где ![]() В ограниченном диапазоне высот до 11 км, важном в практическом отношении (пассажирские самолеты обычно летают на высотах, не превышающих эту высоту, альпинисты поднимаются на вершины, самая высокая из которых Эверест имеет высоту ~ 8,800 м и т.д.) часто аппроксимируют температуру линейной функцией высоты![]() Если принять такое распределение температуры, то давление записывается в виде ![]() p = p0 r gz или р = p0 + r gh, где h глубина жидкости под ее поверхностью, р0 давление на поверхности (рис. 2). Эта формула, известная из школьных учебников, показывает, как давление в жидкости возрастает с ее глубиной. С помощью этой формулы легко рассчитать давление на дно сосуда, заполненного жидкостью. Интересно, что это давление зависит от глубины, но не зависит от формы сосуда. В частности, на рис. 3 давление на дно сосудов 1 и 2 одинаковой площади дна S будет одинаковым или сила, действующая на дно этих сосудов вследствие давления жидкостей, будет одинаковой. Много важных приложений основывается на решениях уравнений гидростатики (закон Архимеда, устойчивость равновесия атмосфер звезд и планет и т.п.). Уравнения гидроаэромеханики для вязких и теплопроводных жидкостей или газов в большинстве очень важных для практики проблем поддаются решению только численными методами. Однако эти уравнения существенно упрощаются в предположении, что для рассматриваемого течения справедливо предположение о его несжимаемости ( r = const). Хотя строго несжимаемых жидкостей или газов в природе не существует, тем не менее, во многих случаях, например, сжимаемый газ можно рассматривать как несжимаемую жидкость, поскольку изменением плотности во многих течениях можно пренебречь. При этом уравнение неразрывности для несжимаемой жидкости принимает вид div = 0. Вместе с уравнением сохранения импульса оно образует замкнутую систему уравнений для определения давления р и скорости V. Два критерия определяют возможность использования модели несжимаемой жидкости для, вообще говоря, сжимаемого газа![]() ![]() где u скорость жидкости вдоль оси х, совпадающей с осью симметрии трубы, а r расстояние от этой оси. Из этой видно, что профиль скорости в трубе является параболическим. На стенках трубы скорость обращается в нуль вследствие прилипания жидкости из-за эффекта вязкости. Такое течение было изучено в середине 19 в. Пуазейлем и Гагеном на примере течений жидкостей в капиллярах и получило название течения Гагена Пуазейля. Очевидно, при постоянном потоке (не зависящем от r) жидкости у входа в трубу и на ее начальном участке профиль скорости не будет совпадать с приведенным решением. Параболический профиль устанавливается лишь на достаточно большом расстоянии от входного участка, именно поэтому для получения решения нужно предположить, что труба достаточно длинная, при этом для таких труб это точное решение хорошо совпадает с экспериментальными данными.Полученное решение описывает стационарное, гладко-слоистое течение, которое обычно называют ламинарным. Однако из практики известно, что в трубах иногда течение бывает нестационарным, с пульсациями скорости, с перемешиванием между слоями, это течение обычно называется турбулентным. Опыты Рейнольдса, проведенные в 1883, показали, что при достаточно больших значениях числа r U L/m, где U средняя по сечению трубы скорость жидкости, параболический профиль становится неустойчивым по отношению к малым возмущениям, а при дальнейшем увеличении этого числа течение в трубе становится турбулентным. Это число получило название числа Рейнольдса (Re), которое играет очень важную роль в различных задачах гидроаэромеханики. В частности оно характеризует отношение инерционных сил (левая часть уравнения) к силам вязкости, при этом часто силами вязкости можно пренебречь и использовать уравнения гидроаэромеханики идеальной жидкоститолько при Re >> 1. Течения идеальных жидкостей и газов. Часто важные в приложениях задачи рассматривают на основе уравнений гидроаэромеханики идеальной жидкости, а не на полных уравнениях. Это связано с тем, что математически уравнения идеальной гидроаэромеханики существенно проще. Если нужно определить подъемную силу крыла самолета при малых дозвуковых скоростях, то вязкие силы пренебрежимо малы и нет необходимости использовать уравнения Навье Стокса. Однако для определения сопротивления такого крыла при движении его в воздухе вязкие силы оказываются определяющими и необходимо использовать более сложный математический аппарат, связанный с уравнениями Навье Стокса. Интеграл Бернулли. При некоторых предположениях уравнения гидромеханики идеальной жидкости можно один раз проинтегрировать, они имеют решения, одним из которых является интеграл Бернулли для стационарных течений (по имени современника Эйлера математика Бернулли, впервые получившего этот интеграл)![]() ![]()
![]() ![]() где ратм атмосферное давление у выходного отверстия. Отсюда легко получить формулу для скорости истечения V2. В частном случае ратм = р1 получаем так называемую формулу Торичелли для истечения жидкости из широкого сосуда с узким выходным отверстием ![]() Уравнения движения идеальной жидкости имеют еще один интеграл для нестационарных течений, который называется интегралом Коши Лагранжа. Он справедлив для течений, в которых отсутствуют вихри. Его часто, например, используют при рассмотрении волновых движений жидкости или газа. Ударные волны как одно из важных проявлений сжимаемости газа. Математически уравнения идеальной гидроаэромеханики допускают разрывные решения, т.е. решения, которые имеют скачки параметров газа (плотности, давления, скорости и температуры). Одним из таких проявлений в природе является образование ударной волны около летящего со сверхзвуковой скоростью тела в плотных слоях атмосферы Земли. Например, образование ударной волны около летающих сверхзвуковых самолетов или ударных волн около метеоритов, вторгающихся в плотные слои атмосферы Земли с большими сверхзвуковыми скоростями. В условиях космического пространства хорошо известны межпланетные ударные волны, которые чаще всего являются результатом активных процессов на Солнце (например, вспышек).Известно, что около пассажирских самолетов, летающих главным образом с большими дозвуковыми, никакие ударные волны не образуются. Пусть есть сферическое тело радиуса R (рис. 6), которое летит в воздухе со сверхзвуковой скоростью. Тогда впереди такого тела образуется ударная волна В, являющаяся границей между областями 1 и 2, которые отличаются значениями параметров газа. В системе координат, связанной с летящим телом. поток газа набегает на покоящееся тело. Пусть ось Оx направлена вдоль скорости потока, а V1, p1, r1 и T1 скорость, давление, плотность и температура, соответственно, в невозмущенном телом потоке газа (до ударной волны). В область 1 возмущения от тела не попадают, поскольку тело движется со сверхзвуковой скоростью. Так как скорость газа в лобовой точке тела А обращается в нуль, то от точки А до точки С на ударной волне есть область дозвуковой скорости газа, которой достигают возмущения воздуха от летящего тела. Физический смысл образования ударной волны и заключается в разделении невозмущенного и возмущенного потоков газа. Если через V2, p2, r2 и T2 обозначить скорость, давление, плотность и температуру газа соответственно сразу же после ударной волны В, то справедливы неравенстваV2 < V1, p2 > p1 , r2 > r1 , T2 > T1. Это означает, что скорость за ударной волной уменьшается, а давление, плотность и температура возрастают. Сильным возрастанием температуры за ударной волной и объясняется оплавление возвращающихся на Землю космических аппаратов и метеоритов, вторгающихся в атмосферу с большими сверхзвуковыми скоростями. Такие ударные волны называются ударными волнами сжатия (плотность газа возрастает). Интересно, что в природе никогда не наблюдались ударные волны разрежения, в которых плотность падает. Математически образование ударных волн разрежения запрещается известной в гидроаэромеханике теоремой Цемплена Соотношения между параметрами с индексами «1» и «2» можно получить из интегральных законов сохранения массы, импульса и энергии, поскольку они справедливы и для разрывных функций. Такие соотношения называются соотношениями Гюгонио и имеют вид (в системе координат, связанной с ударной волной) r1 Vn1 = r2 Vn2 ; r1 Vn1V1 + p1 n = r2 Vn2V2 + p2 n ; [r1 V12/2 + p1 g/(g 1)]Vn1 = [r2 V22/2 + p2g/(g 1)]Vn2.Вместе с уравнением состояния эти соотношения позволяют определить значения параметров газа за ударной волной (индекс «2») по значениям параметров невозмущенного ударной волной потока газа (индекс «1»). Описанный математический аппарат гидроаэромеханики используется во многих областях естественных наук, при этом для корректности использования этого аппарата требуется только выполнение критерия сплошности среды, т.е. для газов, например, длина свободного пробега частиц должна быть много меньше характерных размеров рассматриваемых объектов обтекания. В частности, в условиях космического пространства часто среда очень разрежена. В таких средах, конечно же, длина свободного пробега частиц очень велика, но размеры самих объектов исследования оказываются во многих случаях существенно больше, т.е. методы гидроаэромеханики применимы и к таким объектам. В биомеханике при помощи методов гидромеханики исследуются интересные особенности течений биологических жидкостей по сосудам, а в гидрогеологии исследуются, например, задачи динамики внутренних слоев Земли. Все это свидетельствует о важности науки, которая называется «гидроаэромеханика». Владимир Баранов Чепмен С. и Каулинг Т. Математическая теория неоднородных газов. ИИЛ, М., 1960 Кочин Н.Е., Кибель и Розе. Теоретическая гидромеханика, т.1. Физматгиз, 1963 Кочин Н.Е., Кибель и Розе, Теоретическая гидромеханика, т.2, Физматгиз, 1963 Седов Л.И. Механика сплошной среды. М., Наука, т. 1, 1973 Седов Л.И. Механика сплошной среды. М., Наука, т. 2, 1973 Баранов В.Б. и Краснобаев К.В., Гидродинамическая теория космической плазмы, М., Изд. «Наука», 1977 Эйлер Л. Общие законы движения жидкостей. Известия РАН, сер. МЖГ, 1999, № 6
|