The theory for internal waves differs in the description of interfacial waves and vertically propagating internal waves. These are treated separately below.
Interfacial waves In the simplest case, one considers a two-layer fluid in which a slab of fluid with uniform density \rho_1 overlies a slab of fluid with uniform density \rho_2. Arbitrarily the interface between the two layers is taken to be situated at z=0. The fluid in the upper and lower layers are assumed to be
irrotational. So the velocity in each layer is given by the gradient of a
velocity potential, {\vec{u}=\nabla\phi,} and the potential itself satisfies
Laplace's equation: :\nabla^2\phi=0. Assuming the domain is unbounded and two-dimensional (in the x-z plane), and assuming the wave is
periodic in x with
wavenumber k>0, the equations in each layer reduces to a second-order ordinary differential equation in z. Insisting on bounded solutions the velocity potential in each layer is :\phi_1(x,z,t) = A e^{-kz} \cos(kx - \omega t) and :\phi_2(x,z,t) = A e^{kz} \cos(kx - \omega t), with A the
amplitude of the wave and \omega its
angular frequency. In deriving this structure, matching conditions have been used at the interface requiring continuity of mass and pressure. These conditions also give the
dispersion relation: :\omega^2 = g^\prime k in which the reduced gravity g^\prime is based on the density difference between the upper and lower layers: :g^\prime = \frac{\rho_2-\rho_1}{\rho_2+\rho_1}\, g, with g the
Earth's gravity. Note that the dispersion relation is the same as that for deep water
surface waves by setting g^\prime=g.
Internal waves in uniformly stratified fluid The structure and dispersion relation of internal waves in a uniformly stratified fluid is found through the solution of the linearized conservation of mass, momentum, and internal energy equations assuming the fluid is incompressible and the background density varies by a small amount (the
Boussinesq approximation). Assuming the waves are two dimensional in the x-z plane, the respective equations are :\partial_x u + \partial_z w = 0 :\rho_{00} \partial_t u = - \partial_x p :\rho_{00} \partial_t w = - \partial_z p - \rho g :\partial_t \rho = -w d\rho_0/dz in which \rho is the perturbation density, p is the pressure, and (u,w) is the velocity. The ambient density changes linearly with height as given by \rho_0(z) and \rho_{00}, a constant, is the characteristic ambient density. Solving the four equations in four unknowns for a wave of the form \exp[i(kx+mz-\omega t)] gives the dispersion relation :\omega^2 = N^2 \frac{k^2}{k^2+m^2} = N^2 \cos^2\Theta in which N is the
buoyancy frequency and \Theta=\tan^{-1}(m/k) is the angle of the wavenumber vector to the horizontal, which is also the angle formed by lines of constant phase to the vertical. The
phase velocity and
group velocity found from the dispersion relation predict the unusual property that they are perpendicular and that the vertical components of the phase and group velocities have opposite sign: if a wavepacket moves upward to the right, the crests move downward to the right. ==Internal waves in the ocean==