Thermal convection in differentially rotating systems

Abstract

The onset of convection in a cylindrical fluid annulus is analyzed in the case when the cylindrical walls are rotating differentially, a temperature gradient in the radial direction is applied, and the centrifugal force dominates over gravity. The small gap approximation is used and no-slip conditions on the cylindrical walls are assumed. It is found that over a considerable range of the parameter space either convection rolls aligned with the axis of rotation or rolls in the perpendicular (azimuthal) direction are preferred. It is shown that by a suitable redefinition of parameters, results for finite amplitude Taylor vortices and for convection rolls in the presence of shear can be applied to the present problem. Weakly nonlinear results for transverse rolls in a Couette flow indicate the possibility of subcritical bifurcation for Prandtl numbers P less than 0.82. Heat and momentum transports are derived as functions of P and the problem of interaction between transverse and longitudinal rolls is considered. The relevance of the analysis for problems of convection in planetary and stellar atmospheres is briefly discussed.     

Magnetic tracers,a probe of the solar convective layers

Abstract

A meridional circulation of sunspots has been measured through the digital analysis of the Meudon spectroheliograms from 1978 to 1983. Old and young sunspots follow a zonal meridional circulation, in several bands of latitude, in which two adjacent bands have opposite motions. This meridional circulation pattern is time-dependent. Using the H _α filaments as magnetic field tracers, a large-scale magnetic pattern has been found that was also obtained independently by direct measurement of the magnetic field (Hoeksema, 1988).

The coincidence of a large-scale magnetic pattern with a zonal meridional circulation suggests the existence of azimuthal rolls below the surface, and these azimuthal rolls can explain a number of properties of the solar cycle. New rolls occur with increasing proximity to the Equator, thereby indicating the direction of propagation of the dynamo wave. The occurrence of rolls is very favorable to the emergence of the magnetic regions. The rolls also influence the magnetic complexity of the active regions. They modulate the surface rotation through the Coriolis force, which accelerates or decelerates the fluid particles. They therefore offer a plausible explanation of the torsional oscillation pattern.

There are a number of problems raised by such an unexpected circulation pattern: for example, the coexistence of axisymmeric rolls with hypothetical giant cells, the location of the dynamo source below or within the convective zone, and the coupling of the radiative interior and the convective layers. To resolve these important issues, continuous observational studies are needed of the manifestation of solar activity, as well as of radius and luminosity variations. So, we have aimed our paper at an audience of theoreticians in the hope that they take up the challenges we describe.     

Instability of flows near the onset of convection in a rotating layer with stress-free horizontal boundaries

We investigate instability of convective flows of simple structure (rolls, standing and travelling waves) in a rotating layer with stress-free horizontal boundaries near the onset of convection. We show that the flows are always unstable to perturbations, which are linear combinations of large-scale modes and short-scale modes, whose wave numbers are close to those of the perturbed flows. Depending on asymptotic relations of small parameters α (the difference between the wave number of perturbed flows and the critical wave number for the onset of convection) and ε (ε² being the overcriticality and the perturbed flow amplitude being O(ε)), either small-angle or Eckhaus instability is prevailing. In the case of small-angle instability for rolls the largest growth rate scales as ε^8/5, in agreement with results of Cox and Matthews (Cox, S.M. and Matthews, P.C., Instability of rotating convection. J. Fluid. Mech., 2000, 403, 153–172) obtained for rolls with k = k _c . For waves, the largest growth rate is of the order ε^4/3. In the case of Eckhaus instability the growth rate is of the order of α².     

Topology of Rayleigh–Bénard convection and magnetoconvection in plane layer

ABSTRACT

The present study aims to link the dynamics of geophysical fluid flows with their vortical structures in physical space and to study the transition of these structures due to the control parameters. The simulations are carried in a rectangular box filled with liquid gallium for three different cases, namely, Rayleigh–Bénard convection (RBC), magnetoconvection (MC) and rotating magnetoconvection (RMC). The physical setup and material properties are similar to those considered by Aurnou and Olson in their experimental work. The simulated results are validated with theoretical results of Chandrasekhar and experimental results of Aurnou and Olson. The results are also topologically verified with the help of Euler number given by Ma and Wang. For RBC, the onset is obtained at Ra greater than 1708 and at this Ra, the symmetric rolls are orientated in/along a horizontal axis. As the value of Ra increases further, the width of the horizontal rolls starts to amplify. It is observed that these two-dimensional rolls are nothing but the cross-sections of three-dimensional (3D) cylindrical rolls with wave structures. When the vertically imposed magnetic field is added to RBC, the onset of convection is delayed due to the effect of Lorentz force on the thermal buoyancy force. The presence of 3D rectangular structures is highlighted and analysed. When the magnetically influenced rectangular box rotates about vertical axis at low rotation rates in magnetoconvection model, the onset of convection gets further delayed by magnetic field, which is in general agreement with the theoretical predictions. The critical Ra increases linearly with magnetic field intensity. Coherent thermal oscillations are detected near the onset of convection, at moderate rotation rates.     

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.     

Convective instability when the temperature gradient and rotation vector are oblique to gravity. II. Real fluids with effects of diffusion

Abstract

The linear stability analysis of Hathaway, Gilman and Toomre (1979) (hereafter referred to as Paper I) is repeated for Boussinesq fluids with viscous and thermal diffusion. As in Paper I the fluid is confined between plane parallel boundaries and the rotation vector is oblique to gravity. This tilted rotation vector introduces a preference for roll-like disturbances whose axes are oriented north-south; the preference is particularly strong in the equatorial region. The presence of a latitudinal temperature gradient produces a thermal wind shear which favors axisymmetric convective rolls if the gradient exceeds some critical value. For vanishingly small diffusivities the value of this transition temperature gradient approaches the inviscid value found in Paper I. For larger diffusivities larger gradients are required particularly in the high latitudes. These results are largely independent of the Prandtl number. Diffusion tends to stabilize the large wavenumber rolls with the result that a unique wavenumber can be found at which the growth rate is maximized. These preferred rolls have widths comparable to the depth of the layer and tend to be broader near the equator. The axisymmetric rolls are similar in many respects to the cloud bands on Jupiter provided they extend to a depth of about 15,000 km.     

Effective conductivity of periodic media with cuboid inclusions

This paper presents a numerical solution for the effective conductivity of a periodic binary medium with cuboid inclusions located on an octahedral lattice. The problem is defined by five dimensionless geometric parameters and one dimensionless conductivity contrast parameter. The effective conductivity is determined by considering the flow through the “elementary flow domain” (EFD), which is an octant of the unitary domain of the periodic media. We derive practical bounds of interest for the six-dimensional parameter space of the EFD and numerically compute solutions at regular intervals throughout the entire bounded parameter space. A continuous solution of the effective conductivity within the limits of the simulated parameter space is then obtained via interpolation of the numerical results. Comparison to effective conductivities derived for random heterogeneous media demonstrate similarities and differences in the behavior of the effective conductivity in regular periodic (low entropy) vs. random (high entropy) media. The results define the low entropy bounds of effective conductivity in natural media, which is neither completely random nor completely periodic, over a large range of structural geometries. For aniso-probable inclusion spacing, the absolute bounds of K_eff for isotropic inclusions are the Wiener bounds, not the Hashin-Shtrikman bounds. For isotropic inclusion and isoprobable conditions well below the percolation threshold, the results are in agreement with the self-consistent approach. For anisotropic cuboid inclusions, or at relatively close spacing in at least one direction (p > 0.2) (aniso-probable conditions), the effective conductivity of the periodic media is significantly different from that found in anisotropic random binary or Gaussian media.     

Displacement and rectification of planetary fluids

Abstract

Laboratory experiments on the decay (spin-up) of fluid motion on the β-plane are compared with theory. Under weakly dissipative conditions, some particles conserve potential vorticity during the entire decay. We also study the rectified mean flow which is produced by the lateral Reynolds Stress when a low frequency force is applied to the planetary fluid. The possible connection of effects with oceanic phenomena is briefly discussed.     

A double-porosity model for a fractured aquifer with non-Darcian flow in fractures

Abstract

Non-Darcian flow in a finite fractured confined aquifer is studied. A stream bounds the aquifer at one side and an impervious stratum at the other. The aquifer consists of fractures capable of transmitting water rapidly, and porous blocks which mainly store water. Unsteady flow in the aquifer due to a sudden rise in the stream level is analysed by the double-porosity conceptual model. Governing equations for the flow in fractures and blocks are developed using the continuity equation. The fluid velocity in fractures is often too high for the linear Darcian flow so that the governing equation for fracture flow is modified by Forcheimer's equation, which incorporates a nonlinear term. Governing equations are coupled by an interaction term that controls the quasi-steady-state fracture—block interflow. Governing equations are solved numerically by the Crank-Nicolson implicit scheme. The numerical results are compared to the analytical results for the same problem which assumes Darcian flow in both fractures and blocks. Numerical and analytical solutions give the same results when the Reynolds number is less than 0.1. The effect of nonlinearity on the flow appears when the Reynolds number is greater than 0.1. The higher the rate of flow from the stream to the aquifer, the higher the degree of nonlinearity. The effect of aquifer parameters on the flow is also investigated. The proposed model and its numerical solution provide a useful application of nonlinear flow models to fractured aquifers. It is possible to extend the model to different types of aquifer, as well as boundary conditions at the stream side. Time-dependent flow rates in the analysis of recession hydrographs could also be evaluated by this model.     

Application of Generalized Likelihood Uncertainty Estimation (GLUE) at different temporal scales to reduce the uncertainty level in modelled river flows

ABSTRACT

In this study, the distributed catchment-scale model, DiCaSM, was applied on five catchments across the UK. Given its importance, river flow was selected to study the uncertainty in streamflow prediction using the Generalized Likelihood Uncertainty Estimation (GLUE) methodology at different timescales (daily, monthly, seasonal and annual). The uncertainty analysis showed that the observed river flows were within the predicted bounds/envelope of 5% and 95% percentiles. These predicted river flow bounds contained most of the observed river flows, as expressed by the high containment ratio, CR. In addition to CR, other uncertainty indices – bandwidth B, relative bandwidth RB, degrees of asymmetry S and T, deviation amplitude D, relative deviation amplitude RD and the R factor – also indicated that the predicted river flows have acceptable uncertainty levels. The results show lower uncertainty in predicted river flows when increasing the timescale from daily to monthly to seasonal, with the lowest uncertainty associated with annual flows.     

An approximate method for determining the stability of two-scale flow in the mantle

Summary An approximate method of calculating the stability of two-dimensional convection rolls to cross-roll disturbances is evaluated by comparison with the results from exact analyses and good agreement is obtained. In particular, investigation of the second mode of the disturbance provides a qualitative estimate of the terms excluded in making the approximation. We conclude that this approximate method, sensibly used, gives a good indication of the stability of rolls to cross-roll disturbances.Application of this method is made to convection rolls heated partially or wholly from within, and to rolls in the presence of a longitudinal shear flow. The approximation indicates that, at high Rayleigh Numbers, growth rates are underestimated and so the amount of shear calculated to stabilise longitudinal rolls is a lower bound. Our results suggest that two-dimensional rolls are unlikely to be a stable flow pattern beneath even the fastest moving plates.     

Study of forward–backward solute dispersion profiles in a semi-infinite groundwater system

ABSTRACT

Forward–backward solute dispersion with an intermediate point source in one-dimensional semi-infinite homogeneous porous media is studied in this paper. Solute transport under sorption conditions, first-order decay and zero-order production terms are included. The first type of boundary condition is taken as a constant point source at an intermediate point from where forward and backward solute dispersion is examined. The Laplace transform method is adopted to solve the governing equation analytically. All the analytical results are obtained in graphical form to investigate the forward–backward solute transport in porous media for various hydrological input data. The graphical nature of the analytical solution is compared with numerical data taken from existing literature and similar results are obtained. Also, numerical solution of the governing equation is obtained by the Crank-Nicolson finite difference scheme and validated with the analytical solution, which demonstrates good agreement between them. Accuracy of the solution is also observed by using RMSE.     

Structural uncertainty assessment in a discharge simulation model

Abstract

A major goal in hydrological modelling is to identify and quantify different sources of uncertainty in the modelling process. This paper analyses the structural uncertainty in a streamflow modelling system by investigating a set of models with increasing model structure complexity. The models are applied to two basins: Kielstau in Germany and XitaoXi in China. The results show that the model structure is an important factor affecting model performance. For the Kielstau basin, influences from drainage and wetland are critical for the local runoff generation, while for the XitaoXi basin accurate distributions of precipitation and evapotranspiration are two of the determining factors for the success of the river flow simulations. The derived model uncertainty bounds exhibit appropriate coverage of observations. Both case studies indicate that simulation uncertainty for the low-flow period contributes more to the overall uncertainty than that for the peak-flow period, although the main hydrological features in these two basins differ greatly.

Citation Zhang, X. Y., Hörmann, G., Gao, J. F. & Fohrer, N. (2011) Structural uncertainty assessment in a discharge simulation model. Hydrol. Sci. J. 56(5), 854–869.     

Energy spectra associated with long's model for steady finite-amplitude flow over orography

Abstract

The process of wave steepening in Long's model of steady, two-dimensional stably stratified flow over orography is examined. Under conditions of the long-wave approximation, and constant values of the background static stability and basic flow, Long's equation is cast into the form of a nonlinear advection equation. Spectral properties of this latter equation, which could be useful for the interpretation of data analyses under mountain wave conditions, are presented. The principal features, that apply at the onset of convective instability (density constant with height), are:

i) a power spectrum for available potential energy that exhibits a minus eight-thirds decay, in terms of the vertical wavenumber k _z ^-;

ii) a rate of energy transfer across the spectrum that is inversely proportional to the wavenumber for large k _z ^-;

iii) an equipartition between the kinetic energy of the horizontal motion and the available potential energy, under the longwave approximation, although all the disturbance energy is kinetic at the point where convective instability is initiated. It is also shown that features i) and ii) apply to more general conditions that are appropriate to Long's model, not just the long-wave approximation. Application to fully turbulent flow or to conditions at the onset of shearing instability are not considered to be warranted, since the development only applies to conditions at the onset of convective instability.     

Laboratory experiments in a baroclinic annulus with heating and cooling on the horizontal boundaries

Abstract

Experiments have been performed in a cylindrical annulus with horizontal temperature gradients imposed upon the horizontal boundaries and in which the vertical depth was smaller than the width of the annulus. Qualitative observations were made by the use of small, suspended, reflective flakes in the liquid (water).

Four basic regimes of flow were observed: (1) axisymmetric flow, (2) deep cellular convection, (3) boundary layer convective rolls, and (4) baroclinic waves. In some cases there was a mix of baroclinic and convective instabilities present. As a “mean” interior Richardson number was decreased from a value greater than unity to one less than zero, axisymmetric baroclinic instability of the Solberg type was never observed. Rather, the transition was from non-axisymmetric baroclinic waves, to a mix of baroclinic and convective instability, to irregular cellular convection.     

Indices for assessing the prediction bounds of hydrological models and application by generalised likelihood uncertainty estimation / Indices pour évaluer les bornes de prévision de modèles hydrologiques et mise en œuvre pour une estimation d'incertitude par vraisemblance généralisée

Abstract

To reflect the uncertainties of a hydrological model in simulating and forecasting observed discharges according to rainfall inputs, the estimated result for each time step should not be just a point estimate (a single numerical value), but should be expressed as a prediction interval, i.e. a band defined by the prediction bounds of a particular confidence level α. How best to assess the quality of the prediction bounds thus becomes very important for understanding the modelling uncertainty in a comprehensive and objective way. This paper focuses on seven indices for characterizing the prediction bounds from different perspectives. For the three case-study catchments presented, these indices are calculated for the prediction bounds generated by the generalized likelihood uncertainty estimation (GLUE) method for various threshold values. In addition, the relationships among these indices are investigated, particularly that of the containing ratio (CR) to the other indices. In this context, three main findings are obtained for the prediction bounds estimated by GLUE. Firstly, both the average band-width and the average relative band-width are seen to have very strong linear correlations with the CR index. Secondly, a high CR value, a narrow band-width, and a high degree of symmetry with respect to the observed hydrograph, all of which are clearly desirable properties of the prediction bounds estimated by the uncertainty assessment methods, cannot all be achieved simultaneously. Thirdly, for the prediction bounds considered, the higher CR values and the higher degrees of symmetry with respect to the observed hydrograph are found to be associated with both the larger band-widths and the larger deviation amplitudes. It is recommended that a set of different indices, such as those considered in this study, be employed for assessing and comparing the prediction bounds in a more comprehensive and objective way.     

Note on the flow of a stratified fluid over a stationary obstacle in a channel

Abstract

A new method is introduced to produce a uniform stratified flow over a stationary obstacle in an open channel. The flow is achieved by discharging the flow from the channel through a sink. The details of the sink are unimportant. The flow speed is limited only by the sink capacity. Selective withdrawal at lower densimetric Froude numbers is effectively eliminated through the use of a contraction. The standing, free-surface, long wave arising from the initiation of the flow is also eliminated by the contraction. Experiments are conducted for flow over a sphere for a range of Reynolds numbers from O(10²) to O(10³) and a range of Richardson numbers from O(10^?1) to O(10). Dye and neutrally buoyant droplets are used for quantitative analysis of the wake structure. The wake is also probed by a hot-film anemometer. The frequency of vortex shedding is obtained. Comparison with data from towed experiments is also presented.     

On the generation,propagation and radiation of magneto–acoustic–gravity waves with application to stars

Abstract

The study of the mechanisms controlling the stratification in closed fluid regions is an important branch of geophysical fluid dynamics. Part of this subject can be handled with a simple linear model, consisting of a buoyancy layer at the non-horizontal boundaries of a container and an advective-diffusive interior coupled by volume continuity. The model is valid under the following conditions: firstly, the buoyancy-frequency characterizing the solution must be sufficiently large to give rise to a flow pattern of boundary layer type and, secondly, the non-horizontal walls must not have too large thermal conductivity.

The main purpose of the present paper is to summarise previous work done by the authors in this field and to present some consequences of their theory not previously discussed.

Three important cases are discussed; certain stationary solutions, the decay of a given stratification and the build up of a stratification in a homogeneous fluid. The experimental results concerning the afore-mentioned cases are presented.     

A transient phreatic surface mound,evidenced by a strip of vegetation on an earth dam

Abstract

On the slopes of the embankment of the Al-Khod groundwater recharge–flood protection dam (Oman), a band of scrub vegetation community emerged after torrential rains and temporary filling of the dam reservoir. Species composition differs markedly on both sides of the embankment, with many exotics found on the reservoir side and more typical gravel-desert species on the outside. Hydro-ecologically, the vegetation is interpreted as the footprint of a temporary storage of water, which is a small-sized groundwater mound within the permeable shoulder of the levee. The levee, as an anthropogenic landform, induces a U-turn (gravitational slumping–lateral seepage–transpirational moisture ascent) topology of seepage. The Lembke method of successive variations of steady states is used in modelling the water table dynamics. In the early stage of the mound decay, outflow through a seepage face of the shoulder is modelled by the Barenblatt slumping parabola of the phreatic-zone part of the flow domain, which is perfectly matched with the Youngs exact solution for a purely horizontal flow through a porous wedge. At the stitching cross-section, the flow rates and saturated depths in the two zones coincide. The late stage of mound evolution is characterized by transpiration by the plant roots projected onto a shrinking free surface, with the Barenblatt and Youngs solutions conjugated but without the outcrop of the saturated mound on the levee slope. Ordinary differential equations for the sliding or descending locus of the intersection of the parabola and the triangle hypotenuse are integrated in a closed form or by the Runge-Kutta method. The dwindling saturated volume and the rate of drainage are obtained. They can be used in assessments of the hydro-ecological sustainability of slope-rooted shrubs (vegetation survival between rare rainfall episodes).

Editor D. Koutsoyiannis; Associate editor A. Porporato     

Neural network model for discharge and water-level prediction for Ramganga River catchment of Ganga Basin,India

ABSTRACT

Discharges and water levels are essential components of river hydrodynamics. In unreachable terrains and ungauged locations, it is quite difficult to measure these parameters due to rugged topography. In the present study an artificial neural network model has been developed for the Ramganga River catchment of the Ganga Basin. The modelled network is trained, validated and tested using daily water flow and level data pertaining to 4 years (2010–2013). The network has been optimized using an enumeration technique and a network topology of 4-10-2 with a learning rate set at 0.06, which was found optimum for predicting discharge and water-level values for the considered river. The mean square error values obtained for discharge and water level for the tested data were found to be 0.046 and 0.012, respectively. Thus, monsoon flow patterns can be estimated with an accuracy of about 93.42%.

Editor M.C. Acreman; Associate editor E. Gargouri