Nonlinear waves in rotating magnetohydrodynamics

We consider an electrically conducting fluid in rotating cylindrical coordinates in which the Elsasser and magnetic Reynolds numbers are assumed to be large while the Rossby number is assumed to vanish in an appropriate limit. This may be taken as a simple model for the Earth's outer core. Fully nonlinear waves dominated by the nonlinear Lorentz forces are studied using the method of geometric optics (essentially WKB). These waves are assumed to be of the form of an asymptotic series expanded about ambient magnetic and velocity fields which vanish on the equatorial plane. They take the form of short wave, slowly varying wave trains. The first-order approximation is sinusoidal and basically the same as in the linear problem, with a dispersion relation modified by the appearance of mean terms. These mean terms, as well the undetermined amplitude functions, are found by suppressing secular terms in a “fast” variable in the second-order approximation. The interaction of the mean terms with the dispersion relation is the primary cause of behaviors which differ from the linear case. In particular, new singularities appear in the wave amplitude functions and an initial value problem results in a singularity in one of the mean terms which propagates through the fluid. The singularities corresponding to the linear ones are shown to develop when the corresponding waves propagate toward the equatorial plane.     

A reduced model for nonlinear dispersive waves in a rotating environment

Abstract

The simplest model for geophysical flows is one layer of a constant density fluid with a free surface, where the fluid motions occur on a scale in which the Coriolis force is significant. In the linear shallow water limit, there are non-dispersive Kelvin waves, localized near a boundary or near the equator, and a large family of dispersive waves. We study weakly nonlinear and finite depth corrections to these waves, and derive a reduced system of equations governing the flow. For this system we find approximate solitary Kelvin waves, both for waves traveling along a boundary and along the equator. These waves induce jets perpendicular to their direction of propagation, which may have a role in mixing. We also derive an equivalent reduced system for the evolution of perturbations to a mean geostrophic flow.     

Geostrophic vs magneto-geostrophic adjustment and nonlinear magneto-inertia-gravity waves in rotating shallow water magnetohydrodynamics

The evolution of localised jets and periodic nonlinear waves in rotating shallow water magnetohydrodynamics (rotating SWMHD) and standard rotating shallow water model (RSW) is compared within the framework of translationally-invariant 1.5-dimensional configurations, which are traditionally used in geophysical fluid dynamics for studying geostrophic adjustment and frontogenesis. Such configurations also allow for exact nonlinear wave solutions in both models. A theory of the magneto-geostrophic adjustment, i.e. adjustment of an arbitrary initial configuration to a state of magneto-geostrophic equilibrium in RSWMHD, is developed and confirmed by numerical simulations with a finite-volume well-balanced code. The code is resolving all kinds of waves in the model and corresponding weak solutions equally well. It is benchmarked by reproducing exact solutions – steady essentially nonlinear Alfvèn and mixed magneto-inertia-gravity waves – and used to demonstrate robustness of these solutions with respect to localised along-wave perturbations. It is also shown how the results on adjustment can be extended to the fully 2-dimensional case.     

Hydromagnetic waves in a thin rotating spherical shell

We consider an electrically conducting fluid confined to a thin rotating spherical shell in which the Elsasser and magnetic Reynolds numbers are assumed to be large while the Rossby number is assumed to vanish in an appropriate limit. This may be taken as a simple model for a possible stable layer at the top of the Earth's outer core. It may also be a model for the thin shells which are thought to be a source of the magnetic fields of some planets such as Mercury or Uranus. Linear hydromagnetic waves are studied using a multiple scale asymptotic scheme in which boundary layers and the associated boundary conditions determine the structure of the waves. These waves are assumed to be of the form of an asymptotic series expanded about an ambient magnetic field which vanishes on the equatorial plane and velocity and pressure fields which do not. They take the form of short wave, slowly varying wave trains. The results are compared to the author's previous work on such waves in cylindrical geometry in which the boundary conditions play no role. The approximation obtained is significantly different from that obtained in the previous work in that an essential singularity appears at the equator and nonequatorial wave regions appear.     

Hydromagnetic waves in a differentially rotating,stratified spherical shell

Abstract

Investigations of an earlier paper (Friedlander 1987a) are extended to include the effect of an azimuthal shear flow on the small amplitude oscillations of a rotating, density stratified, electrically conducting, non-dissipative fluid in the geometry of a spherical shell. The basic state mean fields are taken to be arbitrary toroidal axisymmetric functions of space that are consistent with the constraint of the ‘‘magnetic thermal wind equation''. The problem is formulated to emphasize the similarities between the magnetic and the non-magnetic internal wave problem. Energy integrals are constructed and the stabilizing/destabilizing roles of the shears in the basic state functions are examined. Effects of curvature and sphericity are studied for the eigenvalue problem. This is given by a partial differential equation (P.D.E.) of mixed type with, in general, a complex pattern of turning surfaces delineating the hyperbolic and elliptic regimes. Further mathematical complexities arise from a distribution of the magnetic analogue of critical latitudes. The magnetic extension of Laplace's tidal equations are discussed. It is observed that the magnetic analogue of planetary waves may propagate to the east and to the west.     

Nonlinear interactions between gravity waves and tides

In this study, we present the nonlinear interactions between gravity waves (GWs) and tides by using the 2D numerical model for the nonlinear propagation of GWs in the compressible atmosphere. During the propagation in the tidal background, GWs become instable in three regions, that is z = 75―85 km, z = 90―110 km and z = 115―130 km. The vertical wavelength firstly varies gradually from the initial 12 km to 27 km. Then the newly generated longer waves are gradually compressed. The longer and shorter waves occur in the regions where GWs propagate in the reverse and the same direction of the hori-zontal mean wind respectively. In addition, GWs can propagate above the main breaking region (90—110 km). During GWs propagation, not only the mean wind is accelerated, but also the amplitude of tide is amplified. Especially, after GWs become instable, this amplified effect to the tidal amplitude is much obvious.     

Nonlinear internal waves in the upper atmosphere

To investigate the mechanism of mixing in oscillatory doubly diffusive (ODD) convection, we truncate the horizontal modal expansion of the Boussinesq equations to obtain a simplified model of the process. In the astrophysically interesting case with low Prandtl number (traditionally called semiconvection), large-scale shears are generated as in ordinary thermal convection. The interplay between the shear and the oscillatory convection produces intermittent overturning of the fluid with significant mixing. By contrast, in the parameter regime appropriate to sea water, large-scale flows are not generated by the convection. However, if such flows are imposed externally, intermittent overturning with enhanced mixing is observed.     

Trapped internal waves over undular topography in a partially mixed estuary

The flow of a stratified fluid over small-scale topographic features in an estuary may generate significant internal wave activity. Lee waves and upstream influence generated at isolated topographic features have received considerable attention during the past few decades. Field surveys of a partially mixed estuary, the Rotterdam Waterway, in 1987, also showed a plethora of internal wave activity generated by isolated topography, banks and groynes. Additionally it revealed a spectacular series of resonant internal waves trapped above low-amplitude bed waves. The internal waves reached amplitudes of 3–4 m in an estuary with a mean depth of 16 m. The waves were observed during the decreasing flood tide and are thought to make a significant contribution to turbulence production and mixing. However, while stationary linear and finite amplitude theories can be used to explain the presence of these waves, it is important to further investigate their time-dependent and non-linear behaviour. With the development of advanced non-hydrostatic models it now becomes possible to further investigate these waves through numerical experimentation. This is the focus of the work presented here. The non-hydrostatic finite element numerical model FINEL3D developed by Labeur was used in the experiments presented here. The model has been shown to work well in a number of stratified flow investigations. Here, we first show that the model reproduces the field data and for idealised stationary flow scenarios that the results are in agreement with the resonant response predicted by linear theory. Then we explore the effects of non-linearity and time dependence and consider the importance of resonant internal waves for turbulence production in stratified coastal environments.Responsible Editior: Hans Burchard     

An experimental investigation of blocking by partial barriers in a rotating baroclinic annulus

We present a series of experimental investigations in which a differentially-heated annulus was used to investigate the effects of topography on rotating, stratified flows with similarities to the Earth’s atmospheric or oceanic circulation. In particular, we compare and investigate blocking effects via partial mechanical barriers to previous experiments by the authors utilising azimuthally-periodic topography. The mechanical obstacle used was an isolated ridge, forming a partial barrier, employed to study the difference between partially blocked and fully unblocked flow. The topography was found to lead to the formation of bottom-trapped waves, as well as impacting the circulation at a level much higher than the top of the ridge. This produced a unique flow structure when the drifting flow and the topography interacted in the form of an “interference” regime at low Taylor number, but forming an erratic “irregular” regime at higher Taylor number. The results also showed evidence of resonant wave-triads, similar to those noted with periodic wavenumber-3 topography by Marshall and Read (Geophys. Astrophys. Fluid Dyn., 2015, 109), though the component wavenumbers of the wave-triads and their impact on the flow were found to depend on the topography in question. With periodic topography, wave-triads were found to occur between both the baroclinic and barotropic components of the zonal wavenumber-3 mode and the wavenumber-6 baroclinic component, whereas with the partial barrier two nonlinear resonant wave-triads were noted, each sharing a common wavenumber-1 mode.     

Nonlinear wave propagation on a sphere: Interaction between Rossby waves and gravity waves; stability of the Rossby waves

Abstract

An analytical spectral model of the barotropic divergent equations on a sphere is developed using the potential-stream function formulation and the normal modes as basic functions. Explicit expressions of the coefficients of nonlinear interaction are obtained in the asymptotic case of a slowly rotating sphere, i.e. when the normal modes can be expressed as single spherical harmonics.     

Hydromagnetic waves in rapidly rotating spherical shells generated by poloidal decay modes

Abstract

This paper presents the first attempt to examine the stability of a poloidal magnetic field in a rapidly rotating spherical shell of electrically conducting fluid. We find that a steady axisymmetric poloidal magnetic field loses its stability to a non-axisymmetric perturbation when the Elsasser number A based on the maximum strength of the field exceeds a value about 20. Comparing this with observed fields, we find that, for any reasonable estimates of the appropriate parameters in planetary interiors, our theory predicts that all planetary poloidal fields are stable, with the possible exception of Jupiter. The present study therefore provides strong support for the physical relevance of magnetic stability analysis to planetary dynamos. We find that the fluid motions driven by magnetic instabilities are characterized by a nearly two-dimensional columnar structure attempting to satisfy the Proudman-Taylor theorm. This suggests that the most rapidly growing perturbation arranges itself in such a way that the geostrophic condition is satisfied to leading order. A particularly interesting feature is that, for the most unstable mode, contours of the non-axisymmetric azimuthal flow are closely aligned with the basic axisymmetric poloidal magnetic field lines. As a result, the amplitude of the azimuthal component of the instability is smaller than or comparable with that of the poloidal component, in contrast with the instabilities generated by toroidal decay modes (Zhang and Fearn, 1994). It is shown, by examining the same system with and without fluid inertia, that fluid inertia plays a secondary role when the magnetic Taylor number T_m ? 10⁵. We find that the direction of propagation of hydromagnetic waves driven by the instability is influenced strongly by the size of the inner core.     

完整Coriolis力作用下非线性Rossby波的精确解

从包含完整Coriolis力的Boussinesq近似的斜压大气运动方程组出发,利用半地转近似导出β效应和地球旋转水平分量fH=2Ωcosφ共同作用下的大气非线性Rossby波动所满足的KdV方程,求得了椭圆余弦波解和孤立波解.结果分析表明,若扰动与纬度有关,Coriolis参数分量fH将影响波动传播的频率特征,并加强水平散度对斜压Rossby波的作用;如果扰动与纬度无关,则 Coriolis 参数分量fH的影响消失.     

Transverse baroclinic oscillations in Lake Überlingen

An episode of velocity measurements in the epilimnion at a midlake station in Lake Überlingen and taken from the campaign in October 1972 discloses a uni-nodal Poincaré-type baroclinic mode response with a 4 h period. We discuss the data and interpret it in terms of a two-layered linear wave model on the rotating Earth.

---

     

Experimental study of long nonlinear internal waves in rotating fluid

Experiments were performed on the rotating platform 14 m in diameter equipped with a simple internal wave generator. Internal waves were generated for a wide range of Coriolis parameters. When the rotation is very weak, i.e., when the internal Rossby radius of deformation is much larger than the wavelength, then the stable nonlinear waves generated are solitary waves. These have a horizontal crest, as in the nonrotating case. When the rotation is strong, i.e., when the internal Rossby radius is at most comparable with the wavelength, then Sverdrup-like periodic waves can be generated, but no solitary wave can then propagate. For the intermediate case, Ostrovsky waves are generated. Their phase speed increases with increasing amplitude. Then, there are two characteristic wave lengths: one which varies with the inverse square root of the amplitude, as for the KdV wave, and the other, linked with the rotation, which varies as the square root of the amplitude. The experimental results are thus in agreement with most of the conclusions in recent analytical developments.     

Nonlinear stability of baroclinic flows over topography

Abstract

A new nonlinear stability criterion is derived for baroclinic flows over topography in spherical geometry. The stability of a wide class of exact three-dimensional nonlinear steady state solutions subject to arbitrary disturbances is established. The resonance condition, at the highest total wavenumber, for the steady state solutions and the stability criteria for baroclinic flow in the absence of topography provide the boundaries of the regions of stability in the presence of topography. The analogous results for flow on periodic or infinite beta planes incorporating non-orthogonal function large scale flows are also discussed.     

Hydromagnetic waves in a differentially rotating annulus IV. Insulating boundaries

Abstract

The first three papers in this series (Fearn, 1983b, 1984, 1985) have investigated the stability of a strong toroidal magnetic field B_o =B_o(s?)Φ [where (s?. Φ, z?) are cylindrical polars] in a rapidly rotating system. The application is to the cores of the Earth and the planets but a simpler cylindrical geometry was chosen to permit a detailed study of the instabilities present. A further simplification was the use of electrically perfectly conducting boundary conditions. Here, we replace these with the boundary conditions appropriate to an insulating container. As expected, we find the same instabilities as for a perfectly conducting container, with quantitative changes in the critical parameters but no qualitative differences except for some interesting mixing between the ideal (“field gradient”) and resistive modes for azimuthal wavenumber m=1. In addition to these modes, we have also found the “exceptional” slow mode of Roberts and Loper (1979) and we investigate the conditions required for its instability for a variety of fields B_o(s?) Roberts and Loper's analysis was restricted to the case B_o∝s? and they found instability only for m=1 and ?1 <ω<0 [where ω is the frequency non-dimensionalised on the slow timescale τ_x, see (1.5)]. For other fields we found the necessary conditions to be less “exceptional”. One surprising feature of this instability is the importance of inertia for its existence. We show that viscosity is an alternative destabilising agent. The standard (magnetostrophic) approximation of neglecting inertial (and viscous) terms in the equation of motion has the effect of filtering out this instability. The field strength required for this “exceptional” mode to become unstable is found to be very much larger than that thought to be present in the Earth's core, so we conclude that this mode is unlikely to play an important role in the dynamics of the core.     

Seismic response of multi-layered basins with velocity gradients upon incidence of plane shear waves

A boundary integral scheme based on boundary-integral discrete-wave-number approach has been developed to compute the seismic response of two-dimensional irregular-shaped basins with horizontal soil layers. Each layer exhibits a linear gradient of shear wave velocity with depth. The approach combines the boundary-integral representation of the seismic wave field outside the basin with plane wave representation of the seismic wave field inside the basin. The propagation throughout the layers is performed by matrix propagators in which the effect of the vertical variation of the velocity is incorporated by using confluent hyper-geometric functions of the Whittaker type. Our method is tested against otherwell-accepted solutions for the case of a circular basin with excellent agreement. Test of the ground response for a semi-circular basin with radius a shows that stable solutions are obtained if the chosen model parameters satisfy following conditions: (1) the distance from the sources to the interface is greater than 0·1a; (2) the distance between the sources is smaller than a quarter of the incident wavelength; and (3) the discrete wave-number step is smaller than 2π/4a. The computation of ground response of basins with a sharp interface and several horizontal deposits leads to the following main results: (1) the amplification of a basin with velocity gradients is larger than that of a basin with homogeneous layers; (2) the frequencies of the second- and third-order harmonics for a basin with velocity gradients are lower than those of a basin with homogeneous layers; and (3) the response amplitude of the basin with velocity gradients attenuates more slowly in time domain than when layers are homogeneous. Since these results have been obtained for realistic values of basin geometrical and mechanical consideration, they should find some interest in earthquake engineering or seismic microzonation studies. © 1998 John Wiley & Sons, Ltd.