On the onset of convection in rotating spherical shells

Abstract

The linear problem of the onset of convection in rotating spherical shells is analysed numerically in dependence on the Prandtl number. The radius ratio η=r _i/r _o of the inner and outer radii is generally assumed to be 0.4. But other values of η are also considered. The goal of the analysis has been the clarification of the transition between modes drifting in the retrograde azimuthal direction in the low Taylor number regime and modes traveling in the prograde direction at high Taylor numbers. It is shown that for a given value m of the azimuthal wavenumber a single mode describes the onset of convection of fluids of moderate or high Prandtl number. At low Prandtl numbers, however, three different modes for a given m may describe the onset of convection in dependence on the Taylor number. The characteristic properties of the modes are described and the singularities leading to the separation with decreasing Prandtl number are elucidated. Related results for the problem of finite amplitude convection are also reported.     

Convection in rotating annular channels heated from below. Part 1. Linear stability and weakly nonlinear mean flows

Convection in a Boussinesq fluid in an annular channel rotating about a vertical axis with lateral rigid sidewalls, stress-free top and bottom, uniformly heated from below is investigated. The sidewalls are assumed to be either perfectly insulating or conducting. Three different types of convection are identified when the channel is rotating sufficiently fast: (i) global oscillatory convection preferred for small Prandtl numbers in channels with intermediate or large aspect ratios (width to height ratio), (ii) wall-localized oscillatory convection representing the most unstable mode for moderate or large Prandtl numbers in channels with intermediate or large aspect ratios and (iii) global stationary convection preferred in channels with sufficiently small aspect ratios regardless of the size of the Prandtl number. The corresponding weakly nonlinear problem describing differential rotation and meridional circulation is also examined, showing that geostrophic, multiple-peaked (two prograde and two retrograde) differential rotation can be maintained by the Reynolds stresses in wall-localized convective eddies in a rapidly rotating channel.     

Magnetic field generation by convective flows in a plane layer: the dependence on the Prandtl numbers

Investigation of magnetic field generation by convective flows is carried out for three values of kinematic Prandtl number: P = 0.3, 1 and 6.8. We consider Rayleigh–Bénard convection in Boussinesq approximation assuming stress-free boundary conditions on horizontal boundaries and periodicity with the same period in the x and y directions. Convective attractors are modelled for increasing Rayleigh numbers for each value of the kinematic Prandtl number. Linear and non-linear dynamo action of these attractors is studied for magnetic Prandtl numbers P _m ≤ 100. Flows, which can act as magnetic dynamos, have been found for all the three considered values of P, if the Rayleigh number R is large enough. The minimal R, for which of magnetic field generation occurs, increases with P. The minimum (over R) of critical P_m for magnetic field generation in the kinematic regime is admitted for P = 0.3. Thus, our study indicates that smaller values of P are beneficial for magnetic field generation.     

The nonlinear development of disturbances in a zonal shear flow

Abstract

A study is made of the nonlinear stability of a weakly supercritical zonal shear flow in the β-plane approximation. The dynamics of initially small disturbances are examined. The main nonlinear effects are associated with the rearrangement of the critical layer. It is shown that as the wave grows in amplitude, linear regimes of the critical layer (viscous and nonstationary) change over to a nonlinear regime while the exponential law of disturbance growth becomes a power-law.     

Nonlinear modal analysis of penetrative convection

Abstract

The modal expansion procedure has been used to analyze penetrative convection that arises when a thin unstable layer is embedded between two stable regions. The Boussinesq approximation is applied in which the effect of compressibility and stratification are neglected. Various calculations have been made, with one and two modes, for Rayleigh numbers ranging from the critical value to more than 10⁵ times critical. The effect of decreasing the Prandtl number has also been investigated.

It is found that in the nonlinear regime, the convective motions penetrate substantially into the stable regions. The flux of kinetic energy plays a crucial role in such penetration, and its existence puts some requirements on the motions: in the single-mode case, they need to be three-dimensional. The extent of penetration amounts to about half of the thickness of the unstable layer on each side of it when the degree of instability and that of stability are comparable in the two domains; it increases as the stability of the outer region is lowered. The penetration depth appears to be independent of all other parameters defining the problem.     

准周期性磁结构的形成过程及其特征的数值分析

在入流马赫数ＭＡ为０．１－０．４范围内，选取不同的入流，逐一考察在可压缩导电介质中磁场重联的类型随磁雷诺数Ｒｍ的变化．大量的数值结果表明：在不同入流驱动下，磁雷诺数（Ｒｍ）均有临界值存在，如果等离子体系统逐渐趋于稳定的单Ｘ线重联；若间歇性的次级撕裂产生了重复出现的磁结构，而且Ｍａ越大产生间歇性次级撕裂的临界值越高，与ＭＡ之间基本上符合线性关系．此外还发现，当多重Ｘ线重联间歇性地出现时，有关物理量发生准周期性地振荡．     

准周期性磁结构的形成过程及其特征的数值分析

在入流马赫数ＭＡ为０．１－０．４范围内,选取不同的入流,逐一考察在可压缩导电介质中磁场重联的类型随磁雷诺数Ｒｍ的变化．大量的数值结果表明：在不同入流驱动下,磁雷诺数（Ｒｍ）均有临界值存在,如果等离子体系统逐渐趋于稳定的单Ｘ线重联;若间歇性的次级撕裂产生了重复出现的磁结构,而且Ｍａ越大产生间歇性次级撕裂的临界值越高,与ＭＡ之间基本上符合线性关系．此外还发现,当多重Ｘ线重联间歇性地出现时,有关物理量发生准周期性地振荡．     

Simulation of moist mountain waves with an anelastic model

Abstract

A two-dimensional, nonlinear, time-dependent, non-hydrostatic, anelastic, numerical model is used to assess the effect of condensation on the evolution and structure of gravity waves generated by the passage of a stable, moist stream over topography. Precipation is ignored but water phase changes are taken into account explicitly.

The main effect of condensation is to damp the wave intensity and to reduce the wave drag, which can be diminished by as much as 50% compared to its value in dry simulations. This result agrees with some earlier analytical models and some more recent fully compressible numerical models.

This model also confirms that the presence of condensation delays the overturning of isentropes, and the formation of the critical layer that accompanies wave-breaking.     

The Influence of the Prandtl Number on the Style of Vigorous Thermal Convection

The effect of the Prandtl number on convection in a planar three-dimensional geometry is investigated in this study. We have employed a numerical scheme to integrate the governing equations. Differently from previous studies we have chosen stress-free boundaries. Experiments have been performed at a Rayleigh number of Ra = 10 6 for Prandtl numbers (Pr) ranging from 0.025 to 100. We have further conducted one experiment in the limiting case of infinite Prandtl number. Despite the differences in the geometry and the boundary conditions, as compared to other studies, we find a similar transition in the dynamics of the flow when the Prandtl number is increased. While the velocity and the temperature structure show diffusive character at low Pr, sharp thermal boundary layers form at high Pr. The heat transport efficiency increases with Pr until a transition value is reached, from there on Nu behaves almost asymptotically. The transition can not be caused by a change in hierarchies between velocity and thermal boundary layers, as suggested in other studies. Due to the stress-free boundaries, a velocity boundary layer does not exist. We observe that the toroidal part of the flow is strong at low Pr and looses its strength with increasing Pr, thus it is likely to be responsible for the transition. In a further chapter we demonstrate that due to the neglect of the toroidal part in two-dimensional calculations at low Pr results are obtained which are misleading, even in a qualitative sense. Infinite Pr results from 2D calculations closely resemble the dynamics of fully 3D flows.     

The Boussinesq and anelastic liquid approximations for convection in the Earth’s core

Convection in the Earth’s core is usually studied in the Boussinesq approximation in which the compressibility of the liquid is ignored. The density of the Earth’s core varies from ICB to CMB by approximately 20%. The question of whether we need to take this variation into account in core convection and dynamo models is examined. We show that it is in the thermodynamic equations that differences between compressible and Boussinesq models become most apparent. The heat flux conducted down the adiabat is much smaller near the inner core boundary than it is near the core-mantle boundary. In consequence, the heat flux carried by convection is much larger nearer the inner core boundary than it is near the core-mantle boundary. This effect will have an important influence on dynamo models. Boussinesq models also assume implicitly that the rate of working of the gravitational and buoyancy forces, as well as the Ohmic and viscous dissipation, are small compared to the heat flux through the core. These terms are not negligible in the Earth’s core heat budget, and neglecting them makes it difficult to get a thermodynamically consistent picture of core convection. We show that the usual anelastic equations simplify considerably if the anelastic liquid approximation, valid if αT?1, where α is the coefficient of expansion and T a typical core temperature, is used. The resulting set of equations are not significantly more difficult to solve numerically than the usual Boussinesq equations. The relationship of our anelastic liquid equations to the Boussinesq equations is also examined.     

Seasonality of low flows and dominant processes in the Rhine River

Low flow forecasting is crucial for sustainable cooling water supply and planning of river navigation in the Rhine River. The first step in reliable low flow forecasting is to understand the characteristics of low flow. In this study, several methods are applied to understand the low flow characteristics of Rhine River basin. In 108 catchments of the Rhine River, winter and summer low flow regions are determined with the seasonality ratio (SR) index. To understand whether different numbers of processes are acting in generating different low flow regimes in seven major sub-basins (namely, East Alpine, West Alpine, Middle Rhine, Neckar, Main, Mosel and Lower Rhine) aggregated from the 108 catchments, the dominant variable concept is adopted from chaos theory. The number of dominant processes within the seven major sub-basins is determined with the correlation dimension analysis. Results of the correlation dimension analysis show that the minimum and maximum required number of variables to represent the low flow dynamics of the seven major sub-basins, except the Middle Rhine and Mosel, is 4 and 9, respectively. For the Mosel and Middle Rhine, the required minimum number of variables is 2 and 6, and the maximum number of variables is 5 and 13, respectively. These results show that the low flow processes of the major sub-basins of the Rhine could be considered as non-stochastic or chaotic processes. To confirm this conclusion, the rescaled range analysis is applied to verify persistency (i.e. non-randomness) in the processes. The estimated rescaled range statistics (i.e. Hurst exponents) are all above 0.5, indicating that persistent long-term memory characteristics exist in the runoff processes. Finally, the mean values of SR indices are compared with the nonlinear analyses results to find significant relationships. The results show that the minimum and maximum numbers of required variables (i.e. processes) to model the dynamic characteristics for five out of the seven major sub-basins are the same, but the observed low flow regimes are different (winter low flow regime and summer low flow regime). These results support the conclusion that a few interrelated nonlinear variables could yield completely different behaviour (i.e. dominant low flow regime).     

Finite amplitude convection and magnetic field generation in a rotating spherical shell

Abstract

Finite amplitude solutions for convection in a rotating spherical fluid shell with a radius ratio of η=0.4 are obtained numerically by the Galerkin method. The case of the azimuthal wavenumber m=2 is emphasized, but solutions with m=4 are also considered. The pronounced distinction between different modes at low Prandtl numbers found in a preceding linear analysis (Zhang and Busse, 1987) is also found with respect to nonlinear properties. Only the positive-ω-mode exhibits subcritical finite amplitude convection. The stability of the stationary drifting solutions with respect to hydrodynamic disturbances is analyzed and regions of stability are presented. A major part of the paper is concerned with the growth of magnetic disturbances. The critical magnetic Prandtl number for the onset of dynamo action has been determined as function of the Rayleigh and Taylor numbers for the Prandtl numbers P=0.1 and P=1.0. Stationary and oscillatory dynamos with both, dipolar and quadrupolar, symmetries are close competitors in the parameter space of the problem.     

Effect of fracture pressure depletion regimes on the dual-porosity shape factor for flow of compressible fluids in fractured porous media

A precise value of the matrix-fracture transfer shape factor is essential for modeling fluid flow in fractured porous media by a dual-porosity approach. The slightly compressible fluid shape factor has been widely investigated in the literature. In a recent study, we have developed a transfer function for flow of a compressible fluid using a constant fracture pressure boundary condition [Ranjbar E, Hassanzadeh H, Matrix-fracture transfer shape factor for modeling flow of a compressible fluid in dual-porosity media. Adv Water Res 2011;34(5):627-39. doi:10.1016/j.advwatres.2011.02.012]. However, for a compressible fluid, the consequence of a pressure depletion boundary condition on the shape factor has not been investigated in the previous studies. The main purpose of this paper is, therefore, to investigate the effect of the fracture pressure depletion regime on the shape factor for single-phase flow of a compressible fluid. In the current study, a model for evaluation of the shape factor is derived using solutions of a nonlinear diffusivity equation subject to different pressure depletion regimes. A combination of the heat integral method, the method of moments and Duhamel’s theorem is used to solve this nonlinear equation. The developed solution is validated by fine-grid numerical simulations. The presented model can recover the shape factor of slightly compressible fluids reported in the literature. This study demonstrates that in the case of a single-phase flow of compressible fluid, the shape factor is a function of the imposed boundary condition in the fracture and its variability with time. It is shown that such dependence can be described by an exponentially declining fracture pressure with different decline exponents. These findings improve our understanding of fluid flow in fractured porous media.     

Properties of nonlinear dynamo waves

Abstract

Dynamo theory offers the most promising explanation of the generation of the sun's magnetic cycle. Mean field electrodynamics has provided the platform for linear and nonlinear models of solar dynamos. However the nonlinearities included arc (necessarily) arbitrarily imposed in these models. This paper conducts a systematic survey of the role of nonlinearities in the dynamo process, by cousidering the behaviour of dynamo waves in the nonlinear regime. It is demonstrated that only by considering realistic nonlinearities that are non-local in space and time can modulation of the basic dynamo wave be achieved. Moreover this modulation is greatest when there is a large separation of timescales provided by including a low magnetic Prandtl number in the equation for the velocity perturbations.     

On the dynamics of 3-D single thermal plumes at various Prandtl numbers and Rayleigh numbers

Three-dimensional (3-D) numerical simulations of single turbulent thermal plumes in the Boussinesq approximation are used to understand more deeply the interaction of a plume with itself and its environment. In order to do so, we varied the Rayleigh and Prandtl numbers from Ra?～?10⁵ to Ra?～?10⁸ and from Pr?～?0.025 to Pr?～?70. We found that thermal dissipation takes place mostly on the border of the plume. Moreover, the rate of energy dissipation per unit mass ε _T has a critical point around Pr?～?0.7. The reason is that at Pr greater than ～0.7, buoyancy dominates inertia and thermal advection dominates wave formation whereas this trend is reversed at Pr less than ～0.7. We also found that for large enough Prandtl number (Pr?～?70), the velocity field is mostly poloidal although this result was known for Rayleigh–Bénard convection (see Schmalzl et al. [On the validity of two-dimensional numerical approaches to time-dependent thermal convection. Europhys. Lett. 2004, 67, 390--396]). On the other hand, at small Prandtl numbers, the plume has a large helicity at large scale and a non-negligible toroidal part. Finally, as observed recently in details in weakly compressible turbulent thermal plume at Pr?=?0.7 (see Plourde et al. [Direct numerical simulations of a rapidly expanding thermal plume: structure and entrainment interaction. J. Fluid Mech. 2008, 604, 99--123]), we also noticed a two-time cycle in which there is entrainment of some of the external fluid to the plume, this process being most pronounced at the base of the plume. We explain this as a consequence of calculated Richardson number being unity at Pr?=?0.7 when buoyancy balance inertia.     

On Convection Driven Dynamos in Rotating Spherical Shells

Solutions for chaotic dynamos driven by thermal convection in a rotating spherical shell are obtained numerically for different Prandtl numbers. The influence of this parameter which is usually suppressed in the magnetostrophic approximation is emphasized in the present analysis.     

On the upper bounding approach to thermal convection at moderate rayleigh numbers

Abstract

Two upper bounding problems for thermal convection in a layer of fluid contained between perfectly conducting stress-free boundaries are treated numerically. Since the Euler equations resulting from this variational approach are simpler than the Navier-Stokes equations, they allow numerical calculations to be carried out economically to fairly large values of the Rayleigh number. The upper bounding problem formulated by Howard (1963), which yields a Nusselt number independent of Prandtl number, diverges from the correct behavior as the Rayleigh number increases. In hopes of coming closer to results of previous investigations of the Boussinesq equations of motion, a more restrictive upper bounding problem is formulated. For large Prandtl numbers the momentum equation is linearized and is used as an explicit side constraint on the variational problem, thereby forcing the solutions to more closely resemble the solutions of the Boussinesq equations. Numerical calculations at values of the Rayleigh number up to 1.5 × 10⁵ indicate that the additional constraint decreases the upper bound on the Nusselt number; it appears that this upper bound differs by only a multiplicative factor from that calculated from solutions of the full equations of motion and may be a reasonable approximation for large Rayleigh numbers.     

Testing the anelastic nonhydrostatic model EULAG as a prospective dynamical core of a numerical weather prediction model Part I: Dry benchmarks

In this paper, a feasibility of anelastic approach for numerical weather prediction (NWP) is examined. The study concerns the anelastic nonhydrostatic model EULAG as a prospective candidate for the new dynamical core of a high-resolution NWP model. Such an application requires a series of benchmark tests to be performed. The study presents the results of dry idealized two-dimensional linear and non-linear tests. They include evolution of cold and warm density currents in neutrally stratified atmosphere, inertia-gravity waves in short and long channels, as well as mountain gravity waves for a set of different flow regimes. Detailed comparison of the results with the reference solutions, based mainly on the results of compressible models, indicates a high level of conformity for all of the experiments. It verifies the anelastic approach as strongly consistent with the compressible one for a broad class of atmospheric problems. It also corroborates the robustness of EULAG numerics, an essential requirement of dynamical core of NWP model.     

Asymptotic aspects of the Boussinesq approximation for gases and liquids

Abstract

The asymptotic formulation of the Boussinesq approximation relates the pressure of the fluid to a thermodynamical quantity involving the heat capacity c_Po . In this paper we examine the implications of such a scaling, in particular: (i) the singular degeneracy of the equation of state ρ = ρ* (1 ?α* (T?;T*)) of a liquid: this equation of state is valid only for small values of the coefficient α T*; (ii) in which manner the scaling introduces the Mach number of the flow as a small parameter e for a compressible fluid. The equations at order zero with respect to ? are the same equations for gases and for liquids only if the thermodynamics of the medium is described by using the Brunt-Väisälä frequency instead of the temperature.