История гидродинамики
Страницы: 1 2
Непознанная турбулентность
Однако вскоре выяснилось, что в подавляющем большинстве реальных случаев ламинарное течение слабовязкой жидкости становится неустойчивым и рано или поздно переходит в турбулентное [3]:
Турбуле́нтность [...] — явление, заключающееся в том, что при превышении некоторого критического числа Рейнольдса и/или Релея (в частном случае при превышении скорости потока при постоянной плотности и диаметре трубы и/или температуры на внешней границе среды) самопроизвольно образуются многочисленные нелинейные фрактальные волны и обычные, линейные различных размеров, без наличия внешних, случайных, возмущающих среду сил и/или при их присутствии.
Общий критерий возникновения турбулентности установлен Осборном Рейнольдсом (1842–1912) в 1883 году. Он подкрашивал ламинарную струйку тока во входной части стеклянной трубки и следил, когда течение станет турбулентным, фиксируя при этом критическое значение безразмерного определяющего параметра, названного впоследствии в его честь числом Рейнольдса. В Манчестерском университете, где Рейнольдс проводил свои опыты, сохранилась его экспериментальная установка.

Рисунок 1. Порождение турбулентности решеткой. Источник
Представив скорость
в виде суммы средней
и пульсационной
составляющих, а давление в виде
, Рейнольдс получил уравнения для средних величин, носящие его имя:
![]() | (7) |
Здесь подразумевается суммирование по индексу j. По сравнению с уравнением Навье–Стокса (6) это уравнение включает дополнительные напряжения — так называемые напряжения Рейнольдса
. Попытки найти их вид из первых принципов физики оказались безуспешными, поэтому уравнение (7) стало базой для развития эмпирических теорий.
К недоопределенной системе (7) следует добавить уравнения для пульсаций и правила осреднения. Полученную таким образом расширенную систему уравнений следует называть системой Рейнольдса. Решение этой системы кроме среднего и пульсационного слагаемого содержит еще волновой член.
Л.В. Келлер и А.А. Фридман дали аналитическую формулировку проблемы турбулентности, сведя задачу к бесконечной системе уравнений для статистических моментов. Джеффри Тейлор ввёл в рассмотрение корреляционные функции, а также впервые ввёл понятия об однородной и изотропной турбулентности.
Льюис Фри Ричардсон (1881–1953) высказал глубокие соображения о «каскадном процессе» передачи энергии по спектру от крупномасштабных мод к мелкомасштабным. Эта картина развитой турбулентности изображена Ричардсоном в стихотворении, которое входит во многие учебники по турбулентности:
Big whorls have little whorls,
Which feed on their velocity;
Little whorls have smaller whorls,
And so on unto viscosity.
Ричардсон любил поэзию. Свою интуитивную теорию турбулентности он создал, вдохновлённый наблюдениями за эволюцией облаков.
Парадигма турбулентности ознаменовала превращение гидродинамики из чисто теоретической науки, из коллекции оригинальных, но оторванных от жизни решений, в прикладную науку. Разлившаяся река стремилась войти в берега актуальной полезности.
Осборн Рейнольдс (Osborne Reynolds; 1842–1912), английский инженер и физик. Родился в Белфасте в семье священнослужителя. С 18 лет работал в механической мастерской, поступил в Кембриджский университет, где изучал математику и механику. Окончил университет в 1867 г. В 1868–1905 гг. — профессор кафедры строительной механики Манчестерского университета. С 1888 г. возглавлял Витвортовскую инженерную лабораторию. Работы Рейнольдса посвящены механике, гидродинамике, теплоте, электричеству, магнетизму. В 1883 г. Рейнольдс установил, что ламинарное течение переходит в турбулентное, когда введенная им безразмерная величина (число Рейнольдса) превышает критическое значение. Число Рейнольдса широко используется при решении задач гидро- и аэродинамики в случае малых и средних дозвуковых скоростей. Подробнее
Наиболее важным вкладом в теорию турбулентности считаются 4 работы Андрея Николаевича Колмогорова (1903–1987), впервые опубликованные в 1941 году. В совокупности они так и называются — «теория К 41». Колмогоров занимался турбулентностью примерно полгода, до начала войны, а затем, в связи с требованиями военного времени, занялся баллистикой. Ещё 2 работы Колмогорова относятся к 1961–1962 годам и называются «теория К 62».
Были ли предшественники у Рейнольдса? Разумеется, каждый человек видел турбулентность в атмосфере, в ручье, в струе воды. Хаген ещё в 1839 году установил, что характер течения воды в трубе изменяется с увеличением напора. Он считал, что при определённых условиях в потоке возникают внутренние вихри, и это приводит к повышению сопротивления, и, следовательно, к уменьшению расхода воды.
Массивным блоком в основании теории турбулентности лежит теория вероятностей, прекрасные исторические обзоры которой даны Ширяевым (1998) и Гнеденко (2001). В теории турбулентности есть теория устойчивости и перехода ламинарного течения в турбулентное. Рэлею в 1880 году, В. Орру и А. Зоммерфельду в 1907–1908 годах удалось линеаризовать задачу устойчивости в частном случае плоскопараллельного движения. Большое значение имеет теория устойчивости пограничного слоя и ламинарно-турбулентного перехода в пограничном слое (Бэтчелор, 2000). Следует упомянуть далеко продвинутую теорию изотропной турбулентности (Бэтчелор, 1955).
В XIX веке основной гидродинамической моделью являлась модель идеальной жидкости, в XIX–XX веках — модель вязкой жидкости. В XXI веке учёные стали вплотную заниматься построением физических и математических моделей турбулентности. Следует отметить, что парадигма турбулентности в гидродинамике еще не сформулирована, пока видны только её контуры... Построены математические модели случайных процессов в теории броуновского движения и в квантовой механике. Однако построить адекватную модель истории человечества, макроэкономики и турбулентности пока не удалось.
Невозможно провести чёткую границу между гидродинамикой и нелинейной механикой. И там, и там важна нелинейность. Актуально как изучение гидродинамических уравнений (Буссинеска, Кортевега–де Вриза и прочих), так и решение уравнений нелинейной механики (Шрёдингера, Гинзбурга–Ландау, Синус–Гордона и прочих), которые находят важные применения в гидродинамике.
Выдающийся физик Вернер Гейзенберг (1901–1976), серьёзно занимавшийся проблемой турбулентности, признался на смертном одре, что хотел бы задать Господу два вопроса — об основах теории относительности и о причине турбулентности. «Думаю, что Господь ответит мне на первый из них», — галантно заключил он.
Численное моделирование гидродинамических процессов
Усложнение математических моделей привело к тому, что человек перестал справляться с расчетами. Ответом на это явилось создание компьютера. В то время казалось, что все задачи будут немедленно решены, что восторжествовал принцип пифагорейцев: «Всё есть число». Однако конец эмпирической эпохи в науке и технике не наступил, отодвинувшись на следующее столетие. Современный компьютер позволяет решать сложные задачи гидродинамики, физики плазмы, анализа загрязнения воздуха и грунтовых вод, создания новых лекарств, составления карт озонового слоя, сейсмического анализа. Численные расчёты востребованы как в фундаментальных науках, так и в прикладных.
Для учёных и инженеров численный расчёт стал обыденным и мощным оружием наступления на непознанное. Численные расчёты практически вытеснили из употребления такие эмпирические методы, как приближения Бубнова–Галёркина в теории пограничного слоя, методы типа Бетца–Мультхопа в теории крыла большого удлинения, полуэмпирические методы в теории турбулентности и так далее. В настоящее время разработаны надёжные методы расчёта ламинарных течений, однако безэмпирический расчёт турбулентных течений пока ещё невозможен. Вычислительная гидродинамика — основная часть вычислительной математики, ибо в ней, как ни в одной области физики, задействована нелинейность и, как следствие, турбулентность.
Компьютерная модель жидкости представляет собой приближение математической модели. Она представляет собой в некоторой степени «модель модели».
Расскажу вкратце об одном подходе к численному моделированию жидкости: методе конечных элементов. Для того чтобы хранить состояние жидкости в памяти компьютера, имеющей конечный объём, в жидкости выделяют конечное число точек, параметры в которых вычисляются. Считается, что каждая точка непосредственно взаимодействует лишь с несколькими соседними точками. Если соединить взаимодействующие точки отрезками, то вычисляемая область покроется сеткой (см. рисунок 2). Такая сетка называется расчётной сеткой, а вычисляемые точки — узлами этой сетки.
Обычно за основу берутся уравнения Навье–Стокса. Они затем модифицируются путём введения в рассмотрение важных факторов: глубины, гравитации, осадков, силы Кориолисса, ветра, вязкости, турбулентной диффузии, стоков рек,…
Параметры жидкости (давление, компоненты скорости, солёность, температура и другие) в узлах сетки связываются уравнениями, которые получаются интегрированием исходных дифференциальных уравнений по ячейкам сетки. Образуется сеточная система уравнений.
Время тоже дробится на отрезки. Выводится зависимость между значениями параметров жидкости в сетке в текущий момент времени и, например, значениями параметров секунду назад (это называется временным шагом). Таким образом, секунда за секундой, мы можем рассчитать год существования водоёма.
Получается, что параметры каждого узла сетки находятся в зависимости от окружающих его узлов (с которыми он связан рёбрами). Узлы сетки как бы являются бочками с водой, а рёбра — трубами, соединяющими эти бочки. Остаётся только надеяться, что подобная «сеть труб и бочек», будучи достаточно подробной, даст приближение к течению жидкости в реальном водоёме.
После написания компьютерной программы, она запускается на счёт для различных значений параметров. Каждое вычисление карты течений с использованием персонального компьютера может занимать несколько часов и даже дней. Чтобы уменьшить время вычислений создаются параллельные программы, способные работать одновременно на десятках компьютеров, объединяя их вычислительные ресурсы.

Рисунок 3. Результат моделирования. Северный ветер, 5 м/с.
Цветом показано возвышение уровня воды из-за ветра
И вот, наконец, после многократного исправления программы и доработки модели мы получаем результаты, которые считаем итоговыми. Но как узнать, являются ли эти результаты верными? Соответствуют ли они действительности? Для этого используется так называемая верификация модели: модель проверяется на исходных данных, результат для которых известен (измерен или точно вычислен). Это может быть, например, другая, простая форма водоёма с заданными погодными условиями, для которой известна картина течений. И, если модель (программа) дала точный результат на известных данных, считается, что она будет давать приемлемые результаты и на других, практически полезных данных, результат для которых неизвестен. Это похоже на сдачу экзаменов по вождению автомобиля: если ученик знает правила (программа написана без ошибок) и выполнил все упражнения (программа даёт хорошие результаты на тестовых задачах), то считается, что в дальнейшем он будет хорошо ездить на автомобиле по городу.
Однако численное моделирование течений жидкости не так просто и безоблачно, как может показаться. Вот, на мой взгляд, некоторые проблемы вычислительной гидродинамики, над решением которых сейчас упорно работают исследователи:
- Прежде всего, уравнения Навье–Стокса — нелинейные. Возможно существование нескольких решений, численные методы достаточно сложны, трудно добиться гарантированной сходимости.
- Жидкость — сложная субстанция, обладающая различными меняющимися во времени и пространстве характеристиками. Зачастую получается так, что при вычислении одной характеристики мы нарушаем другую. Требуются всякие изощрённые методы, чтобы этого не случалось (например, поправка к давлению).
- Нужно добиться выполнения законов сохранения энергии и импульса. Дифференциальные уравнения гидродинамики содержат эти законы лишь косвенно: если их проинтегрировать по всей области, то получится, что суммарные энергия и импульс в ней сохраняются (не меняются со временем). Численные схемы, даже если они неплохо приближают уравнения (и течения) локально, могут давать существенные расхождения в таких глобальных параметрах, как энергия и импульс.
- Схемная вязкость: каждая ячейка/точка/элемент сетки обладают некоторым набором параметров. Но этот набор — конечный. Например, мы можем хранить средние давление и скорость жидкости в ячейке. Это означает, что на каждом шаге по времени мы будем упорно терять все те особенности течения, которые не могут быть описаны имеющимся в ячейке (и в сетке в целом) набором параметров. В случае хранения средних значений, эти особенности как бы мгновенно «размазываются» по всей ячейке. Это напоминает дополнительную вязкость, которая будет зависеть от размеров и формы ячеек сетки. Моделирование течения невязкой жидкости становится крайне сложным.
- Плохая устойчивость: зачастую приходится решать задачу гидродинамики при очень малых значениях временного шага, так как при больших шагах схема попросту «разваливается» (возникают нарастающие колебания параметров). Большое число временных шагов приводит к медленным вычислениям, накоплению погрешностей, и усиливает схемную вязкость.
- Турбулентность: до сих пор не до конца исследованное явление. Если с расчетом ламинарных потоков учёные научились справляться, то для решения задач с турбулентностью каждый раз приходится думать. Турбулентность возникает резко и на всех масштабах расстояния. Возникают проблемы учёта масштабов турбулентности, сравнимых с размером ячеек и меньших их. Для этого существуют различные модели, так или иначе приближающие турбулентность и её средние характеристики.
- В прикладных задачах зачастую требуется учесть дополнительные факторы, что значительно усложняет моделирование: деформация стенок трубопровода, температура, солёность, рельеф дна, волнение на поверхности (подвижная верхняя граница), ветер, сила Кориоллиса, испарение и осадки, приливы и отливы, подземные источники и прочее, и прочее...
- А ещё бывают сверхзвуковые течения, образование пузырьков (разрыв потока), электромагнитные явления в жидкости при её течении в магнитном поле...
В заключение этого пункта приведу ещё несколько компьютерных иллюстраций численного моделирования течения жидкости:

Рисунок 4. Образование вихрей в вязкой жидкости возле границ из одного начального вихря, смещённого относительно центра области. Показаны изолинии давления: давление внутри вихря всегда ниже, чем давление снаружи
Заключение
В данной работе не затронуты поверхностное натяжение жидкости, фазовые переходы и прочие свойства жидкости, как химической и физической субстанции. Рассмотрено лишь течение (движение) жидкости.
Кроме того, у самой знакомой нам жидкости — воды — есть самые удивительные свойства из всех жидкостей на Земле. Очень рекомендую ознакомиться с книгой [4].
Надеюсь, что я смог заинтересовать читателя удивительной наукой гидродинамикой. История этой науки, очевидно, не закончена. Ещё много столетий она будет занимать умы учёных. Кто знает, возможно, в будущем появятся простые уравнения, описывающие всё многообразие феноменов, наблюдаемых при движении воды. Но может случиться и так, что уравнения (и представления о жидкости) будут усложняться, включая в себя всё новые и новые явления.
Список использованных материалов
После номера стоит процент, показывающий, какая часть текста взята из данного источника.
- 30%. 11th Edition of the Enciclopaedia Britannica. 1911 г. На английском языке. Особенно интересна статья «Hydromechanics».
- 30%. Бетяев С. К. «Пролегомены к метагидродинамике». РХД, 2006 г. 304 стр.
- 20%. Wikipedia, the Free Encyclopedia. На английском и русском языках.
- 0%. Леонид Адольфович Кульский, Воля Васильевна Даль, Людмила Григорьевна Ленчина. «Вода знакомая и загадочная». Издательство «Радянська школа», 1982 г. 53 стр. Ссылка
Остальные 20% — собственный текст.
Страницы: 1 2.
![$$\rho\frac{\partial\bar{u}_j\bar{u}_i}{\partial x_j}=\frac{\partial}{\partial x_j}\left[-\bar{p}\delta_{ij}+\mu\left(\frac{\partial\bar{u}_i}{\partial x_j}+\frac{\partial\bar{u}_j}{\partial x_i}\right)-\rho\overline{u_i^\prime u_j^\prime}\right].$$ \rho\frac{\partial\bar{u}_j\bar{u}_i}{\partial x_j} = \frac{\partial}{\partial x_j}
\left[-\bar{p}\delta_{ij}+\mu\left(\frac{\partial\bar{u}_i}{\partial x_j}+\frac{\partial\bar{u}_j}{\partial x_i}\right)-
\rho\overline{u_i^\prime u_j^\prime}\right].](http://iproc.ru/wp-content/cache_tex/tex_dbbb138b3632b3be4634b3fbc606bc91.png)




Девять отзывов на «историю гидродинамики»
Автор: Ольга. Дата: 5-го февраля 2009 г. Время: 08:39.
Огромное спасибо за «Историю гидродинамики». Моя дочь Вашего возраста в обозримом будущем защищает кандидатскую по одной из тем турбулентности, а моя инженерная специальность от этого далека. Теперь мне будет легче найти общий язык с дочерью.
Блестящая статья, гидродинамикой заинтересовали!
Автор: Борис. Дата: 20-го марта 2010 г. Время: 23:33.
Антон, великолепный сайт. Недавно начал тесно заниматься с CFD-комплексами и математическим и инженерным анализом сложных расчетов — на ваш сайт наткнулся в поиске информации для реферата по истории науки для аспирантуры — очень понравилась статья ваша по гидродинамике, хочу использовать ее в своей работе — в списке источников пренепременно укажу ваш сайт
Автор: Антон. Дата: 20-го марта 2010 г. Время: 23:42.
Статья по истории гидродинамики — это и есть мой реферат по истории науки для аспирантуры. Вот такие совпадения.
Автор: Дима. Дата: 29-го марта 2010 г. Время: 13:17.
Реферат по истории науки для аспирантуры — актуальная тема! Антон, если Вы не возражаете, я тоже воспользуюсь Вашим грамотно обобщенным материалом и тоже пренепременно укажу сайт в списке источников. Рад видеть людей с близкими интересами!
Автор: Евгений. Ростов-на-Дону. Дата: 29-го марта 2011 г. Время: 16:26.
Мне тоже для реферата для аспирантуры на кандидатский нужна история гидродинамики! Классно, спасибо за информацию… Нас потом всех не попрут из аспирантуры за одинаковые рефераты?
Автор: Евгений. Дата: 19-го мая 2011 г. Время: 08:13.
Теперь нас уже трое — тоже реферат, тоже кандидатский!
Автор: Юрий. Дата: 2-го июля 2010 г. Время: 15:42.
Ограничение, что жидкость несжимаема , отбрасывает наличие в жидкости энергии сжатия (коэффициент сжимаемости = 0) .
Но жидкость продолжает иметь текучесть (изменяет форму не изменяя объем) и продолжает течь по трубам .
При истечении жидкости через отверстие из сосуда Торричелли физики могут наблюдать непрерывное преобразование потенциальной энергии накопленной жидкостью в сосуде в поле силы тяготения Земли в кинетическую ( то есть энергию движения) .
Поскольку жидкость несжимаема , то приобретаемая жидкостью кинетическая энергия может быть объяснена только на основе динамики твердого тела .
А именно всем известного второго закона Ньютона.
Его фундамент – изменение скорости тела массой М под действием силы
(то есть влияния на него другого тела).
У нас это влияние выражено в виде воздействия одной части жидкости (одного тела) на другую ее часть (другое тело) .
При выталкивании части жидкости из сосуда (представляем себе твердое тело) , она за конечное время приобретает конечную скорость движения .
А значит и конечное воздействие со стороны жидкости находящейся в сосуде (иначе бы скорость постоянно увеличивалась, что противоречит теореме Торричелли ).
Таким образом происходит не только изменение скорости под действием силы , но что не менее важно — изменение силы воздействия (от максимального значения до нулевого ) при изменении скорости движения вытекаемой из отверстия жидкости.
А уравнение Бернулли – следствие теоремы Торричелли .
Нет сосуда Торричелли – неприменимо уравнение Бернулли .
Движение сжатой жидкости в замкнутой трубе различных сечений не содержит сосуда Торричелли ( жидкость имеет кинетическую энергию и энергию
сжатия ). Как применить уравнение Бернулли в этом случае?
Статья Антона познавательна. Спасибо.
Автор: Виктор. Дата: 15-го мая 2011 г. Время: 12:27.
«История газовой механики изложена Я.Б. Зельдовичем. Ударную волну открыл «на кончике пера» Б. Риман в 1876 году. Однако гипотезу о существовании ударной волны задолго до него, в 1848 году, высказал Стокс, но отказался от нее под влиянием критики В. Томсона и своего ученика лорда Рэлея.»
Риман никак не мог в 1876 году открыть ударную волну (он в 1866 году умер). Это было 1860-м. Исправьте, пожалуйста.
Автор: Виктор. Дата: 15-го мая 2011 г. Время: 12:52.
«Однако гипотезу о существовании ударной волны задолго до него, в 1848 году, высказал Стокс, но отказался от нее под влиянием критики В. Томсона и своего ученика лорда Рэлея.»
Здесь ошибка перевода. В 1848 году Стокс высказал гипотезу о существовании ударной волны, но отказался от своей работы. Однако Лорд Кельвин и Лорд Рэлей позже убедили его что он был всё же прав. Как то так. Источник http://www.potto.org/gasDynamics/node44.html
Ну никак в 1848 году Рэлей никак не мог убедить Стокса (Рэлею тогда было 8 лет)…