On the spin-up of a stably stratified,electrically conducting fluid. Part 1: Spin-up controlled by the hartmann layer

Abstract

The linear spin-up of a stably stratified, electrically conducting fluid within an electrically insulating cylindrical container in the presence of an applied axial magnetic field is analyzed for those cases in which electric currents generated within the steady Hartmann boundary layer control the fluid interior. It is shown how to obtain the known spin-up times for a homogeneous, nonconducting fluid (τ = E ^-½), a stably stratified, nonconducting fluid (τ = (σS/E, E ^?1) and a homogeneous conducting fluid (τ = α^?1 E ^-½) from the present formulation where τ = v/ωt, E = v/ωL ², σS = vN²/κω² and 2α² = σB²/pω. The problem is solved in the parameter range E?α²?1, α²/E?σS using the Laplace transform and two new spin-up times are obtained. Combined into one expression, they are τ = (1 + δ)α^?1E^-½ where δ = σμv. The spin-up mechanism is investigated and it is found that, in contrast to the homogeneous, conducting case, torsional Alfvén waves may be instrumental in the spin-up of a stratified conducting fluid. The effects of viscous and ohmic diffusion on the torsional Alfvén wave fronts are studied and the following regimes are identified: 0 < δ ?E/α², spin-up by meridional circulation of electric current with no Alfvén waves; E^-½/α ? δ ? 1, spin-up by meridional circulation of electric current with transient Alfvén waves; α/E^½ ? δ ? α²/E, spin-up by meridional circulation of current with weak Alfvén waves; 1 ? δ ? α/E^½, spin-up by strong Alfvén waves; α^½/E ? δ ? α²/E, spin-up by viscous diffusion with transient Alfvén waves; α/E ? δ < ∞, spin-up by viscous diffusion with no Alfvén waves.     

MAGNETOTELLURIC RESPONSE ON VERTICALLY INHOMOGENEOUS EARTH HAVING CONDUCTIVITY VARYING LINEARLY WITH DEPTH*

Magnetotelluric response is studied for an inhomogeneous medium having conductivity varying linearly with depth as σ(z) =σ₁+αz. For a medium having conductivity increasing linearly with depth, the phase of the impedance approaches 60° at long periods and the apparent resistivity becomes log (ρ_a) = 2 log (1.36/α^1/3) — 1/3 log (T'). The asymptote of log (ρ_a, T'→∞) when plotted against log (T') has a constant gradient —1/3 and has an intercept on the log (T') axis, which equals 6 log (1.36/α^1/3). When a homogeneous layer with a moderate thickness overlies an inhomogeneous half-space, this layer does not affect the asymptote, but it affects the cut-off period and pushes this toward the long period direction. For a medium having conductivity decreasing linearly with depth, the impedance is equivalent to that of a Cagniard two-layer model; the intercept period related to the thickness is T'₀=σ₁(h₂/2)². Homogeneous multilayer approximations to an inhomogeneous layer are also investigated, and it is shown that the fit to the model variation depends on the number of layers and the layer parameters chosen.     

Spin-up of a two-layer fluid in a rotating cylinder

Abstract

we report the results of experiments on the spin-up of two layers of immiscible fluid with a free upper surface in a rotating cylinder over a wide range of internal Froude numbers. Observations of the evolution of the velocity field by particle tracking indicates that spin-up of the azimuthal velocity in the upper layer take much longer than in a homogeneous fluid. Initially, spin-up occurs at a rate comparable to that of homogeneous fluid but, at high internal Froude number, a second phase follows in which the remaining lative motion decays much more slowly. Quantitative comparison of these measurements to the theory of Pedlosky (1967) shows good agreement.

Visualization of the interface displacement during spin-up detected the presence of transient azimuthal variations in the interface elevation over a wide range of Froude (F), Ekman (E), and Rossby (ε) number. nalysis of the occurrence of the asymmetric variations using the parameter space (Q, F), where Q = E ^1/2/ε, suggested by the baroclinic instability theory and experiments of Hart (1972), showed that the flow was stable for Q > 0.06 with no discernable dependence on F. This result, together with the prediction of Pedlosky's theory that radial gradient of potential vorticity in the two layers have opposite signs, suggests at the baroclinic instability mechanism was responsible for the asymmetries. The location and timing of these instabilities may account for the discrepancies between the observations and the Pedlosky (1967) theory.     

Asymptotic treatment of the stability of a rotating layer of fluid with rigid boundaries

Abstract

Chandrasekhar (1961) has summarized the stability results of Bénard convection in a rotating fluid for the cases where the boundary surfaces are both rigid and free, and for both exchange of stabilities and overstability. His analysis provides very accurate results for a limited range of Taylor number J. Bisshopp and Niiler (1965) presented an asymptotic analysis of the rigid boundary problem for exchange of stabilities which is valid for very large Taylor number. The present paper makes use of modern rotating fluid theory to develop an approximate scheme for evaluating the Rayleigh number and other parameters and variables. Known asymptotic results for the free boundary problem at large J are used and an expansion in powers of E^1/6 (the Ekman number, E = 2J ^?½) yields a sequence of equations and appropriate boundary conditions for the rigid boundary problem. After the algorithm for the calculation is developed, results are given for the problem to second order in the expansion parameter for the case of exchange of stabilities and to first order in the expansion parameters for the overstable case. Ekman boundary layers are important in the development as one might expect. However, an additional, diffusive boundary layer of thickness E^? is necessary to provide the details of the temperature field. This boundary layer is the thermal response in the vertical direction to the horizontal spacing of the cells which is also order E^?. The horizontal spacing of the cells is essentially a series of contiguous, Stewartson (1957) layers of thickness E^?.     

Rotating flow over long shallow ridges

Abstract

The flow of a rotating homogeneous, incompressible fluid past a long ridge is investigated. An analysis is presented for flows in which E ? 1, Ro ～ E^½, H/D ～ E⁰, h/D ～ E^½ and cosα ～ E⁰ where E is the Ekman number, Ro the Rossby number, H/D the fluid depth to ridge width ratio, h/D the ridge height to ridge width ratio and α the angle between the free stream flow and a line perpendicular to the ridge axis. The analysis includes effects of the nonlinear inertial terms. Particular examples of a ridge of triangular cross section and a sinusoidal topography are investigated in some detail. Experiments are presented for a triangular ridge which are in good agreement with the theory.     

Space-Time scales of internal waves

Abstract

We have contrived a model E(αω) α μ^?1ω^?p+1(ω ²?ω _i ²)^?+ for the distribution of internal wave energy in horizontal wavenumber, frequency-space, with wavenumber α extending to some upper limit μ(ω) α ω ^r-1 (ω ²?ω _i ²)^½, and frequency ω extending from the inertial frequency ω _i to the local Väisälä frequency n(y). The spectrum is portrayed as an equivalent continuum to which the modal structure (if it exists) is not vital. We assume horizontal isotropy, E(α, ω) = 2παE(α₁, α₂, ω), with α₁, α₂ designating components of α. Certain moments of E(α₁, α₂, ω) can be derived from observations. (i) Moored (or freely floating) devices measuring horizontal current u(t), vertical displacement η(t),…, yield the frequency spectra F _(u,η,…)(ω) = ∫∫ (U ², Z ²,…)E(α₁, ∞₂, ω) dα₁ dα₂, where U, Z,… are the appropriate wave functions. (ii) Similarly towed measurements give the wavenumber spectrum F _(…)(α1) = ∫∫… dα₂ dω. (iii) Moored measurements horizontally separated by X yield the coherence spectrum R(X, ω) which is related to the horizontal cosine transform ∫∫ E(α1, α2 ω) cos α₁ Xdα₁ dα₁. (iv) Moored measurements vertically separated by Y yield R(Y, ω) and (v) towed measurements vertically separated yield R(Y, α₁), and these are related to similar vertical Fourier transforms. Away from inertial frequencies, our model E(α, ω) α ω ^?p-r for α ≦ μ ω ω ^r, yields F(ω) ∞ ω ^?p, F(α₁) ∞ α₁ ^?q, with q = (p + r ? 1)/r. The observed moored and towed spectra suggest p and q between 5/3 and 2, yielding r between 2/3 and 3/2, inconsistent with a value of r = 2 derived from Webster's measurements of moored vertical coherence. We ascribe Webster's result to the oceanic fine-structure. Our choice (p, q, r) = (2, 2, 1) is then not inconsistent with existing evidence. The spectrum is E(∞, ω) ∞ ω ^?1(ω ²?ω _i ^{2 ?1}, and the α-bandwith μ ∞ (ω ²?ω _i ²)⁺ is equivalent to about 20 modes. Finally, we consider the frequency-of-encounter spectra F([sgrave]) at any towing speed S, approaching F(ω) as S ≦ S _o, and F(α₁) for α₁ = [sgrave]/S as S ≧ S _o, where S _o = 0(1 km/h) is the relevant Doppler velocity scale.     

On hydromagnetic instabilities and the mean electromotive force in a non-uniformly stratified Earth’s core affected by viscosity

Linear magnetoconvection in a model of a non-uniformly stratified horizontal rotating fluid layer with a toroidal magnetic field is investigated for no-slip and finitely electrically conductive boundaries and with very thin stably stratified upper sublayer. The basic parabolic temperature profile is determined by the temperature difference between the boundaries and by the homogeneous heat source distribution in the layer. This results in a density pattern, in which a stably stratified upper sublayer is present. The developed diffusive perturbations (modes) are strongly affected by the complicated coupling of viscous, thermal and magnetic diffusive processes. The calculations were performed for various values of Roberts number (q ≪ 1 and q = O(1)). The mean electromotive force produced by the developed hydromagnetic instabilities is investigated to find the modes, which can be appropriate for creating the α-effect. It was found that the azimuthal part of the EMF is dominant for westward modes when the Elsasser number Λ ≲ O(1).     

Resistive instabilities in rapidly rotating fluids: Linear theory of the convective modes

Abstract

This paper analyzes the linear stability of a rapidly-rotating, stratified sheet pinch in a gravitational field, g, perpendicular to the sheet. The sheet pinch is a layer (O ? z ? d) of inviscid, Boussinesq fluid of electrical conductivity σ, magnetic permeability μ, and almost uniform density ρ _o; z is height. The prevailing magnetic field. B _o(z), is horizontal at each z level, but varies in direction with z. The angular velocity, Ω, is vertical and large (Ω ? V_A/d, where V_A = B₀√(μρ₀) is the Alfvén velocity). The Elsasser number, Λ = σB² ₀/2Ωρ₀, measures σ. A (modified) Rayleigh number, R = gβd²/ρ₀V² _A, measures the buoyancy force, where β is the imposed density gradient, antiparallel to g. A Prandtl number, P_K = μσK, measures the diffusivity, k, of density differences.     

The structure of separated flow past a circular cylinder in a rotating frame

Abstract

This paper examines the detailed E ^1/4-layer structure of separated flow past a circular cylinder in a low-Rossby-number rotating fluid as the Ekman number E tends to zero. This structure is based on an initial proposal by Page (1987) but with some modifications in response to further evidence, outlined both in this paper and elsewhere, on the behaviour of E ^1/4-layer flows in this context. Numerical calculations for flow in an E ^1/4 shear layer along the separated free streamline are described and the mass flux from this layer is then used to calculate the higher-order flow within the separation bubble. The flow structure is found to have two forms, depending on the value of the O(1) parameter λ, and these are compared with results from published “Navier-Stokes” type calculations for the flow at small but finite values of E.     

The equatorial dynamics of a deep homogeneous ocean

Abstract

The flow properties of an homogeneous fluid which is bounded by two concentric spheres and two meridional planes which intersect along a diameter of the spheres are investigated. The spheres rotate about this diameter with slightly different angular velocities. As in the axisymmetric case studied by Proudman (1956) and Stewartson (1966) the viscous terms in the equations of motion are important only in boundary layers on the spheres and on the cylinder C which circumscribes the inner sphere and which has generators parallel to the axis of rotation, provided the Ekman number E is small. In the inviscid region the velocities are independent of the coordinate measuring distance along the axis of rotation and are much weaker, by a factor 0(E ^½), than the velocities in the Ekman layer on the driving surface (outer sphere). (It is assumed that the reference frame is fixed in the slower rotating inner sphere.) If the separation of the spheres is small compared to their radii then the asymmetric circulation inside C is characterized by an intense jet along the western wall. Loss of fluid from this jet sustains the eastward and northward flow in the inviscid interior where motion is driven by the suction of the Ekman layer on the outer sphere. (Geophysical conventions have been adopted.) Outside C an intense current is present on the eastern, not western, wall while motion in the inviscid region is westward, and away from the axis of rotation. Though there is no transport across C in the inviscid region, the meridional transport of the Ekman layer on the outer sphere is continuous across C and increases, through suction, as the equator is approached until it drains into an eastward flowing equatorial current of width 0(E ^1/7). The eastern boundary current outside C and shear layers on C carry this fluid to the intersection of C and the western wall where it feeds the western boundary current inside C.

The relation between this study and the experiments of Baker and Robinson (1970) is discussed.     

An analytical model for ice-edge upwelling

Abstract

This paper presents an analytical, two-dimensional model of the wind-induced homogeneous circulation near the edge of an ice pack floating on the ocean surface. It is shown that a vertical shear layer arises under the ice edge, by which the wind-driven geostrophic motion in the open ocean is matched to the flow region underneath the ice. As in coastal upwelling models, this shear layer consists of a thin E ^1/2-layer inside a thicker E ^1/4-layer (E being the Ekman number). Under certain conditions the shear layer produces a vertical mass flux from the bottom to the surface Ekman layer. Near the surface this upwelling flux is concentrated in the narrow E ^1/2-layer. Comparison with observations of upwelling at the edge of a polar ice pack shows good agreement.     

关于非均匀介质中电导率按深度变化时的势及其视电阻率函数

本文共分六部分内容。在前言中,简述了电测深方面有关研究的情况,尤其叙述了E.S.Sampaio（1976）的结果。 因本文采用的方法不同,从而求解了适合要求的微分方程定解问题,获得了电导率以 （其中α、σ0、σ1均为正实数,β为实数）按深度变化时势的分布。 其次,当β=0时,对电导率σ=（σ0+σ1z）α按深度变化时势的分布,应用Bessel函数的积分表达时,明显地改进和推广了E.S.Sampaio（1976）结果中的重要部分。 最后,结合实际作出了第一类、第二类标准化视电阻率函数。这种视电阻率函数,推广了Wenner和Schlumberger分别作出的视电阻率函数的有关部分,对电阻率法和激发极化法,都有其理论价值和实际意义。     

The electric field produced in the mantle by the dynamo in the core

A rigorous singular perturbation theory is developed to estimate the electric field E produced in the mantle M by the core dynamo when the electrical conductivity σ in M depends only on radius r, and when |r?_rln σ| ? 1 in most of M. It is assumed that σ has only one local minimum in M, either (a) at the Earth's surface ?V, or (b) at a radius b inside the mantle, or (c) at the core-mantle boundary ?K. In all three cases, the region where σ is no more than e times its minimum value constitutes a thin critical layer; in case (a), the radial electric field E_r ≈ 0 there, while in cases (b) and (c), E_r is very large there. Outside the critical layer, E_r ≈ 0 in all three cases. In no case is the tangential electric field E_S small, nearly toroidal, or nearly calculable from the magnetic vector potential A as ??_tA_S. The defects in Muth's (1979) argument which led him to contrary conclusions are identified. Benton (1979) cited Muth's work to argue that the core-fluid velocity u just below ?K can be estimated from measurements on ?V of the magnetic field B and its time derivative $\text{?}{}_{\mathrm{t}}\text{B}$. A simple model for westward drift is discussed which shows that Benton's conclusion is also wrong.In case (a), it is shown that knowledge of σ in M is unnecessary for estimating E_S on ?K with a relative error |r?_r 1nσ|^?1from measurements of E_S on ?V and knowledge of ?_tB in M (calculable from ?_tB on ?V if σ is small). Then, in case (a), u just below ?K can be estimated as $\text{?}\text{r}\text{×}\text{E}{}_{\mathrm{S}}\text{/B}{}_{\mathrm{r}}$. The method is impractical unless the contribution to E_S on ?V from ocean currents can be removed.The perturbation theory appropriate when σ in M is small is considered briefly; smallness of σ and of |r?_r ln σ|^?1 a independent questions. It is found that as σ → 0, B approaches the vacuum field in M but E does not; the explanation lies in the hydromagnetic approximation, which is certainly valid in M but fails as σ → 0. It is also found that the singular perturbation theory for |r?_r ln σ|^?1 is a useful tool in the perturbation calculations for σ when both σ and |r?_r ln σ|^?1 are small.     

α2-Dynamos and taylor's constraint

Abstract

An idealised α²ω-dynamo is considered in which the α-effect is prescribed. The additional ω-effect results from a geostrophic motion whose magnitude is determined indirectly by the Lorentz forces and Ekman suction at the boundary. As the strength of the α-effect is increased, a critical value α? _c is reached at which dynamo activity sets in; α? _c is determined by the solution of the kinematic α²-dynamo problem. In the neighbourhood of the critical value of α? the magnetic field is weak of order E ^1/4(μηρω)^½ due to the control of Ekman suction; E(?1) is the Ekman number. At certain values of α?, viscosity independent solutions are found satisfying Taylor's constraint. They are identified by the bifurcation of a nonlinear eigenvalue problem. Dimensional arguments indicate that following this second bifurcation the magnetic field is strong of order (μηρω)^½. The nature of the transition between the kinematic linear theory and the Taylor state is investigated for various distributions of the α-effect. The character of the transition is found to be strongly model dependent.     

Magnetohydrodynamic disturbances in the earth's core,IV

Summary In Paper III (Mohandis [1]²) we considered the sudden introduction of amagnetic dipole in the earth's core to act as a source of disturbance to the exitation field taken as a poloidal one. A symmetrical case was considered where the dipole axis is placed parallel to the original field and perpendicular to the earth's mantle. In the present work, we consider an unsymmetric case where the axis of themagnetic dipole is placed perpendicular to both the mantle and the exitation field which is taken as a toroidal one. A mathematical study is made for the resulting fluid motion in the core as well as for the generated hydromagnetic perturbations in both the mantle and the earth's fluid core. A more powerful method has been adopted than those used in previous cases.     

A nonlinear,resistive boundary layer in rotating hydromagnetic flow

A new nonlinear boundary layer in rotating hydromagnetic flows is presented. The purpose of this layer is to provide a smooth transition between the magnetic field at a rigid, electrically insulating boundary and that far from the boundary. The boundary layer problem is solved in a cylindrical geometry, assuming that the variables obey the hydromagnetic analog of the von Karman similarity. The primary momentum balance within this boundary layer is between pressure, Coriolis and hydromagnetic forces; the viscous and nonlinear inertia forces are unimportant. Diffusion of the magnetic field is balanced by advection where the advecting fluid flow is induced by the action of the hydromagnetic body forces within the boundary layer. The structure of the layer depends upon a single parameter which measures the strength of the normal magnetic field. Steady solutions are presented and their uniqueness and temporal stability are analyzed. The relevance of this layer to the core of the Earth is discussed. It is estimated that this boundary layer would be at least as thick as one tenth of the core radius if it exists within the core.     

Influence of advection on the characteristics of turbulence over uneven surface in the oasis and the Gobi Desert

Utilizing experimental data of the atmospheric surface layer in the Gobi Oasis of Jinta in a comparative study, we demonstrate that under the condition of unstable stratification, the normalization variances of temperature in the oasis and Gobi Desert meet 1/3(z/Λ)(z/Λ)while normalization variances of both humidity and CO2 in the oasis meet(z)sz\Λ-1/3;s s z Λ the normalization variance of temperature in the oasis is large due to disturbance by advection, whereas variance of CO2 in the Gobi Desert has certain degree of deviation relative to Monin-Obukhov(M-O) scaling, and humidity variance completely deviates from variance M-O scaling. The above result indicates that under the condition of advection, humidity variance meets the relationsm sA sB D and it is determined by relative magnitude of scalar variance of advection transport. Our study reveals that, if the scalar variance of humidity or CO2 transported by advection is much larger than local scalar variance, observation value of scalar variance will deviate from M-O scaling; when scalar variance of advection transport is close to or less than local scalar variance, the observation value of scalar variance approximately meets M-O scaling.     

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.     

Does measuring azimuthal variations in sap flux lead to more reliable stand transpiration estimates?

Stand transpiration (E) estimated using the sap‐flux method includes uncertainty induced by variations in sap flux (F) within a tree (i.e. radial and azimuthal variations) and those between trees. Unlike radial variations, azimuthal variations are not particularly systematic (i.e. higher/lower F is not always recorded for a specific direction). Here, we present a theoretical framework to address the question on how to allocate a limited number of sensors to minimize uncertainty in E estimates. Specifically, we compare uncertainty in E estimates for two cases: (1) measuring F for two or more directions to cover azimuthal variations in F and (2) measuring F for one direction to cover between‐tree variations in F. The framework formulates the variation in the probability density function for E (σ_E) based on F recorded in m different azimuthal directions (e.g. north, east, south and west). This formula allows us to determine the m value that minimizes σ_E. This study applied the framework to F data recorded for a 55‐year‐old Cryptomeria japonica stand. σ_E for m = 1 was found to be less than the values for m = 2, 3 and 4. Our results suggest that measuring F for one azimuthal direction provides more reliable E estimates than measuring F for two or more azimuthal directions for this stand, given a limited number of sensors. Application of this framework to other datasets helps us decide how to allocate sensors most effectively. Copyright © 2016 John Wiley & Sons, Ltd.     

ELECTRICAL POTENTIAL DUE TO A POINT CURRENT SOURCE OVER AN INHOMOGENEOUS ANISOTROPIC EARTH*

A mathematical expression for potential of a direct current point source in an inhomoge-neous anisotropic earth is derived. The coefficient of anisotropy is given by f= (σ_r/σ_z), where σ_r and σ_z are the conductivities parallel and perpendicular to the bedding plane. It is assumed that σ_z varies with depth, whereas σ_r varies transversely. This potential may be useful in interpretation of geoelectrical data in specified geological situations. Master curves for Wenner and Schlumberger configurations are presented