Уравнения мелкой воды

Выходные данные модели воды в ванне с помощью уравнения мелкой воды. Вода испытывает пять брызг, которые создают поверхностные гравитационные волны, которые распространяются от мест брызг и отражаются от стенок ванны.

Уравнения мелкой воды ( SWE ) представляют собой набор гиперболических уравнений в частных производных (или параболических, если рассматривается вязкий сдвиг), которые описывают поток под поверхностью давления в жидкости (иногда, но не обязательно, свободной поверхностью ). Уравнения мелкой воды в однонаправленной форме также называются (де) уравнениями Сен-Венана в честь Адемара Жан-Клода Барре де Сен-Венана (см. соответствующий раздел ниже).

Уравнения получены путем интегрирования по глубине уравнений Навье – Стокса в случае, когда горизонтальный масштаб длины намного больше вертикального масштаба длины. При этом условии сохранение массы означает, что масштаб вертикальной скорости жидкости мал по сравнению с масштабом горизонтальной скорости. Из уравнения количества движения можно показать, что вертикальные градиенты давления почти гидростатические и что горизонтальные градиенты давления возникают из-за смещения поверхности давления, а это означает, что поле горизонтальной скорости постоянно по всей глубине жидкости. Вертикальное интегрирование позволяет исключить вертикальную скорость из уравнений. Таким образом выводятся уравнения мелкой воды.

Хотя член вертикальной скорости не присутствует в уравнениях мелкой воды, обратите внимание, что эта скорость не обязательно равна нулю. Это важное различие, потому что, например, вертикальная скорость не может быть равна нулю, когда дно меняет глубину, и, следовательно, если бы она была равна нулю, в уравнениях мелкой воды можно было бы использовать только плоские днища. Как только решение (т.е. горизонтальные скорости и смещение свободной поверхности) найдено, вертикальная скорость может быть восстановлена ​​с помощью уравнения неразрывности.

В гидродинамике часто встречаются ситуации, когда горизонтальный масштаб длины намного больше вертикального, поэтому уравнения мелкой воды широко применимы. Они используются вместе с силами Кориолиса в моделировании атмосферы и океана как упрощение примитивных уравнений атмосферного потока.

Модели уравнений мелководья имеют только один вертикальный уровень, поэтому они не могут напрямую учитывать какой-либо фактор, изменяющийся с высотой. Однако в тех случаях, когда среднее состояние достаточно простое, вертикальные изменения можно отделить от горизонтальных, и это состояние можно описать несколькими наборами уравнений мелкой воды.

Уравнения

Одномерная диаграмма, представляющая модель мелкой воды.

Консервативная форма

Уравнения мелкой воды выводятся из уравнений сохранения массы и сохранения линейного импульса ( уравнения Навье-Стокса ), которые справедливы даже тогда, когда предположения о мелкой воде не работают, например, при гидравлическом прыжке . В случае горизонтального дна с пренебрежимо малыми силами Кориолиса , трения и силами вязкости уравнения мелкой воды имеют вид: $${\begin{aligned}{\frac {\partial (\rho \eta )}{\partial t}}&+{\frac {\partial (\rho \eta u)}{\partial x}}+{\frac {\partial (\rho \eta v)}{\partial y}}=0,\\[3pt]{\frac {\partial (\rho \eta u)}{\partial t}}&+{\frac {\partial }{\partial x}}\left(\rho \eta u^{2}+{\frac {1}{2}}\rho g\eta ^{2}\right)+{\frac {\partial (\rho \eta uv)}{\partial y}}=0,\\[3pt]{\frac {\partial (\rho \eta v)}{\partial t}}&+{\frac {\partial }{\partial y}}\left(\rho \eta v^{2}+{\frac {1}{2}}\rho g\eta ^{2}\right)+{\frac {\partial (\rho \eta uv)}{\partial x}}=0.\end{aligned}}$$

Здесь η — общая высота столба жидкости (мгновенная глубина жидкости как функция x , y и t ), а двумерный вектор ( u , v горизонтального потока жидкости ) — скорость , усредненная по вертикальному столбу. Далее g — ускорение свободного падения, а ρ — плотность жидкости . Первое уравнение получено из закона сохранения массы, вторые два — из закона сохранения импульса.

Неконсервативная форма

Разлагая производные, приведенные выше, с помощью правила произведения , получаем неконсервативную форму уравнений мелкой воды. Поскольку скорости не подчиняются фундаментальному уравнению сохранения, неконсервативные формы не сохраняются при ударе или гидравлическом прыжке . Также включены соответствующие условия для Кориолиса, сил трения и вязкости, чтобы получить (для постоянной плотности жидкости): $${\begin{aligned}{\frac {\partial h}{\partial t}}&+{\frac {\partial }{\partial x}}{\Bigl (}(H+h)u{\Bigr )}+{\frac {\partial }{\partial y}}{\Bigl (}(H+h)v{\Bigr )}=0,\\[3pt]{\frac {\partial u}{\partial t}}&+u{\frac {\partial u}{\partial x}}+v{\frac {\partial u}{\partial y}}-fv=-g{\frac {\partial h}{\partial x}}-ku+\nu \left({\frac {\partial ^{2}u}{\partial x^{2}}}+{\frac {\partial ^{2}u}{\partial y^{2}}}\right),\\[3pt]{\frac {\partial v}{\partial t}}&+u{\frac {\partial v}{\partial x}}+v{\frac {\partial v}{\partial y}}+fu=-g{\frac {\partial h}{\partial y}}-kv+\nu \left({\frac {\partial ^{2}v}{\partial x^{2}}}+{\frac {\partial ^{2}v}{\partial y^{2}}}\right),\end{aligned}}$$

где

в	это скорость в направлении x или зональная скорость
v	- скорость в направлении y или меридиональная скорость
ЧАС	средняя высота горизонтальной поверхности давления
час	- отклонение высоты горизонтальной поверхности давления от ее средней высоты, где h : η ( x , y , t ) = H ( x , y ) + h ( x , y , t )
б	— топографическая высота от эталона D, где b : H ( x , y ) = D + b ( x , y )
г	это ускорение свободного падения
ж	– коэффициент Кориолиса, связанный с силой Кориолиса . На Земле f равна 2 Ом sin( φ ), где Ω — угловая скорость вращения Земли (π/12 радиан /час), а φ — широта.
к	коэффициент вязкого сопротивления
н	кинематическая вязкость

Анимация линеаризованных уравнений мелкой воды для прямоугольного бассейна без учета трения и силы Кориолиса. Вода испытывает всплеск, который генерирует поверхностные гравитационные волны, которые распространяются от места всплеска и отражаются от стенок бассейна. Анимация создана с использованием точного решения Кэрриера и Йе (2005) для осесимметричных волн.

Часто бывает, что члены, квадратичные по u и v , которые отражают эффект объемной адвекции , малы по сравнению с другими членами. Это называется геострофическим балансом и эквивалентно утверждению, что число Россби мало. Предполагая также, что высота волны очень мала по сравнению со средней высотой ( h ≪ H ), имеем (без учета боковых вязких сил): $${\begin{aligned}{\frac {\partial h}{\partial t}}&+H\left({\frac {\partial u}{\partial x}}+{\frac {\partial v}{\partial y}}\right)=0,\\[3pt]{\frac {\partial u}{\partial t}}&-fv=-g{\frac {\partial h}{\partial x}}-ku,\\[3pt]{\frac {\partial v}{\partial t}}&+fu=-g{\frac {\partial h}{\partial y}}-kv.\end{aligned}}$$

Одномерные уравнения Сен-Венана

Одномерные (1-D) уравнения Сен-Венана были выведены Адемаром Жан-Клодом Барре де Сен-Венаном и обычно используются для моделирования переходного течения в открытом канале и поверхностного стока . Их можно рассматривать как сокращение двумерных (2-D) уравнений мелкой воды, которые также известны как двумерные уравнения Сен-Венана. Одномерные уравнения Сен-Венана в определенной степени содержат основные характеристики формы поперечного сечения канала .

Одномерные уравнения широко используются в компьютерных моделях, таких как TUFLOW , Mascaret (EDF), SIC (Irstea) , HEC-RAS , SWMM5, ИГИЛ, ИнфоВоркс, Разработчик моделей наводнений, SOBEK 1DFlow, MIKE 11 , и МАЙК ШЕ, потому что их значительно легче решить, чем полные уравнения мелкой воды. Общие применения одномерных уравнений Сен-Венана включают маршрут паводков вдоль рек (включая оценку мер по снижению риска наводнений), анализ прорыва плотин, штормовые импульсы в открытом русле, а также ливневые стоки в наземном потоке.

Уравнения

Поперечное сечение открытого канала.

Система уравнений в частных производных , которые описывают одномерное течение несжимаемой жидкости в открытом канале произвольного поперечного сечения , полученная и сформулированная Сен-Венаном в его статье 1871 года (уравнения 19 и 20), — это:

$${\frac {\partial A}{\partial t}}+{\frac {\partial \left(Au\right)}{\partial x}}=0$$ ( 1 )

и

$${\frac {\partial u}{\partial t}}+u\,{\frac {\partial u}{\partial x}}+g\,{\frac {\partial \zeta }{\partial x}}=-{\frac {P}{A}}\,{\frac {\tau }{\rho }},$$ ( 2 )

где x — пространственная координата вдоль оси канала, t — время, A ( x , t поперечного сечения ) — площадь потока в точке x , u ( x , t ) — скорость потока , ζ ( x , t) ) — высота свободной поверхности , а τ( x , t стенки ) — напряжение сдвига вдоль смоченного периметра P ( x , t ) поперечного сечения в точке x . Далее ρ — (постоянная) плотность жидкости , а g — ускорение свободного падения .

Замыкание гиперболической системы уравнений ( 1 )–( 2 ) получается из геометрии поперечных сечений – путем обеспечения функциональной связи между площадью поперечного сечения A и высотой поверхности ζ в каждом положении x . Например, для прямоугольного поперечного сечения с постоянной шириной канала B и высотой русла z _b площадь поперечного сечения равна: A = B (ζ - z _b ) = B h . Мгновенная глубина воды равна h ( x , t ) = ζ( x , t ) − z _b ( x ) , где z _b ( x ) уровень дна (т. е. высота самой низкой точки дна над нулевой отметкой , см. крестик). - рисунок сечения ). Для неподвижных стенок канала площадь поперечного сечения A в уравнении ( 1 ) можно записать как: $$A(x,t)=\int _{0}^{h(x,t)}b(x,h')\,dh',$$где b ( x , h ) эффективная ширина поперечного сечения канала в месте x , когда глубина жидкости равна h – поэтому b ( x , h ) = B ( x ) для прямоугольных каналов.

Напряжение сдвига стенки τ зависит от скорости потока u , их можно связать, используя, например, уравнение Дарси-Вейсбаха , формулу Мэннинга или формулу Шези .

Далее, уравнение ( 1 ) является уравнением неразрывности , выражающим сохранение объема воды для этой несжимаемой однородной жидкости. Уравнение ( 2 ) представляет собой уравнение количества движения , дающее баланс между силами и скоростью изменения количества движения.

Уклон пласта S ( x ), уклон трения S _f ( x , t ) и гидравлический радиус R ( x , t ) определяются как: $$S=-{\frac {\mathrm {d} z_{\mathrm {b} }}{\mathrm {d} x}},$$$$S_{\mathrm {f} }={\frac {\tau }{\rho gR}}$$и $$R={\frac {A}{P}}.$$

Следовательно, уравнение количества движения ( 2 ) можно записать как:

$${\frac {\partial u}{\partial t}}+u\,{\frac {\partial u}{\partial x}}+g\,{\frac {\partial h}{\partial x}}+g\,\left(S_{\mathrm {f} }-S\right)=0.$$ ( 3 )

Сохранение импульса

Уравнению количества движения ( 3 ) также можно придать так называемую форму сохранения посредством некоторых алгебраических манипуляций с уравнениями Сен-Венана ( 1 ) и ( 3 ). По разряду Q = Au :

$${\frac {\partial Q}{\partial t}}+{\frac {\partial }{\partial x}}\left({\frac {Q^{2}}{A}}+g\,I_{1}\right)+g\,A\,\left(S_{f}-S\right)-g\,I_{2}=0,$$ ( 4 )

где A , I ₁ и I ₂ являются функциями геометрии канала, описываемыми через ширину канала B (σ, x ). Здесь σ — высота над самой нижней точкой поперечного сечения в месте x , см. рисунок поперечного сечения . Итак, σ — высота над уровнем дна z _b ( x ) (самой низкой точки поперечного сечения): $${\begin{aligned}A(\sigma ,x)&=\int _{0}^{\sigma }B(\sigma ',x)\;\mathrm {d} \sigma ',\\I_{1}(\sigma ,x)&=\int _{0}^{\sigma }(\sigma -\sigma ')\,B(\sigma ^{\prime },x)\;\mathrm {d} \sigma '\qquad {\text{and}}\\I_{2}(\sigma ,x)&=\int _{0}^{\sigma }(\sigma -\sigma ')\,{\frac {\partial B(\sigma ',x)}{\partial x}}\;\mathrm {d} \sigma '.\end{aligned}}$$

Выше – в уравнении импульса ( 4 ) в форме сохранения – A , I ₁ и I ₂ оцениваются при σ = h ( x , t ) . Член g I ₁ описывает гидростатическую силу в определенном поперечном сечении. А для непризматического канала g I ₂ дает эффект изменения геометрии вдоль оси x канала .

В приложениях, в зависимости от решаемой задачи, часто отдается предпочтение использованию либо уравнения импульса в форме несохранения ( 2 ) или ( 3 ), либо формы сохранения ( 4 ). Например, в случае описания гидравлических прыжков форма сохранения предпочтительна, поскольку поток импульса непрерывен по всему прыжку.

Характеристики

Характеристики, область зависимости и область влияния, связанные с местоположением P = ( x _P , t _P ) в пространстве x и времени t .

Уравнения Сен-Венана ( 1 )–( 2 ) можно проанализировать методом характеристик . Две скорости d x /d t на характеристических кривых: $${\frac {\mathrm {d} x}{\mathrm {d} t}}=u\pm c,$$ с $$c={\sqrt {\frac {gA}{B}}}.$$

Число Фруда Fr = | ты | / c определяет, является ли поток докритическим ( Fr < 1 ) или сверхкритическим ( Fr > 1 ).

Для прямоугольного и призматического канала постоянной ширины B , т.е. с A = B h и c = √ gh , инвариантами Римана являются: $$r_{+}=u+2{\sqrt {gh}}$$ и $$r_{-}=u-2{\sqrt {gh}},$$поэтому уравнения в характеристической форме имеют вид: $${\begin{aligned}&{\frac {\mathrm {d} }{\mathrm {d} t}}\left(u+2{\sqrt {gh}}\right)=g\left(S-S_{f}\right)&&{\text{along}}\quad {\frac {\mathrm {d} x}{\mathrm {d} t}}=u+{\sqrt {gh}}\quad {\text{and}}\\&{\frac {\mathrm {d} }{\mathrm {d} t}}\left(u-2{\sqrt {gh}}\right)=g\left(S-S_{f}\right)&&{\text{along}}\quad {\frac {\mathrm {d} x}{\mathrm {d} t}}=u-{\sqrt {gh}}.\end{aligned}}$$

Инварианты Римана и метод характеристик призматического канала произвольного сечения описаны Диденкуловой и Пелиновским (2011).

Характеристики и инварианты Римана дают важную информацию о поведении течения, а также о том, что их можно использовать в процессе получения (аналитических или численных) решений.

Гамильтонова структура для течения без трения

В случае отсутствия трения и канала прямоугольного призматического сечения уравнения Сен-Венана имеют гамильтонову структуру. Гамильтониан H равен энергии течения на свободной поверхности: $$H=\rho \int \left({\frac {1}{2}}Au^{2}+{\frac {1}{2}}gB\zeta ^{2}\right)\mathrm {d} x,$$с постоянной B - ширина канала и ρ - постоянная плотность жидкости . Тогда уравнения Гамильтона таковы: $${\begin{aligned}&\rho B{\frac {\partial \zeta }{\partial t}}+{\frac {\partial }{\partial x}}\left({\frac {\partial H}{\partial u}}\right)=\rho \left(B{\frac {\partial \zeta }{\partial t}}+{\frac {\partial (Au)}{\partial x}}\right)=\rho \left({\frac {\partial A}{\partial t}}+{\frac {\partial (Au)}{\partial x}}\right)=0,\\&\rho B{\frac {\partial u}{\partial t}}+{\frac {\partial }{\partial x}}\left({\frac {\partial H}{\partial \zeta }}\right)=\rho B\left({\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}+g{\frac {\partial \zeta }{\partial x}}\right)=0,\end{aligned}}$$поскольку ∂ А /∂ ζ знак равно B ) .

Производное моделирование

Динамическая волна

Динамическая волна представляет собой полное одномерное уравнение Сен-Венана. Это сложно решить численно, но оно справедливо для всех сценариев течения в русле. Динамическая волна используется для моделирования переходных штормов в программах моделирования, включая Mascaret (EDF), SIC (Irstea) , HEC-RAS , InfoWorks_ICM. Архивировано 25 октября 2016 г. в Wayback Machine . МАЙК 11 , Мойка 123d и SWMM5 .

В порядке возрастания упрощений, удаляя некоторые члены полных одномерных уравнений Сен-Венана (также известных как динамическое волновое уравнение), мы получаем также классическое диффузионное волновое уравнение и кинематическое волновое уравнение.

Диффузионная волна

Для диффузионной волны предполагается, что инерционные члены меньше, чем члены гравитации, трения и давления. Таким образом, диффузионную волну можно более точно описать как неинерционную волну, и она записывается как: $$g{\frac {\partial h}{\partial x}}+g(S_{f}-S)=0.$$

Диффузная волна действительна, когда инерционное ускорение намного меньше, чем все другие формы ускорения, или, другими словами, когда существует преимущественно докритический поток с низкими значениями Фруда. Модели, использующие предположение о диффузных волнах, включают MIKE SHE. и ЛИСФЛУД-ФП. В программном обеспечении SIC (Irstea) эта опция также доступна, поскольку два члена инерции (или любой из них) могут быть удалены в качестве опции из интерфейса.

Кинематическая волна

Для кинематической волны предполагается, что течение однородно, а наклон трения примерно равен наклону канала. Это упрощает полное уравнение Сен-Венана до кинематической волны: $$S_{f}-S=0.$$

Кинематическая волна действительна, когда изменение высоты волны с расстоянием и скорости с расстоянием и временем незначительно по отношению к уклону дна, например, для мелких потоков на крутых склонах. Кинематическая волна используется в HEC-HMS .

Вывод из уравнений Навье – Стокса.

Одномерное уравнение импульса Сен-Венана может быть получено из уравнений Навье-Стокса, описывающих движение жидкости . Компонент x уравнений Навье – Стокса, выраженный в декартовых координатах в направлении x , можно записать как: $${\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}+v{\frac {\partial u}{\partial y}}+w{\frac {\partial u}{\partial z}}=-{\frac {\partial p}{\partial x}}{\frac {1}{\rho }}+\nu \left({\frac {\partial ^{2}u}{\partial x^{2}}}+{\frac {\partial ^{2}u}{\partial y^{2}}}+{\frac {\partial ^{2}u}{\partial z^{2}}}\right)+f_{x},$$

где u — скорость в направлении x , v — скорость в направлении y , w — скорость в направлении z , t — время, p — давление, ρ — плотность воды, ν — кинематическая вязкость, а f _x — массовая сила в направлении x .

1. Если предположить, что трение учитывается как объемная сила, то $\nu$ можно принять равным нулю, поэтому: $$\nu \left({\frac {\partial ^{2}u}{\partial x^{2}}}+{\frac {\partial ^{2}u}{\partial y^{2}}}+{\frac {\partial ^{2}u}{\partial z^{2}}}\right)=0.$$
2. Предполагая одномерный поток в направлении x , отсюда следует, что: $$v{\frac {\partial u}{\partial y}}+w{\frac {\partial u}{\partial z}}=0$$
3. Предполагая также, что распределение давления приблизительно гидростатическое, отсюда следует, что: $$p=\rho gh$$ или в дифференциальной форме: $$\partial p=\rho g(\partial h).$$ И когда эти предположения применяются к x -компоненте уравнений Навье – Стокса: $$-{\frac {\partial p}{\partial x}}{\frac {1}{\rho }}=-{\frac {1}{\rho }}{\frac {\rho g\left(\partial h\right)}{\partial x}}=-g{\frac {\partial h}{\partial x}}.$$
4. На жидкость в канале действуют две объемные силы: сила тяжести и трение: $$f_{x}=f_{x,g}+f_{x,f}$$ где f _x,g — массовая сила, вызванная гравитацией, а f _x,f — массовая сила, вызванная трением.
5. f _{x , g} можно рассчитать, используя основы физики и тригонометрии: $$F_{g}=\sin(\theta )gM$$где F _g — сила тяжести в направлении x , θ — угол, а M — масса.

   

   Рисунок 1: Схема движения блока по наклонной плоскости.

   Выражение для sin θ можно упростить с помощью тригонометрии: $$\sin \theta ={\frac {\text{opp}}{\text{hyp}}}.$$ Для малых θ (разумных практически для всех потоков) можно предположить, что: $$\sin \theta =\tan \theta ={\frac {\text{opp}}{\text{adj}}}=S$$ и учитывая, что f _x представляет собой силу на единицу массы, выражение принимает вид: $$f_{x,g}=gS.$$
6. Если предположить, что линия уровня энергии не совпадает с уклоном канала, и при постоянном уклоне наблюдаются постоянные потери на трение, отсюда следует, что: $$f_{x,f}=S_{f}g.$$
7. В совокупности все эти предположения приводят к одномерному уравнению Сен-Венана в направлении x : $${\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}+g{\frac {\partial h}{\partial x}}+g(S_{f}-S)=0,$$ $$(a)\quad \ \ (b)\quad \ \ \ (c)\qquad \ \ \ (d)\quad (e)\$$ где (a) — член местного ускорения, (b) — член конвективного ускорения, (c) — член градиента давления, (d) — член трения, и (e) — член гравитации.

Условия

Локальное ускорение (а) также можно рассматривать как «нестационарный член», поскольку оно описывает некоторое изменение скорости с течением времени. Конвективное ускорение (b) — это ускорение, вызванное некоторым изменением скорости в зависимости от положения, например, ускорением или замедлением потока жидкости, попадающей в сужение или отверстие соответственно. Оба эти члена составляют члены инерции одномерного уравнения Сен-Венана.

Термин градиента давления (c) описывает, как давление меняется в зависимости от положения, и, поскольку давление предполагается гидростатическим, это изменение напора относительно положения. Член трения (d) учитывает потери энергии из-за трения, а член силы тяжести (e) представляет собой ускорение из-за наклона пласта.

Моделирование волн с помощью уравнений мелкой воды

Уравнения мелкой воды можно использовать для моделирования волн Россби и Кельвина в атмосфере, реках, озерах и океанах, а также гравитационных волн в меньшей области (например, поверхностных волн в ванне). Чтобы уравнения мелкой воды были действительными, длина волны явления, которое они должны моделировать, должна быть намного больше глубины бассейна, в котором это явление имеет место. С несколько меньшими длинами волн можно справиться, расширив уравнения мелкой воды с помощью приближения Буссинеска, включив в него эффекты дисперсии . Уравнения мелкой воды особенно подходят для моделирования приливов, которые имеют очень большие масштабы длины (более сотен километров). Для приливного движения даже очень глубокий океан можно считать мелким, поскольку его глубина всегда будет намного меньше длины волны прилива.

Генерация и распространение цунами , рассчитанные с помощью уравнений мелкой воды (красная линия; без частотной дисперсии)) и с помощью модели типа Буссинеска (синяя линия; с частотной дисперсией). Обратите внимание, что модель типа Буссинеска (синяя линия) образует солитон с осциллирующим хвостом, остающимся позади. Уравнения мелкой воды (красная линия) образуют крутой фронт, который в дальнейшем приведет к образованию борозды . Глубина воды составляет 100 метров.

Моделирование турбулентности с использованием нелинейных уравнений мелкой воды

Снимок моделирования уравнений мелкой воды, в которых присутствуют ударные волны.

Уравнения мелкой воды в своей нелинейной форме являются очевидным кандидатом для моделирования турбулентности в атмосфере и океанах, то есть геофизической турбулентности . Преимущество этого метода перед квазигеострофическими уравнениями заключается в том, что он допускает такие решения, как гравитационные волны , сохраняя при этом энергию и потенциальную завихренность . Однако с точки зрения геофизических приложений есть и некоторые недостатки - он имеет неквадратичное выражение для полной энергии и склонность волн превращаться в ударные волны . Были предложены некоторые альтернативные модели, предотвращающие образование шока. Одной из альтернатив является изменение «члена давления» в уравнении количества движения, но это приводит к сложному выражению для кинетической энергии . Другой вариант — изменить нелинейные члены во всех уравнениях, что дает квадратичное выражение для кинетической энергии , позволяет избежать образования ударной волны, но сохраняет только линеаризованную потенциальную завихренность .

См. также

• Волны и мелководье

Дальнейшее чтение

• Баттьес, JA ; Лабер, Р.Дж. (2017), Нестационарный поток в открытых каналах , Cambridge University Press, doi : 10.1017/9781316576878 , ISBN  978-1-107-15029-4
• Врейгденхил, CB (1994), Численные методы измерения потока мелкой воды , Kluwer Academic Publishers, ISBN  978-0792331643

Внешние ссылки

• Вывод уравнений мелкой воды из первых принципов (вместо упрощения уравнений Навье – Стокса некоторые аналитические решения)