Taylor–Couette problem and related topics

The O(2)×S¹-equivariant degree is applied to study the Hopf bifurcation for the Taylor–Couette problem. Isotypical crossing numbers and equivariant bifurcation invariants are evaluated and are applied to classify the equivariant Hopf bifurcation for the Taylor–Couette flow.     

The Cauchy–Neumann problem for a multi‐dimensional isentropic hydrodynamic model for semiconductors

We investigate a multi‐dimensional isentropic hydrodynamic (Euler–Poisson) model for semiconductors, where the energy equation is replaced by the pressure–density relation p(n) . We establish the global existence of smooth solutions for the Cauchy–Neumann problem with small perturbed initial data and homogeneous Neumann boundary conditions. We show that, as t→+∞, the solutions converge to the non‐constant stationary solutions of the corresponding drift–diffusion equations. Moreover, we also investigate the existence and uniqueness of the stationary solutions for the corresponding drift–diffusion equations. Copyright © 2005 John Wiley & Sons, Ltd.     

Couette and Poiseuille flows of an Oldroyd 6-constant fluid with magnetic field

In this paper, we develop a set of differential equations describing the steady flow of an Oldroyd 6-constant magnetohydrodynamic fluid. The fluid is electrically conducting in the presence of a uniform transverse magnetic field. The developed non-linear differential equation takes into account the effect of the material constants and the applied magnetic field. We presented the solution for three types of steady flows, namely,

(i)

Couette flow

(ii)

Poiseuille flow and

(iii)

generalized Couette flow.

Homotopy analysis method (HAM) is used to solve the non-linear differential equation analytically. It is found from the present analysis that for steady flow the obtained solutions are strongly dependent on the material constants (non-Newtonian parameters) which is different from the model of Oldroyd 3-constant fluid. Numerical solutions are also given and compared with the solutions by HAM.     

Ridged Domains,Embedding Theorems and Poincaré Inequalities

We study the embeddings E : W(X(Ω), Y(Ω)) ↪ Z(Ω), where X(Ω), Y(Ω) and Z(Ω) are rearrangement–invariant Banach function spaces (BFS) defined on a generalized ridged domain Ω, and W denotes a first–order Sobolev–type space. We obtain two–sided estimates for the measure of non–compactness of E when Z(Ω) = X(Ω) and, in turn, necessary and sufficient conditions for a Poincaré–type inequality to be valid and also for E to be compact. The results are used to analyse the example of a trumpet–shaped domain Ω in Lorentz spaces. We consider the problem of determining the range of possible target spaces Z(Ω), in which case we prove that the problem is equivalent to an analogue on the generalized ridge Γ of Ω. The range of target spaces Z(Ω) is determined amongst a scale of (weighted) Lebesgue spaces for “rooms and passages” and trumpet–shaped domains.     

Comment on: “Fundamental flows with nonlinear slip conditions: exact solutions”, by R. Ellahi,T. Hayat,F. M. Mahomed and A. Zeeshan,Z. Angew. Math. Phys. 61 (2010) 877–888

In their article (Fundamental flows with nonlinear slip conditions: exact solutions, R. Ellahi, T. Hayat, F. M. Mahomed and A. Zeeshan, Z. Angew. Math. Phys. 61 (2010) 877–888.), the authors considered three simple cases of the steady flow of a third grade fluid between parallel plates with slip conditions; namely, Couette flow, Poiseuille flow, and generalized Couette flow. They obtained exact solutions, which were utilized in a way that did not lead to useful results. Their conclusion that the Couette flow cannot be obtained from the generalized Couette flow, by dropping the pressure gradient, is incorrect. Meaningful results based on their solution are herein presented.     

New exact solutions for magnetohydrodynamic flows of an Oldroyd-B fluid

This paper presents the new exact analytical solutions for magnetohydrodynamic (MHD) flows of an Oldroyd-B fluid. The explicit expressions for the velocity field and the associated tangential stress are established by using the Laplace transform method. Three characteristic examples: (i) flow due to impulsive motion of plate, (ii) flow due to uniformly accelerated plate, and (iii) flow due to non-uniformly accelerated plate are considered. The solutions for the hydrodynamic flows are special cases of the presented solutions. Moreover, the similar solutions corresponding to Maxwell and Newtonian fluids in the presence as well as absence of a magnetic field appear as the limiting cases of our solutions. The influences of the exerted magnetic field on the flow are also graphically presented and discussed. In particular, graphical results for the Oldroyd-B fluid are compared with those of a Newtonian fluid.     

Compatible metrics and the diagonalizability of nonlocally bi-Hamiltonian systems of hydrodynamic type

We study bi-Hamiltonian systems of hydrodynamic type with nonsingular (semisimple) nonlocal bi-Hamiltonian structures. We prove that all such systems of hydrodynamic type are diagonalizable and that the metrics of the bi-Hamiltonian structure completely determine the complete set of Riemann invariants constructed for any such system. Moreover, we prove that for an arbitrary nonsingular (semisimple) nonlocally bi-Hamiltonian system of hydrodynamic type, there exist local coordinates (Riemann invariants) such that all matrix differential-geometric objects related to this system, namely, the matrix (affinor) V_jⁱ(u) of this system of hydrodynamic type, the metrics g ₁ ^ij(u) and g ₂ ^ij(u), the affinor υ_jⁱ(u) = g ₁ ^is(u)g _2,sj(u), and also the affinors (w _1,n)_jⁱ(u) and (w _2,n)_jⁱ(u) of the nonsingular nonlocal bi-Hamiltonian structure of this system, are diagonal in these special “diagonalizing” local coordinates (Riemann invariants of the system). The proof is a natural corollary of the general results of our previously developed theories of compatible metrics and of nonlocal bi-Hamiltonian structures; we briefly review the necessary notions and results in those two theories.     

Some rigid Steiner 5‐designs

Hitherto, all known non‐trivial Steiner systems S(5, k, v) have, as a group of automorphisms, either PSL(2, v−1) or PGL(2, (v−2)/2) × C₂. In this article, systems S(5, 6, 72), S(5, 6, 84) and S(5, 6, 108) are constructed that have only the trivial automorphism group. © 2010 Wiley Periodicals, Inc. J Combin Designs 18:392–400, 2010     

Linear stability of radially-heated circular Couette flow with simulated radial gravity

The stability of circular Couette flow between vertical concentric cylinders in the presence of a radial temperature gradient is considered with an effective “radial gravity.” In addition to terrestrial buoyancy − ρg e _z we include the term − ρg _m f(r)e _r where g _m f(r) is the effective gravitational acceleration directed radially inward across the gap. Physically, this body force arises in experiments using ferrofluid in the annular gap of a Taylor–Couette cell whose inner cylinder surrounds a vertical stack of equally spaced disk magnets. The radial dependence f(r) of this force is proportional to the modified Bessel function K ₁(κr), where 2π/κ is the spatial period of the magnetic stack and r is the radial coordinate. Linear stability calculations made to compare with conditions reported by Ali and Weidman (J. Fluid Mech., 220, 1990) show strong destabilization effects, measured by the onset Rayleigh number R, when the inner wall is warmer, and strong stabilization effects when the outer wall is warmer, with increasing values of the dimensionless radial gravity γ = g _m /g. Further calculations presented for the geometry and fluid properties of a terrestrial laboratory experiment reveal a hitherto unappreciated structure of the stability problem for differentially-heated cylinders: multiple wavenumber minima exist in the marginal stability curves. Transitions in global minima among these curves give rise to a competition between differing instabilities of the same spiral mode number, but widely separated axial wavenumbers.     

The hydrodynamic limit for the reaction diffusion equation — An approach in terms of the GPV method

We study the hydrodynamic limit of the reaction diffusion process by means of the GPV technique (Guoet al. ⁽⁴⁾). To this end, we first derivea priori bounds on the moments of the occupation numbers using the local central limit theorem and results of stochastic analysis. The result of De Masi and Presutti⁽²⁾ for the hydrodynamic limit of the reaction diffusion process is generalized here.     

Linear stability of radially-heated circular Couette flow with simulated radial gravity

The stability of circular Couette flow between vertical concentric cylinders in the presence of a radial temperature gradient is considered with an effective “radial gravity.” In addition to terrestrial buoyancy − ρg e _z we include the term − ρg _m f(r)e _r where g _m f(r) is the effective gravitational acceleration directed radially inward across the gap. Physically, this body force arises in experiments using ferrofluid in the annular gap of a Taylor–Couette cell whose inner cylinder surrounds a vertical stack of equally spaced disk magnets. The radial dependence f(r) of this force is proportional to the modified Bessel function K ₁(κr), where 2π/κ is the spatial period of the magnetic stack and r is the radial coordinate. Linear stability calculations made to compare with conditions reported by Ali and Weidman (J. Fluid Mech., 220, 1990) show strong destabilization effects, measured by the onset Rayleigh number R, when the inner wall is warmer, and strong stabilization effects when the outer wall is warmer, with increasing values of the dimensionless radial gravity γ = g _m /g. Further calculations presented for the geometry and fluid properties of a terrestrial laboratory experiment reveal a hitherto unappreciated structure of the stability problem for differentially-heated cylinders: multiple wavenumber minima exist in the marginal stability curves. Transitions in global minima among these curves give rise to a competition between differing instabilities of the same spiral mode number, but widely separated axial wavenumbers.     

Pressure in Classical Statistical Mechanics and Interacting Brownian Particles in Multi-dimensions

. Let P(u) denote the pressure at the density u defined in the Gibbs statistical mechanics determined by a 2 body potential U (_qi - _qj). The function U(x) is supposed rotationally invariant and of finite range but may be unbounded about the origin. We establish a representation of P(u) by means of the law of large numbers for the virial ?_i,j q_i ·? U(q_i-q_j)\sum_{i,j} q_i \cdot {\nabla} U(q_i-q_j), whether or not there occur phase transitions. This result on P(u) is motivated by a study of the hydrodynamic behavior of a system of a large number of interacting Brownian particles moving on a d-dimensional torus (d = 1, 2, ...) in which the interaction is given by binary potential forces of potential U. Employing our representation of P(u), we also show that in the hydrodynamic limit of such a system there arises a non linear evolution equation of the form u_t = ¹/₂ DP(u)u_t = {1\over2} \Delta P(u) under a certain hypothetical postulate concerning concentration of particles.     

Ritz–Galerkin method for solving a parabolic equation with non‐local and time‐dependent boundary conditions

The paper is devoted to the investigation of a parabolic partial differential equation with non‐local and time‐dependent boundary conditions arising from ductal carcinoma in situ model. Approximation solution of the present problem is implemented by the Ritz–Galerkin method, which is a first attempt at tackling parabolic equation with such non‐classical boundary conditions. In the process of dealing with the difficulty caused by integral term in non‐local boundary condition, we use a trick of introducing the transition function G(x,t) to convert non‐local boundary to another non‐classical boundary, which can be handled with the Ritz–Galerkin method. Illustrative examples are included to demonstrate the validity and applicability of the technique in this paper. Copyright © 2015 John Wiley & Sons, Ltd.     

Zero‐sum flows in designs

Let D be a t ‐ ( v, k , λ) design and let N _i (D) , for 1 ≤ i ≤ t , be the higher incidence matrix of D , a ( 0 , 1 )‐matrix of size , where b is the number of blocks of D . A zero‐sum flow of D is a nowhere‐zero real vector in the null space of N ₁ ( D ). A zero‐sum k‐flow of D is a zero‐sum flow with values in { 1 , …, ±( k ? 1 )}. In this article, we show that every non‐symmetric design admits an integral zero‐sum flow, and consequently we conjecture that every non‐symmetric design admits a zero‐sum 5‐flow. Similarly, the definition of zero‐sum flow can be extended to N _i ( D ), 1 ≤ i ≤ t . Let be the complete design. We conjecture that N _t ( D ) admits a zero‐sum 3‐flow and prove this conjecture for t = 2 . © 2011 Wiley Periodicals, Inc. J Combin Designs 19:355‐364, 2011     

Differential-geometric approach to the integrability of hydrodynamic chains: the Haantjes tensor

The integrability of an m-component system of hydrodynamic type, u _t = V(u)u _x , by the generalized hodograph method requires the diagonalizability of the m ×  m matrix V(u). This condition is known to be equivalent to the vanishing of the corresponding Haantjes tensor. We generalize this approach to hydrodynamic chains—infinite-component systems of hydrodynamic type for which the ∞ ×  ∞ matrix V(u) is ‘sufficiently sparse’. For such systems the Haantjes tensor is well-defined, and the calculation of its components involves finite summations only. We illustrate our approach by classifying broad classes of conservative and Hamiltonian hydrodynamic chains with the zero Haantjes tensor. We prove that the vanishing of the Haantjes tensor is a necessary condition for a hydrodynamic chain to possess an infinity of semi-Hamiltonian hydrodynamic reductions, thus providing an easy-to-verify necessary condition for the integrability.     

On global symmetric solutions to the relativistic Vlasov–Poisson equation in three space dimensions

A collisionless plasma is modelled by the Vlasov–Maxwell system. In the presence of very large velocities, relativistic corrections are meaningful. When magnetic effects are ignored this formally becomes the relativistic Vlasov–Poisson equation. The initial datum for the phase space density ƒ₀(x, v) is assumed to be sufficiently smooth, non‐negative and cylindrically symmetric. If the (two‐dimensional) angular momentum is bounded away from zero on the support of ƒ₀(x, v), it is shown that a smooth solution to the Cauchy problem exists for all times. Copyright © 2001 John Wiley & Sons, Ltd.     

Closure for the property of having a hamiltonian prism

We prove that a graph G of order n has a hamiltonian prism if and only if the graph Cl_4n/3–4/3(G) has a hamiltonian prism where Cl_4n/3–4/3(G) is the graph obtained from G by sequential adding edges between non‐adjacent vertices whose degree sum is at least 4n/3–4/3. We show that this cannot be improved to less than 4n/3–5. © 2006 Wiley Periodicals, Inc. J Graph Theory 54: 209–220, 2007     

A Distribution- Theoretic Approach to a Stability Problem of Orr

In a classic paper [1] of 1907, W. M'Farr Orr discovered, among other things, the “infinitesimal” instability of inviscid plane Couette flow. Surprisingly, although Orr's paper remains a standard reference in the field, later investigators [2, 3] have been able to call inviscid plane Couette flow stable without finding it necessary to controvert Orr's result. What has happened is that, at least in problems governed by linear (or linearized) equations with time-independent coefficients, the term “instability” has come to be identified with the presence of solutions exhibiting exponential time-growth. Orr found instability indeed: a class of solutions certain members of which grow in time by more than each preassigned factor. Unlike the exponential instabilities, however, Orr's solutions die away like 1/t after achieving their greatest growth. This ephemerality probably accounts for the discounting of Orr's result. Orr did not look into the general initial value problem. This is done in the sequel, with the result that the situation becomes clear. Under general disturbances, Couette flow turns out to be neither stable nor quasi-asymptotically stable*. The rate of growth depends on the smoothness of the initial data: classical solutions grow no faster than t, but sufficiently rough distribution-valued initial data leads to growth matching any power of t. Before presenting detailed results, we briefly review Orr's fundamental work on the problem.     

Thermodynamics / turbulence analogy modelling dissipating molecular gas flows for high Knudsen numbers

Modelling micro channel flows momentum and heat diffusion / convection are recent parameters modelling the molecule velocity distribution. Macroscopic models describe velocity and energy / enthalpie with integrals of mass increments. Using microscopic models motion and forces of a molecular flow have to be computed by models of physical properties, whose are described by statistical power moments of the molecule velocity. Therefore dilute flows have to be investigated in small channels with a mean free path length of molecules higher than the channel width of the the micro channel itself (λ₀≥H₀). Modelling this process by a continuous flow the boundary conditions have to be modified (e.g. [6]). The present model uses the statistical approximation of the molecule velocity distribution to simulate the behaviour of this discrete flow with a weighted averaged molecule velocity ∼ξ_i, its standard deviation σ and the characterisic molecule collision rate z. The number density N per volume V near one position is used for the weighting factor averaging method describing the mean molecule velocity. The present model is validated computing Poiseuille and Couette flows with different Knudsen numbers. Showing the advantages of the present model the simulation results are compared with simulation results of the wall–distance depending diffusivity model of Lockerby and Reese [4] and BGK results of a Lattice–Boltzmann simulation. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)