Teoria lineare del moto ondoso

La teoria lineare del moto ondoso (Teoria di Airy) permette una trattazione semplificata della cinematica delle particelle fluide in posizione variata rispetto al livello indisturbato giungendo a relazioni per il calcolo della lunghezza e del periodo dell'onda e la sua elevazione dalla superficie libera.

Dimostrazione

Impostazione del problema

Si consideri il sistema di riferimento cartesiano posizionato sulla superficie libera, con asse z verticale diretto verso l'alto e asse y normale al piano. Si definisce la profondità locale del fondale h(x) la distanza tra il fondale e la superficie libera (il "tirante" idrico, costante a meno di variazioni locali qui trascurate) e l'elevazione della superficie libera (η(x,t)) la distanza tra la superficie libera ed il livello indisturbato, concorde con l'asse z. Si introduce quindi la pressione P(x,z,t) e il campo di velocità istantanea dell'acqua V(x,z,t) avente componenti lungo x e z rispettivamente chiamate u(x,z,t) e w(x,z,t).

Equazioni e ipotesi

• Principio di conservazione della massa

Si fa l'ipotesi che il fluido sia incomprimibile a causa delle modeste sovrapressioni presenti, pertanto si ha che $\mathrm{\nabla }\cdot \overrightarrow{V}=0$ ovvero $u_x+w_z=0$ (con il pedice qui si indica la variabile rispetto alla quale si sta operando la derivata parziale della funzione).

• Principio di conservazione della quantità di moto

$\rho \frac{D\overrightarrow{V}}{Dt}=\mathrm{\Sigma }{F}_m+\mathrm{\Sigma }{F}_s$, ovvero la variazione della quantità di moto è pari alla somma delle forze di massa e di superficie su un volume di controllo. Esplicitando si ha che:

* ∑ F_m: Forza di gravità + Forza di Coriolis, qui trascurabile in quanto non consideriamo domini geografici estesi

* ∑ F_s: Forze normali, di pressione, + Forze tangenziali legate alla viscosità. Si trascura in questa analisi l'attrito aria/acqua, acqua/acqua, acqua/fondo, quindi si trascura la perdita di energia col fondo e l'azione tangenziale del vento, allora la seconda ipotesi è quella di considerare il fluido perfetto.

Esplicitando i termini si ha che:

$\{\begin{array}{c}{F}_p=-\mathrm{\nabla }P\\ F_g=\hat{i}0+\hat{k}(-\rho g)\end{array}$

quindi scomponendo si ottiene:

$\{\begin{array}{c}{u}_t+u_{x}u+u_{z}w=-\frac{P_x}{\rho }\\ w_t+w_{x}u+w_{z}w=-\frac{P_z}{\rho }-g\end{array}$.

Si fa ora l'ipotesi di vorticità nulla (valida se l'onda non frange), ovvero - considerando solo l'asse y - $-w_x+u_z=0$. Da ciò ne consegue che il moto è irrotazionale quindi il campo di velocità ammette potenziale scalare φ noto in tutto il dominio che descrive il campo cinematico, tale che ${\varphi }_x=u$ e ${\varphi }_z=w$.

L'equazione di continuità allora si trasforma nel laplaciano ${\mathrm{\nabla }}^2\varphi =0$, mentre riscrivendo le equazioni di bilancio si ha:

$\{\begin{array}{c}{\varphi }_t+\frac{u^2}{2}+\frac{w^2}{2}+\frac{P}{\rho }+gz=C_1(z,t)\\ {\varphi }_t+\frac{u^2}{2}+\frac{w^2}{2}+\frac{P}{\rho }+gz=C_2(x,t)\end{array}$

Il primo termine delle equazioni è identico per entrambe, quindi le costanti dipendono solo dal tempo. Si ottiene allora l'Equazione di Bernoulli:

${\varphi }_t+\frac{u^2}{2}+\frac{w^2}{2}+\frac{P}{\rho }+gz=C(t)$

Condizioni al contorno

• Fondale:

Il fondale è descritto da tutti i punti che verificano $F:z+h=0$. Si impone che sia impermeabile e orizzontale, ovvero che le particelle si muovano solo in direzione orizzontale (altrimenti lascerebbero un vuoto o si creerebbe un vortice). Ne consegue che $\frac{dF}{dt}=0$, ovvero $w+h_{x}u=0$ (condizione cinematica).

• Superficie libera:

La superficie libera è descritta da tutti i punti che verificano $F:z-\eta =0$ e anche qui si pone $\frac{dF}{dt}=0$, ovvero ${\varphi }_z-{\eta }_x{\varphi }_x-{\eta }_t=0$. Essendo in superficie la pressione pari a quella atmosferica (P = P_atm = 0), sostituendo in Bernoulli si ha:

${\varphi }_t+\frac{1}{2}(u^2+w^2)+g\eta =0$ (condizione dinamica).

• Teoria delle onde progressive monocromatiche su fondale costante:

Si considera che le onde mantengono costante la forma e quindi si impone la condizione di periodicità, ovvero ciò che accade su un contorno è identico a ciò che accade sull'altro, quindi $\varphi (x+L,z,T)=\varphi (x,z,t)$.

Semplificazioni

Le equazioni che governano il fenomeno non sono lineari. Si pone quindi la fondamentale ipotesi che le onde abbiano piccola ampiezza rispetto alla profondità e alla lunghezza, ovvero $\frac{a}{L},\frac{a}{h}<<1$ ovvero la ripidità h/L (data anche da η_x) è molto bassa. Dalla condizione sulla superficie libera si può trascurare il prodotto ${\varphi }_x{\eta }_x$ e ${\eta }_t={\varphi }_z$ (acqua va più veloce se onda è più grande a parità di periodo). Dalla condizione dinamica, essendo u e v dipendenti da a, la loro potenza è trascurabile. In z=0 quindi si ha:

$\{\begin{array}{c}{\eta }_t-{\varphi }_z=0\\ {\varphi }_t+g\eta =0\end{array}$

dove però non è conosciuta la η, resa indistinguibile dal livello indisturbato.

Soluzione

La soluzione del problema è:

$\varphi =\frac{ag}{\omega }\frac{\cosh k(h+z)}{\cosh k(h)}\sin kx-\omega t$

con $k=\frac{2\pi }{L}$ e $\omega =\frac{2\pi }{T}$, dove L è la lunghezza dell'onda e T il suo periodo.

Verifiche delle ipotesi

• Laplace:

${\varphi }_{xx}=-\frac{k^{2}ag}{\omega }\frac{\cosh k(h+z)}{\cosh k(h)}\sin kx-\omega t=-{\varphi }_{zz}=\frac{k^{2}ag}{\omega }\frac{\cosh k(h+z)}{\cosh k(h)}\sin kx-\omega t$ ⇒ Verificato

• Fondale:

${\varphi }_z=0$, verificata in quanto $\sinh k(h-h)=0$

• Superficie libera:

Dinamica: ${\varphi }_t+g\eta =0$. Da questa si dimostra che l'elevazione della superficie libera è descritta da una funzione periodica di altezza a: $\eta =a\cos kx-\omega t$

Cinematica: ${\varphi }_x+{\eta }_t=0$. Da questa si ottiene l'importantissima relazione di dispersione ${\omega }^2=gk\tanh kh$ che lega tra loro L,h e T. Esplicitando rispetto a L si ha $L=\frac{g{T}^2}{2\pi }\tanh 2\pi \frac{h}{L}$. Ne consegue che la celerità di fase è data da $c=\frac{gT}{2\pi }\tanh 2\pi \frac{h}{L}$.

Approssimazioni

La relazione della lunghezza è di tipo implicito, risolvibile per tentativi. È possibile però effettuare alcune semplificazioni.

Per h/L → ∝, ovvero in acqua profonda o onde corte, tanh(h/L) → 1 quindi, essendo g ≈ 9,81 m/s², $L_0=1,56T^2$ e $c_0=1,56T$.

Per h/L → 0, ovvero in acqua bassa o onde lunghe, tanh(h/L) → h/L quindi, essendo g ≈ 9,81 m/s², $L_s=\sqrt[2]{gh}T$ e $c=\sqrt[2]{gh}$.

È possibile comunque utilizzare la seguente relazione esplicita per la lunghezza, che fornisce errori dell'ordine del 5%:

$L=1,56T^2\tanh (\sinh \frac{2\sqrt[2]{h}}{T})$