Unsteady convective diffusion in a rotating parallel plate channel

The paper presents an exact analysis of the dispersion of a passive contaminant in a viscous fluid flowing in a parallel plate channel driven by a uniform pressure gradient. The channel rotates about an axis perpendicular to its walls with a uniform angular velocity resulting in a secondary flow. Using a generalized dispersion model which is valid for all time, we evaluate the longitudinal dispersion coefficientsK _i (i=1, 2, ...) as functions of time. It is shown thatK ₁=0 andK ₃,K ₄, ... decay rapidly in comparison withK ₂. ButK ₂ decreases with increasing (the dimensionless rotation parameter) for values of upto approximately =2.2. ThereafterK ₂ increases with further increase in and its value gets saturated for large values of (say, 500) and does not change any further with increase in . A physical explanation of this anomalous behaviour ofK ₂ is given.

---

Instationäre konvektive Diffusion in einem rotierenden Parallelplattenkanal

Zusammenfassung In dieser Untersuchung wird eine exakte Analyse der Ausbreitung eines passiven Kontaminierungsstoffes in einer zähen Flüssigkeit gegeben, die, befördert durch einen gleichförmigen Druckgradienten, in einem Parallelplattenkanal strömt. Der Kanal rotiert mit gleichförmiger Winkelgeschwindigkeit um eine zu seinen Wänden senkrechte Achse, wodurch sich eine Sekundärströmung ausbildet. Unter Verwendung eines generalisierten, für alle Zeiten gültigen Dispersionsmodells werden die longitudinalen DispersionskoeffizientenK _i (i=1, 2, ...) als Funktionen der Zeit ermittelt. Es wird gezeigt, daßK ₁=0 gilt und dieK ₃,K ₄, ... gegenüberK ₂ schnell abnehmen.K ₂ nimmt ab, wenn , der dimensionslose Rotationsparameter, bis etwa zum Wert 2,2 ansteigt. Danach wächstK ₂ mit bis auf einem Endwert an, der etwa ab =500 erreicht wird. Dieses anomale Verhalten vonK ₂ findet eine physikalische Erklärung.

---

List of symbols C solute concentration - D molecular diffusivity - K _i longitudinal dispersion coefficients - 2L depth of the channel - P ₀ dimensionless pressure gradient along main flow - Pe Péclet number - q velocity vector - Q _x,Q _y mass flux along the main flow and the secondary flow directions - dimensionless average velocity along the main flow direction - (x, y, z) Cartesian co-ordinates Greek symbols dimensionless rotation parameter - the inclination of side walls withx-axis - kinematic viscosity - fluid density - dimensionless time - angular velocity of the channel - dimensionless distance along the main flow direction - dimensionless distance along the vertical direction - dimensionless solute concentration - integral of the dispersion coefficientK ₂() over a time interval     

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      

Transport in ordered and disordered porous media I: The cellular average and the use of weighting functions

In this work we consider transport in ordered and disordered porous media using singlephase flow in rigid porous mediaas an example. We defineorder anddisorder in terms of geometrical integrals that arise naturally in the method of volume averaging, and we show that dependent variables for ordered media must generally be defined in terms of thecellular average. The cellular average can be constructed by means of a weighting function, thus transport processes in both ordered and disordered media can be treated with a single theory based on weighted averages. Part I provides some basic ideas associated with ordered and disordered media, weighted averages, and the theory of distributions. In Part II a generalized averaging procedure is presented and in Part III the closure problem is developed and the theory is compared with experiment. Parts IV and V provide some geometrical results for computer generated porous media.Roman Letters A interfacial area of the- interface contained within the macroscopic region, m² - A_e area of entrances and exits for the-phase contained within the macroscopic system, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor, m² - L general characteristic length for volume averaged quantities, m - characteristic length (pore scale) for the-phase - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - N unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p0 reference pressure in the-phase, N/m² - p traditional intrinsic volume averaged pressure, N/m² - r0 radius of a spherical averaging volume, m - r position vector, m - r position vector locating points in the-phase, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - v velocity vector in the-phase, m/s - v traditional superficial volume averaged velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V/V, volume average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m²     

Similar solutions and the asymptotics of filtration equations

In this paper we consider the asymptotic behavior of solutions of the quasilinear equation of filtration as t. We prove that similar solutions of the equation u _t = (u )_xx asymptotically represent solutions of the Cauchy problem for the full equation u _t = [(u)]_xx if (u) is close to u for small u.     

The Thermodynamic Significance of the Local Volume Averaged Temperature

Since the temperature is not an additive function, the traditional thermodynamic point of view suggests that the volume integral of the temperature has no precise physical meaning. This observation conflicts with the customary analysis of non-isothermal catalytic reactors, heat pipes, driers, geothermal processes, etc., in which the volume averaged temperature plays a crucial role. In this paper we identify the thermodynamic significance of the volume averaged temperature in terms of a simple two-phase heat transfer process. Given the internal energy as a function of the point temperature and the density

Shear viscosity behavior of highly concentrated suspensions at low and high shear-rates

A method is proposed for determining the shear viscosity behavior of highly concentrated suspensions at low and high shear-rates through the use of a formulation that is a function of three parameters signifying the effects of particle size distribution. These parameters are the intrinsic viscosity [], a parametern that reflects the level of particle association at the initiation of motion and the maximum packing concentration _m. The formulation reduces to the modified Eilers equation withn = 2 for high shear rates. An analytical method was used for the calculation of maximum packing concentration which was subsequently correlated with the experimental values to account for the surface induced interaction of particles with the fluid. The calculated values of viscosities at low and high shear-rates were found to be in good agreement with various experimental data reported in literature. A brief discussion is also offered on the reliability of the methods of measuring the maximum packing concentration. _r = /₀ relative viscosity of the suspension - volumetric concentration of solids - k _n coefficient which characterizes a specific effect of particle interactions - _m maximum packing concentration - _r,0 relative viscosity at low shear-rates - [] intrinsic viscosity - n, n parameter that reflects the level of particle interactions at low and high shear-rates, respectively - _r, relative viscosity at high shear-rates - (_m)_s, (_m)_i, (_m)_l packing factors for small, intermediate and large diameter classes - v _s, v_i, v_l volume fractions of small, intermediate and large diameter classes, respectively - _si, _sl coefficient to be used in relating a smaller to an intermediate and larger particle group, respectively - _is, _il coefficient to be used in relating an intermediate to a smaller and larger particle group, respectively - _ls, _li coefficient to be used in relating a larger to a smaller and intermediate particle group, respectively - _m0 maximum packing concentration for binary mixtures - _m,e measured maximum packing concentration - _m,c calculated maximum packing concentration     

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)     

Measurements of slip at the wall during flow of high-density polyethylene through a rectangular conduit

A hot-film probe has been used to measure slip of high-density polyethylene flowing through a conduit with a rectangular cross section. A transition from no slip to a stick-slip condition has been observed and associated with irregular extrudate shape. Appreciable extrudate roughness was initiated at the same flow rate as that at which the relationship between Nusselt number and Péclet number for heat transfer from the probe departed from the behavior expected for a no-slip condition at the conduit wall. A ₁ constant defined by eq. (A3) - C dimensionless group used in eq. (7) - C _p heat capacity - D constant in eq. (13) - f u _s/u - f _lin defined by eq. (A6) - G storage modulus - G loss modulus - k thermal conductivity - L length of hot film in thex-direction - L _eff effective length of large probe found from eq. (A3) - Nu _L Nusselt number, defined for a lengthL by eq. (2) - (Nu _L)₀ value ofNu _L atPe = 0 (eq. (A 1)) - Pe Péclet number,uL/ - Pe ₀ Péclet number in slip flow, eq. (6) - Pe ₁ Péclet number in shear flow, eq. (4) - q _L average heat flux over hot film of lengthL - R _i resistances defined by figure 8 - R _AB correlation coefficient defined by eq. (14) for signalsA andB - T temperature - T _s temperature of probe surface - T ambient temperature - T T _s –T - u average velocity - u _s slip velocity - V _b voltage indicated in figure 8 - W probe dimension (figure 6) - x distance in flow direction (figure 1) - y distance perpendicular to flow direction (figure 1) - thermal diffusivity,k/C _p - wall shear rate - _5% thickness of lubricating layer during probe calibration for introduction of an error no greater than 5%; see Appendix I - ^* complex viscosity - density - time - _c critical shear stress, eq. (13) - _w wall shear stress - frequency characterizing extrudate distortion (figures 12 and 13), or frequency of oscillation during rheometric characterization (figures 18–20) - * quantity obtained from normalized Nusselt number, eq. (A1), or complex viscosity ^* - A actual (small) probe (see Appendix I) - M model (large) probe (see Appendix I)     

On the closure problem for Darcy's law

In a previous derivation of Darcy's law, the closure problem was presented in terms of an integro-differential equation for a second-order tensor. In this paper, we show that the closure problem can be transformed to a set of Stokes-like equations and we compare solutions of these equations with experimental data. The computational advantages of the transformed closure problem are considerable.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A _e area of entrances and exits for the-phase contained within the averaging volume, m² - B second-order tensor used to respresent the velocity deviation - b vector used to represent the pressure deviation, m^–1 - C second-order tensor related to the permeability tensor, m^–2 - D second-order tensor used to represent the velocity deviation, m² - d vector used to represent the pressure deviation, m - g gravity vector, m/s² - I unit tensor - K C ^–1,–D, Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - l _i i=1, 2, 3, lattice vectors, m - n unit normal vector pointing from the-phase toward the-phase - n _e outwardly directed unit normal vector at the entrances and exits of the-phase - p pressure in the-phase, N/m ² - p intrinsic phase average pressure, N/m² - p –p , spatial deviation of the pressure in the-phase, N/m² - r position vector locating points in the-phase, m - r ₀ radius of the averaging volume, m - t time, s - v velocity vector in the-phase, m/s - v intrinsic phase average velocity in the-phase, m/s - v phase average or Darcy velocity in the \-phase, m/s - v –v , spatial deviation of the velocity in the-phase m/s - V averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ Greek Letters V /V volume fraction of the-phase - mass density of the-phase, kg/m³ - viscosity of the-phase, Nt/m²     

Lp-estimates for the nonlinear spatially homogeneous Boltzmann equation

This paper studies L^p-estimates for solutions of the nonlinear, spatially homogeneous Boltzmann equation. The molecular forces considered include inverse k^th-power forces with k > 5 and angular cut-off.The main conclusions are the following. Let f be the unique solution of the Boltzmann equation with f(v,t)(1 + ¦v²¦)^(s ₁ ^{+ /p)/2} L¹, when the initial value f ₀ satisfies f ₀(v) 0, f ₀(v) (1 + ¦v¦²)^(s ₁ ^{+ /p)/2} L¹, for some s₁ 2 + /p, and f ₀(v) (1 + ¦v¦²)^s/2 L^p. If s 2/p and 1 < p < , then f(v, t)(1 + ¦v¦²)^{(s s} ₁)^/2 L^p, t > 0. If s >2 and 3/(1+ ) < p < , thenf(v,t) (1 + ¦v¦²)^(s(s ₁ ^{+ 3/p))/2} L^p, t > 0. If s >2 + 2C₀/C₁ and 3/(l + ) < p < , then f(v,t)(1 + ¦v¦²)^s/2 L^p, t > 0. Here 1/p + 1/p = 1, x y = min (x, y), and C₀, C₁, 0 < 1, are positive constants related to the molecular forces under consideration; = (k – 5)/ (k – 1) for k^th-power forces.Some weaker conclusions follow when 1 < p 3/ (1 + ).In the proofs some previously known L-estimates are extended. The results for L^p, 1 < p < , are based on these L-estimates coupled with nonlinear interpolation.     

Parameter identification in a soil with constant diffusivity

We consider infiltration into a soil that is assumed to have hydraulic conductivity of the form K = K = K_se^h and water content of the form = K – _r. Here h denotes capillary pressure head while K_s, , and _r represent soil specific parameters. These assumptions linearize the flow equation and permit a closed form solution that displays the roles of all the parameters appearing in the hydraulic function K and . We assume K_s and _r to be known. A measurement of diffusivity fixes the product of and resulting in a parameter identification problem for one parameter. We show that this parameter identification problem, in some cases, has a unique solution. We also show that, in some cases, this parameter identification problem can have multiple solutions, or no solution. In addition it is shown that solutions to the parameter identification problem can be very sensitive to small changes in the problem data.     

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)     

On the effects of a periodic, spanwise variation of suction upon asymptotic suction flow

Summary The effects of superposing streamwise vorticity, periodic in the lateral direction, upon two-dimensional asymptotic suction flow are analyzed. Such vorticity, generated by prescribing a spanwise variation in the suction velocity, is known to play an important role in unstable and turbulent boundary layers. The flow induced by the variation has been obtained for a freestream velocity which (i) is steady, (ii) oscillates periodically in time, (iii) changes impulsively from rest. For the oscillatory case it is shown that a frequency can exist which maximizes the induced, unsteady wall shear stress for a given spanwise period. For steady flow the heat transfer to, or from a wall at constant temperature has also been computed.Nomenclature (x, y, z) spatial coordinates - (u, v, w) corresponding components of velocity - (, , ) corresponding components of vorticity - t time - stream function for v and w - v _w mean wall suction velocity - nondimensional amplitude of variation in wall suction velocity - characteristic wavenumber for variation in direction of z - T temperature - P pressure - density - coefficient of kinematic viscosity - coefficient of thermal diffusivity - (/v _w)² - frequency of oscillation of freestream velocity - nondimensional amplitude of freestream oscillation - /v _w ² - z z - y –v _w y/ - v _w ² t/4 - /v _w - U ₀ characteristic freestream velocity - u/U ₀ - coefficient of viscosity - _w wall shear stress - Prandtl number (/) - q heat transfer to wall - T _w wall temperature - T (T _w–T)/(T _w–)     

Elongational flow behavior of viscoelastic liquids: Modelling bubble dynamics with viscoelastic constitutive relations

Summary Earlier parts of this series have described a technique based on the collapse of single bubbles in the fluids for studying the elongational rheology of viscoelastic solutions and melts of moderate viscosities ( ₀ > 10²p) at relatively high strain rates . The present paper describes the modelling of bubble collapse with both rate and integral type constitutive relations using a body coordinate system. Predictions of the stress at the bubble wall as a function of time during collapse from a BKZ model and a modified corotational Maxwell model compared favorably with experimental data for two polymer solutions, 1% polyacrylamide in water/glycerine and 2% hydroxypropyl cellulose in water.

---

Zusammenfassung In vorangehenden Veröffentlichungen dieser Reihe wurde eine Methode beschrieben, mit Hilfe derer man aus dem Zerfall von einzelnen Blasen in einer Flüssigkeit auf die Dehn-Rheologie viskoelastischer Lösungen und Schmelzen mittlerer Viskosität ( ₀ > 10² P) bei relativ hohen Dehngeschwindigkeiten schließen kann. Die vorliegende Untersuchung beschreibt Modelle des Blasenzerfalls mit Hilfe von Stoffgleichungen sowohl vom rate- als auch vom Integral-Typ, wobei ein körperfestes Koordinatensystem benutzt wird. Die Voraussagen der Spannung an der Blasenwand als Funktion der Zeit während des Zerfalls bei Verwendung eines BKZ- und eines modifizierten korotatorischen Maxwell-Modells zeigen eine recht gute Übereinstimmung mit experimentellen Werten, die an zwei Polymerlösungen, nämlich einer 1%igen Polyacrylamid-Lösung in einer Wasser-Glycerin-Mischung und einer 2%igen wäßrigen Hydropropylcellulose, erhalten worden sind.

---

Nomenclature a material constant - b material constant - g metric tensor, space coordinates - m material constant - n material constant - p pressure - P _G pressure within bubble - P _R pressure outside bubble at the wall - P pressure far away from the bubble - R bubble radius - dR/dt - R ₀ initial bubble radius - t time - u velocity - U potential function - Y R/R ₀ Greek symbols covariant body metric tensor - surface tension - rate of deformation matrix, II -second invariant of - strain rate - ₀ zero shear rate viscosity - _e elongational viscosity - _ef effective viscosity - ¹, ², ³ coordinates in body system - ^{1 1}/R ₀ ³ - body stress tensor - density - space stress tensor - relaxation time - _ef effective relaxation time - bubble pressure function, defined in eq. [19] - vorticity tensor With 11 figures and 1 table     

Deformation to breaking of deep water gravity waves

The results of laboratory observations of the deformation of deep water gravity waves leading to wave breaking are reported. The specially developed visualization technique which was used is described. A preliminary analysis of the results has led to similar conclusions than recently developed theories. As a main fact, the observed wave breaking appears as the result of, first, a modulational instability which causes the local wave steepness to approach a maximum and, second, a rapidly growing instability leading directly to the breaking.List of symbols L total wave length - H total wave height - crest elevation above still water level - trough depression below still water level - wave steepness =H/L - crest steepness =/L - trough steepness =/L - F ₁ forward horizontal length from zero-upcross point (A) to wave crest - F ₂ backward horizontal length from wave crest to zero-downcross point (B) - crest front steepness =/F ₁ - crest rear steepness =/F ₂ - vertical asymmetry factor=F ₂/F ₁ (describing the wave asymmetry with respect to a vertical axis through the wave crest) - µ horizontal asymmetry factor=/H (describing the wave asymmetry with respect to a horizontal axis: SWL) - T ₀ wavemaker period - L ₀ theoretical wave length of a small amplitude sinusoïdal wave generated at T _inf0 ^sup–1 frequency - ₀ average wave height     

Influence of wall conduction on natural convection in an inclined square enclosure

A numerical study has been made of natural convection in an externally heated inclined enclosure with finitely conducting side walls. Results indicate that conduction along the enclosure walls has a stabilizing influence onthe convective motion in the enclosure and is therefore responsible for reduced heat transfer from the enclosure. The average Nusselt number along the hot and cold surfaces is observed to decrease with decreasing conductivity ratio and increasing wall thickness. The heat flux along the side walls is observed to reverse its direction and therefore, the Nusselt number along the side wall interfaces can be either positive or negative.

---

Einfluß der Wärmeleitung in der Wand auf die freie Konvektion in einem geneigten quadratischen Raum

Zusammenfassung Es wurde eine numerische Studie der freien Konvektion in einem von außen beheizten geneigten Raum mit endlicher Wärmeleitung in den Seitenwänden durchgeführt. Die Ergebnisse zeigen, daß die Leitung längs der Raumwände einen stabilisierenden Einfluß auf die Konvektionsbewegung in dem Raum hat und deshalb für die Verringerung des Wärmeübergangs aus dem Raum verantwortlich ist. Es wurde beobachtet, daß die mittlere Nusselt-Zahl längs der warmen und kalten Oberflächen mit abnehmendem Wärmeleitverhältnis und zunehmender Wandstärke geringer wird. Es wurde auch beobachtet, daß der Wärmestrom längs der Seitenwände seine Richtung umkehrt und deshalb die Nusselt-Zahl an den Oberflächen längs der Seitenwand entweder positiv oder negativ sein kann.

---

Nomenclature g gravitational acceleration - H inside height of the enclosure - k thermal conductivity of fluid - k _w thermal conductivity of enclosure walls - k _r k _w/k/Pr - L inside width of the enclosure - n coordinate direction normal to an interface - Nu Nusselt number - Nusselt number along the cold wall interface and its average value - Nusselt number along the hot wall interface and its average value - Nusselt number along the left wall interface and its average value - Nusselt number along the right wall interface and its average value - p, p*, P thermodynamic, modified and dimensionless pressures - Pr Prandtl number - Ra Rayleigh number,Ra=g (T _h – T _c ) L ³ / - T temperature - T _h ,T _c temperature of the hot and cold surfaces - T ₀ T _h + T _c )/2 - u, U dimensional and dimensionless velocity in thex-direction - U dimensionless vector velocity,i U+jV - , V dimensional and dimensionless velocity in they-direction - x,y dimensional coordinates,x along the heated and cooled walls andy along adiabatic walls - X,Y dimensionless coordinates,X=x/L, Y=y/L Greek symbols thermal diffusivity,k/( c _p ) - thermal expansion coefficient - diffusion coefficient for a general dependent variable - dimensionless fluid thermal conductivity - _s dimensionless solid thermal conductivity - kinematic viscosity - density - ₀ reference density (atT ₀) - general dependent variable - dimensionless temperature in fluid - _s dimensionless temperature in solid - angle of inclination     

The generation of turbulence from displaced cross-members in uniform flow

This paper presents details of an experimental investigation into the nature of turbulence generated in the wake of a single grid node. The latter has been considered as two members placed perpendicular to each other in the geometry of a cross, with square and circular sections representing bars and rods respectively. The effect of member spacing has been examined in an attempt to identify the complex flow phenomena associated with such a configuration, and in this respect a critical gap width has been found.List of symbols C _p pressure coefficient, (p — p ₀)/1/2 U ₀ ² - C _pb pressure coefficient measured on the base centre-line - C _ps pressure coefficient measured at stagnation point - D diameter/section depth of model - L distance between central axis of two cylinders or bars - n vortex shedding frequency - p local pressure on model's surface - p ₀ static pressure - R () autocorrelation coefficient - R _e Reynolds number, DU ₀/v - St Strouhal number, n D/U ₀ - U ₀ mean freestream velocity in X-direction - ovu local mean velocity in X-direction - u velocity fluctuation in X-direction - ovv local mean velocity in Y-direction - ovw local mean velocity in Z-direction - X cartesian co-ordinate in longitudinal direction - Y cartesian co-ordinate perpendicular to wind tunnel floor - Z cartesian co-ordinate in lateral direction - dynamic viscosity of fluid - v kinematic viscosity of fluid, / - _n spectral energy, 4 ₀ R () cos 2 n d - density of fluid - _x longitudinal component of vorticity, /y – /z This paper was presented at the 10th Symposium on Turbulence, University of Missouri-Rolla, Sept. 22–24, 1986     

The finite element solutions of laminar flow and combined convection of air in a staggered or an in-line tube-bank

A finite element method is used to solve the full Navier-Stokes and energy equations for the problems of laminar combined convection from three isothermal heat horizontal cylinders in staggered tube-bank and four isothermal heat horizontal cylinders in in-line tube-bank. The variations of surface shear stress, pressure and Nusselt number are obtained over the entire cylinder surface including the zone beyond the separation point. The predicted values of total, pressure and friction drag coefficients, average Nusselt number and the plots of velocity flow fields and isotherms are also presented.

---

Die Finite-Elemente-Lösung von laminarer Strömung und kombinierter Konvektion von Luft in einer versetzten oder fluchtenden Rohranordnung

Zusammenfassung Eine Methode der finiten Elemente wird zur Lösung der vollständigen Navier-Stokes- und der Energiegleichung für die Probleme der laminaren kombinierten Konvektion an drei isothermen geheizten horizontalen Zylindern in versetzter Rohranordnung sowie für vier isotherme geheizte horizontale Zylinder in fluchtender Anordnung verwendet.Die Veränderung der Wandschubspannung, des Druckes und der Nusselt-Zahl werden für die gesamte Zylinderoberfläche, einschließlich des Bereiches nach dem Ablösepunkt, bestimmt. Die Werte des gesamten Widerstandsbeiwertes aufgrund von Druck und Reibung, die durchschnittliche Nusselt-Zahl und die Diagramme des Geschwindigkeitsfeldes und der Isothermen werden ebenfalls aufgezeigt.

---

Nomenclature C specifie heat - C _D total drag coefficient - C _f friction drag coefficient - C _p pressure drag coefficient - D diameter of cylinder,L=2R ₀ - G, g gravitational acceleration - Gr Grashof number, g(T_w–T )D ³/v ² - h local heat transfer coefficient - K thermal conductivity - L spacing between the centers of cylinder - M _l shape function - N _i shape function - Nu, local and average Nusselt numbers - P dimensionless pressure, p^*/u ² - p ^*,p pressure, free stream pressure - Pe Peclet number,RePr - Pr Prandtl number, c/K - Ra Rayleigh number,Gr Pr - Re Reynolds number,Du /v - R ₀ radius of cylinder - T temperature - T _w temperature on cylinder surface with fixed value - T free stream temperature - v dimensionless x-direction component of velocity,v ^*/u - u ^* x-direction component of velocity - u free stream velocity - v dimensionless Y-direction component of velocity,v ^*/u - v ^* Y-direction component of velocity - X x-direction axis - x dimensionless x-direction coordinate,x ^*/D - x^* x-direction coordinate - Y Y-direction axis - y dimensionless Y-direction coordinate,y ^*/D - y ^* Y-direction coordinate Greek symbols coefficient of volumetric thermal expansion - plane angle - dynamic viscosity - kinematic viscosity, / - density of fluid - _w dimensionless surface shear stress, _* ^w /u ² - skw/* surface shear stress - dimensionless temperature,      

Heat transfer from laminar flow in ducts with elliptic cross-section

Summary The cooling of a hot fluid in laminar Newtonian flow through cooled elliptic tubes has been calculated theoretically. Numerical data have been computed for the two values 1.25 and 4 of the axial ratio of the elliptic cross-section . For =1.25 the influence of non-zero thermal resistance between outmost fluid layer and isothermal surroundings has also been investigated. Special attention has been given to the distribution of heat flux around the perimeter; when increases the flux varies more with the position at the circumference. This positional dependence becomes less pronounced, however, as the (position-independent) thermal resistance of the wall increases.Flattening of the conduit, while maintaining its cross-sectional area constant, improves the cooling. Comparison with rectangular pipes shows that this improvement is not as marked with elliptic as with rectangular pipes.Nomenclature A _k =A _{m, n} coefficients of expansion (6) - a, b half-axes of ellipse, b<a - a _p =a _{r, s} coefficients of representation (V) - D hydraulic diameter, = 4S/P; S = cross-sectional area, P = perimeter - D _e equivalent diameter, according to (13) - n coordinate (outward) normal to the tube wall - T temperature of fluid - T _i temperature of fluid at the inlet - T _s temperature of surroundings - v ₀ mean velocity of fluid - v _z longitudinal velocity of fluid - x, y carthesian coordinates coinciding with axes of ellipse - z coordinate in flow direction - , dimensionless half-axes of ellipse, =a/D and =b/D - _t heat transfer coefficient from fluid at bulk temperature to surroundings; equation (11) - _w heat transfer coefficient at the wall; equation (3) - axial ratio of ellipse, = a/b = / - , , , dimensionless coordinates; =x/D, =y/D, =z/D, =n/D - dimensionless temperature, = (T–T _s)/(T _i–T _s) - ₀ cup-mixing mean value of ; equation (10) - thermal conductivity of fluid - _m,n = _k eigenvalue - c volumetric heat capacity of fluid - _{m, n} = _k = _k eigenfunction; equations (6) and (I) - Nu total Nusselt number, = _t D/ - Nusselt number at large distance from the inlet - Nu _w wall Nusselt number, = _w D/, based on _w - Pé Péclet number, = ₀ Dc/     

Flow visualization of compressible vortex structures using density gradient techniques

Mathematical results are derived for the schlieren and shadowgraph contrast variation due to the refraction of light rays passing through two-dimensional compressible vortices with viscous cores. Both standard and small-disturbance solutions are obtained. It is shown that schlieren and shadowgraph produce substantially different contrast profiles. Further, the shadowgraph contrast variation is shown to be very sensitive to the vortex velocity profile and is also dependent on the location of the peak peripheral velocity (viscous core radius). The computed results are compared to actual contrast measurements made for rotor tip vortices using the shadowgraph flow visualization technique. The work helps to clarify the relationships between the observed contrast and the structure of vortical structures in density gradient based flow visualization experiments.Nomenclature a Unobstructed height of schlieren light source in cutoff plane, m - c Blade chord, m - f Focal length of schlieren focusing mirror, m - C _T Rotor thrust coefficient, T/( ² R ⁴) - I Image screen illumination, Lm/m ² - l Distance from vortex to shadowgraph screen, m - n _b Number of blades - p Pressure,N/m ² - p Ambient pressure, N/m ² - r, , z Cylindrical coordinate system - r _c Vortex core radius, m - Non-dimensional radial coordinate, (r/r _c ) - R Rotor radius, m - Tangential velocity, m/s - Specific heat ratio of air - Circulation (strength of vortex), m ²/s - Non-dimensional quantity, ² 8²p r _c ² - Refractive index of fluid medium - ₀ Refractive index of fluid medium at reference conditions - Gladstone-Dale constant, m ³/kg - Density, kg/m ³ - Density at ambient conditions, kg/m ³ - Non-dimensional density, (/ ) - Rotor solidity, (n _b c/ R) - Rotor rotational frequency, rad/s