A single density contour is drawn over the topography h(x). The isopycnal displacement η is defined as the distance the density contour is displaced from its far upstream value. For waves of depression, η is negative and for waves of elevation it is positive.
A comparison of computation times from the embedded and mapped methods for several grid sizes (N grid ) and waves. Results for a small wave with U 0 = 1.25c j (solid line) and a large wave with U 0 = 1.11c j (dashed line). The embedded method is indicated with circles (o) and the mapped method with stars (*). The results are normalized by longest computation time. N grid = m × n where m and n are the number of grid points in x and z, respectively.
Plots of the density field for trapped depression waves (top) and elevation waves (bottom). On the left (a) and (c), the pycnocline is centered at z 0 = 0.75H and on the right (b) and (d), it is centered at z 0 = 0.25H. In these cases, U 0 = 1.06c j . The thick black line represents the bottom topography.
Diagram of wave properties as a function of the background velocity U 0. (a) Maximum value of |η| (scaled by 1/H) for waves over hole topography with z 0 = 0.6H and hill topography with z 0 = 0.4H. (b) A measure of the wave width for the corresponding cases in panel (a). The wave width (scaled by 1/H) is measured as twice the distance between the wave center and location where the wave induced surface velocities reach half of their extreme value. In both cases the conjugate flow speed is c j = 1.02c lw .
(a) Extreme value of |η| (scaled by 1/H) with varying background speeds U 0 for comparison between non-Boussinesq (solid) and Boussinesq (dashed) results with Δρ = 0.05, z 0 = 0.75H and depression topography. (b) A measurement of wave width (scaled by 1/H) is also included. The background speed is scaled by the conjugate flow speed, c j , as calculated for the non-Boussinesq case. The wave width is measured as in Fig. 4 .
An example of the spatial structure of waves from the non-Boussinesq (solid) and Boussinesq (dashed) solutions in Fig. 5 for U 0 = 1.23c j . (a) Normalized wave-induced velocities at the surface. (b) Wave-induced velocities through x = 0 scaled by the conjugate flow speed.
Comparison of non-Boussinesq (solid) and Boussinesq (dashed) maximum |η| (scaled by 1/H) for several background speeds and stratifications. (a) Δρ = 0.05, (b) Δρ = 0.1, (c) Δρ = 0.15, and (d) Δρ = 0.2.
Density fields after 5 (10) h on the left (right) column, where U 0 = [0.96, 0.91, 0.85, 0.80]c j from top to bottom.
Horizontal velocity fields with three white isopycnals superimposed for the double hill case. The background speed is U 0 = 0.96c j . The results are presented at non-dimensional times (a) t* = 745, (b) t* = 993.3, (c) t* = 1987, and (d) t* = 2483. The scaling time is taken to be t d = H/c j .
(a) An example of an asymmetric solution of the DJL equation, with U 0 = 0.96c j . (b) Wave characteristics for asymmetric DJL solutions over several background speeds. The absolute value of the displacement of the ρ = 1 contour (scaled by 1/H) at x/H = 20 (*) and at x/H = 0 (o) are shown.
Comparison of wave amplitude (left) and wave width (right) for two density stratifications with z 0 = 0.75H; Δρ = 0.05 in (a) and (b) and Δρ = 0.2 in (c) and (d). Three cases are shown: non-Boussinesq solutions with N NB (z) (solid), for the Boussinesq solutions with N NB (z) (dashed), and for the Boussinesq solutions with N B (z). The velocity is scaled by the non-Boussinesq conjugate flow speed in each case.
Local conjugate flow speeds using water column depth above the center of the topography H loc . Upstream of the topography the conjugate flow speed is c j = 1.17c lw .
Article metrics loading...
Full text loading...