On the onset of turbulence in natural convection on inclined plates

The problem of determination of the turbulence onset in natural convection on heated inclined plates in an air environment has been experimentally revisited. The transition has been detected by using hot wire velocity measurements. The onset of turbulence has been considered to take place where velocity fluctuations (measured through turbulence intensity) start to grow. Experiments have shown that the onset depends not only on the Grashof number defined in terms of the temperature difference between the heated plate and the surrounding air. A correlation between dimensionless Grashof and Reynolds numbers has been obtained, fitting quite well the experimental data.

Reentrant and whirling hexagons in non-Boussinesq convection

We review recent computational results for hexagon patterns in non- Boussinesq convection. For sufficiently strong dependence of the fluid parameters on the temperature we find reentrance of steady hexagons, i.e. while near onset the hexagon patterns become unstable to rolls as usually, they become again stable in the strongly nonlinear regime. If the convection apparatus is rotated about a vertical axis the transition from hexagons to rolls is replaced by a Hopf bifurcation to whirling hexagons. For weak non-Boussinesq effects they display defect chaos of the type described by the two-dimensional (2D) complex Ginzburg-andau equation. For stronger non-Boussinesq effects the Hopf bifurcation becomes subcritical and localized bursting of the whirling amplitude is found. In this regime the cou- pling of the whirling amplitude to (small) deformations of the hexagon lattice becomes important. For yet stronger non-Boussinesq effects this coupling breaks up the hexagon lattice and strongly disordered states characterized by whirling and lattice defects are obtained.

Hydrothermal waves and corotating rolls in laterally heated convection in simple liquids.

The stability of a liquid layer with an undeformable interface open to the atmo- sphere, subjected to a horizontal temperature gradient, is theoretically analysed. Buoyancy and surface tension forces give rise to a basic flow for any temperature dif- ference applied on the system. Depending on the liquid depth, this basic flow is desta- bilised either by an oscillatory instability, giving rise to the so-called hydrothermal waves, or by a stationary instability leading to corotating rolls. Oscillatory perturba- tions are driven by the basic flow and therefore one must distinguish between convec- tive and absolute thresholds. The instability mechanisms as well as the di¿erent re- gimes observed in experiments are discussed. The calculations are performed for a fluid used in recent experiments, namely silicone oil of 0.65 cSt ðPr 1?4 10Þ. In partic- ular, it is shown that two branches of absolute instability exist, which may be related to the two types of hydrothermal waves observed experimentally

Hexagonal patterns in a model for rotating convection

We study a model equation that mimics convection under rotation in a fluid with temperature- dependent properties (non-Boussinesq (NB)), high Prandtl number and idealized boundary conditions. It is based on a model equation proposed by Segel [1965] by adding rotation terms that lead to a Kuppers-Lortz instability [Kuppers & Lortz, 1969] and can develop into oscillating hexagons. We perform a weakly nonlinear analysis to find out explicitly the coefficients in the amplitude equation as functions of the rotation rate. These equations describe hexagons and os- cillating hexagons quite well, and include the Busse?Heikes (BH) model [Busse & Heikes, 1980] as a particular case. The sideband instabilities as well as short wavelength instabilities of such hexagonal patterns are discussed and the threshold for oscillating hexagons is determined.

Hexagons and spiral defect chaos in non-Boussinesq convection at low Prandtl numbers

We study the stability and dynamics of non-Boussinesq convection in pure gases ?CO2 and SF6? with Prandtl numbers near Pr? 1 and in a H2-Xe mixture with Pr= 0.17. Focusing on the strongly nonlinear regime we employ Galerkin stability analyses and direct numerical simulations of the Navier-Stokes equations. For Pr ? 1 and intermediate non-Boussinesq effects we find reentrance of stable hexagons as the Rayleigh number is increased. For stronger non-Boussinesq effects the usual, transverse side-band instability is superseded by a longitudinal side-band instability. Moreover, the hexagons do not exhibit any amplitude instability to rolls. Seemingly, this result contradicts the experimentally observed transition from hexagons to rolls. We resolve this discrepancy by including the effect of the lateral walls. Non-Boussinesq effects modify the spiral defect chaos observed for larger Rayleigh numbers. For convection in SF6 we find that non-Boussinesq effects strongly increase the number of small, compact convection cells and with it enhance the cellular character of the patterns. In H2-Xe, closer to threshold, we find instead an enhanced tendency toward roll-like structures. In both cases the number of spirals and of targetlike components is reduced. We quantify these effects using recently developed diagnostics of the geometric properties of the patterns.

Defect chaos and bursts: Hexagonal rotating convection and the complex Ginzburg-Landau equation

We employ numerical computations of the full Navier-Stokes equations to investigate non-Boussinesq convection in a rotating system using water as the working fluid. We identify two regimes. For weak non- Boussinesq effects the Hopf bifurcation from steady to oscillating (whirling) hexagons is supercritical and typical states exhibit defect chaos that is systematically described by the cubic complex Ginzburg-Landau equation. For stronger non-Boussinesq effects the Hopf bifurcation becomes subcritical and the oscil- lations exhibit localized chaotic bursting, which is modeled by a quintic complex Ginzburg-Landau equation.

Reentrant hexagons in non-Boussinesq convection.

While non-Boussinesq hexagonal convection patterns are known to be stable close to threshold (i.e. for Rayleigh numbers R ? Rc ), it has often been assumed that they are always unstable to rolls for slightly higher Rayleigh numbers. Using the incompressible Navier?Stokes equations for parameters corresponding to water as the working fluid, we perform full numerical stability analyses of hexagons in the strongly nonlinear regime ( ? (R ? Rc )/Rc = O(1)). We find ?re-entrant? behaviour of the hexagons, i.e. as is increased they can lose and regain stability. This can occur for values of as low as = 0.2. We identify two factors contributing to the re-entrance: (i) far above threshold there exists a hexagon attractor even in Boussinesq convection as has been shown recently and (ii) the non-Boussinesq effects increase with . Using direct simulations for circular containers we show that the re-entrant hexagons can prevail even for sidewall conditions that favour convection in the form of competing stable rolls. For sufficiently strong non-Boussinesq effects hexagons even become stable over the whole -range considered, 0 6 6 1.5.

A subgrid viscosity Lagrance-Galerkin method for convection-diffusion problems

We present and analyze a subgrid viscosity Lagrange-Galerk in method that combines the subgrid eddy viscosity method proposed in W. Layton, A connection between subgrid scale eddy viscosity and mixed methods. Appl. Math. Comp., 133: 14 7-157, 2002, and a conventional Lagrange-Galerkin method in the framework of P1⊕ cubic bubble finite elements. This results in an efficient and easy to implement stabilized method for convection dominated convection diffusion reaction problems. Numerical experiments support the numerical analysis results and show that the new method is more accurate than the conventional Lagrange-Galerkin one.

Electric current of an electrified jet issuing from a long metallic tube.

This paper presents an analysis of the transport of electric current in a jet of an electrically conducting liquid discharging from a metallic tube into a gas or a vacuum, and subject to an electric field due to a high voltage applied between the tube and a far electrode. The flow, the surface charge and the electric field are computed in the current transfer region of the jet, where conduction current in the liquid becomes surface current due to the convection of electric charge accumulated at its surface. The electric current computed as a function of the flow rate of the liquid injected through the tube increases first as the square root of this flow rate, levels to a nearly constant value when the flow rate is increased and finally sets to a linear increase when the flow rate is further increased. The current increases linearly with the applied voltage at small and moderate values of this variable, and faster than linearly at high voltages. The characteristic length and structure of the current transfer region are determined. Order-of-magnitude estimates for jets which are only weakly stretched by the electric stresses are worked out that qualitatively account for some of the numerical results.

Time Evolution of the Fractal Dimension in Turbulent Plumes

Turbulent mixing is a very important issue in the study of geophysical phenomena because most fluxes arising in geophysics fluids are turbulent. We study turbulent mixing due to convection using a laboratory experimental model with two miscible fluids of different density with an initial top heavy density distribution. The fluids that form the initial unstable stratification are miscible and the turbulence will produce molecular mixing. The denser fluid comes into the lighter fluid layer and it generates several forced plumes which are gravitationally unstable. As the turbulent plumes develop, the denser fluid comes into contact with the lighter fluid layer and the mixing process grows. Their development is caused by the lateral interaction between these plumes at the complex fractal surface between the dense and light fluids

The internal structure of hydrogen-air difussion flames

The Internal Structure of Hydrogen-Air Diffusion Flames. Tho purpose of this paper is to study finite rate chemistry effects in diffusion controlled hydrogenair flames undor conditions appearing in some cases in a supersonic combustor. Since for large reaction rates the flame is close to chemical equilibrium, the reaction takes place in a very thin region, so thata "singular perturbation "treatment" of the problem seems appropriate. It has been shown previously that, within the inner or reaction zone, convection effects may be neglocted, the temperature is constant across the flame, and tho mass fraction distributions are given by ordinary differential equations, whore tho only independent variable involved is tho coordinate normal to the flame surface. Tho solution of the outer problom, which is a pure mixing problem with the additional condition that fuol and oxidizer do not coexist in any zone, provides t h e following information: tho flame position, rates of fuel consumption, temperature, concentrators of species, fluid velocity outside of tho flame, and the boundary conditions required to solve the "inner problem." The main contribution of this paper consists in the introduction of a fairly complicated chemical kinetic scheme representing hydrogen-oxygen reaction. The nonlinear equations expressing the conservation of chemical species are approximately integrated by means of an integral method. It has boen found that, in the case considered of a near-equilibrium diffusion flame, tho role played by the dissociation-recombination reactions is purely marginal, and that somo of the second order "shuffling" reactions are close to equilibrium. The method shown here may be applied to compute the distanco from the injector corresponding to a given separation from equilibrium, say ten to twenty percent. For the casos whore this length is a small fraction of the combustion zone length, the equilibrium treatment describes properly tho flame behavior.

A regular perturbation approach to surface tension driven flows

Surface tension induced convection in a liquid bridge held between two parallel, coaxial, solid disks is considered. The surface tension gradient is produced by a small temperature gradient parallel Co the undisturbed surface. The study is performed by using a mathematical regular perturbation approach based on a small parameter, e, which measures the deviation of the imposed temperature field from its mean value. The first order velocity field is given by a Stokes-type problem (viscous terms are dominant) with relatively simple boundary conditions. The first order temperature field is that imposed from the end disks on a liquid bridge immersed in a non-conductive fluid. Radiative effects are supposed to be negligible. The second order temperature field, which accounts for convective effects, is split into three components, one due to the bulk motion, and the other two to the distortion of the free surface. The relative importance of these components in terms of the heat transfer to or from the end disks is assessed

Some physical and chemical proceddes in fluids

The ability to reproduce reduced gravity conditions for long periods is one of the reasons why the orbiting laboratory is so attractive. In this paper several fluid dynamics problem areas are reviewed in which zero-gravity conditions are of great importance. Although emphasis is placed on space processing, there are some older problems also in which gravity masks the phenomcna, impeding a reasonably simple approach to the solution. Three problems are considered: Thermal convection under reduced gravity. The dumping effect ofsurface gravity waves at the outset of convection induced by surface tractions is discussed in particular. The existence of convection is of concern for some satellite thermal control techniques presently used, and for most of the proposed manufacturing processes. Whereas convection should be normally avoided, problems related to the containerless stirring ofa melt constitute an exception. Secondly, gravity and chemical reactions. Although chemical reactions are independent of gravity because of the small mass of the molecules and atoms involved, in many cases the reaction rate dcpends on the arrival of the species to the reaction zone. When the arrival process is buoyancy-controlled, the net specd of the reaction will be affected by the gravity. Thirdly, two-phase flows under reduced gravity provkle interesting problems from boiling heat transfer to degasslng of melts. This part of the paper deals only with the measurement of sound veiocity in a liquid containing bubbles. It is suggested that such measurements should be mude under reduced gravity to provide reliable residís.

Self-similar motion of laser half-space plasmas. II. Thermal waves and intermediate regimes

The one-dimensional self-similar motion of an initially cold, half-space plasma of electron density n,produced by the (anomalous) absorption of a laser pulse of irradiation

convection is negligible) moving into the undisturbed plasma, from a much thinner isothermal flow expanding into the vacuum. For l€~4'3, a qualitative discussion of how plasma behavior changes with a, is given.

Decomposition driven interface evolution for layers of binary mixtures: I. Model dirivation and base states

Algebraic topology (homology) is used to analyze the state of spiral defect chaos in both laboratory experiments and numerical simulations of Rayleigh-Bénard convection. The analysis reveals topological asymmetries that arise when non-Boussinesq effects are present. The asymmetries are found in different flow fields in the simulations and are robust to substantial alterations to flow visualization conditions in the experiment. However, the asymmetries are not observable using conventional statistical measures. These results suggest homology may provide a new and general approach for connecting spatiotemporal observations of chaotic or turbulent patterns to theoretical models.