Range of validity of weakly non-linear theory in Rayleigh-Bénard convection

In this paper we examine the equilibrium states of periodic finite amplitude flow in a horizontal channel with differential heating between the two rigid boundaries. The solutions to the Navier-Stokes equations are obtained by means of a perturbation method for evaluating the Landau coefficients and through a Newton-Raphson iterative method that results from the Fourier expansion of the solutions that bifurcate above the linear stability threshold of infini- tesimal disturbances. The results obtained from these two different methods of evaluating the convective flow are compared in the neighbourhood of the critical Rayleigh number. We find that for small Prandtl numbers the discrepancy of the two methods is noticeable.

The charge effect on the hindrance factors for diffusion and convection of a solute in pores II

The diffusion and convection of a solute suspended in a fluid across porous membranes are known to be reduced compared to those in a bulk solution, owing to the fluid mechanical interaction between the solute and the pore wall as well as steric restriction. If the solute and the pore wall are electrically charged, the electrostatic interaction between them could affect the hindrance to diffusion and convection. In this study, the transport of charged spherical solutes through charged circular cylindrical pores filled with an electrolyte solution containing small ions was studied numerically by using a fluid mechanical and electrostatic model. Based on a mean field theory, the electrostatic interaction energy between the solute and the pore wall was estimated from the Poisson-Boltzmann equation, and the charge effect on the solute transport was examined for the solute and pore wall of like charge. The results were compared with those obtained from the linearized form of the Poisson-Boltzmann equation, i.e.the Debye-Hückel equation. © 2012 The Japan Society of Fluid Mechanics and IOP Publishing Ltd.

Charge effects on hindrance factors for diffusion and convection of solute in pores I

The transport of a spherical solute through a long circular cylindrical pore filled with an electrolyte solution is studied numerically, in the presence of constant surface charge on the solute and the pore wall. Fluid dynamic analyses were carried out to calculate the flow field around the solute in the pore to evaluate the drag coefficients exerted on the solute. Electrical potentials around the solute in the electrolyte solution were computed based on a mean-field theory to provide the interaction energy between the charged solute and the pore wall. Combining the results of the fluid dynamic and electrostatic analyses, we estimated the rate of the diffusive and convective transport of the solute across the pore. Although the present estimates of the drag coefficients on the solute suggest more than 10% difference from existing studies, depending on the radius ratio of the solute relative to the pore and the radial position of the solute center in the pore, this difference leads to a minor effect on the hindrance factors. It was found that even at rather large ion concentrations, the repulsive electrostatic interaction between the charged solute and the pore wall of like charge could significantly reduce the transport rate of the solute.

Effect of aspect ratio on thermal convection in open and closed systems

Internally heated fluids are found across the nuclear fuel cycle. In certain situations the motion of the fluid is driven by the decay heat (i.e. corium melt pools in severe accidents, the shutdown of liquid metal reactors, molten salt and the passive control of light water reactors) as well as normal operation (i.e. intermediate waste storage and generation IV reactor designs). This can in the long-term affect reactor vessel integrity or lead to localized hot spots and accumulation of solid wastes that may prompt local increases in activity. Two approaches to the modeling of internally heated convection are presented here. These are based on numerical analysis using codes developed in-house and simulations using widely available computational fluid dynamics solvers. Open and closed fluid layers at around the transition between conduction and convection of various aspect ratios are considered. We determine optimum domain aspect ratio (1:7:7 up to 1:24:24 for open systems and 5:5:1, 1:10:10 and 1:20:20 for closed systems), mesh resolutions and turbulence models required to accurately and efficiently capture the convection structures that evolve when perturbing the conductive state of the fluid layer. Note that the open and closed fluid layers we study here are bounded by a conducting surface over an insulating surface. Conclusions will be drawn on the influence of the periodic boundary conditions on the flow patterns observed. We have also examined the stability of the nonlinear solutions that we found with the aim of identifying the bifurcation sequence of these solutions en route to turbulence.

Asymptotic theory of wall-attached convection in a horizontal fluid layer with a vertical magnetic field

A horizontal fluid layer heated from below in the presence of a vertical magnetic field is considered. A simple asymptotic analysis is presented which demonstrates that a convection mode attached to the side walls of the layer sets in at Rayleigh numbers much below those required for the onset of convection in the bulk of the layer. The analysis complements an earlier analysis by Houchens [J. Fluid Mech. 469, 189 (2002)] which derived expressions for the critical Rayleigh number for the onset of convection in a vertical cylinder with an axial magnetic field in the cases of two aspect ratios. © 2008 American Institute of Physics.

Bistable attractors in a model of convection-driven spherical dynamos

The range of existence and the properties of two essentially different chaotic attractors found in a model of nonlinear convection-driven dynamos in rotating spherical shells are investigated. A hysteretic transition between these attractors is established as a function of the rotation parameter t. The width of the basins of attraction is also estimated. © 2012 The Royal Swedish Academy of Sciences.

Bistability and hysteresis of dipolar dynamos generated by turbulent convection in rotating spherical shells

Bistability and hysteresis of magnetohydrodynamic dipolar dynamos generated by turbulent convection in rotating spherical fluid shells is demonstrated. Hysteresis appears as a transition between two distinct regimes of dipolar dynamos with rather different properties including a pronounced difference in the amplitude of the axisymmetric poloidal field component and in the form of the differential rotation. The bistability occurs from the onset of dynamo action up to about 9 times the critical value of the Rayleigh number for onset of convection and over a wide range of values of the ordinary and the magnetic Prandtl numbers including the value unity. Copyright © EPLA, 2009.

Convection induced by instabilities in the presence of a transverse seepage

The transition of laterally heated flows in a vertical layer and in the presence of a streamwise pressure gradient is examined numerically for the case of different values Prandtl number. The stability analysis of the basic flow for the pure hydrodynamic case ( Pr = 0 ) was reported in [1]. We find that in the absence of transverse pumping the previously known critical parameters are recovered [2], while as the strength of the Poiseuille flow component is increased the convective motion is delayed considerably. Following the linear stability analysis for the vertical channel flow our attention is focused on a study of the finite am- plitude secondary travelling-wave (TW) solutions that develop from the perturbations of the transverse roll type imposed on the basic flow and temperature profiles. The linear stability of the secondary TWs against three-dimensional perturbations is also examined and it is shown that the bifurcating tertiary flows are phase-locked to the secondary TWs.