Internal waves in a rotating stratified fluid in an arbitrary gravitational field

Abstract

Small amplitude oscillations of a uniformly rotating, density stratified, Boussinesq, non-dissipative fluid are examined. A mathematical model is constructed to describe timedependent motions which are small deviations from an initial state that is motionless with respect to the rotating frame of reference. The basic stable density distribution is allowed to be an arbitrary prescribed function of the gravitational potential. The problem is considered for a wide class of gravitational fields. General properties of the eigenvalues and eigenfunctions of square integrable oscillations are demonstrated, and a bound is obtained for the magnitude of the frequencies. The modal solutions are classified as to type. The eigenfunctions for the pressure field are shown to satisfy a second-order partial differential equation of mixed type, and the equation is obtained for the critical surfaces which delineate the elliptic and hyperbolic regions. The nature of the problem is examined in detail for certain specific gravitational fields, e.g., a radially symmetric field. Where appropriate, results are compared with those of other investigations of waves in a rotating fluid of spherical configuration and the novel aspects of the present treatment are emphasized. Explicit modal solutions are obtained in the specific example of a fluid contained in a rigid cylinder, stratified in the presence of vertical gravity, with the buoyancy frequency N being an arbitrary prescribed function of depth.     

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.     

Large-scale convective dynamos in a stratified rotating plane layer

We study the effect of stratification on large-scale dynamo action in convecting fluids in the presence of background rotation. The fluid is confined between two horizontal planes and both boundaries are impermeable, stress-free and perfectly conducting. An asymptotic analysis is performed in the limit of rapid rotation (τ???1 where τ is the Taylor number). We analyse asymptotic magnetic dynamo solutions in rapidly rotating systems generalising the results of Soward [A convection-driven dynamo I. The weak field case. Philos. Trans. R. Soc. Lond. A 1974, 275, 611–651] to include the effects of compressibility. We find that in general the presence of stratification delays the efficiency of large-scale dynamo action in this regime, leading to a reduction of the onset of dynamo action and in the nonlinear regime a diminution of the large-scale magnetic energy for flows with the same kinetic energy.     

Penetrative convection in rotating fluids: A model for the base of stellar convection zones

Abstract

To model penetrative convection at the base of a stellar convection zone we consider two plane parallel, co-rotating Boussinesq layers coupled at their fluid interface. The system is such that the upper layer is unstable to convection while the lower is stable. Following the method of Kondo and Unno (1982, 1983) we calculate critical Rayleigh numbers R_c for a wide class of parameters. Here, R_c is typically much less than in the case of a single layer, although the scaling R_c～T^2/3 as T → ∞ still holds, where T is the usual Taylor number. With parameters relevant to the Sun the helicity profile is discontinuous at the interface, and dominated by a large peak in a thin boundary layer beneath the convecting region. In reality the distribution is continuous, but the sharp transition associated with a rapid decline in the effective viscosity in the overshoot region is approximated by a discontinuity here. This source of helicity and its relation to an alpha effect in a mean-field dynamo is especially relevant since it is a generally held view that the overshoot region is the location of magnetic field generation in the Sun.     

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.     

Turbulent mixing in a rotating,stratified fluid

Abstract

An experimental study was carried out to investigate the effect of rotation on turbulent mixing in a stratified fluid when the turbulence in the mixed layer is generated by an oscillating grid. Two types of experiments were carried out: one of them is concerned with the deepening of the upper mixed layer in a stable, two-fluid system, and the other deals with the interaction between a stabilizing buoyancy flux and turbulence.

In the first type of experiments, it was found that rotation suppresses entrainment at larger Rossby numbers. As the Rossby number becomes smaller (Ro 0.1), the entrainment rate increases with rotation—the onset of this phenomenon, however, was found to coincide with the appearance of coherent vortices within the mixed layer. The radiation of energy from the mixed layer to the lower non-turbulent layer was found to occur and the magnitude of the energy flux was found to be increased with the rotational frequency. It was also observed that vortices are generated, rather abruptly, in the lower layer as the mixed layer deepens.

In the second set of experiments a quasi-steady mixed layer was found to develop of which the thickness varies with rotation in a fashion that is consistent with the result of the first experiment. Also the rotation was found to delay the formation of a pycnocline.     

Forces on spheres moving horizontally in a rotating stratified fluid

Abstract

Measurements have been made of the net horizontal force F acting on a sphere moving with horizontal velocity U (Reynolds numbers in the range 10²-10⁴) through a stratified fluid rotating about a vertical axis with uniform angular velocity Ω. In both homogeneous and stratified rotating fluids with small Rossby number R(R = U/Ωa ? 1 where a is the radius of the sphere) the force F is of magnitude 2ΩρUV (where ρ is the density of the fluid and V is the volume of the sphere). In a homogeneous fluid the relative directions of F and U were found to depend on the quantity F = 8Ωa ²/UD (where D is the depth of the fluid in which the object is placed (Mason, 1975)). In a rotating stratified fluid the relative directions of F and U are found to depend on the inverse Froude number k(k = Na/U where N ² = (g/δ)?ρ/?z) provided D > 4aΩ/N. In a homogeneous fluid with F ? 1 the force F is mainly in the U direction (a drag force due to inertial wave radiation) and is ～ ?0.4 |MX 2ΩρUV For F ? 1 a “Taylor column” occurs and the force, in correspondence with theoretical expectations, is ～ - 2Ω |MX UρV In a rotating stratified fluid with N ～2Ω and k ? 1 the force F is mainly in the U direction but is roughly one half of that occurring in the homogeneous situation with F ? 1 (tentatively explained as due to the evanescence of inertia-gravity disturbances). In a rotating stratified fluid with k ? 1 the flow should have no vertical motion (as with F ? 1) and again in correspondence with theoretical expectations the drag is ～ ?2 Ω |MX UρV. In a non-rotating stratified fluid the drag coefficient C _D(C _D = F U/½?ρU ²) was measured in the range k = 0.1 to 10 and had a maximum value ～ 1.2 for k ～ 3.     

The incidence of internal waves onto a thin submerged barrier

Abstract

The problem of oblique incidence of internal ocean waves on a thin submerged ocean barrier is considered when the ocean has exponential density stratification. A Wiener-Hopf approach is used combined with numerical evaluation of series. Results for the reflected energy are obtained and reveal a complex dependence on incidence and barrier height. Application of this model to waves incident on the Mid-Atlantic ridge suggests that the ridge almosts isolates first mode energy on one side of the ocean from the other side. In certain circumstances there, is a surprising appearance of “barrier” waves. These waves are closely confined to the barrier and propagate along it.     

Transitions to chaotic thermal convection in a rapidly rotating spherical fluid shell

Abstract

Numerical simulations of thermal convection in a rapidly rotating spherical fluid shell heated from below and within have been carried out with a nonlinear, three-dimensional, time-dependent pseudospectral code. The investigated phenomena include the sequence of transitions to chaos and the differential mean zonal rotation. At the fixed Taylor number T _a =10⁶ and Prandtl number Pr=1 and with increasing Rayleigh number R, convection undergoes a series of bifurcations from onset of steadily propagating motions SP at R=R _c = 13050, to a periodic state P, and thence to a quasi-periodic state QP and a non-periodic or chaotic state NP. Examples of SP, P, QP, and NP solutions are obtained at R = 1.3R _c , R = 1.7 R _c , R = 2R _c , and R = 5 R _c , respectively. In the SP state, convection rolls propagate at a constant longitudinal phase velocity that is slower than that obtained from the linear calculation at the onset of instability. The P state, characterized by a single frequency and its harmonics, has a two-layer cellular structure in radius. Convection rolls near the upper and lower surfaces of the spherical shell both propagate in a prograde sense with respect to the rotation of the reference frame. The outer convection rolls propagate faster than those near the inner shell. The physical mechanism responsible for the time-periodic oscillations is the differential shear of the convection cells due to the mean zonal flow. Meridional transport of zonal momentum by the convection cells in turn supports the mean zonal differential rotation. In the QP state, the longitudinal wave number m of the convection pattern oscillates among m = 3,4,5, and 6; the convection pattern near the outer shell has larger m than that near the inner shell. Radial motions are very weak in the polar regions. The convection pattern also shifts in m for the NP state at R = 5R _c , whose power spectrum is characterized by broadened peaks and broadband background noise. The convection pattern near the outer shell propagates prograde, while the pattern near the inner shell propagates retrograde with respect to the basic rotation. Convection cells exist in polar regions. There is a large variation in the vigor of individual convection cells. An example of a more vigorously convecting chaotic state is obtained at R = 50R _c . At this Rayleigh number some of the convection rolls have axes perpendicular to the axis of the basic rotation, indicating a partial relaxation of the rotational constraint. There are strong convective motions in the polar regions. The longitudinally averaged mean zonal flow has an equatorial superrotation and a high latitude subrotation for all cases except R = 50R _c , at this highest Rayleigh number, the mean zonal flow pattern is completely reversed, opposite to the solar differential rotation pattern.     

Ringdown of inertial waves in a spheroidal shell of rotating fluid

The role of inertial waves in the dynamies of the Earth's fluid core has been investigated through laboratory experiments on a spheroidal shell of rotating fluid. In these experiments inertial waves of azimuthal wavenumber 1, Ekman number 0(10^?5), Rossby number 0(10^?1) were excited by precession of an inner spheroidal body. Proximity to resonance was achieved by adjusting the ratio of the frequency of precession of the inner body to the rotational speed of the container to be near the eigenfrequency of the inertial wave mode being studied. Once the system was near resonance the perturbation was stopped and ringdown records were obtained. Amplitude, eigenfrequency and decay rate were recovered simultaneously for the principal and neighbouring modes excited using an iterative linearized least squares procedure.The recovery of complex eigenfrequencies for non-axisymmetric inertial waves in this shell geometry has given experimental verification of their existence. For those waves of azimuthal wavenumber one, a significant nonlinear interaction among modes is inferred from the simultaneous recovery of neighbouring modes. Other non linear effects include a mean azimuthal flow which appears to be stable for the low spatial order modes studied. These results contrast with highly unstable mean flow found experimentally in similar experiments carried out in cylindrical geometry.     

Scattering of internal waves in a two-layer fluid flowing through a channel with small undulations

Using two-dimensional linear water wave theory, we consider the problem of normal water wave (internal wave) propagation over small undulations in a channel flow consisting of a two-layer fluid in which the upper layer is bounded by a fixed wall, an approximation to the free surface, and the lower one is bounded by a bottom surface that has small undulations. The effects of surface tension at the surface of separation is neglected. Assuming irrotational motion, a perturbation analysis is employed to calculate the first-order corrections to the velocity potentials in the two-layer fluid by using Green’s integral theorem in a suitable manner and the reflection and transmission coefficients in terms of integrals involving the shape function c(x) representing the bottom undulation. Two special forms of the shape function are considered for which explicit expressions for reflection and transmission coefficients are evaluated. For the specific case of a patch of sinusoidal ripples having the same wave number throughout, the reflection coefficient up to the first order is an oscillatory function in the quotient of twice the interface wave number and the ripple wave number. When this quotient approaches one, the theory predicts a resonant interaction between the bed and the interface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of the incident wave energy occurs if this number is large. Again, when a patch of sinusoidal ripples having two different wave numbers for two consecutive stretches is considered, the interaction between the bed and the interface near resonance attains in the neighborhood of two (singular) points along the x-axis (when the ripple wave number of the bottom undulation become approximately twice as large as the interface wave number). The theoretical observations are presented in graphical form.     

The existence of internal solitary waves in a two-fluid system near the KdV limit

Abstract

Two fluid layers of constant density lying one over the other on top of a rigid horizontal lower boundary with either a free upper surface or a rigid upper boundary can support solitary waves. The existence of a unique branch of such waves emanating from the horizontal flow at a critical speed U _? is demonstrated in both cases by use of the Nash—Moser implicit function theorem. These results complement the global results of Amick and Turner (1986) and are analogous to the work of Friedrichs and Hyers (1954) and Beale (1977) for surface waves. It is also noted that the most obvious variational principle which characterizes these waves as constrained extremals (Benjamin, 1984) is of indefinite type, having a Hessian with infinitely many positive and infinitely many negative eigenvalues.     

Programme of the US.-Swiss conference on internal waves in geophysical contexts

Abstract

Evidence from radio polarization measurements is reviewed that indicates that most galactic magnetic field structures fall into one of two categories: axisymmetric spiral and bisymmetric spiral. The resultant challenges to dynamo theorists is stated. Estimates of the magnetic field strengths based on equipartition of field and cosmic ray energies are given, but deviations from equipartition are inferred. Possible goals for future research are suggested.     

Soret and Dufour effects on thermohaline convection in rotating fluids

Using linear and weakly nonlinear stability theory, the effects of Soret and Dufour parameters are investigated on thermohaline convection in a horizontal layer of rotating fluid, specifically the ocean. Thermohaline circulation is important in mixing processes and contributes to heat and mass transports and hence the earth’s climate. A general conception is that due to the smallness of the Soret and Dufour parameters their effect is negligible. However, it is shown here that the Soret parameter, salinity and rotation stabilise the system, whereas temperature destabilises it and the Dufour parameter has minimal effect on stationary convection. For oscillatory convection, the analysis is difficult as it shows that the Rayleigh number depends on six parameters, the Soret and Dufour parameters, the salinity Rayleigh number, the Lewis number, the Prandtl number, and the Taylor number. We demonstrate the interplay between these parameters and their effects on oscillatory convection in a graphical manner. Furthermore, we find that the Soret parameter enhances oscillatory convection whereas the Dufour parameter, salinity Rayleigh number, the Lewis number, and rotation delay instability. We believe that these results have not been elucidated in this way before for large-scale fluids. Furthermore, we investigate weakly nonlinear stability and the effect of cross diffusive terms on heat and mass transports. We show the existence of new solution bifurcations not previously identified in literature.     

Finite-amplitude mountain waves in a compressible atmosphere including the effects of dissipation

Abstract

Adiabatic, two-dimensional, steady-state finite-amplitude, hydrostatic gravity waves produced by flow over a ridge are considered. Nonlinear self advection steepens the wave until the streamlines attain a vertical slope at a critical height z_c. The height z_c , where this occurs, depends on the ridge crest height and adiabatic expansion of the atmosphere. Dissipation is introduced in order to balance nonlinear self advection, and to maintain a marginal state above z_c. The approach is to assume that the wave is inviscid except in a thin layer, small compared to a vertical wavelength, where dissipation cannot be neglected. The solutions in each region are matched to obtain a continuous solution for the streamline displacement δ. Solutions are presented for different values of the nondimensional dissipation parameter β. Eddy viscosity coefficients and the thickness of the dissipative layer are expressed as functions of β, and their magnitudes are compared to other theoretical evaluations and to values inferred from radar measurements of the stratosphere.

The Fourier spectrum of the solution for z ≫ z_c is shown to decay exponentially at large vertical wave numbers n. In comparison, a spectral decay law n ^?-8/3 characterizes the marginal state of the wave at z = z_c .     

Effects of a sill on selective withdrawal through a point sink in a linearly stratified fluid

Reservoir topography has significant effects on mechanism of selective withdrawal. We study selective withdrawal of a linearly stratified fluid through a point sink in a long rectangular reservoir with a sill at the bed experimentally and analytically. Experiments were conducted for various flow rates and sill heights to evaluate their effects on withdrawal layer in the inertial-buoyancy regime. Transition to the steady state is discussed in terms of dynamics of shear waves over the sill. The results show that in presence of a sill the withdrawal layer is thicker and higher in elevation. For flows controlled at the sill crest, the withdrawal layer thickness in the point-sink flow is almost equal to that in the similar line-sink flow. We propose an analytical relation for the withdrawal layer thickness in presence of a sill and confirm it by experiments.     

On the onset of thermal convection in slowly rotating fluid shells

Abstract

The problem of the removal of the degeneracy of the patterns of convective motion in a spherically symmetric fluid shell by the effects of rotation is considered. It is shown that the axisymmetric solution is preferred in sufficiently thick shells where the minimum Rayleigh number corresponds to degree l = 1 of the spherical harmonics. In all cases with l > 1 the solution described by sectional spherical harmonics Y^l _l (θ,φ) is preferred.     

A quasi-geostrophic numerical model of a rotating internally heated fluid

Abstract

A quasi-geostrophic numerical model of flow in a rotating channel is integrated under conditions typical of laboratory experiments with an internally heated annulus system. Compared to a laboratory experiment, or a full Navier-Stokes simulation, the quasi geostrophic numerical model is a simple system. It includes nonlinear interactions, dissipation via conventional parameterizations of Ekman layers and internal diffusion, and a steady forcing term which represents heating near the centre of the channel and cooling near both sides. Explicit boundary layers, cylindrical geometry effects, horizontal variations in static stability and variations in conductivity and diffusivity with temperature are all absent, and ageostrophic advection is incompletely represented. Nevertheless, over a range of parameters, flows are produced which strongly resemble those seen in the laboratory thus suggesting that the most important physical processes are represented. The numerical model is used to map out a regime diagram which includes examples of steady flows, flows with periodic time dependence (wavenumber vacillations) and flows which are irregularly time dependent.     

The effect of stratification and compressibility on anelastic convection in a rotating plane layer

We study the effect of stratification and compressibility on the threshold of convection and the heat transfer by developed convection in the nonlinear regime in the presence of strong background rotation. We consider fluids both with constant thermal conductivity and constant thermal diffusivity. The fluid is confined between two horizontal planes with both boundaries being impermeable and stress-free. An asymptotic analysis is performed in the limits of weak compressibility of the medium and rapid rotation (τ^?1/12???|θ|???1, where τ is the Taylor number and θ is the dimensionless temperature jump across the fluid layer). We find that the properties of compressible convection differ significantly in the two cases considered. Analytically, the correction to the characteristic Rayleigh number resulting from small compressibility of the medium is positive in the case of constant thermal conductivity of the fluid and negative for constant thermal diffusivity. These results are compared with numerical solutions for arbitrary stratification. Furthermore, by generalizing the nonlinear theory of Julien and Knobloch [Fully nonlinear three-dimensional convection in a rapidly rotating layer. Phys. Fluids 1999, 11, 1469–1483] to include the effects of compressibility, we study the Nusselt number in both cases. In the weakly nonlinear regime we report an increase of efficiency of the heat transfer with the compressibility for fluids with constant thermal diffusivity, whereas if the conductivity is constant, the heat transfer by a compressible medium is more efficient than in the Boussinesq case only if the specific heat ratio γ is larger than two.