Approssimazione di Boussinesq (onde marine)

Da Wikipedia, l'enciclopedia libera.
Vai alla navigazione Vai alla ricerca
Simulazione di onde periodiche al di sopra di un banco di sabbia sottomarino ottenuta con una modellazione dell'approssimazione di Boussinesq. L'esempio combina diversi effetti delle onde in acque basse, come la rifrazione, la diffrazione e una debole non linearità.

In fluidodinamica, l'approssimazione di Boussinesq per le onde marine è un'approssimazione valida per onde debolmente non lineari e abbastanza lunghe. Sono così denominate in onore del francese Joseph Boussinesq, che le derivò nel 1872[1] basandosi sulle osservazioni fatte da John Scott Russell sulle onde di traslazione, note come solitoni.

L'approssimazione di Boussinesq per le onde marine prende in considerazione la struttura verticale della velocità di flusso orizzontale e verticale. Si ottiene un'equazione differenziale alle derivate parziali di tipo non lineare che incorpora la dispersione di frequenza. Nell'ingegneria costiera le equazioni di Boussinesq vengono usate frequentemente nelle modellizzazioni al computer per la simulazione delle onde marine in acque basse e all'interno dei porti.

L'approssimazione di Boussinesq si applica alle onde abbastanza lunghe, cioè quando la lunghezza dell'onda è comparabile con la profondità dell'acqua, mentre la trattazione di Stokes è più appropriata per le onde corte, cioè quando la lunghezza dell'onda è comparabile con la profondità o è più corta.

Formalismo matematico

[modifica | modifica wikitesto]

Derivazione approssimata

[modifica | modifica wikitesto]
Onde periodiche nell'approssimazione di Boussinesq, viste in sezione verticale nella direzione di propagazione dell'onda. Si hanno ventri piatti e creste ripide a causa della non linearità. L'onda qui rappresentata in scala, ha una lunghezza d'onda di 39,1 m, un'altezza di 1,8 m (differenza di elevazione tra ventre e cresta) e una profondità media dell'acqua di 5 m. L'accelerazione di gravità è 9,81 m/s².

Il concetto fondamentale dell'approssimazione di Boussinesq è l'eliminazione della coordinata verticale dalle equazioni del flusso, pur conservando qualche influenza della struttura verticale del flusso al di sotto dell'onda marina. Questo è utile in quanto le onde si propagano nel piano orizzontale e hanno un comportamento diverso (non ondoso) nella direzione verticale. Molto spesso infatti, come nel caso di Boussinesq, l'interesse è soprattutto nella propagazione dell'onda.

L'eliminazione della coordinata verticale fu introdotta per la prima volta Joseph Boussinesq nel 1871, per trovare una soluzione approssimata per l'onda di traslazione. L'anno seguente, 1872, egli derivò le equazioni oggi note come "equazioni di Boussinesq".

I passi dell'approssimazione sono:

A questo punto si applica l'approssimazione di Boussinesq alle rimanenti equazioni di flusso, per eliminare la dipendenza dalla coordinata verticale. La risultante equazione differenziale alle derivate parziali è funzione solo delle coordinate orizzontali e del tempo.

Come esempio si consideri il flusso potenziale su un letto fluido orizzontale nel piano (x,z), cdove x e z sono rispettivamente le coordinate orizzontali e verticali. Il letto è situato a z = −h, dove h è la profondità media dell'acqua. Si fa un'espansione di Taylor del potenziale di velocità φ(x,z,t) attorno al livello del letto z = −h:[2]

dove φb(x,t) è il potenziale di velocità nel letto. Utilizzando l'equazione di Laplace per φ, valida per un flusso incomprimibile, si ottiene:

poiché la velocità verticale φ / ∂z è nulla al letto orizzontale (impermeabile) z = −h. Si può troncare la serie dopo un numero finito di termini.

Derivazione originale

[modifica | modifica wikitesto]

Per un'onda marina su un fluido incomprimibile in un flusso irrotazionale nel piano (x,z), le condizioni al contorno al livello della superficie libera z = η(x,t) sono:[3]

dove:

u è la componente orizzontale della velocità di flusso: u = ∂φ / ∂x,
w è la componente verticale della velocità di flusso: w = ∂φ / ∂z,
g è l' accelerazione di gravità.

L'approssimazione di Boussinesq per il potenziale di velocità φ, come visto sopra, si applica a queste condizioni al contorno. Inoltre nell'equazione risultante vengono mantenuti solo i termini lineari e quadratici rispetto a η e ub (con ub = ∂φb / ∂x velocità orizzontale al letto z = −h). I termini cubici e di ordine superiore si considerano trascurabili.

Si ottengono così i seguenti gruppi di equazioni:

gruppo A – Boussinesq (1872), equazione (25)

Questo gruppo di equazioni è stato derivato per un letto orizzontale piatto, cioè dove la profondità media h è costante indipendentemente dalla posizione di x. Se il secondo termine delle equazioni è posto uguale a zero, esse si riducono alle equazioni per acque basse.

Sotto opportune approssimazioni, ma con lo stesso ordine di accuratezza, il gruppo di equazioni A può essere ridotto a una singola equazione differenziale per l'elevazione della superficie libera η:

gruppo B – Boussinesq (1872), equazione (26)

Dai termini tra parentesi, l'importanza della non linearità delle equazioni può essere espressa in termini del numero di Ursell. In grandezze adimensionali, usando la profondità dell'acqua h e l'accelerazione di gravità g per l'adimensionalità, l'equazione dopo normalizzazione diventa:[4]

con:

: elevazione adimensionale della superficie,
: tempo adimensionale,
: posizione orizzontale adimensionale.

Da tale equazione, imponendo un verso preciso nella propagazione dell'onda, si può ricavare l'equazione di Korteweg-de Vries, che ammette soluzioni solitoniche.

  1. ^ J. Boussinesq, Théorie des ondes et des remous qui se propagent le long d'un canal rectangulaire horizontal, en communiquant au liquide contenu dans ce canal des vitesses sensiblement pareilles de la surface au fond, in Journal de Mathématiques Pures et Appliquées, Deuxième Série, vol. 17, 1872, pp. 55–108.
  2. ^ Dingemans (1997), p. 477.
  3. ^ Dingemans (1997), p. 475.
  4. ^ Johnson (1997), p. 219

Altri progetti

[modifica | modifica wikitesto]