Effect of the combined thermal and mass concentration gradients on the stability of a chemically reacting fluid in a porous medium

Summary The onset of instability, due to the combined effects of thermal and mass concentration gradients, is investigated in the hydrodynamic stability regime. It is found that the critical Rayleigh number, which determines the onset of instability, increases as the chemical reaction rate constant increases hyperbolically over a wide range of values at both moderate and high permeabilities. In addition, the instability grows with increase in porosity. Previous results show that the critical Rayleigh number rises linearly when only the mass concentration gradient is considered.     

Interaction and dynamics of defects in convective roll patterns of anisotropic fluids

We present an overview of the dynamics and interaction of defects in roll patterns of electroconvection in nematic liquid crystals (EHC). For the decay of an Eckhaus-unstable pattern we distinguish three regimes, depending on the width of the system perpendicular to the wavenumber mismatch. Motivated by recent experiments, we examine the annihilation process of defects in patterns with wavenumber near to band center, where the motion of the defects is dominated by the interaction. The comparison with the experiments shows that this process can be described even quantitatively within the framework of Ginzburg-Landau theory.     

Three-dimensional absolute and convective instabilities in mixed convection of a viscoelastic fluid through a porous medium

By using the mathematical formalism of absolute and convective instabilities we study the nature of unstable three-dimensional disturbances of viscoelastic flow convection in a porous medium with horizontal through-flow and vertical temperature gradient. Temporal stability analysis reveals that among three-dimensional (3D) modes the pure down-stream transverse rolls are favored for the onset of convection. In addition, by considering a spatiotemporal stability approach we found that all unstable 3D modes are convectively unstable except the transverse rolls which may experience a transition to absolute instability. The combined influence of through-flow and elastic parameters on the absolute instability threshold, wave number and frequency is then determined, and results are compared to those of a Newtonian fluid.     

Condition for convective instability of dark solitons

Simple derivation of the condition for the transition point from absolute instability of plane dark solitons to their convective instability is suggested. It is shown that unstable wave packet expands with velocity equal to the minimal group velocity of the disturbance waves propagating along a dark soliton. The growth rate of the length of dark solitons generated by the flow of Bose-Einstein condensate past an obstacle is estimated. Analytical theory is confirmed by the results of numerical simulations.     

Phenomenological equations for multicomponent fluids

Ahydrodynamic equation of motion for each component of a multicomponent fluid is derived on the basis of nonequilibrium thermodynamics. Special care has been directed to the choice of state variables. In some limiting cases, this equation leads to customary phenomenological equations, such as the equation for diffusion and the Navier-Stokes equation. The viscosity is a consequence of nonlocal coupling of forces and fluxes. The reciprocity between the linear coefficients is examined closely.     

Kinetic equations for dense fluids

It is demonstrated that the kinetic equation of Davis's effective potential theory follows directly from the application of well-defined approximations to the three-body correlations involved in the second equation of the BBGKY hierarchy. The same, simple mathematical techniques involved in this demonstration are used to derive two other kinetic equations, one of which is a generalization to high densities of the Boltzmann equation. In order to facilitate its application to the calculation of the van Hove and other correlation functions, the kinetic equation of the effective potential theory is Fourier-Laplace transformed: explicit formulae are given for the matrix elements of all operators that occur in this equation.     

Variable amplitude equations for one-dimensional scattering

Nonlinear first-order equations, similar to Calogero's equations, are derived for the forward and backward one-dimensional scattering amplitudes. In particular, the even potential case yields two uncoupled equations for the even and odd parity phase shifts. The present approach provides a fast and accurate means for the numerical solution of one-dimensional scattering problems. It also has many analytic merits, some of which are discussed. The connection between one-dimensional and three-dimensional high-energy scattering is reviewed. It is demonstrated that in the one-dimensional case, a slightly modified WKB wavefunction provides an excellent approximation to the exact wavefunction in the shortwave limit. In this limit, additivity of phase shifts for nonoverlapping static potentials is satisfied.     

Stability analysis on a set of calcium-regulated viscoelastic equations

In recent years some progress has been made in modelling pattern formation and morphogenesis in biological systems in terms of calcium ion regulation of the viscoelastic properties of the cellular cortex. In this paper, linear stability analysis is used on a set of calcium-regulated viscoelastic equations derived by Goodwin and Trainor [5] for the 3-dimensional medium appropriate to regeneration phenomena in the single celled alga Acetabularia mediterranea. The nature of the instabilities is discussed and it is shown how complex patterns arise naturally from the cross-terms linking viscoelastic strain to calcium concentration and concentration gradients.     

Derivation of the amplitude equation for a hydromagnetic convective problem

The amplitude equation for a convective system under a vertical magnetic field is derived. The coefficients in these equations have been numerically calculated for infinite Prandtl number fluids and for boundary conditions both free and rigid top and bottom. The results confirm that, for realistic parameter values only stationary convection can be present and the pattern is made by convective rolls.     

Derivation of the amplitude equation for a hydromagnetic convective problem

The amplitude equation for a convective system under a vertical magnetic field is derived. The coefficients in these equations have been numerically calculated for infinite Prandtl number fluids and for boundary conditions both free and rigid top and bottom. The results confirm that, for realistic parameter values only stationary convection can be present and the pattern is made by convective rolls.     

Derivation of the amplitude equation for a hydromagnetic convective problem

The amplitude equation for a convective system under a vertical magnetic field is derived. The coefficients in these equations have been numerically calculated for infinite Prandtl number fluids and for boundary conditions both free and rigid top and bottom. The results confirm that, for realistic parameter values only stationary convection can be present and the pattern is made by convective rolls.     

On wave equations for viscous fluids

This paper deals with wave equations describing small-amplitude disturbances in horizontally stratified, continuously varying, viscous fluids; gradients of the static pressure and of the coefficient of viscosity are neglected. A set of equations in first-order matrix form, which describes coupled longitudinal and transverse disturbances, is treated by the methods ofClemmow andHeading and ofHeading.The work of this paper could be extended in a number of ways; for example, the effect of a gravitational field could be included, and the coefficient of viscosity could be allowed to vary with position.     

High-order Absorbing Boundary Conditions for anisotropic and convective wave equations

High-order Absorbing Boundary Conditions (ABCs), applied on a rectangular artificial computational boundary that truncates an unbounded domain, are constructed for a general two-dimensional linear scalar time-dependent wave equation which represents acoustic wave propagation in anisotropic and subsonically convective media. They are extensions of the construction of Hagstrom, Givoli and Warburton for the isotropic stationary case. These ABCs are local, and involve only low-order derivatives owing to the use of auxiliary variables on the artificial boundary. The accuracy and well-posedness of these ABCs is analyzed. Special attention is given to the issue of mismatch between the directions of phase and group velocities, which is a potential source of concern. Numerical examples for the anisotropic case are presented, using a finite element scheme.     

Bubble oscillation and inertial cavitation in viscoelastic fluids

Non-linear acoustic oscillations of gas bubbles immersed in viscoelastic fluids are theoretically studied. The problem is formulated by considering a constitutive equation of differential type with an interpolated time derivative. With the aid of this rheological model, fluid elasticity, shear thinning viscosity and extensional viscosity effects may be taken into account. Bubble radius evolution in time is analyzed and it is found that the amplitude of the bubble oscillations grows drastically as the Deborah number (the ratio between the relaxation time of the fluid and the characteristic time of the flow) increases, so that, even for moderate values of the external pressure amplitude, the behavior may become chaotic. The quantitative influence of the rheological fluid properties on the pressure thresholds for inertial cavitation is investigated. Pressure thresholds values in terms of the Deborah number for systems of interest in ultrasonic biomedical applications, are provided. It is found that these critical pressure amplitudes are clearly reduced as the Deborah number is increased.     

Subcritical finite-amplitude solutions for plane Couette flow of viscoelastic fluids

Plane Couette flow of viscoelastic fluids is shown to exhibit a purely elastic subcritical instability at a very small-Reynolds number in spite of being linearly stable. The mechanism of this instability is proposed and the nonlinear stability analysis of plane Couette flow of the Upper-Convected Maxwell fluid is presented. Above a critical Weissenberg number, a small finite-size perturbation is sufficient to create a secondary flow, and the threshold value for the amplitude of the perturbation decreases as the Weissenberg number increases. The results suggest a scenario for weakly turbulent viscoelastic flow which is similar to the one for Newtonian fluids as a function of Reynolds number.     

A Reynolds stress model for turbulent flows of viscoelastic fluids

A second-order closure is developed for predicting turbulent flows of viscoelastic fluids described by a modified generalised Newtonian fluid model incorporating a nonlinear viscosity that depends on a strain-hardening Trouton ratio as a means to handle some of the effects of viscoelasticity upon turbulent flows. Its performance is assessed by comparing its predictions for fully developed turbulent pipe flow with experimental data for four different dilute polymeric solutions and also with two sets of direct numerical simulation data for fluids theoretically described by the finitely extensible nonlinear elastic – Peterlin model. The model is based on a Newtonian Reynolds stress closure to predict Newtonian fluid flows, which incorporates low Reynolds number damping functions to properly deal with wall effects and to provide the capability to handle fluid viscoelasticity more effectively. This new turbulence model was able to capture well the drag reduction of various viscoelastic fluids over a wide range of Reynolds numbers and performed better than previously developed models for the same type of constitutive equation, even if the streamwise and wall-normal turbulence intensities were underpredicted.     

Generalized hydrodynamic equations and correlation functions for thermoviscoelastic fluids

For very viscous liquids a phenomenological theory of thermoviscoelasticity is formulated, in which the retarded reaction of thermal variables, which arises from structural relaxation, is taken into account. The theory describes the effect of the slowing down of the structural relaxation near a glass transition on the fluctuation spectra of density and entropy; in particular, the intensity of the slow relaxational component of the fluctuation spectra, which is frozen in the glass below the glass transition, is derived. Conditions for positive energy dissipation and symmetry relations are obtained in the framework of thermodynamic relaxation theory, and the memory functions occurring in the Mori-Zwanzig projection operator formalism are calculated.Dedicated to K. Dransfeld on the occasion of his 60th birthday     

Droplet detachment and satellite bead formation in viscoelastic fluids

The presence of a very small amount of high molecular weight polymer significantly delays the pinch-off singularity of a drop of water falling from a faucet and leads to the formation of a long-lived cylindrical filament. In this Letter, we present experiments, numerical simulations, and theory which examines the pinch-off process in the presence of polymers. The numerical simulations are found to be in good agreement with experiment. As a test case, we establish the conditions under which a small bead remains on the filament; we find that the presence of a bead is due to the asymmetry induced by the self-similar pinch off of the droplet.     

Exact equations of state for one-dimensional chain fluids

Using the isothermal-isobaric ensemble, exact equations of state are derived for three classical models of one-dimensional chain fluids. Each chain molecule is modeled by a series of linked sites which interact through nearest-neighbor bond potentials. In two of the models, the intramolecular bonds are modeled by infinitely deep square-well potentials, while in the third, the bonds are modeled by a harmonic potential. Intermolecular interactions are modeled by a hard-rod potential. Numerical results are presented for dimer and 8-mer fluids which illustrate the influence of chain length, well width and spring constant on the compressibility factor. The effect of adding an infinitely weak, infinitely long-ranged attractive interaction between the sites is also considered. The attractive tail induces a first-order phase transition of the gas-liquid type in all of the chain models. For certain values of the model parameters, however, two of the models show evidence of a second gas-liquid type transition, which appears to be associated with chain collapse.