Кноидален бран

[десно|мини|300x300пкс| Бомбардерите на американската армија летаат над речиси периодично отекување во плитка вода, блиску до брегот на [Панама (1933). Острите гребени и многу рамните корита се карактеристични за кноидалните бранови.]]
Во хидродинамиката, **кноидалниот бран** е нелинеарно и точно периодично браново решение на равенката Кортевег-де Врис. Овие решенија се во однос на Јакобиевата елиптична функција *cn*, поради што се измислени *cn* кноидалните бранови. Тие се користат за да се опишат површинските гравитациски бранови со прилично долга бранова должина, во споредба со длабочината на водата.

Решенијата на кноидалните бранови биле изведени од Кортевег и де Врис, во нивниот труд од 1895 година во кој тие исто така ја предлагаат нивната дисперзивна равенка со долги бранови, сега позната како равенка Кортевег-де Вриј. Во границата на бесконечна бранова должина, кноидалниот бран станува осамен бран.

Равенката Бенџамин-Бона-Махони го подобрила однесувањето на кратки бранови должини, во споредба со равенката Кортевег-де Вриј, и е уште една еднонасочна бранова равенка со решенија на кноидни бранови. Понатаму, бидејќи равенката Кортевег-де Вриј е приближување на Бусинеск равенките за случајот на еднонасочно ширење на брановите, кноидните бранови се приближни решенија за Бусинеск равенките.

Решенијата на кноидните бранови може да се појават и во други апликации освен површинските гравитациски бранови, на пример за да се опишат јонските акустични бранови во физиката на плазмата.
десно|мини|600x600пкс| Кноиден бран, кој се одликува со поостри гребени и порамни корита отколку во синусен бран. За прикажаниот случај, елиптичниот параметар е ''m''&nbsp;=&nbsp;0,9.
десно|мини|300x300пкс

## Заднина

### Равенките Кортевег-де Врис и Бенџамин-Бона-Махони
[десно|мини|338x338пкс| Валидност на неколку теории за периодични бранови на вода Светло-сината област го дава опсегот на валидност на теоријата на кноидалните бранови; светло-жолта за [Ејрова бранова теорија|теоријата на воздушни бранови ; и испрекинати сини линии се разграничуваат помеѓу потребниот ред во теоријата на бранови на Стоукс. Светло-сивото засенчување го дава проширувањето на опсегот со нумерички апроксимации користејќи теорија на функција-проток од петти ред, за високи бранови ( *H*&nbsp;>&nbsp;¼&nbsp;H кршење ).]]
Равенката Кортевег-де Врис (КдВ равенка) може да се користи за да се опише еднонасочното ширење на слабо нелинеарни и долги бранови - каде што долгиот бран значи: има долги бранови должини во споредба со просечната длабочина на водата - на површинските гравитациски бранови на течниот слој. КдВ равенката е дисперзивна бранова равенка, вклучувајќи ги и дисперзиите на честотата и амплитудните ефекти на дисперзија. Во неговата класична употреба, КдВ равенката е применлива за бранови должини *λ* што надминуваат околу пет пати од просечната длабочина на водата *h*, така што за *λ*&nbsp;>&nbsp;5&nbsp;*ч* ; а за периодот *τ* поголем од \scriptstyle 7 \sqrt{h/g} со *g* јачината на гравитациското забрзување. За да се предвиди позицијата на КдВ равенката во опсегот на апроксимации на класичните бранови, таа се разликува на следниве начини:

* **Равенката Кортевег-де Врис** — го опишува напредното ширење на слабо нелинеарни и дисперзивни бранови, за долги бранови со *λ*&nbsp;>&nbsp;7&nbsp;*ж*.
* **Равенките за плитка вода** - исто така се нелинеарни и имаат амплитудна дисперзија, но без дисперзија на честотата; тие важат за многу долги бранови, *λ*&nbsp;>&nbsp;20&nbsp;*ж*.
* **Равенките** на Бусинеск - го имаат истиот опсег на валидност како и равенката КдВ (во нивната класична форма), но овозможуваат ширење на брановите во произволни насоки, па не само бранови што се шират напред. Недостаток е што равенките на Бусинеск често се потешки за решавање од КдВ равенката; а во многу примени рефлексиите на брановите се мали и може да се занемарат.
* **[[Ејрова бранова теорија|Теорија на воздушни бранови]]** — има целосна честотна дисперзија, толку валидна за произволна длабочина и бранова должина, но е линеарна теорија без амплитудна дисперзија, ограничена на бранови со ниска амплитуда.
* **Стоуксова бранова теорија** — пристап од низа на пертурбации за опис на слабо нелинеарни и дисперзивни бранови, особено успешен во подлабока вода за релативно кратки бранови должини, во споредба со длабочината на водата. Меѓутоа, за долгите бранови често се претпочита пристапот на равенките на Бусинеск - како што се применува и во КдВ равенката. Тоа е затоа што во плитка вода, на серијата на пертурбации на Стоукс им требаат многу термини пред да се приближи кон решението, поради големите врвови и долгите рамни корита на нелинеарните бранови..

КдВ равенката може да се изведе од равенките на Бусинеск, но потребни се дополнителни претпоставки за да може да се одвои ширењето на напредниот бран. За практични примени, равенката Бенџамин-Бона-Махони (ББМ равенка) е претпочитана во однос на КдВ равенката, модел кој се шири напред сличен на КдВ, но со многу подобро однесување на дисперзија на честотата на пократки бранови должини. Дополнителни подобрувања во перформансите на кратки бранови може да се добијат со изведување на еднонасочна бранова равенка од современ подобрен модел, валиден за уште пократки бранови должини.

### Кноидни бранови
Решенијата на кноидните бранови на КдВ равенката биле претставени од Кортевег и де Врис во нивниот труд од 1895 година, кој напис е заснован на докторската теза на Де Врис во 1894 година.  Солитарните бранови решенија за нелинеарни и дисперзивни долги бранови биле пронајдени порано од Бусинеск во 1872 година и Рејли во 1876 година. Потрагата по овие решенија била поттикната од набљудувањата на овој осамен бран од Расел, и во природата и во лабораториските експерименти. Кноидните бранови решенија на КдВ равенката се стабилни во однос на мали пертурбации.

Површинската кота *η* ( *x*, *t* ), во функција на хоризонталната положба *x* и времето *t*, за кноиден бран е дадена со:

: \eta(x,t) = \eta_2 + H\, \operatorname{cn}^2\, \left( \begin{array}{c|c} \displaystyle 2\, K(m)\, \frac{x-c\,t}{\lambda} & m \end{array} \right), 

каде што *H* е висината на бранот, *λ* е брановата должина, *c* е брзината на фазата и *η <sub>2</sub>* е надморската височина. Понатаму, cn е една од елиптичните функции на Јакоби, а *K* (*m* ) е целосниот елиптичен интеграл од првиот вид; и двете се зависни од елиптичниот параметар *m*. Последново, *m*, го одредува обликот на кноидалниот бран. За *m* еднаков на нула, кноидалниот бран станува косинус функција, додека за вредности блиску до еден, ноидалниот бран добива врвови и (многу) рамни корита. За вредности од *m* помали од 0,95, кноидалната функција може да се приближи со тригонометриски функции.

Важен бездимензионален параметар за нелинеарни долги бранови ( *λ*&nbsp;≫&nbsp;*h* ) е Урселовиот број:

: U = \frac{H\, \lambda^2}{h^3} = \frac{H}{h}\, \left( \frac{\lambda}{h} \right)^2.

За мали вредности на *U*, каде *U*&nbsp;\partial_t \eta + \sqrt{gh}\; \partial_x \eta + \tfrac{3}{2}\, \sqrt{\frac{g}{h}}\; \eta\, \partial_x \eta + \tfrac{1}{6}\, h^2\, \sqrt{gh}\; \partial_x^3 \eta = 0, 

; Non-dimensionalisation
All quantities can be made dimensionless using the gravitational acceleration *g* and water depth *h*:

:\tilde{\eta} = \frac{\eta}{h},  \tilde{x} = \frac{x}{h} и  \tilde{t} = \sqrt{\frac{g}{h}}\, t.

Добиената недимензионална форма на КдВ равенката е 

:\partial_\tilde{t} \tilde{\eta} + \partial_\tilde{x} \tilde{\eta} + \tfrac{3}{2}\, \tilde{\eta}\, \partial_\tilde{x} \tilde{\eta} + \tfrac{1}{6}\, \partial_\tilde{x}^3 \tilde{\eta} = 0,

; Врска со стандардна форма
Формата

:\partial_\hat{t} \phi + 6\, \phi\ \partial_\hat{x} \phi + \partial_\hat{x}^3 \phi =0

се добива преку трансформацијата

:\hat{t} = \tfrac16\, t, \,  \hat{x} = x - t\,  и   \phi = \tfrac32\, \eta, \,

но оваа форма повеќе нема да се користи во оваа изведба.

; Бранови кои се шират со фиксна форма
Се бараат решенија за периодични бранови, кои се движат со фазна брзина *c*. Овие постојани бранови треба да бидат од следново: 
:\eta = \eta(\xi)\,  with \xi бранова фаза:  \xi = x - c\, t.\,

Следствено, парцијалните деривати во однос на просторот и времето стануваат:

:\partial_x \eta = \eta'\,  и  \partial_t \eta = -c\, \eta', \,

каде што *η* го означува изводотна *η*(*ξ*) во однос на аргумент *ξ*.

Користејќи ги овие во КдВ равенката, се добива следнава обична диференцијална равенка од трет ред:

:\tfrac16\, \eta** + \left( 1 - c \right)\, \eta' + \tfrac32\, \eta\, \eta' = 0.\,</math>

; Интеграција во обична диференцијална равенка од прв ред
Ова може да се интегрирано еднаш, за да се добие :

:<math>\tfrac16\eta'' + \left( 1 - c \right)\, \eta + \tfrac34\, \eta^2 = \tfrac14 r, \, </math>

со ''r'' и [[константа на интеграција]]. Откако ќе се множи со 4&nbsp;''η**' и ќе се интегрира уште еднаш 

:\tfrac13 \left( \eta' \right)^2 + 2\, \left( 1 - c \right)\, \eta^2 + \eta^3 = r\, \eta + s, \,

со *s* уште една константа на интеграција. Ова е напишано во форма

Кубниот полином *f*(*η*) станува негативен за големи позитивни вредности на *η*, а позитивен за големи негативни вредности на *η*. Бидејќи површинската кота *η* е реално вреднувано, исто така и константите на интеграција *r* и *s* се реални. Полиномот *f* може да се изрази во однос на неговите корени *η<sub>1</sub>*, *η<sub>2</sub>' ' и *η3*:



Бидејќи *f*(*η*) е реално вреднувана, трите корени *η1*, *η2* и * η3* се или сите три реални, или инаку едното е реално, а преостанатите две се пар од [[комплексен конјугат]]. Во вториот случај, со само еден корен со реална вредност, постои само една елевација *η* на која *f*(*η*) е нула. И, следствено, само една кота на која површината [[наклон]] *η* е нула. Како и да е, бараме решенија слични на брановите, со две издигнувања - бранот [[сртот (физика)|сртот]] и [[коритото (физика)]] - каде што наклонот на површината е нула. Заклучокот е дека сите три корени на *f*(*η*) треба да бидат реално вреднувани.

Без губење на општоста, се претпоставува дека трите вистински корени се подредени како: 
:<math>\eta_1 \ge \eta_2 \ge \eta_3.\,</math>

; Решение на обично-диференцијална равенка од прв ред) може да се види дека постојат само реални вредности за наклонот ако *f*(*η*) е позитивно. Ова кореспондира со *η2*&nbsp;≤&nbsp;*η*≤&nbsp;*η1*, што според тоа е опсегот помеѓу висината на површината осцилира, видете го и графикот на *f(*η*). Овој услов е задоволен со следнава претстава за котата *η*(*ξ*) :

во согласност со периодичниот карактер на бараните бранови решенија и со *ψ*(*ξ*) фазата на тригонометриските функции sin и cos. Од оваа форма, следните описи на различни поими во равенки () и () може да се добијат:

:
  \begin{align}
    \eta - \eta_1 &= - \left( \eta_1 - \eta_2 \right)\; \sin^2\, \psi(\xi), \\
    \eta - \eta_2 &= + \left( \eta_1 - \eta_2 \right)\; \cos^2\, \psi(\xi), \\
    \eta - \eta_3 &= \left(\eta_1 - \eta_3 \right) - \left( \eta_1 - \eta_2 \right)\; \sin^2\, \psi(\xi), && \text{and} \\
    \eta' &= -2\, \left( \eta_1 - \eta_2 \right)\; \sin\, \psi(\xi)\; \cos\, \psi(\xi)\;\; \psi'(\xi) 
      && \text{with} \quad \psi'(\xi) = \frac{\text{d}\psi(\xi)}{\text{d}\xi}.
  \end{align}


Using these in equations () и (), се добива следната обична диференцијална равенка која се однесува на *ψ* и *ξ*, по неколку манипулации :

:\frac43\, \left( \frac{\text{d}\psi}{\text{d}\xi} \right)^2 = \left(\eta_1 - \eta_3 \right) - \left( \eta_1 - \eta_2 \right)\; \sin^2\, \psi(\xi), 

со десната страна сè уште позитивна, бидејќи *η<sub>1</sub>*&nbsp;−&nbsp;*η<sub>3</sub>*&nbsp;≥&nbsp;*η<sub >1</sub>*&nbsp;−&nbsp;*η<sub>2</sub>*. Без губење на општоста, можеме да претпоставиме дека *ψ*(*ξ*) е монотона функција, бидејќи *f*(*η*) нема нули во интервалот *η<sub >2</sub>*&nbsp;1*. Така, горната обична диференцијална равенка може да се реши и во однос на тоа што *ξ*(*ψ*) е функција од *ψ*:

:<math>\frac{1}{\Delta}\, \frac{\text{d}\xi}{\text{d}\psi} = \pm\, \frac{1}{\sqrt{1 - m \sin^2\, \psi}},</math>

со:

:<math>\Delta^2 = \frac43\, \frac{1}{\eta_1 - \eta_3}</math>  and  <math>m = \frac{\eta_1 - \eta_2}{\eta_1 - \eta_3},</math>
каде што *m* е таканаречениот елиптичен параметар, задоволува *0*&nbsp;≤&nbsp;*m*&nbsp;≤&nbsp;1 (бидејќи *η3*&nbsp;≤&nbsp;* η2*&nbsp;≤&nbsp;*η1*).

{\sqrt{1 - m\, \sin^2 \hat{\psi}}} = \pm\, \Delta\, F(\psi|m),</math>|}}

со *F*(*ψ*|*m*) [[нецелосен елипсовиден интеграл од првиот вид]]. [[Елиптичните функции Јакоби]] cn и sn се инверзи на *F*(*ψ*|*m*) дадени со

:<math>
  \cos\, \psi = \operatorname{cn} \left( \begin{array}{c|c} \displaystyle \frac{\xi}{\Delta} & m \end{array} \right)
</math>  и <math>
  \sin\, \psi = \operatorname{sn} \left( \begin{array}{c|c} \displaystyle \frac{\xi}{\Delta} & m \end{array} \right).
</math>

Со употреба на равенка (), е пронајдено добиеното кноидно-брановно решение на КдВ равенката 

:<math>
  \eta(\xi) = 
    \eta_2 
    + \left( \eta_1 - \eta_2 \right)\, \operatorname{cn}^2 \left( \begin{array}{c|c} \displaystyle \frac{\xi}{\Delta} & m \end{array} \right).
</math>

What remains, is to determine the parameters: *η1*, *η2*, *Δ* and *m*.

; Relationships between the cnoidal-wave parameters
First, since *η1* is the crest elevation and *η2* is the trough elevation, it is convenient to introduce the [[wave height]], defined as *H*&nbsp;=&nbsp;*η1*&nbsp;−&nbsp;*η2*. Consequently, we find for *m* and for *Δ*:

:<math>m = \frac{H}{\eta_1-\eta_3}</math>  и  <math>\frac{\Delta^2}{m}=\frac{4}{3\,H}</math>  па <math>\Delta=\sqrt{\frac{4}{3}\frac{m}{H}}.</math>

Решението на каноидниот бран може да се запише како:

:<math>\eta(\xi) = \eta_2 + H\, \operatorname{cn}^2 \left( \begin{array}{c|c} \displaystyle \frac{\xi}{\Delta} и m \end{array} \right).</math>

Второ, коритото се наоѓа на *ψ*&nbsp;=&nbsp;½&nbsp;*π*, така што растојанието помеѓу *ξ*&nbsp;=&nbsp;0 и *ξ*&nbsp;=&nbsp;½&nbsp;*λ* е, со *λ* [[бранова должина]], од равенката ( ):

:<math>\tfrac12\, \lambda = \Delta\, F\left( \begin{array}{c|c} \tfrac12\,\pi & m \end{array} \right) = \Delta\, K(m),</math>  давајќи <math>\lambda = 2\, \Delta\, K(m) = \sqrt{\frac{16}{3}\frac{m}{H}}\; K(m),</math>

каде што *K*(*m*) е [[целосен елипсовиден интеграл од првиот вид]]. Трето, бидејќи бранот осцилира околу средната длабочина на водата, просечната вредност на *η*(*ξ*) треба да биде нула. Значи

:<math>
  \begin{align}
    0 &= \int_0^{\lambda} \eta(\xi)\; \text{d}\xi 
       = 2\, \int_0^{\tfrac12\lambda} \left[ \eta_2 + \left( \eta_1 - \eta_2 \right)\, 
                          \operatorname{cn}^2\, \left( \begin{array}{c|c} \displaystyle \frac{\xi}{\Delta} & m \end{array} \right) \right]\; \text{d}\xi
    \\
      &= 2\, \int_0^{\tfrac12\pi} \Bigl[ \eta_2 + \left( \eta_1 - \eta_2 \right)\, \cos^2\, \psi \Bigr]\, \frac{\text{d}\xi}{\text{d}\psi}\; \text{d}\psi
       = 2\, \Delta\, \int_0^{\tfrac12\pi} \frac{\eta_1 - \left( \eta_1 - \eta_2 \right)\, \sin^2\, \psi}{\sqrt{1 - m\, \sin^2\, \psi}}\; \text{d}\psi
    \\
      &= 2\, \Delta\, \int_0^{\tfrac12\pi} \frac{\eta_1 - m\, \left( \eta_1 - \eta_3 \right)\, \sin^2\, \psi}{\sqrt{1 - m\, \sin^2\, \psi}}\; \text{d}\psi
       = 2\, \Delta\, \int_0^{\tfrac12\pi} \left[ \frac{\eta_3}{\sqrt{1 - m\, \sin^2\, \psi}} 
                                                + \left( \eta_1 - \eta_3 \right)\, \sqrt{1 - m\, \sin^2\, \psi} \right]\; \text{d}\psi
    \\
      &= 2\, \Delta\, \Bigl[  \eta_3\, K(m) + \left( \eta_1 - \eta_3 \right)\, E(m) \Bigr]
       = 2\, \Delta\, \Bigl[  \eta_3\, K(m) + \frac{H}{m}\, E(m) \Bigr],
  \end{align}
</math>

каде што *E*(*m*) е [[целосен елипсовиден интеграл од вториот вид]]. Следниве изрази за *η1*, *η2* и *η3* како функција на елиптичниот параметар *m* и висината на бранот *H* резултат:

:<math>\eta_3 = -\, \frac{H}{m}\, \frac{E(m)}{K(m)},</math>  <math>\eta_1 = \frac{H}{m}\, \left( 1 - \frac{E(m)}{K(m)} \right)</math>  and  <math>\eta_2 = \frac{H}{m}\, \left( 1 - m - \frac{E(m)}{K(m)} \right).</math>

Четврто, од равенките () and () a relationship can be established between the [[phase speed]] *c* and the roots *η1*, *η2* and *η3*:

:<math>c = 1 + \tfrac12\, \left( \eta_1 + \eta_2 + \eta_3 \right) = 1 + \frac{H}{m}\, \left( 1 - \frac12\, m - \frac32\, \frac{E(m)}{K(m)} \right). </math>

The relative phase-speed changes are depicted in the figure below. As can be seen, for *m*&nbsp;>&nbsp;0.96 (so for 1&nbsp;−&nbsp;*m*&nbsp;<&nbsp;0.04) the phase speed increases with increasing wave height *H*.  This corresponds with the longer and more nonlinear waves.  The nonlinear change in the phase speed, for fixed *m*, is proportional to the wave height *H*.  Note that the phase speed *c* is related to the wavelength *λ* and [[period (physics)|period]] *τ'' as:

:c = \frac{\lambda}{\tau}.

; Résumé of the solution
All quantities here will be given in their dimensional forms, as valid for surface gravity waves before non-dimensionalisation. 

### Равенка Бенџамин-Бона-Махони
Равенката Бенџамин-Бона-Махони (ББМ равенка), или регуларната равенка на долг бран (RLW) е во димензионална форма дадена со:

: \partial_t \eta + \sqrt{g\,h}\, \partial_x \eta + \tfrac32\, \sqrt{\frac{g}{h}}\, \eta\, \partial_x \eta - \tfrac16\, h^2\, \partial_t\, \partial_x^2 \eta = 0.

Сите величини имаат исто значење како и за КдВ равенката. Равенката BBM често се претпочита пред равенката КдВ бидејќи има подобро однесување на кратки бранови.

## Наводи

Категорија:Бранови