A note on the free convection boundary layer on a vertical surface with prescribed heat flux at small Prandtl number

It is shown that for a particular case of the surface heat flux the equations for small Prandtl number have simple analytical solutions. These are presented and compared with numerical solutions of the general equations

Effects of constant electric fields on the buoyant stability of reaction fronts

The effects that applying constant electric fields have on the buoyant instability of reaction fronts propagating vertically in a Hele-Shaw cell are investigated for a range of electric field strengths and fluid parameters. The reaction produces a decrease in density across the front such that upwards propagating fronts are buoyantly unstable in the field-free situation. The reaction kinetics are modeled by cubic autocatalysis. A linear stability analysis reveals that a positive electric field increases the stability of a reaction front and can stabilize an otherwise unstable front. A negative field has the opposite effect, making the reaction front more unstable. Numerical simulations of the full nonlinear problem confirm these predictions and show the development of cellular fingers on unstable fronts. These simulations show that the electric field effects on the reaction within the front can alter the fluid density so as to give the possibility of destabilizing an otherwise stable downward propagating front

Dispersion curves in the diffusional instability of reaction fronts

A (linear) stability analysis of planar reaction fronts to transverse perturbations is considered for systems based on cubic autocatalysis and a model for the chlorite-tetrathionate reaction. Dispersion curves (plots of the growth rate sigma against a transverse wave-number k) are obtained. In both cases it is seen that there is a nonzero value D-0 of D (the ratio of the diffusion coefficients of autocatalyst and substrate) at which sigma(max), the maximum value of sigma for a given value of D, achieves its largest value, with sigma(max) being less for other values of D and becoming small as D decreases to zero. The existence of the optimum value D-0 for initiating a diffusional instability is confirmed, in the cubic autocatalysis case, by an asymptotic analysis for small wave numbers

Homogenization induced by chaotic mixing and diffusion in an oscillatory chemical reaction

A model for an imperfectly mixed batch reactor with the chlorine dioxide-iodine-malonic acid (CDIMA) reaction, with the mixing being modelled by chaotic advection, is considered. The reactor is assumed to be operating in oscillatory mode and the way in which an initial spatial perturbation becomes homogenized is examined. When the kinetics are such that the only stable homogeneous state is oscillatory then the perturbation is always entrained into these oscillations. The rate at which this occurs is relatively insensitive to the chemical effects, measured by the Damkohler number, and is comparable to the rate of homogenization of a passive contaminant. When both steady and oscillatory states are stable, spatially homogeneous states, two possibilities can occur. For the smaller Damkohler numbers, a localized perturbation at the steady state is homogenized within the background oscillations. For larger Damkohler numbers, regions of both oscillatory and steady behavior can co-exist for relatively long times before the system collapses to having the steady state everywhere. An interpretation of this behavior is provided by the one-dimensional Lagrangian filament model, which is analyzed in detail

Predicting magnetospheric dynamics with a coupled sun-to-Earth model: challenges and first results

Results from the first Sun-to-Earth coupled numerical model developed at the Center for Integrated Space Weather Modeling are presented. The model simulates physical processes occurring in space spanning from the corona of the Sun to the Earth's ionosphere, and it represents the first step toward creating a physics-based numerical tool for predicting space weather conditions in the near-Earth environment. Two 6- to 7-d intervals, representing different heliospheric conditions in terms of the three-dimensional configuration of the heliospheric current sheet, are chosen for simulations. These conditions lead to drastically different responses of the simulated magnetosphere-ionosphere system, emphasizing, on the one hand, challenges one encounters in building such forecasting tools, and on the other hand, emphasizing successes that can already be achieved even at this initial stage of Sun-to-Earth modeling

Flow-distributed spikes for Schnakenberg kinetics

This is the post-print version of the final published paper. The final publication is available at link.springer.com by following the link below. Copyright @ 2011 Springer-Verlag.We study a system of reaction–diffusion–convection equations which combine a reaction–diffusion system with Schnakenberg kinetics and the convective flow equations. It serves as a simple model for flow-distributed pattern formation. We show how the choice of boundary conditions and the size of the flow influence the positions of the emerging spiky patterns and give conditions when they are shifted to the right or to the left. Further, we analyze the shape and prove the stability of the spikes. This paper is the first providing a rigorous analysis of spiky patterns for reaction-diffusion systems coupled with convective flow. The importance of these results for biological applications, in particular the formation of left–right asymmetry in the mouse, is indicated.RGC of Hong Kon