 Solution

Write the density as a hydrostatic reference profile $\hat{\rho }(z)$ plus a small perturbation ${\rho }^{\prime }$, and define buoyancy and buoyancy frequency byThe Boussinesq approximation to the Navier-Stokes equation, together with mass conservation, iswhere $D/Dt$ is the material derivative, $\nu$ is kinematic viscosity, and $\kappa$ is mass diffusivity. Linearizing about rest and eliminating pressure and buoyancy givesFor a plane wave proportional to $e^{i(kx+mz-\omega t)}$, with $K^2=k^2+m^2$, the viscous-diffusive dispersion relation isIn the inviscid limit this becomesFor weak diffusion the two oscillatory roots areso a freely evolving Fourier mode decays.

A single plane wave is also an exact solution of the nonlinear equations. Every field depends only on its phase $\Theta =\mathbf{k}\cdot \mathbf{x}-\omega t$, while incompressible flow gives $\mathbf{u}\cdot \mathbf{k}=0$. Therefore $\mathbf{u}\cdot \nabla$ annihilates both $\mathbf{u}(\Theta )$ and $b(\Theta )$, and all nonlinear advection terms vanish.

For a boundary-forced wave with real $\omega$, nonzero $\nu$ or $\kappa$ instead makes the bulk vertical wavenumber complex, attenuating the propagating beam. Because diffusion raises the spatial order of the equations, additional short vertical-wavenumber roots form viscous and scalar boundary layers; they allow a no-slip velocity condition and a scalar no-flux condition to accompany impermeability. These layers and bulk attenuation become essential near critical internal-wave reflection, where the inviscid reflected wavelength collapses.