The model used by nsevol describes fluid motions in a two-dimensional (2D) plane: the method is two-dimensional. It thus seems strange that it is possible to add effects of a topography, since a topography is essentially a three-dimensional thing. Yet, it is possible since in a 2D-model topography can be represented by an "extra term" in the equations. This is of course only possible if the topography does not alter the fluid motions too much: the assumption of 2D motions must remain valid. Hence, the topography cannot be too high (for mountains) or too deep (for troughs).
In this picture, H(x,y) represents the fluid depth, which by default is H=1 (in dimensionless units). Suppose now that a vortex approaches the ridge: what happens to it?
If the fluid is inviscid there is conservation of potential vorticity:
w + f ------- = constant H where: w(x,y) = relative vorticity distribution f(x,y) = Coriolis parameter (background rotation) H(x,y) = fluid depthWhen the vortex climbs up the ridge, the fluid depth decreases and from this conservation law it then follows that the vortex becomes weaker: w decreases in absolute value (assume that f is constant for a moment). At the same time the "length" of the vortex decreases. Since the fluid is considered to be incompressible, this means that the vortex becomes wider as it climbs the ridge. And when the vortex descends the ridge again, its strength increases and it becomes narrower again.
Since the motions in the numerical model are two-dimensional, the squeezing and stretching of the vortex is absent, but the other effects -- weakening and widening as it climbs the ridge -- are certainly taking place. This is clearly visible in an example of a Lamb dipole crossing a ridge, which is without background rotation.
See also the cases where a Bessel monopole encounters a north-south ridge or a circular mountain.
<=== Numerical simulations of 2D vortex evolution with a Finite Difference Method.
Jos van Geffen -- Home | Site Map | Contact Me