This page describes a numerical simulation of gravity waves due to wintertime
winds blowing over the Scandinavian Mountain ranges in Norway and Sweden.
Computational Domain
The layout of the of the simulation is shown in the figure below. Note that
the boundaries of the computation domain (shown by the thin white line) has
been rotated by 26 degrees with respect to lines of constant latitude.
Just the computational domain is shown in the figure below. In this plot,
the surface has been interpolated to rectangular coordinates.
The mesh is uniform at 500 m spacing from the Norwegian coastline to a
point 50 km east-southeast of Esrange, Sweden. Gentle stretching is used
upstream and downstream of these points. The mesh is also uniformly
spaced at 500 m in the lateral direction over the central regon +/- 50 km
N-NE and S-SW of Esrange. Gentle stretching is once again used from this
region to the outer boundaries. The mesh is uniformly spaced at 250 m for
all altitudes in the z direction. In total the mesh contains
650 x 274 x 576 points and sponge layers are used on all external boundaries
in order to absorb outgoing waves.
Wind and Thermodynamic Profiles
The wind and temperature profiles were based on radiosonde mesurements
taken over Nome Alaska on Jan 5 2017 at 12:00 UTC. The radiosonde
measurements were smoothed somewhat and the data was extended between
30 and 113 km using climatology. Polynomial extrapolation was used above
113 km.
Initial Wind Condition
Here the winds are initially set to zero from the surface to an altitude of
about 20 km. At this point the wind transitions to the profiles shown in the
plots above. Forcing terms then gradually introduce winds near the surface
with the
objective of achieving the profile shown above within a two hour period.
A hyperbolic tangent function is used in order to produce gentle
acceleration of the wind near the beginning and end of the forcing period.
The maximum forcing rate is equivalent to that of a linear ramp with a
duration of 30 minutes.
Movies
Shortened animation of u' in the xz plane at the position y = 0 km
Shortened animation of vorticity magnitude in the xz plane at the position y = 0 km
Animation of u' in the xz plane at the position y = 0 km
Animation of w' in the xz plane at the position y = 0 km
Animation of potential temperature perturbation in the xz plane at the position y = 0 km
Animation of vorticity magnitude in the xz plane at the position y = 0 km
Animation of w' in the xy plane at an altitlde of 35 km
Animation of vorticity magnitude in the xy plane at an altitlde of 35 km
overlayed on the computational domain
Animation of w' in the xy plane at an altitlde of 40 km
Animation of vorticity magnitude in the xy plane at an altitlde of 40 km
overlayed on the computational domain
Animation of w' in the xy plane at an altitlde of 45 km
Animation of vorticity magnitude in the xy plane at an altitlde of 45 km
overlayed on the computational domain
Animation of w' in the xy plane at an altitlde of 60 km
Animation of vorticity magnitude in the xy plane at an altitlde of 60 km
overlayed on the computational domain
Animation of w' in the xy plane at an altitlde of 90 km
Animation of vorticity magnitude in the xy plane at an altitlde of 90 km
overlayed on the computational domain
Animation of temperature fluxuation in the xy plane at an altitlde of 90 km