Dynamics of Brownian particles in a turbulent channel flow

Turbulent channel flows with suspended particles are investigated by means of numerical simulations. The fluid velocity is computed by large eddy simulation. Motion of small graphite particles with diameter of 0.01–10 m, corresponding to the Schmidt number, Sc, of 2.87 × 10²–6.22 × 10⁶ and the particle relaxation time in wall unit, _p⁺, of 9.79 × 10^–5–4.51, is computed by Lagrangian particle tracking. Relation between the particle relaxation time and the computed deposition velocity is found to be in good agreement with an empirical relation. The statistics of the particle motion in the vicinity of the wall are studied. Clear differences are found in dynamical behavior of particles with different sizes. Medium size particles show a strong dependence on the structure of the fluid flow, while small and large particles are considerably less sensitive.     

Landslide generated impulse waves.

Landslide generated impulse waves were investigated in a two-dimensional physical laboratory model based on the generalized Froude similarity. Digital particle image velocimetry (PIV) was applied to the landslide impact and wave generation. Areas of interest up to 0.8 m by 0.8 m were investigated. The challenges posed to the measurement system in an extremely unsteady three-phase flow consisting of granular matter, air, and water were considered. The complex flow phenomena in the first stage of impulse wave initiation are: high-speed granular slide impact, slide deformation and penetration into the fluid, flow separation, hydrodynamic impact crater formation, and wave generation. During this first stage the three phases are separated along sharp interfaces changing significantly within time and space. Digital masking techniques are applied to distinguish between phases thereafter allowing phase separated image processing. PIV provided instantaneous velocity vector fields in a large area of interest and gave insight into the kinematics of the wave generation process. Differential estimates such as vorticity, divergence, elongational, and shear strain were extracted from the velocity vector fields. The fundamental assumption of irrotational flow in the Laplace equation was confirmed experimentally for these non-linear waves. Applicability of PIV at large scale as well as to flows with large velocity gradients is highlighted.List of symbols a wave amplitude (L) - c wave celerity (LT^–1) - d_diff diffraction limited minimum particle image diameter (L) - d_e diffracted particle image diameter (L) - d_g granulate grain diameter (L) - d_p seeding particle diameter (L) - d recorded particle image diameter (L) - f focal length (L) - f_# f number (-) - F slide Froude number (-) - g gravitational acceleration (LT^–2) - h still-water depth (L) - H wave height (L) - l_s slide length (L) - L wavelength (L) - M magnification (-) - m_s slide mass (M) - n refractive index (-) - n_por slide porosity (-) - N_iw number of seeding particles in an interrogation window (-) - N_pair number of detected particle image pairs in window (-) - p interrogation window size p×p pixels; 1 pixel=9 m (L) - P probability (-) - P_il probability of in-plane loss of particle (-) - P_ol probability of out-of-plane loss of particle (-) - s slide thickness (L) - S relative slide thickness (-) - t time after impact (T) - T wave period (T) - v velocity (LT^–1) - v_p particle velocity (LT^–1) - v_px streamwise horizontal component of particle velocity (LT^–1) - v_py crosswise horizontal component of particle velocity (LT^–1) - v_pz vertical component of particle velocity (LT^–1) - v_s slide centroid velocity at impact (LT^–1) - V dimensionless slide volume (-) - V_iw interrogation volume (L³) - V_s slide volume (L³) - x streamwise coordinate (L) - x_ip area of view x dimension in image plane (L) - z vertical coordinate (L) - slide impact angle (°) - bed friction angle (°) - y depth of field (L) - t laser pulse separation (T) - x mean particle image x displacement in interrogation window (L) - _x random displacement x error (L) - _v random velocity v error (LT^–1) - _tot total random velocity v error (LT^–1) - _bias velocity v error due to biased correlation analysis (LT^–1) - _optics velocity v error due to optical imaging errors (LT^–1) - _track velocity v error due to particle flow tracking error (LT^–1) - _xx streamwise horizontal elongational strain component (1/T) - _xz shear strain component (1/T) - _zx shear strain component (1/T) - _zz vertical elongational strain component (1/T) - water surface displacement (L) - wavelength (L) - dynamic viscosity (ML^–1T^–1) - density (ML^–3) - _g granulate density (ML^–3) - _p particle density (ML^–3) - _s mean slide density (ML^–3) - _w water density (ML^–3) - granulate internal friction angle (°) - _y vorticity vector component (out-of-plane) (1/T)     

Behavior of spherical solid particles released in a laminar boundary layer along a flat plate

The behavior of a very small solid spherical particle initially at rest on the surface of a flat plate in a laminar boundary layer along the plate is investigated. The Stokes drag is the only force considered to be acting on the particle. The fluid Reynolds number Re _f is assumed to be large, and the particle Reynolds number Re is assumed to be small. The equations describing the motion of the particle are two simultaneous, second order, nonlinear, ordinary differential equations with one parameter. A complete digital computer solution and analytic limiting solutions for large and small values of a dimensionless time have been obtained. The numerical and the analytic solutions are in close agreement. The results presented are the velocity, trajectory, and time history of the particle and the force acting on the particle. These results show that the particle comes into equilibrium with the fluid very quickly with respect to the spatial coordinates, rising only several radii from the surface in its entire flight.Nomenclature a diameter of the particle - F force acting on the particle - F ₁ lift force - F _M Magnus force - i, j unit vectors parallel and perpendicular to the plate, respectively - Re = u a ^–1 Reynolds number for the particle based on the reference fluid velocity - Re _f = u x ₀ ^–1, Reynolds number for the fluid - Re _p |w _f –w | a ^–1, Reynolds number for the particle - t time - u component of particle velocity parallel to the plate - u _f component of fluid velocity parallel to the plate - u free-stream fluid velocity - U , dimensionless component of particle velocity parallel to the plate - U _f Re _f Re ^–1 u ^–1 U _f, dimensionless component of fluid velocity parallel to the plate - v component of particle velocity perpendicular to the plate - v _f component of fluid velocity perpendicular to the plate - V , dimensionless component of particle velocity perpendicular to the plate - V _f Re _f ^3/2 Re ^–2 u v _f, dimensionless component of fluid velocity perpendicular to the plate - w u i+v j, velocity of the particle - w _f u _f i+v _f j, velocity of the fluid - w dw/dt - x horizontal distance of the center of the particle from the leading edge of the plate - x ₀ initial horizontal distance of the center of the particle from the leading edge of the plate - X xx ₀ ^–1 , dimensionless horizontal distance of the center of the particle from the leading edge of the plate - y vertical distance of the center of the particle from the plate - Y ya ^–1, dimensionless vertical distance of the center of the particle from the plate - ^–1 - 0.332 - Re Re _f ^–3/2 - - - viscosity of the fluid - _f ^–1 , kinematic viscosity of the fluid - _f density of the fluid - _p density of the particle - , dimensionless time - angular velocity of the particle     

Landslide generated impulse waves. 2. Hydrodynamic impact craters

Landslide generated impulse waves were investigated in a two-dimensional physical laboratory model based on the generalized Froude similarity. Digital particle image velocimetry (PIV) was applied to the landslide impact and wave generation. Areas of interest up to 0.8 m by 0.8 m were investigated. PIV provided instantaneous velocity vector fields in a large area of interest and gave insight into the kinematics of the wave generation process. Differential estimates such as vorticity, divergence, and elongational and shear strain were extracted from the velocity vector fields. At high impact velocities flow separation occurred on the slide shoulder resulting in a hydrodynamic impact crater, whereas at low impact velocities no flow detachment was observed. The hydrodynamic impact craters may be distinguished into outward and backward collapsing impact craters. The maximum crater volume, which corresponds to the water displacement volume, exceeded the landslide volume by up to an order of magnitude. The water displacement caused by the landslide generated the first wave crest and the collapse of the air cavity followed by a run-up along the slide ramp issued the second wave crest. The extracted water displacement curves may replace the complex wave generation process in numerical models. The water displacement and displacement rate were described by multiple regressions of the following three dimensionless quantities: the slide Froude number, the relative slide volume, and the relative slide thickness. The slide Froude number was identified as the dominant parameter.List of symbols a wave amplitude (L) - b slide width (L) - c wave celerity (LT^–1) - d _g granulate grain diameter (L) - d _p seeding particle diameter (L) - F slide Froude number - g gravitational acceleration (LT^–2) - h stillwater depth (L) - H wave height (L) - l _s slide length (L) - L wave length (L) - M magnification - m _s slide mass (M) - n _por slide porosity - Q _d water displacement rate (L³) - Q _D maximum water displacement rate (L³) - Q _s maximum slide displacement rate - s slide thickness (L) - S relative slide thickness - t time after impact (T) - t _D time of maximum water displacement volume (L³) - t _qD time of maximum water displacement rate (L³) - t _si slide impact duration (T) - t _sd duration of subaqueous slide motion (T) - T wave period (T) - v velocity (LT^–1) - v _p particle velocity (LT^–1) - v _px streamwise horizontal component of particle velocity (LT^–1) - v _pz vertical component of particle velocity (LT^–1) - v _s slide centroid velocity at impact (LT^–1) - V dimensionless slide volume - V _d water displacement volume (L³) - V _D maximum water displacement volume (L³) - V _s slide volume (L³) - x streamwise coordinate (L) - z vertical coordinate (L) - slide impact angle (°) - bed friction angle (°) - x mean particle image x-displacement in interrogation window (L) - _x random displacement x error (L) - _tot total random velocity v error (LT^–1) - _xx streamwise horizontal elongational strain component (1/T) - _xz shear strain component (1/T) - _zx shear strain component (1/T) - _zz vertical elongational strain component (1/T) - water surface displacement (L) - density (ML^–3) - _g granulate density (ML^–3) - _p particle density (ML^–3) - _s mean slide density (ML^–3) - _w water density (ML^–3) - granulate internal friction angle (°) - _y vorticity vector component (out-of-plane) (1/T)     

Diffusion to a chain of drops (bubbles) at large péclet numbers

The problem of diffusion of a substance, dissolved in a flow, to absorbing drops (bubbles) moving one after another in a viscous incompressible fluid is investigated. An approximate analytic expression is obtained for the differential and integral flows of the substance to the surface of each drop with consideration of the changes of the concentration and velocity fields due to the presence of other drops. A chain of spherical drops of equal radius arranged on the axis of a uniform forward flow is examined. It is shown that if the distance between drops, referred to the radius of the drops, satisfies the inequality 1lP^1/2 (P is the Péclet number), then the integral inflow of the substance to the surface of the second drop of the chain is 2.41 times less than the integral inflow to the first (the drops are enumerated along the flow); the total diffusion flow to the surface of an arbitrary drop with number k is determined by the expression I_k=I₁[k^1/2 – (k–1)^1/2], where I_k is the total flow to the first drop of the chain. The case of diffusion interaction of a solid particle and drop is examined. It is shown that for particles moving one after another with the same velocity in a fluid at rest the presence of a drop before the solid particle leads to a marked decrease of the total diffusion flow of the solid particle [by O(P^1/6) times], whereas the presence of a solid particle before a drop does not affect (in the main approximation with respect to the characteristic diffusion parameter) the total flow of the latter.I _k=I _i[k ^1/2–(k–1)^1/2]Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 59–69, January–February, 1978.     

Measurement of the Coefficient of Transverse Dispersion in Flow Through Packed Beds for a Wide Range of Values of the Schmidt Number

Experimental values of the coefficient of transverse dispersion (D _T) were measured with the system 2-naphthol/water, over a range of temperatures between 293K and 373K, which corresponds to a range of values of viscosity () between 2.83×10^–4 Ns/m² and 1.01×10^–3 Ns/m² and of molecular diffusion coefficient (D _m) between 1.03×10^–9 m²/s and 5.49×10^–9 m²/s. Since the density () of water is close to 10³ kg/m³, the corresponding variation of the Schmidt number (Sc=/D _m) was in the range 1000 – 50. More than 200 experimental values of the transverse dispersion coefficient were obtained using beds of silica sand with average particle sizes (d) of 0.297 and 0.496mm, operated over a range of interstitial liquid velocities (u) between 0.1mm/s and 14mm/s and this gave a variation of the Reynolds number (Re=du/) between 0.01 and 3.5.Plots of the dimensionless coefficient of transverse dispersion (D _T/D _m) vs. the Peclet number (Pe_m=ud/D _m) based on molecular diffusion bring into evidence the influence of Sc on transverse dispersion. As the temperature is increased, the value of Sc decreases and the values of D _T/D _m gradually approach the line corresponding to gas behaviour (i.e. Sc 1), which is known to be well approximated by the equation D _T/D _m=1/+ud/12D _m, where is the tortuosity with regard to diffusion.     

Effect of longitudinal riblets on the aerodynamic characteristics of a straight wing

The results of balance aerodynamic tests on model straight wings with smooth and ribbed surfaces at an angle of attack =–4°–12°, Mach number M=0.15–0.63, and Reynolds number Re=2.4·10⁶–3.5·10⁶ are discussed. The nondimensional riblet spacings ⁺, which determines the effect of the riblets on the turbulent friction drag, and the effect of riblets on the upper and/or lower surface of a straight wing on its drag, lift, and moment characteristics are estimated.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 33–38, March–April, 1995.     

Diffusion of moisture in deformable porous media. Anisotropic fibrous materials

This paper presents a study on the deformation of anisotropic fibrous porous media subjected to moistening by water in the liquid phase. The deformation of the medium is studied by applying the concept of effective stress. Given the structure of the medium, the displacement of the solid matrix is not taken into account with respect to the displacement of the liquid phase. The transport equations are derived from the model proposed by Narasimhan. The transport coefficients and the relation between the variation in apparent density and effective stress are obtained by test measurements. A numerical model has been established and applied for studying drip moistening of mineral wool samples capable or incapable of deformation.Nomenclature D mass diffusion coefficient [L²t^–1] - e void fraction - g gravity acceleration [Lt^–2] - J mass transfer density [ML^–2t^–1] - K hydraulic conductivity [Lt^–1] - K _s hydraulic conductivity of the solid phase [Lt^–1] - K ^* hydraulic conductivity of the deformable porous medium [Lt^–1] - P pressure of moistening liquid [ML^–1 t^–2] - S degree of saturation - t time [t] - V speed [Lt^–1] - X horizontal coordinate [L] - Z vertical coordinate measured from the bottom of porous medium [L] - z z-coordinate [L] Greek Letters porosity - ₁ total hydric potential [L] - _g gas density [ML^–3] - ₁ liquid density [ML^–3] - ₀ apparent density [ML^–3] - _s density of the solid phase [ML^–3] - density of the moist porous medium [ML^–3] - external load [ML^–1t^–2] - effective stress [ML^–1t^–2] - bishop's parameter - matrix potential or capillary suction [L] Indices g gas - 1 moistening liquid - p direction perpendicular to fiber planes - s solid matrix - t direction parallel to fiber planes - v pore Exponent * movement of solid particles taken into account     

Development of a full-field planar Mie scattering technique for evaluating swirling mixers

A full-field planar optical diagnostic technique for studying mixing in swirling flows is described. Results were obtained using this technique to provide planar mixing information by seeding a simulated fuel stream with aluminum oxide particles, then inferring concentration from Mie scattering intensity distributions. This facility and measurement technique are unique for several reasons. First, they allow spatial variations in laser sheet energy to be corrected for on a shot-to-shot basis. Second, they allow experiments to be performed for swirlers with practical fuel and oxidizer flow rates, i.e. on the order of 150 g/s (0.33 lbm/s). Finally, they allow full size swirler models to be evaluated, with the entire exit plane imaged simultaneously. Representative results are presented as false color images of the planar mixing fields. These images allow rapid assessment of the mixing process and its changes with variations in operating conditions or swirler geometry.List of Symbols C seed particle concentration, m^–3 - mean component of seed particle concentration, m^–3 - C fluctuating component of seed particle concentration, m^–3 - C ^* time averaged ratio of rms particle concentration fluctuations to average particle concentration, dimensionless - d _p particle diameter, m - I laser energy after passing through the flow, J/m² - mean laser energy, J/m² - I ₀ laser energy before passing through the flow, J/m² - L _v eddy length scale, m - l laser beam path length, m - U _v eddy velocity scale, m/s - V diode voltage reading after passing through the flow, V - mean diode voltage, V - V ₀ diode voltage reading before passing through the flow, V - absorptivity, m² - _rel relative equivalence ratio, dimensionless - fluid viscosity, Ns - _p particle density, kg/m³ - Stokes number= _p / _f , dimensionless - _f flow time scale, s - _p particle response time, s     

Rheological properties of skim milk gels at various temperatures; interrelation between the dynamic moduli and the relaxation modulus

The rheological properties of rennet-induced skim milk gels were determined by two methods, i.e., via stress relaxation and dynamic tests. The stress relaxation modulusG _c (t) was calculated from the dynamic moduliG andG by using a simple approximation formula and by means of a more complex procedure, via calculation of the relaxation spectrum. Either calculation method gave the same results forG _c (t). The magnitude of the relaxation modulus obtained from the stress relaxation experiments was 10% to 20% lower than that calculated from the dynamic tests.Rennet-induced skim milk gels did not show an equilibrium modulus. An increase in temperature in the range from 20° to 35 °C resulted in lower moduli at a given time scale and faster relaxation. Dynamic measurements were also performed on acid-induced skim milk gels at various temperatures andG _c (t) was calculated. The moduli of the acid-induced gels were higher than those of the rennet-induced gels and a kind of permanent network seemed to exist, also at higher temperatures. G storage shear modulus,N·m^–2; - G loss shear modulus,N·m^–2; - G _c calculated storage shear modulus,N·m^–2; - G _c calculated loss shear modulus,N·m^–2; - G _e equilibrium shear modulus,N·m^–2; - G _ec calculated equilibrium shear modulus,N·m^–2; - G(t) relaxation shear modulus,N·m^–2; - G _c (t) calculated relaxation shear modulus,N·m^–2; - G ^*(t) pseudo relaxation shear modulus,N·m^–2; - H relaxation spectrum,N·m^–2; - t time,s; - relaxation time,s; - angular frequency, rad·s^–1. Partly presented at the Conference on Rheology of Food, Pharmaceutical and Biological Materials, Warwick, UK, September 13–15, 1989 [33].     

Diffusion of14C in dense saturated bentonite under steady-state conditions

Diffusion coefficients are critical parameters for predicting migration rates and fluxes of contaminants through clay-based barrier materials used in many waste containment strategies. Cabon-14 is present in high-level nuclear fuel waste and also in many low-level wastes such as those generated from some medical research activities. Diffusion coefficients were measured for¹⁴C (in the form of carbonate) in bentonite compacted to a series of dry bulk densities, _b, ranging from about 0.9 to 1.6 Mg/m³. The clay was saturated with a Na-Ca-Cl-dominated groundwater solution typical of those found deep in plutonic rock on the Canadian Shield. Both effective,D _e, and apparent,D _a, diffusion coefficients were determined.D _e is defined asD _{0 a} n _e, where D₀ is the diffusion coefficient in pure bulk water, _a the apparent tortuosity factor, andn _e the effective porosity available for diffusion; andD _a is defined asD _{0 a} n _e/(n _e + _b K _d ), where K_d is the solid/liquid distribution coefficient. BothD _e andD _a decrease with increasing _b:D _e values range from about 10×10^–12 m²/s at _b0.9 Mg/m³ to 0.6×10^–12 m²/s at 1.6 Mg/m³, andD _a values vary from approximately 40×10^–12 to 4×10^–12 m²/s over the same density range. The decrease inD _e andD _a is attributed to a decrease in both _a andn _e as _b increases. The data indicate thatn _e is <10% of the total solution-filled porosity of the clay at all densities.K _d values for¹⁴C with the clay range from about 0.3 to <0.1 m³/Mg; this indicates there is a small amount of¹⁴C sorbed on the clay and/or some¹⁴C is isotopically exchanged with¹²C in carbonate phases present in the clay. Finally, theD _e values for¹⁴C are lower than those of other diffusants — I^–, Cl^–, TcO₄ ^–, and Cs⁺ — that have been measured in this clay and pore-water solution. This is attributed to lower values for bothn _e andD ₀ for¹⁴C species relative to those of the other diffusants.     

Shear stress and normal stress measurements of aqueous polymer solutions in a cone-and-plate rheogoniometer

Summary The rheological behaviour of aqueous solutions of Separan AP-30 and Polyox WSR-301 in a concentration range of 10–10000 wppm is investigated by means of a cone-and-plate rheogoniometer. The relation between the shear stress and the shear rate is for lower shear rates characterized by a timet ₀, which is concentration dependent. Both polymers show for 4000 s^–1 < < 10000 s^–1 a behaviour similar to that of a Bingham material, characterized by a dynamic viscosity ₀ and an apparent yield stress ₀, which also depend on the concentration. The inertial forces are measured for water and some other Newtonian liquids. An explanation is given why the theoretical model developed for these forces does not match the experimental values; the shape of the liquid surface is shear rate dependent. To obtain the first normal stress difference, we have to correct for these inertial forces, the surface tension and the buoyancy. The normal forces, measured for Separan AP-30, appear to be a linear function of the shear rate for 350 s^–1 < < 3300 s^–1.

---

Zusammenfassung Das rheologische Verhalten wäßriger Polymerlösungen von Separan AP-30 und Polyox WSR-301 wird in einem Konzentrationsgebiet von 10–10000 wppm in einem Kegel-Platte-Rheogoniometer untersucht. Der Zusammenhang zwischen Schubspannung und Schergeschwindigkeit wird für niedrige Schergeschwindigkeiten durch eine konzentrationsabhängige Zeitt ₀ gekennzeichnet. Für Schergeschwindigkeiten 4000 s^–1 < < 10000 s^–1 zeigen beide Polymere ein genähert binghamsches Verhalten, gekennzeichnet durch eine dynamische Viskosität ₀ und eine scheinbare Fließgrenze ₀, welche ebenfalls konzentrationsabhängig sind. Die Trägheitskräfte werden für Wasser und einige newtonsche Öle bestimmt. Die Abweichung der experimentellen Ergebnisse vom theoretischen Modell wird durch die Abhängigkeit der Gestalt der Flüssigkeitsoberfläche von der Schergeschwindigkeit erklärt. Um die Werte der ersten Normalspannungsdifferenz zu erhalten, muß man bezüglich der Trägheitskräfte, der Oberflächenspannung und der Auftriebskräfte korrigieren. Die Normalspannungen für Separan AP-30, gemessen für 350 s^–1 < < 3300 s^–1, zeigen eine lineare Abhängigkeit von der Schergeschwindigkeit.

---

c concentration (wppm) - g acceleration of gravity (ms^–2) - K force (N) - K _b buoyant force (N) - K _c force, acting on the cone (N) - K ₀ dimensional constant def. by eq. [24] (N) - K _s force, def. by eq. [22] (N) - M dimensional constant def. by eq. [24] (Ns) - P _s pressure def. by eq. [17] (Nm^–2) - P ₀ average pressure in the liquid atr = 0 (Nm^–2) - P _R average pressure in the liquid atr = R (Nm^–2) - r ₁,r ₂ radii of curved liquid surface (m) - R platen radius (m) - R _w radius of wetted platen area (m) - S _x standard deviation ofx - t ₀ characteristic time def. by eq. [1] (s) - T temperature (°C) - V volume of the submerged part of the cone (m³) - v tangential velocity of liquid (ms^–1) - x distance (m) - angle (rad) - ₀ cone angle (rad) - calibration constant (Nm^–3) - shear rate (s^–1) - dynamic viscosity (mPa · s) - ₀ viscosity def. by eq. [1] (mPa · s) - contact angle (rad) - density (kgm^–3) - static surface tension (Nm^–1) - shear stress (Nm^–2) - ₀ yield stress def. by eq. [1] (Nm^–2) - _c, _p angular velocity (c = cone,p = plate) (s^–1) With 8 figures and 3 tables     

Simultaneous PIV and pulsed shadow technique in slug flow: a solution for optical problems

A recent technique of simultaneous particle image velocimetry (PIV) and pulsed shadow technique (PST) measurements, using only one black and white CCD camera, is successfully applied to the study of slug flow. The experimental facility and the operating principle are described. The technique is applied to study the liquid flow pattern around individual Taylor bubbles rising in an aqueous solution of glycerol with a dynamic viscosity of 113×10^–3 Pa s. With this technique the optical perturbations found in PIV measurements at the bubble interface are completely solved in the nose and in annular liquid film regions as well as in the rear of the bubble for cases in which the bottom is flat. However, for Taylor bubbles with concave oblate bottoms, some optical distortions appear and are discussed. The measurements achieved a spatial resolution of 0.0022 tube diameters. The results reported show high precision and are in agreement with theoretical and experimental published data.Symbols D internal column diameter (m) - g acceleration due to gravity (m s^–2) - l _w wake length (m) - Q _v liquid volumetric flow rate (m³ s^–1) - r radial position (m) - r * radial position of the wake boundary (m) - R internal column radius (m) - U _s Taylor bubble velocity (m s^–1) - u _z axial component of the velocity (m s^–1) - u _r radial component of the velocity (m s^–1) - z distance from the Taylor bubble nose (m) - Z * distance from the Taylor bubble nose for which the annular liquid film stabilizes (m) Dimensionless groups Re Reynolds number ( ) - N _f inverse viscosity number ( ) Greek letters liquid film thickness (m) - liquid kinematic viscosity (m² s^–1) - liquid dynamic viscosity (Pa s) - liquid density (kg m^–3)     

Two-phase geothermal flows with conduction and the connection with Buckley-Leverett theory

Two-phase flows of boiling water and steam in geothermal reservoirs satisfy a pair of conservation equations for mass and energy which can be combined to yield a hyperbolic wave equation for liquid saturation changes. Recent work has established that in the absence of conduction, the geothermal saturation equation is, under certain conditions, asymptotically identical with the Buckley-Leverett equation of oil recovery theory. Here we summarise this work and show that it may be extended to include conduction. In addition we show that the geothermal saturation wave speed is under all conditions formally identical with the Buckley-Leverett wave speed when the latter is written as the saturation derivative of a volumetric flow.Roman Letters C(P, S,q) geothermal saturation wave speed [ms^–1] (14) - c _t (P, S) two-phase compressibility [Pa^–1] (10) - D(P, S) diffusivity [m s^–2] (8) - E(P, S) energy density accumulation [J m^–3] (3) - g gravitational acceleration (positive downwards) [ms^–2] - h _w (P),h _w (P) specific enthalpies [J kg^–1] - J _M (P, S,P) mass flow [kg m^–2 s^–1] (5) - J _E (P, S,P) energy flow [J m^–2s^–1] (5) - k absolute permeability (constant) [m²] - k _w (S),k _s (S) relative permeabilities of liquid and vapour phases - K formation thermal conductivity (constant) [Wm^–1 K^–1] - L lower sheetC<0 in flow plane - m, c gradient and intercept - M(P, S) mass density accumulation [kg m^–3] (3) - O flow plane origin - P(x,t) pressure (primary dependent variable) [Pa] - q volume flow [ms^–1] (6) - S(x, t) liquid saturation (primary dependent variable) - S ^*(x,t) normalised saturation (Appendix) - t time (primary independent variable) [s] - T temperature (degrees Kelvin) [K] - T _sat(P) saturation line temperature [K] - TdT _sat/dP saturation line temperature derivative [K Pa^–1] (4) - T _c ,T _D convective and diffusive time constants [s] - u _w (P),u _s (P),u _r (P) specific internal energies [J kg^–1] - U upper sheetC > 0 in flow plane - U(x,t) shock velocity [m s^–1] - x spatial position (primary independent variable) [m] - X representative length - x, y flow plane coordinates - z depth variable (+z vertically downwards) [m] Greek Letters _P , _S remainder terms [Pa s^–1], [s^–1] - double-valued saturation region in the flow plane - h =h _s –h _w latent heat [J kg^–1] - = _w – _s density difference [kg m^–3] - line envelope - =D _K /D ₀ diffusivity ratio - porosity (constant) - _w (P), _s (P), _t (P, S) dynamic viscosities [Pa s] - v _w (P),v _s (P) kinematic viscosities [m²s^–1] - v ₀ =kh/KT kinematic viscosity constant [m² s^–1] - ₀ =v ₀ dynamic viscosity constant [m² s^–1] - _w (P), _s (P) density [kg m^–3] Suffixes r rock matrix - s steam (vapour) - w water (liquid) - t total - av average - 0 without conduction - K with conduction     

Rheologische Messungen an Kunststoff-Schmelzen mit einem neuen Rotationsrheometer

Zusammenfassung In der vorliegenden Arbeit wird ein neues Rotationsrheometer vorgestellt und über Messungen an zwei Polymethylmethacrylat-Formmassen berichtet. Bei dem Rheometer handelt es sich um ein Couette-Rheometer mit feststehendem Innenzylinder als Meßkörper. Der Meßkörper ist beidseitig eingespannt. In dem geschlossenen Meßraum können die Schmelzen bis zu einem Druck von 500 bar belastet werden.Der zeitliche Verlauf der Schubspannung in den Schmelzen wird in Abhängigkeit von Temperatur und Druck aufgezeichnet.

---

Summary A new type of rotational rheometer is described, and results for two samples of polymethylmethacrylate are reported. The rheometer consists of a Couette system with fixed inner cylinder, supported at both ends for torque measurements. Pressure may be varied up to 500 bar. Shear stresses have been recorded as a function of time, temperature and pressure.

---

Nomenklatur C [kp cm^–2 s^–1] Steigung der Anlaufkurve im Nullpunkt - D [kp cm rad^–1] Direktionsmoment - E ₀ [kcal mol^–1] Aktivierungsenergie der Newtonschen Viskosität - G [kp cm^–2] Schubmodul - G [—] Griffith-Zahl - l [mm] Länge des Meßkörpers - p [kp cm^–2] Druck - R _i [mm] Radius des Innenzylinders - R _a [mm] Radius des Außenzylinders - t _max [s] Zeit, bei der das Maximum in der Anlaufkurve auftritt - T [°C] Temperatur - ₀ [cm² kp^–1] Druckkoeffizient der Newtonschen Viskosität - [s^–1] Schergeschwindigkeit - ₀ [kp s cm^–2] Newtonsche Viskosität - (g cm²] Trägheitsmoment des Meßkörpers - v ₀ [s^–1] Eigenfrequenz des Meßsystems - _max [kp cm^–2] maximale Schubspannung - _st [kp cm^–2] stationäre Schubspannung Mit 7 Abbildungen und 1 Tabelle     

Influence of Mach and Reynolds numbers on heat transfer to the leeward surface of a half-cone at hypersonic speed

Results are given of an investigation of heat transfer on the flat surface of a blunted half-cone, washed at zero angle of attack by a helium flow at high Mach number (up to 23.5). A correlation is given for the experimental data obtained over a wide range of Mach numbers (M = 3–23.5) and Reynolds numbers (Re_a = 10⁴–3.5·⁵, wherea is the nose radius).Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 105–109, September–October, 1976.     

The wake flow control behind a circular cylinder using ion wind

Active and passive flow control methods have been studied for decades, but there have been only a few studies of flow control methods using ion wind, which is the bulk motion of neutral molecules driven by locally ionized air from a corona discharge. This paper describes an experimental study of ion wind wake control behind a circular cylinder. The experimental conditions consisted of a range of electrohydrodynamic numbers—the ratio of an electrical body force to a fluid inertial force—from 0 to 2 and a range of Reynolds numbers from 4×10³ to 8×10³. Pressure distributions over the cylinder surface were measured and flow visualizations were carried out using a smoke-wire method. The flow visualizations confirmed that ion wind significantly affects the wake structure behind a circular cylinder, and that the pressure drag can be dramatically reduced by superimposing ion wind.List of symbols BR blockage ratio - C _d coefficient of the pressure drag - C _p coefficient of the surface pressure, 2(p–p ₀)/(U ₀ ²) - C _pb coefficient of the base surface pressure, 2(p _b–p ₀)/(U ₀ ²) - D diameter of the cylinder - D _P pressure drag - d _p diameter of particle - E the electric field - F _e Coulombian force (qE) - F _v viscous force - H wire-to-cylinder spacing - I total electric current (A) - L the axial length of cylinder (m) - N _EHD electrohydrodynamic number - p _b base pressure of cylinder at =180° - p ₀ reference static pressure at 10D upstream - q the charge on the particle - R radius of the cylinder - V applied voltage (kV) - U ₀ mean flow velocity (m/s) - ion mobility in air (m²/(s V)) - ₀ permittivity of free space - viscosity of fluid (kg/ms) - density of fluid (kg/m³) - installation angle of a wire electrode (°)     

Numerical analysis of simulation accuracy for hypersonic heat transfer in subsonic jets of dissociated nitrogen

One-dimensional problems of the flow in a boundary layer of finite thickness on the end face of a model and in a thin viscous shock layer on a sphere are solved numerically for three regimes of subsonic flow past a model with a flat blunt face exposed to subsonic jets of pure dissociated nitrogen in an induction plasmatron [1] (for stagnation pressures of (10⁴–3·10⁴) N/m² and an enthalpy of 2.1·10⁷ m²/sec²) and three regimes of hypersonic flow past spheres with parameters related by the local heat transfer simulation conditions [2, 3]. It is established that given equality of the stagnation pressures, enthalpies and velocity gradients on the outer edges of the boundary layers at the stagnation points on the sphere and the model, for a model of radius R_m=1.5·10^–2 m in a subsonic jet the accuracy of reproduction of the heat transfer to the highly catalytic surface of a sphere in a uniform hypersonic flow is about 3%. For surfaces with a low level of catalytic activity the accuracy of simulation of the nonequilibrium heat transfer is determined by the deviations of the temperatures at the outer edges of the boundary layers on the body and the model. For this case the simulation conditions have the form: dU_e/dx=idem, p₀=idem, T_e=idem. At stagnation pressuresP ₀2·10⁴ N/m² irrespective of the catalycity of the surface the heat flux at the stagnation point and the structure of the boundary layer near the axis of symmetry of models with a flat blunt face of radius R_m1.5·10^–2 m exposed to subsonic nitrogen jets in a plasmatron with a discharge channel radius R_c=3·10^–2 m correspond closely to the case of spheres in hypersonic flows with parameters determined by the simulation conditions [2, 3].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 135–143, March–April, 1990.     

Injection moulding of plastics

Summary In continuation of a previous investigation a simple analytical expression is derived in closed form for the thickness distribution of the freeze-off layer which is vitrified at the (flat) wall of an oblong rectangular cavity. As has been pointed out previously, this layer is marked for amorphous polymers by the molecular orientation (birefringence pattern) in the moulded sample. One can show that a more detailed study with the aid of the coupled equations of energy and of motion will not furnish essential improvements. Problems of polymer physics like glass transition or crystallization kinetics at extreme rates of cooling and shearing must be solved first.

---

Zusammenfassung In Fortsetzung einer früheren Untersuchung wurde ein einfacher analytischer Ausdruck in geschlossener Form für die Dickenverteilung der eingefrorenen Schicht abgeleitet, die an der (flachen) Wand eines langgestreckten rechteckigen Formnestes während des Einspritzvorgangs glasig erstarrt. Wie früher auseinandergesetzt wurde, wird diese Schicht bei amorphen Polymeren durch die Molekülorientierung (Doppelbrechungsmuster) im gespritzten Formteil markiert. Man kann zeigen, daß eine eingehendere Studie mit Hilfe der gekoppelten Energie- und Impulsgleichungen keine essentiellen Verbesserungen bringt. Probleme der Polymerphysik, wie Glasübergang oder Kristallisationskinetik bei extremen Abkühlungs- und Schergeschwindigkeiten, müssen erst gelöst werden.

---

List of Symbols a heat diffusivity of polymer melt (averaged overT) [m²s^–1] - B breadth of mould cavity [m] - Br Brinkman number ( ) - c heat capacity of polymer melt (averaged overT) [J kg^–1 K^–1] - F ₀ Fourier number (at _i/4H ²) - h heat transfer coefficient by melt flow [J K^–1 s^–1 m^–2] - h heat transfer coefficient by layer growth [J K^–1 s^–1 m^–2] - H half height of mould cavity [m] - L length of mould cavity [m] - n exponent in eq. [18] (= 0.417) - Nu Nußelt number (2Hh/) - P pressure gradientdP/dz in mould [N m^–3] - t time [s] - t _i injection time [s] - T _g glass transition temperature of polymer [K] - T _i injection temperature of polymer melt [K] - T _l stagnation temperature [K] - T _m mould wall temperature [K] - speed of flow front during mould filling [m s^–1] - x coordinate perpendicular to mould wall [m] - z coordinate in the injection direction [m] - thickness of stagnant layer (atT _l) [m] - ₀ optically detectable freeze-off thickness [m] - ⁺ apparent layer thickness (atT _i) [m] - dimensionless freeze-off thickness (= ₀/2H) - dimensionless distance from entrance (=z/L) - _m dimensionless coordinate of layer maximum - _g dimensionless temperature (= (T _i –T _l)/(T _g –T _m)) - _i dimensionless temperature (= (T _i –T _l)/(T _i –T _m)) - _l dimensionless temperature (= (T _i –T _l)/(T _l –T _m)) - _i viscosity of polymer atT _i [N s m^–3] - _l viscosity of polymer atT _l [N s m^–3] - heat conductivity of polymer melt (averaged) [J K^–1 s^–1 m^–1] - density of polymer melt (averaged) [kg m^–3] - dimensionless time (eq. [11]) - ⁺ dimensionless parameter (eqs. [19a] and [19b]) - dimensionless layer thickness (eq. [12]) - ⁺ dimensionless parameter (eq. [20a]) - dimensionless parameter (eqs. [11a] and [11b]) Formerly at Delft University of Technology, Delft (The Netherlands).Paper presented at the Conference on Chemical Engineering Rheology, Annual Meeting of the Deutsche Rheologische Gesellschaft in Aachen, March 5–7, 1979.With 3 figures and 1 table     

Effect of solids concentration in dense suspensions of silica-iron(III) oxide and water

The rheological properties of dense suspensions, of silica, iron (III) oxide and water, were studied over a range of solids concentrations using a viscometer, which was modified so as to prevent settling of the solid components. Over the conditions studied, the material behaved according to power—law flow relationships. As the concentrations of silica and iron(III) oxide were increased, an entropy term in the flow equation was identified which had a silica dependent and an iron (III) oxide dependent component. This was attributed to a tendency to order into some form of structural regularity. A, A, B, C pre-exponential functions (K Pa^–n s^–1) - C _ox volume fraction iron (III) oxide - Q activation energy (kJ mol^–1) - R gas constant (kJ mol^–1 K^–1) - R _v silica/water volume ratio - T temperature (K) - n power-law index - H enthalpy (kJ mol^–1) - S entropy change (kJ mol^–1 K^–1) - shear strain rate (s^–1) - shear stress (Pa)