Boussinesq approximation (water waves)

Short description: Approximation valid for weakly non-linear and fairly long waves

Simulation of periodic waves over an underwater shoal with a Boussinesq-type model. The waves propagate over an elliptic-shaped underwater shoal on a plane beach. This example combines several effects of waves and shallow water, including refraction, diffraction, shoaling and weak non-linearity.

In fluid dynamics, the Boussinesq approximation for water waves is an approximation valid for weakly non-linear and fairly long waves. The approximation is named after Joseph Boussinesq, who first derived them in response to the observation by John Scott Russell of the wave of translation (also known as solitary wave or soliton). The 1872 paper of Boussinesq introduces the equations now known as the Boussinesq equations.[1]

The Boussinesq approximation for water waves takes into account the vertical structure of the horizontal and vertical flow velocity. This results in non-linear partial differential equations, called Boussinesq-type equations, which incorporate frequency dispersion (as opposite to the shallow water equations, which are not frequency-dispersive). In coastal engineering, Boussinesq-type equations are frequently used in computer models for the simulation of water waves in shallow seas and harbours.

While the Boussinesq approximation is applicable to fairly long waves – that is, when the wavelength is large compared to the water depth – the Stokes expansion is more appropriate for short waves (when the wavelength is of the same order as the water depth, or shorter).

## Contents

* 1 Boussinesq approximation
* 2 Original Boussinesq equations
* 2.1 Derivation
* 2.2 Linear frequency dispersion
* 3 Boussinesq-type equations and extensions
* 4 Further approximations for one-way wave propagation
* 5 Numerical models
* 6 Notes
* 7 References

## Boussinesq approximation

Periodic waves in the Boussinesq approximation, shown in a vertical cross section in the wave propagation direction. Notice the flat troughs and sharp crests, due to the wave nonlinearity. This case (drawn on scale) shows a wave with the wavelength equal to 39.1 m, the wave height is 1.8 m (i.e. the difference between crest and trough elevation), and the mean water depth is 5 m, while the gravitational acceleration is 9.81 m/s2.

The essential idea in the Boussinesq approximation is the elimination of the vertical coordinate from the flow equations, while retaining some of the influences of the vertical structure of the flow under water waves. This is useful because the waves propagate in the horizontal plane and have a different (not wave-like) behaviour in the vertical direction. Often, as in Boussinesq's case, the interest is primarily in the wave propagation.

This elimination of the vertical coordinate was first done by Joseph Boussinesq in 1871, to construct an approximate solution for the solitary wave (or wave of translation). Subsequently, in 1872, Boussinesq derived the equations known nowadays as the Boussinesq equations.

The steps in the Boussinesq approximation are:

* a Taylor expansion is made of the horizontal and vertical flow velocity (or velocity potential) around a certain elevation,
* this Taylor expansion is truncated to a finite number of terms,
* the conservation of mass (see continuity equation) for an incompressible flow and the zero-curl condition for an irrotational flow are used, to replace vertical partial derivatives of quantities in the Taylor expansion with horizontal partial derivatives.

Thereafter, the Boussinesq approximation is applied to the remaining flow equations, in order to eliminate the dependence on the vertical coordinate. As a result, the resulting partial differential equations are in terms of functions of the horizontal coordinates (and time).

As an example, consider potential flow over a horizontal bed in the (x,z) plane, with x the horizontal and z the vertical coordinate. The bed is located at z=−h, where h is the mean water depth. A Taylor expansion is made of the velocity potential φ(x,z,t) around the bed level z=−h:[2]

φ=φb+(z+h)[∂φ∂z]z=−h+12(z+h)2[∂2φ∂z2]z=−h+16(z+h)3[∂3φ∂z3]z=−h+124(z+h)4[∂4φ∂z4]z=−h+⋯,

where φb(x,t) is the velocity potential at the bed. Invoking Laplace's equation for φ, as valid for incompressible flow, gives:

φ={φb−12(z+h)2∂2φb∂x2+124(z+h)4∂4φb∂x4+⋯}+{(z+h)[∂φ∂z]z=−h−16(z+h)3∂2∂x2[∂φ∂z]z=−h+⋯}={φb−12(z+h)2∂2φb∂x2+124(z+h)4∂4φb∂x4+⋯},

since the vertical velocity ∂φ/∂z is zero at the – impermeable – horizontal bed z=−h. This series may subsequently be truncated to a finite number of terms.

## Original Boussinesq equations

### Derivation

For water waves on an incompressible fluid and irrotational flow in the (x,z) plane, the boundary conditions at the free surface elevation z=η(x,t) are:[3]

∂η∂t+u∂η∂x−w=0∂φ∂t+12(u2+w2)+gη=0,

where:

* u is the horizontal flow velocity component: u=∂φ/∂x,
* w is the vertical flow velocity component: w=∂φ/∂z,
* g is the acceleration by gravity.

Now the Boussinesq approximation for the velocity potential φ, as given above, is applied in these boundary conditions. Further, in the resulting equations only the linear and quadratic terms with respect to η and ub are retained (with ub=∂φb/∂x the horizontal velocity at the bed z=−h). The cubic and higher order terms are assumed to be negligible. Then, the following partial differential equations are obtained:

set A – Boussinesq (1872), equation (25)

∂η∂t+∂∂x[(h+η)ub]=16h3∂3ub∂x3,∂ub∂t+ub∂ub∂x+g∂η∂x=12h2∂3ub∂t∂x2.

This set of equations has been derived for a flat horizontal bed, i.e. the mean depth h is a constant independent of position x. When the right-hand sides of the above equations are set to zero, they reduce to the shallow water equations.

Under some additional approximations, but at the same order of accuracy, the above set A can be reduced to a single partial differential equation for the free surface elevation η:

set B – Boussinesq (1872), equation (26)

∂2η∂t2−gh∂2η∂x2−gh∂2∂x2(32η2h+13h2∂2η∂x2)=0.

From the terms between brackets, the importance of nonlinearity of the equation can be expressed in terms of the Ursell number. In dimensionless quantities, using the water depth h and gravitational acceleration g for non-dimensionalization, this equation reads, after normalization:[4]

∂2ψ∂τ2−∂2ψ∂ξ2−∂2∂ξ2(3ψ2+∂2ψ∂ξ2)=0,

with:

ψ=12ηh | : the dimensionless surface elevation,

|

τ=3tgh | : the dimensionless time, and
ξ=3xh | : the dimensionless horizontal position.
Linear phase speed squared c2/(gh) as a function of relative wave number kh.
A = Boussinesq (1872), equation (25),
B = Boussinesq (1872), equation (26),
C = full linear wave theory, see dispersion (water waves)

### Linear frequency dispersion

Water waves of different wave lengths travel with different phase speeds, a phenomenon known as frequency dispersion. For the case of infinitesimal wave amplitude, the terminology is linear frequency dispersion. The frequency dispersion characteristics of a Boussinesq-type of equation can be used to determine the range of wave lengths, for which it is a valid approximation.

The linear frequency dispersion characteristics for the above set A of equations are:[5]

c2=gh1+16k2h21+12k2h2,

with:

* c the phase speed,
* k the wave number (k=2π/λ, with λ the wave length).

The relative error in the phase speed c for set A, as compared with linear theory for water waves, is less than 4% for a relative wave number kh<π/2. So, in engineering applications, set A is valid for wavelengths λ larger than 4 times the water depth h.

The linear frequency dispersion characteristics of equation B are:[5]

c2=gh(1−13k2h2).

The relative error in the phase speed for equation B is less than 4% for kh<2π/7, equivalent to wave lengths λ longer than 7 times the water depth h, called fairly long waves.[6]

For short waves with k2h2>3 equation B become physically meaningless, because there are no longer real-valued solutions of the phase speed. The original set of two partial differential equations (Boussinesq, 1872, equation 25, see set A above) does not have this shortcoming.

The shallow water equations have a relative error in the phase speed less than 4% for wave lengths λ in excess of 13 times the water depth h.

## Boussinesq-type equations and extensions

There are an overwhelming number of mathematical models which are referred to as Boussinesq equations. This may easily lead to confusion, since often they are loosely referenced to as the Boussinesq equations, while in fact a variant thereof is considered. So it is more appropriate to call them Boussinesq-type equations. Strictly speaking, the Boussinesq equations is the above-mentioned set B, since it is used in the analysis in the remainder of his 1872 paper.

Some directions, into which the Boussinesq equations have been extended, are:

* varying bathymetry,
* improved frequency dispersion,
* improved non-linear behavior,
* making a Taylor expansion around different vertical elevations,
* dividing the fluid domain in layers, and applying the Boussinesq approximation in each layer separately,
* inclusion of wave breaking,
* inclusion of surface tension,
* extension to internal waves on an interface between fluid domains of different mass density,
* derivation from a variational principle.

## Further approximations for one-way wave propagation

While the Boussinesq equations allow for waves traveling simultaneously in opposing directions, it is often advantageous to only consider waves traveling in one direction. Under small additional assumptions, the Boussinesq equations reduce to:

* the Korteweg–de Vries equation for wave propagation in one horizontal dimension,
* the Kadomtsev–Petviashvili equation for (near uni-directional) wave propagation in two horizontal dimensions,
* the nonlinear Schrödinger equation (NLS equation) for the complex-valued amplitude of narrowband waves (slowly modulated waves).

Besides solitary wave solutions, the Korteweg–de Vries equation also has periodic and exact solutions, called cnoidal waves. These are approximate solutions of the Boussinesq equation.

## Numerical models

A simulation with a Boussinesq-type wave model of nearshore waves travelling towards a harbour entrance. The simulation is with the BOUSS-2D module of SMS. Faster than real-time simulation with the Boussinesq module of Celeris, showing wave breaking and refraction near the beach. The model provides an interactive environment.

For the simulation of wave motion near coasts and harbours, numerical models – both commercial and academic – employing Boussinesq-type equations exist. Some commercial examples are the Boussinesq-type wave modules in MIKE 21 and SMS. Some of the free Boussinesq models are Celeris,[7] COULWAVE,[8] and FUNWAVE.[9] Most numerical models employ finite-difference, finite-volume or finite element techniques for the discretization of the model equations. Scientific reviews and intercomparisons of several Boussinesq-type equations, their numerical approximation and performance are e.g. (Kirby 2003), (Dingemans 1997) and (Hamm Madsen).

## Notes

1. ↑ This paper (Boussinesq, 1872) starts with: "Tous les ingénieurs connaissent les belles expériences de J. Scott Russell et M. Basin sur la production et la propagation des ondes solitaires" ("All engineers know the beautiful experiments of J. Scott Russell and M. Basin on the generation and propagation of solitary waves").
2. ↑ Dingemans (1997), p. 477.
3. ↑ Dingemans (1997), p. 475.
4. ↑ Johnson (1997), p. 219
5. ↑ 5.0 5.1 Dingemans (1997), p. 521.
6. ↑ Dingemans (1997), p. 473 & 516.
7. ↑ "Celeria.org - Celeris Boussinesq Wave Model". http://www.celeria.org.
8. ↑ "ISEC - Models". http://isec.nacse.org/models/coulwave_download.php.
9. ↑ "James T. Kirby, Funwave program". https://www1.udel.edu/kirby/programs/funwave/funwave.html.

## References

* Boussinesq, J. (1871). "Théorie de l'intumescence liquide, applelée onde solitaire ou de translation, se propageant dans un canal rectangulaire". Comptes Rendus de l'Académie des Sciences 72: 755–759. http://gallica.bnf.fr/ark:/12148/bpt6k3029x/f759.
* Boussinesq, J. (1872). "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". Journal de Mathématiques Pures et Appliquées. Deuxième Série 17: 55–108. http://visualiseur.bnf.fr/ConsulterElementNum?O=NUMM-16416&Deb=63&Fin=116&E=PDF.
* Dingemans, M.W. (1997). Wave propagation over uneven bottoms. Advanced Series on Ocean Engineering 13. World Scientific, Singapore. ISBN 978-981-02-0427-3. http://www.worldscibooks.com/engineering/1241.html. Retrieved 2008-01-21. See Part 2, Chapter 5.
* Hamm, L.; Madsen, P.A.; Peregrine, D.H. (1993). "Wave transformation in the nearshore zone: A review". Coastal Engineering 21 (1–3): 5–39. doi:10.1016/0378-3839(93)90044-9.
* Johnson, R.S. (1997). A modern introduction to the mathematical theory of water waves. Cambridge Texts in Applied Mathematics. 19. Cambridge University Press. ISBN 0-521-59832-X.
* Kirby, J.T. (2003). "Boussinesq models and applications to nearshore wave propagation, surfzone processes and wave-induced currents". in Lakhan, V.C.. Advances in Coastal Modeling. Elsevier Oceanography Series. 67. Elsevier. pp. 1–41. ISBN 0-444-51149-0.
* Peregrine, D.H. (1967). "Long waves on a beach". Journal of Fluid Mechanics 27 (4): 815–827. doi:10.1017/S0022112067002605. Bibcode: 1967JFM....27..815P.
* Peregrine, D.H. (1972). "Equations for water waves and the approximations behind them". in Meyer, R.E.. Waves on Beaches and Resulting Sediment Transport. Academic Press. pp. 95–122. ISBN 0-12-493250-9.

* v
* t
* e

Physical oceanography

Waves|

* Airy wave theory
* Ballantine scale
* Benjamin–Feir instability
* Boussinesq approximation
* Breaking wave
* Clapotis
* Cnoidal wave
* Cross sea
* Dispersion
* Edge wave
* Equatorial waves
* Fetch
* Gravity wave
* Green's law
* Infragravity wave
* Internal wave
* Iribarren number
* Kelvin wave
* Kinematic wave
* Longshore drift
* Luke's variational principle
* Mild-slope equation
* Radiation stress
* Rogue wave
* Rossby wave
* Rossby-gravity waves
* Sea state
* Seiche
* Significant wave height
* Soliton
* Stokes boundary layer
* Stokes drift
* Stokes wave
* Swell
* Trochoidal wave
* Tsunami
* megatsunami
* Undertow
* Ursell number
* Wave action
* Wave base
* Wave height
* Wave power
* Wave radar
* Wave setup
* Wave shoaling
* Wave turbulence
* Wave–current interaction
* Waves and shallow water
* one-dimensional Saint-Venant equations
* shallow water equations
* Wind wave
* model

Circulation|

* Atmospheric circulation
* Baroclinity
* Boundary current
* Coriolis force
* Coriolis–Stokes force
* Craik–Leibovich vortex force
* Downwelling
* Eddy
* Ekman layer
* Ekman spiral
* Ekman transport
* El Niño–Southern Oscillation
* General circulation model
* Geochemical Ocean Sections Study
* Geostrophic current
* Global Ocean Data Analysis Project
* Gulf Stream
* Halothermal circulation
* Humboldt Current
* Hydrothermal circulation
* Langmuir circulation
* Longshore drift
* Loop Current
* Modular Ocean Model
* Ocean dynamics
* Ocean gyre
* Princeton ocean model
* Rip current
* Subsurface currents
* Sverdrup balance
* Thermohaline circulation
* shutdown
* Upwelling
* Whirlpool
* World Ocean Circulation Experiment

Tides|

* Amphidromic point
* Earth tide
* Head of tide
* Internal tide
* Lunitidal interval
* Perigean spring tide
* Rip tide
* Rule of twelfths
* Slack water
* Tidal bore
* Tidal force
* Tidal power
* Tidal race
* Tidal range
* Tidal resonance
* Tide gauge
* Tideline
* Theory of tides

Landforms|

* Abyssal fan
* Abyssal plain
* Atoll
* Bathymetric chart
* Coastal geography
* Cold seep
* Continental margin
* Continental rise
* Continental shelf
* Contourite
* Guyot
* Hydrography
* Oceanic basin
* Oceanic plateau
* Oceanic trench
* Passive margin
* Seabed
* Seamount
* Submarine canyon
* Submarine volcano

Plate
tectonics|

* Convergent boundary
* Divergent boundary
* Fracture zone
* Hydrothermal vent
* Marine geology
* Mid-ocean ridge
* Mohorovičić discontinuity
* Vine–Matthews–Morley hypothesis
* Oceanic crust
* Outer trench swell
* Ridge push
* Seafloor spreading
* Slab pull
* Slab suction
* Slab window
* Subduction
* Transform fault
* Volcanic arc

Ocean zones|

* Benthic
* Deep ocean water
* Deep sea
* Littoral
* Mesopelagic
* Oceanic
* Pelagic
* Photic
* Surf
* Swash

Sea level|

* Deep-ocean Assessment and Reporting of Tsunamis
* Future sea level
* Global Sea Level Observing System
* North West Shelf Operational Oceanographic System
* Sea-level curve
* Sea level rise
* World Geodetic System

Acoustics|

* Deep scattering layer
* Hydroacoustics
* Ocean acoustic tomography
* Sofar bomb
* SOFAR channel
* Underwater acoustics

Satellites|

* Jason-1
* Jason-2 (Ocean Surface Topography Mission)
* Jason-3

Related|

* Argo
* Benthic lander
* Color of water
* DSV Alvin
* Marginal sea
* Marine energy
* Marine pollution
* Mooring
* National Oceanographic Data Center
* Ocean
* Ocean exploration
* Ocean observations
* Ocean reanalysis
* Ocean surface topography
* Ocean thermal energy conversion
* Oceanography
* Pelagic sediment
* Sea surface microlayer
* Sea surface temperature
* Seawater
* Science On a Sphere
* Thermocline
* Underwater glider
* Water column
* World Ocean Atlas

* Category
* Commons

0.00

(0 votes)

Original source: https://en.wikipedia.org/wiki/Boussinesq approximation (water waves). Read more

|

Retrieved from "https://handwiki.org/wiki/index.php?title=Boussinesq_approximation_(water_waves)&oldid=235678"
*[v]: View this template
*[t]: Discuss this template
*[e]: Edit this template