Linear stability of co-flowing liquid–gas jets

2001 ◽  
Vol 448 ◽  
pp. 23-51 ◽  
Author(s):  
J. M. GORDILLO ◽  
M. PÉREZ-SABORID ◽  
A. M. GAÑÁN-CALVO

A temporal, inviscid, linear stability analysis of a liquid jet and the co-flowing gas stream surrounding the jet has been performed. The basic liquid and gas velocity profiles have been computed self-consistently by solving numerically the appropriate set of coupled Navier–Stokes equations reduced using the slenderness approximation. The analysis in the case of a uniform liquid velocity profile recovers the classical Rayleigh and Weber non-viscous results as limiting cases for well-developed and very thin gas boundary layers respectively, but the consideration of realistic liquid velocity profiles brings to light new families of modes which are essential to explain atomization experiments at large enough Weber numbers, and which do not appear in the classical stability analyses of non-viscous parallel streams. In fact, in atomization experiments with Weber numbers around 20, we observe a change in the breakup pattern from axisymmetric to helicoidal modes which are predicted and explained by our theory as having an hydrodynamic origin related to the structure of the liquid-jet basic velocity profile. This work has been motivated by the recent discovery by Gañán-Calvo (1998) of a new atomization technique based on the acceleration to large velocities of coaxial liquid and gas jets by means of a favourable pressure gradient and which are of emerging interest in microfluidic applications (high-quality atomization, micro-fibre production, biomedical applications, etc.).

Author(s):  
James Q. Feng

The steady axisymmetric behaviour of a relatively small bubble moving with a flowing liquid in a straight round tube is studied by computationally solving the nonlinear Navier–Stokes equations, using a Galerkin finite-element method with boundary-fitted mesh, for wide ranges of capillary number C a and Reynolds number R e . Here a bubble is considered relatively small when its volume-equivalent radius is less than that of the tube. At small values of R e , the velocity of a bubble increases with bubble size for large values of C a but decreases with bubble size for small values of C a . At large values of R e , however, a bubble of large size appears to move at a slower velocity for any given value of C a . When R e is large (e.g. R e = 100) and C a > 0.1, a bubble of radius greater than half of the tube radius moves at a velocity that seems to be independent of bubble size. The strong inertial effect at large R e makes a small bubble of radius greater than a quarter of the tube radius to deform into a noticeable oblate shape as C a increases from very small value, and then to be elongated into a bullet shape with further increasing C a after C a reaches an intermediate value. Even very small bubbles (e.g. of radius equal to one-tenth of the tube radius) can still be significantly deformed provided that the value of C a is adequately large. Despite significant shape deformations that may still occur, bubbles of radius less than a quarter of that of the tube almost always move at the same velocity as that of the local liquid flow at the tube centreline (i.e. twice that of the average liquid velocity), regardless the values of R e and C a . This fact suggests that very small bubbles are basically carried by the local liquid flow.


2020 ◽  
Author(s):  
RAJDEEP TAH ◽  
SARBAJIT MAZUMDAR ◽  
Krishna Kant Parida

The shape of the liquid surface for a fluid present in a uniformly rotating cylinder is generally determined by making a Tangential velocity gradient along the radius of the rotating cylindrical container. A very similar principle can be applied if the direction of the produced velocity gradient is reversed, for which the source of rotation will be present at the central axis of the cylindrical vessel in which the liquid is present. Now if the described system is completely closed, the angular velocity will decrease as a function of time. But when the surface of the rotating fluid is kept free, then the Tangential velocity profile would be similar to that of the Taylor-Couette Flow, with a modification that; due to formation of a curvature at the surface, the Navier-Stokes law is to be modified. Now the final equation may not seem to have a proper general solution, but can be approximated to certain solvable expressions for specific cases of angular velocity.


Author(s):  
V.N. Popov ◽  
A.N. Cherepanov

The purpose of the research was to numerically simulate the processes when melting drops fall on a substrate. The paper deals with the solidification on the metal surface of a binary aluminum alloy modified by activated refractory nanosized particles, which are the centers of crystalline phase nucleation. We formulated a mathematical model which describes the thermo- and hydrodynamic phenomena in the drop upon interaction with a solid substrate, heterogeneous nucleation during melt cooling, and subsequent crystallization. The flow in a liquid is described by the Navier --- Stokes equations in the Boussinesq approximation. The position of the free boundary of the melt is fixed by marker particles moving with the local liquid velocity. On the melt --- substrate contact surface, thermal resistance is taken into account. The hydrodynamic problem is considered under conditions of crystallization of molten metal. The temperature conditions and the kinetics of the growth of the solid phase in the solidifying aluminum alloy are described for various sizes of formed splats. Satisfactory agreement was found between the shape of the splat obtained by the results of numerical simulation and the available experimental data. The adequacy of the crystallization model in the presence of ultradisperse refractory particles in a binary alloy is confirmed. It was determined that, regardless of the size of the drop, bulk crystallization of the metal takes place. It was found that at a high rate of collision of a drop with a substrate during the melt spreading, a small fraction of the solid phase can be formed.


1996 ◽  
Vol 06 (02) ◽  
pp. 409-413
Author(s):  
A. J. CONLEY

The flow of an incompressible viscous fluid between parallel plates becomes unstable when the plates are tumbled. As the tumbling rate increases, the flow restabilizes. This phenomenon is elucidated by path-following techniques. The solution of the Navier-Stokes equations is approximated by spectral techniques. The linear stability of these solutions is studied.


1996 ◽  
Vol 314 ◽  
pp. 227-246 ◽  
Author(s):  
Philip A. J. Mees ◽  
K. Nandakumar ◽  
J. H. Masliyah

Steady developing flow of an incompressible Newtonian fluid in a curved duct of square cross-section (the Dean problem) is investigated both experimentally and numerically. This study is a continuation of the work by Bara, Nandakumar & Masliyah (1992) and is focused on flow rates between Dn = 200 and Dn = 600 (Dn = Re/(R/a)1/2, where Re is the Reynolds number, R is the radius of curvature of the duct and a is the duct dimension; the curvature ratio, R/a, is 15.1).Numerical simulations based on the steady three-dimensional Navier – Stokes equations predict the development of a 6-cell secondary flow pattern above a Dean number of 350. The 6-cell state consists of two large Ekman vortices and two pairs of small Dean vortices near the outer wall that result from the primary instability that is of centrifugal nature. The 6-cell flow state develops near θ = 80° and breaks down symmetrically into a 2-cell flow pattern.The apparatus used to verify the simulations had a duct dimension of 1.27 cm and a streamwise length of 270°. At a Dean number of 453, different velocity profiles of the 6-cell flow state at θ = 90° and spanwise profiles of the streamwise velocity at every 20° were measured using a laser-Doppler anemometer. All measured velocity profiles, as well as flow visualization of secondary flow patterns, are in very good agreement with the simulations, indicating that the parabolized Navier – Stokes equations give an accurate description of the flow.Based on the similarity with boundary layer flow over a concave wall (the Görtler problem), it is suggested that the transition to the 6-cell flow state is the result of a decreasing spanwise wavelength of the Dean vortices with increasing flow rate. A numerical stability analysis shows that the 6-cell flow state is unconditionally unstable. This is the first time that detailed experiments and simulations of the development of a 6-cell flow state are reported.


Author(s):  
Ehsan Gholamalizadeh ◽  
Farzad Pahlevanzadeh ◽  
Kamal Ghani ◽  
Arash Karimipour ◽  
Truong Khang Nguyen ◽  
...  

Purpose This study aims to numerically study the forced convection effects on a two-dimensional microchannel filled with a porous material containing the water/FMWCNT nanofluid. The upper and lower microchannel walls were fully insulated thermally along 15 per cent of their lengths at each end of the microchannel, with the in-between length being exposed to a constant temperature. The slip velocity boundary condition was applied along the microchannel walls. Design/methodology/approach The Navier–Stokes equations were discretized before being solved numerically via a FORTRAN computer code. The following ranges were considered for the studied parameters: slip factor (B) equal to 0.001, 0.01 and 0.1; Reynolds number (Re) between 10 and 100; solid nanoparticle mass fraction (ϕ) between 0.0012 and 0.0025; Darcy number (Da) between 0.001 and 0.1; and porosity factor (ε) between 0.4 and 0.9. Findings Increasing the Da caused a greater increase in the velocity profile than increasing Re, whereas increasing porosity did not affect the velocity profile growth at all. Originality/value This paper is the continuation of the authors’ previous studies. Using the water/FMWCNT nanofluid as the working fluid in microchannels is among the achievements of this study.


1974 ◽  
Vol 96 (4) ◽  
pp. 333-340 ◽  
Author(s):  
J. E. R. Coney ◽  
M. A. I. El-Shaarawi

The boundary layer simplification of the Navier-Stokes equations for hydrodynamically developing laminar flow with constant physical properties in the entrance region of concentric annuli with rotating inner walls have been numerically solved using a simple linearized finite-difference scheme. Additional results to those existing in the literature by Martin and Payne [1–2] will be presented here. An advantage of the analysis used in this paper is that it does not solve for the stream function and vorticity, but predicts the development of tangential, axial and radial velocity profiles directly, thus avoiding numerical differentiation. Results for the development of these velocity profiles, pressure drop and friction factor are presented for five annuli radii ratios (0.3, 0.5, 0.674, 0.727 and 0.90) at various values of the parameter Re2/Ta. The paper may be considered as a direct comparison between the boundary layer solution and the complete solution of the Navier-Stokes equations [1–2] for that special case.


1992 ◽  
Vol 114 (4) ◽  
pp. 504-511 ◽  
Author(s):  
X. Y. Xu ◽  
M. W. Collins ◽  
C. J. H. Jones

Three-dimensional flows through canine femoral bifurcation models were predicted under physiological flow conditions by solving numerically the time-dependent threedimensional Navier-stokes equations. In the calculations, two models were assumed for the blood, those of (a) a Newtonian fluid, and (b) a non-Newtonian fluid obeying the power law. The blood vessel wall was assumed to be rigid this being the only approximation to the prediction model. The numerical procedure utilized a finite volume approach on a finite element mesh to discretize the equations, and the code used (ASTEC) incorporated the SIMPLE velocity-pressure algorithm in performing the calculations. The predicted velocity profiles were in good qualitative agreement with the in vivo measurements recently obtained by Jones et al. [1]. The non-Newtonian effects on the bifurcation flow field were also investigated, and no great differences in velocity profiles were observed. This indicated that the non-Newtonian characteristics of the blood might not be an important factor in determining the general flow patterns for these bifurcations, but could have local significance. Current work involves modeling wall distensibility in an empirically valid manner. Predictions accommodating these will permit a true quantitative comparison with experiment.


2019 ◽  
Vol 72 (4) ◽  
pp. 545-562 ◽  
Author(s):  
C V Tran ◽  
X Yu

Summary In this article, we present a new kind of regularity criteria for the global well-posedness problem of the three-dimensional Navier–Stokes equations in the whole space. The novelty of the new results is that they involve only the profiles of the magnitude of the velocity. One particular consequence of our theorem is as follows. If for every fixed $t\in (0,T)$, the ‘large velocity’ region $\Omega:=\{(x,t)\mid |u(x,t)|>C(q)\left|\mkern-2mu\left|{u}\right|\mkern-2mu\right|_{L^{3q-6}}\}$, for some $C(q)$ appropriately defined, shrinks fast enough as $q\nearrow \infty$, then the solution remains regular beyond $T$. We examine and discuss velocity profiles satisfying our criterion. It remains to be seen whether these profiles are typical of general Navier–Stokes flows.


A method is described of predicting the growth of a tangential velocity profile in fully developed laminar axial flow through a concentric annulus when the inner surface is rotated at speeds which are insufficient to generate Taylor vortices. The treatment, which is based on simplification and subsequent solution of the Navier-Stokes equations, as Fourier-Bessel series, appears preferable to momentum-integral techniques through greater simplicity of expression and in requiring fewer assumptions about the developing tangential profile. The validity of the predictions is best at high axial Reynolds number.


Sign in / Sign up

Export Citation Format

Share Document