Earth
Maree in un estuario
Guarda la marea della sera risalire un estuario a imbuto oltre il faro: l’alveo si restringe verso il fiume, così l’onda cresce avanzando e l’escursione alla testata è molto maggiore che alla foce, il flusso dura meno del riflusso, e l’acqua salata del mare penetra decine di chilometri nell’entroterra contro il fiume. Alza la portata del fiume e il sale è respinto verso il mare; accentua l’imbuto o raddrizzalo; chiudi la testata e l’alveo risuona come un quarto d’onda, con la marea della testata a / cos(kL) di quella del mare. Acqua chiusa e calma resta esattamente calma, un estuario salato senza fiume resta esattamente salato, e acqua e sale si bilanciano con ciò che ha attraversato foce e testata.

Cosa succede
Un canale di 60 km profondo 8 m la cui larghezza si restringe da 3 km come exp(−x/25 km), in 120 celle su una griglia sfalsata: la superficie ai centri, la velocità alle facce. Passi avanti-indietro di 20 s aggiornano la superficie dalle portate alle facce, poi la velocità dalla nuova pendenza, con avvezione upwind e attrito di Manning implicito. La marea M2 del mare di 1,5 m entra dalla foce nell’arco di due periodi; il fiume porta 100 m³/s alla testata. Il sale viaggia upwind con le stesse portate con una dispersione mareale di 150 m²/s, così un estuario salato senza fiume resta salato e l’acqua chiusa e calma resta calma. Un canale prismatico chiuso verifica il quarto d’onda a/cos(kL).
Da provare
Guarda l’escursione crescere risalendo l’imbuto e il flusso durare meno del riflusso; alza il fiume e spingi fuori il sale; raddrizza o accentua l’imbuto; chiudi la testata e confronta con il quarto d’onda.
Il modello e i suoi limiti
Un canale unidimensionale ben rimescolato di profondità uniforme, niente prosciugamento delle piane di marea, salinità senza effetti di densità, una sola componente di marea. Le escursioni appartengono a questo estuario idealizzato.
Equazioni del modulo
Equazioni b·∂η/∂t + ∂(b·H·u)/∂x = 0, b = b0·e^(−x/L_b) ∂u/∂t + u·∂u/∂x = −g·∂η/∂x − g·n²·|u|·u/H^(4/3) ∂(b·H·S)/∂t + ∂(b·H·u·S − b·H·K·∂S/∂x)/∂x = 0 closed prismatic channel: a_head = a/cos(kL), k = ω/√(g·h)


