Knoidal to'lqin - Cnoidal wave

AQSh armiyasi bombardimonchi samolyotlar vaqti-vaqti bilan uchib yurishadi shishiradi ga yaqin bo'lgan sayoz suvda Panama qirg'oq (1933). Knoidal to'lqinlar uchun keskin tepaliklar va juda tekis oluklar xarakterlidir.

Yilda suyuqlik dinamikasi, a knoidal to'lqin a chiziqli emas va aniq davriy to'lqin ning echimi Korteweg – de Fris tenglamasi. Ushbu echimlar Jakobi elliptik funktsiyasi cn, shuning uchun ular o'ylab topilgan cnoidal to'lqinlar. Ular tasvirlash uchun ishlatiladi sirt tortishish to'lqinlari juda uzoq to'lqin uzunligi, suv chuqurligi bilan taqqoslaganda.

Knoidal to'lqinli eritmalar tomonidan olingan Korteweg va de Fris, ularning 1895 yilgi qog'ozida ular o'zlarini taklif qilishadi tarqoq uzoq to'lqinli tenglama, endi Korteweg – de Vriz tenglamasi deb nomlanadi. Chegarasida cheksiz to'lqin uzunligi, knoidal to'lqin a ga aylanadi yolg'iz to'lqin.

The Benjamin - Bona - Maaxoni tenglamasi qisqa vaqt ichida yaxshilandito'lqin uzunligi Korteweg-de-Vriz tenglamasiga taqqoslaganda, bu yana bir yo'nalishli to'lqin tenglamasi bo'lib, knoidal to'lqin echimlariga ega. Bundan tashqari, Korteweg-de-Vriz tenglamasi ga yaqinlashgani uchun Bussinesq tenglamalari bir tomonlama ish uchun to'lqin tarqalishi, knoidal to'lqinlar - Bussinesq tenglamalarining taxminiy echimlari.

Masalan, tasvirlash uchun knoidal to'lqinli eritmalar sirt tortishish to'lqinlaridan tashqari, boshqa dasturlarda ham paydo bo'lishi mumkin ionli akustik to'lqinlar yilda plazma fizikasi.[1]

Knoidal to'lqin, o'tkirroq xarakterlanadi tepaliklar va xushomadgo'ylik oluklar a ga qaraganda sinus to'lqin. Ko'rsatilgan holat uchun elliptik parametr m = 0.9.
Kesib o'tish shishiradi, yaqin knoidal to'lqinli poezdlardan iborat. Surat g'arbiy qismida joylashgan Phares des Baleines (Whale Lighthouse) dan olingan Dele de Ré (Rhe oroli), Frantsiya, yilda Atlantika okeani.

Fon

Korteweg-de-Fris va Benjamin-Bona-Maxoni tenglamalari

Le Mexauté (1976) ga ko'ra davriy suv to'lqinlari uchun bir necha nazariyalarning amal qilish muddati.[2] Ochiq ko'k maydon knoidal to'lqinlar nazariyasining amal qilish doirasini beradi; uchun ochiq sariq Havo to'lqinlari nazariyasi; va kesilgan ko'k chiziqlar kerakli tartib bilan chegaralanadi Stoksning to'lqinlar nazariyasi. Ochiq kulrang soyalar beshinchi darajadan foydalangan holda raqamli yaqinlashuvlar oralig'ida kengayishni beradi oqim funktsiyasi nazariya, yuqori to'lqinlar uchun (H > ¼ Hbuzish).

The Korteweg – de Fris tenglamasi (KdV tenglamasi) zaif chiziqli va uzun to'lqinlarning bir yo'nalishda tarqalishini tavsiflash uchun ishlatilishi mumkin, bu erda uzun to'lqin degani: o'rtacha qatlam chuqurligi bilan solishtirganda uzoq to'lqin uzunliklariga ega - suyuqlik qatlamidagi sirt tortishish to'lqinlari. KdV tenglamasi a tarqoq ikkalasini ham o'z ichiga olgan to'lqin tenglamasi chastota tarqalishi va amplituda dispersiya effektlari. KdV tenglamasi klassik qo'llanilishida to'lqin uzunliklari uchun amal qiladi λ taxminan besh baravaridan ortiq o'rtacha suv chuqurligi h, shuning uchun λ > 5 h; va uchun davr τ dan katta bilan g ning kuchi tortishish tezlashishi.[3] KdV tenglamasining klassik to'lqinlar yaqinlashuvi doirasidagi holatini tasavvur qilish uchun u o'zini quyidagicha ajratib turadi:

  • Korteweg – de Fris tenglamasi - kuchsiz chiziqli va dispersiv to'lqinlarning oldinga tarqalishini tavsiflaydi, uzun to'lqinlar uchun λ > 7 h.
  • Sayoz suv tenglamalari - shuningdek, chiziqli emas va amplitudali dispersiyaga ega, ammo chastotali dispersiyasiz; ular juda uzun to'lqinlar uchun amal qiladi, λ > 20 h.
  • Bussinesq tenglamalari - KdV tenglamasi bilan bir xil kuchga ega (ularning klassik shaklida), lekin o'zboshimchalik yo'nalishlarida to'lqin tarqalishiga imkon beradi, shuning uchun nafaqat oldinga tarqaladigan to'lqinlar. Kamchilik shundaki, Bussinesq tenglamalarini echish ko'pincha KdV tenglamasidan ko'ra qiyinroq; va ko'plab dasturlarda to'lqin aks etishi kichik va ularni e'tiborsiz qoldirish mumkin.
  • Havo to'lqinlari nazariyasi - to'liq chastotali dispersiyaga ega, shuning uchun ixtiyoriy chuqurlik va to'lqin uzunligi uchun amal qiladi, ammo amplituda dispersiyasiz, past amplituda to'lqinlar bilan cheklangan chiziqli nazariya.
  • Stoksning to'lqinlar nazariyasi - zaif chiziqli va dispersiv to'lqinlarni tavsiflashda bezovtalanish seriyali yondashuv, ayniqsa suvning chuqurligi bilan taqqoslaganda nisbatan qisqa to'lqin uzunliklari uchun chuqurroq suvda muvaffaqiyatli. Biroq, uzoq to'lqinlar uchun Kuss tenglamasida qo'llaniladigan Bussinesq yondashuvi ko'pincha tanlanadi. Buning sababi shundaki, sayoz suvda, Stoksning bezovtalanish seriyasi, eng yuqori darajaga etganligi sababli, eritma tomon yaqinlashguncha juda ko'p atamalarga muhtoj tepaliklar va uzun tekis oluklar chiziqli bo'lmagan to'lqinlarning KdV yoki Boussinesq modellari bu uzoq chiziqli bo'lmagan to'lqinlar uchun yaxshi taxminlarni beradi.

KdV tenglamasini Bussinesq tenglamasidan olish mumkin, ammo oldinga to'lqin tarqalishini ajratish uchun qo'shimcha taxminlar zarur. Amaliy qo'llanmalar uchun Benjamin - Bona - Maaxoni tenglamasi (BBM tenglamasi) KdV tenglamasidan afzalroq, KdV ga o'xshash oldinga tarqaladigan model, ammo qisqa to'lqin uzunliklarida chastotali-dispersiyaviy xatti-harakatlar yaxshi. Qisqa to'lqinli ishlashni yanada takomillashtirishni zamonaviy takomillashtirilgan Boussinesq modelidan to'lqin uzunliklari uchun ham amal qiladigan bir tomonlama to'lqinli tenglamani chiqarishni boshlash orqali olish mumkin.[4]

Knoidal to'lqinlar

Elliptik parametrning uchta qiymati uchun knoidal to'lqin rejimlari m.
ko'k : m = 0,
qizil : m = 0,9 va
qora : m = 0.99999.

KdV tenglamasining knoidal to'lqinli echimlari Korteveg va de Vriz tomonidan 1895 yilda nashr etilgan, bu maqola 1894 yilda de Vriz tomonidan doktorlik dissertatsiyasiga asoslangan.[5] Lineer bo'lmagan va dispersiv uzun to'lqinlar uchun yagona to'lqinli eritmalar ilgari topilgan Bussinesq 1872 yilda va Reyli 1876 ​​yilda. Ushbu echimlarni izlash bunga oid kuzatuvlar bilan boshlangan yolg'iz to'lqin (yoki "tarjima to'lqini") tomonidan Rassel, ham tabiatda, ham laboratoriya tajribalarida.[4] KdV tenglamasining knoidal to'lqinli echimlari kichik buzilishlarga nisbatan barqarordir.[6]

Sirt balandligi η(x,t), gorizontal holat funktsiyasi sifatida x va vaqt t, knoidal to'lqin uchun quyidagilar beriladi:[7]

qayerda H bo'ladi to'lqin balandligi, λ bo'ladi to'lqin uzunligi, v bo'ladi o'zgarishlar tezligi va η2 bo'ladi truba balandlik. Bundan tashqari, cn ulardan biri Jakobi elliptik funktsiyalari va K(m) bo'ladi birinchi turdagi to'liq elliptik integral; ikkalasi ham elliptik parametrga bog'liq m. Keyingisi, m, knoidal to'lqin shaklini aniqlaydi. Uchun m nolga teng knoidal to'lqin a ga aylanadi kosinus funktsiyasini bajaradi, biriga yaqin qiymatlar uchun esa knoidal to'lqin eng yuqori darajaga ko'tariladi tepaliklar va (juda) tekis oluklar. Ning qiymatlari uchun m 0,95 dan kam bo'lsa, knoidal funktsiyani trigonometrik funktsiyalar bilan taxmin qilish mumkin.[8]

Lineer bo'lmagan uzun to'lqinlar uchun muhim o'lchovsiz parametr (λ ≫ h) bo'ladi Ursell parametri:

Ning kichik qiymatlari uchun U, demoq U < 5,[9] chiziqli nazariyadan foydalanish mumkin, va yuqori qiymatlarda knoidal to'lqin nazariyasi kabi chiziqli bo'lmagan nazariyalardan foydalanish kerak. Uchinchi yoki beshinchi tartib - Stoks va knoidal to'lqinlar nazariyalari orasidagi demarkatsiya zonasi Ursell parametrining 10-25 oralig'ida.[10] Ursell parametri formulasidan ko'rinib turibdiki, ma'lum bir to'lqin balandligi uchun H/h Ursell parametri va shu bilan bir qatorda nochiziqlik ham nisbiy to'lqin uzunligining oshishi bilan tez o'sib boradi λ/h.

Ichidagi tortishish to'lqinlarining to'la chiziqli bo'lmagan muammosini tahlil qilish asosida potentsial oqim nazariyasiga ko'ra, yuqoridagi knoidal to'lqinlarni bezovtalanish qatoridagi eng past tartibli muddat deb hisoblash mumkin. Yuqori darajali knoidal to'lqinlar nazariyalari qisqa va ko'proq chiziqli to'lqinlar uchun amal qiladi. Beshinchi darajali knoidal to'lqin nazariyasi 1979 yilda Fenton tomonidan ishlab chiqilgan.[11] Beshinchi darajali Stoklar va beshinchi darajali knoidal to'lqinlar nazariyalarining batafsil tavsifi va taqqoslanishi Fentonning sharh maqolasida keltirilgan.[12]

Renoidalizatsiya orqali knoidal to'lqinlarning tavsiflari, shuningdek chuqur suvdagi to'lqinlarga, hatto cheksiz suv chuqurligiga juda mos keladi; Klamond tomonidan topilganidek.[13][14] Haqiqiy dengizlarda topilgan sayoz suvdagi knoidal to'lqinlarning o'zaro ta'sirining tavsifi 1994 yilda Osborne tomonidan taqdim etilgan.[15]

Yuzaki taranglik

Agar sirt tarangligi effektlari (shuningdek) muhim bo'lsa, ularni uzun to'lqinlar uchun knoidal to'lqin eritmalariga kiritish mumkin.[16]

Davriy to'lqinli eritmalar

Korteweg – de Fris tenglamasi

The Korteweg – de Fris tenglamasi (KdV tenglamasi), suv to'lqinlari uchun va o'lchovli shaklda ishlatiladigan:[17]

qayerda

η: sirt balandlik, funktsiyasi x va t, ijobiy yo'nalish yuqoriga qarab (qarama-qarshi tortishish kuchi),
x: gorizontal koordinata,
t: vaqt,
g: qiymati Yerning tortishish kuchi,
h: the anglatadi suv chuqurligi va
x va ∂t: qisman lotin nisbatan operatorlar x va t.
  Xulosa tafsilotlari
O'lchovsizlashtirish

Barcha miqdorlarni bajarish mumkin o'lchovsiz tortishish tezlanishidan foydalangan holda g va suv chuqurligi h:

    va  

Natijada hosil bo'lgan KdV tenglamasining o'lchovsiz shakli[17]

Qolgan qismida tillar yozuvlar qulayligi uchun tashlanadi.

Standart shaklga aloqadorlik

Shakl

transformatsiya orqali olinadi

    va  

ammo ushbu shakl bundan keyin ushbu hosilada ishlatilmaydi.

Ruxsat etilgan shaklda tarqaladigan to'lqinlar

Vaqti-vaqti bilan to'lqinli echimlar, bilan sayohat qilish o'zgarishlar tezligi v, qidirilmoqda. Ushbu doimiy to'lqinlar quyidagilardan iborat bo'lishi kerak:

  bilan The to'lqin fazasi:  

Natijada, makon va vaqtga nisbatan qisman hosilalar quyidagicha bo'ladi:

  va  

qayerda η ' belgisini bildiradi oddiy lotin ning η(ξ) ga nisbatan dalil ξ.

KdV tenglamasida ulardan foydalanib, quyidagi uchinchi tartib oddiy differentsial tenglama olinadi:[18]

Birinchi tartibli oddiy differentsial tenglamaga integratsiya

Bu bo'lishi mumkin birlashtirilgan bir marta olish uchun:[18]

bilan r an integratsiya doimiysi. 4 bilan ko'paytirilgandan so'ngη 'va yana bir bor integratsiya[18]

Kubik polinom f (η) ning davriy to'lqinli eritmalarida uchraydigan kabi Korteweg – de Fris tenglamasi va Benjamin - Bona - Maaxoni tenglamasi.

bilan s boshqa integratsiya doimiysi. Bu shaklda yozilgan

  bilan  

 

 

 

 

(A)

Kubik polinom f(η) ning katta ijobiy qiymatlari uchun salbiy bo'ladi η, ning katta salbiy qiymatlari uchun ijobiy η. Sirt balandligidan η bu haqiqiy qadrlanadi, shuningdek, integratsiya barqarorlari r va s haqiqiydir. Polinom f bilan ifodalanishi mumkin ildizlar η1, η2 va η3:[7]

 

 

 

 

(B)

Chunki f(η) haqiqiy qadrlanadi, uchta ildiz η1, η2 va η3 yoki uchalasi ham haqiqiydir, aks holda bittasi haqiqiy, qolgan ikkitasi esa juftlik murakkab konjugatlar. Ikkinchi holatda, faqat bitta haqiqiy qiymatga ega ildiz bilan bitta balandlik mavjud η unda f(η) nolga teng. Binobarin, sirt faqat bitta balandlik Nishab η ' nolga teng. Biroq, biz to'lqinlarni ikkita balandlik bilan to'lqinlar kabi qidirmoqdamiz tepalik va truba (fizika) - bu erda sirt qiyaligi nolga teng. Xulosa shuki, uchta ildiz ham f(η) haqiqiy qiymatga ega bo'lishi kerak.

Umumiylikni yo'qotmasdan, uchta haqiqiy ildiz quyidagicha tartiblangan deb taxmin qilinadi.

Birinchi tartibli oddiy-differentsial tenglamaning echimi

Endi, tenglamadan (A) faqat qiyalik uchun haqiqiy qiymatlar mavjudligini ko'rish mumkin, agar f(η) ijobiy. Bu bilan mos keladi η2 ≤ η≤ η1, shuning uchun bu sirt balandligi tebranadigan oraliqdir, shuningdek, ning grafigiga qarang f (η). Ushbu holat balandlikning quyidagi tasviri bilan qondiriladi η(ξ):[7]

 

 

 

 

(C)

izlanayotgan to'lqin echimlarining davriy xarakteriga muvofiq va ψ(ξ) bosqichi trigonometrik funktsiyalar gunoh va cos. Ushbu shakldan, turli xil atamalarning quyidagi tavsiflari tenglamalarda (A) va (B) olish mumkin:

Bularni tenglamalarda ishlatish (A) va (B), quyidagi oddiy differentsial tenglama ψ va ξ ba'zi manipulyatsiyalardan so'ng olinadi:[7]

chunki o'ng qo'li hali ham ijobiy, chunki η1 − η3 ≥ η1 − η2. Umumiylikni yo'qotmasdan, biz buni taxmin qilishimiz mumkin ψ(ξ) - bu monoton funktsiya, chunki f(η) oralig'ida nolga ega emas η2 < η < η1. Demak, yuqoridagi oddiy differentsial tenglama ham nuqtai nazaridan echilishi mumkin ξ(ψ) ning funktsiyasi bo'lish ψ:[7]

bilan:

  va  

qayerda m - elliptik parametr,[19][20] qoniqarli 0 ≤ m ≤ 1 (chunki η3 ≤ η2 ≤ η1). Agar ξ = 0 to'lqin tepasida tanlanadi η(0) = η1 integratsiya beradi[7]

 

 

 

 

(D.)

bilan F(ψ|m) birinchi turdagi to'liq bo'lmagan elliptik integral. The Jakobi elliptik funktsiyalari cn va sn ning teskari tomonlari F(ψ|m) tomonidan berilgan

  va  

Tenglama yordamida (C), natijada KdV tenglamasining knoidal to'lqinli eritmasi topildi[7]

Qolgan narsa, parametrlarni aniqlash: η1, η2, Δ va m.

Knoidal-to'lqin parametrlari o'rtasidagi munosabatlar

Birinchidan, beri η1 tepalik balandligi va η2 truba balandligi, uni joriy etish qulay to'lqin balandligi sifatida belgilanadi H = η1 − η2. Binobarin, biz buni topamiz m va uchun Δ:

  va     shunday  

Knoidal to'lqin eritmasi quyidagicha yozilishi mumkin:

Ikkinchidan, truba joylashgan ψ = ½ π, shuning uchun orasidagi masofa ξ = 0 va ξ = ½ λ bilan, bilan λ The to'lqin uzunligi, tenglamadan (D.):

  berib  

qayerda K(m) bo'ladi birinchi turdagi to'liq elliptik integral. Uchinchidan, to'lqin suvning o'rtacha chuqurligi atrofida tebranishi sababli o'rtacha qiymati η(ξ) nolga teng bo'lishi kerak. Shunday qilib[7]

qayerda E(m) bo'ladi ikkinchi turdagi to'liq elliptik integral. Uchun quyidagi iboralar η1, η2 va η3 elliptik parametr funktsiyasi sifatida m va to'lqin balandligi H natija:[7]

    va  

To'rtinchidan, tenglamalardan (A) va (B) o'rtasida munosabatlar o'rnatilishi mumkin o'zgarishlar tezligi v va ildizlar η1, η2 va η3:[7]

Nisbatan faza tezligining o'zgarishi quyidagi rasmda tasvirlangan. Ko'rinib turibdiki, uchun m > 0,96 (shuning uchun 1 - uchunm <0.04) faza tezligi to'lqin balandligi oshishi bilan ortadi H. Bu uzunroq va ko'proq chiziqsiz to'lqinlarga mos keladi. Faza tezligining chiziqli bo'lmagan o'zgarishi, sobit uchun m, to'lqin balandligi bilan mutanosib H. Faza tezligini unutmang v to'lqin uzunligi bilan bog'liq λ va davr τ kabi:

Eritmaning xulosasi

Bu erdagi barcha miqdorlar, ularning amal qilishlari uchun o'lchovli shakllarida beriladi sirt tortishish to'lqinlari oldin o'lchovsizlashtirish.

Nisbiy o'zgarishlar tezligi uchun knoidal to'lqin eritmalarining ko'payishi Korteweg – de Fris tenglamasi 1− funktsiyasi sifatidam, bilan m elliptik parametr.
Gorizontal o'qi a ga to'g'ri keladi logaritmik o'lchov, 10 dan−6 10 ga0=1.
Shakl o'lchovsiz miqdorlar uchun, ya'ni o'zgarishlar tezligi v sayoz suv fazasi tezligi bilan o'lchovsiz qilingan va to'lqin balandligi H suvning o'rtacha chuqurligi bilan o'lchovsiz qilingan h.

KdV tenglamasining knoidal to'lqinli echimi:[7]

bilan H The to'lqin balandligi - orasidagi farq tepalik va truba balandlik, η2 truba balandligi, m elliptik parametr, v The o'zgarishlar tezligi va ulardan biri Jakobi elliptik funktsiyalari. Chuqurlik darajasi η2 va kenglik parametri Δ bilan ifodalanishi mumkin H, h va m:[7]

  va  

bilan K(m) birinchi turdagi to'liq elliptik integral va E(m) ikkinchi turdagi to'liq elliptik integral. Yozib oling K(m) va E(m) bu erda elliptik parametr funktsiyasi sifatida belgilanadi m va elliptik modulning funktsiyasi sifatida emas k, bilan m = k2.

The to'lqin uzunligi λ, faza tezligi v va to'lqin davr τ bilan bog'liq H, h va m tomonidan:[7]

    va  

bilan g The Yerning tortishish kuchi.

Ko'pincha ma'lum to'lqin parametrlari to'lqin balandligi H, suv chuqurligi degani h, tortishish tezlashishi g, yoki to'lqin uzunligi λ yoki boshqa davr τ. Keyin yuqoridagi munosabatlar λ, v va τ elliptik parametrni topish uchun ishlatiladi m. Bu talab qiladi raqamli echim kimdir tomonidan takroriy usul.[3]

Benjamin - Bona - Maaxoni tenglamasi

The Benjamin - Bona - Maaxoni tenglamasi (BBM tenglamasi) yoki muntazam uzun to'lqin (RLW) tenglamasi o'lchovli shaklda berilgan:[21]

Barcha kattaliklar KdV tenglamasi bilan bir xil ma'noga ega. BBM tenglamasi ko'pincha KdV tenglamasidan afzalroq, chunki u qisqa to'lqinli harakatga ega.[21]

  Xulosa tafsilotlari
Hosil qilish

Chiqish KdV tenglamasi uchun o'xshashdir.[22] O'lchamsiz BBM tenglamasi o'rtacha suv chuqurligidan foydalanib o'lchovsiz h va tortishish tezlashishi g:[21]

Bu standart shaklga keltirilishi mumkin

o'zgartirish orqali:

    va  

ammo bu standart shakl bu erda ishlatilmaydi.

KdV tenglamasi uchun knoidal to'lqin eritmasining chiqarilish analogi, davriy to'lqin eritmalari η(ξ) bilan ξ = xct Bu holda BBM tenglamasi uchinchi darajali oddiy differentsial tenglamaga aylanadi, uni ikki marta integrallash mumkin:

  bilan  

Bu faqat KdV tenglamasi uchun tenglamadan faktor orqali farq qiladi v ni oldida (η ′)2 chap tomonda. Koordinatali transformatsiya orqali β = ξ /  omil v olib tashlanishi mumkin, natijada KdV va BBM tenglamalari uchun bir xil birinchi darajali oddiy differentsial tenglama bo'ladi. Biroq, bu erda oldingi tenglamada berilgan shakl ishlatiladi. Buning uchun boshqa formulalar paydo bo'ladi Δ KdV tenglamasi uchun topilgan:

To'lqin uzunligining aloqasi λfunktsiyasi sifatida H va m, o'zgarishi ta'sir qiladi

Qolganlari uchun, KdV tenglamasi uchun o'xshashdir va bu erda takrorlanmaydi.

Rezyume; qayta boshlash

Natijalar chuqurlikdagi suyuq qatlamda suv to'lqinlari uchun o'lchovli shaklda keltirilgan h.

BBM tenglamasining knoidal to'lqinli eritmasi va parametrlari uchun bog'liq munosabatlar quyidagicha:[22]

KdV tenglamasining knoidal to'lqinli eritmasi bilan farq faqat uchun tenglamada to'lqin uzunligi λ.[22] Amaliy qo'llanmalar uchun odatda suv chuqurligi h, to'lqin balandligi H, tortishish tezlashishi gva ham to'lqin uzunligi λ, yoki - ko'pincha - davr (fizika) τ taqdim etiladi. Keyin elliptik parametr m uchun yuqoridagi munosabatlardan aniqlanishi kerak λ, v va τ ba'zilari orqali takroriy usul.[3]

Misol

Korteweg-de-Vriz tenglamasining knoidal to'lqinli echimlari uchun parametr munosabatlari. −log ko'rsatilgan10 (1−m) bilan m ning elliptik parametri to'liq elliptik integrallar,[20] o'lchovsiz funktsiya sifatida davr τ g/h va nisbiy to'lqin balandligi H / h. Kontur chiziqlari bo'yicha qiymatlar −log10 (1−m), shuning uchun 1 qiymati mos keladi m = 1 − 10−1 = 0,9 va qiymati 40 bilan m = 1 − 10−40.

Ushbu misolda Korteweg-de-Vriz (KdV) tenglamasiga muvofiq knoidal to'lqin ko'rib chiqilgan. To'lqinning quyidagi parametrlari berilgan:

Davr o'rniga τ, boshqa holatlarda to'lqin uzunligi λ oldindan ma'lum bo'lgan miqdor sifatida yuzaga kelishi mumkin.

Birinchidan, o'lchovsiz davr hisoblanadi:

bu yettidan kattaroq, shuning uchun knoidal nazariya amal qilishi uchun etarli. Asosiy noma'lum - bu elliptik parametr m. Bu shunday to'lqin davri aniqlanishi kerak τ, KdV tenglamasi uchun knoidal to'lqin nazariyasidan kelib chiqqan holda:

    va  

ning berilgan qiymatiga mos keladi τ; Bu yerga λ to'lqin uzunligi va v bo'ladi o'zgarishlar tezligi to'lqinning Bundan tashqari, K(m) va E(m) bor to'liq elliptik integrallar navbati bilan birinchi va ikkinchi turdagi. Elliptik parametrni qidirish m tomonidan amalga oshirilishi mumkin sinov va xato yoki raqam yordamida ildiz topish algoritmi. Bunday holda, dastlabki taxminlardan boshlab minit = 0.99, sinov va xato bilan javob

topildi. Jarayon davomida to'lqin uzunligi λ va o'zgarishlar tezligi v hisoblab chiqilgan:

  • to'lqin uzunligi λ = 50,8 m (167 fut) va
  • o'zgarishlar tezligi v = 7,26 m / s (23,8 fut / s).

Faza tezligi v qiymati bilan solishtirish mumkin ga ko'ra sayoz suv tenglamalari:

chiziqli bo'lmagan ta'sir tufayli 3,8% o'sishni ko'rsatmoqda amplituda tarqalish, bu holda faza tezligini kamaytirishdan yutadi chastota tarqalish.

Endi to'lqin uzunligi ma'lum, Ursell raqami hisoblash mumkin:

bu kichik emas, shuning uchun chiziqli to'lqinlar nazariyasi amal qilmaydi, lekin knoidal to'lqin nazariyasi amal qiladi. Va nihoyat, to'lqin uzunligining chuqurlikka nisbati λ / h = 10.2> 7, yana bu to'lqinning knoidal to'lqin deb hisoblash uchun etarlicha uzunligini ko'rsatmoqda.

Yagona to'lqin chegarasi

Parametr bilan juda uzun chiziqli bo'lmagan to'lqinlar uchun m biriga yaqin, m → 1 ga binoan, Jacobi elliptik funktsiyasi cn ga yaqinlashishi mumkin[23]

  bilan  

Bu erda sinx, cosh, tanh va sech mavjud giperbolik funktsiyalar. Chegarada m = 1:

sech bilan (z) = 1 / cosh (z).

Bundan tashqari, xuddi shu chegara uchun m → 1, birinchi turdagi to'liq elliptik integral K(m) cheksizlikka boradi, shu bilan birga ikkinchi turdagi to'liq elliptik integral E(m) biriga boradi.[24] Bu faza tezligining chegara qiymatlarini nazarda tutadi v va minimal daraja η2 bo'lish:[25]

  va  

Binobarin, kenglik parametri bo'yicha Δ, yolg'iz to'lqin ikkala KdV va BBM tenglamasining echimi:[25]

Kenglik parametri, knoidal to'lqinlar uchun topilgan va endi chegarada m → 1, KdV va BBM tenglamasi uchun farq qiladi:[25]

 : KdV tenglamasi va
 : BBM tenglamasi.

Ammo ikkala tenglamadagi yakka to'lqinning faza tezligi bir xil, balandlikning ma'lum birikmasi uchun H va chuqurlik h.

Infinitesimal to'lqin balandligining chegarasi

Uchun cheksiz to'lqinlar balandligi knoidal to'lqinlar nazariyasining natijalari bilan yaqinlashishi kutilmoqda Havo to'lqinlari nazariyasi uzoq to'lqinlar chegarasi uchun λ ≫ h. Avval cheksiz kichik to'lqin balandligi uchun knoidal to'lqinlarning sirt balandligi va undan keyin faza tezligi ko'rib chiqiladi.

Yuzaki balandlik

  Xulosa tafsilotlari

The Jakobi elliptik funktsiyasi cn ni a ga kengaytirish mumkin Fourier seriyasi[26]

K '(m) xayoliy chorak davri sifatida tanilgan, while K(m) Jakobi elliptik funktsiyasining haqiqiy chorak davri deb ham ataladi. Ular quyidagilar bilan bog'liq: K '(m) = K(1−m)[27]

Bu erda qiziqish kichik parametrga mos keladigan kichik to'lqin balandligida m ≪ 1, ni ko'rib chiqish qulay Maklaurin seriyasi tegishli parametrlar uchun bilan boshlash kerak to'liq elliptik integrallar K va E:[28][29]

Keyin Furye qatorida paydo bo'ladigan giperbolik-kosinus atamalari kichraytirishi mumkin m ≪ 1 quyidagicha:[26]

  nome bilan q tomonidan berilgan  

Nome q kichik uchun quyidagi xatti-harakatlarga ega m:[30]

Binobarin, amplitudalar Furye seriyasidagi birinchi atamalardan:

:  
:  
:  

Shunday qilib, uchun m ≪ 1 Jakobi elliptik funktsiyasida birinchi Furye seriyasining birinchi shartlari mavjud:

  bilan  

Va uning maydoni

Erkin sirt η(x,t) knoidal to'lqinning elliptik parametrining kichik qiymatlari uchun uning Furye qatorida ifodalanadi m. Birinchidan, cn funktsiyasining argumenti ekanligini unutmang ξ/Δva bu to'lqin uzunligi λ = 2 Δ K(m), shuning uchun:

 

Bundan tashqari, o'rtacha erkin sirt balandligi nolga teng. Shuning uchun kichik amplituda to'lqinlarning sirt balandligi

Shuningdek, to'lqin uzunligi λ elliptik parametrning Maclaurin seriyasiga kengaytirilishi mumkin m, KdV va BBM tenglamasi uchun boshqacha, ammo bu hozirgi maqsad uchun zarur emas.

Eslatma: Nolga cheklovchi xatti-harakatlar mTo'lqinning cheksiz kichik balandligida - quyidagilarni ham ko'rish mumkin.[31]
 
lekin yuqori tartibli muddat mutanosib m bu taxminiy tarkibida a mavjud dunyoviy atama, cn davrining mos kelmasligi tufayli (z|m), bu 4 ga tengK(m) va davr 2π kosinus uchun cos (z). Yuqoridagi Fourier seriyasi kichik uchun m bu kamchilikka ega emas va yordamida topilgan shakllarga mos keladi Lindstedt-Puankare usuli yilda bezovtalanish nazariyasi.

Uchun cheksiz to'lqin balandligi, chegarada m → 0 bo'lsa, erkin sirt balandligi quyidagicha bo'ladi:

  bilan  

Shunday qilib to'lqin amplituda ½ dirH, yarmi to'lqin balandligi. Bu o'rganilgan shaklga o'xshashdir Havo to'lqinlari nazariyasi, ammo knoidal to'lqinlar nazariyasi faqat to'lqin uzunligi o'rtacha suv chuqurligidan ancha uzun bo'lgan uzun to'lqinlar uchun amal qiladi.

Faza tezligi

  Xulosa tafsilotlari

Knoidal to'lqinning fazaviy tezligi, ham KdV, ham BBM tenglamasi uchun quyidagicha berilgan.[7][22]

Ushbu formulada fazaning tezligi quyidagicha funktsiya hisoblanadi to'lqin balandligi H va parametr m. Shu bilan birga, cheksiz kichik balandlikdagi to'lqinlar uchun to'lqin tarqalishini aniqlash uchun doimiy tezlikda faza tezligini harakatini aniqlash kerak. to'lqin uzunligi λ parametr chegarasida m nolga yaqinlashadi. Buni KdV va BBM tenglamalari uchun farq qiladigan to'lqin uzunligi tenglamasidan foydalanish orqali amalga oshirish mumkin:[7][22]

KdV:  
BBM :  

Introducing the relative gulchambar κh:

and using the above equations for the phase speed and wavelength, the factor H / m in the phase speed can be replaced by κh va m. The resulting phase speeds are:

KdV :  
BBM :
 

The limiting behaviour for small m can be analysed through the use of the Maklaurin seriyasi uchun K(m) va E(m),[28] resulting in the following expression for the common factor in both formulas for v:

so in the limit m → 0, the factor γ → −16. The limiting value of the phase speed for m ≪ 1 directly results.

The o'zgarishlar tezligi for infinitesimal wave height, according to the cnoidal wave theories for the KdV equation and BBM equation, are[32]

KdV:  
BBM:

bilan κ = 2π / λ The gulchambar va κh the relative wavenumber. These phase speeds are in full agreement with the result obtained by directly searching for sine-wave solutions of the linearised KdV and BBM equations. As is evident from these equations, the linearised BBM equation has a positive phase speed for all κh. On the other hand, the phase speed of the linearised KdV equation changes sign for short waves with κh > . This is in conflict with the derivation of the KdV equation as a one-way wave equation.

Direct derivation from the full inviscid-flow equations

Undular bore and yordam beradi og'ziga yaqin Araguari daryosi in north-eastern Brazil. View is oblique toward mouth from airplane at approximately 100 ft (30 m) altitude.[33]

Cnoidal waves can be derived directly from the noaniq, irrotatsion va siqilmaydigan flow equations, and expressed in terms of three invariants of the flow, as shown by Benjamin & Lighthill (1954) in their research on undular bores. A ma'lumotnoma doirasi bilan harakatlanuvchi o'zgarishlar tezligi, in which reference frame the flow becomes a barqaror oqim, the cnoidal wave solutions can directly be related to the ommaviy oqim, momentum flux va energy head oqimning. Keyingi Benjamin & Lighthill (1954) —using a oqim funktsiyasi description of this incompressible flow—the horizontal and vertical components of the flow velocity are the spatial derivatives of the stream function Ψ(ξ,z): +zΨ va -ξΨ, ichida ξ va z direction respectively (ξ = xct). The vertical coordinate z is positive in the upward direction, opposite to the direction of the gravitational acceleration, and the zero level of z is at the impermeable lower boundary of the fluid domain. While the free surface is at z = ζ(ξ); yozib oling ζ is the local water depth, related to the surface elevation η(ξ) kabi ζ = h + η bilan h the mean water depth.

In this steady flow, the tushirish Q through each vertical cross section is a constant independent of ξ, and because of the horizontal bed also the horizontal momentum flux S, ga bo'lingan zichlik r, through each vertical cross section is conserved. Further, for this inviscid and irrotational flow, Bernulli printsipi can be applied and has the same Bernoulli constant R oqim domenidagi hamma joyda. Ular quyidagicha ta'riflanadi:[34]

For fairly long waves, assuming the water depth ζ is small compared to the wavelength λ, the following relation is obtained between the water depth ζ(ξ) and the three invariants Q, R va S:[34]

 

 

 

 

(E)

This nonlinear and first-order oddiy differentsial tenglama has cnoidal wave solutions.

For very long waves of cheksiz amplituda on a fluid of depth h and with a uniform flow velocity v, the flow constants are according to the sayoz suv tenglamalari:[34]

    va  

Tenglama (E) can be brought into o'lchovsiz form by use of the discharge Q and gravitational acceleration g, and defining the critical depth hv:

bilan bog'liq critical flow orasidagi chegarani belgilash subkritik oqim va superkritik oqim (Shuningdek qarang Froude number ). Consequently, the non-dimensional form of the equation is

bilan

      va  

Hosil qilish

First eliminate the pressure p from the momentum flux S by use of the Bernoulli equation:

The streamfunction Ψ is expanded as a Maklaurin seriyasi around the bed at z = 0, and using that the impermeable bed is a streamline and the irrotationality of the flow: Ψ = 0 and ∂z2Ψ = 0 ot z = 0:[34]

bilan sizb the horizontal velocity at the bed z = 0. Because the waves are long, h ≫ λ, only terms up to z3 va ζ3 are retained in the approximations to Q va S. The momentum flux S keyin bo'ladi:[34]

The discharge Q becomes, since it is the value of the streamfunction Ψ at the free surface z = ζ:

As can be seen, the discharge Q is an O(ζ) miqdor. From this, the bed velocity is seen to be[34]

Yozib oling Q / ζ is an order one quantity. This relation will be used to replace the bed velocity sizb tomonidan Q va ζ in the momentum flux S. The following terms can be derived from it:

Consequently, the momentum flux S becomes, again retaining only terms up to proportional to ζ3:[34]

Which can directly be recast in the form of equation (E).

Potentsial energiya

The potential energy density

bilan r the fluid zichlik, is one of the infinite number of invariantlar of the KdV equation.[35] This can be seen by multiplying the KdV equation with the surface elevation η(x,t); after repeated use of the zanjir qoidasi the result is:

which is in conservation form, and is an invariant after integration over the interval of periodicity—the wavelength for a cnoidal wave. The potential energy is not an invariant of the BBM equation, but ½rg [η2 + 16 h2 (x η)2] is.[36]

Birinchidan dispersiya of the surface elevation in a cnoidal wave is computed. Yozib oling η2 = −(1/λ0λ H cn2(ξ/Δ|m) dx, cn(ξ/Δ|m) = cos ψ(ξ) va λ = 2 Δ K(m), shuning uchun[37]

The potential energy, both for the KdV and the BBM equation, is subsequently found to be[37]

The infinitesimal wave-height limit (m → 0) of the potential energy is Eqozon = 116 r g H2, which is in agreement with Havo to'lqinlari nazariyasi.[37] The wave height is twice the amplitude, H = 2a, in the infinitesimal wave limit.

Shuningdek qarang

Izohlar va ma'lumotnomalar

Izohlar

  1. ^ Nezlin, M.V. (1993), Physics of intense beams in plasmas, CRC Press, p. 205, ISBN  978-0-7503-0186-2
  2. ^ Le Méhauté, B. (1976), An introduction to hydrodynamics and water waves, Springer, ISBN  978-0-387-07232-6
  3. ^ a b v Dingemans (1997) pp. 718–721.
  4. ^ a b Dingemans (1997) pp. 689–691.
  5. ^ de Jager, E.M. (2006). "On the origin of the Korteweg–de Vries equation". arXiv:math/0602661v1.
  6. ^ Drazin, P.G. (1977), "On the stability of cnoidal waves", Mexanika va amaliy matematikaning har choraklik jurnali, 30 (1): 91–105, doi:10.1093/qjmam/30.1.91
  7. ^ a b v d e f g h men j k l m n o Dingemans (1997) 708-715 betlar.
  8. ^ Yunfeng Xu; Xiaohe Xia; Jianhua Vang (2012), "Knoidal to'lqin nazariyasida knoidal funktsiyani hisoblash va yaqinlashtirish", Kompyuterlar va suyuqliklar, 68: 244–247, doi:10.1016 / j.compfluid.2012.07.012
  9. ^ U normallashtirilganligi sababli Ursell parametri chiziqli nazariyani qachon qo'llash mumkinligini ko'rsatadi U ≪ 32 π2 / 3 ≈ 100.
  10. ^ Sorensen, R.M. (1993), Asosiy to'lqin mexanikasi: qirg'oq va okean muhandislari uchun, Wiley-Interscience, ISBN  978-0-471-55165-2, p. 61.
  11. ^ Fenton, JD (1979), "Yuqori darajali knoidal to'lqinlar nazariyasi", Suyuqlik mexanikasi jurnali, 94 (1): 129–161, Bibcode:1979JFM .... 94..129F, doi:10.1017 / S0022112079000975
  12. ^ Fenton, JD (1990), "Lineer bo'lmagan to'lqin nazariyalari", Le Mehauté, B.; Xans, D.M. (tahr.), Okean muhandisligi fanlari, Dengiz, 9A, Wiley Interscience, 3-25 betlar
  13. ^ Klamond, D. (1999), "Ixtiyoriy chuqurlikdagi gorizontal dengiz tubidagi barqaror sonli amplituda to'lqinlar", Suyuqlik mexanikasi jurnali, 398 (1): 45–60, Bibcode:1999 JFM ... 398 ... 45C, doi:10.1017 / S0022112099006151
  14. ^ Klamond, D. (2003), "Chuqur suvdagi knoidal tipdagi sirt to'lqinlari", Suyuqlik mexanikasi jurnali, 489: 101–120, Bibcode:2003JFM ... 489..101C, CiteSeerX  10.1.1.573.3434, doi:10.1017 / S0022112003005111
  15. ^ Osborne, A.R. (1994), "Sayoz suvning knoidal to'lqinlarning o'zaro ta'siri" (PDF), Geofizikadagi chiziqli bo'lmagan jarayonlar, 1 (4): 241–251, doi:10.5194 / npg-1-241-1994
  16. ^ Vanden-Bruk, J.-M.; Shen, M.C. (1983), "Sirt tarangligi bilan yakka va knoidal to'lqinlar to'g'risida eslatma", Zeitschrift für angewandte Mathematik und Physik, 34: 112–117, doi:10.1007 / BF00962619, S2CID  119997409
  17. ^ a b Dingemans (1997) 692-693 betlar.
  18. ^ a b v Dingemans (1997) p. 701.
  19. ^ Abramovits va Stegun (1965) p. 590.
  20. ^ a b Elliptik parametr m elliptik moduldan ajralib turadi k: m = k2. Abramowitz & Stegun (1965) betiga qarang. 590.
  21. ^ a b v Dingemans (1997) p. 694-696.
  22. ^ a b v d e Dingemans (1997) p. 715.
  23. ^ Abramovits va Stegun (1965) tenglama. 16.15.2, p. 574.
  24. ^ Abramovits va Stegun (1965) Shakllar 17.1 va 17.2, p. 592.
  25. ^ a b v Dingemans (1997) 702-704 betlar.
  26. ^ a b Abramovits va Stegun (1965) tenglama. 16.23.2, p. 575.
  27. ^ Abramovits va Stegun (1965) tenglama. 17.3.5, p. 590.
  28. ^ a b Dingemans (1997) p. 784.
  29. ^ Abramovits va Stegun (1965). 17.3.9 & 17.3.10, p. 591.
  30. ^ Abramovits va Stegun (1965) 17.3.21, p. 591.
  31. ^ Abramovits va Stegun (1965) tenglama. 16.13.2, p. 573.
  32. ^ Dingemans (1997) p. 695
  33. ^ 5-rasm: Syuzan Bartsch-Vinkler; Devid K. Linch (1988), Dunyo bo'ylab to'lqinlarning paydo bo'lishi va xususiyatlari katalogi (Doiraviy 1022), U. S. Geologik tadqiqotlar
  34. ^ a b v d e f g Benjamin va Lighthill (1954)
  35. ^ Dingemans (1997) 730-733 betlar.
  36. ^ Benjamin, Bona va Maoni (1972)
  37. ^ a b v Dingemans (1997) 791-794 betlar.

Adabiyotlar

Qo'shimcha o'qish