Theoretical derivation of tangential velocity profiles in a flat vortex chamber-influence of turbulence and wall friction

Summary As part of a study on the hydrodynamics of a cyclone separator, a theoretical investigation of the flow pattern in a flat box cyclone (vortex chamber) has been carried out. Expressions have been derived for the tangential velocity profile as influenced by internal friction (eddy viscosity) and wall friction. The most important parameter controlling the tangential velocity profile is = –u ₀ R/(v+ ), where u ₀ is the radial velocity at the outer radius R of the cyclone, the kinematic liquid viscosity and is the kinematic eddy viscosity. For values of greater than about 10 the tangential velocity profile is nearly hyperbolic, for smaller than 1 the tangential velocity even decreases towards the centre. It is shown how and also the wall friction coefficient may be obtained from experimental velocity profiles with the aid of suitable graphs. Because of the close relation between eddy viscosity and eddy diffusion, measurements of velocity profiles in flat box cyclones will also provide information on the eddy motion of particles in a cyclone, a motion reducing its separation efficiency.List of symbols A cross-sectional area of cyclone inlet - h height of cyclone - p static pressure in cyclone - p static pressure difference in cyclone between two points on different radius - r radius in cyclone - r ₁ radius of cyclone outlet - R radius of cyclone circumference - u radial velocity in cyclone - u ₀ radial velocity at circumference of flat box cyclone - v tangential velocity - v ₀ tangential velocity at circumference of flat box cyclone - w axial velocity - z axial co-ordinate in cyclone - friction coefficient in flat box cyclone (for definition see § 5) - ₁ value of friction coefficient for ₁<< ₂ - ₂ value of friction coefficient for ₂<<1 - = - ₁ value of for ₁<< ₂ - ₂ value of for ₂<<1 - thickness of laminar boundary layer - =/h - turbulent kinematic viscosity - ratio of z to h - k ratio of height of cyclone to radius R of cyclone - parameter describing velocity profile in cyclone =–u ₀ R/(+) - kinematic viscosity of fluid - density of fluid - ratio of r to R - ₁ value of at outlet of cyclone - ₂ value of at inner radius of cyclone inlet - _w shear stress at cyclone wall - angular momentum in cyclone/angular momentum in cyclone inlet - ₁ value of at = ₁ - ₂ value of at = ₂     

Critical heat fluxes in boiling of liquid nitrogen underheated to saturation point in forced-flow conditions

This paper gives the results of experimental determinations of the critical heat fluxes in the boiling of Liquid nitrogen in forced-flow conditions in the mass velocity range 2 · 10³-40 · 10³ kg/m² · sec, pressure range 29 · 10⁴–245 · 10⁴ N/m², and at underheatings corresponding to the onset of normal boiling crises.Notation q₀ critical heat flux - r heat of vaporization - i enthalpy of flow corresponding to saturation point - i enthalpy of flow corresponding to liquid temperature - surface tension - density of liquid - density of saturated vapor - C _f friction factor - W_g mass velocity - Fr_* Froude number - g acceleration due to gravity     

Quantitative liquid jet instability measurement system using asymmetric magnification and digital image processing

A new technique for measuring the growth of instabilities on the surface of liquid jets flowing into gas is demonstrated. A collimated beam of white light illuminates the jet from behind, forming a shadow image. A pair of cylindrical lenses are arranged to provide different magnifications in the streamwise and cross-stream directions. A number of streamwise diameters and one cross-stream diameter are thus captured with maximum resolution in a single image on a charge-coupled device (CCD) electronic camera. A short-duration spark is used to freeze the jet motion. A mask representing the theoretical edge-response of the imaging system is digitally convolved with the cross-stream gray scale data to obtain sub-pixel resolution of the jet edge profile. The method is demonstrated using the well-known capillary jet instability and a ratio of streamwise to cross-stream magnifications of 40. Well-resolved single images show the development of the instability from small perturbations through the formation of the first drop. The system forms an accurate automated method of measuring the development of liquid jet instabilities. It can readily be applied to practical problems including liquid jet atomization.List of symbols a undisturbed jet radius - k nondimensional wavenumber (= 2a/) - Q gas-to-liquid density ratio - r ₀ mean jet radius, from initial region of image - R Reynolds number (= 2Ua/) - U mean jet velocity - We Weber number - z streamwise coordinate, origin at jet orifice - temporal growth rate - _s measured spatial growth rate - nondimensional temporal growth rate - r absolute value of height of peaks or troughs relative to r ₀ - r ₁ height of first extremum in a particular record - instability wavelength - liquid viscosity - liquid density - surface tension of liquid-gas interface     

The rheology of polymeric liquid crystals

The molecular theory of Doi has been used as a framework to characterize the rheological behavior of polymeric liquid crystals at the low deformation rates for which it was derived, and an appropriate extension for high deformation rates is presented. The essential physics behind the Doi formulation has, however, been retained in its entirety. The resulting four-parameter equation enables prediction of the shearing behavior at low and high deformation rates, of the stress in extensional flows, of the isotropic-anisotropic phase transition and of the molecular orientation. Extensional data over nearly three decades of elongation rate (10^–2–10¹) and shearing data over six decades of shear rate (10^–2–10⁴) have been correlated using this analysis. Experimental data are presented for both homogeneous and inhomogeneous shearing stress fields. For the latter, a 20-fold range of capillary tube diameters has been employed and no effects of system geometry or the inhomogeneity of the flow-field are observed. Such an independence of the rheological properties from these effects does not occur for low molecular weight liquid crystals and this is, perhaps, the first time this has been reported for polymeric lyotropic liquid crystals; the physical basis for this major difference is discussed briefly. A Semi-empirical constant in eq. (18), N/m² - c rod concentration, rods/m³ - c ^* critical rod concentration at which the isotropic phase becomes unstable, rods/m³ - C interaction potential in the Doi theory defined in eq. (3) - d rod diameter, m - D semi-empirical constant in eq. (19), s^–1 - D _r lumped rotational diffusivity defined in eq. (4), s^–1 - rotational diffusivity of rods in a concentrated (liquid crystalline) system, s^–1 - D _ro rotational diffusivity of a dilute solution of rods, s^–1 - f distribution function defining rod orientation - F tensorial term in the Doi theory defined in eq. (7) (or eq. (19)), s^–1 - G tensorial term in the Doi theory defined in eq. (8) - K _B Boltzmann constant, 1.38 × 10^–23 J/K-molecule - L rod length, m - S scalar order parameter - S tensor order parameter defined in eq. (5) - t time, s - T absolute temperature, K - u unit vector describing the orientation of an individual rod - rate of change ofu due to macroscopic flow, s^–1 - v fluid velocity vector, m/s - v velocity gradient tensor defined in eq. (9), s^–1 - V mean field (aligning) potential defined in eq. (2) - x coordinate direction, m - Kronecker delta (= 0 if = 1 if = ) - _r ratio of viscosity of suspension to that of the solvent at the same shear stress - _s solvent viscosity, Pa · s - ^* viscosity at the critical concentrationc ^*, Pa · s - v ₁, v₂ numerical factors in eqs. (3) and (4), respectively - deviatoric stress tensor, N/m² - volume fraction of rods - ⁰ constant in eq. (16) - ^* volume fraction of rods at the critical concentrationc ^* - average over the distribution functionf(u, t) (= d ²u f(u, t)) - gradient operator - d ²u integral over the surface of the sphere (|u| = 1)     

Thermal free convection from a solid sphere

Summary Thermal free convection from a sphere has been studied by melting solid benzene spheres in excess liquid benzene (Pr=8,3; 10⁸<Gr<10⁹). Overall heat transfer as well as local heat transfer were investigated. For the effect of cold liquid produced by the melting a correction has been applied. Results are compared with those obtained by other workers who used alternative experimental methods.Nomenclature coefficient of heat transfer - d characteristic length, here diameter of sphere - thermal conductivity - g acceleration of free fall - cubic expansion coefficient - T temperature difference between wall and fluid at infinity - kinematic viscosity - density - c specific heat capacity - a thermal diffusivity (=/c) - D diffusion coefficient - Nu dimensionless Nusselt number (=d/) - Nu* the analogous number for mass transfer (=kd/D) - mean value of Nusselt number - Gr dimensionless Grashof number (=gd ³T/ ²) - Gr* the analogous number for mass transfer (=gd ³x/ ²) - Pr dimensionless Prandtl number (=/a) - Sc dimensionless Schmidt number (=/D)     

A visualization study of the interaction of a free vortex with the wake behind an airfoil

The two-dimensional interaction of a single vortex with a thin symmetrical airfoil and its vortex wake has been investigated in a low turbulence wind tunnel having velocity of about 2 m/s in the measuring section. The flow Reynolds number based on the airfoil chord length was 4.5 × 10³. The investigation was carried out using a smoke-wire visualization technique with some support of standard hot-wire measurements. The experiment has proved that under certain conditions the vortex-airfoil-wake interaction leads to the formation of new vortices from the part of the wake positioned closely to the vortex. After the formation, the vortices rotate in the direction opposite to that of the incident vortex.List of symbols c test airfoil chord - C vortex generator airfoil chord - TA test airfoil - TE test airfoil trailing edge - TE _G vortex generator airfoil trailing edge - t dimensionless time-interval measured from the vortex passage by the test airfoil trailing edge: gDt=(T-T- _TE)·U/c - T time-interval measured from the start of VGA rotation - U free stream velocity - U vortex induced velocity fluctuation - VGA vortex generator airfoil - y distance in which the vortex passes the test airfoil - Z vortex circulation coefficient: Z=/(U · c/2) - vortex generator airfoil inclination angle - vortex circulation - vortex strength: =/2     

Flow of second-order fluids over an enclosed rotating disc

Summary This paper is devoted to a study of the flow of a second-order fluid (flowing with a small mass rate of symmetrical radial outflow m, taken negative for a net radial inflow) over a finite rotating disc enclosed within a coaxial cylinderical casing. The effects of the second-order terms are observed to depend upon two dimensionless parameters ₁ and ₂. Maximum values ₁ and ₂ of the dimensionless radial distances at which there is no recirculation, for the cases of net radial outflow (m>0) and net radial inflow (m<0) respectively, decrease with an increase in the second-order effects [represented by T(=₁+₂)]. The velocities at ₁ and ₂ as well as at some other fixed radii have been calculated for different T and the associated phenomena of no-recirculation/recirculation discussed. The change in flow phenomena due to a reversal of the direction of net radial flow has also been studied. The moment on the rotating disc increases with T.Nomenclature , , z coordinates in a cylindrical polar system - z ₀ distance between rotor and stator (gap length) - =/z ₀, dimensionless radial distance - =z/z ₀, dimensionless axial distance - _s = _s/z₀, dimensionless disc radius - V =(u, v, w), velocity vector - dimensionless velocity components - uniform angular velocity of the rotor - , p fluid density and pressure - P =p/(₂ z ₀₂ ² , dimensionless pressure - ₁, ₂, ₃ kinematic coefficients of Newtonian viscosity, elastico-viscosity and cross-viscosity respectively - ₁, _{2 2}/z ₀ ² , resp. ₃/z ₀ ² , dimensionless parameters representing the ratio of second-order and inertial effects - m = , mass rate of symmetrical radial outflow - l a number associated with induced circulatory flow - Rm =m/(z ₀₁), Reynolds number of radial outflow - R _l =l/(z ₀₁), Reynolds number of induced circulatory flow - Rz =z ₀ ² /₁, Reynolds number based on the gap - ₁, ₂ maximum radii at which there is no recirculation for the cases Rm>0 and Rm<0 respectively - ₁(T), ₂(T) ₁ and ₂ for different T - U _1(T) ⁽⁺⁾ = dimensionless radial velocity, Rm>0 - V _1(T) ⁽⁺⁾ = , dimensionless transverse velocity, Rm>0 - U _2(T) ^(–) = , dimensionless radial velocity, Rm=–Rn<0, m=–n - V _2(T) ^(–) = , dimensionless transverse velocity, Rm<0 - C _m moment coefficient     

A nonequilibrium theory of a slurry

A nonequilibrium theory of a slurry is developed and its practical use is illustrated by a simple stability analysis. Here a slurry is defined as a deformable continuum consisting of a liquid phase containing in suspension a large number of small solid particles which have formed by solidification from the liquid. The liquid is assumed to consist of two components and the solid to contain only one of the two. Consequently, the process of change of phase requires redistribution of material on the scale of the solid particles. This process is assumed to take a finite amount of time, requiring a nonequilibrium macroscopic theory. This theory contains four thermodynamic variables, three to represent the equilibrium state of the binary system and a fourth measuring the departure from thermodynamic equilibrium. The process of microscale diffusion of material is parameterized in the macroscale theory, leading to a Landau-type relaxation term in the equation of evolution of the fourth variable. The theory is simplified to yield a Boussinesq-like set of governing equations. Their practical use is illustrated by analyzing the stability of a simple steady solution of the equations and the effects of a non-zero relaxation time are discussed. A novel instability mechanism involving sedimentation of particles, previously found to occur in the equilibrium case, is found to persist in nonequilibrium, but disappears in the limit of no change of phase.Key to symbols a, b, c thermodynamic coefficients; see (3.36)–(3.38) - sedimentation coefficient; see (5.18) - C _p specific heat; see (3.24) - C _p ^de specific heat of the slurry; see (3.28) and (3.30) - c radius of solid particle (in §4) - D, D diffusive coefficients; see (3.40) and (3.41) - material diffusivity in liquid phase - D ^* modified diffusion coefficient; see (5.15) - d thermodynamic coefficient; see (3.39) - E specific internal energy - f, g, h thermodynamic coefficients; see (3.36)–(3.38) - g acceleration of gravity - reduced gravity; see (5.10) - i total diffusive flux vector of constituent 1 - i diffusive flux vector of constituent 1 in the liquid phase - j diffusive flux vector of solid phase - k thermal conductivity - k entropy flux vector - k _T, k_T thermodiffusion coefficients; see (3.40) and (3.41) - L latent heat of solidification per unit mass; see (3.7) and (3.24) - m wave number - m ^s rate of creation of mass of solid per unit volume through solidification - m ₁ ^s rate of creation of mass of solid constituent 1 per unit volume through solidification - mass rate of freezing per unit area per unit time - N number of solid particles per unit volume - p pressure - p _H hydrostatic component of pressure - p _m mechanical pressure - p ₁ dynamic component of pressure - q heat flux vector - Q _D rate of regeneration of heat through diffusive fluxes - Q _M rate of regeneration of heat through phase-change processes - Q _v rate of regeneration of heat through viscosity - Q vector defined by (3.16) - r heat externally supplied per unit mass (in §3); spherical radial coordinate (in §4) - S specific entropy of slurry - change of specific entropy with mass fraction of constituent 1; also change of chemical potential of liquid phase with temperature barring change of phase - change of chemical potential of liquid phase with temperature in phase equilibrium; see (3.28) and (3.30) - T temperature - t time - t ₀ relaxation time; see (5.30) - u barycentric velocity - u _H horizontal perturbation velocity - V sedimentation speed - w _a upward speed of simple state; see (6.5) and (6.12) - z upward vertical coordinate - upward unit vector - thermal expansion coefficient barring change of phase; see (3.23) - > ^* thermal expansion coefficient in phase equilibrium; see (3.27) and (3.30) - modified thermal expansion coefficient; see (5.1) and (5.4) - isothermal compressibility of slurry barring change of phase; see (3.23) - ^* isothermal compressibility of slurry in phase equilibrium; see (3.27) and (3.30) - dimensionless measure of departure from liquidus equilibrium; see (5.2) - _a deviation from phase equilibrium in simple state; see (6.6) and (6.13) - vertical wave number - volume expansion per unit mass upon melting; see (3.6) - change of chemical potential of liquid phase with pressure; see (3.25) - change of chemical potential of liquid phase with pressure for slurry; see (3.29) and (3.30) - compositional gradient in the static state; see (6.15) - vector defined by (3.35) - constant of integration; see (6.7) and (6.8) - coefficient defined by (6.23) - nonequilibrium expansion coefficient; see (5.1) and (5.4) - thermal diffusivity; =k/C _p - modified thermal diffusivity; see (5.33) - relaxation rate to phase equilibrium; see (2.2) - ₁ relaxation rate to solid-composition equilibrium; see (2.3) - sedimentation coefficient; see (4.29) - horizontal wave number vector - sedimentation coefficient; see (4.30) - ^L , ^s chemical potential of constituent 1 relative to constituent 2 in liquid and solid phase per unit mass; see (2.6) - change of chemical potential of liquid with liquid composition; see (3.8) - coefficient defined by (3.10) - kinematic shear viscosity - total mass fraction of constituent 1 (i.e., solute) - ^L, ^s mass fraction of constituent 1 in liquid and solid phases - density of slurry - _s density of solid phase - - - , growth rate of disturbance - stress tensor - deviatoric stress tensor - dimensionless temperature; see (5,3) - _a constant of integration; see (6.7) - mass fraction of solid phase in slurry - _b vertical gradient of mass fraction of solid; see (6.1) - dimensionless measure of _b; see (6.22) - _c temporal gradient of mass fraction of solid; see (6.1) - specific Gibbs free energy; see (3.13) - _L,_s specific Gibbs free energy of liquid and solid phases; see (2.12) - measure of departure from liquidus equilibrium; see (2.14) - measure of departure from solidus equilibrium; see (2.5) - spherical polar coordinate (in §4); see (4.20); wave angle (in §6); see (6.38)     

The swirling radial jet

The similarity solution of the radial turbulent jet with weak swirl is discussed and a new solution of the radial turbulent jet with swirl is proposed without restrictions assumed in the weak swirl solution.Nomenclature e swirl parameter - k experimental constant - l non-negative constant - M, M , N, P integral invariants - q velocity component in -direction - q _max maximum velocity component in -direction - u radial velocity component - u _max maximum radial velocity component - v axial velocity component - w peripheral velocity component - w _max maximum peripheral velocity component - x radial coordinate - y transverse coordinate - angle introduced in (28) - characteristic width of a jet - (x, y) similarity variable (scaled x and y coordinate) - molecular kinematic viscosity - _T eddy kinematic viscosity - tangential coordinate - fluid density - turbulent shear stress in -direction - _xy , _y components of turbulent shear stress tensor - (x, y) stream function     

Parameter determination for thermodynamical models of the pseudoelastic behaviour of shape memory alloy by resistivity measurements and infrared thermography

Two thermodynamical models of pseudoelastic behaviour of shape memory alloys have been formulated. The first corresponds to the ideal reversible case. The second takes into account the hysteresis loop characteristic of this shape memory alloys.Two totally independent techniques are used during a loading-unloading tensile test to determine the whole set of model parameters, namely resistivity and infrared thermography measurements. In the ideal case, there is no difficulty in identifying parameters.Infrared thermography measurements are well adapted for observing the phase transformation thermal effects.Notations 1 austenite 2 martensite - ₍₎ Macroscopic infinitesimal strain tensor of phase - ₍₂₎ ^f Traceless strain tensor associated with the formation of martensite phase - Macroscopic infiniesimal strain tensor - Macroscopic infinitesimal strain tensor deviator - _f Trace - Equivalent strain - ^pe Macroscopic pseudoelastic strain tensor - x Distortion due to parent (austenite =1)product (martensite =2) phase transformation (traceless symmetric second order tensor) - M Total mass of a system - M() Total mass of phase - V Total volume of a system - V() Total volume of phase - z=M(2)/M Weight fraction of martensite - 1-z=M(1)/M Weight fraction of austenite - u ₀ ^* () Specific internal energy of phase (=1,2) - s ₀ ^* () Specific internal entropy of phase - Specific configurational energy - Specific configurational entropy - ₀ ^f (T) Driving force for temperature-induced martensitic transformation at stress free state ( ₀ ^f T) = T ^*–Ts ^*) - Kirchhoff stress tensor - Kirchhoff stress tensor deviator - Equivalent stress - Cauchy stress tensor - Mass density - K Bulk moduli (K ₀=K) - L Elastic moduli tensor (order 4) - E Young modulus - Energetic shear (₀ = ) - Poisson coefficient - M _s ^o (M _F ^o ) Martensite start (finish) temperature at stress free state - A _s ^o (A _F ^o ) Austenite start (finish) temperature at stress free state - C _v Specific heat at constant volume - k Conductivity - Pseudoelastic strain obtained in tensile test after complete phase transformation (AM) (unidimensional test) - ₀ Thermal expansion tensor - r Resistivity - 1MPa 10⁶ N/m ² - ⁽⁾ Specific free energy of phase - _n Specific free energy at non equilibrium (R model) - _n ^eq Specific free energy at equilibrium (R model) - _n ^v Volumic part of ^eq - Specific free energy at non equilibrium (R _L model) - _conf Specific coherency energy (R _L model) - _c Specific free energy at constrained equilibria (R _L model) - _it (T) Coherency term (R _L model)     

The contribution of internal degrees of freedom to the non-Newtonian rheology of model polymer fluids

We report non-equilibrium molecular dynamics simulations of rigid and non-rigid dumbbell fluids to determine the contribution of internal degrees of freedom to strain-rate-dependent shear viscosity. The model adopted for non-rigid molecules is a modification of the finitely extensible nonlinear elastic (FENE) dumbbell commonly used in kinetic theories of polymer solutions. We consider model polymer melts — that is, fluids composed of rigid dumbbells and of FENE dumbbells. We report the steady-state stress tensor and the transient stress response to an applied Couerte strain field for several strain rates. We find that the rheological properties of the rigid and FENE dumbbells are qualitatively and quantitatively similar. (The only exception to this is the zero strain rate shear viscosity.) Except at high strain rates, the average conformation of the FENE dumbbells in a Couette strain field is found to be very similar to that of FENE dumbbells in the absence of strain. The theological properties of the two dumbbell fluids are compared to those of a corresponding fluid of spheres which is shown to be the most non-Newtonian of the three fluids considered.Symbol Definition b dimensionless time constant relating vibration to other forms of motion - F force on center of mass of dumbbell - F _i force on bead i of dumbbell - F force between center of masses of dumbbells and - F _ij force between beads i and j - h vector connecting bead to center of mass of dumbbell - H dimensionless spring constant for dumbbells, in units of / ² - I moment of inertia of dumbbell - J general current induced by applied field - k _B Boltzmann's constant - L angular momentum - m mass of bead, (= m/2) - M mass of dumbbell, g - N number of dumbbells in simulation cell - P translational momentum of center of mass of dumbbell - P pressure tensor - P _xy xy component of pressure tensor - Q separation of beads in dumbbell - Q _eq equilibrium extension of FENE dumbbell and fixed extension of rigid dumbbell - Q ₀ maximum extension of dumbbell - r _ij vector connecting beads i and j - r position vector of center of mass dumbbell - R vector connecting centers of mass of two dumbbells - t time - t ^* dimensionless time, in units of m/ - T ^* dimensionless temperature, in units of /k - u potential energy - u velocity vector of flow field - u _x x component of velocity vector - V volume of simulation cell - X general applied field - strain rate, s^–1 - ^* dimensionless shear rate, in units of /m ² - general transport property - Lennard-Jones potential well depth - friction factor for Gaussian thermostat - shear viscosity, g/cms - ^* dimensionless shear viscosity, in units of m/ ² - ^* dimensionless number density, in units of ^–3 - Lennard-Jones separation of minimum energy - relaxation time of a fluid - angular velocity of dumbbell - orientation angle of dumbbell      

Control of low-speed turbulent separated flow using jet vortex generators

A parametric study has been performed with jet vortex generators to determine their effectiveness in controlling flow separation associated with low-speed turbulent flow over a two-dimensional rearward-facing ramp. Results indicate that flow-separation control can be accomplished, with the level of control achieved being a function of jet speed, jet orientation (with respect to the free-stream direction), and jet location (distance from the separation region in the free-stream direction). Compared to slot blowing, jet vortex generators can provide an equivalent level of flow control over a larger spanwise region (for constant jet flow area and speed).Nomenclature C _p pressure coefficient, 2(P-P_)/V ² - C _Q total flow coefficient, Q/ v - D ₀ jet orifice diameter - Q total volumetric flow rate - R Reynolds number based on momentum thickness - u fluctuating velocity component in the free-stream (x) direction - V free-stream flow speed - VR ratio of jet speed to free-stream flow speed - x coordinate along the wall in the free-stream direction - jet inclination angle (angle between the jet axis and the wall) - jet azimuthal angle (angle between the jet axis and the free-stream direction in a horizontal plane) - boundary-layer thickness - momentum thickness - lateral distance between jet orifices A version of this paper was presented at the 12th Symposium on Turbulence, University of Missouri-Rolla, 24–26 Sept. 1990     

Supersonic air flow past a sphere with account for vibrational relaxation

A numerical solution is obtained for the problem of air flow past a sphere under conditions when nonequilibrium excitation of the vibrational degrees of freedom of the molecular components takes place in the shock layer. The problem is solved using the method of [1]. In calculating the relaxation rates account was taken of two processes: 1) transition of the molecular translational energy into vibrational energy during collision; 2) exchange of vibrational energy between the air components. Expressions for the relaxation rates were computed in [2]. The solution indicates that in the state far from equilibrium a relaxation layer is formed near the sphere surface. A comparison is made of the calculated values of the shock standoff with the experimental data of [3].Notation uV_max, vV_max velocity components normal and tangential to the sphere surface - V_max maximal velocity - P V _max ² pressure - density - TT temperature - e_viRT vibrational energy of the i-th component per mole (i=–O₂, N₂) - =rb–1 shock wave shape - a _f the frozen speed of sound - HRT/m gas total enthalpy     

Investigations of mean and turbulent flow characteristics of a two dimensional plane diffuser

This paper reports the investigation of mean and turbulent flow characteristics of a two-dimensional plane diffuser. Both experimental and theoretical details are considered. The experimental investigation consists of the measurement of mean velocity profiles, wall static pressure and turbulence stresses. Theoretical study involves the prediction of downstream velocity profiles and the distribution of turbulence kinetic energy using a well tested finite difference procedure. Two models, viz., Prandtl's mixing length hypothesis and k- model of turbulence, have been used and compared. The nondimensional static pressure distribution, the longitudinal pressure gradient, the pressure recovery coefficient, percentage recovery of static pressure, the variation of U _max/U _bar along the length of the diffuser and the blockage factor have been valuated from the predicted results and compared with the experimental data. Further, the predicted and the measured value of kinetic energy of turbulence have also been compared. It is seen that for the prediction of mean flow characteristics and to evaluate the performance of the diffuser, a simple turbulence model like Prandtl's mixing length hypothesis is quite adequate.List of symbols C ₁ , C ₂ ,C turbulence model constants - F _x body force - k kinetic energy of turbulence - l _m mixing length - L length of the diffuser - u, v, w rms value of the fluctuating velocity - u, v, w turbulent component of the velocity - mean velocity in the x direction - _A average velocity at inlet - U _bar average velocity in any cross section - U _max maximum velocity in any cross section - V mean velocity in the y direction - W local width of the diffuser at any cross section - x, y coordinates - dissipation rate of turbulence - _m eddy diffusivity - Von Karman constant - mixing length constant - _l laminar viscosity - _eff effective viscosity - v kinematic viscosity - density - _k effective Schmidt number for k - effective Schmidt number for - stream function - non dimensional stream function     

An experimental study of gas-liquid slug flow

Experimental measurements were carried out for upward gas-liquid slug flow in a 50.8 mm diameter pipe. Parallel conductance wires were used to distinguish the Taylor bubbles and liquid slugs and to determine translation velocities and lengths, an electrochemical probe provided the magnitude and direction of the wall shear stress and a radio-frequency local probe was used for the axial and radial distribution of voidage in the liquid slugs. Data are reported over wide range of flow conditions covering slug flow and into the churn flow pattern. Comparison with the Fernandes model predictions are presented. Numerical simulation of slug flow provided information on the structure of flow in a liquid slug and, in particular, on the process of mixing behind a Taylor bubble.List of symbols D pipe diameter - f Taylor bubble frequency - F _Gi (x) gas existence function for i-th liquid slug - g gravitational acceleration - l _A distance for the wall shear stress reversal in a liquid slug - l _B distance for the wall shear stress reversal in a Taylor bubble region - l _LS length of a liquid slug - l _TB length of a Taylor bubble - n number of samples in an ensemble - u axial velocity - U _M superficial mixture velocity (U _SG + U_SL) - U _N translation velocity of the leading Taylor bubble - U _NLS average translation velocity of liquid slugs - U _NTB average translation velocity of Taylor bubbles - U _OT overtaking velocity of the trailing Taylor bubble - U _SG superficial gas velocity - U _SL superficial liquid velocity - v radial velocity - w (y) velocity profile at the inlet to a liquid slug - x axial coordinate - y radial coordinate - void fraction - _LS void fraction in a liquid slug - =l _TB /(l_TB + l_LS) - density - surface tension - shear stress - saturation ratio, = _w / g h - ensemble average     

Transient flow towards a well in an aquifer including the effect of fluid inertia

Transient propagation of weak pressure perturbations in a homogeneous, isotropic, fluid saturated aquifer has been studied. A damped wave equation for the pressure in the aquifer is derived using the macroscopic, volume averaged, mass conservation and momentum equations. The equation is applied to the case of a well in a closed aquifer and analytical solutions are obtained to two different flow cases. It is shown that the radius of influence propagates with a finite velocity. The results show that the effect of fluid inertia could be of importance where transient flow in porous media is studied.List of symbols b Thickness of the aquifer, m - c ₀ Wave velocity, m/s - k Permeability of the porous medium, m² - n Porosity of the porous medium - p( ,t) Pressure, N/m² - Q Volume flux, m³/s - r Radial coordinate, m - r _w Radius of the well, m - s Transform variable - S Storativity of the aquifer - S _d(r, t) Drawdown, m - t Time, s - T Transmissivity of the aquifer, m²/s - ( ,t) Velocity of the fluid, m/s - Coordinate vector, m - z Vertical coordinate, m - Coefficient of compressibility, m²/N - Coefficient of fluid compressibility, m²/N - Relaxation time, s - (r, t) Hydraulic potential, m - Dynamic viscosity of the fluid, Ns/m² - Dimensionless radius - Density of the fluid, Ns²/m⁴ - (, ) Dimensionless drawdown - Dimensionless time - , x Dummy variables - ₀, ₁ Auxilary functions     

Streamwise pseudo-vortical structures and associated vorticity in the near-wall region of a wall-bounded turbulent shear flow

Streamwise pseudo-vortical motions near the wall in a fully-developed two-dimensional turbulent channel flow are clearly visualized in the plane perpendicular to the flow direction by a sophisticated hydrogen-bubble technique. This technique utilizes partially insulated fine wires, which generate hydrogen-bubble clusters at several distances from the wall. These flow visualizations also supply quantitative data on two instantaneous velocity components, and w, as well as the streamwise vorticity, _x . The vorticity field thus obtained shows quasi-periodicity in the spanwise direction and also a double-layer structure near the wall, both of which are qualitatively in good agreement with a pseudo-vortical motion model of the viscous wall-region.List of symbols C _i ,c _i ,d _i constants in Eqs. (2), (3) and (4) - H channel width (m) - Re _H Reynolds number (= U _c H/) - Re Reynolds number (= U _c /) - T period (s) - t time (s) - U mean streamwise velocity (m/s) - U _c center-line velocity (m/s) - u friction velocity (m/s) - u, , w velocity fluctuations (m/s) - x, y, z coordinates (m) - ^* displacement thickness (m) - momentum thickness (m) - mean low-speed streak spacing (m) - kinematic viscosity (m²/s) - phase difference - _x streamwise vorticity fluctuation (1/s) - ( )⁺ normalized by u and - (^–) root mean square value - (^–) statistical average This paper was presented at the Ninth Symposium on Turbulence, University of Missouri-Rolla, October 1–3, 1984